The UNC-Wisconsin Rhesus Macaque Neurodevelopment Database: A Structural MRI and DTI Database of Early Postnatal Development

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.

nBEST: Deep-learning-based non-human primates Brain Extraction and Segmentation Toolbox across ages, sites and species

Accurate processing and analysis of non-human primate (NHP) brain magnetic resonance imaging (MRI) serves an indispensable role in understanding brain evolution, development, aging, and diseases. Despite the accumulation of diverse NHP brain MRI datasets at various developmental stages and from various imaging sites/scanners, existing computational tools designed for human MRI typically perform poor on NHP data, due to huge differences in brain sizes, morphologies, and imaging appearances across species, sites, and ages, highlighting the imperative for NHP-specialized MRI processing tools. To address this issue, in this paper, we present a robust, generic, and fully automated computational pipeline, called non-human primates Brain Extraction and Segmentation Toolbox (nBEST), whose main functionality includes brain extraction, non-cerebrum removal, and tissue segmentation. Building on cutting-edge deep learning techniques by employing lifelong learning to flexibly integrate data from diverse NHP populations and innovatively constructing 3D U-NeXt architecture, nBEST can well handle structural NHP brain MR images from multi-species, multi-site, and multi-developmental-stage (from neonates to the elderly). We extensively validated nBEST based on, to our knowledge, the largest assemblage dataset in NHP brain studies, encompassing 1,469 scans with 11 species (e.g., rhesus macaques, cynomolgus macaques, chimpanzees, marmosets, squirrel monkeys, etc.) from 23 independent datasets. Compared to alternative tools, nBEST outperforms in precision, applicability, robustness, comprehensiveness, and generalizability, greatly benefiting downstream longitudinal, cross-sectional, and cross-species quantitative analyses. We have made nBEST an open-source toolbox (https://github.com/TaoZhong11/nBEST) and we are committed to its continual refinement through lifelong learning with incoming data to greatly contribute to the research field.

Effects of multiple anesthetic exposures on rhesus macaque brain development: a longitudinal structural MRI analysis

Concerns about the potential neurotoxic effects of anesthetics on developing brain exist. When making clinical decisions, the timing and dosage of anesthetic exposure are critical factors to consider due to their associated risks. In our study, we investigated the impact of repeated anesthetic exposures on the brain development trajectory of a cohort of rhesus monkeys (n = 26) over their first 2 yr of life, utilizing longitudinal magnetic resonance imaging data. We hypothesized that early or high-dose anesthesia exposure could negatively influence structural brain development. By employing the generalized additive mixed model, we traced the longitudinal trajectories of brain volume, cortical thickness, and white matter integrity. The interaction analysis revealed that age and cumulative anesthetic dose were variably linked to white matter integrity but not to morphometric measures. Early high-dose exposure was associated with increased mean, axial, and radial diffusivities across all white matter regions, compared to late-low-dose exposure. Our findings indicate that early or high-dose anesthesia exposure during infancy disrupts structural brain development in rhesus monkeys. Consequently, the timing of elective surgeries and procedures that require anesthesia for children and pregnant women should be strategically planned to account for the cumulative dose of volatile anesthetics, aiming to minimize the potential risks to brain development.

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

Postnatal growth and spatial conformity of the cranium, brain, eyeballs and masseter muscles in the macaque (Macaca mulatta)

Spatial growth constraints in the head region can lead to coordinated patterns of morphological variation that pleiotropically modify genetically defined phenotypes as the tissues compete for space. Here we test for such architectural modifications during rhesus macaque (Macaca mulatta) postnatal ontogeny. We captured cranium and brain shape from 153 MRI datasets spanning 13 to 1090 postnatal days and tested for patterns of covariation with measurements of relative brain, eyeball, and masseter muscle size as well as callosal tract length. We find that the shape of the infant (<365 days) macaque cranium was most closely aligned to masseter muscle and brain size measured relative to face size. Infant brain and juvenile (365-1090 days) cranium shape were more closely linked with brain size relative to basicranium and face size. Meanwhile, the juvenile macaque brain shape was dominated by the size of the brain relative to that of the basicranium. Associations with relative eyeball size and commissural tract lengths were weaker. Our results are consistent with a spatial-packing regime operating during postnatal macaque ontogeny, in which relative growth of the masseter, face and basicranium have a greater influence than brain growth on the overall shape of the cranium and brain.

Longitudinal mapping of the development of cortical thickness and surface area in rhesus macaques during the first three years

Studying dynamic spatiotemporal patterns of early brain development in macaque monkeys is critical for understanding the cortical organization and evolution in humans, given the phylogenetic closeness between humans and macaques. However, due to huge challenges in the analysis of early brain Magnetic Resonance Imaging (MRI) data typically with extremely low contrast and dynamic imaging appearances, our knowledge of the early macaque cortical development remains scarce. To fill this critical gap, this paper characterizes the early developmental patterns of cortical thickness and surface area in rhesus macaques by leveraging advanced computing tools tailored for early developing brains based on a densely sampled longitudinal dataset with 140 rhesus macaque MRI scans seamlessly covering from birth to 36 mo of age. The average cortical thickness exhibits an inverted U-shaped trajectory with peak thickness at around 4.3 mo of age, which is remarkably in line with the age of peak thickness at 14 mo in humans, considering the around 3:1 age ratio of human to macaque. The total cortical surface area in macaques increases monotonically but with relatively lower expansions than in humans. The spatial distributions of thicker and thinner regions are quite consistent during development, with gyri having a thicker cortex than sulci. By 4 mo of age, over 81% of cortical vertices have reached their peaks in thickness, except for the insula and medial temporal cortices, while most cortical vertices keep expanding in surface area, except for the occipital cortex. These findings provide important insights into early brain development and evolution in primates.

Effects of anesthesia exposure on postnatal maturation of white matter in rhesus monkeys

There is increasing concern about the potential effects of anesthesia exposure on the developing brain. The effects of relatively brief anesthesia exposures used repeatedly to acquire serial magnetic resonance imaging scans could be examined prospectively in rhesus macaques. We analyzed magnetic resonance diffusion tensor imaging (DTI) of 32 rhesus macaques (14 females, 18 males) aged 2 weeks to 36 months to assess postnatal white matter (WM) maturation. We investigated the longitudinal relationships between each DTI property and anesthesia exposure, taking age, sex, and weight of the monkeys into consideration. Quantification of anesthesia exposure was normalized to account for variation in exposures. Segmented linear regression with two knots provided the best model for quantifying WM DTI properties across brain development as well as the summative effect of anesthesia exposure. The resulting model revealed statistically significant age and anesthesia effects in most WM tracts. Our analysis indicated there were major effects on WM associated with low levels of anesthesia even when repeated as few as three times. Fractional anisotropy values were reduced across several WM tracts in the brain, indicating that anesthesia exposure may delay WM maturation, and highlight the potential clinical concerns with even a few exposures in young children.

Gyral peaks: Novel gyral landmarks in developing macaque brains

Cerebral cortex development undergoes a variety of processes, which provide valuable information for the study of the developmental mechanism of cortical folding as well as its relationship to brain structural architectures and brain functions. Despite the variability in the anatomy-function relationship on the higher-order cortex, recent studies have succeeded in identifying typical cortical landmarks, such as sulcal pits, that bestow specific functional and cognitive patterns and remain invariant across subjects and ages with their invariance being related to a gene-mediated proto-map. Inspired by the success of these studies, we aim in this study at defining and identifying novel cortical landmarks, termed gyral peaks, which are the local highest foci on gyri. By analyzing data from 156 MRI scans of 32 macaque monkeys with the age spanned from 0 to 36 months, we identified 39 and 37 gyral peaks on the left and right hemispheres, respectively. Our investigation suggests that these gyral peaks are spatially consistent across individuals and relatively stable within the age range of this dataset. Moreover, compared with other gyri, gyral peaks have a thicker cortex, higher mean curvature, more pronounced hub-like features in structural connective networks, and are closer to the borders of structural connectivity-based cortical parcellations. The spatial distribution of gyral peaks was shown to correlate with that of other cortical landmarks, including sulcal pits. These results provide insights into the spatial arrangement and temporal development of gyral peaks as well as their relation to brain structure and function.

Imaging the fetal nonhuman primate brain with SV2A positron emission tomography (PET)

PURPOSE

Exploring synaptic density changes during brain growth is crucial to understanding brain development. Previous studies in nonhuman primates report a rapid increase in synapse number between the late gestational period and the early neonatal period, such that synaptic density approaches adult levels by birth. Prenatal synaptic development may have an enduring impact on postnatal brain development, but precisely how synaptic density changes in utero are unknown because current methods to quantify synaptic density are invasive and require post-mortem brain tissue.

METHODS

We used synaptic vesicle glycoprotein 2A (SV2A) positron emission tomography (PET) radioligands [¹¹C]UCB-J and [¹⁸F]Syn-VesT-1 to conduct the first assessment of synaptic density in the developing fetal brain in gravid rhesus monkeys. Eight pregnant monkeys were scanned twice during the third trimester at two imaging sites. Fetal post-mortem samples were collected near term in a subset of subjects to quantify SV2A density by Western blot.

RESULTS

Image-derived fetal brain SV2A measures increased during the third trimester. SV2A concentrations were greater in subcortical regions than in cortical regions at both gestational ages. Near term, SV2A density was higher in primary motor and visual areas than respective associative regions. Post-mortem quantification of SV2A density was significantly correlated with regional SV2A PET measures.

CONCLUSION

While further study is needed to determine the exact relationship of SV2A and synaptic density, the imaging paradigm developed in the current study allows for the effective in vivo study of SV2A development in the fetal brain.

Tracing Modification to Cortical Circuits in Human and Nonhuman Primates from High-Resolution Tractography, Transcription, and Temporal Dimensions

The neural circuits that support human cognition are a topic of enduring interest. Yet, there are limited tools available to map brain circuits in the human and nonhuman primate brain. We harnessed high-resolution diffusion MR tractography, anatomic, and transcriptomic data from individuals of either sex to investigate the evolution and development of frontal cortex circuitry. We applied machine learning to RNA sequencing data to find corresponding ages between humans and macaques and to compare the development of circuits across species. We transcriptionally defined neural circuits by testing for associations between gene expression and white matter maturation. We then considered transcriptional and structural growth to test whether frontal cortex circuit maturation is unusually extended in humans relative to other species. We also considered gene expression and high-resolution diffusion MR tractography of adult brains to test for cross-species variation in frontal cortex circuits. We found that frontal cortex circuitry development is extended in primates, and concomitant with an expansion in corticocortical pathways compared with mice in adulthood. Importantly, we found that these parameters varied relatively little across humans and studied primates. These data identify a surprising collection of conserved features in frontal cortex circuits across humans and Old World monkeys. Our work demonstrates that integrating transcriptional and structural data across temporal dimensions is a robust approach to trace the evolution of brain pathways in primates.SIGNIFICANCE STATEMENT Diffusion MR tractography is an exciting method to explore pathways, but there are uncertainties in the accuracy of reconstructed tracts. We broaden the repertoire of toolkits to enhance our ability to trace human brain pathways from diffusion MR tractography. Our integrative approach finds corresponding ages across species and transcriptionally defines neural circuits. We used this information to test for variation in circuit maturation across species and found a surprising constellation of similar features in frontal cortex neural circuits across humans and primates. Integrating across scales of biological organization expands the repertoire of tools available to study pathways in primates, which opens new avenues to study pathways in health and diseases of the human brain.

Cortical development coupling between surface area and sulcal depth on macaque brains

Postnatal development of cerebral cortex is associated with a variety of neuronal processes and is thus critical to development of brain function and cognition. Longitudinal changes of cortical morphology and topology, such as postnatal cortical thinning and flattening have been widely studied. However, thorough and systematic investigation of such cortical change, including how to quantify it from multiple spatial directions and how to relate it to surface topology, is rarely found. In this work, based on a longitudinal macaque neuroimaging dataset, we quantified local changes in gyral white matter's surface area and sulcal depth during early development. We also investigated how these two metrics are coupled and how this coupling is linked to cortical surface topology, underlying white matter, and positions of functional areas. Semi-parametric generalized additive models were adopted to quantify the longitudinal changes of surface area (A) and sulcal depth (D), and the coupling patterns between them. This resulted in four classes of regions, according to how they change compared with global change throughout early development: slower surface area change and slower sulcal depth change (slowA_slowD), slower surface area change and faster sulcal depth change (slowA_fastD), faster surface area change and slower sulcal depth change (fastA_slowD), and faster surface area change and faster sulcal depth change (fastA_fastD). We found that cortex-related metrics, including folding pattern and cortical thickness, vary along slowA_fastD-fastA_slowD axis, and structural connection-related metrics vary along fastA_fastD-slowA_slowD axis, with which brain functional sites align better. It is also found that cortical landmarks, including sulcal pits and gyral hinges, spatially reside on the borders of the four patterns. These findings shed new lights on the relationship between cortex development, surface topology, axonal wiring pattern and brain functions.

Scan Once, Analyse Many: Using Large Open-Access Neuroimaging Datasets to Understand the Brain

We are now in a time of readily available brain imaging data. Not only are researchers now sharing data more than ever before, but additionally large-scale data collecting initiatives are underway with the vision that many future researchers will use the data for secondary analyses. Here I provide an overview of available datasets and some example use cases. Example use cases include examining individual differences, more robust findings, reproducibility-both in public input data and availability as a replication sample, and methods development. I further discuss a variety of considerations associated with using existing data and the opportunities associated with large datasets. Suggestions for further readings on general neuroimaging and topic-specific discussions are also provided.

Comparative connectomics of the primate social brain

Social interaction is thought to provide a selection pressure for human intelligence, yet little is known about its neurobiological basis and evolution throughout the primate lineage. Recent advances in neuroimaging have enabled whole brain investigation of brain structure, function, and connectivity in humans and non-human primates (NHPs), leading to a nascent field of comparative connectomics. However, linking social behavior to brain organization across the primates remains challenging. Here, we review the current understanding of the macroscale neural mechanisms of social behaviors from the viewpoint of system neuroscience. We first demonstrate an association between the number of cortical neurons and the size of social groups across primates, suggesting a link between neural information-processing capacity and social capabilities. Moreover, by capitalizing on recent advances in species-harmonized functional MRI, we demonstrate that portions of the mirror neuron system and default-mode networks, which are thought to be important for representation of the other's actions and sense of self, respectively, exhibit similarities in functional organization in macaque monkeys and humans, suggesting possible homologies. With respect to these two networks, we describe recent developments in the neurobiology of social perception, joint attention, personality and social complexity. Together, the Human Connectome Project (HCP)-style comparative neuroimaging, hyperscanning, behavioral, and other multi-modal investigations are expected to yield important insights into the evolutionary foundations of human social behavior.

Computational models of cortical folding: A review of common approaches

The process of gyrification, by which the brain develops the intricate pattern of gyral hills and sulcal valleys, is the result of interactions between biological and mechanical processes during brain development. Researchers have developed a vast array of computational models in order to investigate cortical folding. This review aims to summarize these studies, focusing on five essential elements of the brain that affect development and gyrification and how they are represented in computational models: (i) the constraints of skull, meninges, and cerebrospinal fluid; (ii) heterogeneity of cortical layers and regions; (iii) anisotropic behavior of subcortical fiber tracts; (iv) material properties of brain tissue; and (v) the complex geometry of the brain. Finally, we highlight areas of need for future simulations of brain development.

Strengths and challenges of longitudinal non-human primate neuroimaging

Longitudinal non-human primate neuroimaging has the potential to greatly enhance our understanding of primate brain structure and function. Here we describe its specific strengths, compared to both cross-sectional non-human primate neuroimaging and longitudinal human neuroimaging, but also its associated challenges. We elaborate on factors guiding the use of different analytical tools, subject-specific versus age-specific templates for analyses, and issues related to statistical power.

CIVET-Macaque: An automated pipeline for MRI-based cortical surface generation and cortical thickness in macaques

The MNI CIVET pipeline for automated extraction of cortical surfaces and evaluation of cortical thickness from in-vivo human MRI has been extended for processing macaque brains. Processing is performed based on the NIMH Macaque Template (NMT), as the reference template, with the anatomical parcellation of the surface following the D99 and CHARM atlases. The modifications needed to adapt CIVET to the macaque brain are detailed. Results have been obtained using CIVET-macaque to process the anatomical scans of the 31 macaques used to generate the NMT and another 95 macaques from the PRIME-DE initiative. It is anticipated that the open usage of CIVET-macaque will promote collaborative efforts in data collection and processing, sharing, and automated analyses from which the non-human primate brain imaging field will advance.

General anaesthesia during infancy reduces white matter micro-organisation in developing rhesus monkeys

BACKGROUND

Non-human primates are commonly used in neuroimaging research for which general anaesthesia or sedation is typically required for data acquisition. In this analysis, the cumulative effects of exposure to ketamine, Telazol® (tiletamine and zolazepam), and the inhaled anaesthetic isoflurane on early brain development were evaluated in two independent cohorts of typically developing rhesus macaques.

METHODS

Diffusion MRI scans were analysed from 43 rhesus macaques (20 females and 23 males) at either 12 or 18 months of age from two separate primate colonies.

RESULTS

Significant, widespread reductions in fractional anisotropy with corresponding increased axial, mean, and radial diffusivity were observed across the brain as a result of repeated anaesthesia exposures. These effects were dose dependent and remained after accounting for age and sex at time of exposure in a generalised linear model. Decreases of up to 40% in fractional anisotropy were detected in some brain regions.

CONCLUSIONS

Multiple exposures to commonly used anaesthetics were associated with marked changes in white matter microstructure. This study is amongst the first to examine clinically relevant anaesthesia exposures on the developing primate brain. It will be important to examine if, or to what degree, the maturing brain can recover from these white matter changes.

Impact of gonadectomy on maturational changes in brain volume in adolescent macaques

Adolescence is a transitional period between childhood and adulthood characterized by significant changes in global and regional brain tissue volumes. It is also a period of increasing vulnerability to psychiatric illness. The relationship between these patterns and increased levels of circulating sex steroids during adolescence remains unclear. The objective of the current study was to determine whether gonadectomy, prior to puberty, alters adolescent brain development in male rhesus macaques. Ninety-six structural MRI scans were acquired from 12 male rhesus macaques (8 time points per animal over a two-year period). Six animals underwent gonadectomy and 6 animals underwent a sham operation at 29 months of age. Mixed-effects models were used to determine whether gonadectomy altered developmental trajectories of global and regional brain tissue volumes. We observed a significant effect of gonadectomy on the developmental trajectory of prefrontal gray matter (GM), with intact males showing peak volumes around 3.5 years of age with a subsequent decline. In contrast, prefrontal GM volumes continued to increase in gonadectomized males until the end of the study. We did not observe a significant effect of gonadectomy on prefrontal white matter or on any other global or regional brain tissue volumes, though we cannot rule out that effects might be detected in a larger sample. Results suggest that the prefrontal cortex is more vulnerable to gonadectomy than other brain regions.

Normative Analysis of Individual Brain Differences Based on a Population MRI-Based Atlas of Cynomolgus Macaques

The developmental trajectory of the primate brain varies substantially with aging across subjects. However, this ubiquitous variability between individuals in brain structure is difficult to quantify and has thus essentially been ignored. Based on a large-scale structural magnetic resonance imaging dataset acquired from 162 cynomolgus macaques, we create a species-specific 3D template atlas of the macaque brain, and deploy normative modeling to characterize individual variations of cortical thickness (CT) and regional gray matter volume (GMV). We observed an overall decrease in total GMV and mean CT, and an increase in white matter volume from juvenile to early adult. Specifically, CT and regional GMV were greater in prefrontal and temporal cortices relative to early unimodal areas. Age-dependent trajectories of thickness and volume for each cortical region revealed an increase in the medial temporal lobe, and decreases in all other regions. A low percentage of highly individualized deviations of CT and GMV were identified (0.0021%, 0.0043%, respectively, P < 0.05, false discovery rate [FDR]-corrected). Our approach provides a natural framework to parse individual neuroanatomical differences for use as a reference standard in macaque brain research, potentially enabling inferences regarding the degree to which behavioral or symptomatic variables map onto brain structure in future disease studies.

PREEMACS: Pipeline for preprocessing and extraction of the macaque brain surface

Accurate extraction of the cortical brain surface is critical for cortical thickness estimation and a key element to perform multimodal imaging analysis, where different metrics are integrated and compared in a common space. While brain surface extraction has become widespread practice in human studies, several challenges unique to neuroimaging of non-human primates (NHP) have hindered its adoption for the study of macaques. Although, some of these difficulties can be addressed at the acquisition stage, several common artifacts can be minimized through image preprocessing. Likewise, there are several image analysis pipelines for human MRIs, but very few automated methods for extraction of cortical surfaces have been reported for NHPs and none have been tested on data from diverse sources. We present PREEMACS, a pipeline that standardizes the preprocessing of structural MRI images (T1- and T2-weighted) and carries out an automatic surface extraction of the macaque brain. Building upon and extending pre-existing tools, the first module performs volume orientation, image cropping, intensity non-uniformity correction, and volume averaging, before skull-stripping through a convolutional neural network. The second module performs quality control using an adaptation of MRIqc method to extract objective quality metrics that are then used to determine the likelihood of accurate brain surface estimation. The third and final module estimates the white matter (wm) and pial surfaces from the T1-weighted volume (T1w) using an NHP customized version of FreeSurfer aided by the T2-weighted volumes (T2w). To evaluate the generalizability of PREEMACS, we tested the pipeline using 57 T1w/T2w NHP volumes acquired at 11 different sites from the PRIME-DE public dataset. Results showed an accurate and robust automatic brain surface extraction from images that passed the quality control segment of our pipeline. This work offers a robust, efficient and generalizable pipeline for the automatic standardization of MRI surface analysis on NHP.

DIKA-Nets: Domain-invariant knowledge-guided attention networks for brain skull stripping of early developing macaques

As non-human primates, macaques have a close phylogenetic relationship to human beings and have been proven to be a valuable and widely used animal model in human neuroscience research. Accurate skull stripping (aka. brain extraction) of brain magnetic resonance imaging (MRI) is a crucial prerequisite in neuroimaging analysis of macaques. Most of the current skull stripping methods can achieve satisfactory results for human brains, but when applied to macaque brains, especially during early brain development, the results are often unsatisfactory. In fact, the early dynamic, regionally-heterogeneous development of macaque brains, accompanied by poor and age-related contrast between different anatomical structures, poses significant challenges for accurate skull stripping. To overcome these challenges, we propose a fully-automated framework to effectively fuse the age-specific intensity information and domain-invariant prior knowledge as important guiding information for robust skull stripping of developing macaques from 0 to 36 months of age. Specifically, we generate Signed Distance Map (SDM) and Center of Gravity Distance Map (CGDM) based on the intermediate segmentation results as guidance. Instead of using local convolution, we fuse all information using the Dual Self-Attention Module (DSAM), which can capture global spatial and channel-dependent information of feature maps. To extensively evaluate the performance, we adopt two relatively-large challenging MRI datasets from rhesus macaques and cynomolgus macaques, respectively, with a total of 361 scans from two different scanners with different imaging protocols. We perform cross-validation by using one dataset for training and the other one for testing. Our method outperforms five popular brain extraction tools and three deep-learning-based methods on cross-source MRI datasets without any transfer learning.

The macaque brain ONPRC18 template with combined gray and white matter labelmap for multimodal neuroimaging studies of Nonhuman Primates

Macaques are the most common nonhuman primate (NHP) species used in neuroscience research. With the advancement of many neuroimaging techniques, new studies are beginning to apply multiple types of in vivo mag netic resonance imaging (MRI), such as structural imaging (sMRI) with T1 and T2 weighted contrasts alongside diffusion weighed (DW) imaging. In studies involving rhesus macaques, this approach can be used to better under stand micro-structural changes that occur during development, in various disease states or with normative aging. However, many of the available rhesus brain atlases have been designed for only one imaging modality, making it difficult to consistently define the same brain regions across multiple imaging modalities in the same subject. To address this, we created a brain atlas from 18 adult rhesus macaques that includes co-registered templates constructed from images frequently used to characterize macroscopic brain structure (T2/SPACE and T1/MP-RAGE), and a diffusion tensor imaging (DTI) template. The DTI template was up-sampled from 1 mm isotropic resolution to resolution match to the T1 and T2-weighted images (0.5 mm isotropic), and the parameter maps were derived for FA, AD, RD and MD. The labelmap volumes delineate 57 gray matter regions of interest (ROIs; 36 cortical regions and 21 subcortical structures), as well as 74 white matter tracts. Importantly, the labelmap overlays both the structural and diffusion templates, enabling the same regions to be consistently identified across imaging modalities. A specialized condensed version of the labelmap ROIs are also included to further extend the usefulness of this tool for imaging data with lower spatial resolution, such as functional MRI (fMRI) or positron emission tomography (PET).

Diffusion Tensor Imaging-Based Studies at the Group-Level Applied to Animal Models of Neurodegenerative Diseases

The understanding of human and non-human microstructural brain alterations in the course of neurodegenerative diseases has substantially improved by the non-invasive magnetic resonance imaging (MRI) technique of diffusion tensor imaging (DTI). Animal models (including disease or knockout models) allow for a variety of experimental manipulations, which are not applicable to humans. Thus, the DTI approach provides a promising tool for cross-species cross-sectional and longitudinal investigations of the neurobiological targets and mechanisms of neurodegeneration. This overview with a systematic review focuses on the principles of DTI analysis as used in studies at the group level in living preclinical models of neurodegeneration. The translational aspect from in-vivo animal models toward (clinical) applications in humans is covered as well as the DTI-based research of the non-human brains' microstructure, the methodological aspects in data processing and analysis, and data interpretation at different abstraction levels. The aim of integrating DTI in multiparametric or multimodal imaging protocols will allow the interrogation of DTI data in terms of directional flow of information and may identify the microstructural underpinnings of neurodegeneration-related patterns.

Rhesus Macaque Brain Developmental Trajectory: A Longitudinal Analysis Using Tensor-Based Structural Morphometry and Diffusion Tensor Imaging

The typical developmental trajectory of brain structure among nonhuman primates (NHPs) remains poorly understood. In this study, we characterized the normative trajectory of developmental change among a cohort of rhesus monkeys (n = 28), ranging in age from 2 to 22 months, using structural MRI datasets that were longitudinally acquired every 3-4 months. We hypothesized that NHP-specific transient intracranial volume decreases reported during late infancy would be part of the typical developmental process, which is driven by volumetric contraction of gray matter in primary functional areas. To this end, we performed multiscale analyses from the whole brain to voxel level, characterizing regional heterogeneity, hemispheric asymmetry, and sexual dimorphism in developmental patterns. The longitudinal trajectory of brain development was explained by three different regional volumetric growth patterns (exponentially decreasing, undulating, and linearly increasing), which resulted in developmental brain volume curves with transient brain volumetric decreases. White matter (WM) fractional anisotropy increased with age, corresponding to WM volume increases, while mean diffusivity (MD) showed biphasic patterns. The longitudinal trajectory of brain development in young rhesus monkeys follows typical maturation patterns seen in humans, but regional volumetric and MD changes are more dynamic in rhesus monkeys compared with humans, with marked decreases followed by "rebound-like" increases.

High Angular Resolution Diffusion MRI Reveals Conserved and Deviant Programs in the Paths that Guide Human Cortical Circuitry

Diffusion magnetic resonance (MR) tractography represents a novel opportunity to investigate conserved and deviant developmental programs between humans and other species such as mice. To that end, we acquired high angular resolution diffusion MR scans of mice [embryonic day (E) 10.5 to postnatal week 4] and human brains [gestational week (GW) 17-30] at successive stages of fetal development to investigate potential evolutionary changes in radial organization and emerging pathways between humans and mice. We compare radial glial development as well as commissural development (e.g., corpus callosum), primarily because our findings can be integrated with previous work. We also compare corpus callosal growth trajectories across primates (i.e., humans and rhesus macaques) and rodents (i.e., mice). One major finding is that the developing cortex of humans is predominated by pathways likely associated with a radial glial organization at GW 17-20, which is not as evident in age-matched mice (E 16.5, 17.5). Another finding is that, early in development, the corpus callosum follows a similar developmental timetable in primates (i.e., macaques and humans) as in mice. However, the corpus callosum grows for an extended period of time in primates compared with rodents. Taken together, these findings highlight deviant developmental programs underlying the emergence of cortical pathways in the human brain.

Mapping hemispheric asymmetries of the macaque cerebral cortex during early brain development

Studying cortical hemispheric asymmetries during the dynamic early postnatal stages in macaque monkeys (with close phylogenetic relationship to humans) would increase our limited understanding on the possible origins, developmental trajectories, and evolutional mechanisms of brain asymmetries in nonhuman primates, but remains a blind spot to the community. Via cortical surface-based morphometry, we comprehensively analyze hemispheric structural asymmetries in 134 longitudinal MRI scans from birth to 20 months of age from 32 healthy macaque monkeys. We reveal that most clusters of hemispheric asymmetries of cortical properties, such as surface area, cortical thickness, sulcal depth, and vertex positions, expand in the first 4 months of life, and evolve only moderately thereafter. Prominent hemispheric asymmetries are found at the inferior frontal gyrus, precentral gyrus, posterior temporal cortex, superior temporal gyrus (STG), superior temporal sulcus (STS), and cingulate cortex. Specifically, the left planum temporale and left STG consistently have larger area and thicker cortices than those on the right hemisphere, while the right STS, right cingulate cortex, and right anterior insula are consistently deeper than the left ones, partially consistent with the findings in human infants and adults. Our results thus provide a valuable reference in studying early brain development and evolution.

Large-scale comparative neuroimaging: Where are we and what do we need?

Neuroimaging has a lot to offer comparative neuroscience. Although invasive "gold standard" techniques have a better spatial resolution, neuroimaging allows fast, whole-brain, repeatable, and multi-modal measurements of structure and function in living animals and post-mortem tissue. In the past years, comparative neuroimaging has increased in popularity. However, we argue that its most significant potential lies in its ability to collect large-scale datasets of many species to investigate principles of variability in brain organisation across whole orders of species-an ambition that is presently unfulfilled but achievable. We briefly review the current state of the field and explore what the current obstacles to such an approach are. We propose some calls to action.

atlasBREX: Automated template-derived brain extraction in animal MRI

We proposed a generic template-derived approach for (semi-) automated brain extraction in animal MRI studies and evaluated our implementation with different animal models (macaque, marmoset, rodent) and MRI protocols (T1, T2). While conventional MR-neuroimaging studies perform brain extraction as an initial step priming subsequent image-registration from subject to template, our proposed approach propagates an anatomical template to (whole-head) individual subjects in reverse order, which is challenging due to the surrounding extracranial tissue, greater differences in contrast pattern and larger areas with field inhomogeneity. As a novel approach, the herein introduced brain extraction algorithm derives whole-brain segmentation using rigid and non-rigid deformation based on unbiased anatomical atlas building with a priori estimates from study-cohort and an initial approximate brain extraction. We evaluated our proposed method in comparison to several other technical approaches including “Marker based watershed scalper”, “Brain-Extraction-Tool”, “3dSkullStrip”, “Primatologist-Toolbox”, “Rapid Automatic Tissue Segmentation” and “Robust automatic rodent brain extraction using 3D pulse-coupled neural networks” with manual skull-stripping as reference standard. ABX demonstrated best performance with accurate (≥92%) and consistent results throughout datasets and across species, age and MRI protocols. ABX was made available to the public with documentation, templates and sample material (https://www.github.com/jlohmeier/atlasBREX).

CHARTING DEVELOPMENT-BASED JOINT PARCELLATION MAPS OF HUMAN AND MACAQUE BRAINS DURING INFANCY

Comparative characterization of early brain development between human and macaque using neuroimaging data is crucial to understand the mechanisms of brain development and evolution. To this end, joint cortical parcellation maps of human and macaque infant brains with corresponding regions are highly desirable, since they provide basic cortical parcels for both region-based and network-based studies of two closely-related species. To address this issue, we propose to leverage developmental patterns of cortical properties of both human and macaque infants for creating joint parcellation maps with inter-species comparability. The motivation is that the developmental patterns of cortical properties indicate underlying rapid changes of microstructures, which determine the molecular and functional principles of the cortex. Thus, developmental patterns are well suitable for defining distinct cortical regions in both structures and functions. To comprehensively capture the similarities of developmental patterns of vertices on cortical surfaces, for each species, we first construct two complementary similarity matrices: a low-order matrix and a high-order matrix. Then, we non-linearly fuse these four matrices together as a single matrix in a hierarchical manner, thus capturing the common and complementary information of both human and macaque infants. Finally, based on the fused similarity matrix, we apply the spectral clustering to derive the joint parcellation maps. By applying our method to 210 longitudinal human infant MRI scans and 140 longitudinal macaque infant MRI scans, we generate the first biologically-meaningful joint parcellation maps of human and macaque infants.

Individual variation in longitudinal postnatal development of the primate brain

Quantifying individual variation in postnatal brain development can provide insight into cognitive diversity within a population and the aetiology of common neuropsychiatric and neurodevelopmental disorders. Non-invasive studies of the non-human primate can aid understanding of human brain development, facilitating longitudinal analysis during early postnatal development when comparative human populations are difficult to sample. In this study, we perform analysis of a longitudinal MRI dataset of 32 macaques, each with up to five magnetic resonance imaging (MRI) scans acquired between 3 and 36 months of age. Using nonlinear mixed effects model we derive growth trajectories for whole brain, cortical and subcortical grey matter, cerebral white matter and cerebellar volume. We then test the association between individual variation in postnatal tissue volumes and birth weight. We report nonlinear growth models for all tissue compartments, as well as significant variation in total intracranial volume between individuals. We also demonstrate that regional subcortical grey matter varies both in total volume and rate of change between individuals and is associated with differences in birth weight. This supports evidence that birth weight may act as a marker of subsequent brain development and highlights the importance of longitudinal MRI analysis in developmental studies.

Optimization of macaque brain DMRI connectome by neuron tracing and myelin stain data

Accurate assessment of connectional anatomy of primate brains can be an important avenue to better understand the structural and functional organization of brains. To this end, numerous connectome projects have been initiated to create a comprehensive map of the connectional anatomy over a large spatial expanse. Tractography based on diffusion MRI (dMRI) data has been used as a tool by many connectome projects in that it is widely used to visualize axonal pathways and reveal microstructural features on living brains. However, the measures obtained from dMRI are indirect inference of microstructures. This intrinsic limitation reduces the reliability of dMRI in constructing connectomes for brains. In this work, we proposed a framework to increase the accuracy of constructing a dMRI-based connectome on macaque brains by integrating meso-scale connective information from tract-tracing data and micro-scale axonal orientation information from myelin stain data. Our results suggest that this integrative framework could advance the mapping accuracy of dMRI based connections and axonal pathways, and demonstrate the prospect of the proposed framework in constructing a large-scale connectome on living primate brains.

An Open Resource for Non-human Primate Imaging

Non-human primate neuroimaging is a rapidly growing area of research that promises to transform and scale translational and cross-species comparative neuroscience. Unfortunately, the technological and methodological advances of the past two decades have outpaced the accrual of data, which is particularly challenging given the relatively few centers that have the necessary facilities and capabilities. The PRIMatE Data Exchange (PRIME-DE) addresses this challenge by aggregating independently acquired non-human primate magnetic resonance imaging (MRI) datasets and openly sharing them via the International Neuroimaging Data-sharing Initiative (INDI). Here, we present the rationale, design, and procedures for the PRIME-DE consortium, as well as the initial release, consisting of 25 independent data collections aggregated across 22 sites (total = 217 non-human primates). We also outline the unique pitfalls and challenges that should be considered in the analysis of non-human primate MRI datasets, including providing automated quality assessment of the contributed datasets.

CONSTRUCTION OF SPATIOTEMPORAL INFANT CORTICAL SURFACE ATLAS OF RHESUS MACAQUE

As a widely used animal model in MR imaging studies, rhesus macaque helps to better understand both normal and abnormal neural development in the human brain. However, the available adult macaque brain atlases are not well suitable for study of brain development at the early postnatal stage, since this stage undergoes dramatic changes in brain appearances and structures. Building age matched atlases for this critical period is thus highly desirable yet still lacking. In this paper, we construct the first spatiotemporal (4D) cortical surface atlases for rhesus macaques from 2 weeks to 24 months, using 138 longitudinal MRI scans from 32 healthy rhesus monkeys. Specifically, we first perform intra-subject cortical surface registration to obtain within-subject mean cortical surfaces. Then, we perform inter-subject registration of within-subject mean surfaces to obtain unbiased and longitudinally-consistent 4D cortical surface atlases. Based on our 4D rhesus monkey atlases, we further chart the first developmental-trajectories-based parcellation maps using the local surface area and spectral clustering algorithm. Our 4D macaque surface atlases and parcellation maps will greatly facilitate early brain development studies of macaques.

Morphometric Similarity Networks Detect Microscale Cortical Organization and Predict Inter-Individual Cognitive Variation

Macroscopic cortical networks are important for cognitive function, but it remains challenging to construct anatomically plausible individual structural connectomes from human neuroimaging. We introduce a new technique for cortical network mapping based on inter-regional similarity of multiple morphometric parameters measured using multimodal MRI. In three cohorts (two human, one macaque), we find that the resulting morphometric similarity networks (MSNs) have a complex topological organization comprising modules and high-degree hubs. Human MSN modules recapitulate known cortical cytoarchitectonic divisions, and greater inter-regional morphometric similarity was associated with stronger inter-regional co-expression of genes enriched for neuronal terms. Comparing macaque MSNs with tract-tracing data confirmed that morphometric similarity was related to axonal connectivity. Finally, variation in the degree of human MSN nodes accounted for about 40% of between-subject variability in IQ. Morphometric similarity mapping provides a novel, robust, and biologically plausible approach to understanding how human cortical networks underpin individual differences in psychological functions.

Primatologist: A modular segmentation pipeline for macaque brain morphometry

Because they bridge the genetic gap between rodents and humans, non-human primates (NHPs) play a major role in therapy development and evaluation for neurological disorders. However, translational research success from NHPs to patients requires an accurate phenotyping of the models. In patients, magnetic resonance imaging (MRI) combined with automated segmentation methods has offered the unique opportunity to assess in vivo brain morphological changes. Meanwhile, specific challenges caused by brain size and high field contrasts make existing algorithms hard to use routinely in NHPs. To tackle this issue, we propose a complete pipeline, Primatologist, for multi-region segmentation. Tissue segmentation is based on a modular statistical model that includes random field regularization, bias correction and denoising and is optimized by expectation-maximization. To deal with the broad variety of structures with different relaxing times at 7 T, images are segmented into 17 anatomical classes, including subcortical regions. Pre-processing steps insure a good initialization of the parameters and thus the robustness of the pipeline. It is validated on 10 T2-weighted MRIs of healthy macaque brains. Classification scores are compared with those of a non-linear atlas registration, and the impact of each module on classification scores is thoroughly evaluated.

Advances in nonhuman primate models of autism: Integrating neuroscience and behavior

Given the prevalence and societal impact of autism spectrum disorders (ASD), there is an urgent need to develop innovative preventative strategies and treatments to reduce the alarming number of cases and improve core symptoms for afflicted individuals. Translational efforts between clinical and preclinical research are needed to (i) identify and evaluate putative causes of ASD, (ii) determine the underlying neurobiological mechanisms, (iii) develop and test novel therapeutic approaches and (iv) ultimately translate basic research into safe and effective clinical practices. However, modeling a uniquely human brain disorder, such as ASD, will require sophisticated animal models that capitalize on unique advantages of diverse species including drosophila, zebra fish, mice, rats, and ultimately, species more closely related to humans, such as the nonhuman primate. Here we discuss the unique contributions of the rhesus monkey (Macaca mulatta) model to ongoing efforts to understand the neurobiology of the disorder, focusing on the convergence of brain and behavior outcome measures that parallel features of human ASD.

A population MRI brain template and analysis tools for the macaque

The use of standard anatomical templates is common in human neuroimaging, as it facilitates data analysis and comparison across subjects and studies. For non-human primates, previous in vivo templates have lacked sufficient contrast to reliably validate known anatomical brain regions and have not provided tools for automated single-subject processing. Here we present the "National Institute of Mental Health Macaque Template", or NMT for short. The NMT is a high-resolution in vivo MRI template of the average macaque brain generated from 31 subjects, as well as a neuroimaging tool for improved data analysis and visualization. From the NMT volume, we generated maps of tissue segmentation and cortical thickness. Surface reconstructions and transformations to previously published digital brain atlases are also provided. We further provide an analysis pipeline using the NMT that automates and standardizes the time-consuming processes of brain extraction, tissue segmentation, and morphometric feature estimation for anatomical scans of individual subjects. The NMT and associated tools thus provide a common platform for precise single-subject data analysis and for characterizations of neuroimaging results across subjects and studies.