Cingulate cortex shapes early postnatal development of social vocalizations.

The social dynamics of vocal behavior has major implications for social development in humans. We asked whether early life damage to the anterior cingulate cortex (ACC), which is closely associated with socioemotional regulation more broadly, impacts the normal development of vocal expression. The common marmoset provides a unique opportunity to study the developmental trajectory of vocal behavior, and to track the consequences of early brain damage on aspects of social vocalizations. We created ACC lesions in neonatal marmosets and compared their pattern of vocalization to that of age-matched controls throughout the first 6 weeks of life. We found that while early life ACC lesions had little influence on the production of vocal calls, developmental changes to the quality of social contact calls and their associated syntactical and acoustic characteristics were compromised. These animals made fewer social contact calls, and when they did, they were short, loud and monotonic. We further determined that damage to ACC in infancy results in a permanent alteration in downstream brain areas known to be involved in social vocalizations, such as the amygdala and periaqueductal gray. Namely, in the adult, these structures exhibited diminished GABA-immunoreactivity relative to control animals, likely reflecting disruption of the normal inhibitory balance following ACC deafferentation. Together, these data indicate that the normal development of social vocal behavior depends on the ACC and its interaction with other areas in the vocal network during early life.

Multimodal anatomical mapping of subcortical regions in marmoset monkeys using high-resolution MRI and matched histology with multiple stains.

Subcortical nuclei and other deep brain structures play essential roles in regulating the central and peripheral nervous systems. However, many of these nuclei and their subregions are challenging to identify and delineate in conventional MRI due to their small size, hidden location, and often subtle contrasts compared to neighboring regions. To address these limitations, we scanned the whole brain of the marmoset monkeys in ex vivo using a clinically feasible diffusion MRI method, called the mean apparent propagator (MAP)-MRI, along with T2W and MTR (T1-like contrast) images acquired at 7 Tesla. Additionally, we registered these multimodal MRI volumes to the high-resolution images of matched whole-brain histology sections with seven different stains obtained from the same brain specimens. At high spatial resolution, the microstructural parameters and fiber orientation distribution functions derived with MAP-MRI can distinguish the subregions of many subcortical and deep brain structures, including fiber tracts of different sizes and orientations. The good correlation with multiple but distinct histological stains from the same brain serves as a thorough validation of the structures identified with MAP-MRI and other MRI parameters. Moreover, the anatomical details of deep brain structures found in the volumes of MAP-MRI parameters are not visible in conventional T1W or T2W images. The high-resolution mapping using novel MRI contrasts, combined and correlated with histology, can elucidate structures that were previously invisible radiologically. Thus, this multimodal approach offers a roadmap toward identifying salient brain areas in vivo in future neuroradiological studies. It also provides a useful anatomical standard reference for the region definition of subcortical targets and the generation of a 3D digital template atlas for the marmoset brain research (Saleem et al., 2023). Additionally, we conducted a cross-species comparison between marmoset and macaque monkeys using results from our previous studies (Saleem et al., 2021). We found that the two species had distinct patterns of iron distribution in subregions of the basal ganglia, red nucleus, and deep cerebellar nuclei, confirmed with T2W MRI and histology.

Multimodal anatomical mapping of subcortical regions in Marmoset monkeys using high-resolution MRI and matched histology with multiple stains.

