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

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

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.

A comprehensive macaque fMRI pipeline and hierarchical atlas.

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

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.

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

Digitized neuroanatomical atlases that can be overlaid onto functional data are crucial for localizing brain structures and analyzing functional networks identified by neuroimaging techniques. To aid in functional and structural data analysis, we have created a comprehensive parcellation of the rhesus macaque subcortex using a high-resolution ex vivo structural imaging scan. This anatomical scan and its parcellation were warped to the updated NIMH Macaque Template (NMT v2), an in vivo population template, where the parcellation was refined to produce the Subcortical Atlas of the Rhesus Macaque (SARM) with 210 primary regions-of-interest (ROIs). The subcortical parcellation and nomenclature reflect those of the 4th edition of the Rhesus Monkey Brain in Stereotaxic Coordinates (Paxinos et al., in preparation), rather than proposing yet another novel atlas. The primary ROIs are organized across six spatial hierarchical scales from small, fine-grained ROIs to broader composites of multiple ROIs, making the SARM suitable for analysis at different resolutions and allowing broader labeling of functional signals when more accurate localization is not possible. As an example application of this atlas, we have included a functional localizer for the dorsal lateral geniculate (DLG) nucleus in three macaques using a visual flickering checkerboard stimulus, identifying and quantifying significant fMRI activation in this atlas region. The SARM has been made openly available to the neuroimaging community and can easily be used with common MRI data processing software, such as AFNI, where the atlas has been embedded into the software alongside cortical macaque atlases.

LayNii: A software suite for layer-fMRI.

High-resolution fMRI in the sub-millimeter regime allows researchers to resolve brain activity across cortical layers and columns non-invasively. While these high-resolution data make it possible to address novel questions of directional information flow within and across brain circuits, the corresponding data analyses are challenged by MRI artifacts, including image blurring, image distortions, low SNR, and restricted coverage. These challenges often result in insufficient spatial accuracy of conventional analysis pipelines. Here we introduce a new software suite that is specifically designed for layer-specific functional MRI: LayNii. This toolbox is a collection of command-line executable programs written in C/C++ and is distributed opensource and as pre-compiled binaries for Linux, Windows, and macOS. LayNii is designed for layer-fMRI data that suffer from SNR and coverage constraints and thus cannot be straightforwardly analyzed in alternative software packages. Some of the most popular programs of LayNii contain 'layerification' and columnarization in the native voxel space of functional data as well as many other layer-fMRI specific analysis tasks: layer-specific smoothing, model-based vein mitigation of GE-BOLD data, quality assessment of artifact dominated sub-millimeter fMRI, as well as analyses of VASO data.

A collaborative resource platform for non-human primate neuroimaging.

Neuroimaging non-human primates (NHPs) is a growing, yet highly specialized field of neuroscience. Resources that were primarily developed for human neuroimaging often need to be significantly adapted for use with NHPs or other animals, which has led to an abundance of custom, in-house solutions. In recent years, the global NHP neuroimaging community has made significant efforts to transform the field towards more open and collaborative practices. Here we present the PRIMatE Resource Exchange (PRIME-RE), a new collaborative online platform for NHP neuroimaging. PRIME-RE is a dynamic community-driven hub for the exchange of practical knowledge, specialized analytical tools, and open data repositories, specifically related to NHP neuroimaging. PRIME-RE caters to both researchers and developers who are either new to the field, looking to stay abreast of the latest developments, or seeking to collaboratively advance the field .

The Haskins pediatric atlas: a magnetic-resonance-imaging-based pediatric template and atlas.

BACKGROUND: Spatial normalization plays an essential role in multi-subject MRI and functional MRI (fMRI) experiments by facilitating a common space in which group analyses are performed. Although many prominent adult templates are available, their use for pediatric data is problematic. Generalized templates for pediatric populations are limited or constructed using older methods that result in less ideal normalization. OBJECTIVE: The Haskins pediatric templates and atlases aim to provide superior registration and more precise accuracy in labeling of anatomical and functional regions essential for all fMRI studies involving pediatric populations. MATERIALS AND METHODS: The Haskins pediatric templates and atlases were generated with nonlinear methods using structural MRI from 72 children (age range 7-14 years, median 10 years), allowing for a detailed template with corresponding parcellations of labeled atlas regions. The accuracy of these templates and atlases was assessed using multiple metrics of deformation distance and overlap. RESULTS: When comparing the deformation distances from normalizing pediatric data between this template and both the adult templates and other pediatric templates, we found significantly less deformation distance for the Haskins pediatric template (P<0.0001). Further, the correct atlas classification was higher using the Haskins pediatric template in 74% of regions (P<0.0001). CONCLUSION: The Haskins pediatric template results in more accurate correspondence across subjects because of lower deformation distances. This correspondence also provides better accuracy in atlas locations to benefit structural and functional imaging analyses of pediatric populations.

Sub-millimeter fMRI reveals multiple topographical digit representations that form action maps in human motor cortex.

