Comparative validation of automated presurgical tractography based on constrained spherical deconvolution and diffusion tensor imaging with direct electrical stimulation

OBJECTIVES

Accurate presurgical brain mapping enables preoperative risk assessment and intraoperative guidance. This cross-sectional study investigated whether constrained spherical deconvolution (CSD) methods were more accurate than diffusion tensor imaging (DTI)-based methods for presurgical white matter mapping using intraoperative direct electrical stimulation (DES) as the ground truth.

METHODS

Five different tractography methods were compared (three DTI-based and two CSD-based) in 22 preoperative neurosurgical patients undergoing surgery with DES mapping. The corticospinal tract (CST, N = 20) and arcuate fasciculus (AF, N = 7) bundles were reconstructed, then minimum distances between tractograms and DES coordinates were compared between tractography methods. Receiver-operating characteristic (ROC) curves were used for both bundles. For the CST, binary agreement, linear modeling, and posthoc testing were used to compare tractography methods while correcting for relative lesion and bundle volumes.

RESULTS

Distance measures between 154 positive (functional response, pDES) and negative (no response, nDES) coordinates, and 134 tractograms resulted in 860 data points. Higher agreement was found between pDES coordinates and CSD-based compared to DTI-based tractograms. ROC curves showed overall higher sensitivity at shorter distance cutoffs for CSD (8.5 mm) compared to DTI (14.5 mm). CSD-based CST tractograms showed significantly higher agreement with pDES, which was confirmed by linear modeling and posthoc tests (PFWE < .05).

CONCLUSIONS

CSD-based CST tractograms were more accurate than DTI-based ones when validated using DES-based assessment of motor and sensory function. This demonstrates the potential benefits of structural mapping using CSD in clinical practice.

Fiber Orientation Estimation from X-ray Dark Field Images of Fiber Reinforced Polymers Using Constrained Spherical Deconvolution

The properties of fiber reinforced polymers are strongly related to the length and orientation of the fibers within the polymer matrix, the latter of which can be studied using X-ray computed tomography (XCT). Unfortunately, resolving individual fibers is challenging because they are small compared to the XCT voxel resolution and because of the low attenuation contrast between the fibers and the surrounding resin. To alleviate both problems, anisotropic dark field tomography via grating based interferometry (GBI) has been proposed. Here, the fiber orientations are extracted by applying a Funk-Radon transform (FRT) to the local scatter function. However, the FRT suffers from a low angular resolution, which complicates estimating fiber orientations for small fiber crossing angles. We propose constrained spherical deconvolution (CSD) as an alternative to the FRT to resolve fiber orientations. Instead of GBI, edge illumination phase contrast imaging is used because estimating fiber orientations with this technique has not yet been explored. Dark field images are generated by a Monte Carlo simulation framework. It is shown that the FRT cannot estimate the fiber orientation accurately for crossing angles smaller than 70∘, while CSD performs well down to a crossing angle of 50∘. In general, CSD outperforms the FRT in estimating fiber orientations.

Alterations in white matter fiber density associated with structural MRI and metabolic PET lesions following multimodal therapy in glioma patients

BACKGROUND

In glioma patients, multimodality therapy and recurrent tumor can lead to structural brain tissue damage characterized by pathologic findings in MR and PET imaging. However, little is known about the impact of different types of damage on the fiber architecture of the affected white matter.

PATIENTS AND METHODS

This study included 121 pretreated patients (median age, 52 years; ECOG performance score, 0 in 48%, 1-2 in 51%) with histomolecularly characterized glioma (WHO grade IV glioblastoma, n=81; WHO grade III anaplastic astrocytoma, n=28; WHO grade III anaplastic oligodendroglioma, n=12), who had a resection, radiotherapy, alkylating chemotherapy, or combinations thereof. After a median follow-up time of 14 months (range, 1-214 months), anatomic MR and O-(2-[18F]fluoroethyl)-L-tyrosine (FET) PET images were acquired on a 3T hybrid PET/MR scanner. Post-therapeutic findings comprised resection cavities, regions with contrast enhancement or increased FET uptake and T2/FLAIR hyperintensities. Local fiber density was determined from high angular-resolution diffusion-weighted imaging and advanced tractography methods. A cohort of 121 healthy subjects selected from the 1000BRAINS study matched for age, gender and education served as a control group.

RESULTS

Lesion types differed in both affected tissue volumes and relative fiber densities compared to control values (resection cavities: median volume 20.9 mL, fiber density 16% of controls; contrast-enhanced lesions: 7.9 mL, 43%; FET uptake areas: 30.3 mL, 49%; T2/FLAIR hyperintensities: 53.4 mL, 57%, p<0.001). In T2/FLAIR-hyperintense lesions caused by peritumoral edema due to recurrent glioma (n=27), relative fiber density was as low as in lesions associated with radiation-induced gliosis (n=13, 48% vs. 53%, p=0.17). In regions with pathologically increased FET uptake, local fiber density was inversely related (p=0.005) to the extent of uptake. Total fiber loss associated with contrast-enhanced lesions (p=0.006) and T2/FLAIR hyperintense lesions (p=0.013) had a significant impact on overall ECOG score.

