An ex vivo imaging pipeline for producing high-quality and high-resolution diffusion-weighted imaging datasets

Beyond the surface: how ex-vivo diffusion-weighted imaging reveals large animal brain microstructure and connectivity

Diffusion-weighted Imaging (DWI) is an effective and state-of-the-art neuroimaging method that non-invasively reveals the microstructure and connectivity of tissues. Recently, novel applications of the DWI technique in studying large brains through ex-vivo imaging enabled researchers to gain insights into the complex neural architecture in different species such as those of Perissodactyla (e.g., horses and rhinos), Artiodactyla (e.g., bovids, swines, and cetaceans), and Carnivora (e.g., felids, canids, and pinnipeds). Classical in-vivo tract-tracing methods are usually considered unsuitable for ethical and practical reasons, in large animals or protected species. Ex-vivo DWI-based tractography offers the chance to examine the microstructure and connectivity of formalin-fixed tissues with scan times and precision that is not feasible in-vivo. This paper explores DWI's application to ex-vivo brains of large animals, highlighting the unique insights it offers into the structure of sometimes phylogenetically different neural networks, the connectivity of white matter tracts, and comparative evolutionary adaptations. Here, we also summarize the challenges, concerns, and perspectives of ex-vivo DWI that will shape the future of the field in large brains.

Exploring the Potential of Machine Learning Algorithms to Improve Diffusion Nuclear Magnetic Resonance Imaging Models Analysis

Purpose

This paper explores different machine learning (ML) algorithms for analyzing diffusion nuclear magnetic resonance imaging (dMRI) models when analytical fitting shows restrictions. It reviews various ML techniques for dMRI analysis and evaluates their performance on different b-values range datasets, comparing them with analytical methods.

Materials and Methods

After standard fitting for reference, four sets of diffusion-weighted nuclear magnetic resonance images were used to train/test various ML algorithms for prediction of diffusion coefficient (D), pseudo-diffusion coefficient (D*), perfusion fraction (f), and kurtosis (K). ML classification algorithms, including extra-tree classifier (ETC), logistic regression, C-support vector, extra-gradient boost, and multilayer perceptron (MLP), were used to determine the existence of diffusion parameters (D, D*, f, and K) within single voxels. Regression algorithms, including linear regression, polynomial regression, ridge, lasso, random forest (RF), elastic-net, and support-vector machines, were used to estimate the value of the diffusion parameters. Performance was evaluated using accuracy (ACC), area under the curve (AUC) tests, and cross-validation root mean square error (RMSECV). Computational timing was also assessed.

Results

ETC and MLP were the best classifiers, with 94.1% and 91.7%, respectively, for the ACC test and 98.7% and 96.3% for the AUC test. For parameter estimation, RF algorithm yielded the most accurate results The RMSECV percentages were: 8.39% for D, 3.57% for D*, 4.52% for f, and 3.53% for K. After the training phase, the ML methods demonstrated a substantial decrease in computational time, being approximately 232 times faster than the conventional methods.

Conclusions

The findings suggest that ML algorithms can enhance the efficiency of dMRI model analysis and offer new perspectives on the microstructural and functional organization of biological tissues. This paper also discusses the limitations and future directions of ML-based dMRI analysis.

Eddy current-induced artifact correction in high b-value ex vivo human brain diffusion MRI with dynamic field monitoring

PURPOSE

To investigate whether spatiotemporal magnetic field monitoring can correct pronounced eddy current-induced artifacts incurred by strong diffusion-sensitizing gradients up to 300 mT/m used in high b-value diffusion-weighted (DW) EPI.

METHODS

A dynamic field camera equipped with 16 1 H NMR field probes was first used to characterize field perturbations caused by residual eddy currents from diffusion gradients waveforms in a 3D multi-shot EPI sequence on a 3T Connectom scanner for different gradient strengths (up to 300 mT/m), diffusion directions, and shots. The efficacy of dynamic field monitoring-based image reconstruction was demonstrated on high-gradient strength, submillimeter resolution whole-brain ex vivo diffusion MRI. A 3D multi-shot image reconstruction framework was developed that incorporated the nonlinear phase evolution measured with the dynamic field camera.

RESULTS

Phase perturbations in the readout induced by residual eddy currents from strong diffusion gradients are highly nonlinear in space and time, vary among diffusion directions, and interfere significantly with the image encoding gradients, changing the k-space trajectory. During the readout, phase modulations between odd and even EPI echoes become non-static and diffusion encoding direction-dependent. Superior reduction of ghosting and geometric distortion was achieved with dynamic field monitoring compared to ghosting reduction approaches such as navigator- and structured low-rank-based methods or MUSE followed by image-based distortion correction with the FSL tool "eddy."

CONCLUSION

Strong eddy current artifacts characteristic of high-gradient strength DW-EPI can be well corrected with dynamic field monitoring-based image reconstruction.

Mapping the individual human cortex using multidimensional MRI and unsupervised learning

Human evolution has seen the development of higher-order cognitive and social capabilities in conjunction with the unique laminar cytoarchitecture of the human cortex. Moreover, early-life cortical maldevelopment has been associated with various neurodevelopmental diseases. Despite these connections, there is currently no noninvasive technique available for imaging the detailed cortical laminar structure. This study aims to address this scientific and clinical gap by introducing an approach for imaging human cortical lamina. This method combines diffusion-relaxation multidimensional MRI with a tailored unsupervised machine learning approach that introduces enhanced microstructural sensitivity. This new imaging method simultaneously encodes the microstructure, the local chemical composition and importantly their correlation within complex and heterogenous tissue. To validate our approach, we compared the intra-cortical layers obtained using our ex vivo MRI-based method with those derived from Nissl staining of postmortem human brain specimens. The integration of unsupervised learning with diffusion-relaxation correlation MRI generated maps that demonstrate sensitivity to areal differences in cytoarchitectonic features observed in histology. Significantly, our observations revealed layer-specific diffusion-relaxation signatures, showing reductions in both relaxation times and diffusivities at the deeper cortical levels. These findings suggest a radial decrease in myelin content and changes in cell size and anisotropy, reflecting variations in both cytoarchitecture and myeloarchitecture. Additionally, we demonstrated that 1D relaxation and high-order diffusion MRI scalar indices, even when aggregated and used jointly in a multimodal fashion, cannot disentangle the cortical layers. Looking ahead, our technique holds the potential to open new avenues of research in human neurodevelopment and the vast array of disorders caused by disruptions in neurodevelopment.

versaFlow: a versatile pipeline for resolution adapted diffusion MRI processing and its application to studying the variability of the PRIME-DE database

The lack of "gold standards" in Diffusion Weighted Imaging (DWI) makes validation cumbersome. To tackle this task, studies use translational analysis where results in humans are benchmarked against findings in other species. Non-Human Primates (NHP) are particularly interesting for this, as their cytoarchitecture is closely related to humans. However, tools used for processing and analysis must be adapted and finely tuned to work well on NHP images. Here, we propose versaFlow, a modular pipeline implemented in Nextflow, designed for robustness and scalability. The pipeline is tailored to in vivo NHP DWI at any spatial resolution; it allows for maintainability and customization. Processes and workflows are implemented using cutting-edge and state-of-the-art Magnetic Resonance Imaging (MRI) processing technologies and diffusion modeling algorithms, namely Diffusion Tensor Imaging (DTI), Constrained Spherical Deconvolution (CSD), and DIstribution of Anisotropic MicrOstructural eNvironments in Diffusion-compartment imaging (DIAMOND). Using versaFlow, we provide an in-depth study of the variability of diffusion metrics computed on 32 subjects from 3 sites of the Primate Data Exchange (PRIME-DE), which contains anatomical T1-weighted (T1w) and T2-weighted (T2w) images, functional MRI (fMRI), and DWI of NHP brains. This dataset includes images acquired over a range of resolutions, using single and multi-shell gradient samplings, on multiple scanner vendors. We perform a reproducibility study of the processing of versaFlow using the Aix-Marseilles site's data, to ensure that our implementation has minimal impact on the variability observed in subsequent analyses. We report very high reproducibility for the majority of metrics; only gamma distribution parameters of DIAMOND display less reproducible behaviors, due to the absence of a mechanism to enforce a random number seed in the software we used. This should be taken into consideration when future applications are performed. We show that the PRIME-DE diffusion data exhibits a great level of variability, similar or greater than results obtained in human studies. Its usage should be done carefully to prevent instilling uncertainty in statistical analyses. This hints at a need for sufficient harmonization in acquisition protocols and for the development of robust algorithms capable of managing the variability induced in imaging due to differences in scanner models and/or vendors.

CACTUS: a computational framework for generating realistic white matter microstructure substrates