The human brain coordinates a wide variety of motor activities. On a large scale, the cortical motor system is topographically organized such that neighboring body parts are represented by neighboring brain areas. This homunculus-like somatotopic organization along the central sulcus has been observed using neuroimaging for large body parts such as the face, hands and feet. However, on a finer scale, invasive electrical stimulation studies show deviations from this somatotopic organization that suggest an organizing principle based on motor actions rather than body part moved. It has not been clear how the action-map organization principle of the motor cortex in the mesoscopic (sub-millimeter) regime integrates into a body map organization principle on a macroscopic scale (cm). Here we developed and applied advanced mesoscopic (sub-millimeter) fMRI and analysis methodology to non-invasively investigate the functional organization topography across columnar and laminar structures in humans. Compared to previous methods, in this study, we could capture locally specific blood volume changes across entire brain regions along the cortical curvature. We find that individual fingers have multiple mirrored representations in the primary motor cortex depending on the movements they are involved in. We find that individual digits have cortical representations up to 3 â€‹mm apart from each other arranged in a column-like fashion. These representations are differentially engaged depending on whether the digits' muscles are used for different motor actions such as flexion movements, like grasping a ball or retraction movements like releasing a ball. This research provides a starting point for non-invasive investigation of mesoscale topography across layers and columns of the human cortex and bridges the gap between invasive electrophysiological investigations and large coverage non-invasive neuroimaging.

A series of five population-specific Indian brain templates and atlases spanning ages 6-60 years.

Anatomical brain templates are commonly used as references in neurological MRI studies, for bringing data into a common space for group-level statistics and coordinate reporting. Given the inherent variability in brain morphology across age and geography, it is important to have templates that are as representative as possible for both age and population. A representative-template increases the accuracy of alignment, decreases distortions as well as potential biases in final coordinate reports. In this study, we developed and validated a new set of T1w Indian brain templates (IBT) from a large number of brain scans (total n = 466) acquired across different locations and multiple 3T MRI scanners in India. A new tool in AFNI, make_template_dask.py, was created to efficiently make five age-specific IBTs (ages 6-60 years) as well as maximum probability map (MPM) atlases for each template; for each age-group's template-atlas pair, there is both a "population-average" and a "typical" version. Validation experiments on an independent Indian structural and functional-MRI dataset show the appropriateness of IBTs for spatial normalization of Indian brains. The results indicate significant structural differences when comparing the IBTs and MNI template, with these differences being maximal along the Anterior-Posterior and Inferior-Superior axes, but minimal Left-Right. For each age-group, the MPM brain atlases provide reasonably good representation of the native-space volumes in the IBT space, except in a few regions with high intersubject variability. These findings provide evidence to support the use of age and population-specific templates in human brain mapping studies.

A tail of two sides: Artificially doubled false positive rates in neuroimaging due to the sidedness choice with t-tests.

One-sided t-tests are widely used in neuroimaging data analysis. While such a test may be applicable when investigating specific regions and prior information about directionality is present, we argue here that it is often mis-applied, with severe consequences for false positive rate (FPR) control. Conceptually, a pair of one-sided t-tests conducted in tandem (e.g., to test separately for both positive and negative effects), effectively amounts to a two-sided t-test. However, replacing the two-sided test with a pair of one-sided tests without multiple comparisons correction essentially doubles the intended FPR of statements made about the same study; that is, the actual family-wise error (FWE) of results at the whole brain level would be 10% instead of the 5% intended by the researcher. Therefore, we strongly recommend that, unless otherwise explicitly justified, two-sided t-tests be applied instead of two simultaneous one-sided t-tests.

An integrative Bayesian approach to matrix-based analysis in neuroimaging.

Understanding the correlation structure associated with brain regions is a central goal in neuroscience, as it informs about interregional relationships and network organization. Correlation structure can be conveniently captured in a matrix that indicates the relationships among brain regions, which could involve electroencephalogram sensors, electrophysiology recordings, calcium imaging data, or functional magnetic resonance imaging (FMRI) data-We call this type of analysis matrix-based analysis, or MBA. Although different methods have been developed to summarize such matrices across subjects, including univariate general linear models (GLMs), the available modeling strategies tend to disregard the interrelationships among the regions, leading to "inefficient" statistical inference. Here, we develop a Bayesian multilevel (BML) modeling framework that simultaneously integrates the analyses of all regions, region pairs (RPs), and subjects. In this approach, the intricate relationships across regions as well as across RPs are quantitatively characterized. The adoption of the Bayesian framework allows us to achieve three goals: (a) dissolve the multiple testing issue typically associated with seeking evidence for the effect of each RP under the conventional univariate GLM; (b) make inferences on effects that would be treated as "random" under the conventional linear mixed-effects framework; and (c) estimate the effect of each brain region in a manner that indexes their relative "importance". We demonstrate the BML methodology with an FMRI dataset involving a cognitive-emotional task and compare it to the conventional GLM approach in terms of model efficiency, performance, and inferences. The associated program MBA is available as part of the AFNI suite for general use.

