Diffusion-MRI in neurodegenerative disorders

Limbic-predominant age-related TDP-43 encephalopathy neuropathological change (LATE-NC) is associated with abnormalities in white matter structural integrity and connectivity: An ex-vivo diffusion MRI and pathology investigation

Limbic predominant age-related TDP-43 encephalopathy neuropathological change (LATE-NC) is common in older adults and is associated with neurodegeneration, cognitive decline and dementia. In this MRI and pathology investigation we tested the hypothesis that LATE-NC is associated with abnormalities in white matter structural integrity and connectivity of a network of brain regions typically harboring TDP-43 inclusions in LATE, referred to here as the "LATE-NC network". Ex-vivo diffusion MRI and detailed neuropathological data were collected on 184 community-based older adults. Linear regression revealed an independent association of higher LATE-NC stage with lower diffusion anisotropy in a set of white matter connections forming a pattern of connectivity that is consistent with the stereotypical spread of this pathology in the brain. Graph theory analysis revealed an association of higher LATE-NC stage with weaker integration and segregation in the LATE-NC network. Abnormalities were significant in stage 3, suggesting that they are detectable in later stages of the disease. Finally, LATE-NC network abnormalities were associated with faster cognitive decline, specifically in episodic and semantic memory.

Dynamic parallel transmit diffusion MRI at 7T

Diffusion MRI (dMRI) is inherently limited by SNR. Scanning at 7 T increases intrinsic SNR but 7 T MRI scans suffer from regions of signal dropout, especially in the temporal lobes and cerebellum. We applied dynamic parallel transmit (pTx) to allow whole-brain 7 T dMRI and compared with circularly polarized (CP) pulses in 6 subjects. Subject-specific 2-spoke dynamic pTx pulses were designed offline for 8 slabs covering the brain. We used vendor-provided B0 and B1+ mapping. Spokes positions were set using the Fourier difference approach, and RF coefficients optimized with a Jacobi-matrix high-flip-angle optimizer. Diffusion data were analyzed with FSL. Comparing whole-brain averages for pTx against CP scans: mean flip angle error improved by 15% for excitation (2-spoke-VERSE 15.7° vs CP 18.4°, P = 0.012) and improved by 14% for refocusing (2-spoke-VERSE 39.7° vs CP 46.2°, P = 0.008). Computed spin-echo signal standard deviation improved by 14% (2-spoke-VERSE 0.185 vs 0.214 CP, P = 0.025). Temporal SNR increased by 5.4% (2-spoke-VERSE 8.47 vs CP 8.04, P = 0.004) especially in the inferior temporal lobes. Diffusion fitting uncertainty decreased by 6.2% for first fibers (2-spoke VERSE 0.0655 vs CP 0.0703, P < 0.001) and 1.3% for second fibers (2-spoke VERSE 0.139 vs CP 0.141, P = 0.01). In conclusion, dynamic parallel transmit improves the uniformity of 7 T diffusion-weighted imaging. In future, less restrictive SAR limits for parallel transmit scans are expected to allow further improvements.

Integration of diffusion tensor imaging parameters with mesh morphing for in-depth analysis of brain white matter fibre tracts

Averaging is commonly used for data reduction/aggregation to analyse high-dimensional MRI data, but this often leads to information loss. To address this issue, we developed a novel technique that integrates diffusion tensor metrics along the whole volume of the fibre bundle using a 3D mesh-morphing technique coupled with principal component analysis for delineating case and control groups. Brain diffusion tensor MRI scans of high school rugby union players (n = 30, age 16-18) were acquired on a 3 T MRI before and after the sports season. A non-contact sport athlete cohort with matching demographics (n = 12) was also scanned. The utility of the new method in detecting differences in diffusion tensor metrics of the right corticospinal tract between contact and non-contact sport athletes was explored. The first step was to run automated tractography on each subject's native space. A template model of the right corticospinal tract was generated and morphed into each subject's native shape and space, matching individual geometry and diffusion metric distributions with minimal information loss. The common dimension of the 20 480 diffusion metrics allowed further data aggregation using principal component analysis to cluster the case and control groups as well as visualization of diffusion metric statistics (mean, ±2 SD). Our approach of analysing the whole volume of white matter tracts led to a clear delineation between the rugby and control cohort, which was not possible with the traditional averaging method. Moreover, our approach accounts for the individual subject's variations in diffusion tensor metrics to visualize group differences in quantitative MR data. This approach may benefit future prediction models based on other quantitative MRI methods.

Practical considerations of diffusion-weighted MRS with ultra-strong diffusion gradients

Introduction

Diffusion-weighted magnetic resonance spectroscopy (DW-MRS) offers improved cellular specificity to microstructure-compared to water-based methods alone-but spatial resolution and SNR is severely reduced and slow-diffusing metabolites necessitate higher b-values to accurately characterize their diffusion properties. Ultra-strong gradients allow access to higher b-values per-unit time, higher SNR for a given b-value, and shorter diffusion times, but introduce additional challenges such as eddy-current artefacts, gradient non-uniformity, and mechanical vibrations.

Methods

In this work, we present initial DW-MRS data acquired on a 3T Siemens Connectom scanner equipped with ultra-strong (300 mT/m) gradients. We explore the practical issues associated with this manner of acquisition, the steps that may be taken to mitigate their impact on the data, and the potential benefits of ultra-strong gradients for DW-MRS. An in-house DW-PRESS sequence and data processing pipeline were developed to mitigate the impact of these confounds. The interaction of TE, b-value, and maximum gradient amplitude was investigated using simulations and pilot data, whereby maximum gradient amplitude was restricted. Furthermore, two DW-MRS voxels in grey and white matter were acquired using ultra-strong gradients and high b-values.

Results

Simulations suggest T2-based SNR gains that are experimentally confirmed. Ultra-strong gradient acquisitions exhibit similar artefact profiles to those of lower gradient amplitude, suggesting adequate performance of artefact mitigation strategies. Gradient field non-uniformity influenced ADC estimates by up to 4% when left uncorrected. ADC and Kurtosis estimates for tNAA, tCho, and tCr align with previously published literature.

Discussion

In conclusion, we successfully implemented acquisition and data processing strategies for ultra-strong gradient DW-MRS and results indicate that confounding effects of the strong gradient system can be ameliorated, while achieving shorter diffusion times and improved metabolite SNR.

Automated motion artifact detection in early pediatric diffusion MRI using a convolutional neural network

Diffusion MRI (dMRI) is a widely used method to investigate the microstructure of the brain. Quality control (QC) of dMRI data is an important processing step that is performed prior to analysis using models such as diffusion tensor imaging (DTI) or neurite orientation dispersion and density imaging (NODDI). When processing dMRI data from infants and young children, where intra-scan motion is common, the identification and removal of motion artifacts is of the utmost importance. Manual QC of dMRI data is (1) time-consuming due to the large number of diffusion directions, (2) expensive, and (3) prone to subjective errors and observer variability. Prior techniques for automated dMRI QC have mostly been limited to adults or school-age children. Here, we propose a deep learning-based motion artifact detection tool for dMRI data acquired from infants and toddlers. The proposed framework uses a simple three-dimensional convolutional neural network (3DCNN) trained and tested on an early pediatric dataset of 2,276 dMRI volumes from 121 exams acquired at 1 month and 24 months of age. An average classification accuracy of 95% was achieved following four-fold cross-validation. A second dataset with different acquisition parameters and ages ranging from 2-36 months (consisting of 2,349 dMRI volumes from 26 exams) was used to test network generalizability, achieving 98% classification accuracy. Finally, to demonstrate the importance of motion artifact volume removal in a dMRI processing pipeline, the dMRI data were fit to the DTI and NODDI models and the parameter maps were compared with and without motion artifact removal.

Cerebellar structural, astrocytic, and neuronal abnormalities in the SMNΔ7 mouse model of spinal muscular atrophy

Spinalmuscular atrophy (SMA) is a neuromuscular disease that affects as many as 1 in 6000 individuals at birth, making it the leading genetic cause of infant mortality. A growing number of studies indicate that SMA is a multi-system disease. The cerebellum has received little attention even though it plays an important role in motor function and widespread pathology has been reported in the cerebella of SMA patients. In this study, we assessed SMA pathology in the cerebellum using structural and diffusion magnetic resonance imaging, immunohistochemistry, and electrophysiology with the SMNΔ7 mouse model. We found a significant disproportionate loss in cerebellar volume, decrease in afferent cerebellar tracts, selective lobule-specific degeneration of Purkinje cells, abnormal lobule foliation and astrocyte integrity, and a decrease in spontaneous firing of cerebellar output neurons in the SMA mice compared to controls. Our data suggest that defects in cerebellar structure and function due to decreased survival motor neuron (SMN) levels impair the functional cerebellar output affecting motor control, and that cerebellar pathology should be addressed to achieve comprehensive treatment and therapy for SMA patients.

