High-resolution mapping and digital atlas of subcortical regions in the macaque monkey based on matched MAP-MRI and histology

Deep brain stimulation induces white matter remodeling and functional changes to brain-wide networks

Deep brain stimulation (DBS) is an emerging therapeutic option for treatment resistant neurological and psychiatric disorders, most notably depression. Despite this, little is known about the anatomical and functional mechanisms that underlie this therapy. Here we targeted stimulation to the white matter adjacent to the subcallosal anterior cingulate cortex (SCC-DBS) in macaques, modeling the location in the brain proven effective for depression. We demonstrate that SCC-DBS has a selective effect on white matter macro- and micro-structure in the cingulum bundle distant to where stimulation was delivered. SCC-DBS also decreased functional connectivity between subcallosal and posterior cingulate cortex, two areas linked by the cingulum bundle and implicated in depression. Our data reveal that white matter remodeling as well as functional effects contribute to DBS's therapeutic efficacy.

The neural basis of resting-state fMRI functional connectivity in fronto-limbic circuits revealed by chemogenetic manipulation

Measures of fMRI resting-state functional connectivity (rs-FC) are an essential tool for basic and clinical investigations of fronto-limbic circuits. Understanding the relationship between rs-FC and the underlying patterns of neural activity in these circuits is therefore vital. Here we introduced inhibitory designer receptors exclusively activated by designer drugs (DREADDs) into the amygdala of two male macaques. We evaluated the causal effect of activating the DREADD receptors on rs-FC and neural activity within circuits connecting amygdala and frontal cortex. Activating the inhibitory DREADD increased rs-FC between amygdala and ventrolateral prefrontal cortex. Neurophysiological recordings revealed that the DREADD-induced increase in fMRI rs-FC was associated with increased local field potential coherency in the alpha band (6.5-14.5 Hz) between amygdala and ventrolateral prefrontal cortex. Thus, our multi-modal approach reveals the specific signature of neuronal activity that underlies rs-FC in fronto-limbic circuits.

A complementary approach for neocortical cytoarchitecture inspection with cellular resolution imaging at whole brain scale

Cytoarchitecture, the organization of cells within organs and tissues, serves as a crucial anatomical foundation for the delineation of various regions. It enables the segmentation of the cortex into distinct areas with unique structural and functional characteristics. While traditional 2D atlases have focused on cytoarchitectonic mapping of cortical regions through individual sections, the intricate cortical gyri and sulci demands a 3D perspective for unambiguous interpretation. In this study, we employed fluorescent micro-optical sectioning tomography to acquire architectural datasets of the entire macaque brain at a resolution of 0.65 μm × 0.65 μm × 3 μm. With these volumetric data, the cortical laminar textures were remarkably presented in appropriate view planes. Additionally, we established a stereo coordinate system to represent the cytoarchitectonic information as surface-based tomograms. Utilizing these cytoarchitectonic features, we were able to three-dimensionally parcel the macaque cortex into multiple regions exhibiting contrasting architectural patterns. The whole-brain analysis was also conducted on mice that clearly revealed the presence of barrel cortex and reflected biological reasonability of this method. Leveraging these high-resolution continuous datasets, our method offers a robust tool for exploring the organizational logic and pathological mechanisms of the brain's 3D anatomical structure.

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.

Macaque Brainnetome Atlas: A multifaceted brain map with parcellation, connection, and histology

The rhesus macaque (Macaca mulatta) is a crucial experimental animal that shares many genetic, brain organizational, and behavioral characteristics with humans. A macaque brain atlas is fundamental to biomedical and evolutionary research. However, even though connectivity is vital for understanding brain functions, a connectivity-based whole-brain atlas of the macaque has not previously been made. In this study, we created a new whole-brain map, the Macaque Brainnetome Atlas (MacBNA), based on the anatomical connectivity profiles provided by high angular and spatial resolution ex vivo diffusion MRI data. The new atlas consists of 248 cortical and 56 subcortical regions as well as their structural and functional connections. The parcellation and the diffusion-based tractography were evaluated with invasive neuronal-tracing and Nissl-stained images. As a demonstrative application, the structural connectivity divergence between macaque and human brains was mapped using the Brainnetome atlases of those two species to uncover the genetic underpinnings of the evolutionary changes in brain structure. The resulting resource includes: (1) the thoroughly delineated Macaque Brainnetome Atlas (MacBNA), (2) regional connectivity profiles, (3) the postmortem high-resolution macaque diffusion and T2-weighted MRI dataset (Brainnetome-8), and (4) multi-contrast MRI, neuronal-tracing, and histological images collected from a single macaque. MacBNA can serve as a common reference frame for mapping multifaceted features across modalities and spatial scales and for integrative investigation and characterization of brain organization and function. Therefore, it will enrich the collaborative resource platform for nonhuman primates and facilitate translational and comparative neuroscience research.

