Mapping vascular network architecture in primate brain using ferumoxytol-weighted laminar MRI.

Mapping the vascular organization of the brain is of great importance across various domains of basic neuroimaging research, diagnostic radiology, and neurology. However, the intricate task of precisely mapping vasculature across brain regions and cortical layers presents formidable challenges, resulting in a limited understanding of neurometabolic factors influencing the brain's microvasculature. Addressing this gap, our study investigates whole-brain vascular volume using ferumoxytol-weighted laminar-resolution multi-echo gradient-echo imaging in macaque monkeys. We validate the results with published data for vascular densities and compare them with cytoarchitecture, neuron and synaptic densities. The ferumoxytol-induced change in transverse relaxation rate (ΔR2*), an indirect proxy measure of cerebral blood volume (CBV), was mapped onto twelve equivolumetric laminar cortical surfaces. Our findings reveal that CBV varies 3-fold across the brain, with the highest vascular volume observed in the inferior colliculus and lowest in the corpus callosum. In the cerebral cortex, CBV is notably high in early primary sensory areas and low in association areas responsible for higher cognitive functions. Classification of CBV into distinct groups unveils extensive replication of translaminar vascular network motifs, suggesting distinct computational energy supply requirements in areas with varying cytoarchitecture types. Regionally, baseline R2* and CBV exhibit positive correlations with neuron density and negative correlations with receptor densities. Adjusting image resolution based on the critical sampling frequency of penetrating cortical vessels, allows us to delineate approximately 30% of the arterial-venous vessels. Collectively, these results mark significant methodological and conceptual advancements, contributing to the refinement of cerebrovascular MRI. Furthermore, our study establishes a linkage between neurometabolic factors and the vascular network architecture in the primate brain.

Cortical adaptation of the night monkey to a nocturnal niche environment: a comparative non-invasive T1w/T2w myelin study.

Night monkeys (Aotus) are the only genus of monkeys within the Simian lineage that successfully occupy a nocturnal environmental niche. Their behavior is supported by their sensory organs' distinctive morphological features; however, little is known about their evolutionary adaptations in sensory regions of the cerebral cortex. Here, we investigate this question by exploring the cortical organization of night monkeys using high-resolution in-vivo brain MRI and comparative cortical-surface T1w/T2w myeloarchitectonic mapping. Our results show that the night monkey cerebral cortex has a qualitatively similar but quantitatively different pattern of cortical myelin compared to the diurnal macaque and marmoset monkeys. T1w/T2w myelin and its gradient allowed us to parcellate high myelin areas, including the middle temporal complex (MT +) and auditory cortex, and a low-myelin area, Brodmann area 7 (BA7) in the three species, despite species differences in cortical convolutions. Relative to the total cortical-surface area, those of MT + and the auditory cortex are significantly larger in night monkeys than diurnal monkeys, whereas area BA7 occupies a similar fraction of the cortical sheet in all three species. We propose that the selective expansion of sensory areas dedicated to visual motion and auditory processing in night monkeys may reflect cortical adaptations to a nocturnal environment.

Empirical transmit field bias correction of T1w/T2w myelin maps.

T1-weighted divided by T2-weighted (T1w/T2w) myelin maps were initially developed for neuroanatomical analyses such as identifying cortical areas, but they are increasingly used in statistical comparisons across individuals and groups with other variables of interest. Existing T1w/T2w myelin maps contain radiofrequency transmit field (B1+) biases, which may be correlated with these variables of interest, leading to potentially spurious results. Here we propose two empirical methods for correcting these transmit field biases using either explicit measures of the transmit field or alternatively a 'pseudo-transmit' approach that is highly correlated with the transmit field at 3T. We find that the resulting corrected T1w/T2w myelin maps are both better neuroanatomical measures (e.g., for use in cross-species comparisons), and more appropriate for statistical comparisons of relative T1w/T2w differences across individuals and groups (e.g., sex, age, or body-mass-index) within a consistently acquired study at 3T. We recommend that investigators who use the T1w/T2w approach for mapping cortical myelin use these B1+ transmit field corrected myelin maps going forward.

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.

Comparative connectomics of the primate social brain.

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

Minimal specifications for non-human primate MRI: Challenges in standardizing and harmonizing data collection.