Monte-Carlo diffusion simulations are a powerful tool for validating tissue microstructure models by generating synthetic diffusion-weighted magnetic resonance images (DW-MRI) in controlled environments. This is fundamental for understanding the link between micrometre-scale tissue properties and DW-MRI signals measured at the millimetre-scale, optimizing acquisition protocols to target microstructure properties of interest, and exploring the robustness and accuracy of estimation methods. However, accurate simulations require substrates that reflect the main microstructural features of the studied tissue. To address this challenge, we introduce a novel computational workflow, CACTUS (Computational Axonal Configurator for Tailored and Ultradense Substrates), for generating synthetic white matter substrates. Our approach allows constructing substrates with higher packing density than existing methods, up to 95% intra-axonal volume fraction, and larger voxel sizes of up to 500μm3 with rich fibre complexity. CACTUS generates bundles with angular dispersion, bundle crossings, and variations along the fibres of their inner and outer radii and g-ratio. We achieve this by introducing a novel global cost function and a fibre radial growth approach that allows substrates to match predefined targeted characteristics and mirror those reported in histological studies. CACTUS improves the development of complex synthetic substrates, paving the way for future applications in microstructure imaging.

An open resource combining multi-contrast MRI and microscopy in the macaque brain

Understanding brain structure and function often requires combining data across different modalities and scales to link microscale cellular structures to macroscale features of whole brain organisation. Here we introduce the BigMac dataset, a resource combining in vivo MRI, extensive postmortem MRI and multi-contrast microscopy for multimodal characterisation of a single whole macaque brain. The data spans modalities (MRI and microscopy), tissue states (in vivo and postmortem), and four orders of spatial magnitude, from microscopy images with micrometre or sub-micrometre resolution, to MRI signals on the order of millimetres. Crucially, the MRI and microscopy images are carefully co-registered together to facilitate quantitative multimodal analyses. Here we detail the acquisition, curation, and first release of the data, that together make BigMac a unique, openly-disseminated resource available to researchers worldwide. Further, we demonstrate example analyses and opportunities afforded by the data, including improvement of connectivity estimates from ultra-high angular resolution diffusion MRI, neuroanatomical insight provided by polarised light imaging and myelin-stained histology, and the joint analysis of MRI and microscopy data for reconstruction of the microscopy-inspired connectome. All data and code are made openly available.

Linear multi-scale modeling of diffusion MRI data: A framework for characterization of oriented structures across length scales

Diffusion-weighted magnetic resonance imaging (DW-MRI) has evolved to provide increasingly sophisticated investigations of the human brain's structural connectome in vivo. Restriction spectrum imaging (RSI) is a method that reconstructs the orientation distribution of diffusion within tissues over a range of length scales. In its original formulation, RSI represented the signal as consisting of a spectrum of Gaussian diffusion response functions. Recent technological advances have enabled the use of ultra-high b-values on human MRI scanners, providing higher sensitivity to intracellular water diffusion in the living human brain. To capture the complex diffusion time dependence of the signal within restricted water compartments, we expand upon the RSI approach to represent restricted water compartments with non-Gaussian response functions, in an extended analysis framework called linear multi-scale modeling (LMM). The LMM approach is designed to resolve length scale and orientation-specific information with greater specificity to tissue microstructure in the restricted and hindered compartments, while retaining the advantages of the RSI approach in its implementation as a linear inverse problem. Using multi-shell, multi-diffusion time DW-MRI data acquired with a state-of-the-art 3 T MRI scanner equipped with 300 mT/m gradients, we demonstrate the ability of the LMM approach to distinguish different anatomical structures in the human brain and the potential to advance mapping of the human connectome through joint estimation of the fiber orientation distributions and compartment size characteristics.

Tissue optimization strategies for high-quality ex vivo diffusion imaging

Ex vivo diffusion imaging can be used to study healthy and pathological tissue microstructure in the rodent brain with high resolution, providing a link between in vivo MRI and ex vivo microscopy techniques. Major challenges for the successful acquisition of ex vivo diffusion imaging data however are changes in the relaxivity and diffusivity of brain tissue following perfusion fixation. In this study we address this question by examining the combined effects of tissue preparation factors that influence signal-to-noise ratio (SNR) and consequently image quality, including fixative concentration, contrast agent concentration and tissue rehydration time. We present an optimization strategy combining these factors to manipulate theT 1 andT 2 of fixed tissue and maximize SNR efficiency. We apply this strategy in the rat brain, for a diffusion-weighted spin echo protocol with TE = 27 ms on a 9.4 T scanner with a 39 mm volume coil and 660 mT/m 114 mm gradient insert. We used a reduced fixative concentration of 2% paraformaldehyde (PFA), rehydration time more than 20 days, 15 mM Gd-DTPA in perfusate and TR 250 ms. This resulted in a doubling of SNR and an increase in SNR per unit time of 135% in cortical grey matter and 88% in white matter compared with 4% PFA and no contrast agent. This improved SNR efficiency enabled the acquisition of excellent-quality high-resolution (78 μ m isotropic voxel size) diffusion data with b = 4000 s/mm2 , 30 diffusion directions and a field of view of 40 × 13 × 18 mm3 in less than 4 days. It was also possible to achieve comparable data quality for a standard resolution (150 μ m) diffusion dataset in 2 1 4 h. In conclusion, the tissue optimization strategy presented here may be used to improve SNR, increase spatial resolution and/or allow faster acquisitions in preclinical ex vivo diffusion MRI experiments.

Multimodal 3D Mouse Brain Atlas Framework with the Skull-Derived Coordinate System

Magnetic resonance imaging (MRI) and light-sheet fluorescence microscopy (LSFM) are technologies that enable non-disruptive 3-dimensional imaging of whole mouse brains. A combination of complementary information from both modalities is desirable for studying neuroscience in general, disease progression and drug efficacy. Although both technologies rely on atlas mapping for quantitative analyses, the translation of LSFM recorded data to MRI templates has been complicated by the morphological changes inflicted by tissue clearing and the enormous size of the raw data sets. Consequently, there is an unmet need for tools that will facilitate fast and accurate translation of LSFM recorded brains to in vivo, non-distorted templates. In this study, we have developed a bidirectional multimodal atlas framework that includes brain templates based on both imaging modalities, region delineations from the Allen's Common Coordinate Framework, and a skull-derived stereotaxic coordinate system. The framework also provides algorithms for bidirectional transformation of results obtained using either MR or LSFM (iDISCO cleared) mouse brain imaging while the coordinate system enables users to easily assign in vivo coordinates across the different brain templates.

Eddy current-induced artifacts correction in high gradient strength diffusion MRI with dynamic field monitoring: demonstration in ex vivo human brain imaging

Purpose

To demonstrate the advantages of spatiotemporal magnetic field monitoring to correct eddy current-induced artifacts (ghosting and geometric distortions) in high gradient strength diffusion MRI (dMRI).

Methods

A dynamic field camera with 16 NMR field probes was used to characterize eddy current fields induced from diffusion gradients for different gradients strengths (up to 300 mT/m), diffusion directions, and shots in a 3D multi-shot EPI sequence on a 3T Connectom scanner. The efficacy of dynamic field monitoring-based image reconstruction was demonstrated on high-resolution whole brain ex vivo dMRI. A 3D multi-shot image reconstruction framework was informed with the actual nonlinear phase evolution measured with the dynamic field camera, thereby accounting for high-order eddy currents fields on top of the image encoding gradients in the image formation model.

Results

Eddy current fields from diffusion gradients at high gradient strength in a 3T Connectom scanner are highly nonlinear in space and time, inducing high-order spatial phase modulations between odd/even echoes and shots that are not static during the readout. Superior reduction of ghosting and geometric distortion was achieved with dynamic field monitoring compared to ghosting approaches such as navigator- and structured low-rank-based methods or MUSE, followed by image-based distortion correction with eddy. Improved dMRI analysis is demonstrated with diffusion tensor imaging and high-angular resolution diffusion imaging.

Conclusion

Strong eddy current artifacts characteristic of high gradient strength dMRI can be well corrected with dynamic field monitoring-based image reconstruction, unlike the two-step approach consisting of ghosting correction followed by geometric distortion reduction with eddy.

High-fidelity, high-spatial-resolution diffusion magnetic resonance imaging of ex vivo whole human brain at ultra-high gradient strength with structured low-rank echo-planar imaging ghost correction