Subcortical nuclei and other deep brain structures play essential roles in regulating the central and peripheral nervous systems. However, many of these nuclei and their subregions are challenging to identify and delineate in conventional MRI due to their small size, hidden location, and often subtle contrasts compared to neighboring regions. To address these limitations, we scanned the whole brain of the marmoset monkeys in ex vivo using a clinically feasible diffusion MRI method, called the mean apparent propagator (MAP)-MRI, along with T2W and MTR (T1-like contrast) images acquired at 7 Tesla. Additionally, we registered these multimodal MRI volumes to the high-resolution images of matched whole-brain histology sections with seven different stains obtained from the same brain specimens. At high spatial resolution, the microstructural parameters and fiber orientation distribution functions derived with MAP-MRI can distinguish the subregions of many subcortical and deep brain structures, including fiber tracts of different sizes and orientations. The good correlation with multiple but distinct histological stains from the same brain serves as a thorough validation of the structures identified with MAP-MRI and other MRI parameters. Moreover, the anatomical details of deep brain structures found in the volumes of MAP-MRI parameters are not visible in conventional T1W or T2W images. The high-resolution mapping using novel MRI contrasts, combined and correlated with histology, can elucidate structures that were previously invisible radiologically. Thus, this multimodal approach offers a roadmap toward identifying salient brain areas in vivo in future neuroradiological studies. It also provides a useful anatomical standard reference for the region definition of subcortical targets and the generation of a 3D digital template atlas for the marmoset brain research (Saleem et al., 2023). Additionally, we conducted a cross-species comparison between marmoset and macaque monkeys using results from our previous studies (Saleem et al., 2021). We found that the two species had distinct patterns of iron distribution in subregions of the basal ganglia, red nucleus, and deep cerebellar nuclei, confirmed with T2W MRI and histology.

In vivo MRI is sensitive to remyelination in a nonhuman primate model of multiple sclerosis.

Remyelination is crucial to recover from inflammatory demyelination in multiple sclerosis (MS). Investigating remyelination in vivo using magnetic resonance imaging (MRI) is difficult in MS, where collecting serial short-interval scans is challenging. Using experimental autoimmune encephalomyelitis (EAE) in common marmosets, a model of MS that recapitulates focal cerebral inflammatory demyelinating lesions, we investigated whether MRI is sensitive to, and can characterize, remyelination. In six animals followed with multisequence 7 T MRI, 31 focal lesions, predicted to be demyelinated or remyelinated based on signal intensity on proton density-weighted images, were subsequently assessed with histopathology. Remyelination occurred in four of six marmosets and 45% of lesions. Radiological-pathological comparison showed that MRI had high statistical sensitivity (100%) and specificity (90%) for detecting remyelination. This study demonstrates the prevalence of spontaneous remyelination in marmoset EAE and the ability of in vivo MRI to detect it, with implications for preclinical testing of pro-remyelinating agents.

An integrated resource for functional and structural connectivity of the marmoset brain.

Comprehensive integration of structural and functional connectivity data is required to model brain functions accurately. While resources for studying the structural connectivity of non-human primate brains already exist, their integration with functional connectivity data has remained unavailable. Here we present a comprehensive resource that integrates the most extensive awake marmoset resting-state fMRI data available to date (39 marmoset monkeys, 710 runs, 12117 mins) with previously published cellular-level neuronal tracing data (52 marmoset monkeys, 143 injections) and multi-resolution diffusion MRI datasets. The combination of these data allowed us to (1) map the fine-detailed functional brain networks and cortical parcellations, (2) develop a deep-learning-based parcellation generator that preserves the topographical organization of functional connectivity and reflects individual variabilities, and (3) investigate the structural basis underlying functional connectivity by computational modeling. This resource will enable modeling structure-function relationships and facilitate future comparative and translational studies of primate brains.

High-resolution cortical MAP-MRI reveals areal borders and laminar substructures observed with histological staining.

