Kurt G. Schilling

Comparison of 3D orientation distribution functions measured with confocal microscopy and diffusion MRI.

The ability of diffusion MRI (dMRI) fiber tractography to non-invasively map three-dimensional (3D) anatomical networks in the human brain has made it a valuable tool in both clinical and research settings. However, there are many assumptions inherent to any tractography algorithm that can limit the accuracy of the reconstructed fiber tracts. Among them is the assumption that the diffusion-weighted images accurately reflect the underlying fiber orientation distribution (FOD) in the MRI voxel. Consequently, validating dMRIs ability to assess the underlying fiber orientation in each voxel is critical for its use as a biomedical tool. Here, using post-mortem histology and confocal microscopy, we present a method to perform histological validation of orientation functions in 3D, which has previously been limited to two-dimensional analysis of tissue sections. We demonstrate the ability to extract the 3D FOD from confocal z-stacks, and quantify the agreement between the MRI estimates of orientation information obtained using constrained spherical deconvolution (CSD) and the true geometry of the fibers. We find an orientation error of approximately 6° in voxels containing nearly parallel fibers, and 10-11° in crossing fiber regions, and note that CSD was unable to resolve fibers crossing at angles below 60° in our dataset. This is the first time that the 3D white matter orientation distribution is calculated from histology and compared to dMRI. Thus, this technique serves as a gold standard for dMRI validation studies - providing the ability to determine the extent to which the dMRI signal is consistent with the histological FOD, and to establish how well different dMRI models can predict the ground truth FOD.

Confirmation of a gyral bias in diffusion MRI fiber tractography

Diffusion MRI fiber tractography has been increasingly used to map the structural connectivity of the human brain. However, this technique is not without limitations; for example, there is a growing concern over anatomically correlated bias in tractography findings. In this study, we demonstrate that there is a bias for fiber tracking algorithms to terminate preferentially on gyral crowns, rather than the banks of sulci. We investigate this issue by comparing diffusion MRI (dMRI) tractography with equivalent measures made on myelin‐stained histological sections. We begin by investigating the orientation and trajectories of axons near the white matter/gray matter boundary, and the density of axons entering the cortex at different locations along gyral blades. These results are compared with dMRI orientations and tract densities at the same locations, where we find a significant gyral bias in many gyral blades across the brain. This effect is shown for a range of tracking algorithms, both deterministic and probabilistic, and multiple diffusion models, including the diffusion tensor and a high angular resolution diffusion imaging technique. Additionally, the gyral bias occurs for a range of diffusion weightings, and even for very high‐resolution datasets. The bias could significantly affect connectivity results using the current generation of tracking algorithms.

Histological validation of diffusion MRI fiber orientation distributions and dispersion

ABSTRACT Diffusion magnetic resonance imaging (dMRI) is widely used to probe tissue microstructure, and is currently the only non‐invasive way to measure the brains fiber architecture. While a large number of approaches to recover the intra‐voxel fiber structure have been utilized in the scientific community, a direct, 3D, quantitative validation of these methods against relevant histological fiber geometries is lacking. In this study, we investigate how well different high angular resolution diffusion imaging (HARDI) models and reconstruction methods predict the ground‐truth histologically defined fiber orientation distribution (FOD), as well as investigate their behavior over a range of physical and experimental conditions. The dMRI methods tested include constrained spherical deconvolution (CSD), Q‐ball imaging (QBI), diffusion orientation transform (DOT), persistent angular structure (PAS), and neurite orientation dispersion and density imaging (NODDI) methods. Evaluation criteria focus on overall agreement in FOD shape, correct assessment of the number of fiber populations, and angular accuracy in orientation. In addition, we make comparisons of the histological orientation dispersion with the fiber spread determined from the dMRI methods. As a general result, no HARDI method outperformed others in all quality criteria, with many showing tradeoffs in reconstruction accuracy. All reconstruction techniques describe the overall continuous angular structure of the histological FOD quite well, with good to moderate correlation (median angular correlation coefficient > 0.70) in both single‐ and multiple‐fiber voxels. However, no method is consistently successful at extracting discrete measures of the number and orientations of FOD peaks. The major inaccuracies of all techniques tend to be in extracting local maxima of the FOD, resulting in either false positive or false negative peaks. Median angular errors are ˜10° for the primary fiber direction and ˜20° for the secondary fiber, if present. For most methods, these results did not vary strongly over a wide range of acquisition parameters (number of diffusion weighting directions and b value). Regardless of acquisition parameters, all methods show improved successes at resolving multiple fiber compartments in a voxel when fiber populations cross at near‐orthogonal angles, with no method adequately capturing low to moderate angle (<60°) crossing fibers. Finally, most methods are limited in their ability to capture orientation dispersion, resulting in low to moderate, yet statistically significant, correlation with histologically‐derived dispersion with both HARDI and NODDI methodologies. Together, these results provide quantitative measures of the reliability and limitations of dMRI reconstruction methods and can be used to identify relative advantages of competing approaches as well as potential strategies for improving accuracy. Highlights3D histological validation of diffusion MRI measures of fiber orientation.All methods capture the overall structure of the FOD quite well.Most inaccuracies occur when extracting discrete peaks from the FOD.No method consistently resolves fibers crossing at low to moderate angles.Measures of dispersion show modest correlation with histological measures.