Diffusion magnetic resonance imaging (dMRI) of whole ex vivo human brain specimens enables three-dimensional (3D) mapping of structural connectivity at the mesoscopic scale, providing detailed evaluation of fiber architecture and tissue microstructure at a spatial resolution that is difficult to access in vivo. To account for the short T2 and low diffusivity of fixed tissue, ex vivo dMRI is often acquired using strong diffusion-sensitizing gradients and multishot/segmented 3D echo-planar imaging (EPI) sequences to achieve high spatial resolution. However, the combination of strong diffusion-sensitizing gradients and multishot/segmented EPI readout can result in pronounced ghosting artifacts incurred by nonlinear spatiotemporal variations in the magnetic field produced by eddy currents. Such ghosting artifacts cannot be corrected with conventional correction solutions and pose a significant roadblock to leveraging human MRI scanners with ultrahigh gradients for ex vivo whole-brain dMRI. Here, we show that ghosting-correction approaches that correct for either polarity-related ghosting or shot-to-shot variations in a separate manner are suboptimal for 3D multishot diffusion-weighted EPI experiments in fixed human brain specimens using strong diffusion-sensitizing gradients on the 3-T Connectom MRI scanner, resulting in orientationally biased dMRI estimates. We apply a recently developed advanced k-space reconstruction method based on structured low-rank matrix (SLM) modeling that handles both polarity-related ghosting and shot-to-shot variation simultaneously, to mitigate artifacts in high-angular resolution multishot dMRI data acquired in several fixed human brain specimens at 0.7-0.8-mm isotropic spatial resolution using b-values up to 10,000 s/mm2 and gradient strengths up to 280 mT/m. We demonstrate the improved mapping of diffusion tensor imaging and fiber orientation distribution functions in key neuroanatomical areas distributed across the whole brain using SLM-based EPI ghost correction compared with alternative techniques.

Diffusion radiomics for subtyping and clustering in autism spectrum disorder: A preclinical study

Autism spectrum disorder (ASD) is a highly prevalent, heterogenous neurodevelopmental disorder. Neuroimaging methods such as functional, structural, and diffusion MRI have been used to identify candidate imaging biomarkers for ASD, but current findings remain non-specific and likely arise from the heterogeneity present in ASD. To account for this, efforts to subtype ASD have emerged as a potential strategy for both the study of ASD and advancement of tailored behavioral therapies and therapeutics. Towards these ends, to improve upon current neuroimaging methods, we propose combining biologically sensitive neurite orientation dispersion and density index (NODDI) diffusion MR imaging with radiomics image processing to create a new methodological approach that, we hypothesize, can sensitively and specifically capture neurobiology. We demonstrate this method can sensitively distinguish differences between four genetically distinct rat models of ASD (Fmr1, Pten, Nrxn1, Disc1). Further, we demonstrate diffusion radiomic analyses hold promise for subtyping in ASD as we show unsupervised clustering of NODDI radiomic data generates clusters specific to the underlying genetic differences between the animal models. Taken together, our findings suggest the unique application of radiomic analysis on NODDI diffusion MRI may have the capacity to sensitively and specifically disambiguate the neurobiological heterogeneity present in the ASD population.

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.

Magic angle effect on diffusion tensor imaging in ligament and brain

PURPOSE

To evaluate the magic angle effect on diffusion tensor imaging (DTI) measurements in rat ligaments and mouse brains.

METHODS

Three rat knee joints and three mouse brains were scanned at 9.4 T using a modified 3D diffusion-weighted spin echo pulse sequence with the isotropic spatial resolution of 45 μm. The b value was 1000 s/mm2 for rat knee and 4000 s/mm2 for mouse brain. DTI model was used to investigate the quantitative metrics at different orientations with respect to the main magnetic field. The collagen fiber structure of the ligament was validated with polarized light microscopy (PLM) imaging.

RESULTS

The signal intensity, signal-to-noise ratio (SNR), and DTI metrics in the ligament were strongly dependent on the collagen fiber orientation with respect to the main magnetic field from both simulation and actual MRI scans. The variation of fractional anisotropy (FA) was about ~32%, and the variation of mean diffusivity (MD) was ~11%. These findings were further validated with the numerical simulation at different SNRs (~10.0 to 86.0). Compared to the ligament, the DTI metrics showed little orientation dependence in mouse brains.

CONCLUSION

Magic angle effect plays an important role in DTI measurements in the highly ordered collagen-rich tissues, while MD showed less orientation dependence than FA.

Dose-response relationship between the variables of unilateral optogenetic stimulation and transcallosal evoked responses in rat motor cortex

Efficient interhemispheric integration of neural activity between left and right primary motor cortex (M1) is critical for inter-limb motor control. We employed optogenetic stimulation to establish a framework for probing transcallosal M1–M1 interactions in rats. We performed optogenetic stimulation of excitatory neurons in right M1 of male Sprague-Dawley rats. We recorded the transcallosal evoked potential in contralateral left M1 via chronically implanted electrodes. Recordings were performed under anesthesia combination of dexmedetomidine and a low concentration of isoflurane. We systematically varied the stimulation intensity and duration to characterize the relationship between stimulation parameters in right M1 and the characteristics of the evoked intracortical potentials in left M1. Optogenetic stimulation of right M1 consistently evoked a transcallosal response in left M1 with a consistent negative peak (N1) that sometimes was preceded by a smaller positive peak (P1). Higher stimulation intensity or longer stimulation duration gradually increased N1 amplitude and reduced N1 variability across trials. A combination of stimulation intensities of 5–10 mW with stimulus durations of 1–10 ms were generally sufficient to elicit a robust transcallosal response in most animal, with our optic fiber setup. Optogenetically stimulated excitatory neurons in M1 can reliably evoke a transcallosal response in anesthetized rats. Characterizing the relationship between “stimulation dose” and “response magnitude” (i.e., the gain function) of transcallosal M1-to-M1 excitatory connections can be used to optimize the variables of optogenetic stimulation and ensure stimulation efficacy.

Insights from the IronTract challenge: Optimal methods for mapping brain pathways from multi-shell diffusion MRI

Limitations in the accuracy of brain pathways reconstructed by diffusion MRI (dMRI) tractography have received considerable attention. While the technical advances spearheaded by the Human Connectome Project (HCP) led to significant improvements in dMRI data quality, it remains unclear how these data should be analyzed to maximize tractography accuracy. Over a period of two years, we have engaged the dMRI community in the IronTract Challenge, which aims to answer this question by leveraging a unique dataset. Macaque brains that have received both tracer injections and ex vivo dMRI at high spatial and angular resolution allow a comprehensive, quantitative assessment of tractography accuracy on state-of-the-art dMRI acquisition schemes. We find that, when analysis methods are carefully optimized, the HCP scheme can achieve similar accuracy as a more time-consuming, Cartesian-grid scheme. Importantly, we show that simple pre- and post-processing strategies can improve the accuracy and robustness of many tractography methods. Finally, we find that fiber configurations that go beyond crossing (e.g., fanning, branching) are the most challenging for tractography. The IronTract Challenge remains open and we hope that it can serve as a valuable validation tool for both users and developers of dMRI analysis methods.

Post mortem mapping of connectional anatomy for the validation of diffusion MRI

Diffusion MRI (dMRI) is a unique tool for the study of brain circuitry, as it allows us to image both the macroscopic trajectories and the microstructural properties of axon bundles in vivo. The Human Connectome Project ushered in an era of impressive advances in dMRI acquisition and analysis. As a result of these efforts, the quality of dMRI data that could be acquired in vivo improved substantially, and large collections of such data became widely available. Despite this progress, the main limitation of dMRI remains: it does not image axons directly, but only provides indirect measurements based on the diffusion of water molecules. Thus, it must be validated by methods that allow direct visualization of axons but that can only be performed in post mortem brain tissue. In this review, we discuss methods for validating the various features of connectional anatomy that are extracted from dMRI, both at the macro-scale (trajectories of axon bundles), and at micro-scale (axonal orientations and other microstructural properties). We present a range of validation tools, including anatomic tracer studies, Klingler's dissection, myelin stains, label-free optical imaging techniques, and others. We provide an overview of the basic principles of each technique, its limitations, and what it has taught us so far about the accuracy of different dMRI acquisition and analysis approaches.

Intra-Areal Visual Topography in Primate Brains Mapped with Probabilistic Tractography of Diffusion-Weighted Imaging

Noninvasive diffusion-weighted magnetic resonance imaging (dMRI) can be used to map the neural connectivity between distinct areas in the intact brain, but the standard resolution achieved fundamentally limits the sensitivity of such maps. We investigated the sensitivity and specificity of high-resolution postmortem dMRI and probabilistic tractography in rhesus macaque brains to produce retinotopic maps of the lateral geniculate nucleus (LGN) and extrastriate cortical visual area V5/MT based on their topographic connections with the previously established functional retinotopic map of primary visual cortex (V1). We also replicated the differential connectivity of magnocellular and parvocellular LGN compartments with V1 across visual field positions. Predicted topographic maps based on dMRI data largely matched the established retinotopy of both LGN and V5/MT. Furthermore, tractography based on in vivo dMRI data from the same macaque brains acquired at standard field strength (3T) yielded comparable topographic maps in many cases. We conclude that tractography based on dMRI is sensitive enough to reveal the intrinsic organization of ordered connections between topographically organized neural structures and their resultant functional organization.