The variations in cellular composition and tissue architecture measured with histology provide the biological basis for partitioning the brain into distinct cytoarchitectonic areas and for characterizing neuropathological tissue alterations. Clearly, there is an urgent need to develop whole-brain neuroradiological methods that can assess cortical cyto- and myeloarchitectonic features non-invasively. Mean apparent propagator (MAP) MRI is a clinically feasible diffusion MRI method that quantifies efficiently and comprehensively the net microscopic displacements of water molecules diffusing in tissues. We investigate the sensitivity of high-resolution MAP-MRI to detecting areal and laminar variations in cortical cytoarchitecture and compare our results with observations from corresponding histological sections in the entire brain of a rhesus macaque monkey. High-resolution images of MAP-derived parameters, in particular the propagator anisotropy (PA), non-gaussianity (NG), and the return-to-axis probability (RTAP) reveal cortical area-specific lamination patterns in good agreement with the corresponding histological stained sections. In a few regions, the MAP parameters provide superior contrast to the five histological stains used in this study, delineating more clearly boundaries and transition regions between cortical areas and laminar substructures. Throughout the cortex, various MAP parameters can be used to delineate transition regions between specific cortical areas observed with histology and to refine areal boundaries estimated using atlas registration-based cortical parcellation. Using surface-based analysis of MAP parameters we quantify the cortical depth dependence of diffusion propagators in multiple regions-of-interest in a consistent and rigorous manner that is largely independent of the cortical folding geometry. The ability to assess cortical cytoarchitectonic features efficiently and non-invasively, its clinical feasibility, and translatability make high-resolution MAP-MRI a promising 3D imaging tool for studying whole-brain cortical organization, characterizing abnormal cortical development, improving early diagnosis of neurodegenerative diseases, identifying targets for biopsies, and complementing neuropathological investigations.

An open access resource for functional brain connectivity from fully awake marmosets.

The common marmoset (Callithrix jacchus) is quickly gaining traction as a premier neuroscientific model. However, considerable progress is still needed in understanding the functional and structural organization of the marmoset brain to rival that documented in longstanding preclinical model species, like mice, rats, and Old World primates. To accelerate such progress, we present the Marmoset Functional Brain Connectivity Resource (marmosetbrainconnectome.org), currently consisting of over 70 h of resting-state fMRI (RS-fMRI) data acquired at 500 µm isotropic resolution from 31 fully awake marmosets in a common stereotactic space. Three-dimensional functional connectivity (FC) maps for every cortical and subcortical gray matter voxel are stored online. Users can instantaneously view, manipulate, and download any whole-brain functional connectivity (FC) topology (at the subject- or group-level) along with the raw datasets and preprocessing code. Importantly, researchers can use this resource to test hypotheses about FC directly - with no additional analyses required - yielding whole-brain correlations for any gray matter voxel on demand. We demonstrate the resource's utility for presurgical planning and comparison with tracer-based neuronal connectivity as proof of concept. Complementing existing structural connectivity resources for the marmoset brain, the Marmoset Functional Brain Connectivity Resource affords users the distinct advantage of exploring the connectivity of any voxel in the marmoset brain, not limited to injection sites nor constrained by regional atlases. With the entire raw database (RS-fMRI and structural images) and preprocessing code openly available for download and use, we expect this resource to be broadly valuable to test novel hypotheses about the functional organization of the marmoset brain.

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.

Ultrahigh-resolution MRI Reveals Extensive Cortical Demyelination in a Nonhuman Primate Model of Multiple Sclerosis.

Cortical lesions are a primary driver of disability in multiple sclerosis (MS). However, noninvasive detection of cortical lesions with in vivo magnetic resonance imaging (MRI) remains challenging. Experimental autoimmune encephalomyelitis (EAE) in the common marmoset is a relevant animal model of MS for investigating the pathophysiological mechanisms leading to brain damage. This study aimed to characterize cortical lesions in marmosets with EAE using ultrahigh-field (7 T) MRI and histological analysis. Tissue preparation was optimized to enable the acquisition of high-spatial resolution (50-µm isotropic) T2*-weighted images. A total of 14 animals were scanned in this study, and 70% of the diseased animals presented at least one cortical lesion on postmortem imaging. Cortical lesions identified on MRI were verified with myelin proteolipid protein immunostaining. An optimized T2*-weighted sequence was developed for in vivo imaging and shown to capture 65% of cortical lesions detected postmortem. Immunostaining confirmed extensive demyelination with preserved neuronal somata in several cortical areas of EAE animals. Overall, this study demonstrates the relevance and feasibility of the marmoset EAE model to study cortical lesions, among the most important yet least understood features of MS.