CONCLUSIONS

These results suggest that apart from resection cavities, reduction in local fiber density is greatest in contrast-enhancing recurrent tumors, but total fiber loss induced by edema or gliosis has an equal detrimental effect on the patients' performance status due to the larger volume affected.

Effects of tDCS on Language Recovery in Post-Stroke Aphasia: A Pilot Study Investigating Clinical Parameters and White Matter Change with Diffusion Imaging

Objectives: In this pilot study we investigated the effects of transcranial direct current stimulation (tDCS) on language recovery in the subacute stage of post-stroke aphasia using clinical parameters and diffusion imaging with constrained spherical deconvolution-based tractography. Methods: The study included 21 patients with subacute post-stroke aphasia. Patients were randomly classified into two groups with a ratio of 2:1 to receive real tDCS or sham tDCS as placebo control. Patients received 10 sessions (5/week) bi-hemispheric tDCS treatments over the left affected Broca's area (anodal electrode) and over the right unaffected Broca's area (cathodal stimulation). Aphasia score was assessed clinically using the language section of the Hemispheric Stroke Scale (HSS) before and after treatment sessions. Diffusion imaging and tractography were performed for seven patients of the real group, both before and after the 10th session. Dissection of language-related white matter tracts was achieved, and diffusion measures were extracted. A paired Student's t-test was used to compare the clinical recovery and diffusion measures of the dissected tracts both pre- and post- treatment. The partial correlation between changes in diffusion measures and the language improvements was calculated. Results: At baseline assessment, there were no significant differences between groups in demographic and clinical HSS language score. No significant clinical recovery in HSS was evident in the sham group. However, significant improvements in the different components of HSS were only observed in patients receiving real tDCS. Associated significant increase in the fractional anisotropy of the right uncinate fasciculus and a significant reduction in the mean diffusivity of the right frontal aslant tract were reported. A significant positive correlation was found between the changes in the right uncinate fasciculus and fluency improvement. Conclusions: Aphasia recovery after bi-hemispheric transcranial direct current stimulation was associated with contralesional right-sided white matter changes at the subacute stage. These changes probably reflect neuroplasticity that could contribute to the recovery. Both the right uncinate fasciculus and right frontal aslant tract seem to be involved in aphasia recovery.

Tractographic and Microstructural Analysis of the Dentato-Rubro-Thalamo-Cortical Tracts in Children Using Diffusion MRI

The dentato-rubro-thalamo-cortical tract (DRTC) is the main outflow pathway of the cerebellum, contributing to a finely balanced corticocerebellar loop involved in cognitive and sensorimotor functions. Damage to the DRTC has been implicated in cerebellar mutism syndrome seen in up to 25% of children after cerebellar tumor resection. Multi-shell diffusion MRI (dMRI) combined with quantitative constrained spherical deconvolution tractography and multi-compartment spherical mean technique modeling was used to explore the frontocerebellar connections and microstructural signature of the DRTC in 30 healthy children. The highest density of DRTC connections were to the precentral (M1) and superior frontal gyri (F1), and from cerebellar lobules I-IV and IX. The first evidence of a topographic organization of anterograde projections to the frontal cortex at the level of the superior cerebellar peduncle (SCP) is demonstrated, with streamlines terminating in F1 lying dorsomedially in the SCP compared to those terminating in M1. The orientation dispersion entropy of DRTC regions appears to exhibit greater contrast than that shown by fractional anisotropy. Analysis of a separate reproducibility cohort demonstrates good consistency in the dMRI metrics described. These novel anatomical insights into this well-studied pathway may prove to be of clinical relevance in the surgical resection of cerebellar tumors.

The effect of gradient nonlinearities on fiber orientation estimates from spherical deconvolution of diffusion magnetic resonance imaging data

Gradient nonlinearities in magnetic resonance imaging (MRI) cause spatially varying mismatches between the imposed and the effective gradients and can cause significant biases in rotationally invariant diffusion MRI measures derived from, for example, diffusion tensor imaging. The estimation of the orientational organization of fibrous tissue, which is nowadays frequently performed with spherical deconvolution techniques ideally using higher diffusion weightings, can likewise be biased by gradient nonlinearities. We explore the sensitivity of two established spherical deconvolution approaches to gradient nonlinearities, namely constrained spherical deconvolution (CSD) and damped Richardson-Lucy (dRL). Additionally, we propose an extension of dRL to take into account gradient imperfections, without the need of data interpolation. Simulations show that using the effective b-matrix can improve dRL fiber orientation estimation and reduces angular deviations, while CSD can be more robust to gradient nonlinearity depending on the implementation. Angular errors depend on a complex interplay of many factors, including the direction and magnitude of gradient deviations, underlying microstructure, SNR, anisotropy of the effective response function, and diffusion weighting. Notably, angular deviations can also be observed at lower b-values in contrast to the perhaps common assumption that only high b-value data are affected. In in vivo Human Connectome Project data and acquisitions from an ultrastrong gradient (300 mT/m) scanner, angular differences are observed between applying and not applying the effective gradients in dRL estimation. As even small angular differences can lead to error propagation during tractography and as such impact connectivity analyses, incorporating gradient deviations into the estimation of fiber orientations should make such analyses more reliable.