Mapping brain-wide excitatory projectome of primate prefrontal cortex at submicron resolution and comparison with diffusion tractography

Resolving trajectories of axonal pathways in the primate prefrontal cortex remains crucial to gain insights into higher-order processes of cognition and emotion, which requires a comprehensive map of axonal projections linking demarcated subdivisions of prefrontal cortex and the rest of brain. Here, we report a mesoscale excitatory projectome issued from the ventrolateral prefrontal cortex (vlPFC) to the entire macaque brain by using viral-based genetic axonal tracing in tandem with high-throughput serial two-photon tomography, which demonstrated prominent monosynaptic projections to other prefrontal areas, temporal, limbic, and subcortical areas, relatively weak projections to parietal and insular regions but no projections directly to the occipital lobe. In a common 3D space, we quantitatively validated an atlas of diffusion tractography-derived vlPFC connections with correlative green fluorescent protein-labeled axonal tracing, and observed generally good agreement except a major difference in the posterior projections of inferior fronto-occipital fasciculus. These findings raise an intriguing question as to how neural information passes along long-range association fiber bundles in macaque brains, and call for the caution of using diffusion tractography to map the wiring diagram of brain circuits.

In the brain is a web of interconnected nerve cells that send messages to one another via spindly projections called axons. These axons join together at junctions called synapses to create circuits of nerve cells which connect neighboring or distant brain regions. Notably, long-range neural connections underpin higher-order cognitive skills (such as planning and emotion regulation) which make humans distinct from our primate relatives. Only by untangling these far-reaching networks can researchers begin to delineate what sets the human brain apart from other species.

Researchers deploy a range of imaging techniques to map neural networks: scanning entire brains using MRI machines, or imaging thin slices of fluorescently labelled brain tissue using powerful microscopes. However, tracing long-range axons at a high resolution is challenging, and has stirred up debate about whether some neural tracts, such as the inferior fronto-occipital fasciculus, are present in all primates or only humans.

To address these discrepancies, Yan, Yu et al. employed a two-pronged approach to map neural circuits in the brains of macaques. First, two techniques – called viral tracing and two-photon microscopy – were used to create a three-dimensional, fine-grain map showing how the ventrolateral prefrontal cortex (vlPFC), which regulates complex behaviors, connects to the rest of the brain. This revealed prominent axons from the vlPFC projecting via a single synapse to distant brain regions involved in higher-order functions, such as encoding memories and processing emotion. However, there were no direct, monosynaptic connections between the vlPFC and the occipital lobe, the brain’s visual processing center at the back of the head.

Next, Yan, Yu et al. used a specialized MRI scanner to create an atlas of neural circuits connected to the vlPFC, and compared these results to a technique tracing axons stained with a fluorescent dye. In general, there was good agreement between the two methods, except for major differences in the rear-end projections that typically form the inferior fronto-occipital fasciculus. This suggests that this long-range neural pathway exists in monkeys, but it connects via multiple synapses instead of a single junction as was previously thought.

The findings of Yan, Yu et al. provide new insights on the far-reaching neural pathways connecting distant parts of the macaque brain. It also suggests that atlases of neural circuits from whole brain scans should be taken with caution and validated using neural tracing experiments.

Resolution and b value dependent Structural Connectome in ex vivo Mouse Brain

Diffusion magnetic resonance imaging has been widely used in both clinical and preclinical studies to characterize tissue microstructure and structural connectivity. The diffusion MRI protocol for the Human Connectome Project (HCP) has been developed and optimized to obtain high-quality, high-resolution diffusion MRI (dMRI) datasets. However, such efforts have not been fully explored in preclinical studies, especially for rodents. In this study, high quality dMRI datasets of mouse brains were acquired at 9.4T system from two vendors. In particular, we acquired a high-spatial resolution dMRI dataset (25 μm isotropic with 126 diffusion encoding directions), which we believe to be the highest spatial resolution yet obtained; and a high-angular resolution dMRI dataset (50 μm isotropic with 384 diffusion encoding directions), which we believe to be the highest angular resolution compared to the dMRI datasets at the microscopic resolution. We systematically investigated the effects of three important parameters that affect the final outcome of the connectome: b value (1000s/mm² to 8000 s/mm²), angular resolution (10 to 126), and spatial resolution (25 μm to 200 μm). The stability of tractography and connectome increase with the angular resolution, where more than 50 angles is necessary to achieve consistent results. The connectome and quantitative parameters derived from graph theory exhibit a linear relationship to the b value (R² > 0.99); a single-shell acquisition with b value of 3000 s/mm² shows comparable results to the multi-shell high angular resolution dataset. The dice coefficient decreases and both false positive rate and false negative rate gradually increase with coarser spatial resolution. Our study provides guidelines and foundations for exploration of tradeoffs among acquisition parameters for the structural connectome in ex vivo mouse brain.

A novel male Japanese quail structural connectivity atlas using ultra-high field diffusion MRI at 11.7 T

The structural connectivity of animal brains can be revealed using post-mortem diffusion-weighted magnetic resonance imaging (MRI). Despite the existence of several structural atlases of avian brains, few of them address the bird’s structural connectivity. In this study, a novel atlas of the structural connectivity is proposed for the male Japanese quail (Coturnix japonica), aiming at investigating two lines divergent on their emotionality trait: the short tonic immobility (STI) and the long tonic immobility (LTI) lines. The STI line presents a low emotionality trait, while the LTI line expresses a high emotionality trait. 21 male Japanese quail brains from both lines were scanned post-mortem for this study, using a preclinical Bruker 11.7 T MRI scanner. Diffusion-weighted MRI was performed using a 3D segmented echo planar imaging (EPI) pulsed gradient spin-echo (PGSE) sequence with a 200 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\upmu$$\end{document}μm isotropic resolution, 75 diffusion-encoding directions and a b-value fixed at 4500 s/mm². Anatomical MRI was likewise performed using a 2D anatomical T₂-weighted spin-echo (SE) sequence with a 150 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\upmu$$\end{document}μm isotropic resolution. This very first anatomical connectivity atlas of the male Japanese quail reveals 34 labeled fiber tracts and the existence of structural differences between the connectivity patterns characterizing the two lines. Thus, the link between the male Japanese quail’s connectivity and its underlying anatomical structures has reached a better understanding.

The Digital Brain Bank, an open access platform for post-mortem imaging datasets

Post-mortem magnetic resonance imaging (MRI) provides the opportunity to acquire high-resolution datasets to investigate neuroanatomy and validate the origins of image contrast through microscopy comparisons. We introduce the Digital Brain Bank (open.win.ox.ac.uk/DigitalBrainBank), a data release platform providing open access to curated, multimodal post-mortem neuroimaging datasets. Datasets span three themes-Digital Neuroanatomist: datasets for detailed neuroanatomical investigations; Digital Brain Zoo: datasets for comparative neuroanatomy; and Digital Pathologist: datasets for neuropathology investigations. The first Digital Brain Bank data release includes 21 distinctive whole-brain diffusion MRI datasets for structural connectivity investigations, alongside microscopy and complementary MRI modalities. This includes one of the highest-resolution whole-brain human diffusion MRI datasets ever acquired, whole-brain diffusion MRI in fourteen nonhuman primate species, and one of the largest post-mortem whole-brain cohort imaging studies in neurodegeneration. The Digital Brain Bank is the culmination of our lab's investment into post-mortem MRI methodology and MRI-microscopy analysis techniques. This manuscript provides a detailed overview of our work with post-mortem imaging to date, including the development of diffusion MRI methods to image large post-mortem samples, including whole, human brains. Taken together, the Digital Brain Bank provides cross-scale, cross-species datasets facilitating the incorporation of post-mortem data into neuroimaging studies.

A time-course study of actively stained mouse brains: Diffusion tensor imaging parameters and connectomic stability over 1 year

While the application of diffusion tensor imaging (DTI), tractography, and connectomics to fixed tissue is a common practice today, there have been limited studies examining the effects of fixation on brain microstructure over extended periods. This mouse model time-course study reports the changes of regional brain volumes and diffusion scalar parameters, such as fractional anisotropy, across 12 representative brain regions as measures of brain structural stability. The scalar DTI parameters and regional volumes were highly variable over the first 2 weeks after fixation. The same parameters were consistent over a 2-8-week window after fixation, which means confounds from tissue stability over that scanning window were minimal. Quantitative connectomes were analyzed over the same time with extension out to 1 year. While there was some change in the scalar metrics at 1 year after fixation, these changes were sufficiently small, particularly in white matter, to support reproducible connectomes over a period ranging from 2-weeks to 1-year post-fixation. These findings delineate a scanning period, during which brain volumes, diffusion scalar metrics, and connectomes are remarkably consistent.