Recent methodological advances in MRI have enabled substantial growth in neuroimaging studies of non-human primates (NHPs), while open data-sharing through the PRIME-DE initiative has increased the availability of NHP MRI data and the need for robust multi-subject multi-center analyses. Streamlined acquisition and analysis protocols would accelerate and improve these efforts. However, consensus on minimal standards for data acquisition protocols and analysis pipelines for NHP imaging remains to be established, particularly for multi-center studies. Here, we draw parallels between NHP and human neuroimaging and provide minimal guidelines for harmonizing and standardizing data acquisition. We advocate robust translation of widely used open-access toolkits that are well established for analyzing human data. We also encourage the use of validated, automated pre-processing tools for analyzing NHP data sets. These guidelines aim to refine methodological and analytical strategies for small and large-scale NHP neuroimaging data. This will improve reproducibility of results, and accelerate the convergence between NHP and human neuroimaging strategies which will ultimately benefit fundamental and translational brain science.

The nonhuman primate neuroimaging and neuroanatomy project.

Multi-modal neuroimaging projects such as the Human Connectome Project (HCP) and UK Biobank are advancing our understanding of human brain architecture, function, connectivity, and their variability across individuals using high-quality non-invasive data from many subjects. Such efforts depend upon the accuracy of non-invasive brain imaging measures. However, 'ground truth' validation of connectivity using invasive tracers is not feasible in humans. Studies using nonhuman primates (NHPs) enable comparisons between invasive and non-invasive measures, including exploration of how "functional connectivity" from fMRI and "tractographic connectivity" from diffusion MRI compare with long-distance connections measured using tract tracing. Our NonHuman Primate Neuroimaging & Neuroanatomy Project (NHP_NNP) is an international effort (6 laboratories in 5 countries) to: (i) acquire and analyze high-quality multi-modal brain imaging data of macaque and marmoset monkeys using protocols and methods adapted from the HCP; (ii) acquire quantitative invasive tract-tracing data for cortical and subcortical projections to cortical areas; and (iii) map the distributions of different brain cell types with immunocytochemical stains to better define brain areal boundaries. We are acquiring high-resolution structural, functional, and diffusion MRI data together with behavioral measures from over 100 individual macaques and marmosets in order to generate non-invasive measures of brain architecture such as myelin and cortical thickness maps, as well as functional and diffusion tractography-based connectomes. We are using classical and next-generation anatomical tracers to generate quantitative connectivity maps based on brain-wide counting of labeled cortical and subcortical neurons, providing ground truth measures of connectivity. Advanced statistical modeling techniques address the consistency of both kinds of data across individuals, allowing comparison of tracer-based and non-invasive MRI-based connectivity measures. We aim to develop improved cortical and subcortical areal atlases by combining histological and imaging methods. Finally, we are collecting genetic and sociality-associated behavioral data in all animals in an effort to understand how genetic variation shapes the connectome and behavior.

Towards HCP-Style macaque connectomes: 24-Channel 3T multi-array coil, MRI sequences and preprocessing.

Macaque monkeys are an important animal model where invasive investigations can lead to a better understanding of the cortical organization of primates including humans. However, the tools and methods for noninvasive image acquisition (e.g. MRI RF coils and pulse sequence protocols) and image data preprocessing have lagged behind those developed for humans. To resolve the structural and functional characteristics of the smaller macaque brain, high spatial, temporal, and angular resolutions combined with high signal-to-noise ratio are required to ensure good image quality. To address these challenges, we developed a macaque 24-channel receive coil for 3-T MRI with parallel imaging capabilities. This coil enables adaptation of the Human Connectome Project (HCP) image acquisition protocols to the in-vivo macaque brain. In addition, we adapted HCP preprocessing methods to the macaque brain, including spatial minimal preprocessing of structural, functional MRI (fMRI), and diffusion MRI (dMRI). The coil provides the necessary high signal-to-noise ratio and high efficiency in data acquisition, allowing four- and five-fold accelerations for dMRI and fMRI. Automated FreeSurfer segmentation of cortex, reconstruction of cortical surface, removal of artefacts and nuisance signals in fMRI, and distortion correction of dMRI all performed well, and the overall quality of basic neurobiological measures was comparable with those for the HCP. Analyses of functional connectivity in fMRI revealed high sensitivity as compared with those from publicly shared datasets. Tractography-based connectivity estimates correlated with tracer connectivity similarly to that achieved using ex-vivo dMRI. The resulting HCP-style in vivo macaque MRI data show considerable promise for analyzing cortical architecture and functional and structural connectivity using advanced methods that have previously only been available in studies of the human brain.

