Histologically derived fiber response functions for diffusion MRI vary across white matter fibers-An ex vivo validation study in the squirrel monkey brain

Validation of Structure Tensor Analysis for Orientation Estimation in Brain Tissue Microscopy

Accurate localization of white matter pathways using diffusion MRI is critical to investigating brain connectivity, but the accuracy of current methods is not thoroughly understood. A fruitful approach to validating accuracy is to consider microscopy data that have been co-registered with MRI of post mortem samples. In this setting, structure tensor analysis is a standard approach to computing local orientations for validation. However, structure tensor analysis itself has not been well-validated and is subject to uncertainty in its angular resolution, and selectivity to specific spatial scales. In this work, we conducted a simulation study to investigate the accuracy of using structure tensors to estimate the orientations of fibers arranged in configurations with and without crossings. We examined a range of simulated conditions, with a focus on investigating the method's behavior on images with anisotropic resolution, which is particularly common in microscopy data acquisition. We also analyzed 2D and 3D optical microscopy data. Our results show that parameter choice in structure tensor analysis has relatively little effect on accuracy for estimating single orientations, although accuracy decreases with anisotropy. On the other hand, when estimating the orientations of crossing fibers, the choice of parameters becomes critical, and poor choices result in orientation estimates that are essentially random. This work provides a set of recommendations for researchers seeking to apply structure tensor analysis effectively in the study of axonal orientations in brain imaging data and quantifies the method's limitations, particularly in the case of anisotropic data.

White Matter Tract Crossing and Bottleneck Regions in the Fetal Brain

There is a growing interest in using diffusion MRI to study the white matter tracts and structural connectivity of the fetal brain. Recent progress in data acquisition and processing suggests that this imaging modality has a unique role in elucidating the normal and abnormal patterns of neurodevelopment in utero. However, there have been no efforts to quantify the prevalence of crossing tracts and bottleneck regions, important issues that have been investigated for adult brains. In this work, we determined the brain regions with crossing tracts and bottlenecks between 23 and 36 gestational weeks. We performed probabilistic tractography on 62 fetal brain scans and extracted a set of 51 distinct white matter tracts, which we grouped into 10 major tract bundle groups. We analyzed the results to determine the patterns of tract crossings and bottlenecks. Our results showed that 20%-25% of the white matter voxels included two or three crossing tracts. Bottlenecks were more prevalent. Between 75% and 80% of the voxels were characterized as bottlenecks, with more than 40% of the voxels involving four or more tracts. These results highlight the relevance of these regions to key developmental processes, specifically, the dispersion of projection fibers, the protracted growth of commissural pathways, and the emergence of association tracts that contribute to the formation of complex intersection regions. These developmental interactions lead to a high prevalence of crossing fibers and bottleneck areas, reflecting the intricate organization required for establishing structural and functional connectivity. Additionally, our results highlight the challenge of fetal brain tractography and structural connectivity assessment and call for innovative image acquisition and analysis methods to mitigate these problems.

Diffusion MRI approaches for investigating microstructural complexity in a rat model of traumatic brain injury

Our study explores the potential of conventional and advanced diffusion MRI techniques including diffusion tensor imaging (DTI), and single-shell 3-tissue constrained spherical deconvolution (SS3T-CSD) to investigate complex microstructural changes following severe traumatic brain injury in rats at a chronic phase. Rat brains after sham-operation or lateral fluid percussion (LFP) injury were scanned ex vivo in a 9.4 T scanner. Our region-of-interest-based approach of tensor-, and SS3T-CSD derived fixel-, 3-tissue signal fraction maps were sensitive to changes in both white matter (WM) and grey matter (GM) areas. Tensor-based measures, such as fractional anisotropy (FA) and radial diffusivity (RD), detected more changes in WM and GM areas as compared to fixel-based measures including apparent fiber density (AFD), peak FOD amplitude and primary fiber bundle density, while 3-tissue signal fraction maps revealed distinct changes in WM, GM, and phosphate-buffered saline (PBS) fractions highlighting the complex tissue microstructural alterations post-trauma. Track-weighted imaging demonstrated changes in track morphology including reduced curvature and average pathlength distal from the primary lesion in severe TBI rats. In histological analysis, changes in the diffusion MRI measures could be associated to decreased myelin density, loss of myelinated axons, and increased cellularity, revealing progressive microstructural alterations in these brain areas five months after injury. Overall, this study highlights the use of combined conventional and advanced diffusion MRI measures to obtain more precise insights into the complex tissue microstructural alterations in chronic phase of severe brain injury.

