Meyer's loop tractography for image-guided surgery depends on imaging protocol and hardware

Introduction Surgical resection is an effective treatment for temporal lobe epilepsy but can result in visual field defects. This could be minimized if surgeons knew the exact location of the anterior part of the optic radiation (OR), the Meyer's loop. To this end, there is increasing prevalence of image-guided surgery using diffusion MRI tractography. Despite considerable effort in developing analysis methods, a wide discrepancy in Meyer's loop reconstructions is observed in the literature. Moreover, the impact of differences in image acquisition on Meyer's loop tractography remains unclear. Here, while employing the same state-of-the-art analysis protocol, we explored the extent to which variance in data acquisition leads to variance in OR reconstruction. Methods Diffusion MRI data were acquired for the same thirteen healthy subjects using standard and state-of-the-art protocols on three scanners with different maximum gradient amplitudes (MGA): Siemens Connectom (MGA = 300 mT/m); Siemens Prisma (MGA = 80 mT/m) and GE Excite-HD (MGA = 40 mT/m). Meyer's loop was reconstructed on all subjects and its distance to the temporal pole (ML-TP) was compared across protocols. Results A significant effect of data acquisition on the ML-TP distance was observed between protocols (p < .01 to 0.0001). The biggest inter-acquisition discrepancy for the same subject across different protocols was 16.5 mm (mean: 9.4 mm, range: 3.7–16.5 mm). Conclusion We showed that variance in data acquisition leads to substantive variance in OR tractography. This has direct implications for neurosurgical planning, where part of the OR is at risk due to an under-estimation of its location using conventional acquisition protocols

Scanner invariant representations for diffusion MRI harmonization

Purpose In the present work, we describe the correction of diffusion‐weighted MRI for site and scanner biases using a novel method based on invariant representation. Theory and Methods Pooled imaging data from multiple sources are subject to variation between the sources. Correcting for these biases has become very important as imaging studies increase in size and multi‐site cases become more common. We propose learning an intermediate representation invariant to site/protocol variables, a technique adapted from information theory‐based algorithmic fairness; by leveraging the data processing inequality, such a representation can then be used to create an image reconstruction that is uninformative of its original source, yet still faithful to underlying structures. To implement this, we use a deep learning method based on variational auto‐encoders (VAE) to construct scanner invariant encodings of the imaging data. Results To evaluate our method, we use training data from the 2018 MICCAI Computational Diffusion MRI (CDMRI) Challenge Harmonization dataset. Our proposed method shows improvements on independent test data relative to a recently published baseline method on each subtask, mapping data from three different scanning contexts to and from one separate target scanning context. Conclusions As imaging studies continue to grow, the use of pooled multi‐site imaging will similarly increase. Invariant representation presents a strong candidate for the harmonization of these data

Automated characterization of noise distributions in diffusion MRI data

Knowledge of the noise distribution in diffusion MRI is the centerpiece to quantify uncertainties arising from the acquisition process. Accurate estimation beyond textbook distributions often requires information about the acquisition process, which is usually not available. We introduce two new automated methods using the moments and maximum likelihood equations of the Gamma distribution to estimate all unknown parameters using only the magnitude data. A rejection step is used to make the framework automatic and robust to artifacts. Simulations were created for two diffusion weightings with parallel imaging. Furthermore, MRI data of a water phantom with different combinations of parallel imaging were acquired. Finally, experiments on freely available datasets are used to assess reproducibility when limited information about the acquisition protocol is available. Additionally, we demonstrated the applicability of the proposed methods for a bias correction and denoising task on an in vivo dataset. A generalized version of the bias correction framework for non integer degrees of freedom is also introduced. The proposed framework is compared with three other algorithms with datasets from three vendors, employing different reconstruction methods. Simulations showed that assuming a Rician distribution can lead to misestimation of the noise distribution in parallel imaging. Results showed that signal leakage in multiband can also lead to a misestimation of the noise distribution. Repeated acquisitions of in vivo datasets show that the estimated parameters are stable and have lower variability than compared methods. Results show that the proposed methods reduce the appearance of noise at high b-value. The proposed algorithms herein can estimate both parameters of the noise distribution automatically, are robust to signal leakage artifacts and perform best when used on acquired noise maps.Comment: v3: Peer reviewed version v2: Manuscript as submitted to Medical image analysis v1: Manuscript as submitted to Magnetic resonance in medicin

The importance of correcting for signal drift in diffusion MRI