Diffusion Tensor Model links to Neurite Orientation Dispersion and Density Imaging at high b-value in Cerebral Cortical Gray Matter.

Diffusion tensor imaging (DTI) and neurite orientation dispersion and density imaging (NODDI) are widely used models to infer microstructural features in the brain from diffusion-weighted MRI. Several studies have recently applied both models to increase sensitivity to biological changes, however, it remains uncertain how these measures are associated. Here we show that cortical distributions of DTI and NODDI are associated depending on the choice of b-value, a factor reflecting strength of diffusion weighting gradient. We analyzed a combination of high, intermediate and low b-value data of multi-shell diffusion-weighted MRI (dMRI) in healthy 456 subjects of the Human Connectome Project using NODDI, DTI and a mathematical conversion from DTI to NODDI. Cortical distributions of DTI and DTI-derived NODDI metrics were remarkably associated with those in NODDI, particularly when applied highly diffusion-weighted data (b-value = 3000 sec/mm2). This was supported by simulation analysis, which revealed that DTI-derived parameters with lower b-value datasets suffered from errors due to heterogeneity of cerebrospinal fluid fraction and partial volume. These findings suggest that high b-value DTI redundantly parallels with NODDI-based cortical neurite measures, but the conventional low b-value DTI is hard to reasonably characterize cortical microarchitecture.

Mitochondrial 2,4-dienoyl-CoA reductase (Decr) deficiency and impairment of thermogenesis in mouse brown adipose tissue.

A large number of studies have demonstrated significance of polyunsaturated fatty acids (PUFAs) for human health. However, many aspects on signals translating PUFA-sensing into body homeostasis have remained enigmatic. To shed light on PUFA physiology, we have generated a mouse line defective in mitochondrial dienoyl-CoA reductase (Decr), which is a key enzyme required for ß-oxidation of PUFAs. Previously, we have shown that these mice, whose oxidation of saturated fatty acid is intact but break-down of unsaturated fatty acids is blunted, develop severe hypoglycemia during metabolic stresses and fatal hypothermia upon acute cold challenge. In the current work, indirect calorimetry and thermography suggested that cold intolerance of Decr-/- mice is due to failure in maintaining appropriate heat production at least partly due to failure of brown adipose tissue (BAT) thermogenesis. Magnetic resonance imaging, electron microscopy, mass spectrometry and biochemical analysis showed attenuation in activation of lipolysis despite of functional NE-signaling and inappropriate expression of genes contributing to thermogenesis in iBAT when the Decr-/- mice were exposed to cold. We hypothesize that the failure in turning on BAT thermogenesis occurs due to accumulation of unsaturated long-chain fatty acids or their metabolites in Decr-/- mice BAT suppressing down-stream propagation of NE-signaling.

A novel Tungsten-based fiducial marker for multi-modal brain imaging.

BACKGROUND: Multi-modal brain image registration is a prerequisite for accurate mapping of brain structure and function in neuroscience. Image registration is commonly performed using automated software; however, its accuracy decreases when images differ in modality, contrast, uniformity, and resolution. This limitation could be overcome by using an external reference point; however, high-contrast agents in multi-modal imaging have not been previously reported. NEW METHODS: Here, we propose a novel multi-modal fiducial marker that contains Tungsten solution and provides high contrast in magnetic resonance imaging (MRI), computed tomography (CT), and positron emission tomography (PET). The basic characteristics of this multi-modal marker were investigated by assessing major sources of image contrast in the following modalities: density and T1-, T2-relaxivity in comparison with conventional contrast agents. RESULTS: Tungsten solution had lower T1- and T2-relaxivity and high solubility, and showed high contrast in T1- and T2-weighted MR and CT images at a high-density concentration (˜3.0 g/mL), whereas other conventional solutions did not show sufficient contrast in either CT or MRI. COMPARISON WITH EXISTING METHODS: The use of this Tungsten-based multi-modal marker allowed more accurate registration than a software-only method in phantom and animal experiments. Application of this method demonstrated accurate cortical surface mapping of neurotransmitter function (dopamine transporter, DAT) using PET and MRI, and provided a neurobiologically relevant cortical distribution consistent with previous literature on histology-based DAT immunoreactivity. CONCLUSIONS: The Tungsten-based multi-modal fiducial marker is non-radioactive, easy to handle, and aids precise registration across different modalities of brain imaging.