Evolutionary continuity and divergence of auditory dorsal and ventral pathways in primates revealed by ultra-high field diffusion MRI

Auditory dorsal and ventral pathways in the human brain play important roles in supporting speech and language processing. However, the evolutionary root of the dual auditory pathways in the primate brain is unclear. By parcellating the auditory cortex of marmosets (a New World monkey species), macaques (an Old World monkey species), and humans using the same individual-based analysis method and tracking the pathways from the auditory cortex based on multi-shell diffusion-weighted MRI (dMRI), homologous auditory dorsal and ventral fiber tracks were identified in these primate species. The ventral pathway was found to be well conserved in all three primate species analyzed but extend to more anterior temporal regions in humans. In contrast, the dorsal pathway showed a divergence between monkey and human brains. First, frontal regions in the human brain have stronger connections to the higher-level auditory regions than to the lower-level auditory regions along the dorsal pathway, while frontal regions in the monkey brain show opposite connection patterns along the dorsal pathway. Second, the left lateralization of the dorsal pathway is only found in humans. Moreover, the connectivity strength of the dorsal pathway in marmosets is more similar to that of humans than macaques. These results demonstrate the continuity and divergence of the dual auditory pathways in the primate brains along the evolutionary path, suggesting that the putative neural networks supporting human speech and language processing might have emerged early in primate evolution.

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

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

Neuronal microstructural changes in the human brain are associated with neurocognitive aging

Gray matter (GM) alterations play a role in aging-related disorders like Alzheimer's disease and related dementias, yet MRI studies mainly focus on macroscopic changes. Although reliable indicators of atrophy, morphological metrics like cortical thickness lack the sensitivity to detect early changes preceding visible atrophy. Our study aimed at exploring the potential of diffusion MRI in unveiling sensitive markers of cortical and subcortical age-related microstructural changes and assessing their associations with cognitive and behavioral deficits. We leveraged the Human Connectome Project-Aging cohort that included 707 unimpaired participants (394 female; median age = 58, range = 36-90 years) and applied the powerful mean apparent diffusion propagator model to measure microstructural parameters, along with comprehensive behavioral and cognitive test scores. Both macro- and microstructural GM characteristics were strongly associated with age, with widespread significant microstructural correlations reflective of cellular morphological changes, reduced cellular density, increased extracellular volume, and increased membrane permeability. Importantly, when correlating MRI and cognitive test scores, our findings revealed no link between macrostructural volumetric changes and neurobehavioral performance. However, we found that cellular and extracellular alterations in cortical and subcortical GM regions were associated with neurobehavioral performance. Based on these findings, it is hypothesized that increased microstructural heterogeneity and decreased neurite orientation dispersion precede macrostructural changes, and that they play an important role in subsequent cognitive decline. These alterations are suggested to be early markers of neurocognitive performance that may distinctly aid in identifying the mechanisms underlying phenotypic aging and subsequent age-related functional decline.

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

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

Model incorporating multiple diffusion MRI features: development and validation of a radiomics-based model to predict adult-type diffuse gliomas grade

OBJECTIVES

To develop and validate a radiomics-based model (ADGGIP) for predicting adult-type diffuse gliomas (ADG) grade by combining multiple diffusion modalities and clinical and imaging morphologic features.

METHODS

In this prospective study, we recruited 103 participants diagnosed with ADG and collected their preoperative conventional MRI and multiple diffusion imaging (diffusion tensor imaging, diffusion kurtosis imaging, neurite orientation dispersion and density imaging, and mean apparent propagator diffusion-MRI) data in our hospital, as well as clinical information. Radiomic features of the diffusion images and clinical information and morphological data from the radiological reports were extracted, and multiple pipelines were used to construct the optimal model. Model validation was performed through a time-independent validation cohort. ROC curves were used to evaluate model performance. The clinical benefit was determined by decision curve analysis.

RESULTS