Tests of cortical parcellation based on white matter connectivity using diffusion tensor imaging

ABSTRACT The cerebral cortex is conventionally divided into a number of domains based on cytoarchitectural features. Diffusion tensor imaging (DTI) enables noninvasive parcellation of the cortex based on white matter connectivity patterns. However, the correspondence between DTI‐connectivity‐based and cytoarchitectural parcellation has not been systematically established. In this study, we compared histological parcellation of New World monkey neocortex to DTI‐ connectivity‐based classification and clustering in the same brains. First, we used supervised classification to parcellate parieto‐frontal cortex based on DTI tractograms and the cytoarchitectural prior (obtained using Nissl staining). We performed both within and across sample classification, showing reasonable classification performance in both conditions. Second, we used unsupervised clustering to parcellate the cortex and compared the clusters to the cytoarchitectonic standard. We then explored the similarities and differences with several post‐hoc analyses, highlighting underlying principles that drive the DTI‐connectivity‐based parcellation. The differences in parcellation between DTI‐connectivity and Nissl histology probably represent both DTIs bias toward easily‐tracked bundles and true differences between cytoarchitectural and connectivity defined domains. DTI tractograms appear to cluster more according to functional networks, rather than mapping directly onto cytoarchitectonic domains. Our results show that caution should be used when DTI‐tractography classification, based on data from another brain, is used as a surrogate for cytoarchitectural parcellation. HIGHLIGHTSDTI‐connectivity‐based parcellation and Nissl histology are compared.Intra‐subject supervised DTI parcellation has up to 87% agreement with Nissl histology.Inter‐subject supervised DTI parcellation has up to 71% agreement with Nissl histology.Unsupervised DTI parcellation has 39% agreement with Nissl histology.Differences are from DTI errors and true differences between connectivity‐defined and cytoarchitectural domains.

Reproducibility and variation of diffusion measures in the squirrel monkey brain, in vivo and ex vivo

PURPOSE Animal models are needed to better understand the relationship between diffusion MRI (dMRI) and the underlying tissue microstructure. One promising model for validation studies is the common squirrel monkey, Saimiri sciureus. This study aims to determine (1) the reproducibility of in vivo diffusion measures both within and between subjects; (2) the agreement between in vivo and ex vivo data acquired from the same specimen and (3) normal diffusion values and their variation across brain regions. METHODS Data were acquired from three healthy squirrel monkeys, each imaged twice in vivo and once ex vivo. Reproducibility of fractional anisotropy (FA), mean diffusivity (MD), and principal eigenvector (PEV) was assessed, and normal values were determined both in vivo and ex vivo. RESULTS The calculated coefficients of variation (CVs) for both intra-subject and inter-subject MD were below 10% (low variability) while FA had a wider range of CVs, 2-14% intra-subject (moderate variability), and 3-31% inter-subject (high variability). MD in ex vivo tissue was lower than in vivo (30%-50% decrease), while FA values increased in all regions (30-39% increase). The mode of angular differences between in vivo and ex vivo PEVs was 12 degrees. CONCLUSION This study characterizes the diffusion properties of the squirrel monkey brain and serves as the groundwork for using the squirrel monkey, both in vivo and ex vivo, as a model for diffusion MRI studies.

The VALiDATe29 MRI Based Multi-Channel Atlas of the Squirrel Monkey Brain

We describe the development of the first digital atlas of the normal squirrel monkey brain and present the resulting product, VALiDATe29. The VALiDATe29 atlas is based on multiple types of magnetic resonance imaging (MRI) contrast acquired on 29 squirrel monkeys, and is created using unbiased, nonlinear registration techniques, resulting in a population-averaged stereotaxic coordinate system. The atlas consists of multiple anatomical templates (proton density, T1, and T2* weighted), diffusion MRI templates (fractional anisotropy and mean diffusivity), and ex vivo templates (fractional anisotropy and a structural MRI). In addition, the templates are combined with histologically defined cortical labels, and diffusion tractography defined white matter labels. The combination of intensity templates and image segmentations make this atlas suitable for the fundamental atlas applications of spatial normalization and label propagation. Together, this atlas facilitates 3D anatomical localization and region of interest delineation, and enables comparisons of experimental data across different subjects or across different experimental conditions. This article describes the atlas creation and its contents, and demonstrates the use of the VALiDATe29 atlas in typical applications. The atlas is freely available to the scientific community.

A brain MRI atlas of the common squirrel monkey, Saimiri sciureus