Marmoset Brain Mapping V3: Population multi-modal standard volumetric and surface-based templates.

The standard anatomical brain template provides a common space and coordinate system for visualizing and analyzing neuroimaging data from large cohorts of subjects. Previous templates and atlases for the common marmoset brain were either based on data from a single individual or lacked essential functionalities for neuroimaging analysis. Here, we present new population-based in-vivo standard templates and tools derived from multi-modal data of 27 marmosets, including multiple types of T1w and T2w contrast images, DTI contrasts, and large field-of-view MRI and CT images. We performed multi-atlas labeling of anatomical structures on the new templates and constructed highly accurate tissue-type segmentation maps to facilitate volumetric studies. We built fully featured brain surfaces and cortical flat maps to facilitate 3D visualization and surface-based analyses, which are compatible with most surface analyzing tools, including FreeSurfer, AFNI/SUMA, and the Connectome Workbench. Analysis of the MRI and CT datasets revealed significant variations in brain shapes, sizes, and regional volumes of brain structures, highlighting substantial individual variabilities in the marmoset population. Thus, our population-based template and associated tools provide a versatile analysis platform and standard coordinate system for a wide range of MRI and connectome studies of common marmosets. These new template tools comprise version 3 of our Marmoset Brain Mapping Project and are publicly available via marmosetbrainmapping.org/v3.html.

Long-distance aberrant heterotopic connectivity in a mouse strain with a high incidence of callosal anomalies.

Corpus callosum dysgenesis (CCD) is a developmental brain condition in which some white matter fibers fail to find their natural course across the midplane, reorganizing instead to form new aberrant pathways. This type of white matter reorganization is known as long-distance plasticity (LDP). The present work aimed to characterize the Balb/c mouse strain as a model of CCD. We employed high-resolution anatomical MRI in 81 Balb/c and 27 C57bl6 mice to show that the Balb/c mouse strain presents a variance in the size of the CC that is 3.9 times higher than the variance of normotypical C57bl6. We also performed high-resolution diffusion-weighted imaging (DWI) in 8 Balb/c and found that the Balb/c strain shows aberrant white matter bundles, such as the Probst (5/8 animals) and the Sigmoid bundles (7/8 animals), which are similar to those found in humans with CCD. Using a histological tracer technique, we confirmed the existence of these aberrant bundles in the Balb/c strain. Interestingly, we also identified sigmoid-like fibers in the C57bl6 strain, thought to a lesser degree. Next, we used a connectome approach and found widespread brain connectivity differences between Balb/c and C57bl6 strains. The Balb/c strain also exhibited increased variability of global connectivity. These findings suggest that the Balb/c strain presents local and global changes in brain structural connectivity. This strain often presents with callosal abnormalities, along with the Probst and the Sigmoid bundles, making it is an attractive animal model for CCD and LDP in general. Our results also show that even the C57bl6 strain, which typically serves as a normotypical control animal in a myriad of studies, presents sigmoid-fashion pattern fibers laid out in the brain. These results suggest that these aberrant fiber pathways may not necessarily be a pathological hallmark, but instead an alternative roadmap for misguided axons. Such findings offer new insights for interpreting the significance of CCD-associated LDP in humans.

A resource for the detailed 3D mapping of white matter pathways in the marmoset brain.

While the fundamental importance of the white matter in supporting neuronal communication is well known, existing publications of primate brains do not feature a detailed description of its complex anatomy. The main barrier to achieving this is that existing primate neuroimaging data have insufficient spatial resolution to resolve white matter pathways fully. Here we present a resource that allows detailed descriptions of white matter structures and trajectories of fiber pathways in the marmoset brain. The resource includes: (1) the highest-resolution diffusion-weighted MRI data available to date, which reveal white matter features not previously described; (2) a comprehensive three-dimensional white matter atlas depicting fiber pathways that were either omitted or misidentified in previous atlases; and (3) comprehensive fiber pathway maps of cortical connections combining diffusion-weighted MRI tractography and neuronal tracing data. The resource, which can be downloaded from marmosetbrainmapping.org, will facilitate studies of brain connectivity and the development of tractography algorithms in the primate brain.