Axonal T2 estimation using the spherical variance of the strongly diffusion-weighted MRI signal

In magnetic resonance imaging, the application of a strong diffusion weighting suppresses the signal contributions from the less diffusion-restricted constituents of the brain's white matter, thus enabling the estimation of the transverse relaxation time T₂ that arises from the more diffusion-restricted constituents such as the axons. However, the presence of cell nuclei and vacuoles can confound the estimation of the axonal T₂, as diffusion within those structures is also restricted, causing the corresponding signal to survive the strong diffusion weighting. We devise an estimator of the axonal T₂ based on the directional spherical variance of the strongly diffusion-weighted signal. The spherical variance T₂ estimates are insensitive to the presence of isotropic contributions to the signal like those provided by cell nuclei and vacuoles. We show that with a strong diffusion weighting these estimates differ from those obtained using the directional spherical mean of the signal which contains both axonal and isotropically-restricted contributions. Our findings hint at the presence of an MRI-visible isotropically-restricted contribution to the signal in the white matter ex vivo fixed tissue (monkey) at 7T, and do not allow us to discard such a possibility also for in vivo human data collected with a clinical 3T system.

Data-driven separation of MRI signal components for tissue characterization

PURPOSE

MRI can be utilized for quantitative characterization of tissue. To assess e.g. water fractions or diffusion coefficients for compartments in the brain, a decomposition of the signal is necessary. Imposing standard models carries the risk of estimating biased parameters if model assumptions are violated. This work introduces a data-driven multicomponent analysis, the monotonous slope non-negative matrix factorization (msNMF), tailored to extract data features expected in MR signals.

METHODS

The msNMF was implemented by extending the standard NMF with monotonicity constraints on the signal profiles and their first derivatives. The method was validated using simulated data, and subsequently applied to both ex vivo DWI data and in vivo relaxometry data. Reproducibility of the method was tested using the latter.

RESULTS

The msNMF recovered the multi-exponential signals in the simulated data and showed superiority to standard NMF (based on the explained variance, area under the ROC curve, and coefficient of variation). Diffusion components extracted from the DWI data reflected the cell density of the underlying tissue. The relaxometry analysis resulted in estimates of edema water fractions (EWF) highly correlated with published results, and demonstrated acceptable reproducibility.

CONCLUSION

The msNMF can robustly separate MR signals into components with relation to the underlying tissue composition, and may potentially be useful for e.g. tumor tissue characterization.

Does powder averaging remove dispersion bias in diffusion MRI diameter estimates within real 3D axonal architectures?

Noninvasive estimation of axon diameter with diffusion MRI holds the potential to investigate the dynamic properties of the brain network and pathology of neurodegenerative diseases. Recent studies use powder averaging to account for complex white matter architectures, but these have not been validated for real axonal geometries from regions that contain fibre crossings. Here, we present 120-304μm long segmented axons from X-ray nano-holotomography volumes of a splenium and crossing fibre region of a vervet monkey brain. We show that the axons in the complex crossing fibre region, which contains callosal, association, and corticospinal connections, are larger and exhibit a wider distribution than those of the splenium region. To accurately estimate the axon diameter in these regions, therefore, sensitivity to a wide range of diameters is required. We demonstrate how the q-value, b-value, signal-to-noise ratio and the assumed intra-axonal parallel diffusivity influence the range of measurable diameters with powder average approaches. Furthermore, we show how Gaussian distributed noise results in a wider range of measurable diameter at high b-values than Rician distributed noise, even at high signal-to-noise ratios of 100. The number of gradient directions is also shown to impose a lower bound on measurable diameter. Our results indicate that axon diameter estimation can be performed with only few b-shells, and that additional shells do not improve the accuracy of the estimate. For strong gradients available on human Connectom and preclinical scanners, Monte Carlo simulations of diffusion confirm that powder averaging techniques succeed in providing accurate estimates of axon diameter across a range of sequence parameters and diffusion times, even in complex white matter architectures. At relatively low b-values, the diameter estimate becomes sensitive to axonal microdispersion and the intra-axonal parallel diffusivity shows time dependency at both in vivo and ex vivo intrinsic diffusivities.

Diffusion MRI and anatomic tracing in the same brain reveal common failure modes of tractography

Anatomic tracing is recognized as a critical source of knowledge on brain circuitry that can be used to assess the accuracy of diffusion MRI (dMRI) tractography. However, most prior studies that have performed such assessments have used dMRI and tracer data from different brains and/or have been limited in the scope of dMRI analysis methods allowed by the data. In this work, we perform a quantitative, voxel-wise comparison of dMRI tractography and anatomic tracing data in the same macaque brain. An ex vivo dMRI acquisition with high angular resolution and high maximum b-value allows us to compare a range of q-space sampling, orientation reconstruction, and tractography strategies. The availability of tracing in the same brain allows us to localize the sources of tractography errors and to identify axonal configurations that lead to such errors consistently, across dMRI acquisition and analysis strategies. We find that these common failure modes involve geometries such as branching or turning, which cannot be modeled well by crossing fibers. We also find that the default thresholds that are commonly used in tractography correspond to rather conservative, low-sensitivity operating points. While deterministic tractography tends to have higher sensitivity than probabilistic tractography in that very conservative threshold regime, the latter outperforms the former as the threshold is relaxed to avoid missing true anatomical connections. On the other hand, the q-space sampling scheme and maximum b-value have less of an impact on accuracy. Finally, using scans from a set of additional macaque brains, we show that there is enough inter-individual variability to warrant caution when dMRI and tracer data come from different animals, as is often the case in the tractography validation literature. Taken together, our results provide insights on the limitations of current tractography methods and on the critical role that anatomic tracing can play in identifying potential avenues for improvement.

A high-resolution interactive atlas of the human brainstem using magnetic resonance imaging

Conventional atlases of the human brainstem are limited by the inflexible, sparsely-sampled, two-dimensional nature of histology, or the low spatial resolution of conventional magnetic resonance imaging (MRI). Postmortem high-resolution MRI circumvents the challenges associated with both modalities. A single human brainstem specimen extending from the rostral diencephalon through the caudal medulla was prepared for imaging after the brain was removed from a 65-year-old male within 24 h of death. The specimen was formalin-fixed for two weeks, then rehydrated and placed in a custom-made MRI compatible tube and immersed in liquid fluorocarbon. MRI was performed in a 7-Tesla scanner with 120 unique diffusion directions. Acquisition time for anatomic and diffusion images were 14 h and 208 h, respectively. Segmentation was performed manually. Deterministic fiber tractography was done using strategically chosen regions of interest and avoidance, with manual editing using expert knowledge of human neuroanatomy. Anatomic and diffusion images were rendered with isotropic resolutions of 50 μm and 200 μm, respectively. Ninety different structures were segmented and labeled, and 11 different fiber bundles were rendered with tractography. The complete atlas is available online for interactive use at https://www.civmvoxport.vm.duke.edu/voxbase/login.php?return_url=%2Fvoxbase%2F. This atlas presents multiple contrasting datasets and selected tract reconstruction with unprecedented resolution for MR imaging of the human brainstem. There are immediate applications in neuroanatomical education, with the potential to serve future applications for neuroanatomical research and enhanced neurosurgical planning through "safe" zones of entry into the human brainstem.

White matter changes in the trigeminal spinal tract in chronic migraineurs: an ex vivo study combining ultra-high field diffusion tensor imaging and polarized light imaging microscopy

ABSTRACT

Chronic migraine (CM) is a disabling neurologic disorder that affects approximately 2% of the general population. Neuroimaging studies show functional involvement of trigeminal structures, such as the trigeminal spinal nucleus (Sp5) in migraine. However, structural changes in the Sp5 and the afferent trigeminal spinal tract (sp5) have never been found. The aim of this study was to test the hypothesis that white matter changes in the sp5 are a key feature of brain alterations in CM patients. We used diffusion Magnetic Resonance Imaging (dMRI) and polarized light imaging (PLI) of post mortem brainstem specimens from healthy controls (n = 5) and CM patients (n = 5) to study white matter alterations in the sp5. Within the sp5, dMRI metrics included fractional anisotropy (FA), mean diffusivity (MD), axial diffusivity (AD) and radial diffusivity (RD) values. PLI was used to assess myelin density by measure of the retardance values in the sp5. The present study provides histological evidence that structural alterations occur in the sp5 in CM patients as compared to healthy controls. Myelin-density, as assessed by retardance values, showed to be higher and a corresponding increase in FA-values was observed. In addition, accompanying decreases in MD-, AD- and RD-values were observed. This study shows that the sp5 undergoes neuroplastic changes, a feature which substantiates evidence for the hyperactivity of the Sp5 in migraine patients. More insights are needed to observe whether these changes only occur in CM patients.