Prevalence of white matter pathways coming into a single white matter voxel orientation: The bottleneck issue in tractography

Characterizing and understanding the limitations of diffusion MRI fiber tractography is a prerequisite for methodological advances and innovations which will allow these techniques to accurately map the connections of the human brain. The so-called "crossing fiber problem" has received tremendous attention and has continuously triggered the community to develop novel approaches for disentangling distinctly oriented fiber populations. Perhaps an even greater challenge occurs when multiple white matter bundles converge within a single voxel, or throughout a single brain region, and share the same parallel orientation, before diverging and continuing towards their final cortical or sub-cortical terminations. These so-called "bottleneck" regions contribute to the ill-posed nature of the tractography process, and lead to both false positive and false negative estimated connections. Yet, as opposed to the extent of crossing fibers, a thorough characterization of bottleneck regions has not been performed. The aim of this study is to quantify the prevalence of bottleneck regions. To do this, we use diffusion tractography to segment known white matter bundles of the brain, and assign each bundle to voxels they pass through and to specific orientations within those voxels (i.e. fixels). We demonstrate that bottlenecks occur in greater than 50-70% of fixels in the white matter of the human brain. We find that all projection, association, and commissural fibers contribute to, and are affected by, this phenomenon, and show that even regions traditionally considered "single fiber voxels" often contain multiple fiber populations. Together, this study shows that a majority of white matter presents bottlenecks for tractography which may lead to incorrect or erroneous estimates of brain connectivity or quantitative tractography (i.e., tractometry), and underscores the need for a paradigm shift in the process of tractography and bundle segmentation for studying the fiber pathways of the human brain.

Tractography methods and findings in brain tumors and traumatic brain injury

White matter fiber tracking using diffusion magnetic resonance imaging (dMRI) provides a noninvasive approach to map brain connections, but improving anatomical accuracy has been a significant challenge since the birth of tractography methods. Utilizing tractography in brain studies therefore requires understanding of its technical limitations to avoid shortcomings and pitfalls. This review explores tractography limitations and how different white matter pathways pose different challenges to fiber tracking methodologies. We summarize the pros and cons of commonly-used methods, aiming to inform how tractography and its related analysis may lead to questionable results. Extending these experiences, we review the clinical utilization of tractography in patients with brain tumors and traumatic brain injury, starting from tensor-based tractography to more advanced methods. We discuss current limitations and highlight novel approaches in the context of these two conditions to inform future tractography developments.

Design and characterization of a 3D-printed axon-mimetic phantom for diffusion MRI

PURPOSE

To introduce and characterize inexpensive and easily produced 3D-printed axon-mimetic diffusion MRI phantoms in terms of pore geometry and diffusion kurtosis imaging metrics.

METHODS

Phantoms were 3D-printed with a composite printing material that, after the dissolution of the polyvinyl alcohol, exhibits microscopic fibrous pores. Confocal microscopy and synchrotron phase-contrast micro-CT imaging were performed to visualize and assess the pore sizes. Diffusion MRI scans of four identical phantoms and phantoms with varying print parameters in water were performed at 9.4 T. Diffusion kurtosis imaging was fit to both data sets and used to assess the reproducibility between phantoms and effects of print parameters on diffusion kurtosis imaging metrics. Identical scans were performed 25 and 76 days later, to test their stability.

RESULTS