Dynamic Interhemispheric Desynchronization in Marmosets and Humans With Disorders of the Corpus Callosum.

The corpus callosum, the principal structural avenue for interhemispheric neuronal communication, controls the brain's lateralization. Developmental malformations of the corpus callosum (CCD) can lead to learning and intellectual disabilities. Currently, there is no clear explanation for these symptoms. Here, we used resting-state functional MRI (rsfMRI) to evaluate the dynamic resting-state functional connectivity (rsFC) in both the cingulate cortex (CG) and the sensory areas (S1, S2, A1) in three marmosets (Callithrix jacchus) with spontaneous CCD. We also performed rsfMRI in 10 CCD human subjects (six hypoplasic and four agenesic). We observed no differences in the strength of rsFC between homotopic CG and sensory areas in both species when comparing them to healthy controls. However, in CCD marmosets, we found lower strength of quasi-periodic patterns (QPP) correlation in the posterior interhemispheric sensory areas. We also found a significant lag of interhemispheric communication in the medial CG, suggesting asynchrony between the two hemispheres. Correspondingly, in human subjects, we found that the CG of acallosal subjects had a higher QPP correlation than controls. In comparison, hypoplasic subjects had a lower QPP correlation and a delay of 1.6 s in the sensory regions. These results show that CCD affects the interhemispheric synchrony of both CG and sensory areas and that, in both species, its impact on cortical communication varies along the CC development gradient. Our study shines a light on how CCD misconnects homotopic regions and opens a line of research to explain the causes of the symptoms exhibited by CCD patients and how to mitigate them.

Potential role of iron in repair of inflammatory demyelinating lesions.

Inflammatory destruction of iron-rich myelin is characteristic of multiple sclerosis (MS). Although iron is needed for oligodendrocytes to produce myelin during development, its deposition has also been linked to neurodegeneration and inflammation, including in MS. We report perivascular iron deposition in multiple sclerosis lesions that was mirrored in 72 lesions from 13 marmosets with experimental autoimmune encephalomyelitis. Iron accumulated mainly inside microglia/macrophages from 6 weeks after demyelination. Consistently, expression of transferrin receptor, the brain's main iron-influx protein, increased as lesions aged. Iron was uncorrelated with inflammation and postdated initial demyelination, suggesting that iron is not directly pathogenic. Iron homeostasis was at least partially restored in remyelinated, but not persistently demyelinated, lesions. Taken together, our results suggest that iron accumulation in the weeks after inflammatory demyelination may contribute to lesion repair rather than inflammatory demyelination per se.

Anatomical and functional investigation of the marmoset default mode network.

The default mode network (DMN) is associated with a wide range of brain functions. In humans, the DMN is marked by strong functional connectivity among three core regions: medial prefrontal cortex (mPFC), posterior parietal cortex (PPC), and the medial parietal and posterior cingulate cortex (PCC). Neuroimaging studies have shown that the DMN also exists in non-human primates, suggesting that it may be a conserved feature of the primate brain. Here, we found that, in common marmosets, the dorsolateral prefrontal cortex (dlPFC; peak at A8aD) has robust fMRI functional connectivity and reciprocal anatomical connections with the posterior DMN core regions (PPC and PCC), while the mPFC has weak connections with the posterior DMN core regions. This strong dlPFC but weak mPFC connectivity in marmoset differs markedly from the stereotypical DMN in humans. The mPFC may be involved in brain functions that are further developed in humans than in other primates.

Herpesvirus trigger accelerates neuroinflammation in a nonhuman primate model of multiple sclerosis.