Optimized rectification of fiber orientation density function

PURPOSE

To demonstrate an optimized rectification strategy for fiber orientation density functions (fODFs).

THEORY AND METHODS

In white matter, fODFs can be estimated with diffusion MRI. However, because of signal noise, imaging artifacts and other factors, experimentally determined fODFs may take on unphysical negative values in some directions. Here, we show how to rectify such fODFs to eliminate all negative values while minimizing the mean square difference between the original and rectified fODFs. The method is demonstrated for a mathematical model and for fODFs estimated from experimental human data using both constrained spherical deconvolution and fiber ball imaging. Comparison with an alternative nonoptimized rectification approach is also provided.

RESULTS

For the mathematical model, it is found that the optimized rectification procedure removes negative fODF values while at the same time reducing the mean square error. Relative to the alternative rectification approach, the optimized fODFs are substantially more accurate. For the experimental data, the optimized fODFs have a lower average fractional anisotropy axonal and often fewer small peaks than the original, unrectified fODFs. The calculation of optimized fODFs is straightforward where the main step is the finding of the root to an equation in one variable, as may be efficiently accomplished with the bisection method.

CONCLUSION

Unphysical negative fODF values can be easily eliminated in a manner that minimizes the mean square difference between the original and rectified fODFs. Optimized fODF rectification may be useful in applications for which negative values are problematic.

Measurement of white matter fiber-bundle cross-section in multiple sclerosis using diffusion-weighted imaging

BACKGROUND

When investigating white matter (WM) microstructure, the axonal fiber orientation should be considered. Constrained spherical deconvolution (CSD) is a diffusion-weighted imaging (DWI) method that estimates distribution of fibers within each imaging voxel.

OBJECTIVE

To study fiber-bundle cross-section (FC) as measured by CSD in multiple sclerosis (MS) patients versus healthy controls (HCs).

METHODS

DWI and three-dimensional (3D) T₁-weighted magnetic resonance imaging (MRI) were obtained from 45 MS patients and 45 HCs. We applied fixel-based morphometry analysis to assess differences of FC in MS against HCs and voxel-based analysis of fractional anisotropy (FA).

RESULTS

We found a significant widespread reduction of WM FC in MS compared to HCs. The decrease in FA was less extensive, mainly located in regions with high lesion occurrence such as the periventricular WM and the corpus callosum. Progressive MS patients showed a significant FC reduction in the right anterior cingulum, bilateral cerebellum, and in several mesencephalic and diencephalic regions compared to relapsing-remitting MS patients.

CONCLUSION

The CSD method can be applied in MS for a fiber-specific study of WM microstructure and quantification of FC. Fixel-based findings offered greater anatomical specificity and biological interpretability by identifying tract-specific differences and allowed substantial abnormalities to be detected.

Impact of b-value on estimates of apparent fibre density

Recent advances in diffusion magnetic resonance imaging (dMRI) analysis techniques have improved our understanding of fibre‐specific variations in white matter microstructure. Increasingly, studies are adopting multi‐shell dMRI acquisitions to improve the robustness of dMRI‐based inferences. However, the impact of b‐value choice on the estimation of dMRI measures such as apparent fibre density (AFD) derived from spherical deconvolution is not known. Here, we investigate the impact of b‐value sampling scheme on estimates of AFD. First, we performed simulations to assess the correspondence between AFD and simulated intra‐axonal signal fraction across multiple b‐value sampling schemes. We then studied the impact of sampling scheme on the relationship between AFD and age in a developmental population (n = 78) aged 8–18 (mean = 12.4, SD = 2.9 years) using hierarchical clustering and whole brain fixel‐based analyses. Multi‐shell dMRI data were collected at 3.0T using ultra‐strong gradients (300 mT/m), using 6 diffusion‐weighted shells ranging from b = 0 to 6,000 s/mm². Simulations revealed that the correspondence between estimated AFD and simulated intra‐axonal signal fraction was improved with high b‐value shells due to increased suppression of the extra‐axonal signal. These results were supported by in vivo data, as sensitivity to developmental age‐relationships was improved with increasing b‐value (b = 6,000 s/mm², median R² = .34; b = 4,000 s/mm², median R² = .29; b = 2,400 s/mm², median R² = .21; b = 1,200 s/mm², median R² = .17) in a tract‐specific fashion. Overall, estimates of AFD and age‐related microstructural development were better characterised at high diffusion‐weightings due to improved correspondence with intra‐axonal properties.

Fast learning of fiber orientation distribution function for MR tractography using convolutional neural network

PURPOSE