From June 2018 to May 2021, 72 participants were recruited for the training cohort. Between June 2021 and February 2022, 31 participants were enrolled in the prospective validation cohort. In the training cohort (AUC 0.958), internal validation cohort (0.942), and prospective validation cohort (0.880), ADGGIP had good accuracy in predicting ADG grade. ADGGIP was also significantly better than the single-modality prediction model (AUC 0.860) and clinical imaging morphology model (0.841) (all p < .01) in the prospective validation cohort. When the threshold probability was greater than 5%, ADGGIP provided the greatest net benefit.

CONCLUSION

ADGGIP, which is based on advanced diffusion modalities, can predict the grade of ADG with high accuracy and robustness and can help improve clinical decision-making.

CLINICAL RELEVANCE STATEMENT

Integrated multi-modal predictive modeling is beneficial for early detection and treatment planning of adult-type diffuse gliomas, as well as for investigating the genuine clinical significance of biomarkers.

KEY POINTS

• Integrated model exhibits the highest performance and stability. • When the threshold is greater than 5%, the integrated model has the greatest net benefit. • The advanced diffusion models do not demonstrate better performance than the simple technology.

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.

The neural basis of resting-state fMRI functional connectivity in fronto-limbic circuits revealed by chemogenetic manipulation

Measures of fMRI resting-state functional connectivity (rs-FC) are an essential tool for basic and clinical investigations of fronto-limbic circuits. Understanding the relationship between rs-FC and neural activity in these circuits is therefore vital. Here we introduced inhibitory designer receptors exclusively activated by designer drugs (DREADDs) into the macaque amygdala and activated them with a highly selective and potent DREADD agonist, deschloroclozapine. We evaluated the causal effect of activating the DREADD receptors on rs-FC and neural activity within circuits connecting amygdala and frontal cortex. Interestingly, activating the inhibitory DREADD increased rs-FC between amygdala and ventrolateral prefrontal cortex. Neurophysiological recordings revealed that the DREADD-induced increase in fMRI rs-FC was associated with increased local field potential coherency in the alpha band (6.5-14.5Hz) between amygdala and ventrolateral prefrontal cortex. Thus, our multi-disciplinary approach reveals the specific signature of neuronal activity that underlies rs-FC in fronto-limbic circuits.

Opposing functions of glutamatergic inputs between the globus pallidus external segment and substantia nigra pars reticulata

The indirect pathway of the basal ganglia, including the subthalamic nucleus (STN) and globus pallidus external segment (GPe), is believed to play a crucial role in suppressing involuntary movements. However, recent evidence suggests the STN and GPe also facilitate voluntary movements. This study hypothesized that excitatory inputs from the STN to the GPe contribute to this facilitation, and that excitatory projections to the substantia nigra pars reticulata (SNr) are involved in the inhibition. To disrupt the STN-GPe or STN-SNr projections in monkeys during choice and fixation tasks, glutamate receptor inhibitors were injected into the GPe or SNr, which induced delayed saccade latencies toward good choices in the choice task (GPe) and caused frequent reflexive saccades to objects in the fixation task (SNr). Our findings suggest excitatory inputs to the GPe and SNr work in opposing manners, providing new insights that redefine our understanding of the functions of basal ganglia pathways.

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

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

Adult lifespan maturation and degeneration patterns in gray and white matter: A mean apparent propagator (MAP) MRI study

The relationship between brain microstructure and aging has been the subject of intense study, with diffusion MRI perhaps the most effective modality for elucidating these associations. Here, we used the mean apparent propagator (MAP)-MRI framework, which is suitable to characterize complex microstructure, to investigate age-related cerebral differences in a cohort of cognitively unimpaired participants and compared the results to those derived using diffusion tensor imaging. We studied MAP-MRI metrics, among them the non-Gaussianity (NG) and propagator anisotropy (PA), and established an opposing pattern in white matter of higher NG alongside lower PA among older adults, likely indicative of axonal degradation. In gray matter, however, these two indices were consistent with one another, and exhibited regional pattern heterogeneity compared to other microstructural parameters, which could indicate fewer neuronal projections across cortical layers along with an increased glial concentration. In addition, we report regional variations in the magnitude of age-related microstructural differences consistent with the posterior-anterior shift in aging paradigm. These results encourage further investigations in cognitive impairments and neurodegeneration.

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.

High-throughput sequencing of macaque basolateral amygdala projections reveals dissociable connectional motifs with frontal cortex