Pathogens, particularly human herpesviruses (HHVs), are implicated as triggers of disease onset/progression in multiple sclerosis (MS) and other neuroinflammatory disorders. However, the time between viral acquisition in childhood and disease onset in adulthood complicates the study of this association. Using nonhuman primates, we demonstrate that intranasal inoculations with HHV-6A and HHV-6B accelerate an MS-like neuroinflammatory disease, experimental autoimmune encephalomyelitis (EAE). Although animals inoculated intranasally with HHV-6 (virus/EAE marmosets) were asymptomatic, they exhibited significantly accelerated clinical EAE compared with control animals. Expansion of a proinflammatory CD8 subset correlated with post-EAE survival in virus/EAE marmosets, suggesting that a peripheral (viral?) antigen-driven expansion may have occurred post-EAE induction. HHV-6 viral antigen in virus/EAE marmosets was markedly elevated and concentrated in brain lesions, similar to previously reported localizations of HHV-6 in MS brain lesions. Collectively, we demonstrate that asymptomatic intranasal viral acquisition accelerates subsequent neuroinflammation in a nonhuman primate model of MS.

Neuroprotective Effects of MAGL (Monoacylglycerol Lipase) Inhibitors in Experimental Ischemic Stroke.

BACKGROUND AND PURPOSE: MAGL (monoacylglycerol lipase) is an enzyme that hydrolyzes the endocannabinoid 2-arachidonoylglycerol and regulates the production of arachidonic acid and prostaglandins-substances that mediate tissue inflammatory response. Here, we have studied the effects of the selective MAGL inhibitors JZL184 and MJN110 and their underlying molecular mechanisms on 3 different experimental models of focal cerebral ischemia. METHODS: SHR (spontaneously hypertensive rats) and normotensive WKY (Wistar Kyoto) rats were subject to an intracortical injection of the potent vasoconstrictor endothelin-1, permanent occlusion of a distal segment of the middle cerebral artery via craniectomy, or transient occlusion of the middle cerebral artery by the intraluminal suture method. JZL184 or MJN110 was administered 60 minutes after focal cerebral ischemia. Infarct volumes, hemispheric swelling, and functional outcomes were assessed between days 1 to 28 by magnetic resonance imaging, histology, and behavioral tests. RESULTS: Pharmacological inhibition of MAGL significantly attenuated infarct volume and hemispheric swelling. MAGL inhibition also ameliorated sensorimotor deficits, suppressed inflammatory response, and decreased the number of degenerating neurons. These beneficial effects of MAGL inhibition were not fully abrogated by selective antagonists of cannabinoid receptors, indicating that the anti-inflammatory effects are caused by inhibition of eicosanoid production rather than by activation of cannabinoid receptors. CONCLUSIONS: Our results suggest that MAGL may contribute to the pathophysiology of focal cerebral ischemia and is thus a promising therapeutic target for the treatment of ischemic stroke.

Investigation of the BOLD and CBV fMRI responses to somatosensory stimulation in awake marmosets (Callithrix jacchus).

Understanding the spatiotemporal features of the hemodynamic response function (HRF) to brain stimulation is essential for the correct application of neuroimaging methods to study brain function. Here, we investigated the spatiotemporal evolution of the blood oxygen level-dependent (BOLD) and cerebral blood volume (CBV) HRF in conscious, awake marmosets (Callithrix jacchus), a New World non-human primate with a lissencephalic brain and with growing use in biomedical research. The marmosets were acclimatized to head fixation and placed in a 7-T magnetic resonance imaging (MRI) scanner. Somatosensory stimulation (333-µs pulses; amplitude, 2 mA; 64 Hz) was delivered bilaterally via pairs of contact electrodes. A block design paradigm was used in which the stimulus duration increased in pseudo-random order from a single pulse up to 256 electrical pulses (4 s). For CBV measurements, 30 mg/kg of ultrasmall superparamagnetic ironoxide particles (USPIO) injected intravenously, were used. Robust BOLD and CBV HRFs were obtained in the primary somatosensory cortex (S1), secondary somatosensory cortex (S2) and caudate at all stimulus conditions. In particular, BOLD and CBV responses to a single 333-µs-long stimulus were reliably measured, and the CBV HRF presented shorter onset time and time to peak than the BOLD HRF. Both the size of the regions of activation and the peak amplitude of the HRFs grew quickly with increasing stimulus duration, and saturated for stimulus durations greater than 1 s. Onset times in S1 and S2 were faster than in caudate. Finally, the fine spatiotemporal features of the HRF in awake marmosets were similar to those obtained in humans, indicating that the continued refinement of awake non-human primate models is essential to maximize the applicability of animal functional MRI studies to the investigation of human brain function.