In diffusion-weighted magnetic resonance imaging (DW-MRI), the fiber orientation distribution function (fODF) is of great importance for solving complex fiber configurations to achieve reliable tractography throughout the brain, which ultimately facilitates the understanding of brain connectivity and exploration of neurological dysfunction. Recently, multi-shell multi-tissue constrained spherical deconvolution (MSMT-CSD) method has been explored for reconstructing full fODFs. To achieve a reliable fitting, similar to other model-based approaches, a large number of diffusion measurements is typically required for MSMT-CSD method. The prolonged acquisition is, however, not feasible in practical clinical routine and is prone to motion artifacts. To accelerate the acquisition, we proposed a method to reconstruct the fODF from downsampled diffusion-weighted images (DWIs) by leveraging the strong inference ability of the deep convolutional neural network (CNN).

METHODS

The method treats spherical harmonics (SH)-represented DWI signals and fODF coefficients as inputs and outputs, respectively. To compensate for the reduced gradient directions with reduced number of DWIs in acquisition in each voxel, its surrounding voxels are incorporated by the network for exploiting their spatial continuity. The resulting fODF coefficients are fitted with applying the CNN in a multi-target regression model. The network is composed of two convolutional layers and three fully connected layers. To obtain an initial evaluation of the method, we quantitatively measured its performance on a simulated dataset. Then, for in vivo tests, we employed data from 24 subjects from the Human Connectome Project (HCP) as training set and six subjects as test set. The performance of the proposed method was primarily compared to the super-resolved MSMT-CSD with the decreasing number of DWIs. The fODFs reconstructed by MSMT-CSD from all available 288 DWIs were used as training labels and the reference standard. The performance was quantitatively measured by the angular correlation coefficient (ACC) and the mean angular error (MAE).

RESULTS

For the simulated dataset, the proposed method exhibited the potential advantage over the model reconstruction. For the in vivo dataset, it achieved superior results over the MSMT-CSD in all the investigated cases, with its advantage more obvious when a limited number of DWIs were used. As the number of DWIs was reduced from 95 to 25, the median ACC ranged from 0.96 to 0.91 for the CNN, but 0.93 to 0.77 for the MSMT-CSD (with perfect score of 1). The angular error in the typical regions of interest (ROIs) was also much lower, especially in multi-fiber regions. The average MAE for the CNN method in regions containing one, two, three fibers was, respectively, 1.09°, 2.75°, and 8.35° smaller than the MSMT-CSD method. The visual inception of the fODF further confirmed this superiority. Moreover, the tractography results validated the effectiveness of the learned fODF, in preserving known major branching fibers with only 25 DWIs.

CONCLUSION

Experiments on HCP datasets demonstrated the feasibility of the proposed method in recovering fODFs from up to 11-fold reduced number of DWIs. The proposed method offers a new streamlined reconstruction procedure and exhibits promising potential in acquisition acceleration for the reconstruction of fODFs with good accuracy.

Crossing muscle fibers of the human tongue resolved in vivo using constrained spherical deconvolution

BACKGROUND

Surgical resection of tongue cancer may impair swallowing and speech. Knowledge of tongue muscle architecture affected by the resection could aid in patient counseling. Diffusion tensor imaging (DTI) enables reconstructions of muscle architecture in vivo. Reconstructing crossing fibers in the tongue requires a higher-order diffusion model.

PURPOSE

To develop a clinically feasible diffusion imaging protocol, which facilitates both DTI and constrained spherical deconvolution (CSD) reconstructions of tongue muscle architecture in vivo.

STUDY TYPE

Cross-sectional study.

SUBJECTS/SPECIMEN

One ex vivo bovine tongue resected en bloc from mandible to hyoid bone. Ten healthy volunteers (mean age 25.5 years; range 21-34 years; four female).

FIELD STRENGTH/SEQUENCE

Diffusion-weighted echo planar imaging at 3 T using a high-angular resolution diffusion imaging scheme acquired twice with opposing phase-encoding for B₀ -field inhomogeneity correction. The scan of the healthy volunteers was divided into four parts, in between which the volunteers were allowed to swallow, resulting in a total acquisition time of 10 minutes.

ASSESSMENT

The ability of resolving crossing muscle fibers using CSD was determined on the bovine tongue specimen. A reproducible response function was estimated and the optimal peak threshold was determined for the in vivo tongue. The quality of tractography of the in vivo tongue was graded by three experts.

STATISTICAL TESTS

The within-subject coefficient of variance was calculated for the response function. The qualitative results of the grading of DTI and CSD tractography were analyzed using a multilevel proportional odds model.

RESULTS

Fiber orientation distributions in the bovine tongue specimen showed that CSD was able to resolve crossing muscle fibers. The response function could be determined reproducibly in vivo. CSD tractography displayed significantly improved tractography compared with DTI tractography (P = 0.015).

DATA CONCLUSION

The 10-minute diffusion imaging protocol facilitates CSD fiber tracking with improved reconstructions of crossing tongue muscle fibers compared with DTI.

LEVEL OF EVIDENCE

2 Technical Efficacy: Stage 1 J. Magn. Reson. Imaging 2019;50:96-105.

White matter degeneration profile in the cognitive cortico-subcortical tracts in Parkinson's disease

BACKGROUND

In Parkinson's disease cognitive impairment is an early nonmotor feature, but it is still unclear why some patients are able to maintain their cognitive performance at normal levels, as quantified by neuropsychological tests, whereas others cannot. The objectives of this study were to perform a cross-sectional study and analyze the white matter changes in the cognitive and motor bundles in patients with Parkinson's disease.