The basolateral amygdala (BLA) projects widely across the macaque frontal cortex1-4, and amygdalo-frontal projections are critical for optimal emotional responding5 and decision-making6. Yet, little is known about the single-neuron architecture of these projections: namely, whether single BLA neurons project to multiple parts of the frontal cortex. Here, we use MAPseq7 to determine the projection patterns of over 3000 macaque BLA neurons. We found that one-third of BLA neurons have two or more distinct targets in parts of frontal cortex and of subcortical structures. Further, we reveal non-random structure within these branching patterns such that neurons with four targets are more frequently observed than those with two or three, indicative of widespread networks. Consequently, these multi-target single neurons form distinct networks within medial and ventral frontal cortex consistent with their known functions in regulating mood and decision-making. Additionally, we show that branching patterns of single neurons shape functional networks in the brain as assessed by fMRI-based functional connectivity. These results provide a neuroanatomical basis for the role of the BLA in coordinating brain-wide responses to valent stimuli8 and highlight the importance of high-resolution neuroanatomical data for understanding functional networks in the brain.

A multimodal imaging-guided software for access to primate brains

Background

Imaging-guided access to the brain has become a routine procedure for various research and clinical applications, including drug administration, neurophysiological recording, and sampling tissue. Therefore, open-source software is required to handle such datasets in these specific applications.

New methods

Here, we proposed an open-source tool utilizing different imaging modalities for automating the steps to access the brain. This tool provides means for easily calculating the coordination of the area of interest concerning a specific point of entry. The source and documentation are available at this link.

Results

We have used this software for three different applications: electrophysiological recording, drug infusion in the nonhuman primate brain, and guided biopsy procedure in the human brain. We performed a neural recording of two monkeys' prefrontal cortex and inferior temporal cortex using this software in submillimeter resolution. We also applied our procedure for infusion in the putamen and caudate nuclei in both hemispheres of another group of rhesus monkeys with histological proof in one animal. More so, we validated this software in the human subjects that underwent biopsy surgery with the commercial software used in human biopsy surgery.

Comparison with existing methods

Our software uses different imaging modalities by co-registering them. This will provide structural details of the skull and brain tissue. We can calculate each brain region's coordination at the point of entry by re-slicing the images. Atlas-based image segmentation were implemented in our software. Three mentioned applications of our software in neuroscience will be further discussed in this paper.

Conclusion

In our procedure, working with different imaging modalities provides a precise estimation of the specific region in the brain related to the location of implants or stereotaxic frames. There is no limitation to using metal implants in this procedure.

Structural Connectivity Changes After Fornix Transection in Macaques Using Probabilistic Diffusion Tractography

The fornix, the limbic system's white matter tract connecting the extended hippocampal system to subcortical structures of the medial diencephalon, is strongly associated with learning and memory in humans and nonhuman primates (NHPs). Here, we sought to investigate alterations in structural connectivity across key cortical and subcortical regions after fornix transection in NHPs. We collected diffusion-weighted MRI (dMRI) data from three macaque monkeys that underwent bilateral fornix transection during neurosurgery and from four age- and cohort-matched control macaques that underwent surgery to implant a head-post but remained neurologically intact. dMRI data were collected from both groups at two time points, before and after the surgeries, and scans took place at around the same time for the two groups. We used probabilistic tractography and employed the number of tracking streamlines to quantify connectivity across our regions of interest (ROIs), in all dMRI sessions. In the neurologically intact monkeys, we observed high connectivity across certain ROIs, including the CA3 hippocampal subfield with the retrosplenial cortex (RSC), the anterior thalamus with the RSC, and the RSC with the anterior cingulate cortex (ACC). However, we found that, compared to the control group, the fornix-transected monkeys showed marked, significant, connectivity changes including increases between the anterior thalamus and the ACC and between the CA3 and the ACC, as well as decreases between the CA3 and the RSC. Our results highlight cortical and subcortical network changes after fornix transection and identify candidate indirect connectivity routes that may support memory functions after damage and/or neurodegeneration.

COnstrained Reference frame diffusion TEnsor Correlation Spectroscopic (CORTECS) MRI: A practical framework for high-resolution diffusion tensor distribution imaging