Identifying perfusion deficits on CT perfusion images using temporal similarity perfusion (TSP) mapping.

OBJECTIVES: Deconvolution-derived maps of CT perfusion (CTP) data may be confounded by transit delays. We propose temporal similarity perfusion (TSP) analysis to decrease CTP maps' dependence on transit times and investigate its sensitivity to detect perfusion deficits. METHODS: CTP data of acute stroke patients obtained within 9 h of symptom onset was analyzed using a delay-insensitive singular value decomposition method and with TSP. The TSP method applies an iterative process whereby a pixel's highest Pearson's R value is obtained through comparison of a pixel's time-shifted signal density time-series curve and the average whole brain signal density time-series curve. Our evaluation included a qualitative and quantitative rating of deconvolution maps (MTT, CBV, and TTP), of TSP maps, and of follow-up CT. RESULTS: Sixty-five patients (mean 68 (SD 13) years, 34 male) were included. A perfusion deficit was identified in 90%, 86%, 65%, and 84% of MTT, TTP, CBV, and TSP maps. The agreement of MTT, TTP, and TSP with CT follow-up was comparable but noticeably lower for CBV. CBV had the best relationship with final infarct volume (R2 = 0.77, p < 0.001), followed by TSP (R2 = 0.63, p < 0.001). Intra-rater agreement of an inexperienced reader was higher for TSP than for CBV/MTT maps (kappa's of 0.79-0.84 and 0.63-0.7). Inter-rater agreement for experienced readers was comparable across maps. CONCLUSIONS: TSP maps are easier to interpret for inexperienced readers. Perfusion deficits detected by TSP are smaller which may suggest less dependence on transit delays although more investigation is required. KEY POINTS: • Temporal similarity perfusion mapping assesses CTP data based on similarities in signal time-curves. • TSP maps are comparable in perfusion deficit detection to deconvolution maps. • TSP maps are easier to interpret for inexperienced readers.

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.

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.

Statistical power comparisons at 3T and 7T with a GO / NOGO task.

The field of cognitive neuroscience is weighing evidence about whether to move from standard field strength to ultra-high field (UHF). The present study contributes to the evidence by comparing a cognitive neuroscience paradigm at 3 Tesla (3T) and 7 Tesla (7T). The goal was to test and demonstrate the practical effects of field strength on a standard GO/NOGO task using accessible preprocessing and analysis tools. Two independent matched healthy samples (N = 31 each) were analyzed at 3T and 7T. Results show gains at 7T in statistical strength, the detection of smaller effects and group-level power. With an increased availability of UHF scanners, these gains may be exploited by cognitive neuroscientists and other neuroimaging researchers to develop more efficient or comprehensive experimental designs and, given the same sample size, achieve greater statistical power at 7T.

Three-Dimensional Digital Template Atlas of the Macaque Brain.

We present a new 3D template atlas of the anatomical subdivisions of the macaque brain, which is based on and aligned to the magnetic resonance imaging (MRI) data set and histological sections of the Saleem and Logothetis atlas. We describe the creation and validation of the atlas that, when registered with macaque structural or functional MRI scans, provides a straightforward means to estimate the boundaries between architectonic areas, either in a 3D volume with different planes of sections, or on an inflated brain surface (cortical flat map). As such, this new template atlas is intended for use as a reference standard for macaque brain research. Atlases and templates are available as both volumes and surfaces in standard NIFTI and GIFTI formats.

Untangling the relatedness among correlations, part I: Nonparametric approaches to inter-subject correlation analysis at the group level.

FMRI data acquisition under naturalistic and continuous stimuli (e.g., watching a video or listening to music) has become popular recently due to the fact that it entails less manipulation and more realistic/complex contexts involved in the task, compared to the conventional task-based experimental designs. The synchronization or response similarities among subjects are typically measured through inter-subject correlation (ISC) between any pair of subjects. At the group level, summarizing the collection of ISC values is complicated by their intercorrelations, which necessarily lead to the violation of independence assumed in typical parametric approaches such as Student's t-test. Nonparametric methods, such as bootstrapping and permutation testing, have previously been adopted for testing purposes by resampling the time series of each subject, but the quantitative validity of these specific approaches in terms of controllability of false positive rate (FPR) has never been explored before. Here we survey the methods of ISC group analysis that have been employed in the literature, and discuss the issues involved in those methods. We then propose less computationally intensive nonparametric methods that can be performed at the group level (for both one- and two-sample analyses), as compared to the popular method of circularly shifting the EPI time series at the individual level. As part of the new approaches, subject-wise (SW) resampling is adopted instead of element-wise (EW) resampling, so that exchangeability and independence assumptions are satisfied, and the patterned correlation structure among the ISC values can be more accurately captured. We examine the FPR controllability and power achievement of all the methods through simulations, as well as their performance when applied to a real experimental dataset.