Dynamic and inherent B0 correction for DTI using stimulated echo spiral imaging

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.

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.

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.

Improved dynamic distortion correction for fMRI using single-echo EPI and a readout-reversed first image (REFILL)

The boundaries between tissues with different magnetic susceptibilities generate inhomogeneities in the main magnetic field which change over time due to motion, respiration and system instabilities. The dynamically changing field can be measured from the phase of the fMRI data and corrected. However, methods for doing so need multi-echo data, time-consuming reference scans and/or involve error-prone processing steps, such as phase unwrapping, which are difficult to implement robustly on the MRI host. The improved dynamic distortion correction method we propose is based on the phase of the single-echo EPI data acquired for fMRI, phase offsets calculated from a triple-echo, bipolar reference scan of circa 3-10 s duration using a method which avoids the need for phase unwrapping and an additional correction derived from one EPI volume in which the readout direction is reversed. This Reverse-Encoded First Image and Low resoLution reference scan (REFILL) approach is shown to accurately measure B0 as it changes due to shim, motion and respiration, even with large dynamic changes to the field at 7 T, where it led to a > 20% increase in time-series signal to noise ratio compared to data corrected with the classic static approach. fMRI results from REFILL-corrected data were free of stimulus-correlated distortion artefacts seen when data were corrected with static field mapping. The method is insensitive to shim changes and eddy current differences between the reference scan and the fMRI time series, and employs calculation steps that are simple and robust, allowing most data processing to be performed in real time on the scanner image reconstruction computer. These improvements make it feasible to routinely perform dynamic distortion correction in fMRI.

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.

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.

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 T₁-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.

Whole-Brain Imaging of Subvoxel T1-Diffusion Correlation Spectra in Human Subjects

T1 relaxation and water mobility generate eloquent MRI tissue contrasts with great diagnostic value in many neuroradiological applications. However, conventional methods do not adequately quantify the microscopic heterogeneity of these important biophysical properties within a voxel, and therefore have limited biological specificity. We describe a new correlation spectroscopic (CS) MRI method for measuring how T1 and mean diffusivity (MD) co-vary in microscopic tissue environments. We develop a clinical pulse sequence that combines inversion recovery (IR) with single-shot isotropic diffusion encoding (IDE) to efficiently acquire whole-brain MRIs with a wide range of joint T1-MD weightings. Unlike conventional diffusion encoding, the IDE preparation ensures that all subvoxel water pools are weighted by their MDs regardless of the sizes, shapes, and orientations of their corresponding microscopic diffusion tensors. Accordingly, IR-IDE measurements are well-suited for model-free, quantitative spectroscopic analysis of microscopic water pools. Using numerical simulations, phantom experiments, and data from healthy volunteers we demonstrate how IR-IDE MRIs can be processed to reconstruct maps of two-dimensional joint probability density functions, i.e., correlation spectra, of subvoxel T1-MD values. In vivo T1-MD spectra show distinct cerebrospinal fluid and parenchymal tissue components specific to white matter, cortical gray matter, basal ganglia, and myelinated fiber pathways, suggesting the potential for improved biological specificity. The one-dimensional marginal distributions derived from the T1-MD correlation spectra agree well with results from other relaxation spectroscopic and quantitative MRI studies, validating the T1-MD contrast encoding and the spectral reconstruction. Mapping subvoxel T1-diffusion correlations in patient populations may provide a more nuanced, comprehensive, sensitive, and specific neuroradiological assessment of the non-specific changes seen on fluid-attenuated inversion recovery (FLAIR) and diffusion-weighted MRIs (DWIs) in cancer, ischemic stroke, or brain injury.

Repeatability and reproducibility of prospective motion- and shim corrected 2D glycoCEST MRI

Background

Repeated glycoCEST MRI measurements on the same subject should produce similar results under the same environmental and experimental conditions. However, fluctuations in the static B₀ field, which may occur between and within measurements due to heating of the shim iron or subject motion, may alter results and affect reproducibility. Here we investigate the repeatability and reproducibility of glycoCEST measurements and examine the effectiveness of a real-time shim- and motion navigated chemical exchange saturation transfer (CEST) sequence to improve reproducibility.

Methods

In nine subjects, double volumetric navigated (DvNav)-CEST acquisitions in the calf muscle were repeated five times in each of two sessions-the first without correction, and the second with real-time shim- and motion correction applied. In both sessions a dynamically changing field was introduced by running a 5-minute gradient intensive diffusion sequence. We evaluated the effect of the introduced B₀ inhomogeneity on the reproducibility of glycoCEST, where the small chemical shift difference between the hydroxyl and bulk water protons at 3 T makes CEST quantification extremely sensitive to magnetic field inhomogeneities.

Results

With real-time shim- and motion correction, glycoCEST results were relatively consistent with mean coefficient of variation (CoV) 2.7%±1.4% across all subjects, whereas without correction the results were less consistent with CoV 84%±71%.

Conclusions

Our results demonstrate that real-time shim- and motion correction can mitigate effects of B0 field fluctuations and improve reproducibility of glycoCEST data. This is important when conducting longitudinal studies or when using glycoCEST MRI to assess treatment or physiological responses over time.

The quest for high spatial resolution diffusion-weighted imaging of the human brain in vivo

Diffusion-weighted imaging, a contrast unique to MRI, is used for assessment of tissue microstructure in vivo. However, this exquisite sensitivity to finer scales far above imaging resolution comes at the cost of vulnerability to errors caused by sources of motion other than diffusion motion. Addressing the issue of motion has traditionally limited diffusion-weighted imaging to a few acquisition techniques and, as a consequence, to poorer spatial resolution than other MRI applications. Advances in MRI imaging methodology have allowed diffusion-weighted MRI to push to ever higher spatial resolution. In this review we focus on the pulse sequences and associated techniques under development that have pushed the limits of image quality and spatial resolution in diffusion-weighted MRI.

Magnetic resonance microdynamic imaging reveals distinct tissue microenvironments

Magnetic resonance imaging (MRI) provides a powerful set of tools with which to investigate biological tissues noninvasively and in vivo. Tissues are heterogeneous in nature; an imaging voxel contains an ensemble of different cells and extracellular matrix components. A long-standing challenge has been to infer the content of and interactions among these microscopic tissue components within a macroscopic imaging voxel. Spatially resolved multidimensional relaxation-diffusion correlation (REDCO) spectroscopy holds the potential to deliver such microdynamic information. However, to date, vast data requirements have mostly relegated these type of measurements to nuclear magnetic resonance applications and prevented them from being widely and successfully used in conjunction with imaging. By using a novel data acquisition and processing strategy in this study, spatially resolved REDCO could be performed in reasonable scanning times with excellent prospects for clinical applications. This new MR imaging framework-which we term "magnetic resonance microdynamic imaging (MRMI)"-permits the simultaneous noninvasive and model-free quantification of multiple subcellular, cellular, and interstitial tissue microenvironments within a voxel. MRMI is demonstrated with a fixed spinal cord specimen, enabling the quantification of microscopic tissue components with unprecedented specificity. Tissue components, such as axons, neuronal and glial soma, and myelin were identified on the basis of their multispectral signature within individual imaging voxels. These tissue elements could then be composed into images and be correlated with immunohistochemistry findings. MRMI provides novel image contrasts of tissue components and a new family of microdynamic biomarkers that could lead to new diagnostic imaging approaches to probe biological tissue alterations accompanied by pathological or developmental changes.