The common squirrel monkey, Saimiri sciureus, is a New World monkey with functional and microstructural organization of central nervous system similar to that of humans. It is one of the most commonly used South American primates in biomedical research. Unlike its Old World macaque cousins, no digital atlases have described the organization of the squirrel monkey brain. Here, we present a multi-modal magnetic resonance imaging (MRI) atlas constructed from the brain of an adult female squirrel monkey. In vivo MRI acquisitions include high resolution T2 structural imaging and low resolution diffusion tensor imaging. Ex vivo MRI acquisitions include high resolution T2 structural imaging and high resolution diffusion tensor imaging. Cortical regions were manually annotated on the co-registered volumes based on published histological sections.

Can increased spatial resolution solve the crossing fiber problem for diffusion MRI

It is now widely recognized that voxels with crossing fibers or complex geometrical configurations present a challenge for diffusion MRI (dMRI) reconstruction and fiber tracking, as well as microstructural modeling of brain tissues. This “crossing fiber” problem has been estimated to affect anywhere from 30% to as many as 90% of white matter voxels, and it is often assumed that increasing spatial resolution will decrease the prevalence of voxels containing multiple fiber populations. The aim of this study is to estimate the extent of the crossing fiber problem as we progressively increase the spatial resolution, with the goal of determining whether it is possible to mitigate this problem with higher resolution spatial sampling. This is accomplished using ex vivo MRI data of the macaque brain, followed by histological analysis of the same specimen to validate these measurements, as well as to extend this analysis to resolutions not yet achievable in practice with MRI. In both dMRI and histology, we find unexpected results: the prevalence of crossing fibers increases as we increase spatial resolution. The problem of crossing fibers appears to be a fundamental limitation of dMRI associated with the complexity of brain tissue, rather than a technical problem that can be overcome with advances such as higher fields and stronger gradients.

Empirical consideration of the effects of acquisition parameters and analysis model on clinically feasible q-ball imaging

Q-ball imaging (QBI) is a popular high angular resolution diffusion imaging (HARDI) technique used to study brain architecture in vivo. Simulation and phantom-based studies suggest that QBI results are affected by the b-value, the number of diffusion weighting directions, and the signal-to-noise ratio (SNR). However, optimal acquisition schemes for QBI in clinical settings are largely undetermined given empirical (observed) imaging considerations. In this study, we acquire a HARDI dataset at five b-values with 11 repetitions on a single subject to investigate the effects of acquisition scheme and subsequent analysis models on the accuracy and precision of measures of tissue composition and fiber orientation derived from clinically feasible QBI at 3T. Clinical feasibility entails short scan protocols - less than 5minutes in the current study - resulting in lower SNR, lower b-values, and fewer diffusion directions than are typical in most QBI protocols with research applications, where time constraints are less prevalent. In agreement with previous studies, we find that the b-value and number of diffusion directions impact the magnitude and variation of QBI indices in both white matter and gray matter regions; however, QBI indices are most heavily dependent on the maximum order of the spherical harmonic (SH) series used to represent the diffusion orientation distribution function (ODF). Specifically, to ensure numerical stability and reduce the occurrence of false peaks and inflated anisotropy, we recommend oversampling by at least 8-12 more diffusion directions than the number of estimated coefficients for a given SH order. In addition, in an equal scan time comparison of QBI accuracy, we find that increasing the directional resolution of the HARDI dataset is preferable to repeating observations; however, our results indicate that as few as 32 directions at a low b-value (1000s/mm2) captures most of the angular information in the q-ball ODF. Our findings provide guidance for determining an optimal acquisition scheme for QBI in the low SNR and low scan time regime, and suggest that care must be taken when choosing the basis functions used to represent the QBI ODF.

Integrating histology and MRI in the first digital brain of common squirrel monkey, Saimiri sciureus

This effort is a continuation of development of a digital brain atlas of the common squirrel monkey, Saimiri sciureus, a New World monkey with functional and microstructural organization of central nervous system similar to that of humans. Here, we present the integration of histology with multi-modal magnetic resonance imaging (MRI) atlas constructed from the brain of an adult female squirrel monkey. The central concept of this work is to use block face photography to establish an intermediate common space in coordinate system which preserves the high resolution in-plane resolution of histology while enabling 3-D correspondence with MRI. In vivo MRI acquisitions include high resolution T2 structural imaging (300 μm isotropic) and low resolution diffusion tensor imaging (600 um isotropic). Ex vivo MRI acquisitions include high resolution T2 structural imaging and high resolution diffusion tensor imaging (both 300 μm isotropic). Cortical regions were manually annotated on the co-registered volumes based on published histological sections in-plane. We describe mapping of histology and MRI based data of the common squirrel monkey and construction of a viewing tool that enable online viewing of these datasets. The previously descried atlas MRI is used for its deformation to provide accurate conformation to the MRI, thus adding information at the histological level to the MRI volume. This paper presents the mapping of single 2D image slice in block face as a proof of concept and this can be extended to map the atlas space in 3D coordinate system as part of the future work and can be loaded to an XNAT system for further use.