Age-related white matter changes in high b-value q-space diffusion-weighted imaging

A practical evaluation of measures derived from compressed sensing diffusion spectrum imaging

Diffusion Spectrum Imaging (DSI) using dense Cartesian sampling of q-space has been shown to provide important advantages for modeling complex white matter architecture. However, its adoption has been limited by the lengthy acquisition time required. Sparser sampling of q-space combined with compressed sensing (CS) reconstruction techniques has been proposed as a way to reduce the scan time of DSI acquisitions. However prior studies have mainly evaluated CS-DSI in post-mortem or non-human data. At present, the capacity for CS-DSI to provide accurate and reliable measures of white matter anatomy and microstructure in the living human brain remains unclear. We evaluated the accuracy and inter-scan reliability of 6 different CS-DSI schemes that provided up to 80% reductions in scan time compared to a full DSI scheme. We capitalized on a dataset of 26 participants who were scanned over eight independent sessions using a full DSI scheme. From this full DSI scheme, we subsampled images to create a range of CS-DSI images. This allowed us to compare the accuracy and inter-scan reliability of derived measures of white matter structure (bundle segmentation, voxel-wise scalar maps) produced by the CS-DSI and the full DSI schemes. We found that CS-DSI estimates of both bundle segmentations and voxel-wise scalars were nearly as accurate and reliable as those generated by the full DSI scheme. Moreover, we found that the accuracy and reliability of CS-DSI was higher in white matter bundles that were more reliably segmented by the full DSI scheme. As a final step, we replicated the accuracy of CS-DSI in a prospectively acquired dataset (n = 20, scanned once). Together, these results illustrate the utility of CS-DSI for reliably delineating in vivo white matter architecture in a fraction of the scan time, underscoring its promise for both clinical and research applications.

Predicting peritumoral edema development after gamma knife radiosurgery of meningiomas using machine learning methods: a multicenter study

OBJECTIVES

Edema is a complication of gamma knife radiosurgery (GKS) in meningioma patients that leads to a variety of consequences. The aim of this study is to construct radiomics-based machine learning models to predict post-GKS edema development.

METHODS

In total, 445 meningioma patients who underwent GKS in our institution were enrolled and partitioned into training and internal validation datasets (8:2). A total of 150 cases from multicenter data were included as the external validation dataset. In each case, 1132 radiomics features were extracted from each pre-treatment MRI sequence (contrast-enhanced T1WI, T2WI, and ADC maps). Nine clinical features and eight semantic features were also generated. Nineteen random survival forest (RSF) and nineteen neural network (DeepSurv) models with different combinations of radiomics, clinical, and semantic features were developed with the training dataset, and evaluated with internal and external validation. A nomogram was derived from the model achieving the highest C-index in external validation.

RESULTS

All the models were successfully validated on both validation datasets. The RSF model incorporating clinical, semantic, and ADC radiomics features achieved the best performance with a C-index of 0.861 (95% CI: 0.748-0.975) in internal validation, and 0.780 (95% CI: 0.673-0.887) in external validation. It stratifies high-risk and low-risk cases effectively. The nomogram based on the predicted risks provided personalized prediction with a C-index of 0.962 (95%CI: 0.951-0.973) and satisfactory calibration.

CONCLUSION

This RSF model with a nomogram could represent a non-invasive and cost-effective tool to predict post-GKS edema risk, thus facilitating personalized decision-making in meningioma treatment.

CLINICAL RELEVANCE STATEMENT

The RSF model with a nomogram built in this study represents a handy, non-invasive, and cost-effective tool for meningioma patients to assist in better counselling on the risks, appropriate individual treatment decisions, and customized follow-up plans.

KEY POINTS

• Machine learning models were built to predict post-GKS edema in meningioma. The random survival forest model with clinical, semantic, and ADC radiomics features achieved excellent performance. • The nomogram based on the predicted risks provides personalized prediction with a C-index of 0.962 (95%CI: 0.951-0.973) and satisfactory calibration and shows the potential to assist in better counselling, appropriate treatment decisions, and customized follow-up plans. • Given the excellent performance and convenient acquisition of the conventional sequence, we envision that this non-invasive and cost-effective tool will facilitate personalized medicine in meningioma treatment.

High-fidelity approximation of grid- and shell-based sampling schemes from undersampled DSI using compressed sensing: Post mortem validation

While many useful microstructural indices, as well as orientation distribution functions, can be obtained from multi-shell dMRI data, there is growing interest in exploring the richer set of microstructural features that can be extracted from the full ensemble average propagator (EAP). The EAP can be readily computed from diffusion spectrum imaging (DSI) data, at the cost of a very lengthy acquisition. Compressed sensing (CS) has been used to make DSI more practical by reducing its acquisition time. CS applied to DSI (CS-DSI) attempts to reconstruct the EAP from significantly undersampled q-space data. We present a post mortem validation study where we evaluate the ability of CS-DSI to approximate not only fully sampled DSI but also multi-shell acquisitions with high fidelity. Human brain samples are imaged with high-resolution DSI at 9.4T and with polarization-sensitive optical coherence tomography (PSOCT). The latter provides direct measurements of axonal orientations at microscopic resolutions, allowing us to evaluate the mesoscopic orientation estimates obtained from diffusion MRI, in terms of their angular error and the presence of spurious peaks. We test two fast, dictionary-based, L2-regularized algorithms for CS-DSI reconstruction. We find that, for a CS acceleration factor of R=3, i.e., an acquisition with 171 gradient directions, one of these methods is able to achieve both low angular error and low number of spurious peaks. With a scan length similar to that of high angular resolution multi-shell acquisition schemes, this CS-DSI approach is able to approximate both fully sampled DSI and multi-shell data with high accuracy. Thus it is suitable for orientation reconstruction and microstructural modeling techniques that require either grid- or shell-based acquisitions. We find that the signal-to-noise ratio (SNR) of the training data used to construct the dictionary can have an impact on the accuracy of CS-DSI, but that there is substantial robustness to loss of SNR in the test data. Finally, we show that, as the CS acceleration factor increases beyond R=3, the accuracy of these reconstruction methods degrade, either in terms of the angular error, or in terms of the number of spurious peaks. Our results provide useful benchmarks for the future development of even more efficient q-space acceleration techniques.