The Complex Hodological Architecture of the Macaque Dorsal Intraparietal Areas as Emerging from Neural Tracers and DW-MRI Tractography

In macaque monkeys, dorsal intraparietal areas are involved in several daily visuomotor actions. However, their border and sources of cortical afferents remain loosely defined. Combining retrograde histologic tracing and MRI diffusion-based tractography, we found a complex hodology of the dorsal bank of the intraparietal sulcus (db-IPS), which can be subdivided into a rostral intraparietal area PEip, projecting to the spinal cord, and a caudal medial intraparietal area MIP lacking such projections. Both include an anterior and a posterior sector, emerging from their ipsilateral, gradient-like connectivity profiles. As tractography estimations, we used the cross-sectional area of the white matter bundles connecting each area with other parietal and frontal regions, after selecting regions of interest (ROIs) corresponding to the injection sites of neural tracers. For most connections, we found a significant correlation between the proportions of cells projecting to all sectors of PEip and MIP along the continuum of the db-IPS and tractography. The latter also revealed “false positive” but plausible connections awaiting histologic validation.

Investigation of the magnetic susceptibility properties of fresh and fixed mouse heart, liver, skeletal muscle and brain tissue

PURPOSE

Several magnetic resonance imaging (MRI) techniques exploit the difference in magnetic susceptibilities between tissues, but systematic measurements of tissue susceptibility are lacking. Furthermore, there is the question as to whether chemical fixation that is used for ex vivo MRI studies, affects the magnetic properties of the tissue. Here, we determined the magnetic susceptibility and water content of fresh and chemically fixed mouse tissue.

METHODS

Mass susceptibility of brain, heart, liver and skeletal muscle samples were determined on a vibrating sample magnetometer at room temperature. Measurements at 50, 125, 200 and 295 K were performed to assess the temperature dependence of susceptibility. Moreover, we measured water content of fresh and fixed samples.

RESULTS

All samples show mass susceptibilities between -0.068 and -1.929 × 10^-8 m³/kg, compared to -9.338 × 10^-9 m³/kg of double distilled water. Heart tissue has a more diamagnetic susceptibility than the other tissues. Compared to fresh tissue, fixed tissue has a less diamagnetic susceptibility. Fixed tissue was not different in water content to fresh tissue and showed no consistent dependence of susceptibility with temperature, whereas fresh tissue shows a decrease to at least 125 K, indicative of a paramagnetic component.

CONCLUSIONS

Biological tissues are diamagnetic in comparison to water, where the heart is more diamagnetic than the other tissues, with paramagnetic contributions. Fixation rendered tissue less diamagnetic compared to fresh tissue. Our measurements revealed differences in tissue susceptibility between VSM and QSM, inviting more research to compare susceptibility-based MRI methods with physical measurements of tissue susceptibility.

Evaluation of early microstructural changes in the R6/1 mouse model of Huntington's disease by ultra-high field diffusion MR imaging

Diffusion MRI (dMRI) has been able to detect early structural changes related to neurological symptoms present in Huntington's disease (HD). However, there is still a knowledge gap to interpret the biological significance at early neuropathological stages. The purpose of this study is two-fold: (i) establish if the combination of Ultra-High Field Diffusion MRI (UHFD-MRI) techniques can add a more comprehensive analysis of the early microstructural changes observed in HD, and (ii) evaluate if early changes in dMRI microstructural parameters can be linked to cellular biomarkers of neuroinflammation. Ultra-high field magnet (16.7T), diffusion tensor imaging (DTI), and neurite orientation dispersion and density imaging (NODDI) techniques were applied to fixed ex-vivo brains of a preclinical model of HD (R6/1 mice). Fractional anisotropy (FA) was decreased in deep and superficial grey matter (GM) as well as white matter (WM) brain regions with well-known early HD microstructure and connectivity pathology. NODDI parameters associated with the intracellular and extracellular compartment, such as intracellular ventricular fraction (ICVF), orientation dispersion index (ODI), and isotropic volume fractions (IsoVF) were altered in R6/1 mice GM. Further, histological studies in these areas showed that glia cell markers associated with neuroinflammation (GFAP & Iba1) were consistent with the dMRI findings. dMRI can be used to extract non-invasive information of neuropathological events present in the early stages of HD. The combination of multiple imaging techniques represents a better approach to understand the neuropathological process allowing the early diagnosis and neuromonitoring of patients affected by HD.

In vivo human whole-brain Connectom diffusion MRI dataset at 760 µm isotropic resolution

We present a whole-brain in vivo diffusion MRI (dMRI) dataset acquired at 760 μm isotropic resolution and sampled at 1260 q-space points across 9 two-hour sessions on a single healthy participant. The creation of this benchmark dataset is possible through the synergistic use of advanced acquisition hardware and software including the high-gradient-strength Connectom scanner, a custom-built 64-channel phased-array coil, a personalized motion-robust head stabilizer, a recently developed SNR-efficient dMRI acquisition method, and parallel imaging reconstruction with advanced ghost reduction algorithm. With its unprecedented resolution, SNR and image quality, we envision that this dataset will have a broad range of investigational, educational, and clinical applications that will advance the understanding of human brain structures and connectivity. This comprehensive dataset can also be used as a test bed for new modeling, sub-sampling strategies, denoising and processing algorithms, potentially providing a common testing platform for further development of in vivo high resolution dMRI techniques. Whole brain anatomical T1-weighted and T2-weighted images at submillimeter scale along with field maps are also made available.

Dedicated container for postmortem human brain ultra-high field magnetic resonance imaging

BACKGROUND

The emerging field of ultra-high field MRI (UHF-MRI, 7 Tesla and higher) provides the opportunity to image human brains at a higher resolution and with higher signal-to-noise ratios compared to the more widely available 1.5 and 3T scanners. Scanning postmortem tissue additionally allows for greatly increased scan times and fewer movement issues leading to improvements in image quality. However, typical postmortem neuroimaging routines involve placing the tissue within plastic bags that leave room for susceptibility artifacts from tissue-air interfaces, inadequate submersion, and leakage issues. To address these challenges in postmortem imaging, a custom-built nonferromagnetic container was developed that allows whole brain hemispheres to be scanned at sub-millimeter resolution within typical head-coils.

METHOD

The custom-built polymethylmethacrylaat container consists of a cylinder with a hemispheric side and a lid with valves on the adjacent side. This shape fits within common MR head-coils and allows whole hemispheres to be submerged and vacuum sealed within it reducing imaging artifacts that would otherwise arise at air-tissue boundaries. Two hemisphere samples were scanned on a Siemens 9.4T Magnetom MRI scanner. High resolution T2* weighted data was obtained with a custom 3D gradient echo (GRE) sequence and diffusion-weighted imaging (DWI) scans were obtained with a 3D kT-dSTEAM sequence along 48 directions.

RESULTS

The custom-built container proved to submerge and contain tissue samples effectively and showed no interferences with MR scanning acquisition. The 3D GRE sequence provided high resolution isotropic T2* weighted data at 250 μm which showed a clear visualization of gray and white matter structures. DWI scans allowed for dense reconstruction of structural white matter connections via tractography.

CONCLUSION

Using this custom-built container worked towards achieving high quality MR images of postmortem brain material. This procedure can have advantages over traditional schemes including utilization of a standardized protocol and the reduced likelihood of leakage. This methodology could be adjusted and used to improve typical postmortem imaging routines.

Comparison of diffusion MRI and CLARITY fiber orientation estimates in both gray and white matter regions of human and primate brain

Diffusion MRI (dMRI) represents one of the few methods for mapping brain fiber orientations non-invasively. Unfortunately, dMRI fiber mapping is an indirect method that relies on inference from measured diffusion patterns. Comparing dMRI results with other modalities is a way to improve the interpretation of dMRI data and help advance dMRI technologies. Here, we present methods for comparing dMRI fiber orientation estimates with optical imaging of fluorescently labeled neurofilaments and vasculature in 3D human and primate brain tissue cuboids cleared using CLARITY. The recent advancements in tissue clearing provide a new opportunity to histologically map fibers projecting in 3D, which represents a captivating complement to dMRI measurements. In this work, we demonstrate the capability to directly compare dMRI and CLARITY in the same human brain tissue and assess multiple approaches for extracting fiber orientation estimates from CLARITY data. We estimate the three-dimensional neuronal fiber and vasculature orientations from neurofilament and vasculature stained CLARITY images by calculating the tertiary eigenvector of structure tensors. We then extend CLARITY orientation estimates to an orientation distribution function (ODF) formalism by summing multiple sub-voxel structure tensor orientation estimates. In a sample containing part of the human thalamus, there is a mean angular difference of 19o±15o between the primary eigenvectors of the dMRI tensors and the tertiary eigenvectors from the CLARITY neurofilament stain. We also demonstrate evidence that vascular compartments do not affect the dMRI orientation estimates by showing an apparent lack of correspondence (mean angular difference = 49o±23o) between the orientation of the dMRI tensors and the structure tensors in the vasculature stained CLARITY images. In a macaque brain dataset, we examine how the CLARITY feature extraction depends on the chosen feature extraction parameters. By varying the volume of tissue over which the structure tensor estimates are derived, we show that orientation estimates are noisier with more spurious ODF peaks for sub-voxels below 30 µm3 and that, for our data, the optimal gray matter sub-voxel size is between 62.5 µm3 and 125 µm3. The example experiments presented here represent an important advancement towards robust multi-modal MRI-CLARITY comparisons.