METHODS

Sixteen Parkinson's disease patients with normal cognitive performance, 19 with mild cognitive impairment (based on their performance of 1.5 standard deviations below the healthy population mean), and 16 healthy controls were compared with respect to their tractography patterns between the cortical cognitive / motor regions and subcortical structures, using high angular resolution diffusion imaging and constrained spherical deconvolution computation.

RESULTS

Motor bundles showed decreased apparent fiber density in both PD groups, associated with a significant increase in diffusivity metrics, number of reconstructed streamlines, and track volumes, compared with healthy controls. By contrast, in the cognitive bundles, decreased fiber density in both Parkinson's groups was compounded by the absence of changes in diffusivity in patients with normal cognition, whereas patients with cognitive impairment had increased diffusivity metrics, lower numbers of reconstructed streamlines, and lower track volumes.

CONCLUSIONS

Both PD groups showed similar patterns of white matter neurodegeneration in the motor bundles, whereas cognitive bundles showed a distinct pattern: Parkinson's patients with normal cognition had white matter diffusivity metrics similar to healthy controls, whereas in patients with cognitive impairment white matter showed a neurodegeneration pattern. © 2018 International Parkinson and Movement Disorder Society.

Comparison of Diffusion-Weighted MRI Reconstruction Methods for Visualization of Cranial Nerves in Posterior Fossa Surgery

Diffusion-weighted imaging (DWI)-based tractography has gained increasing popularity as a method for detailed visualization of white matter (WM) tracts. Different imaging techniques, and more novel, advanced imaging methods provide significant WM structural detail. While there has been greater focus on improving tract visualization for larger WM pathways, the relative value of each method for cranial nerve reconstruction and how this methodology can assist surgical decision-making is still understudied. Images from 10 patients with posterior fossa tumors (4 male, mean age: 63.5), affecting either the trigeminal nerve (CN V) or the facial/vestibular complex (CN VII/VIII), were employed. Three distinct reconstruction methods [two tensor-based methods: single diffusion tensor tractography (SDT) (3D Slicer), eXtended streamline tractography (XST), and one fiber orientation distribution (FOD)-based method: streamline tractography using constrained spherical deconvolution (CSD)-derived estimates (MRtrix3)], were compared to determine which of these was best suited for use in a neurosurgical setting in terms of processing speed, anatomical accuracy, and accurate depiction of the relationship between the tumor and affected CN. Computation of the tensor map was faster when compared to the implementation of CSD to provide estimates of FOD. Both XST and CSD-based reconstruction methods tended to give more detailed representations of the projections of CN V and CN VII/VIII compared to SDT. These reconstruction methods were able to more accurately delineate the course of CN V and CN VII/VIII, differentiate CN V from the cerebellar peduncle, and delineate compression of CN VII/VIII in situations where SDT could not. However, CSD-based reconstruction methods tended to generate more invalid streamlines. XST offers the best combination of anatomical accuracy and speed of reconstruction of cranial nerves within this patient population. Given the possible anatomical limitations of single tensor models, supplementation with more advanced tensor-based reconstruction methods might be beneficial.

A Connectomic Analysis of the Human Basal Ganglia Network

The current model of basal ganglia circuits has been introduced almost two decades ago and has settled the basis for our understanding of basal ganglia physiology and movement disorders. Although many questions are yet to be answered, several efforts have been recently made to shed new light on basal ganglia function. The traditional concept of “direct” and “indirect” pathways, obtained from axonal tracing studies in non-human primates and post-mortem fiber dissection in the human brain, still retains a remarkable appeal but is somehow obsolete. Therefore, a better comprehension of human structural basal ganglia connectivity in vivo, in humans, is of uttermost importance given the involvement of these deep brain structures in many motor and non-motor functions as well as in the pathophysiology of several movement disorders. By using diffusion magnetic resonance imaging and tractography, we have recently challenged the traditional model of basal ganglia network by showing the possible existence, in the human brain, of cortico-pallidal, cortico-nigral projections, which could be mono- or polysynaptic, and an extensive subcortical network connecting the cerebellum and basal ganglia. Herein, we aimed at reconstructing the basal ganglia connectome providing a quantitative connectivity analysis of the reconstructed pathways. The present findings reinforce the idea of an intricate, not yet unraveled, network involving the cerebral cortex, basal ganglia, and cerebellum. Our findings may pave the way for a more comprehensive and holistic pathophysiological model of basal ganglia circuits.

Lipofibromatous hamartoma of the median nerve: 3T MRI evaluation by constrained spherical deconvolution analysis

In this study we described a case of lipofibromatous hamartoma involving the median nerve. We adopted diffusion tensor imaging and constrained spherical deconvolution-based tractography to reconstruct the affected median nerve. Moreover, we extracted diffusion-based parameters reflecting axonal integrity loss of median nerve fibres. Our data showed that constrained spherical deconvolution-based tractography outperformed the diffusion tensor imaging-based method, allowing the detection of the entire median nerve, including its branches, thus offering a robust method to investigate the involvement of the median nerve in pathological conditions. All clinical and technical implications are extensively described.