Segmentation of pores in three microscopy images yielded a mean, median, and SD of equivalent pore diameters of 7.57 μm, 3.51 μm, and 12.13 μm, respectively. Phantoms had T₁ /T₂ = 2 seconds/180 ms, and those with identical parameters showed a low coefficient of variation (~10%) in mean diffusivity (1.38 × 10^-3 mm² /s) and kurtosis (0.52) metrics and radial diffusivity (1.01 × 10^-3 mm² /s) and kurtosis (1.13) metrics. Printing temperature and speed had a small effect on diffusion kurtosis imaging metrics (< 16%), whereas infill density had a larger and more variable effect (> 16%). The stability analysis showed small changes over 2.5 months (< 7%).

CONCLUSION

Three-dimension-printed axon-mimetic phantoms can mimic the fibrous structure of axon bundles on a microscopic scale, serving as complex, anisotropic diffusion MRI phantoms.

Effective connectivity and spatial selectivity-dependent fMRI changes elicited by microstimulation of pulvinar and LIP

The thalamic pulvinar and the lateral intraparietal area (LIP) share reciprocal anatomical connections and are part of an extensive cortical and subcortical network involved in spatial attention and oculomotor processing. The goal of this study was to compare the effective connectivity of dorsal pulvinar (dPul) and LIP and to probe the dependency of microstimulation effects on task demands and spatial tuning properties of a given brain region. To this end, we applied unilateral electrical microstimulation in the dPul (mainly medial pulvinar) and LIP in combination with event-related BOLD fMRI in monkeys performing fixation and memory-guided saccade tasks. Microstimulation in both dPul and LIP enhanced task-related activity in monosynaptically-connected fronto-parietal cortex and along the superior temporal sulcus (STS) including putative face patch locations, as well as in extrastriate cortex. LIP microstimulation elicited strong activity in the opposite homotopic LIP while no homotopic activation was found with dPul stimulation. Both dPul and LIP stimulation also elicited activity in several heterotopic cortical areas in the opposite hemisphere, implying polysynaptic propagation of excitation. Despite extensive activation along the intraparietal sulcus evoked by LIP stimulation, there was a difference in frontal and occipital connectivity elicited by posterior and anterior LIP stimulation sites. Comparison of dPul stimulation with the adjacent but functionally dissimilar ventral pulvinar also showed distinct connectivity. On the level of single trial timecourses within each region of interest (ROI), most ROIs did not show task-dependence of stimulation-elicited response modulation. Across ROIs, however, there was an interaction between task and stimulation, and task-specific correlations between the initial spatial selectivity and the magnitude of stimulation effect were observed. Consequently, stimulation-elicited modulation of task-related activity was best fitted by an additive model scaled down by the initial response amplitude. In summary, we identified overlapping and distinct patterns of thalamocortical and corticocortical connectivity of pulvinar and LIP, highlighting the dorsal bank and fundus of STS as a prominent node of shared circuitry. Spatial task-specific and partly polysynaptic modulations of cue and saccade planning delay period activity in both hemispheres exerted by unilateral pulvinar and parietal stimulation provide insight into the distributed interhemispheric processing underlying spatial behavior.

Fiber orientation distribution from diffusion MRI: Effects of inaccurate response function calibration

Background and Purpose

Diffusion MRI of the brain enables to quantify white matter fiber orientations noninvasively. Several approaches have been proposed to estimate such characteristics from diffusion MRI data with spherical deconvolution being one of the most widely used methods. Spherical deconvolution requires to define––or derive from the data––a response function, which is used to compute the fiber orientation distribution (FOD). Different characteristics of the response function are expected to affect the FOD computation and the subsequent fiber tracking.

Methods

In this work, we explored the effects of inaccuracies in the shape factors of the response function on the FOD characteristics.

Results

With simulations, we show that the apparent fiber density could be doubled in the presence of underestimated shape factors in the response functions, whereas the overestimation of the shape factor will cause more spurious peaks in the FOD, especially when the signal‐to‐noise ratio is below 15. Moreover, crossing fiber populations with a separation angle smaller than 60° were more sensitive to inaccuracies in the response function than fiber populations with more orthogonal separation angles. Results with in vivo data demonstrate angular deviations in the FODs and spurious peaks as a result of modified shape factors of the response function, while the reconstruction of the main parts of fiber bundles is well preserved.