A digital 3D atlas of the marmoset brain based on multi-modal MRI.

The common marmoset (Callithrix jacchus) is a New-World monkey of growing interest in neuroscience. Magnetic resonance imaging (MRI) is an essential tool to unveil the anatomical and functional organization of the marmoset brain. To facilitate identification of regions of interest, it is desirable to register MR images to an atlas of the brain. However, currently available atlases of the marmoset brain are mainly based on 2D histological data, which are difficult to apply to 3D imaging techniques. Here, we constructed a 3D digital atlas based on high-resolution ex-vivo MRI images, including magnetization transfer ratio (a T1-like contrast), T2w images, and multi-shell diffusion MRI. Based on the multi-modal MRI images, we manually delineated 54 cortical areas and 16 subcortical regions on one hemisphere of the brain (the core version). The 54 cortical areas were merged into 13 larger cortical regions according to their locations to yield a coarse version of the atlas, and also parcellated into 106 sub-regions using a connectivity-based parcellation method to produce a refined atlas. Finally, we compared the new atlas set with existing histology atlases and demonstrated its applications in connectome studies, and in resting state and stimulus-based fMRI. The atlas set has been integrated into the widely-distributed neuroimaging data analysis software AFNI and SUMA, providing a readily usable multi-modal template space with multi-level anatomical labels (including labels from the Paxinos atlas) that can facilitate various neuroimaging studies of marmosets.

Investigating the spatiotemporal characteristics of the deoxyhemoglobin-related and deoxyhemoglobin-unrelated functional hemodynamic response across cortical layers in awake marmosets.

Blood oxygenation level dependent (BOLD) functional magnetic resonance imaging (fMRI) has become a major tool to map neural activity. However, the spatiotemporal characteristics of the BOLD functional hemodynamic response across the cortical layers remain poorly understood. While human fMRI studies suffer from low spatiotemporal resolution, the use of anesthesia in animal models introduces confounding factors. Additionally, inflow contributions to the fMRI signal become non-negligible when short repetition times (TRs) are used. In the present work, we mapped the BOLD fMRI response to somatosensory stimulation in awake marmosets. To address the above technical concerns, we used a dual-echo gradient-recalled echo planar imaging (GR-EPI) sequence to separate the deoxyhemoglobin-related response (absolute T2* differences) from the deoxyhemoglobin-unrelated response (relative S0 changes). We employed a spatial saturation pulse to saturate incoming arterial spins and reduce inflow effects. Functional GR-EPI images were obtained from a single coronal slice with two different echo times (13.5 and 40.5ms) and TR=0.2s. BOLD, T2*, and S0 images were calculated and their functional responses were detected in both hemispheres of primary somatosensory cortex, from which five laminar regions (L1+2, L3, L4, L5, and L6) were derived. The spatiotemporal distribution of the BOLD response across the cortical layers was heterogeneous, with the middle layers having the highest BOLD amplitudes and shortest onset times. ΔT2* also showed a similar trend. However, functional S0 changes were detected only in L1+2, with a fast onset time. Because inflow effects were minimized, the source of S0 functional changes in L1+2 could be attributed to a reduction of cerebrospinal fluid volume fraction due to the functional increase in cerebral blood volume and to unmodeled T2* changes in the extra- and intra-venous compartments. Caution should be exercised when interpreting laminar BOLD fMRI changes in superficial layers as surrogates of underlying neural activity.