Viability of AMURA biomarkers from single-shell diffusion MRI in clinical studies

Diffusion Tensor Imaging (DTI) is the most employed method to assess white matter properties using quantitative parameters derived from diffusion MRI, but it presents known limitations that restrict the evaluation of complex structures. The objective of this study was to validate the reliability and robustness of complementary diffusion measures extracted with a novel approach, Apparent Measures Using Reduced Acquisitions (AMURA), with a typical diffusion MRI acquisition from a clinical context in comparison with DTI with application to clinical studies. Fifty healthy controls, 51 episodic migraine and 56 chronic migraine patients underwent single-shell diffusion MRI. Four DTI-based and eight AMURA-based parameters were compared between groups with tract-based spatial statistics to establish reference results. On the other hand, following a region-based analysis, the measures were assessed for multiple subsamples with diverse reduced sample sizes and their stability was evaluated with the coefficient of quartile variation. To assess the discrimination power of the diffusion measures, we repeated the statistical comparisons with a region-based analysis employing reduced sample sizes with diverse subsets, decreasing 10 subjects per group for consecutive reductions, and using 5,001 different random subsamples. For each sample size, the stability of the diffusion descriptors was evaluated with the coefficient of quartile variation. AMURA measures showed a greater number of statistically significant differences in the reference comparisons between episodic migraine patients and controls compared to DTI. In contrast, a higher number of differences was found with DTI parameters compared to AMURA in the comparisons between both migraine groups. Regarding the assessments reducing the sample size, the AMURA parameters showed a more stable behavior than DTI, showing a lower decrease for each reduced sample size or a higher number of regions with significant differences. However, most AMURA parameters showed lower stability in relation to higher coefficient of quartile variation values than the DTI descriptors, although two AMURA measures showed similar values to DTI. For the synthetic signals, there were AMURA measures with similar quantification to DTI, while other showed similar behavior. These findings suggest that AMURA presents favorable characteristics to identify differences of specific microstructural properties between clinical groups in regions with complex fiber architecture and lower dependency on the sample size or assessing technique than DTI.

Potential Pitfalls of Using Fractional Anisotropy, Axial Diffusivity, and Radial Diffusivity as Biomarkers of Cerebral White Matter Microstructure

Fractional anisotropy (FA), axial diffusivity (AD), and radial diffusivity (RD) are commonly used as MRI biomarkers of white matter microstructure in diffusion MRI studies of neurodevelopment, brain aging, and neurologic injury/disease. Some of the more frequent practices include performing voxel-wise or region-based analyses of these measures to cross-sectionally compare individuals or groups, longitudinally assess individuals or groups, and/or correlate with demographic, behavioral or clinical variables. However, it is now widely recognized that the majority of cerebral white matter voxels contain multiple fiber populations with different trajectories, which renders these metrics highly sensitive to the relative volume fractions of the various fiber populations, the microstructural integrity of each constituent fiber population, and the interaction between these factors. Many diffusion imaging experts are aware of these limitations and now generally avoid using FA, AD or RD (at least in isolation) to draw strong reverse inferences about white matter microstructure, but based on the continued application and interpretation of these metrics in the broader biomedical/neuroscience literature, it appears that this has perhaps not yet become common knowledge among diffusion imaging end-users. Therefore, this paper will briefly discuss the complex biophysical underpinnings of these measures in the context of crossing fibers, provide some intuitive “thought experiments” to highlight how conventional interpretations can lead to incorrect conclusions, and suggest that future studies refrain from using (over-interpreting) FA, AD, and RD values as standalone biomarkers of cerebral white matter microstructure.

Microstructural changes in prodromal dementia with Lewy bodies compared to normal aging: multiparametric quantitative MRI evidences

Dementia with Lewy bodies (DLB) patients show few significant macroscopic structural changes, especially at the early stages of the disease, making quantitative MRI especially interesting to explore more subtle changes that are not detectable by conventional volumetric techniques. Microstructural alterations have been reported in DLB at the dementia stage, but no study to date was conducted in prodromal patients. Here, quantitative MRI data were collected from 46 DLB prodromal patients and 20 healthy elderly subjects, who also underwent a detailed clinical examination including the Mayo Clinic Fluctuation Scale. We conducted voxel-wise between-group comparisons in diffusion tensor imaging (DTI) metrics and in R2* mapping, along with a multivariate analysis combining the two modalities. We highlighted multiple grey matter and white matter microstructural changes in DLB patients at the prodromal stage, compared to control subjects. Our multivariate analysis identified three distinct regional patterns of DTI and R2* changes (anterior, anteromedial, posterior) in DLB patients, that could reflect different neuropathological processes across brain regions. We also observed an association between R2* alterations in the thalamus, and the severity of fluctuations, in the DLB group. These preliminary findings are promising and require future investigations to better understand the biological underpinnings of microstructural alterations.

Same Brain, Different Look?-The Impact of Scanner, Sequence and Preprocessing on Diffusion Imaging Outcome Parameters

In clinical diagnostics and longitudinal studies, the reproducibility of MRI assessments is of high importance in order to detect pathological changes, but developments in MRI hard- and software often outrun extended periods of data acquisition and analysis. This could potentially introduce artefactual changes or mask pathological alterations. However, if and how changes of MRI hardware, scanning protocols or preprocessing software affect complex neuroimaging outcomes from, e.g., diffusion weighted imaging (DWI) remains largely understudied. We therefore compared DWI outcomes and artefact severity of 121 healthy participants (age range 19–54 years) who underwent two matched DWI protocols (Siemens product and Center for Magnetic Resonance Research sequence) at two sites (Siemens 3T Magnetom Verio and Skyra^fit). After different preprocessing steps, fractional anisotropy (FA) and mean diffusivity (MD) maps, obtained by tensor fitting, were processed with tract-based spatial statistics (TBSS). Inter-scanner and inter-sequence variability of skeletonised FA values reached up to 5% and differed largely in magnitude and direction across the brain. Skeletonised MD values differed up to 14% between scanners. We here demonstrate that DTI outcome measures strongly depend on imaging site and software, and that these biases vary between brain regions. These regionally inhomogeneous biases may exceed and considerably confound physiological effects such as ageing, highlighting the need to harmonise data acquisition and analysis. Future studies thus need to implement novel strategies to augment neuroimaging data reliability and replicability.

Cross-scanner reproducibility and harmonization of a diffusion MRI structural brain network: A traveling subject study of multi-b acquisition