Neurite imaging reveals microstructural variations in human cerebral cortical gray matter.

We present distinct patterns of neurite distribution in the human cerebral cortex using diffusion magnetic resonance imaging (MRI). We analyzed both high-resolution structural (T1w and T2w images) and diffusion MRI data in 505 subjects from the Human Connectome Project. Neurite distributions were evaluated using the neurite orientation dispersion and density imaging (NODDI) model, optimized for gray matter, and mapped onto the cortical surface using a method weighted towards the cortical mid-thickness to reduce partial volume effects. The estimated neurite density was high in both somatosensory and motor areas, early visual and auditory areas, and middle temporal area (MT), showing a strikingly similar distribution to myelin maps estimated from the T1w/T2w ratio. The estimated neurite orientation dispersion was particularly high in early sensory areas, which are known for dense tangential fibers and are classified as granular cortex by classical anatomists. Spatial gradients of these cortical neurite properties revealed transitions that colocalize with some areal boundaries in a recent multi-modal parcellation of the human cerebral cortex, providing mutually supportive evidence. Our findings indicate that analyzing the cortical gray matter neurite morphology using diffusion MRI and NODDI provides valuable information regarding cortical microstructure that is related to but complementary to myeloarchitecture.

Parenchymal spin-lock fMRI signals associated with cortical spreading depression.

We found novel types of parenchymal functional magnetic resonance imaging (fMRI) signals in the rat brain during large increases in metabolism. Cortical spreading depression (CSD), a self-propagating wave of cellular activation, is associated with several pathologic conditions such as migraine and stroke. It was used as a paradigm to evoke transient neuronal depolarization leading to enhanced energy consumption. Activation of CSD was investigated using spin-lock (SL), diffusion, blood oxygenation level-dependent and cerebral blood volume fMRI techniques. Our results show that the SL-fMRI signal is generated by endogenous parenchymal mechanisms during CSD propagation, and these mechanisms are not associated with hemodynamic changes or cellular swelling. Protein phantoms suggest that pH change alone does not explain the observed SL-fMRI signal changes. However, increased amounts of inorganic phosphates released from high-energy phosphates combined with pH changes may produce SL- power-dependent longitudinal relaxation in the rotating frame (R1ρ) changes in protein phantoms that are similar to those observed during CSD, as seen before in acute ischemia under our experimental conditions. This links SL-fMRI changes intimately to energy metabolism and supports the use of the SL technique as a new, promising functional approach for noninvasive imaging of metabolic transitions in the active or pathologic brain.

High b-value diffusion-weighted fMRI in a rat forepaw electrostimulation model at 7 T.

Spin-echo diffusion-weighted functional MRI (DW-fMRI) was performed on a rat forepaw electrostimulation model at 7 T. This small animal model used electric (rather than visual) stimulation and allowed DW-fMRI experiments to be performed over a broader range of acquisition parameters than previous work on humans and cats. Resting state experiments with injections of ultra-small superparamagnetic iron oxide (USPIO) were also used to investigate the effects of gradient coupling on the signal change. The experiments were performed over five b-values (0, 200, 800, 1400 and 2000s/mm(2)) and three echo-times (30, 60 and 90 ms). Alterations to the stimulation-induced response with respect to TE and b-value were evaluated in two intervals: the positive stimulus-correlated response (5-20s after stimulus onset) and the post-stimulus undershoot (27-40s). There was no strong dependence of the signal change on b-value for any of the intervals or TEs. Similarly, changes to the apparent transverse relaxation rate showed no clear dependence on b-value. In contrast to previous DW-fMRI studies, the simplest explanation for the observed data is a single-compartment signal model with the functional signal changes probably corresponding to extravascular SE-BOLD. Experiments with USPIO suggested that at 7 T and within the range of parameters used, the influence of gradient coupling may be sufficient to explain minor DW-fMRI signal changes.