High-resolution imaging studies have consistently shown that in cortical tissue water diffuses preferentially along radial and tangential orientations with respect to the cortical surface, in agreement with histology. These dominant orientations do not change significantly even if the relative contributions from microscopic water pools to the net voxel signal vary across experiments that use different diffusion times, b-values, TEs, and TRs. With this in mind, we propose a practical new framework for imaging non-parametric diffusion tensor distributions (DTDs) by constraining the microscopic diffusion tensors of the DTD to be diagonalized using the same orthonormal reference frame of the mesoscopic voxel. In each voxel, the constrained DTD (cDTD) is completely determined by the correlation spectrum of the microscopic principal diffusivities associated with the axes of the voxel reference frame. Consequently, all cDTDs are inherently limited to the domain of positive definite tensors and can be reconstructed efficiently using Inverse Laplace Transform methods. Moreover, the cDTD reconstruction can be performed using only data acquired efficiently with single diffusion encoding, although it also supports datasets with multiple diffusion encoding. In tissues with a well-defined architecture, such as the cortex, we can further constrain the cDTD to contain only cylindrically symmetric diffusion tensors and measure the 2D correlation spectra of principal diffusivities along the radial and tangential orientation with respect to the cortical surface. To demonstrate this framework, we perform numerical simulations and analyze high-resolution dMRI data from a fixed macaque monkey brain. We estimate 2D cDTDs in the cortex and derive, in each voxel, the marginal distributions of the microscopic principal diffusivities, the corresponding distributions of the microscopic fractional anisotropies and mean diffusivities along with their 2D correlation spectra to quantify the cDTD shape-size characteristics. Signal components corresponding to specific bands in these cDTD-derived spectra show high specificity to cortical laminar structures observed with histology. Our framework drastically simplifies the measurement of non-parametric DTDs in high-resolution datasets with mesoscopic voxel sizes much smaller than the radius of curvature of the underlying anatomy, e.g., cortical surface, and can be applied retrospectively to analyze existing diffusion MRI data from fixed cortical tissues.

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.

Fusion of quantitative susceptibility maps and T1-weighted images improve brain tissue contrast in primates

Recent progress in quantitative susceptibility mapping (QSM) has enabled the accurate delineation of submillimeter-scale subcortical brain structures in humans. However, the simultaneous visualization of cortical, subcortical, and white matter structure remains challenging, utilizing QSM data solely. Here we present TQ-SILiCON, a fusion method that enhances the contrast of cortex and subcortical structures and provides an excellent white matter delineation by combining QSM and conventional T1-weighted (T1w) images. In this study, we first applied QSM in the macaque monkey to map iron-rich subcortical structures. Implementing the same QSM acquisition and analysis methods allowed a similar accurate delineation of subcortical structures in humans. However, the QSM contrast of white and cortical gray matter was not sufficient for appropriate segmentation. Applying automatic brain tissue segmentation to TQ-SILiCON images of the macaque improved the classification of subcortical brain structures as compared to the single T1 contrast by maintaining an excellent white to cortical gray matter contrast. Furthermore, we validated our dual-contrast fusion approach in humans and similarly demonstrated improvements in automated segmentation of the cortex and subcortical structures. We believe the proposed contrast will facilitate translational studies in nonhuman primates to investigate the pathophysiology of neurodegenerative diseases that affect subcortical structures such as the basal ganglia in humans.

Coexistence of chronic hyperalgesia and multilevel neuroinflammatory responses after experimental SCI: a systematic approach to profiling neuropathic pain

Background

People with spinal cord injury (SCI) frequently develop neuropathic pain (NP) that worsens disability and diminishes rehabilitation efficacy. Chronic NP is presently incurable due to poor understanding of underlying mechanisms. We hypothesized that multilocus neuroinflammation (NIF) might be a driver of SCI NP, and tested it by investigating whether NP coexisted with central NIF, neurotransmission (NTM), neuromodulation (NML) and neuroplasticity (NPL) changes post-SCI.

Methods

Female Sprague–Dawley rats (230–250 g) with T10 compression or laminectomy were evaluated for physical conditions, coordinated hindlimb functions, neurological reflexes, and mechanical/thermal sensitivity thresholds at 1 day post-injury (p.i.) and weekly thereafter. Eight weeks p.i., central nervous system tissues were histochemically and immunohistochemically characterized for parameters/markers of histopathology and NIF/NTM/NML/NPL. Also analyzed was the correlative relationship between levels of selected biomarkers and thermosensitivity thresholds via statistical linear regression.

Results