Purpose To investigate previously unreported effects of signal drift as a result of temporal scanner instability on diffusion MRI data analysis and to propose a method to correct this signal drift. Methods We investigated the signal magnitude of non-diffusion-weighted EPI volumes in a series of diffusion-weighted imaging experiments to determine whether signal magnitude changes over time. Different scan protocols and scanners from multiple vendors were used to verify this on phantom data, and the effects on diffusion kurtosis tensor estimation in phantom and in vivo data were quantified. Scalar metrics (eigenvalues, fractional anisotropy, mean diffusivity, mean kurtosis) and directional information (first eigenvectors and tractography) were investigated. Results Signal drift, a global signal decrease with subsequently acquired images in the scan, was observed in phantom data on all three scanners, with varying magnitudes up to 5% in a 15-min scan. The signal drift has a noticeable effect on the estimation of diffusion parameters. All investigated quantitative parameters as well as tractography were affected by this artifactual signal decrease during the scan. Conclusion By interspersing the non-diffusion-weighted images throughout the session, the signal decrease can be estimated and compensated for before data analysis; minimizing the detrimental effects on subsequent MRI analyses. Magn Reson Med 77:285–299, 2017. © 2016 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine

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 0 – 6000 s/mm2. 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=6000 s/mm2, median R2 = .34; b=4000 s/mm2, median R2 = .29; b=2400 s/mm2, median R2 = .21; b=1200 s/mm2, median R2 = .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

Transferring principles of solid-state and Laplace NMR to the field of in vivo brain MRI

Magnetic resonance imaging (MRI) is the primary method for non-invasive investigations of the human brain in health, disease, and development, but yields data that are difficult to interpret whenever the millimeter scale voxels contain multiple microscopic tissue environments with different chemical and structural properties. We propose a novel MRI framework to quantify the microscopic heterogeneity of the living human brain as spatially resolved five-dimensional relaxation-diffusion distributions by augmenting a conventional diffusion-weighted imaging sequence with signal encoding principles from multidimensional solid-state nuclear magnetic resonance (NMR) spectroscopy, relaxation-diffusion correlation methods from Laplace NMR of porous media, and Monte Carlo data inversion. The high dimensionality of the distribution space allows resolution of multiple microscopic environments within each heterogeneous voxel as well as their individual characterization with novel statistical measures that combine the chemical sensitivity of the relaxation rates with the link between microstructure and the anisotropic diffusivity of tissue water. The proposed framework is demonstrated on a healthy volunteer using both exhaustive and clinically viable acquisition protocols

The effect of gradient nonlinearities on fiber orientation estimates from spherical deconvolution of diffusion MRI 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

Magnetic resonance imaging of T2 - and diffusion snisotropy using a tiltable receive coil

The anisotropic microstructure of white matter is reflected in various MRI contrasts. Transverse relaxation rates can be probed as a function of fibre-orientation with respect to the main magnetic field, while diffusion properties are probed as a function of fibre-orientation with respect to an encoding gradient. While the latter is easy to obtain by varying the orientation of the gradient, as the magnetic field is fixed, obtaining the former requires re-orienting the head. In this work we deployed a tiltable RF-coil to study T2 - and diffusional anisotropy of the brain white matter simultaneously in diffusion- T2 correlation experiments

Structural magnetic resonance imaging in dystonia: A systematic review of methodological approaches and findings

Background and purpose: Structural magnetic resonance techniques have been widely applied in neurological disorders to better understand tissue changes, probing characteristics such as volume, iron deposition and diffusion. Dystonia is a hyperkinetic movement disorder, resulting in abnormal postures and pain. Its pathophysiology is poorly understood, with normal routine clinical imaging in idiopathic forms. More advanced tools provide an opportunity to identify smaller scale structural changes which may underpin pathophysiology. This review aims to provide an overview of methodological approaches undertaken in structural brain imaging of dystonia cohorts, and to identify commonly identified pathways, networks or regions that are implicated in pathogenesis. Methods: Structural magnetic resonance imaging studies of idiopathic and genetic forms of dystonia were systematically reviewed. Adhering to strict inclusion and exclusion criteria, PubMed and Embase databases were searched up to January 2022, with studies reviewed for methodological quality and key findings. Results: Seventy‐seven studies were included, involving 1945 participants. The majority of studies employed diffusion tensor imaging (DTI) (n = 45) or volumetric analyses (n = 37), with frequently implicated areas of abnormality in the brainstem, cerebellum, basal ganglia and sensorimotor cortex and their interconnecting white matter pathways. Genotypic and motor phenotypic variation emerged, for example fewer cerebello‐thalamic tractography streamlines in genetic forms than idiopathic and higher grey matter volumes in task‐specific than non‐task‐specific dystonias. Discussion: Work to date suggests microstructural brain changes in those diagnosed with dystonia, although the underlying nature of these changes remains undetermined. Employment of techniques such as multiple diffusion weightings or multi‐exponential relaxometry has the potential to enhance understanding of these differences