Original experimental rat model of blast-induced mild traumatic brain injury: a pilot study

OBJECTIVE

Diagnosing blast-induced mild traumatic brain injury (mTBI) is difficult due to minimal imaging findings. This study aimed to establish a rat model of behavioral abnormality caused by blast-induced mTBI and detect new findings for therapeutic intervention.

METHODS

We used a bench-top blast wave generator with the blast wave exiting through a 20-mm I.D. nozzle aimed at the focused target. The blast wave was directed at the head of male Wistar rats under general anesthesia positioned prone 2.5 cm below the nozzle. Peak shock wave pressure was 646.2 ± 70.3 kPa.

RESULTS

After blast injury, mTBI rats did not show the findings of brain hemorrhage or contusion macroscopically and on hematoxylin-eosin-stained frozen sections but did show anorexia and weight loss in the early post-injury phase. Behavioral experiments revealed short-term memory impairment at 2 weeks and depression-like behavior at 2 and 6 weeks. Diffusion-weighted ex vivo MRI showed high-intensity areas in layers of the bilateral hippocampus. Immunohistochemical analysis revealed accumulation of reactive microglia and GFAP-positive astrocytes in the same region and loss of NeuN-positive neurons in the hippocampal pyramidal cell layer.

CONCLUSIONS

This model can reflect the pathophysiology of blast-induced mTBI and could potentially be used to develop therapeutic interventions in the future.

Towards in vivo g-ratio mapping using MRI: Unifying myelin and diffusion imaging

BACKGROUND

The g-ratio, quantifying the comparative thickness of the myelin sheath encasing an axon, is a geometrical invariant that has high functional relevance because of its importance in determining neuronal conduction velocity. Advances in MRI data acquisition and signal modelling have put in vivo mapping of the g-ratio, across the entire white matter, within our reach. This capacity would greatly increase our knowledge of the nervous system: how it functions, and how it is impacted by disease.

NEW METHOD

This is the second review on the topic of g-ratio mapping using MRI.

RESULTS

This review summarizes the most recent developments in the field, while also providing methodological background pertinent to aggregate g-ratio weighted mapping, and discussing pitfalls associated with these approaches.

COMPARISON WITH EXISTING METHODS

Using simulations based on recently published data, this review reveals caveats to the state-of-the-art calibration methods that have been used for in vivo g-ratio mapping. It highlights the need to estimate both the slope and offset of the relationship between these MRI-based markers and the true myelin volume fraction if we are really to achieve the goal of precise, high sensitivity g-ratio mapping in vivo. Other challenges discussed in this review further evidence the need for gold standard measurements of human brain tissue from ex vivo histology.

CONCLUSIONS

We conclude that the quest to find the most appropriate MRI biomarkers to enable in vivo g-ratio mapping is ongoing, with the full potential of many novel techniques yet to be investigated.

Diffusion MRI revealed altered inter-hippocampal projections in the mouse brain after intrauterine inflammation

Diffusion MRI (dMRI) is commonly used to map large axonal pathways in the white matter. Recent technical advances have also enabled dMRI to resolve the small and complex axonal and dendritic projections in the gray matter. This study investigated whether high-resolution dMRI can resolve the hippocampal neuronal projections and detect abnormal connections due to neurological injury. We performed 3D high spatial and angular resolution dMRI of the mouse brains of the offspring survivors from a model of intrauterine (UI) inflammation, who had known functional deficiency in the hippocampus. We used a novel hippocampal connection mapping method to quantify the intra- and inter-hippocampal projections among 34 automatically segmented hippocampal sub-regions. The results demonstrated wide-spread intra-hippocampal projections, but rather specific intra-hippocampal projections that primarily connected through the CA3 region. Compared with the control group (n = 9), UI-injured mice (n = 11) exhibited significantly reduced inter-hippocampal projection strength (p < 0.01), which correlated well with the neurobehavioral assessments (R² = 0.47). Furthermore, using a whole-brain fixel-based analysis, we identified reduced fiber-density in the CA3 and the ventral hippocampal commissure of the UI-injured mice, which may explain the reduced inter-hippocampal projections. Histological findings also indicated reduced commissural fibers due to the UI-injury. Our study suggested that the dMRI-based connectivity mapping technique can potentially characterize abnormal hippocampal projections in neurological disorders.

The sensitivity of diffusion MRI to microstructural properties and experimental factors

Diffusion MRI is a non-invasive technique to study brain microstructure. Differences in the microstructural properties of tissue, including size and anisotropy, can be represented in the signal if the appropriate method of acquisition is used. However, to depict the underlying properties, special care must be taken when designing the acquisition protocol as any changes in the procedure might impact on quantitative measurements. This work reviews state-of-the-art methods for studying brain microstructure using diffusion MRI and their sensitivity to microstructural differences and various experimental factors. Microstructural properties of the tissue at a micrometer scale can be linked to the diffusion signal at a millimeter-scale using modeling. In this paper, we first give an introduction to diffusion MRI and different encoding schemes. Then, signal representation-based methods and multi-compartment models are explained briefly. The sensitivity of the diffusion MRI signal to the microstructural components and the effects of curvedness of axonal trajectories on the diffusion signal are reviewed. Factors that impact on the quality (accuracy and precision) of derived metrics are then reviewed, including the impact of random noise, and variations in the acquisition parameters (i.e., number of sampled signals, b-value and number of acquisition shells). Finally, yet importantly, typical approaches to deal with experimental factors are depicted, including unbiased measures and harmonization. We conclude the review with some future directions and recommendations on this topic.

A comprehensive atlas of white matter tracts in the chimpanzee

Chimpanzees (Pan troglodytes) are, along with bonobos, humans’ closest living relatives. The advent of diffusion MRI tractography in recent years has allowed a resurgence of comparative neuroanatomical studies in humans and other primate species. Here we offer, in comparative perspective, the first chimpanzee white matter atlas, constructed from in vivo chimpanzee diffusion-weighted scans. Comparative white matter atlases provide a useful tool for identifying neuroanatomical differences and similarities between humans and other primate species. Until now, comprehensive fascicular atlases have been created for humans (Homo sapiens), rhesus macaques (Macaca mulatta), and several other nonhuman primate species, but never in a nonhuman ape. Information on chimpanzee neuroanatomy is essential for understanding the anatomical specializations of white matter organization that are unique to the human lineage.

Diffusion MRI tractography reveals the first complete atlas of white matter of the chimpanzee, with the potential to help understand differences between the organization of human and chimpanzee brains.

Axon morphology is modulated by the local environment and impacts the noninvasive investigation of its structure-function relationship

Axonal conduction velocity, which ensures efficient function of the brain network, is related to axon diameter. Noninvasive, in vivo axon diameter estimates can be made with diffusion magnetic resonance imaging, but the technique requires three-dimensional (3D) validation. Here, high-resolution, 3D synchrotron X-ray nano-holotomography images of white matter samples from the corpus callosum of a monkey brain reveal that blood vessels, cells, and vacuoles affect axonal diameter and trajectory. Within single axons, we find that the variation in diameter and conduction velocity correlates with the mean diameter, contesting the value of precise diameter determination in larger axons. These complex 3D axon morphologies drive previously reported 2D trends in axon diameter and g-ratio. Furthermore, we find that these morphologies bias the estimates of axon diameter with diffusion magnetic resonance imaging and, ultimately, impact the investigation and formulation of the axon structure-function relationship.

Axon diameter index estimation independent of fiber orientation distribution using high-gradient diffusion MRI