SCI impaired sensorimotor functions, altered reflexes, and produced spontaneous pain signs and hypersensitivity to evoked nociceptive, mechanical, and thermal inputs. Only injured spinal cords exhibited neural lesion, microglia/astrocyte activation, and abnormal expression of proinflammatory cytokines, as well as NIF/NTM/NML/NPL markers. Brains of SCI animals displayed similar pathophysiological signs in the gracile and parabrachial nuclei (GrN and PBN: sensory relay), raphe magnus nucleus and periaqueduct gray (RMN and PAG: pain modulation), basolateral amygdala (BLA: emotional-affective dimension of pain), and hippocampus (HPC: memory/mood/neurogenesis). SCI augmented sensory NTM/NPL (GrN and PBN); increased GAD67 (PAG) level; reduced serotonin (RMN) and fear-off neuronal NTR2 (BLA) expressions; and perturbed neurogenesis (HPC).

Conclusion

T10 compression caused chronic hyperalgesia that coexisted with NIF/NTM/NML/NPL responses at multilevel neuroaxis centers. The data have provided multidimensional biomarkers as new mechanistic leads to profile SCI NP for therapeutic/therapy development.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12974-022-02628-2.

Corresponding anatomical of the macaque superior parietal lobule areas 5 (PE) subdivision reveal similar connectivity patterns with humans

Using the animal brain as a cross-species tool for human brain research based on imaging features can provide more potential to reveal comprehensive human brain analysis. Previous studies have shown that human Brodmann area 5 (BA5) and macaque PE are homologous regions. They are both involved in processes depth and direction information during the touch process in the arm movement. However, recent studies show that both BA5 and PE are not homogeneous. According to the cytoarchitecture, BA5 is subdivided into three different subregions, and PE can be subdivided into PEl, PEla, and PEm. The species homologous relationship among the subregions is not clear between BA5 and PE. At the same time, the subdivision of PE based on the anatomical connection of white matter fiber bundles needs more verification. This research subdivided the PE of macaques based on the anatomical connection of white matter fiber bundles. Two PE subregions are defined based on probabilistic fiber tracking, one on the anterior side and the other on the dorsal side. Finally, the research draws connectivity fingerprints with predefined homologous target areas for the BA5 and PE subregions to reveal the characteristics of structure and functions and gives the homologous correspondence identified.

Cumulant expansion with localization: A new representation of the diffusion MRI signal

Diffusion MR is sensitive to the microstructural features of a sample. Fine-scale characteristics can be probed by employing strong diffusion gradients while the low b-value regime is determined by the cumulants of the distribution of particle displacements. A signal representation based on the cumulants, however, suffers from a finite convergence radius and cannot represent the 'localization regime' characterized by a stretched exponential decay that emerges at large gradient strengths. Here, we propose a new representation for the diffusion MR signal. Our method provides not only a robust estimate of the first three cumulants but also a meaningful extrapolation of the entire signal decay.

Anatomical variability, multi-modal coordinate systems, and precision targeting in the marmoset brain

Localising accurate brain regions needs careful evaluation in each experimental species due to their individual variability. However, the function and connectivity of brain areas is commonly studied using a single-subject cranial landmark-based stereotactic atlas in animal neuroscience. Here, we address this issue in a small primate, the common marmoset, which is increasingly widely used in systems neuroscience. We developed a non-invasive multi-modal neuroimaging-based targeting pipeline, which accounts for intersubject anatomical variability in cranial and cortical landmarks in marmosets. This methodology allowed creation of multi-modal templates (MarmosetRIKEN20) including head CT and brain MR images, embedded in coordinate systems of anterior and posterior commissures (AC-PC) and CIFTI grayordinates. We found that the horizontal plane of the stereotactic coordinate was significantly rotated in pitch relative to the AC-PC coordinate system (10 degrees, frontal downwards), and had a significant bias and uncertainty due to positioning procedures. We also found that many common cranial and brain landmarks (e.g., bregma, intraparietal sulcus) vary in location across subjects and are substantial relative to average marmoset cortical area dimensions. Combining the neuroimaging-based targeting pipeline with robot-guided surgery enabled proof-of-concept targeting of deep brain structures with an accuracy of 0.2 mm. Altogether, our findings demonstrate substantial intersubject variability in marmoset brain and cranial landmarks, implying that subject-specific neuroimaging-based localization is needed for precision targeting in marmosets. The population-based templates and atlases in grayordinates, created for the first time in marmoset monkeys, should help bridging between macroscale and microscale analyses.