Characterization of brain networks by diffusion MRI (dMRI) has rapidly evolved, and there are ongoing movements toward data sharing and multi-center studies. To extract meaningful information from multi-center data, methods to correct for the bias caused by scanner differences, that is, harmonization, are urgently needed. In this work, we report the cross-scanner differences in structural network analyses using data from nine traveling subjects (four males and five females, 21-49 years-old) who underwent scanning using four 3T scanners (public database available from the Brain/MINDS Beyond Human Brain MRI project (http://mriportal.umin.jp/)). The reliability and reproducibility were compared to those of data from another set of four subjects (all males, 29-42 years-old) who underwent scan-rescan (interval, 105-147 days) with the same scanner as well as scan-rescan data from the Human Connectome Project database. The results demonstrated that the reliability of the edge weights and graph theory metrics was lower for data including different scanners, compared to the scan-rescan with the same scanner. Besides, systematic differences between scanners were observed, indicating the risk of bias in comparing networks obtained from different scanners directly. We further demonstrate that it is feasible to reduce inter-scanner variabilities while preserving the inter-subject differences among healthy individuals by modeling the scanner effects at the level of network matrices, when traveling-subject data are available for calibration between scanners. The present data and results are expected to serve as a basis for developing and evaluating novel harmonization methods.

Conductance-Based Structural Brain Connectivity in Aging and Dementia

Background: Structural brain connectivity has been shown to be sensitive to the changes that the brain undergoes during Alzheimer's disease (AD) progression. Methods: In this work, we used our recently proposed structural connectivity quantification measure derived from diffusion magnetic resonance imaging, which accounts for both direct and indirect pathways, to quantify brain connectivity in dementia. We analyzed data from the second phase of Alzheimer's Disease Neuroimaging Initiative and third release in the Open Access Series of Imaging Studies data sets to derive relevant information for the study of the changes that the brain undergoes in AD. We also compared these data sets to the Human Connectome Project data set, as a reference, and eventually validated externally on two cohorts of the European DTI Study in Dementia database. Results: Our analysis shows expected trends of mean conductance with respect to age and cognitive scores, significant age prediction values in aging data, and regional effects centered among subcortical regions, and cingulate and temporal cortices. Discussion: Results indicate that the conductance measure has prediction potential, especially for age, that age and cognitive scores largely overlap, and that this measure could be used to study effects such as anticorrelation in structural connections. Impact statement This work presents a methodology and a set of analyses that open new possibilities in the study of healthy and pathological aging. The methodology used here is sensitive to direct and indirect pathways in deriving brain connectivity measures from diffusion-weighted magnetic resonance imaging, and therefore provides information that many state-of-the-art methods do not account for. As a result, this technique may provide the research community with ways to detect subtle effects of healthy aging and Alzheimer's disease.

A multi-inversion-recovery magnetic resonance fingerprinting for multi-compartment water mapping

PURPOSE

This study aimed at introducing short-T1/T2 compartment to MR fingerprinting (MRF) at 3 T. Water that is bound to myelin macromolecules have significantly shorter T1 and T2 than free water and can be distinguished from free water by multi-compartment analysis.

METHODS

We developed a new multi-inversion-recovery (mIR) water mapping-MRF based on an unbalanced steady-state coherent sequence (FISP). mIR pulses with an interval of 400 or 500 repetition times (TRs) were inserted into the conventional FISP MRF sequence. Data from our proposed mIR MRF was used to quantify different compartments, including myelin water, gray matter free water, and white matter free water, of brain water by virtue of the iterative non-negative least square (NNLS) with reweighting. Three healthy volunteers were scanned with mIR MRF on a clinical 3 T MRI.

RESULTS

Using an extended phase graph simulation, we found that our proposed mIR scheme with four IR pulses allowed differentiation between short and long T1/T2 components. For in vivo experiments, we achieved the quantification of myelin water, gray matter water, and white matter water at an image resolution of 1.17 × 1.17 × 5 mm³/pixel. As compared to the conventional MRF technique with single IR, our proposed mIR improved the detection of myelin water content. In addition, mIR MRF using spiral-in/out trajectory provided a higher signal level compared with that with spiral-out trajectory. Myelin water quantification using mIR MRF with 4 IR and 5 IR pulses were qualitatively similar. Meanwhile, 5 IR MRF showed fewer artifacts in myelin water detection.

CONCLUSION

We developed a new mIR MRF sequence for the rapid quantification of brain water compartments.

Comparison of Neurite Orientation Dispersion and Density Imaging and Two-Compartment Spherical Mean Technique Parameter Maps in Multiple Sclerosis

Background: Neurite orientation dispersion and density imaging (NODDI) and the spherical mean technique (SMT) are diffusion MRI methods providing metrics with sensitivity to similar characteristics of white matter microstructure. There has been limited comparison of changes in NODDI and SMT parameters due to multiple sclerosis (MS) pathology in clinical settings.

Purpose: To compare group-wise differences between healthy controls and MS patients in NODDI and SMT metrics, investigating associations with disability and correlations with diffusion tensor imaging (DTI) metrics.

Methods: Sixty three relapsing-remitting MS patients were compared to 28 healthy controls. NODDI and SMT metrics corresponding to intracellular volume fraction (v_in), orientation dispersion (ODI and ODE), diffusivity (D) (SMT only) and isotropic volume fraction (v_iso) (NODDI only) were calculated from diffusion MRI data, alongside DTI metrics (fractional anisotropy, FA; axial/mean/radial diffusivity, AD/MD/RD). Correlations between all pairs of MRI metrics were calculated in normal-appearing white matter (NAWM). Associations with expanded disability status scale (EDSS), controlling for age and gender, were evaluated. Patient-control differences were assessed voxel-by-voxel in MNI space controlling for age and gender at the 5% significance level, correcting for multiple comparisons. Spatial overlap of areas showing significant differences were compared using Dice coefficients.

Results: NODDI and SMT show significant associations with EDSS (standardised beta coefficient −0.34 in NAWM and −0.37 in lesions for NODDI v_in; 0.38 and −0.31 for SMT ODE and v_in in lesions; p < 0.05). Significant correlations in NAWM are observed between DTI and NODDI/SMT metrics. NODDI v_in and SMT v_in strongly correlated (r = 0.72, p < 0.05), likewise NODDI ODI and SMT ODE (r = −0.80, p < 0.05). All DTI, NODDI and SMT metrics detect widespread differences between patients and controls in NAWM (12.57% and 11.90% of MNI brain mask for SMT and NODDI v_in, Dice overlap of 0.42).

Data Conclusion: SMT and NODDI detect significant differences in white matter microstructure between MS patients and controls, concurring on the direction of these changes, providing consistent descriptors of tissue microstructure that correlate with disability and show alterations beyond focal damage. Our study suggests that NODDI and SMT may play a role in monitoring MS in clinical trials and practice.

On the dependence of quantitative diffusion-weighted imaging on scanner system characteristics and acquisition parameters: A large multicenter and multiparametric phantom study with unsupervised clustering analysis

PURPOSE

The purpose of this multicenter phantom study was to exploit an innovative approach, based on an extensive acquisition protocol and unsupervised clustering analysis, in order to assess any potential bias in apparent diffusion coefficient (ADC) estimation due to different scanner characteristics. Moreover, we aimed at assessing, for the first time, any effect of acquisition plan/phase encoding direction on ADC estimation.

METHODS

Water phantom acquisitions were carried out on 39 scanners. DWI acquisitions (b-value = 0-200-400-600-800-1000 s/mm²) with different acquisition plans (axial, coronal, sagittal) and phase encoding directions (anterior/posterior and right/left, for the axial acquisition plan), for 3 orthogonal diffusion weighting gradient directions, were performed. For each acquisition setup, ADC values were measured in-center and off-center (6 different positions), resulting in an entire dataset of 84 × 39 = 3276 ADC values. Spatial uniformity of ADC maps was assessed by means of the percentage difference between off-center and in-center ADC values (Δ).

RESULTS

No significant dependence of in-center ADC values on acquisition plan/phase encoding direction was found. Ward unsupervised clustering analysis showed 3 distinct clusters of scanners and an association between Δ-values and manufacturer/model, whereas no association between Δ-values and maximum gradient strength, slew rate or static magnetic field strength was revealed. Several acquisition setups showed significant differences among groups, indicating the introduction of different biases in ADC estimation.

CONCLUSIONS

Unsupervised clustering analysis of DWI data, obtained from several scanners using an extensive acquisition protocol, allows to reveal an association between measured ADC values and manufacturer/model of scanner, as well as to identify suboptimal DWI acquisition setups for accurate ADC estimation.

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

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

A Role of PET/MR Imaging in Dementia?

Since many years, magnetic resonance imaging (MRI) and positron emission tomography (PET) have a prominent role in neurodegenerative disorders and dementia, not only in a research setting but also in a clinical setting. For several decades, information from both modalities is combined ranging from individual visual assessments to fully integrating all images. Several tools are available to coregister images from MRI and PET and to covisualize these images. When studying neurodegenerative disorders with PET it is important to perform a partial volume correction and this can be done using the structural information obtained by MRI. With the advent of PET/MR, the question arises in how far this hybrid imaging modality is an added value compared to combining PET and MRI data from two separate modalities. One issue in PET/MR is still not yet completely settled, that is, the attenuation correction. This is of less importance for visual assessments but it can become an issue when combining data from PET/CT and PET/MR scanners in multicenter studies or when using cut-off values to classify patients. Simultaneous imaging has clearly some advantages: for the patient it is beneficial to have only one scan session instead of two but also in cases in which PET data are related to functional of physiological data acquired with MRI (such as functional MRI or arterial spin labeling). However, the most important benefit is currently the more integrated use of PET and MRI. This is also possible with separate measurements but requires more streamlining of the whole process. In that case coregistration of images is mandatory. It needs to be determined in which cases simultaneous PET/MRI leads to new insights or improved diagnosis compared to multimodal imaging using dedicated scanners.

Brain White-Matter Degeneration Due to Aging and Parkinson Disease as Revealed by Double Diffusion Encoding

Microstructure imaging by means of multidimensional diffusion encoding is increasingly applied in clinical research, with expectations that it yields a parameter that better correlates with clinical disability than current methods based on single diffusion encoding. Under the assumption that diffusion within a voxel can be well described by a collection of diffusion tensors, several parameters of this diffusion tensor distribution can be derived, including mean size, variance of sizes, orientational dispersion, and microscopic anisotropy. The information provided by multidimensional diffusion encoding also enables us to decompose the sources of the conventional fractional anisotropy and mean kurtosis. In this study, we explored the utility of the diffusion tensor distribution approach for characterizing white-matter degeneration in aging and in Parkinson disease by using double diffusion encoding. Data from 23 healthy older subjects and 27 patients with Parkinson disease were analyzed. Advanced age was associated with greater mean size and size variances, as well as smaller microscopic anisotropy. By analyzing the parameters underlying diffusion kurtosis, we found that the reductions of kurtosis in aging and Parkinson disease reported in the literature are likely driven by the reduction in microscopic anisotropy. Furthermore, microscopic anisotropy correlated with the severity of motor impairment in the patients with Parkinson disease. The present results support the use of multidimensional diffusion encoding in clinical studies and are encouraging for its future clinical implementation.

Estimation error bound for GRAPPA diffusion-weighted MRI

In recent years, diffusion weight magnetic resonance imaging (DW-MRI) has become one of the most important MRI imaging modalities. The importance of the DW-MRI grew thanks to the combination of parallel magnetic resonance imaging (pMRI) techniques with the echo-planar imaging (EPI), which minimize scan time and lead to reduced distortion, allowing the DW-MRI to become a routine clinical exam. Additionally, this has brought various new parameters that influence image quality and biomarkers used in DW-MRI. This work aims to investigate the effects of these parameters on the estimation quality, by using the Cramér-Rao bound tool, which gives analytical expressions of the lower limit on the estimation error variance of different DW-MRI variables when using the pMRI technique. In particular, these bounds will be used to study and optimize the impact of different factors of generalized autocalibrating partially parallel acquisition (GRAPPA) technique and system parameters on the estimation quality of the desired clinical metrics. Moreover, the obtained results of this study can be exploited and adapted in all human body DW-MRI clinical routines, further improving disease diagnosis, and tractography studies.

Diffusion Tensor Imaging-Based Studies at the Group-Level Applied to Animal Models of Neurodegenerative Diseases

The understanding of human and non-human microstructural brain alterations in the course of neurodegenerative diseases has substantially improved by the non-invasive magnetic resonance imaging (MRI) technique of diffusion tensor imaging (DTI). Animal models (including disease or knockout models) allow for a variety of experimental manipulations, which are not applicable to humans. Thus, the DTI approach provides a promising tool for cross-species cross-sectional and longitudinal investigations of the neurobiological targets and mechanisms of neurodegeneration. This overview with a systematic review focuses on the principles of DTI analysis as used in studies at the group level in living preclinical models of neurodegeneration. The translational aspect from in-vivo animal models toward (clinical) applications in humans is covered as well as the DTI-based research of the non-human brains' microstructure, the methodological aspects in data processing and analysis, and data interpretation at different abstraction levels. The aim of integrating DTI in multiparametric or multimodal imaging protocols will allow the interrogation of DTI data in terms of directional flow of information and may identify the microstructural underpinnings of neurodegeneration-related patterns.

NODDI in clinical research

Diffusion MRI (dMRI) has proven to be a useful imaging approach for both clinical diagnosis and research investigating the microstructures of nervous tissues, and it has helped us to better understand the neurophysiological mechanisms of many diseases. Though diffusion tensor imaging (DTI) has long been the default tool to analyze dMRI data in clinical research, acquisition with stronger diffusion weightings beyond the DTI regimen is now possible with modern clinical scanners, potentially enabling even more detailed characterization of tissue microstructures. To take advantage of such data, neurite orientation dispersion and density imaging (NODDI) has been proposed as a way to relate the dMRI signal to tissue features via biophysically inspired modeling. The number of reports demonstrating the potential clinical utility of NODDI is rapidly increasing. At the same time, the pitfalls and limitations of NODDI, and general challenges in microstructure modeling, are becoming increasingly recognized by clinicians. dMRI microstructure modeling is a rapidly evolving field with great promise, where people from different scientific backgrounds, such as physics, medicine, biology, neuroscience, and statistics, are collaborating to build novel tools that contribute to improving human healthcare. Here, we review the applications of NODDI in clinical research and discuss future perspectives for investigations toward the implementation of dMRI microstructure imaging in clinical practice.

Towards accurate and unbiased imaging-based differentiation of Parkinson's disease, progressive supranuclear palsy and corticobasal syndrome

The early and accurate differential diagnosis of parkinsonian disorders is still a significant challenge for clinicians. In recent years, a number of studies have used magnetic resonance imaging data combined with machine learning and statistical classifiers to successfully differentiate between different forms of Parkinsonism. However, several questions and methodological issues remain, to minimize bias and artefact-driven classification. In this study, we compared different approaches for feature selection, as well as different magnetic resonance imaging modalities, with well-matched patient groups and tightly controlling for data quality issues related to patient motion. Our sample was drawn from a cohort of 69 healthy controls, and patients with idiopathic Parkinson’s disease (n = 35), progressive supranuclear palsy Richardson’s syndrome (n = 52) and corticobasal syndrome (n = 36). Participants underwent standardized T1-weighted and diffusion-weighted magnetic resonance imaging. Strict data quality control and group matching reduced the control and patient numbers to 43, 32, 33 and 26, respectively. We compared two different methods for feature selection and dimensionality reduction: whole-brain principal components analysis, and an anatomical region-of-interest based approach. In both cases, support vector machines were used to construct a statistical model for pairwise classification of healthy controls and patients. The accuracy of each model was estimated using a leave-two-out cross-validation approach, as well as an independent validation using a different set of subjects. Our cross-validation results suggest that using principal components analysis for feature extraction provides higher classification accuracies when compared to a region-of-interest based approach. However, the differences between the two feature extraction methods were significantly reduced when an independent sample was used for validation, suggesting that the principal components analysis approach may be more vulnerable to overfitting with cross-validation. Both T1-weighted and diffusion magnetic resonance imaging data could be used to successfully differentiate between subject groups, with neither modality outperforming the other across all pairwise comparisons in the cross-validation analysis. However, features obtained from diffusion magnetic resonance imaging data resulted in significantly higher classification accuracies when an independent validation cohort was used. Overall, our results support the use of statistical classification approaches for differential diagnosis of parkinsonian disorders. However, classification accuracy can be affected by group size, age, sex and movement artefacts. With appropriate controls and out-of-sample cross validation, diagnostic biomarker evaluation including magnetic resonance imaging based classifiers may be an important adjunct to clinical evaluation.

Reproducibility, reliability and variability of FA and MD in the older healthy population: A test-retest multiparametric analysis

•

In older healthy subjects, FA and MD show overall good test-retest reliability & reproducibility.

•

MD is sistematically more reproducible than FA across the entire brain anatomy.

•

FA is more reliable than MD in subcortical white matter regions.

•

In high reliability & low reproducibility regions, variability between subjects is high and statistical power is low.

•

In low reliability & high reproducibility regions, variability between subjects is low and statistical power is high.

Reduced betweenness centrality of a sensory-motor vestibular network in subclinical agoraphobia

Agoraphobic patients feel dizzy in crowded open spaces and respond to this symptom with excessive fear and avoidance. These clinical features show great similitude with the newly defined syndrome of persistent postural perceptual dizziness (PPPD). Patients with PPPD show decreased activity and connectivity in regions of the vestibular cortex. Due to the great overlap between these two conditions, we hypothesized that individuals with sub-clinical agoraphobia would show reduction in the connectivity features of these regions. We selected a group of healthy individuals from the Human Connectome Project that self-reported agoraphobia episodes, and compared it with a control group. We accurately matched the two groups for psychological measures and personality traits in order to study the neural correlates of vestibular symptoms independently of possible psychiatric vulnerabilities. We found that the agoraphobia group showed reduced betweenness centrality of a network encompassing key regions of the vestibular cortex. Dysfunctions of the vestibular cortex may explain the dizziness symptom for a disorder previously labelled as psychogenic.

MR Biomarkers of Degenerative Brain Disorders Derived From Diffusion Imaging

The incidence of neurodegenerative diseases has shown an increasing trend. These conditions typically cause progressive functional disability. Identification of robust biomarkers of neurodegenerative diseases is a key imperative to facilitate early identification of the pathological features and to foster a better understanding of the pathogenetic mechanisms of individual diseases. Diffusion tensor imaging (DTI) is the most widely used diffusion MRI technique for assessment of neurodegenerative diseases. The DTI parameters are promising biomarkers for evaluation of microstructural changes; however, some limitations of DTI restrict its wider clinical use. New diffusion MRI techniques, such as diffusion kurtosis imaging (DKI), bi-tensor DTI, and neurite orientation density and dispersion imaging (NODDI) have been demonstrated to provide value addition to DTI for evaluation of neurodegenerative diseases. In this review article, we summarize the key technical aspects and provide an overview of the current state of knowledge regarding the role of DKI, bi-tensor DTI, and NODDI as biomarkers of microstructural changes in representative neurodegenerative diseases including Alzheimer's disease, Parkinson's disease, amyotrophic lateral sclerosis, and Huntington's disease. LEVEL OF EVIDENCE: 5 TECHNICAL EFFICACY STAGE: 2 J. MAGN. RESON. IMAGING 2020;52:1620-1636.

Multishell diffusion imaging reveals sex-specific trajectories of early white matter degeneration in normal aging

During aging, human white matter (WM) is subject to dynamic structural changes which have a deep impact on healthy and pathological evolution of the brain through the lifespan; characterizing this pattern is of key importance for understanding brain development, maturation, and aging as well as for studying its pathological alterations. Diffusion magnetic resonance imaging (MRI) can provide a quantitative assessment of the white-matter microstructural organization that characterizes these trajectories. Here, we use both conventional and advanced diffusion MRI in a cohort of 91 individuals (age range: 13-62 years) to study region- and sex-specific features of WM microstructural integrity in healthy aging. We focus on the age at which microstructural imaging parameters invert their development trend as the time point which marks the onset of microstructural decline in WM. Importantly, our results indicate that age-related brain changes begin earlier in males than females and affect more frontal regions-in accordance with evolutionary theories and numerous evidences across non-MRI domains. Advanced diffusion MRI reveals age-related WM modification patterns which cannot be detected using conventional diffusion tensor imaging.

Correction for fast pseudo-diffusive fluid motion contaminations in diffusion tensor imaging

In this prospective study, we quantified the fast pseudo-diffusion contamination by blood perfusion or cerebrospinal fluid (CSF) intravoxel incoherent movements on the measurement of the diffusion tensor metrics in healthy brain tissue. Diffusion-weighted imaging (TR/TE = 4100 ms/90 ms; b-values: 0, 5, 10, 20, 35, 55, 80, 110, 150, 200, 300, 500, 750, 1000, 1300 s/mm², 20 diffusion-encoding directions) was performed on a cohort of five healthy volunteers at 3 Tesla. The projections of the diffusion tensor along each diffusion-encoding direction were computed using a two b-value approach (2b), by fitting the signal to a monoexponential curve (mono), and by correcting for fast pseudo-diffusion compartments using the biexponential intravoxel incoherent motion model (IVIM) (bi). Fractional anisotropy (FA) and mean diffusivity (MD) of the diffusion tensor were quantified in regions of interest drawn over white matter areas, gray matter areas, and the ventricles. A significant dependence of the MD from the evaluation method was found in all selected regions. A lower MD was computed when accounting for the fast-diffusion compartments. A larger dependence was found in the nucleus caudatus (bi: median 0.86 10^-3 mm²/s, Δ2b: -11.2%, Δmono: -14.4%; p = 0.007), in the anterior horn (bi: median 2.04 10^-3 mm²/s, Δ2b: -9.4%, Δmono: -11.5%, p = 0.007) and in the posterior horn of the lateral ventricles (bi: median 2.47 10^-3 mm²/s, Δ2b: -5.5%, Δmono: -11.7%; p = 0.007). Also for the FA, the signal modeling affected the computation of the anisotropy metrics. The deviation depended on the evaluated region with significant differences mainly in the nucleus caudatus (bi: median 0.15, Δ2b: +39.3%, Δmono: +14.7%; p = 0.022) and putamen (bi: median 0.19, Δ2b: +3.1%, Δmono: +17.3%; p = 0.015). Fast pseudo-diffusive regimes locally affect diffusion tensor imaging (DTI) metrics in the brain. Here, we propose the use of an IVIM-based method for correction of signal contaminations through CSF or perfusion.

Handedness Side and Magnetization Transfer Ratio in the Primary Sensorimotor Cortex Central Sulcus

Left-handers show lower asymmetry in manual ability when compared to right-handers. Unlike right-handers, left-handers do not show larger deactivation of the ipsilateral primary sensorimotor (SM1) cortex on functional magnetic resonance imaging when moving their dominant than their nondominant hand. However, it should be noted that morphometric MRI studies have reported no differences between right-handers and left-handers in volume, thickness, or surface area of the SM1 cortex. In this regard, magnetization transfer (MT) imaging is a technique with the potential to provide information on microstructural organization and macromolecular content of tissue. In particular, MT ratio index of the brain gray matter is assumed to reflect the variable content of afferent or efferent myelinated fibers, with higher MT ratio values being associated with increased fibers number or degree of myelination. The aim of this study was hence to assess, for the first time, through quantitative MT ratio measurements, potential differences in microstructural organization/characteristics of SM1 cortex between left- and right-handers, which could underlay handedness side. Nine left-handed and 9 right-handed healthy subjects, as determined by the Edinburgh handedness inventory, were examined with T1-weighted and MT-weighted imaging on a 3 T scanner. The hands of subjects were kept still during all acquisitions. Using FreeSurfer suite and the automatic anatomical labeling parcellations defined by the Destrieux atlas, we measured MT ratio, as well as cortical thickness, in three regions of interests corresponding to the precentral gyrus, the central sulcus, and the postcentral gyrus in the bilateral SM1 cortex. No significant difference between left- and right-handers was revealed in the thickness of the three partitions of the SM1 cortex. However, left-handers showed a significantly (p = 0.007) lower MT ratio (31.92%  ± 0.96%) in the right SM1 central sulcus (i.e., the hand representation area for left-handers) as compared to right-handers (33.28%  ± 0.94%). The results of this preliminary study indicate that quantitative MT imaging, unlike conventional morphometric MRI measurements, can be a useful tool to reveal, in SM1 cortex, potential microstructural correlates of handedness side.

Longer Repetition Time Proton MR Spectroscopy Shows Increasing Hippocampal and Parahippocampal Metabolite Concentrations with Aging

BACKGROUND AND PURPOSE

Previous magnetic resonance spectroscopy (MRS) studies have concluded that hippocampal and parahippocampal metabolite concentrations remain stable during healthy adult aging. However, these studies used short repetition times (TR ≤ 2 seconds), which lead to incomplete longitudinal magnetization recovery, and thus, heavily T₁ -weighted measurements. It is important to accurately characterize brain metabolites changes with age to enable appropriate interpretations of MRS findings in the context of neurodegenerative diseases. Our goal was to assess hippocampal brain metabolite concentrations in a large cohort of diversely aged healthy volunteers using a longer TR of 4 seconds.

METHODS

Left hippocampal MR spectra were collected from 38 healthy volunteers at 3T. Absolute metabolite concentrations were determined for total N-acetyl-aspartate (tNAA), total creatine (tCr), total choline (tCho), glutamate and glutamine (Glx), and myoinositol (mI). Individual partial correlations between each metabolite with age were assessed using demographic information and voxel compartmentation as confounders.

RESULTS

Hippocampal tNAA, tCr, tCho, and mI all increased with age (NAA: R² = .17, P = .041; tCr: R² = .45, P = .0002; tCho: R² = .37, P = .001; mI: R² = .44, P = .0003). There were no relationships between age and signal to noise ratio, linewidth, or scan date, indicating the correlations were not confounded by spectral quality. Furthermore, we did not observe a trend with age in the voxel tissue compartmentations.

CONCLUSIONS

We observed increases in hippocampal/parahippocampal metabolite concentrations with age, a finding that is in contrast to previous literature. Our findings illustrate the importance of using a sufficiently long TR in MRS to avoid T₁ -relaxation effects influencing the measurement of absolute metabolite concentrations.

Disruption of brain network organization in primary open angle glaucoma

Resting-state functional magnetic resonance imaging (rs-fMRI) is commonly employed to study changes in functional brain connectivity. Recently, the hypothesis of a brain involvement in primary open angle glaucoma has sprung interest for neuroimaging studies in this pathology. The purpose of this study is to evaluate a putative reorganization of brain networks in glaucomatous patients through graph-theoretical measures of integration, segregation and centrality by exploiting a multivariate networks association measure and a recently introduced global and local brain network disruption index. Nineteen glaucoma patients and sixteen healthy control subjects (age: 50 - 76, mean 61 years) underwent rs-fMRI examination at 3T. After preprocessing, rs-fMRI time series were parcellated into 116 regions (AAL atlas), adjacency matrices were computed based on partial correlations and graph-theoretical measures of integration, segregation and centrality as well as group-wise and subject-wise disruption index estimates were generated for all subjects. We found that the group-wise disruption index was negative and statistically different from 0 in for all graph theoretical metrics. Additionally, statistically significant group-wise differences in subject-wise disruption indexes were found in all local metrics. The differences in local network measures highlight cerebral reorganization of brain networks in glaucoma patients, supporting the interpretation of glaucoma as central nervous system disease, likely part of the heterogeneous group of recently described disconnection syndromes.

Recurrent neural networks for reconstructing complex directed brain connectivity

While Granger Causality(GC)-based approaches have been widely employed in a vast number of problems in network science, the vast majority of GC applications are based on linear multivariate autoregressive (MVAR) models. However, it is well known that real-life system (and biological networks in particular) exhibit notable nonlinear behavior, hence undermining that validity of MVAR-based approaches to estimating GC (MVAR-GC). In this paper, we define a novel approach to estimating nonlinear, directed within-network interactions based on a specific class of recurrent neural networks (RNN) termed echo-state networks (ESN). We reformulate the classical GC framework in terms of ESN-based models for multivariate signals generated by arbitrarily complex networks, and characterize the ability of our ESN-based Granger Causality (ES-GC) to capture nonlinear causal relations by simulating multivariate coupling in a network of nonlinearly interacting, noisy Duffing oscillators operating in a chaotic regime. Synthetic validation shows a net advantage of ES-GC over all other estimators in detecting nonlinear, causal links. We then explore the structure of EC-GC networks in the human brain in functional MRI data from 1003 healthy subjects scanned at rest at 3T, discovering previously unknown between-network interactions. In summary, ES-GC performs significantly better than commonly used and recently developed GC detection tools, making it a superior tool for the analysis of e.g. multivariate biological networks.

Image correction for diffusion tensor imaging of Rhesus monkey thoracic spinal cord

BACKGROUND

The relatively tiny spinal cord of non-human primate (NHP) causes increased challenge in diffusion tensor imaging (DTI) post-processing. This study aimed to establish a reliable correction strategy applied to clinical DTI images of NHP.

METHODS

Six normal and partial spinal cord injury (SCI) rhesus monkeys underwent 3T MR scanning. A correction strategy combining multiple iterations and non-rigid deformation was used for DTI image post-processing. Quantitative evaluations were then conducted to investigate effects of distortion correction.

RESULTS

After correction, longitudinal geometric distortion, global distortion, and residual distance errors were all significantly decreased (P < 0.05). Fractional anisotropy at the injured site was remarkably lower than that at the contralateral site (P = 0.0488) and was substantially lower than those at the adjacent superior (P = 0.0157) and inferior (P = 0.0128) areas at the same side.

CONCLUSIONS

Our image correction strategy can improve the quality of the DTI images of NHP thoracic cords, contributing to the development of SCI preclinical research.

White matter changes in the perforant path area in patients with amyotrophic lateral sclerosis

OBJECTIVE

The aim of this study was to test the hypothesis that white matter degeneration of the perforant path - as part of the Papez circuit - is a key feature of amyotrophic lateral sclerosis (ALS), even in the absence of frontotemporal dementia (FTD) or deposition of pTDP-43 inclusions in hippocampal granule cells.

METHODS

We used diffusion Magnetic Resonance Imaging (dMRI), polarized light imaging (PLI) and immunohistochemical analysis of post mortem hippocampus specimens from controls (n = 5) and ALS patients (n = 14) to study white matter degeneration in the perforant path.

RESULTS

diffusion Magnetic Resonance Imaging demonstrated a decrease in fractional anisotropy (P = 0.01) and an increase in mean diffusivity (P = 0.01) in the perforant path in ALS compared to controls. PLI-myelin density was lower in ALS (P = 0.05) and correlated with fractional anisotropy (r = 0.52, P = 0.03). These results were confirmed by immunohistochemistry; both myelin (proteolipid protein, P = 0.03) and neurofilaments (SMI-312, P = 0.02) were lower in ALS. Two out of the fourteen ALS cases showed pTDP-43 pathology in the dentate gyrus, but with comparable myelination levels in the perforant path to other ALS cases.

CONCLUSION

We conclude that degeneration of the perforant path occurs in ALS patients and that this may occur before, or independent of, pTDP-43 aggregation in the dentate gyrus of the hippocampus. Future research should focus on correlating the degree of cognitive decline to the amount of white matter atrophy in the perforant path.

Diffusion tensor distribution imaging

Conventional diffusion MRI yields voxel-averaged parameters that suffer from ambiguities for heterogeneous anisotropic materials such as brain tissue. Using principles from solid-state NMR spectroscopy, we have previously introduced the shape of the diffusion encoding tensor as a separate acquisition dimension that disentangles isotropic and anisotropic contributions to the observed diffusivities, thereby allowing for unconstrained data inversion into diffusion tensor distributions with "size," "shape," and orientation dimensions. Here we combine our recent non-parametric data inversion algorithm and data acquisition protocol with an imaging pulse sequence to demonstrate spatial mapping of diffusion tensor distributions using a previously developed composite phantom with multiple isotropic and anisotropic components. We propose a compact format for visualizing two-dimensional arrays of the distributions, new scalar parameters quantifying intra-voxel heterogeneity, and a binning procedure giving maps of all relevant parameters for each of the components resolved in the multidimensional distribution space.

Diffusional kurtosis imaging in head and neck cancer: On the use of trace-weighted images to estimate indices of non-Gaussian water diffusion

PURPOSE

While previous studies have demonstrated the feasibility and potential usefulness of quantitative non-Gaussian diffusional kurtosis imaging (DKI) of the brain, more recent research has focused on oncological application of DKI in various body regions such as prostate, breast, and head and neck (HN). Given the need to minimize scan time during most routine magnetic resonance imaging (MRI) acquisitions of body regions, diffusion-weighted imaging (DWI) with only three orthogonal diffusion weighting directions (x, y, z) is usually performed. Moreover, as water diffusion within malignant tumors is generically thought to be almost isotropic, DWI with only three diffusion weighting directions is considered sufficient for oncological application and it represents the de facto standard in body DKI. In this context, since the kurtosis tensor and diffusion tensor cannot be obtained, the averages of the three directional (K_x , K_y , K_z ) and (D_x , D_y , D_z ) - namely K and D, respectively - represent the best-possible surrogates of directionless DKI-derived indices of kurtosis and diffusivity, respectively. This would require fitting the DKI model to the diffusion-weighted images acquired along each direction (x, y, z) prior to averaging. However, there is a growing tendency to perform only a single fit of the DKI model to the geometric means of the images acquired with diffusion-sensitizing gradient along (x, y, z), referred to as trace-weighted (TW) images. To the best of our knowledge, no in vivo studies have evaluated how TW images affect estimates of DKI-derived indices of K and D. Thus, the aim of this study was to assess the potential bias and error introduced in estimated K and D by fitting the DKI model to the TW images in HN cancer patients.

METHODS

Eighteen patients with histologically proven malignant tumors of the HN were enrolled in the study. They underwent pretreatment 3 T MRI, including DWI (b-values: 0, 500, 1000, 1500, 2000 s/mm² ). Some patients had multiple lesions, and thus a total of 34 lesions were analyzed. DKI-derived indices were estimated, voxel-by-voxel, using single diffusion-weighted images along (x, y, z) as well as TW images. A comparison between the two estimation methods was performed by calculating the percentage error in D (D_err ) and K (K_err ). Also, diffusivity anisotropy (D_anis ) and diffusional kurtosis anisotropy (K_anis ) were estimated. Agreements between the two estimation methods were assessed by Bland-Altman plots. The Spearman rank correlation test was used to study the correlations between K_err /D_err and D_anis /K_anis. RESULTS: The median (95% confidence interval) K_err and D_err were 5.1% (0.8%, 32.6%) and 1.7% (-2.5%, 5.3%), respectively. A significant relationship was observed between K_err  and D_anis (correlation coefficient R = 0.694, P < 0.0001), as well as between K_err and K_anis (R = 0.848, P < 0.0001).

CONCLUSIONS

In HN cancer, the fit of the DKI model to TW images can introduce bias and error in the estimation of K and D, which may be non-negligible for single lesions, and should hence be adopted with caution.

Dependence of apparent diffusion coefficient measurement on diffusion gradient direction and spatial position - A quality assurance intercomparison study of forty-four scanners for quantitative diffusion-weighted imaging

PURPOSE

To propose an MRI quality assurance procedure that can be used for routine controls and multi-centre comparison of different MR-scanners for quantitative diffusion-weighted imaging (DWI).

MATERIALS AND METHODS

44 MR-scanners with different field strengths (1 T, 1.5 T and 3 T) were included in the study. DWI acquisitions (b-value range 0-1000 s/mm²), with three different orthogonal diffusion gradient directions, were performed for each MR-scanner. All DWI acquisitions were performed by using a standard spherical plastic doped water phantom. Phantom solution ADC value and its dependence with temperature was measured using a DOSY sequence on a 600 MHz NMR spectrometer. Apparent diffusion coefficient (ADC) along each diffusion gradient direction and mean ADC were estimated, both at magnet isocentre and in six different position 50 mm away from isocentre, along positive and negative AP, RL and HF directions.

RESULTS

A good agreement was found between the nominal and measured mean ADC at isocentre: more than 90% of mean ADC measurements were within 5% from the nominal value, and the highest deviation was 11.3%. Away from isocentre, the effect of the diffusion gradient direction on ADC estimation was larger than 5% in 47% of included scanners and a spatial non uniformity larger than 5% was reported in 13% of centres.

CONCLUSION

ADC accuracy and spatial uniformity can vary appreciably depending on MR scanner model, sequence implementation (i.e. gradient diffusion direction) and hardware characteristics. The DWI quality assurance protocol proposed in this study can be employed in order to assess the accuracy and spatial uniformity of estimated ADC values, in single- as well as multi-centre studies.

Alterations of the nigrostriatal pathway in a 6-OHDA rat model of Parkinson's disease evaluated with multimodal MRI

Parkinson's disease is characterized by neurodegeneration of the dopaminergic neurons in the substantia nigra pars compacta. The 6-hydroxydopamine (6-OHDA) rat model has been used to study neurodegeneration in the nigro-striatal dopaminergic system. The goal of this study was to evaluate the reliability of diffusion MRI and resting-state functional MRI biomarkers in monitoring neurodegeneration in the 6-OHDA rat model assessed by quantitative histology. We performed a unilateral injection of 6-OHDA in the striatum of Sprague Dawley rats to produce retrograde degeneration of the dopamine neurons in the substantia nigra pars compacta. We carried out a longitudinal study with a multi-modal approach combining structural and functional MRI together with quantitative histological validation to follow the effects of the lesion. Functional and structural connectivity were assessed in the brain of 6-OHDA rats and sham rats (NaCl injection) at 3 and 6 weeks post-lesioning using resting-state functional MRI and diffusion-weighted. Our results showed (i) increased functional connectivity in ipsi- and contra-lesioned regions of the cortico-basal ganglia network pathway including the motor cortex, the globus pallidus, and the striatum regions at 3 weeks; (ii) increased fractional anisotropy (FA) in the ipsi- and contralateral striatum of the 6-OHDA group at 3 weeks, and increased axial diffusivity (AD) and mean diffusivity in the ipsilateral striatum at 6 weeks; (iii) a trend for increased FA in both substantia nigra of the 6-OHDA group at 3 weeks. Optical density measurements of tyrosine-hydroxylase (TH) staining of the striatum showed good correlations with the FA and AD measurements in the striatum. No correlations were found between the number of TH-stained dopaminergic neurons and MRI measurements in the substantia nigra. This study suggested that (i) FA and AD were reliable biomarkers to evaluate neurodegeneration in the cortico-basal ganglia network of the 6-OHDA model, (ii) diffusion MRI and resting-state functional MRI (rsfMRI) were not sensitive enough to detect changes in the substantia nigra in this model.

Early axonal damage in normal appearing white matter in multiple sclerosis: Novel insights from multi-shell diffusion MRI

Conventional diffusion-weighted MR imaging techniques provide limited specificity in disentangling disease-related microstructural alterations involving changes in both axonal density and myelination. By simultaneously probing multiple diffusion regimens, multi-shell diffusion MRI is capable of increasing specificity to different tissue sub-compartments and hence separate different contributions to changes in diffusion-weighted signal attenuation. Advanced multi-shell diffusion models impose significant requirements on the amount of diffusion weighting (i.e. gradient coil performance) and angular resolution (i.e. in-scanner subject time), which commonly limits their applicability in a clinical setting. In this paper, we apply a high-b-value, high angular resolution multi-shell diffusion MRI protocol to a population of early multiple sclerosis (MS) patients and healthy controls. Through the Composite Hindered and Restricted Model of Diffusion (CHARMED) model, we extract indices for axonal density as well as parameters sensitive to myelin. We demonstrate increased sensitivity to microstructural changes in normal appearing white matter and in lesions in MS as compared to traditional models like DTI. These changes appear to be predominantly in axonal density, pointing towards the existence of axonal damage mechanisms in early MS.

Histogram analysis of DTI-derived indices reveals pontocerebellar degeneration and its progression in SCA2

PURPOSE

To assess the potential of histogram metrics of diffusion-tensor imaging (DTI)-derived indices in revealing neurodegeneration and its progression in spinocerebellar ataxia type 2 (SCA2).

MATERIALS AND METHODS

Nine SCA2 patients and 16 age-matched healthy controls, were examined twice (SCA2 patients 3.6±0.7 years and controls 3.3±1.0 years apart) on the same 1.5T scanner by acquiring T1-weighted and diffusion-weighted (b-value = 1000 s/mm2) images. Cerebrum and brainstem-cerebellum regions were segmented using FreeSurfer suite. Histogram analysis of DTI-derived indices, including mean diffusivity (MD), fractional anisotropy (FA), axial (AD) / radial (RD) diffusivity and mode of anisotropy (MO), was performed.

RESULTS

At baseline, significant differences between SCA2 patients and controls were confined to brainstem-cerebellum. Median values of MD/AD/RD and FA/MO were significantly (p<0.001) higher and lower, respectively, in SCA2 patients (1.11/1.30/1.03×10(-3) mm2/s and 0.14/0.19) than in controls (0.80/1.00/0.70×10(-3) mm2/s and 0.20/0.41). Also, peak location values of MD/AD/RD and FA were significantly (p<0.001) higher and lower, respectively, in SCA2 patients (0.91/1.11/0.81×10(-3) mm2/s and 0.12) than in controls (0.71/0.91/0.63×10(-3) mm2/s and 0.18). Peak height values of FA and MD/AD/RD/MO were significantly (p<0.001) higher and lower, respectively, in SCA2 patients (0.20 and 0.07/0.06/0.07×10(-3) mm2/s/year /0.07) than in controls (0.15 and 0.14/0.11/0.12/×10(-3) mm2/s/year /0.09). The rate of change of MD median values was significantly (p<0.001) higher (i.e., increased) in SCA2 patients (0.010×10(-3) mm2/s/year) than in controls (-0.003×10(-3) mm2/s/year) in the brainstem-cerebellum, whereas no significant difference was found for other indices and in the cerebrum.

CONCLUSION

Histogram analysis of DTI-derived indices is a relatively straightforward approach which reveals microstructural changes associated with pontocerebellar degeneration in SCA2 and the median value of MD is capable to track its progression.

Diffusional kurtosis imaging and white matter microstructure modeling in a clinical study of major depressive disorder

Major depressive disorder (MDD) is a globally prevalent psychiatric disorder that results from disruption of multiple neural circuits involved in emotional regulation. Although previous studies using diffusion tensor imaging (DTI) found smaller values of fractional anisotropy (FA) in the white matter, predominantly in the frontal lobe, of patients with MDD, studies using diffusion kurtosis imaging (DKI) are scarce. Here, we used DKI whole-brain analysis with tract-based spatial statistics (TBSS) to investigate the brain microstructural abnormalities in MDD. Twenty-six patients with MDD and 42 age- and sex-matched control subjects were enrolled. To investigate the microstructural pathology underlying the observations in DKI, a compartment model analysis was conducted focusing on the corpus callosum. In TBSS, the patients with MDD showed significantly smaller values of FA in the genu and frontal portion of the body of the corpus callosum. The patients also had smaller values of mean kurtosis (MK) and radial kurtosis (RK), but MK and RK abnormalities were distributed more widely compared with FA, predominantly in the frontal lobe but also in the parietal, occipital, and temporal lobes. Within the callosum, the regions with smaller MK and RK were located more posteriorly than the region with smaller FA. Model analysis suggested significantly smaller values of intra-neurite signal fraction in the body of the callosum and greater fiber dispersion in the genu, which were compatible with the existing literature of white matter pathology in MDD. Our results show that DKI is capable of demonstrating microstructural alterations in the brains of patients with MDD that cannot be fully depicted by conventional DTI. Though the issues of model validation and parameter estimation still remain, it is suggested that diffusion MRI combined with a biophysical model is a promising approach for investigation of the pathophysiology of MDD.

Longitudinal white matter microstructural change in Parkinson's disease

Postmortem studies of Parkinson's disease (PD) suggest that Lewy body pathology accumulates in a predictable topographical sequence, beginning in the olfactory bulb, followed by caudal brainstem, substantia nigra, limbic cortex, and neocortex. Diffusion-weighted imaging (DWI) is sensitive, if not specific, to early disease-related white matter (WM) change in a variety of traumatic and degenerative brain diseases. Although numerous cross-sectional studies have reported DWI differences in cerebral WM in PD, only a few longitudinal studies have investigated whether DWI change exceeds that of normal aging or coincides with regional Lewy body accumulation. This study mapped regional differences in the rate of DWI-based microstructural change between 29 PD patients and 43 age-matched controls over 18 months. Iterative within- and between-subject tensor-based registration was completed on motion- and eddy current-corrected DWI images, then baseline versus follow-up difference maps of fractional anisotropy, mean, radial, and axial diffusivity were analyzed in the Biological Parametric Mapping toolbox for MATLAB. This analysis showed that PD patients had a greater decline in WM integrity in the rostral brainstem, caudal subcortical WM, and cerebellar peduncles, compared with controls. In addition, patients with unilateral clinical signs at baseline experienced a greater rate of WM change over the 18-month study than patients with bilateral signs. These findings suggest that rate of WM microstructural change in PD exceeds that of normal aging and is maximal during early stage disease. In addition, the neuroanatomic locations (rostral brainstem and subcortical WM) of accelerated WM change fit with current theories of topographic disease progression.

Fractional cable equation for general geometry: A model of axons with swellings and anomalous diffusion

Different experimental studies have reported anomalous diffusion in brain tissues and notably this anomalous diffusion is expressed through fractional derivatives. Axons are important to understand neurodegenerative diseases such as multiple sclerosis, Alzheimer's disease, and Parkinson's disease. Indeed, abnormal accumulation of proteins and organelles in axons is a hallmark of these diseases. The diffusion in the axons can become anomalous as a result of this abnormality. In this case the voltage propagation in axons is affected. Another hallmark of different neurodegenerative diseases is given by discrete swellings along the axon. In order to model the voltage propagation in axons with anomalous diffusion and swellings, in this paper we propose a fractional cable equation for a general geometry. This generalized equation depends on fractional parameters and geometric quantities such as the curvature and torsion of the cable. For a cable with a constant radius we show that the voltage decreases when the fractional effect increases. In cables with swellings we find that when the fractional effect or the swelling radius increases, the voltage decreases. Similar behavior is obtained when the number of swellings and the fractional effect increase. Moreover, we find that when the radius swelling (or the number of swellings) and the fractional effect increase at the same time, the voltage dramatically decreases.

Retained executive abilities in mild cognitive impairment are associated with increased white matter network connectivity

PURPOSE

To describe structural network differences in individuals with mild cognitive impairment (MCI) with high versus low executive abilities, as reflected by measures of white matter connectivity using diffusion tensor imaging (DTI).

MATERIALS AND METHODS

This was a retrospective, cross-sectional study. Of the 128 participants from the Alzheimer's Disease Neuroimaging Initiative database who had both a DTI scan as well as a diagnosis of MCI, we used an executive function score to classify the top 15 scoring patients as high executive ability, and the bottom-scoring 16 patients as low executive ability. Using a regions-of-interest-based analysis, we constructed networks and calculated graph theory measures on the constructed networks. We used automated tractography in order to compare differences in major white matter tracts.

RESULTS

The high executive ability group yielded greater network size, density and clustering coefficient. The high executive ability group reflected greater fractional anisotropy bilaterally in the inferior and superior longitudinal fasciculi.

CONCLUSIONS

The network measures of the high executive ability group demonstrated greater white matter integrity. This suggests that white matter reserve may confer greater protection of executive abilities. Loss of this reserve may lead to greater impairment in the progression to Alzheimer's disease dementia.

KEY POINTS

• The MCI high executive ability group yielded a larger network. • The MCI high executive ability group had greater FA in numerous tracts. • White matter reserve may confer greater protection of executive abilities. • Loss of executive reserve may lead to greater impairment in AD dementia.

Sensitivity of diffusion MRI to perilesional reactive astrogliosis in focal ischemia

Reactive astrogliosis is a response to injury in the central nervous system that plays an essential role in inflammation and tissue repair. It is characterized by hypertrophy of astrocytes, alterations in astrocyte gene expression and astrocyte proliferation. Reactive astrogliosis occurs in multiple neuropathologies, including stroke, traumatic brain injury and Alzheimer's disease, and it has been proposed as a possible source of the changes in diffusion magnetic resonance imaging (dMRI) metrics observed with these diseases. In this study, the sensitivity of dMRI to reactive astrogliosis was tested in an animal model of focal acute and subacute ischemia induced by the vasoconstricting peptide, endothelin-1. Reactive astrogliosis in perilesional cortex was quantified by calculating the astrocyte surface density as determined with a glial fibrillary acidic protein (GFAP) antibody, whereas perilesional diffusion changes were measured in vivo with diffusional kurtosis imaging. We found substantial changes in the surface density of GFAP-positive astrocyte processes and modest changes in dMRI metrics in the perilesional motor cortex following stroke. Although there are time point-specific correlations between dMRI and histological measures, there is no definitive evidence for a causal relationship.

Distribution-aware estimation of the minimum achievable uncertainty in diffusion-tensor imaging (DTI)

Diffusion tensor imaging (DTI) provides exquisite sensitivity to structural and microstructural characteristics of brain tissue, and is routinely employed in advanced neuroimaging applications. DTI is commonly performed using intrinsically noisy echo-planar imaging techniques and poses high demands both on scanner performance and on in-scanner subject time, which in turn is directly related to the number of diffusion-weighting direction one requires. While DTI-derived indices such as fractional anisotropy (FA), diffusion tensor trace and anisotropy mode have proven extremely useful in characterizing disease-related aberrations, their estimation is commonly performed using fitting routines that do not properly take into account MRI noise distribution. In this paper, we present a distribution-aware maximum likelihood tensor estimation framework which also allows, for the first time, separate local noise estimation in both diffusion weighted and reference images. We validate our framework using multiple water phantom diffusion weighted acquisitions, and demonstrate its feasibility in human data. We then employ our framework within Monte Carlo simulations to show how the minimum achievable uncertainty attainable in DTI depends on signal-to-noise ratio (SNR) and number of diffusion gradient directions, demonstrating that these dependencies could be recast into simple power laws which may serve as guidelines for application-specific DTI protocol design.

Analytical performance bounds for multi-tensor diffusion-MRI

PURPOSE

To examine the effects of MR acquisition parameters on brain white matter fiber orientation estimation and parameter of clinical interest in crossing fiber areas based on the Multi-Tensor Model (MTM).

MATERIAL AND METHODS

We compute the Cramér-Rao Bound (CRB) for the MTM and the parameter of clinical interest such as the Fractional Anisotropy (FA) and the dominant fiber orientations, assuming that the diffusion MRI data are recorded by a multi-coil, multi-shell acquisition system. Considering the sum-of-squares method for the reconstructed magnitude image, we introduce an approximate closed-form formula for Fisher Information Matrix that has the simplicity and easy interpretation advantages. In addition, we propose to generalize the FA and the mean diffusivity to the multi-tensor model.

RESULTS

We show the application of the CRB to reduce the scan time while preserving a good estimation precision. We provide results showing how the increase of the number of acquisition coils compensates the decrease of the number of diffusion gradient directions. We analyze the impact of the b-value and the Signal-to-Noise Ratio (SNR). The analysis shows that the estimation error variance decreases with a quadratic rate with the SNR, and that the optimum b-values are not unique but depend on the target parameter, the context, and eventually the target cost function.

CONCLUSION

In this study we highlight the importance of choosing the appropriate acquisition parameters especially when dealing with crossing fiber areas. We also provide a methodology for the optimal tuning of these parameters using the CRB.

Diffusion alterations associated with Parkinson's disease symptomatology: A review of the literature

Parkinson's disease (PD) is a heterogeneous neurological disorder with a variety of motor and non-motor symptoms. The underlying mechanisms of these symptoms are not fully understood. An increased interest in structural connectivity analyses using diffusion tensor imaging (DTI) in PD has led to an expansion of our understanding of the impact of abnormalities in diffusivity on phenotype. This review outlines the contribution of these abnormalities to symptoms of PD including bradykinesia, tremor and non-tremor phenotypes, freezing of gait, cognitive impairment, mood, sleep disturbances, visual hallucinations and olfactory dysfunction. Studies have shown that impairments in cognitive functioning are related to diffusion abnormalities in frontal and parietal regions, as well as in the corpus callosum and major fibres connecting midbrain and subcortical structures with the neocortex. However, the impact of diffusion alterations on motor, mood and other symptoms of PD are less well understood. The findings presented here highlight the challenges faced and the potential areas of future research avenues where DTI may be beneficial. Larger cohort studies and standardized imaging protocols are required to investigate current promising preliminary findings.