Axon diameter mapping using high-gradient diffusion MRI has generated great interest as a noninvasive tool for studying trends in axonal size in the human brain. One of the main barriers to mapping axon diameter across the whole brain is accounting for complex white matter fiber configurations (e.g., crossings and fanning), which are prevalent throughout the brain. Here, we present a framework for generalizing axon diameter index estimation to the whole brain independent of the underlying fiber orientation distribution using the spherical mean technique (SMT). This approach is shown to significantly benefit from the use of real-valued diffusion data with Gaussian noise, which reduces the systematic bias in the estimated parameters resulting from the elevation of the noise floor when using magnitude data with Rician noise. We demonstrate the feasibility of obtaining whole-brain orientationally invariant estimates of axon diameter index and relative volume fractions in six healthy human volunteers using real-valued diffusion data acquired on a dedicated high-gradient 3-Tesla human MRI scanner with 300 mT/m maximum gradient strength. The trends in axon diameter index are consistent with known variations in axon diameter from histology and demonstrate the potential of this generalized framework for revealing coherent patterns in axonal structure throughout the living human brain. The use of real-valued diffusion data provides a viable solution for eliminating the Rician noise floor and should be considered for all spherical mean approaches to microstructural parameter estimation.

Reproducing macaque lateral grasping and oculomotor networks using resting state functional connectivity and diffusion tractography

Cortico-cortical networks involved in motor control have been well defined in the macaque using a range of invasive techniques. The advent of neuroimaging has enabled non-invasive study of these large-scale functionally specialized networks in the human brain; however, assessing its accuracy in reproducing genuine anatomy is more challenging. We set out to assess the similarities and differences between connections of macaque motor control networks defined using axonal tracing and those reproduced using structural and functional connectivity techniques. We processed a cohort of macaques scanned in vivo that were made available by the open access PRIME-DE resource, to evaluate connectivity using diffusion imaging tractography and resting state functional connectivity (rs-FC). Sectors of the lateral grasping and exploratory oculomotor networks were defined anatomically on structural images, and connections were reproduced using different structural and functional approaches (probabilistic and deterministic whole-brain and seed-based tractography; group template and native space functional connectivity analysis). The results showed that parieto-frontal connections were best reproduced using both structural and functional connectivity techniques. Tractography showed lower sensitivity but better specificity in reproducing connections identified by tracer data. Functional connectivity analysis performed in native space had higher sensitivity but lower specificity and was better at identifying connections between intrasulcal ROIs than group-level analysis. Connections of AIP were most consistently reproduced, although those connected with prefrontal sectors were not identified. We finally compared diffusion MR modelling with histology based on an injection in AIP and speculate on anatomical bases for the observed false negatives. Our results highlight the utility of precise ex vivo techniques to support the accuracy of neuroimaging in reproducing connections, which is relevant also for human studies.

Mesoscale diffusion magnetic resonance imaging of the ex vivo human hippocampus

Mesoscale diffusion magnetic resonance imaging (MRI) endeavors to bridge the gap between macroscopic white matter tractography and microscopic studies investigating the cytoarchitecture of human brain tissue. To ensure a robust measurement of diffusion at the mesoscale, acquisition parameters were arrayed to investigate their effects on scalar indices (mean, radial, axial diffusivity, and fractional anisotropy) and streamlines (i.e., graphical representation of axonal tracts) in hippocampal layers. A mesoscale resolution afforded segementation of the pyramidal cell layer (CA1-4), the dentate gyrus, as well as stratum moleculare, radiatum, and oriens. Using ex vivo samples, surgically excised from patients with intractable epilepsy (n = 3), we found that shorter diffusion times (23.7 ms) with a b-value of 4,000 s/mm2 were advantageous at the mesoscale, providing a compromise between mean diffusivity and fractional anisotropy measurements. Spatial resolution and sample orientation exerted a major effect on tractography, whereas the number of diffusion gradient encoding directions minimally affected scalar indices and streamline density. A sample temperature of 15°C provided a compromise between increasing signal-to-noise ratio and increasing the diffusion properties of the tissue. Optimization of the acquisition afforded a system's view of intra- and extra-hippocampal connections. Tractography reflected histological boundaries of hippocampal layers. Individual layer connectivity was visualized, as well as streamlines emanating from individual sub-fields. The perforant path, subiculum and angular bundle demonstrated extra-hippocampal connections. Histology of the samples confirmed individual cell layers corresponding to ROIs defined on MR images. We anticipate that this ex vivo mesoscale imaging will yield novel insights into human hippocampal connectivity.

Diffusion Tensor Imaging Tractography of White Matter Tracts in the Equine Brain

Tractography, a noninvasive technique tracing brain pathways from diffusion tensor magnetic resonance imaging (DTI) data, is increasingly being used for brain investigation of domestic mammals. In the equine species, such a technique could be useful to improve our knowledge about structural connectivity or to assess structural changes of white matter tracts potentially associated with neurodegenerative diseases. The goals of the present study were to establish the feasibility of DTI tractography in the equine brain and to provide a morphologic description of the most representative tracts in this species. Postmortem DTI and susceptibility-weighted imaging (SWI) of an equine brain were acquired with a 3-T system using a head coil. Association, commissural, and projection fibers, the three fiber groups typically investigated in tractography studies, were successfully reconstructed and overlaid on SWI or fractional anisotropy maps. The fibers derived from DTI correlate well with their description in anatomical textbooks. Our results demonstrate the feasibility of using postmortem DTI data to reconstruct the main white matter tracts of the equine brain. Further DTI acquisitions and corresponding dissections of equine brains will be necessary to validate these findings and create an equine stereotaxic white matter atlas that could be used in future neuroimaging research.

On the cortical connectivity in the macaque brain: A comparison of diffusion tractography and histological tracing data

Diffusion-weighted magnetic resonance imaging (DW-MRI) tractography is a non-invasive tool to probe neural connections and the structure of the white matter. It has been applied successfully in studies of neurological disorders and normal connectivity. Recent work has revealed that tractography produces a high incidence of false-positive connections, often from "bottleneck" white matter configurations. The rich literature in histological connectivity analysis studies in the macaque monkey enables quantitative evaluation of the performance of tractography algorithms. In this study, we use the intricate connections of frontal, cingulate, and parietal areas, well established by the anatomical literature, to derive a symmetrical histological connectivity matrix composed of 59 cortical areas. We evaluate the performance of fifteen diffusion tractography algorithms, including global, deterministic, and probabilistic state-of-the-art methods for the connectivity predictions of 1711 distinct pairs of areas, among which 680 are reported connected by the literature. The diffusion connectivity analysis was performed on a different ex-vivo macaque brain, acquired using multi-shell DW-MRI protocol, at high spatial and angular resolutions. Across all tested algorithms, the true-positive and true-negative connections were dominant over false-positive and false-negative connections, respectively. Moreover, three-quarters of streamlines had endpoints location in agreement with histological data, on average. Furthermore, probabilistic streamline tractography algorithms show the best performances in predicting which areas are connected. Altogether, we propose a method for quantitative evaluation of tractography algorithms, which aims at improving the sensitivity and the specificity of diffusion-based connectivity analysis. Overall, those results confirm the usefulness of tractography in predicting connectivity, although errors are produced. Many of the errors result from bottleneck white matter configurations near the cortical grey matter and should be the target of future implementation of methods.

Histopathology of diffusion-weighted imaging-positive lesions in cerebral amyloid angiopathy

Small subclinical hyperintense lesions are frequently encountered on brain diffusion-weighted imaging (DWI) scans of patients with cerebral amyloid angiopathy (CAA). Interpretation of these DWI+ lesions, however, has been limited by absence of histopathological examination. We aimed to determine whether DWI+ lesions represent acute microinfarcts on histopathology in brains with advanced CAA, using a combined in vivo MRI-ex vivo MRI-histopathology approach. We first investigated the histopathology of a punctate cortical DWI+ lesion observed on clinical in vivo MRI 7 days prior to death in a CAA case. Subsequently, we assessed the use of ex vivo DWI to identify similar punctate cortical lesions post-mortem. Intact formalin-fixed hemispheres of 12 consecutive cases with CAA and three non-CAA controls were subjected to high-resolution 3 T ex vivo DWI and T2 imaging. Small cortical lesions were classified as either DWI+/T2+ or DWI-/T2+. A representative subset of lesions from three CAA cases was selected for detailed histopathological examination. The DWI+ lesion observed on in vivo MRI could be matched to an area with evidence of recent ischemia on histopathology. Ex vivo MRI of the intact hemispheres revealed a total of 130 DWI+/T2+ lesions in 10/12 CAA cases, but none in controls (p = 0.022). DWI+/T2+ lesions examined histopathologically proved to be acute microinfarcts (classification accuracy 100%), characterized by presence of eosinophilic neurons on hematoxylin and eosin and absence of reactive astrocytes on glial fibrillary acidic protein-stained sections. In conclusion, we suggest that small DWI+ lesions in CAA represent acute microinfarcts. Furthermore, our findings support the use of ex vivo DWI as a method to detect acute microinfarcts post-mortem, which may benefit future histopathological investigations on the etiology of microinfarcts.