Improvement in White Matter Tract Reconstruction with Constrained Spherical Deconvolution and Track Density Mapping in Low Angular Resolution Data: A Pediatric Study and Literature Review

INTRODUCTION

Diffusion-weighted magnetic resonance imaging (DW-MRI) allows noninvasive investigation of brain structure in vivo. Diffusion tensor imaging (DTI) is a frequently used application of DW-MRI that assumes a single main diffusion direction per voxel, and is therefore not well suited for reconstructing crossing fiber tracts. Among the solutions developed to overcome this problem, constrained spherical deconvolution with probabilistic tractography (CSD-PT) has provided superior quality results in clinical settings on adult subjects; however, it requires particular acquisition parameters and long sequences, which may limit clinical usage in the pediatric age group. The aim of this work was to compare the results of DTI with those of track density imaging (TDI) maps and CSD-PT on data from neonates and children, acquired with low angular resolution and low b-value diffusion sequences commonly used in pediatric clinical MRI examinations.

MATERIALS AND METHODS

We analyzed DW-MRI studies of 50 children (eight neonates aged 3-28 days, 20 infants aged 1-8 months, and 22 children aged 2-17 years) acquired on a 1.5 T Philips scanner using 34 gradient directions and a b-value of 1,000 s/mm². Other sequence parameters included 60 axial slices; acquisition matrix, 128 × 128; average scan time, 5:34 min; voxel size, 1.75 mm × 1.75 mm × 2 mm; one b = 0 image. For each subject, we computed principal eigenvector (EV) maps and directionally encoded color TDI maps (DEC-TDI maps) from whole-brain tractograms obtained with CSD-PT; the cerebellar-thalamic, corticopontocerebellar, and corticospinal tracts were reconstructed using both CSD-PT and DTI. Results were compared by two neuroradiologists using a 5-point qualitative score.

RESULTS

The DEC-TDI maps obtained presented higher anatomical detail than EV maps, as assessed by visual inspection. In all subjects, white matter (WM) tracts were successfully reconstructed using both tractography methodologies. The mean qualitative scores of all tracts obtained with CSD-PT were significantly higher than those obtained with DTI (p-value < 0.05 for all comparisons).

CONCLUSION

CSD-PT can be successfully applied to DW-MRI studies acquired at 1.5 T with acquisition parameters adapted for pediatric subjects, thus providing TDI maps with greater anatomical detail. This methodology yields satisfactory results for clinical purposes in the pediatric age group.

Probabilistic versus deterministic tractography for delineation of the cortico-subthalamic hyperdirect pathway in patients with Parkinson disease selected for deep brain stimulation

OBJECTIVE Diffusion-weighted MRI (DWI) and tractography allows noninvasive mapping of the structural connections of the brain, and may provide important information for neurosurgical planning. The hyperdirect pathway, connecting the subthalamic nucleus (STN) with the motor cortex, is assumed to play a key role in mediating the effects of deep brain stimulation (DBS), which is an effective but poorly understood treatment for Parkinson disease. This study aimed to apply recent methodological advances in DWI acquisition and analysis to the delineation of the hyperdirect pathway in patients with Parkinson disease selected for surgery. METHODS High spatial and angular resolution DWI data were acquired preoperatively from 5 patients with Parkinson disease undergoing DBS. The authors compared the delineated hyperdirect pathways and associated STN target maps generated by 2 different tractography methods: a tensor-based deterministic method, typically available in clinical settings, and an advanced probabilistic method based on constrained spherical deconvolution. In addition, 10 high-resolution data sets with the same scanning parameters were acquired from a healthy control participant to assess the robustness of the tractography results. RESULTS Both tractography approaches identified connections between the ipsilateral motor cortex and the STN. However, the 2 methods provided substantially different target regions in the STN, with the target center of gravity differing by > 1.4 mm on average. The probabilistic method (based on constrained spherical deconvolution) plausibly reconstructed a continuous set of connections from the motor cortex, terminating in the dorsolateral region of the STN. In contrast, the tensor-based method reconstructed a comparatively sparser and more variable subset of connections. Furthermore, across the control scans, the probabilistic method identified considerably more consistent targeting regions within the STN compared with the deterministic tensor-based method, which demonstrated a 1.9-2.4 times higher variation. CONCLUSIONS These data provide a strong impetus for the use of a robust probabilistic tractography framework based on constrained spherical deconvolution, or similar advanced DWI models, in clinical settings. The inherent limitations and demonstrated inaccuracy of the tensor-based method leave it questionable for use in high-precision stereotactic DBS surgery. The authors have also described a straightforward method for importing tractography-derived information into any clinical neuronavigation system, based on the generation of track-density images.

Isotropic non-white matter partial volume effects in constrained spherical deconvolution