Generalized diffusion spectrum magnetic resonance imaging (GDSI) for model-free reconstruction of the ensemble average propagator

Diffusion spectrum MRI (DSI) provides model-free estimation of the diffusion ensemble average propagator (EAP) and orientation distribution function (ODF) but requires the diffusion data to be acquired on a Cartesian q-space grid. Multi-shell diffusion acquisitions are more flexible and more commonly acquired but have, thus far, only been compatible with model-based analysis methods. Here, we propose a generalized DSI (GDSI) framework to recover the EAP from multi-shell diffusion MRI data. The proposed GDSI approach corrects for q-space sampling density non-uniformity using a fast geometrical approach. The EAP is directly calculated in a preferable coordinate system by multiplying the sampling density corrected q-space signals by a discrete Fourier transform matrix, without any need for gridding. The EAP is demonstrated as a way to map diffusion patterns in brain regions such as the thalamus, cortex and brainstem where the tissue microstructure is not as well characterized as in white matter. Scalar metrics such as the zero displacement probability and displacement distances at different fractions of the zero displacement probability were computed from the recovered EAP to characterize the diffusion pattern within each voxel. The probability averaged across directions at a specific displacement distance provides a diffusion property based image contrast that clearly differentiates tissue types. The displacement distance at the first zero crossing of the EAP averaged across directions orthogonal to the primary fiber orientation in the corpus callosum is found to be larger in the body (5.65 ± 0.09 μm) than in the genu (5.55 ± 0.15 μm) and splenium (5.4 ± 0.15 μm) of the corpus callosum, which corresponds well to prior histological studies. The EAP also provides model-free representations of angular structure such as the diffusion ODF, which allows estimation and comparison of fiber orientations from both the model-free and model-based methods on the same multi-shell data. For the model-free methods, detection of crossing fibers is found to be strongly dependent on the maximum b-value and less sensitive compared to the model-based methods. In conclusion, our study provides a generalized DSI approach that allows flexible reconstruction of the diffusion EAP and ODF from multi-shell diffusion data and data acquired with other sampling patterns.

Diffusion MRI of the spinal cord: from structural studies to pathology

Diffusion MRI is extensively used to study brain microarchitecture and pathologies, and water diffusion appears highly anisotropic in the white matter (WM) of the spinal cord (SC). Despite these facts, the use of diffusion MRI to study the SC, which has increased in recent years, is much less common than that in the brain. In the present review, after a brief outline of early studies of diffusion MRI (DWI) and diffusion tensor MRI (DTI) of the SC, we provide a short survey on DTI and on diffusion MRI methods beyond the tensor that have been used to study SC microstructure and pathologies. After introducing the porous view of WM and describing the q-space approach and q-space diffusion MRI (QSI), we describe other methodologies that can be applied to study the SC. Selected applications of the use of DTI, QSI, and other more advanced diffusion MRI methods to study SC microstructure and pathologies are presented, with some emphasis on the use of less conventional diffusion methodologies. Because of length constraints, we concentrate on structural studies and on a few selected pathologies. Examples of the use of diffusion MRI to study dysmyelination, demyelination as in experimental autoimmune encephalomyelitis and multiple sclerosis, amyotrophic lateral sclerosis, and traumatic SC injury are presented. We conclude with a brief summary and a discussion of challenges and future directions for diffusion MRI of the SC. Copyright © 2016 John Wiley & Sons, Ltd.

Associations among q-space MRI, diffusion-weighted MRI and histopathological parameters in meningiomas

OBJECTIVES

The purposes of this MR-based study were to calculate q-space imaging (QSI)-derived mean displacement (MDP) in meningiomas, to evaluate the correlation of MDP values with apparent diffusion coefficient (ADC) and to investigate the relationships among these diffusion parameters, tumour cell count (TCC) and MIB-1 labelling index (LI).

METHODS

MRI, including QSI and conventional diffusion-weighted imaging (DWI), was performed in 44 meningioma patients (52 lesions). ADC and MDP maps were acquired from post-processing of the data. Quantitative analyses of these maps were performed by applying regions of interest. Pearson correlation coefficients were calculated for ADC and MDP in all lesions and for ADC and TCC, MDP and TCC, ADC and MIB-1 LI, and MDP and MIB-1 LI in 17 patients who underwent subsequent surgery.

RESULTS

ADC and MDP values were found to have a strong correlation: r = 0.78 (P = <0.0001). Both ADC and MDP values had a significant negative association with TCC: r = -0.53 (p = 0.02) and -0.48 (P = 0.04), respectively. MIB-1 LI was not, however, found to have a significant association with these diffusion parameters.

CONCLUSION

In meningiomas, both ADC and MDP may be representative of cell density.

KEY POINTS

• Diffusion-weighted MRI offers possibilities to assess the aggressiveness of meningiomas. • The q-space imaging-derived mean displacement correlates strongly with apparent diffusion coefficients. • Both diffusion parameters showed a strong negative association with tumour cell counts. • Derived mean displacement may help assess the aggressiveness of meningiomas preoperatively.