Conclusions

This work sheds light on how specific aspects of the response function shape can affect the estimated FODs, and highlights the importance of a proper calibration/definition of the response function.

The squirrel monkey model in clinical neuroscience

Clinical neuroscience research relying on animal models brought valuable translational insights into the function and pathologies of the human brain. The anatomical, physiological, and behavioural similarities between humans and mammals have prompted researchers to study cerebral mechanisms at different levels to develop and test new treatments. The vast majority of biomedical research uses rodent models, which are easily manipulable and have a broadly resembling organisation to the human nervous system but cannot satisfactorily mimic some disorders. For these disorders, macaque monkeys have been used as they have a more comparable central nervous system. Still, this research has been hampered by limitations, including high costs and reduced samples. This review argues that a squirrel monkey model might bridge the gap by complementing translational research from rodents, macaque, and humans. With the advent of promising new methods such as ultrasound imaging, tool miniaturisation, and a shift towards open science, the squirrel monkey model represents a window of opportunity that will potentially fuel new translational discoveries in the diagnosis and treatment of brain pathologies.

Assessment of the structural complexity of diffusion MRI voxels using 3D electron microscopy in the rat brain

Validation and interpretation of diffusion magnetic resonance imaging (dMRI) requires detailed understanding of the actual microstructure restricting the diffusion of water molecules. In this study, we used serial block-face scanning electron microscopy (SBEM), a three-dimensional electron microscopy (3D-EM) technique, to image seven white and grey matter volumes in the rat brain. SBEM shows excellent contrast of cellular membranes, which are the major components restricting the diffusion of water in tissue. Additionally, we performed 3D structure tensor (3D-ST) analysis on the SBEM volumes and parameterised the resulting orientation distributions using Watson and angular central Gaussian (ACG) probability distributions as well as spherical harmonic (SH) decomposition. We analysed how these parameterisations described the underlying orientation distributions and compared their orientation and dispersion with corresponding parameters from two dMRI methods, neurite orientation dispersion and density imaging (NODDI) and constrained spherical deconvolution (CSD). Watson and ACG parameterisations and SH decomposition captured well the 3D-ST orientation distributions, but ACG and SH better represented the distributions due to its ability to model asymmetric dispersion. The dMRI parameters corresponded well with the 3D-ST parameters in the white matter volumes, but the correspondence was less evident in the more complex grey matter. SBEM imaging and 3D-ST analysis also revealed that the orientation distributions were often not axially symmetric, a property neatly captured by the ACG distribution. Overall, the ability of SBEM to image diffusion barriers in intricate detail, combined with 3D-ST analysis and parameterisation, provides a step forward toward interpreting and validating the dMRI signals in complex brain tissue microstructure.

Joint modelling of diffusion MRI and microscopy

The combination of diffusion MRI (dMRI) with microscopy provides unique opportunities to study microstructural features of tissue, particularly when acquired in the same sample. Microscopy is frequently used to validate dMRI microstructure models, addressing the indirect nature of dMRI signals. Typically, these modalities are analysed separately, and microscopy is taken as a gold standard against which dMRI-derived parameters are validated. Here we propose an alternative approach in which we combine dMRI and microscopy data obtained from the same tissue sample to drive a single, joint model. This simultaneous analysis allows us to take advantage of the breadth of information provided by complementary data acquired from different modalities. By applying this framework to a spherical-deconvolution analysis, we are able to overcome a known degeneracy between fibre dispersion and radial diffusion. Spherical-deconvolution based approaches typically estimate a global fibre response function to determine the fibre orientation distribution in each voxel. However, the assumption of a 'brain-wide' fibre response function may be challenged if the diffusion characteristics of white matter vary across the brain. Using a generative joint dMRI-histology model, we demonstrate that the fibre response function is dependent on local anatomy, and that current spherical-deconvolution based models may be overestimating dispersion and underestimating the number of distinct fibre populations per voxel.