Diffusion-weighted (DW) magnetic resonance imaging (MRI) is a non-invasive imaging method, which can be used to investigate neural tracts in the white matter (WM) of the brain. Significant partial volume effects (PVEs) are present in the DW signal due to relatively large voxel sizes. These PVEs can be caused by both non-WM tissue, such as gray matter (GM) and cerebrospinal fluid (CSF), and by multiple non-parallel WM fiber populations. High angular resolution diffusion imaging (HARDI) methods have been developed to correctly characterize complex WM fiber configurations, but to date, many of the HARDI methods do not account for non-WM PVEs. In this work, we investigated the isotropic PVEs caused by non-WM tissue in WM voxels on fiber orientations extracted with constrained spherical deconvolution (CSD). Experiments were performed on simulated and real DW-MRI data. In particular, simulations were performed to demonstrate the effects of varying the diffusion weightings, signal-to-noise ratios (SNRs), fiber configurations, and tissue fractions. Our results show that the presence of non-WM tissue signal causes a decrease in the precision of the detected fiber orientations and an increase in the detection of false peaks in CSD. We estimated 35-50% of WM voxels to be affected by non-WM PVEs. For HARDI sequences, which typically have a relatively high degree of diffusion weighting, these adverse effects are most pronounced in voxels with GM PVEs. The non-WM PVEs become severe with 50% GM volume for maximum spherical harmonics orders of 8 and below, and already with 25% GM volume for higher orders. In addition, a low diffusion weighting or SNR increases the effects. The non-WM PVEs may cause problems in connectomics, where reliable fiber tracking at the WM-GM interface is especially important. We suggest acquiring data with high diffusion-weighting 2500-3000 s/mm(2), reasonable SNR (~30) and using lower SH orders in GM contaminated regions to minimize the non-WM PVEs in CSD.

Abnormal functional connectivity during visuospatial processing is associated with disrupted organisation of white matter in autism

Disruption of structural and functional neural connectivity has been widely reported in Autism Spectrum Disorder (ASD) but there is a striking lack of research attempting to integrate analysis of functional and structural connectivity in the same study population, an approach that may provide key insights into the specific neurobiological underpinnings of altered functional connectivity in autism. The aims of this study were (1) to determine whether functional connectivity abnormalities were associated with structural abnormalities of white matter (WM) in ASD and (2) to examine the relationships between aberrant neural connectivity and behavior in ASD. Twenty-two individuals with ASD and 22 age, IQ-matched controls completed a high-angular-resolution diffusion MRI scan. Structural connectivity was analysed using constrained spherical deconvolution (CSD) based tractography. Regions for tractography were generated from the results of a previous study, in which 10 pairs of brain regions showed abnormal functional connectivity during visuospatial processing in ASD. WM tracts directly connected 5 of the 10 region pairs that showed abnormal functional connectivity; linking a region in the left occipital lobe (left BA19) and five paired regions: left caudate head, left caudate body, left uncus, left thalamus, and left cuneus. Measures of WM microstructural organization were extracted from these tracts. Fractional anisotropy (FA) reductions in the ASD group relative to controls were significant for WM connecting left BA19 to left caudate head and left BA19 to left thalamus. Using a multimodal imaging approach, this study has revealed aberrant WM microstructure in tracts that directly connect brain regions that are abnormally functionally connected in ASD. These results provide novel evidence to suggest that structural brain pathology may contribute (1) to abnormal functional connectivity and (2) to atypical visuospatial processing in ASD.

Auditory tracts identified with combined fMRI and diffusion tractography

The auditory tracts in the human brain connect the inferior colliculus (IC) and medial geniculate body (MGB) to various components of the auditory cortex (AC). While in non-human primates and in humans, the auditory system is differentiated in core, belt and parabelt areas, the correspondence between these areas and anatomical landmarks on the human superior temporal gyri is not straightforward, and at present not completely understood. However it is not controversial that there is a hierarchical organization of auditory stimuli processing in the auditory system. The aims of this study were to demonstrate that it is possible to non-invasively and robustly identify auditory projections between the auditory thalamus/brainstem and different functional levels of auditory analysis in the cortex of human subjects in vivo combining functional magnetic resonance imaging (fMRI) with diffusion MRI, and to investigate the possibility of differentiating between different components of the auditory pathways (e.g. projections to areas responsible for sound, pitch and melody processing). We hypothesized that the major limitation in the identification of the auditory pathways is the known problem of crossing fibres and addressed this issue acquiring DTI with b-values higher than commonly used and adopting a multi-fibre ball-and-stick analysis model combined with probabilistic tractography. Fourteen healthy subjects were studied. Auditory areas were localized functionally using an established hierarchical pitch processing fMRI paradigm. Together fMRI and diffusion MRI allowed the successful identification of tracts connecting IC with AC in 64 to 86% of hemispheres and left sound areas with homologous areas in the right hemisphere in 86% of hemispheres. The identified tracts corresponded closely with a three-dimensional stereotaxic atlas based on postmortem data. The findings have both neuroscientific and clinical implications for delineation of the human auditory system in vivo.

Cortical and subcortical connections of the human claustrum revealed in vivo by constrained spherical deconvolution tractography

The claustrum is a thin layer of gray matter that is at the center of an active scientific debate. Recently, Constrained Spherical Deconvolution (CSD) tractography has proved to be an extraordinary tool allowing to track white matter fibers from cortex to cortical and subcortical targets with subvoxel resolution. The aim of this study was to evaluate claustral connectivity in the human brain. Ten normal brains were analyzed by using the High Angular Resolution Diffusion Imaging CSD-based technique. Tractography revealed 4 groups of white matter fibers connecting the claustrum with the brain cortex: Anterior, posterior, superior, and lateral. The anterior and posterior cortico-claustral tracts connected the claustrum to prefrontal cortex and visual areas. The superior tract linked the claustrum with sensory-motor areas, while the lateral pathway connected the claustrum to the auditory cortex. In addition, we demonstrated a claustral medial pathway connecting the claustrum with the basal ganglia, specifically with caudate nucleus, putamen, and globus pallidus. An interesting and exciting new finding was the demonstration of a bilateral connection between claustrum and contralateral cortical areas and a well-represented interclaustral communication with interconnection bundles interspersed within the bulk of the trunk of the corpus callosum. The physiological and pathophysiological relevance of these findings are discussed.

White matter and visuospatial processing in autism: a constrained spherical deconvolution tractography study

Autism spectrum disorders (ASDs) are associated with a marked disturbance of neural functional connectivity, which may arise from disrupted organization of white matter. The aim of this study was to use constrained spherical deconvolution (CSD)-based tractography to isolate and characterize major intrahemispheric white matter tracts that are important in visuospatial processing. CSD-based tractography avoids a number of critical confounds that are associated with diffusion tensor tractography, and to our knowledge, this is the first time that this advanced diffusion tractography method has been used in autism research. Twenty-five participants with ASD and aged 25, intelligence quotient-matched controls completed a high angular resolution diffusion imaging scan. The inferior fronto-occipital fasciculus (IFOF) and arcuate fasciculus were isolated using CSD-based tractography. Quantitative diffusion measures of white matter microstructural organization were compared between groups and associated with visuospatial processing performance. Significant alteration of white matter organization was present in the right IFOF in individuals with ASD. In addition, poorer visuospatial processing was associated in individuals with ASD with disrupted white matter in the right IFOF. Using a novel, advanced tractography method to isolate major intrahemispheric white matter tracts in autism, this research has demonstrated that there are significant alterations in the microstructural organization of white matter in the right IFOF in ASD. This alteration was associated with poorer visuospatial processing performance in the ASD group. This study provides an insight into structural brain abnormalities that may influence atypical visuospatial processing in autism.

Reliability of two clinically relevant fiber pathways reconstructed with constrained spherical deconvolution

The single diffusion tensor model is inadequate for the reconstruction of fiber pathways in brain regions with multiple fiber orientations. To overcome this limitation, constrained spherical deconvolution has been proposed. A high reliability of constrained spherical deconvolution is, however, a pre-requisite for its use in clinical applications. Reliability of reconstructed fiber pathways can be assessed in terms of architectural (addressing their spatial configuration) and microstructural (addressing diffusion-derived measures along the fibers) reproducibility. We assess the reliability for two clinically relevant fiber pathways: the corticospinal tract and arcuate fasciculus. The fiber pathways were reconstructed using constrained spherical deconvolution in 11 healthy subjects who were scanned on three occasions. Coefficients of variations of diffusion-derived measures were used to assess the microstructural reproducibility. Image correlation and fiber overlap were used to assess the architectural reproducibility. The mean correlation between sessions was 72% for both the corticospinal tract and arcuate fasciculus. The mean overlap between sessions was 63% for the corticospinal tract and 58% for the arcuate fasciculus. Coefficients of variations of diffusion-derived measures showed very low variation (all measures <3.1%). These results are comparable with reliability results based on the diffusion tensor model, which is commonly used in clinical settings. The reliability results found here are, therefore, promising to further investigate the use of constrained spherical deconvolution in clinical practice.

Investigating the prevalence of complex fiber configurations in white matter tissue with diffusion magnetic resonance imaging

It has long been recognized that the diffusion tensor model is inappropriate to characterize complex fiber architecture, causing tensor-derived measures such as the primary eigenvector and fractional anisotropy to be unreliable or misleading in these regions. There is however still debate about the impact of this problem in practice. A recent study using a Bayesian automatic relevance detection (ARD) multicompartment model suggested that a third of white matter (WM) voxels contain crossing fibers, a value that, whilst already significant, is likely to be an underestimate. The aim of this study is to provide more robust estimates of the proportion of affected voxels, the number of fiber orientations within each WM voxel, and the impact on tensor-derived analyses, using large, high-quality diffusion-weighted data sets, with reconstruction parameters optimized specifically for this task. Two reconstruction algorithms were used: constrained spherical deconvolution (CSD), and the ARD method used in the previous study. We estimate the proportion of WM voxels containing crossing fibers to be ~90% (using CSD) and 63% (using ARD). Both these values are much higher than previously reported, strongly suggesting that the diffusion tensor model is inadequate in the vast majority of WM regions. This has serious implications for downstream processing applications that depend on this model, particularly tractography, and the interpretation of anisotropy and radial/axial diffusivity measures.