Optimal experimental design for diffusion kurtosis imaging

Conditional generative diffusion deep learning for accelerated diffusion tensor and kurtosis imaging

PURPOSE

The purpose of this study was to develop DiffDL, a generative diffusion probabilistic model designed to produce high-quality diffusion tensor imaging (DTI) and diffusion kurtosis imaging (DKI) metrics from a reduced set of diffusion-weighted images (DWIs). This model addresses the challenge of prolonged data acquisition times in diffusion MRI while preserving metric accuracy.

METHODS

DiffDL was trained using data from the Human Connectome Project, including 300 training/validation subjects and 50 testing subjects. High-quality DTI and DKI metrics were generated using many DWIs and combined with subsets of DWIs to form training pairs. A UNet architecture was used for denoising, trained over 500 epochs with a linear noise schedule. Performance was evaluated against conventional DTI/DKI modeling and a reference UNet model using normalized mean absolute error (NMAE), peak signal-to-noise ratio (PSNR), and Pearson correlation coefficient (PCC).

RESULTS

DiffDL showed significant improvements in the quality and accuracy of fractional anisotropy (FA) and mean diffusivity (MD) maps compared to conventional methods and the baseline UNet model. For DKI metrics, DiffDL outperformed conventional DKI modeling and the UNet model across various acceleration scenarios. Quantitative analysis demonstrated superior NMAE, PSNR, and PCC values for DiffDL, capturing the full dynamic range of DTI and DKI metrics. The generative nature of DiffDL allowed for multiple predictions, enabling uncertainty quantification and enhancing performance.

CONCLUSION

The DiffDL framework demonstrated the potential to significantly reduce data acquisition times in diffusion MRI while maintaining high metric quality. Future research should focus on optimizing computational demands and validating the model with clinical cohorts and standard MRI scanners.

Surface-based Analyses of Diffusional Kurtosis Imaging in Amyotrophic Lateral Sclerosis: Relationship with Onset Subtypes

PURPOSE

Here, we aimed to characterize the cortical and subcortical microstructural alterations in the brains of patients with amyotrophic lateral sclerosis (ALS). In particular, we compared these features between bulbar-onset ALS (b-ALS) and limb-onset ALS (l-ALS).

METHODS

Diffusion MRI data (b = 0, 700, 2000 ms/mm2, 1.7-mm isotropic voxel) from 28 patients with ALS (9 b-ALS and 19 l-ALS) and 17 healthy control subjects (HCs) were analyzed. Diffusional kurtosis imaging (DKI) metrics were sampled at the mid-cortical and subcortical surfaces. We used permutation testing with a nonparametric combination of mean diffusivity (MD), fractional anisotropy (FA), and mean kurtosis (MK) to assess intergroup differences over the cerebrum. We also carried out an atlas-based analysis focusing on Brodmann Area 4 and 6 (primary motor and premotor areas) and investigated the correlation between MRI metrics and clinical parameters.

RESULTS

At both the mid-cortical and subcortical surfaces, b-ALS was associated with significantly greater MD, smaller FA, and smaller MK in the motor and premotor areas than HC. In contrast, the patients with l-ALS showed relatively moderate differences relative to HCs. The ALS Functional Rating Scale-Revised bulbar subscore was significantly correlated with the diffusion metrics in Brodmann Area 4.

CONCLUSION

The distribution of abnormalities over the cerebral hemispheres and the more severe microstructural alteration in b-ALS compared to l-ALS were in good agreement with findings from postmortem histology. Our results suggest the feasibility of surface-based DKI analyses for exploring brain microstructural pathologies in ALS. The observed differences between b-ALS and l-ALS and their correlations with functional bulbar impairment support the clinical relevance of DKI measurement in the cortical and juxtacortical regions of patients with ALS.

White matter and latency of visual evoked potentials during maturation: A miniature pig model of adolescent development

BACKGROUND

Continuous myelination of cerebral white matter (WM) during adolescence overlaps with the formation of higher cognitive skills and the onset of many neuropsychiatric disorders. We developed a miniature-pig model of adolescent brain development for neuroimaging and neurophysiological assessment during this critical period. Minipigs have gyroencephalic brains with a large cerebral WM compartment and a well-defined adolescence period.

METHODS

Eight Sinclair™ minipigs (Sus scrofa domestica) were evaluated four times during weeks 14-28 (40, 28 and 28 days apart) of adolescence using monocular visual stimulation (1 Hz)-evoked potentials and diffusion MRI (dMRI) of WM. The latency for the pre-positive 30 ms (PP30), positive 30 ms (P30) and negative 50 ms (N50) components of the flash visual evoked potentials (fVEPs) and their interhemispheric latency (IL) were recorded in the frontal, central and occipital areas during ten 60-second stimulations for each eye. The dMRI imaging protocol consisted of fifteen b-shells (b = 0-3500 s/mm2) with 32 directions/shell, providing measurements that included fractional anisotropy (FA), radial kurtosis, kurtosis anisotropy (KA), axonal water fraction (AWF), and the permeability-diffusivity index (PDI).

RESULTS

Significant reductions (p < 0.05) in the latency and IL of fVEP measurements paralleled significant rises in FA, KA, AWF and PDI over the same period. The longitudinal latency changes in fVEPs were primarily associated with whole-brain changes in diffusion parameters, while fVEP IL changes were related to maturation of the corpus callosum.

CONCLUSIONS

Good agreement between reduction in the latency of fVEPs and maturation of cerebral WM was interpreted as evidence for ongoing myelination and confirmation of the minipig as a viable research platform. Adolescent development in minipigs can be studied using human neuroimaging and neurophysiological protocols and followed up with more invasive assays to investigate key neurodevelopmental hypotheses in psychiatry.

Robust, fast and accurate mapping of diffusional mean kurtosis

Diffusional kurtosis imaging (DKI) is a methodology for measuring the extent of non-Gaussian diffusion in biological tissue, which has shown great promise in clinical diagnosis, treatment planning, and monitoring of many neurological diseases and disorders. However, robust, fast, and accurate estimation of kurtosis from clinically feasible data acquisitions remains a challenge. In this study, we first outline a new accurate approach of estimating mean kurtosis via the sub-diffusion mathematical framework. Crucially, this extension of the conventional DKI overcomes the limitation on the maximum b-value of the latter. Kurtosis and diffusivity can now be simply computed as functions of the sub-diffusion model parameters. Second, we propose a new fast and robust fitting procedure to estimate the sub-diffusion model parameters using two diffusion times without increasing acquisition time as for the conventional DKI. Third, our sub-diffusion-based kurtosis mapping method is evaluated using both simulations and the Connectome 1.0 human brain data. Exquisite tissue contrast is achieved even when the diffusion encoded data is collected in only minutes. In summary, our findings suggest robust, fast, and accurate estimation of mean kurtosis can be realised within a clinically feasible diffusion-weighted magnetic resonance imaging data acquisition time.

Impact of truncating diffusion MRI scans on diffusional kurtosis imaging

OBJECTIVE

Diffusional kurtosis imaging (DKI) extends diffusion tensor imaging (DTI), characterizing non-Gaussian diffusion effects but requires longer acquisition times. To ensure the robustness of DKI parameters, data acquisition ordering should be optimized allowing for scan interruptions or shortening. Three methodologies were used to examine how reduced diffusion MRI scans impact DKI histogram-metrics: 1) the electrostatic repulsion model (OptEEM); 2) spherical codes (OptSC); 3) random (RandomTRUNC).

MATERIALS AND METHODS

Pre-acquired diffusion multi-shell data from 14 female healthy volunteers (29±5 years) were used to generate reordered data. For each strategy, subsets containing different amounts of the full dataset were generated. The subsampling effects were assessed on histogram-based DKI metrics from tract-based spatial statistics (TBSS) skeletonized maps. To evaluate each subsampling method on simulated data at different SNRs and the influence of subsampling on in vivo data, we used a 3-way and 2-way repeated measures ANOVA, respectively.

RESULTS

Simulations showed that subsampling had different effects depending on DKI parameter, with fractional anisotropy the most stable (up to 5% error) and radial kurtosis the least stable (up to 26% error). RandomTRUNC performed the worst while the others showed comparable results. Furthermore, the impact of subsampling varied across distinct histogram characteristics, the peak value the least affected (OptEEM: up to 5% error; OptSC: up to 7% error) and peak height (OptEEM: up to 8% error; OptSC: up to 11% error) the most affected.

CONCLUSION

The impact of truncation depends on specific histogram-based DKI metrics. The use of a strategy for optimizing the acquisition order is advisable to improve DKI robustness to exam interruptions.

Progression/Recurrence of Meningioma: An Imaging Review Based on Magnetic Resonance Imaging

Meningiomas are the most common primary central nervous system tumors. The preferred treatment is maximum safe resection, and the heterogeneity of meningiomas results in a variable prognosis. Progression/recurrence (P/R) can occur at any grade of meningioma and is a common adverse outcome after surgical treatment and a major cause of postoperative rehospitalization, secondary surgery, and mortality. Early prediction of P/R plays an important role in postoperative management, further adjuvant therapy, and follow-up of patients. Therefore, it is essential to thoroughly analyze the heterogeneity of meningiomas and predict postoperative P/R with the aid of noninvasive preoperative imaging. In recent years, the development of advanced magnetic resonance imaging technology and machine learning has provided new insights into noninvasive preoperative prediction of meningioma P/R, which helps to achieve accurate prediction of meningioma P/R. This narrative review summarizes the current research on conventional magnetic resonance imaging, functional magnetic resonance imaging, and machine learning in predicting meningioma P/R. We further explore the significance of tumor microenvironment in meningioma P/R, linking imaging features with tumor microenvironment to comprehensively reveal tumor heterogeneity and provide new ideas for future research.

Spatiotemporal Patterns of White Matter Maturation after Pre-Adolescence: A Diffusion Kurtosis Imaging Study

Diffusion tensor imaging (DTI) enables the assessment of changes in brain tissue microstructure during maturation and ageing. In general, patterns of cerebral maturation and decline render non-monotonic lifespan trajectories of DTI metrics with age, and, importantly, the rate of microstructural changes is heterochronous for various white matter fibres. Recent studies have demonstrated that diffusion kurtosis imaging (DKI) metrics are more sensitive to microstructural changes during ageing compared to those of DTI. In a previous work, we demonstrated that the Cohen's d of mean diffusional kurtosis (dMK) represents a useful biomarker for quantifying maturation heterochronicity. However, some inferences on the maturation grades of different fibre types, such as association, projection, and commissural, were of a preliminary nature due to the insufficient number of fibres considered. Hence, the purpose of this follow-up work was to further explore the heterochronicity of microstructural maturation between pre-adolescence and middle adulthood based on DTI and DKI metrics. Using the effect size of the between-group parametric changes and Cohen's d, we observed that all commissural fibres achieved the highest level of maturity, followed by the majority of projection fibres, while the majority of association fibres were the least matured. We also demonstrated that dMK strongly correlates with the maxima or minima of the lifespan curves of DTI metrics. Furthermore, our results provide substantial evidence for the existence of spatial gradients in the timing of white matter maturation. In conclusion, our data suggest that DKI provides useful biomarkers for the investigation of maturation spatial heterogeneity and heterochronicity.

Maximizing SNR per unit time in diffusion MRI with multiband T-Hex spirals

PURPOSE

The characterization of tissue microstructure using diffusion MRI (dMRI) signals is rapidly evolving, with increasing sophistication of signal representations and microstructure models. However, this progress often requires signals to be acquired with very high b-values (e.g., b > 30 ms/μm2 ), along many directions, and using multiple b-values, leading to long scan times and extremely low SNR in dMRI images. The purpose of this work is to boost the SNR efficiency of dMRI by combining three particularly efficient spatial encoding techniques and utilizing a high-performance gradient system (Gmax  ≤ 300 mT/m) for efficient diffusion encoding.

METHODS

Spiral readouts, multiband imaging, and sampling on tilted hexagonal grids (T-Hex) are combined and implemented on a 3T MRI system with ultra-strong gradients. Image reconstruction is performed through an iterative cg-SENSE algorithm incorporating static off-resonance distributions and field dynamics as measured with an NMR field camera. Additionally, T-Hex multiband is combined with a more conventional EPI-readout and compared with state-of-the-art blipped-CAIPIRINHA sampling. The advantage of the proposed approach is furthermore investigated for clinically available gradient performance and diffusion kurtosis imaging.

RESULTS

High fidelity in vivo images with b-values up to 40 ms/μm2 are obtained. The approach provides superior SNR efficiency over other state-of-the-art multiband diffusion readout schemes.

CONCLUSION

The demonstrated gains hold promise for the widespread dissemination of advanced microstructural scans, especially in clinical populations.

Value of diffusion kurtosis MR imaging and conventional diffusion weighed imaging for evaluating response to first-line chemotherapy in unresectable pancreatic cancer

OBJECTIVE

To investigate the diagnostic value of diffusion kurtosis magnetic resonance imaging (DKI) and conventional diffusion-weighted imaging (DWI) for evaluating the response to first-line chemotherapy in unresectable pancreatic cancer.

MATERIALS AND METHODS

We retrospectively analyzed 21 patients with clinically and pathologically confirmed unresected pancreatic cancer who received palliative chemotherapy. Three-tesla MRI examinations containing DWI sequences with b values of 0, 100, 700, 1400, and 2100 s/mm2 were performed before and after chemotherapy. Parameters included the apparent diffusion coefficient (ADC), mean diffusion coefficient (MD), and mean diffusional kurtosis (MK). The performances of the DWI and DKI parameters in distinguishing the response to chemotherapy were evaluated by the area under the curve (AUC) of the receiver operating characteristic (ROC) curve. Overall survival (OS) was calculated from the date of first treatment to the date of death or the latest follow-up date.

RESULTS

The ADCchange and MDchange were significantly higher in the responding group (PR group) than in the nonresponding group (non-PR group) (ADCchange: 0.21 ± 0.05 vs. 0.11 ± 0.09, P = 0.02; MDchange: 0.37 ± 0.24 vs. 0.10 ± 0.12, P = 0.002). No statistical significance was shown when comparing ADCpre, ADCpost, MKpre, MKpost, MKchange, MDpre, and MDpost between the PR and non-PR groups. The ROC curve analysis indicated that MDchange (AUC = 0.898, cutoff value = 0.7143) performed better than ADCchange (AUC = 0.806, cutoff value = 0.1369) in predicting the response to chemotherapy.

CONCLUSION

The ADCchange and MDchange demonstrated strong potential for evaluating the response to chemotherapy in unresectable pancreatic cancer. The MDchange showed higher specificity in the classification of PR and non-PR than the ADCchange. Other parameters, including ADCpre, ADCpost, MKpre, MKpost, MKchange, MDpre, and MDpost, are not suitable for response evaluation. The combined model SUMchange demonstrated superior performance compared to the individual DWI and DKI models. Further experiments are needed to evaluate the potential of DWI and DKI parameters in predicting the prognosis of patients with unresectable pancreatic cancer.

Diffusion tensor and kurtosis imaging reveal microstructural changes in the trigeminal nerves of patients with trigeminal neuralgia

OBJECTIVES

To evaluate the use of diffusion tensor imaging (DTI) and diffusion kurtosis imaging (DKI) for detection of microstructural changes in the trigeminal nerves of trigeminal neuralgia (TN) patients.

METHODS

Forty TN patients and 40 healthy controls were examined using 3 T magnetic resonance imaging (MRI) to evaluate DTI and DKI parameters in trigeminal nerves. One-way ANOVA was used to test the differences in age, sex, and DTI and DKI parameters between the TN-affected sides, TN-unaffected sides, and controls. For parameters with inter-group differences, pairwise comparisons were performed. Then, the difference ratios (DRs) of the parameters with statistical differences were calculated and used for the receiver operating characteristic (ROC) analysis to assess their diagnostic performance. In addition, for the DTI and DKI parameter values with differences, we used pure DTI and DKI values to perform the ROC analysis.

RESULTS

Compared to the unaffected sides and controls, the FA, MK, and Kr of the affected sides of TN patients were significantly reduced, while ADC was significantly increased (p < 0.05). The diagnostic efficiency of the FA DRs (AUC: 0.974; cutoff value: 0.038; sensitivity: 100%; specificity: 95.0%) was the highest among all DTI and DKI parameters. The DRs of FA and MK more efficiently diagnosed TN than pure FA and MK values.

CONCLUSIONS

DTI and DKI allowed detection of microstructural changes in TN-affected trigeminal nerves. FA DR was the best independent predictor of microstructural changes in TN.

CLINICAL RELEVANCE STATEMENT

Both DTI and DKI can be used for noninvasive quantitative evaluation of the changes in the microstructure of the cisternal segment of the cranial nerves in clinical practice.

KEY POINTS

• Diffusion tensor imaging (DTI) can be used to evaluate the in vivo integrity of white matter and nerve fiber pathway. • Diffusion kurtosis imaging (DKI) has been shown to be considerable sensitive to microstructural changes. • DTI combined with DKI can comprehensively evaluate the status of the TN-affected trigeminal nerve.

Mutual-information based optimal experimental design for hyperpolarized [Formula: see text]C-pyruvate MRI

A key parameter of interest recovered from hyperpolarized (HP) MRI measurements is the apparent pyruvate-to-lactate exchange rate, [Formula: see text], for measuring tumor metabolism. This manuscript presents an information-theory-based optimal experimental design approach that minimizes the uncertainty in the rate parameter, [Formula: see text], recovered from HP-MRI measurements. Mutual information is employed to measure the information content of the HP measurements with respect to the first-order exchange kinetics of the pyruvate conversion to lactate. Flip angles of the pulse sequence acquisition are optimized with respect to the mutual information. A time-varying flip angle scheme leads to a higher parameter optimization that can further improve the quantitative value of mutual information over a constant flip angle scheme. However, the constant flip angle scheme, 35 and 28 degrees for pyruvate and lactate measurements, leads to an accuracy and precision comparable to the variable flip angle schemes obtained from our method. Combining the comparable performance and practical implementation, optimized pyruvate and lactate flip angles of 35 and 28 degrees, respectively, are recommended.

Optimal experimental design and estimation for q-space trajectory imaging

Tensor-valued diffusion encoding facilitates data analysis by q-space trajectory imaging. By modeling the diffusion signal of heterogeneous tissues with a diffusion tensor distribution (DTD) and modulating the encoding tensor shape, this novel approach allows disentangling variations in diffusivity from microscopic anisotropy, orientation dispersion, and mixtures of multiple isotropic diffusivities. To facilitate the estimation of the DTD parameters, a parsimonious acquisition scheme coupled with an accurate and precise estimation of the DTD is needed. In this work, we create two precision-optimized acquisition schemes: one that maximizes the precision of the raw DTD parameters, and another that maximizes the precision of the scalar measures derived from the DTD. The improved precision of these schemes compared to a naïve sampling scheme is demonstrated in both simulations and real data. Furthermore, we show that the weighted linear least squares (WLLS) estimator that uses the squared reciprocal of the noisy signal as weights can be biased, whereas the iteratively WLLS estimator with the squared reciprocal of the predicted signal as weights outperforms the conventional unweighted linear LS and nonlinear LS estimators in terms of accuracy and precision. Finally, we show that the use of appropriate constraints can considerably increase the precision of the estimator with only a limited decrease in accuracy.

Mean kurtosis-Curve (MK-Curve) correction improves the test-retest reproducibility of diffusion kurtosis imaging at 3 T

Diffusion kurtosis imaging (DKI) is applied to gain insights into the microstructural organization of brain tissues. However, the reproducibility of DKI outside brain white matter, particularly in combination with advanced estimation to remedy its noise sensitivity, remains poorly characterized. Therefore, in this study, we investigated the variability and reliability of DKI metrics while correcting implausible values with a fit method called mean kurtosis (MK)-Curve. A total of 10 volunteers (four women; age: 41.4 ± 9.6 years) were included and underwent two MRI examinations of the brain. The images were acquired on a clinical 3-T scanner and included a T1-weighted image and a diffusion sequence with multiple diffusion weightings suitable for DKI. Region of interest analysis of common kurtosis and tensor metrics derived with the MK-Curve DKI fit was performed, including intraclass correlation (ICC) and Bland-Altman (BA) plot statistics. A p value of less than 0.05 was considered statistically significant. The analyses showed good to excellent agreement of both kurtosis tensor- and diffusion tensor-derived MK-Curve-corrected metrics (ICC values: 0.77-0.98 and 0.87-0.98, respectively), with the exception of two DKI-derived metrics (axial kurtosis in the cortex: ICC = 0.68, and radial kurtosis in deep gray matter: ICC = 0.544). Non-MK-Curve-corrected kurtosis tensor-derived metrics ranged from 0.01 to 0.52 and diffusion tensor-derived metrics from 0.06 to 0.66, indicating poor to moderate reliability. No structural bias was observed in the BA plots for any of the diffusion metrics. In conclusion, MK-Curve-corrected DKI metrics of the human brain can be reliably acquired in white and gray matter at 3 T and DKI metrics have good to excellent agreement in a test-retest setting.

Tuned bipolar oscillating gradients for mapping frequency dispersion of diffusion kurtosis in the human brain

PURPOSE

Oscillating gradient spin-echo (OGSE) sequences have demonstrated an ability to probe time-dependent microstructural features, although they often suffer from low SNR due to increased TEs. In this work we introduce frequency-tuned bipolar (FTB) gradients as a variation of oscillating gradients with reduced TE and demonstrate their utility by mapping the frequency dispersion of kurtosis in human subjects.

METHODS

An FTB oscillating gradient waveform is presented that provides encoding of 1.5 net oscillation periods, thereby reducing the TE of the acquisition. Simulations were performed to determine an optimal protocol based on the SNR of kurtosis frequency dispersion-defined as the difference in kurtosis between pulsed and oscillating gradient acquisitions. Healthy human subjects were scanned at 7T using pulsed gradient and an optimized 23 Hz FTB protocol, which featured a maximum b-value of 2500 s/mm² . In addition, to directly compare existing methods, measurements using traditional cosine OGSE were also acquired.

RESULTS

FTB oscillating gradients demonstrated equivalent frequency-dependent diffusion measurements compared with cosine-modulated OGSE while enabling a significant reduction in TE. Optimization and in vivo results suggest that FTB gradients provide increased SNR of kurtosis dispersion maps compared with traditional cosine OGSE. The optimized FTB gradient protocol demonstrated consistent reductions in apparent kurtosis values and increased diffusivity in generated frequency dispersion maps.

CONCLUSIONS

This work presents an alternative to traditional cosine OGSE sequences, enabling more time-efficient acquisitions of frequency-dependent diffusion quantities as demonstrated through in vivo kurtosis frequency dispersion maps.

A correlational study between microstructural, macrostructural and functional age-related changes in the human visual cortex

Age-related changes in the human brain can be investigated from either structural or functional perspectives. Analysis of structural and functional age-related changes throughout the lifespan may help to understand the normal brain development process and monitor the structural and functional pathology of the brain. This study, combining dedicated electroencephalography (EEG) and magnetic resonance imaging (MRI) approaches in adults (20-78 years), highlights the complex relationship between micro/macrostructural properties and the functional responses to visual stimuli. Here, we aimed to relate age-related changes of the latency of visual evoked potentials (VEPs) to micro/macrostructural indexes and find any correlation between micro/macrostructural features, as well. We studied age-related structural changes in the brain, by using the MRI and diffusion-weighted imaging (DWI) as preferred imaging methods for extracting brain macrostructural parameters such as the cortical thickness, surface area, folding and curvature index, gray matter volume, and microstructural parameters such as mean diffusivity (MD), radial diffusivity (RD) and axial diffusivity (AD). All the mentioned features were significantly correlated with age in V1 and V2 regions of the visual cortex. Furthermore, we highlighted, negative correlations between structural features extracted from T1-weighted images and DWI. The latency and amplitude of the three dominants peaks (C1, P1, N1) of the VEP were considered as the brain functional features to be examined for correlation with age and structural features of the corresponding age. We observed significant correlations between mean C1 latency and GM volume averaged in V1 and V2. In hierarchical regression analysis, the structural index did not contribute to significant variance in the C1 latency after regressing out the effect of age. However, the age explained significant variance in the model after regressing out the effect of structural feature.

T2* quantification using multi-echo gradient echo sequences: a comparative study of different readout gradients

To quantify T2*, multiple echoes are typically acquired with a multi-echo gradient echo sequence using either monopolar or bipolar readout gradients. The use of bipolar readout gradients achieves a shorter echo spacing time, enabling the acquisition of a larger number of echoes in the same scan time. However, despite their relative time efficiency and the potential for more accurate quantification, a comparative investigation of these readout gradients has not yet been addressed. This work aims to compare the performance of monopolar and bipolar readout gradients for T2* quantification. The differences in readout gradients were theoretically investigated with a Cramér-Rao lower bound and validated with computer simulations with respect to the various imaging parameters (e.g., flip angle, TR, TE, TE range, and BW). The readout gradients were then compared at 3 T using phantom and in vivo experiments. The bipolar readout gradients provided higher precision than monopolar readout gradients in both computer simulations and experimental results. The difference between the two readout gradients increased for a lower SNR and smaller TE range, consistent with the prediction made using Cramér-Rao lower bound. The use of bipolar readout gradients is advantageous for regions or situations where a lower SNR is expected or a shorter acquisition time is required.

Diffusion MRI with pulsed and free gradient waveforms: Effects of restricted diffusion and exchange

Monitoring time dependence with diffusion MRI yields observables sensitive to compartment sizes (restricted diffusion) and membrane permeability (water exchange). However, restricted diffusion and exchange have opposite effects on the diffusion-weighted signal, which can lead to errors in parameter estimates. In this work, we propose a signal representation that incorporates the effects of both restricted diffusion and exchange up to second order in b-value and is compatible with gradient waveforms of arbitrary shape. The representation features mappings from a gradient waveform to two scalars that separately control the sensitivity to restriction and exchange. We demonstrate that these scalars span a two-dimensional space that can be used to choose waveforms that selectively probe restricted diffusion or exchange, eliminating the correlation between the two phenomena. We found that waveforms with specific but unconventional shapes provide an advantage over conventional pulsed and oscillating gradient acquisitions. We also show that parametrization of waveforms into a two-dimensional space can be used to understand protocols from other approaches that probe restricted diffusion and exchange. For example, we found that the variation of mixing time in filter-exchange imaging corresponds to variation of our exchange-weighting scalar at a fixed value of the restriction-weighting scalar. The proposed signal representation was evaluated using Monte Carlo simulations in identical parallel cylinders with hexagonal and random packing as well as parallel cylinders with gamma-distributed radii. Results showed that the approach is sensitive to sizes in the interval 4-12 μm and exchange rates in the simulated range of 0 to 20 s - 1 , but also that there is a sensitivity to the extracellular geometry. The presented theory constitutes a simple and intuitive description of how restricted diffusion and exchange influence the signal as well as a guide to protocol design capable of separating the two effects.

Optimization of quasi-diffusion magnetic resonance imaging for quantitative accuracy and time-efficient acquisition

PURPOSE

Quasi-diffusion MRI (QDI) is a novel quantitative technique based on the continuous time random walk model of diffusion dynamics. QDI provides estimates of the diffusion coefficient, <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mspace/> <mml:msub><mml:mi>D</mml:mi> <mml:mrow><mml:mn>1</mml:mn> <mml:mo>,</mml:mo> <mml:mn>2</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ {D}_{1,2} $$</mml:annotation></mml:semantics> </mml:math> in mm²  s^-1 and a fractional exponent, <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mi>α</mml:mi></mml:mrow> <mml:annotation>$$ \upalpha $$</mml:annotation></mml:semantics> </mml:math> , defining the non-Gaussianity of the diffusion signal decay. Here, the b-value selection for rapid clinical acquisition of QDI tensor imaging (QDTI) data is optimized.

METHODS

Clinically appropriate QDTI acquisitions were optimized in healthy volunteers with respect to a multi-b-value reference (MbR) dataset comprising 29 diffusion-sensitized images arrayed between <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mi>b</mml:mi> <mml:mo>=</mml:mo> <mml:mn>0</mml:mn></mml:mrow> <mml:annotation>$$ b=0 $$</mml:annotation></mml:semantics> </mml:math> and 5000 s mm^-2 . The effects of varying maximum b-value ( <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msub><mml:mi>b</mml:mi> <mml:mi>max</mml:mi></mml:msub> </mml:mrow> <mml:annotation>$$ {b}_{\mathrm{max}} $$</mml:annotation></mml:semantics> </mml:math> ), number of b-value shells, and the effects of Rician noise were investigated.

RESULTS

QDTI measures showed <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msub><mml:mi>b</mml:mi> <mml:mi>max</mml:mi></mml:msub> </mml:mrow> <mml:annotation>$$ {b}_{\mathrm{max}} $$</mml:annotation></mml:semantics> </mml:math> dependence, most significantly for <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mi>α</mml:mi></mml:mrow> <mml:annotation>$$ \upalpha $$</mml:annotation></mml:semantics> </mml:math> in white matter, which monotonically decreased with higher <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msub><mml:mi>b</mml:mi> <mml:mi>max</mml:mi></mml:msub> </mml:mrow> <mml:annotation>$$ {b}_{\mathrm{max}} $$</mml:annotation></mml:semantics> </mml:math> leading to improved tissue contrast. Optimized 2 b-value shell acquisitions showed small systematic differences in QDTI measures relative to MbR values, with overestimation of <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mspace/> <mml:mspace/> <mml:msub><mml:mi>D</mml:mi> <mml:mrow><mml:mn>1</mml:mn> <mml:mo>,</mml:mo> <mml:mn>2</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ \kern0.50em {D}_{1,2} $$</mml:annotation></mml:semantics> </mml:math> and underestimation of <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mi>α</mml:mi></mml:mrow> <mml:annotation>$$ \upalpha $$</mml:annotation></mml:semantics> </mml:math> in white matter, and overestimation of <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msub><mml:mi>D</mml:mi> <mml:mrow><mml:mn>1</mml:mn> <mml:mo>,</mml:mo> <mml:mn>2</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ {D}_{1,2} $$</mml:annotation></mml:semantics> </mml:math> and <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mi>α</mml:mi></mml:mrow> <mml:annotation>$$ \upalpha $$</mml:annotation></mml:semantics> </mml:math> anisotropies in gray and white matter. Additional shells improved the accuracy, precision, and reliability of QDTI estimates with 3 and 4 shells at <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msub><mml:mi>b</mml:mi> <mml:mi>max</mml:mi></mml:msub> <mml:mo>=</mml:mo> <mml:mn>5000</mml:mn></mml:mrow> <mml:annotation>$$ {b}_{\mathrm{max}}=5000 $$</mml:annotation></mml:semantics> </mml:math>  s mm^-2 , and 4 b-value shells at <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msub><mml:mi>b</mml:mi> <mml:mi>max</mml:mi></mml:msub> <mml:mo>=</mml:mo> <mml:mn>3960</mml:mn></mml:mrow> <mml:annotation>$$ {b}_{\mathrm{max}}=3960 $$</mml:annotation></mml:semantics> </mml:math>  s mm^-2 , providing minimal bias in <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msub><mml:mi>D</mml:mi> <mml:mrow><mml:mn>1</mml:mn> <mml:mo>,</mml:mo> <mml:mn>2</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ {D}_{1,2} $$</mml:annotation></mml:semantics> </mml:math> and <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mi>α</mml:mi></mml:mrow> <mml:annotation>$$ \upalpha $$</mml:annotation></mml:semantics> </mml:math> compared to the MbR.

CONCLUSION

A highly detailed optimization of non-Gaussian dMRI for in vivo brain imaging was performed. QDI provided robust parameterization of non-Gaussian diffusion signal decay in clinically feasible imaging times with high reliability, accuracy, and precision of QDTI measures.

Optimization of scan parameters to reduce acquisition time for RESOLVE-based diffusion kurtosis imaging (DKI) in nasopharyngeal carcinoma (NPC)

OBJECTIVE

To shorten acquisition time of readout segmentation of long variable echo trains (RESOLVE)-based diffusion kurtosis imaging (DKI) via Readout Partial Fourier (RPF) and b-value combinations.

METHODS

The RESOLVE-based DKI images of 38 patients with nasopharyngeal carcinoma (NPC) were prospectively enrolled. For RESOLVE-based DKI images with 5/8 RPF and without RPF, objective and subjective evaluations of image quality were performed. A total of nine groups with different b-value combinations were simulated, and the influence of different b-value combinations for RESOLVE-RPF-based DKI sequences was assessed using the intraclass correlation coefficient (ICC).

RESULTS

The mean values of signal-to-noise ratio (SNR) and contrast-to-noise ratio (CNR) in DKI images without RPF were higher than those with 5/8 RPF (252.9 ± 77.7 vs 247.3 ± 85.5 and 5.8 ± 2.8 vs 5.4 ± 2.3, respectively), but not significantly (p = 0.460 and p = 0.180, respectively). In comparing the ICCs between nine groups of different b-value combinations in RESOLVE-RPF-based DKI, group (200, 800, 2000 s/mm2), group (200, 400, 800, 2000 s/mm2) and group (200, 800, 1500, 2000 s/mm2) were not significantly different (p > 0.001) and showed excellent agreement (0.81-1.00) with that of group (200, 400, 800, 1500, 2000 s/mm2). Using b-value optimization and RPF technology, the group with RPF (200, 400, 800, 2000 s/mm2) showed a 56% reduced scanning compared with the group without RPF (200, 400, 800, 1500, 2000 s/mm2; 3 min 46 s vs 8 min 31 s, respectively).

CONCLUSION

DKI with RPF did not significantly affect image quality, but both RPF and different b-value combinations can affect the scanning time. The combination of RPF and b-value optimization can ensure the stability of DKI parameters and reduce the scanning time by 56%.

ADVANCES IN KNOWLEDGE

This work is to optimize scan parameters, e.g. RPF and b-value combinations, to reduce acquisition time for RESOLVE-based DKI in NPC. To our knowledge, the effect of RESOLVE-RPF and b-value combinations on DKI has not been reported.

The COMPLETE trial: HolistiC early respOnse assessMent for oroPharyngeaL cancEr paTiEnts; Protocol for an observational study

INTRODUCTION

The locoregional failure (LRF) rate in human papilloma virus (HPV)-negative oropharyngeal squamous cell carcinoma (OPSCC) remains disappointingly high and toxicity is substantial. Response prediction prior to or early during treatment would provide opportunities for personalised treatment. Currently, there are no accurate predictive models available for correct OPSCC patient selection. Apparently, the pivotal driving forces that determine how a OPSCC responds to treatment, have yet to be elucidated. Therefore, the holistiC early respOnse assessMent for oroPharyngeaL cancer paTiEnts study focuses on a holistic approach to gain insight in novel potential prognostic biomarkers, acquired before and early during treatment, to predict response to treatment in HPV-negative patients with OPSCC.

METHODS AND ANALYSIS

This single-centre prospective observational study investigates 60 HPV-negative patients with OPSCC scheduled for primary radiotherapy (RT) with cisplatin or cetuximab, according to current clinical practice. A holistic approach will be used that aims to map the macroscopic (with Intra Voxel Incoherent Motion Diffusion Kurtosis Imaging (IVIM-DKI); before, during, and 3 months after RT), microscopic (with biopsies of the primary tumour acquired before treatment and irradiated ex vivo to assess radiosensitivity), and molecular landscape (with circulating tumour DNA (ctDNA) analysed before, during and 3 months after treatment). The main end point is locoregional control (LRC) 2 years after treatment. The primary objective is to determine whether a relative change in the mean of the diffusion coefficient D (an IVIM-DKI parameter) in the primary tumour early during treatment, improves the performance of a predictive model consisting of tumour volume only, for 2 years LRC after treatment. The secondary objectives investigate the potential of other IVIM-DKI parameters, ex vivo sensitivity characteristics, ctDNA, and combinations thereof as potential novel prognostic markers.

ETHICS AND DISSEMINATION

The study was approved by the Medical Ethical Committee of Erasmus Medical Center. The main results of the trial will be presented in international meetings and medical journals.

TRIAL REGISTRATION NUMBER

NL8458.

Improved diffusion parameter estimation by incorporating T2 relaxation properties into the DKI-FWE model

The free water elimination (FWE) model and its kurtosis variant (DKI-FWE) can separate tissue and free water signal contributions, thus providing tissue-specific diffusional information. However, a downside of these models is that the associated parameter estimation problem is ill-conditioned, necessitating the use of advanced estimation techniques that can potentially bias the parameter estimates. In this work, we propose the T₂-DKI-FWE model that exploits the T₂ relaxation properties of both compartments, thereby better conditioning the parameter estimation problem and providing, at the same time, an additional potential biomarker (the T₂ of tissue). In our approach, the T₂ of tissue is estimated as an unknown parameter, whereas the T₂ of free water is assumed known a priori and fixed to a literature value (1573 ms). First, the error propagation of an erroneous assumption on the T₂ of free water is studied. Next, the improved conditioning of T₂-DKI-FWE compared to DKI-FWE is illustrated using the Cramér-Rao lower bound matrix. Finally, the performance of the T₂-DKI-FWE model is compared to that of the DKI-FWE and T₂-DKI models on both simulated and real datasets. The error due to a biased approximation of the T₂ of free water was found to be relatively small in various diffusion metrics and for a broad range of erroneous assumptions on its underlying ground truth value. Compared to DKI-FWE, using the T₂-DKI-FWE model is beneficial for the identifiability of the model parameters. Our results suggest that the T₂-DKI-FWE model can achieve precise and accurate diffusion parameter estimates, through effective reduction of free water partial volume effects and by using a standard nonlinear least squares approach. In conclusion, incorporating T₂ relaxation properties into the DKI-FWE model improves the conditioning of the model fitting, while only requiring an acquisition scheme with at least two different echo times.

Generalisation of continuous time random walk to anomalous diffusion MRI models with an age-related evaluation of human corpus callosum

Diffusion MRI measures of the human brain provide key insight into microstructural variations across individuals and into the impact of central nervous system diseases and disorders. One approach to extract information from diffusion signals has been to use biologically relevant analytical models to link millimetre scale diffusion MRI measures with microscale influences. The other approach has been to represent diffusion as an anomalous transport process and infer microstructural information from the different anomalous diffusion equation parameters. In this study, we investigated how parameters of various anomalous diffusion models vary with age in the human brain white matter, particularly focusing on the corpus callosum. We first unified several established anomalous diffusion models (the super-diffusion, sub-diffusion, quasi-diffusion and fractional Bloch-Torrey models) under the continuous time random walk modelling framework. This unification allows a consistent parameter fitting strategy to be applied from which meaningful model parameter comparisons can be made. We then provided a novel way to derive the diffusional kurtosis imaging (DKI) model, which is shown to be a degree two approximation of the sub-diffusion model. This link between the DKI and sub-diffusion models led to a new robust technique for generating maps of kurtosis and diffusivity using the sub-diffusion parameters β_SUB and D_SUB. Superior tissue contrast is achieved in kurtosis maps based on the sub-diffusion model. 7T diffusion weighted MRI data for 65 healthy participants in the age range 19-78 years was used in this study. Results revealed that anomalous diffusion model parameters α and β have shown consistent positive correlation with age in the corpus callosum, indicating α and β are sensitive to tissue microstructural changes in ageing.

A Joint Framework for Denoising and Estimating Diffusion Kurtosis Tensors Using Multiple Prior Information

Diffusion kurtosis imaging (DKI) has been shown to be valuable in a wide range of neuroscientific and clinical applications. However, reliable estimation of DKI tensors is often compromised by noise, especially for the kurtosis tensor (KT). Here, we propose a joint denoising and estimating framework that integrates multiple sources of prior information, including nonlocal structural self-similarity (NSS), local spatial smoothness (LSS), physical relevance (PR) of the DKI model, and noise characteristics of magnitude diffusion MRI (dMRI) images for improved estimation of DKI tensors. The local and nonlocal spatial smoothing constraints are complementary to each other, making the proposed framework highly effective in reducing the noise fluctuations on DKI tensors, especially KT. As an additional refinement, we propose to impose a physically relevant constraint within our joint denoising and estimation framework. We further adopt the first-moment noise-corrected fitting model (M¹NCM) to remove the noncentral χ -distribution noise bias. The effectiveness of integrating multiple sources of priors into the joint framework is verified by comparing the proposed M¹NCM-NSS-LSS-PR method with various versions of M¹NCM-based estimators and two state-of-the-art methods. Results show that the proposed method outperformed the compared methods in simulations and in-vivo dMRI datasets of both spatially stationary and nonstationary noise distributions. The in-vivo experiments also show that the proposed M¹NCM-NSS-LSS-PR method was robust to the number of diffusion directions.

Spatiotemporal characterisation of ischaemic lesions in transient stroke animal models using diffusion free water elimination and mapping MRI with echo time dependence

BACKGROUND AND PURPOSE

The excess fluid as a result of vasogenic oedema and the subsequent tissue cavitation obscure the microstructural characterisation of ischaemic tissue by conventional diffusion and relaxometry MRI. They lead to a pseudo-normalisation of the water diffusivity and transverse relaxation time maps in the subacute and chronic phases of stroke. Within the context of diffusion MRI, the free water elimination and mapping method (FWE) with echo time dependence has been proposed as a promising approach to measure the amount of free fluid in brain tissue robustly and to eliminate its biasing effect on other biomarkers. In this longitudinal study of transient middle cerebral artery occlusion (MCAo) in the rat brain, we investigated the use of FWE MRI with echo time dependence for the characterisation of the tissue microstructure and explored the potential of the free water fraction as a novel biomarker of ischaemic tissue condition.

METHODS

Adult rats received a transient MCAo. Diffusion- and transverse relaxation-weighted MRI experiments were performed longitudinally, pre-occlusion and on days 1, 3, 4, 5, 6, 7 and 10 after MCAo on four rats. Histology was performed for non-stroke and 1, 3 and 10 days after MCAo on three different rats at each time point.

RESULTS

The free water fraction was homogeneously increased in the ischaemic cortex one day after stroke. Between three and ten days after stroke, the core of the ischaemic tissue showed a progressive normalisation in the amount of free water, whereas the inner and outer border zones of the ischaemic cortex depicted a large, monotonous increase with time. The specific lesions in brain sections were verified by H&E and immunostaining. The tissue-specific diffusion and relaxometry MRI metrics in the ischaemic cortex were significantly different compared to their conventional counterpart.

CONCLUSIONS

Our results demonstrate that the free water fraction in FWE MRI with echo time dependence is a valuable biomarker, sensitive to the progressive degeneration in ischaemic tissue. We showed that part of the heterogeneity previously observed in conventional parameter maps can be accounted for by a heterogeneous distribution of free water in the tissue. Our results suggest that the temporal evolution of the free fluid fraction map at the core and inner border zone can be associated with the pathological changes linked to the evolution of vasogenic oedema. Namely, the homogeneous increase in free water one day after stroke and its tendency to normalise in the core of the ischaemic cortex starting three days after stroke, followed by a progressive increase in free water at the inner border zone from three to ten days after stroke. Finally, the monotonous increase in free fluid in the outer border zone of the cortex reflects the formation of fluid-filled cysts.

Fisher information for optimal planning of X-ray diffraction experiments

Fisher information is a powerful mathematical tool suitable for quantification of data `informativity' and optimization of the experimental setup and measurement conditions. Here, it is applied to X-ray diffraction and an informational approach to choosing the optimal measurement configuration is proposed. The core idea is maximization of the information which can be extracted from the measured data set by the selected analysis technique, over the sets of accessible reflections and measurement geometries. The developed approach is applied to high-resolution X-ray diffraction measurements and microstructure analysis of multilayer samples, and its efficiency and consistency are demonstrated with the results of more straightforward Monte Carlo simulations.

Task-driven assessment of experimental designs in diffusion MRI: A computational framework

This paper proposes a task-driven computational framework for assessing diffusion MRI experimental designs which, rather than relying on parameter-estimation metrics, directly measures quantitative task performance. Traditional computational experimental design (CED) methods may be ill-suited to experimental tasks, such as clinical classification, where outcome does not depend on parameter-estimation accuracy or precision alone. Current assessment metrics evaluate experiments' ability to faithfully recover microstructural parameters rather than their task performance. The method we propose addresses this shortcoming. For a given MRI experimental design (protocol, parameter-estimation method, model, etc.), experiments are simulated start-to-finish and task performance is computed from receiver operating characteristic (ROC) curves and associated summary metrics (e.g. area under the curve (AUC)). Two experiments were performed: first, a validation of the pipeline's task performance predictions against clinical results, comparing in-silico predictions to real-world ROC/AUC; and second, a demonstration of the pipeline's advantages over traditional CED approaches, using two simulated clinical classification tasks. Comparison with clinical datasets validates our method's predictions of (a) the qualitative form of ROC curves, (b) the relative task performance of different experimental designs, and (c) the absolute performance (AUC) of each experimental design. Furthermore, we show that our method outperforms traditional task-agnostic assessment methods, enabling improved, more useful experimental design. Our pipeline produces accurate, quantitative predictions of real-world task performance. Compared to current approaches, such task-driven assessment is more likely to identify experimental designs that perform well in practice. Our method is not limited to diffusion MRI; the pipeline generalises to any task-based quantitative MRI application, and provides the foundation for developing future task-driven end-to end CED frameworks.

Brain Structural and Functional Connectivity: A Review of Combined Works of Diffusion Magnetic Resonance Imaging and Electro-Encephalography

Implications of structural connections within and between brain regions for their functional counterpart are timely points of discussion. White matter microstructural organization and functional activity can be assessed in unison. At first glance, however, the corresponding findings appear variable, both in the healthy brain and in numerous neuro-pathologies. To identify consistent associations between structural and functional connectivity and possible impacts for the clinic, we reviewed the literature of combined recordings of electro-encephalography (EEG) and diffusion-based magnetic resonance imaging (MRI). It appears that the strength of event-related EEG activity increases with increased integrity of structural connectivity, while latency drops. This agrees with a simple mechanistic perspective: the nature of microstructural white matter influences the transfer of activity. The EEG, however, is often assessed for its spectral content. Spectral power shows associations with structural connectivity that can be negative or positive often dependent on the frequencies under study. Functional connectivity shows even more variations, which are difficult to rank. This might be caused by the diversity of paradigms being investigated, from sleep and resting state to cognitive and motor tasks, from healthy participants to patients. More challenging, though, is the potential dependency of findings on the kind of analysis applied. While this does not diminish the principal capacity of EEG and diffusion-based MRI co-registration, it highlights the urgency to standardize especially EEG analysis.

Diffusion-Weighted Imaging: Recent Advances and Applications

Quantitative diffusion imaging techniques enable the characterization of tissue microstructural properties of the human brain "in vivo", and are widely used in neuroscientific and clinical contexts. In this review, we present the basic physical principles behind diffusion imaging and provide an overview of the current diffusion techniques, including standard and advanced techniques as well as their main clinical applications. Standard diffusion tensor imaging (DTI) offers sensitivity to changes in microstructure due to diseases and enables the characterization of single fiber distributions within a voxel as well as diffusion anisotropy. Nonetheless, its inability to represent complex intravoxel fiber topologies and the limited biological specificity of its metrics motivated the development of several advanced diffusion MRI techniques. For example, high-angular resolution diffusion imaging (HARDI) techniques enabled the characterization of fiber crossing areas and other complex fiber topologies in a single voxel and supported the development of higher-order signal representations aiming to decompose the diffusion MRI signal into distinct microstructure compartments. Biophysical models, often known by their acronym (e.g., CHARMED, WMTI, NODDI, DBSI, DIAMOND) contributed to capture the diffusion properties from each of such tissue compartments, enabling the computation of voxel-wise maps of axonal density and/or morphology that hold promise as clinically viable biomarkers in several neurological and neuroscientific applications; for example, to quantify tissue alterations due to disease or healthy processes. Current challenges and limitations of state-of-the-art models are discussed, including validation efforts. Finally, novel diffusion encoding approaches (e.g., b-tensor or double diffusion encoding) may increase the biological specificity of diffusion metrics towards intra-voxel diffusion heterogeneity in clinical settings, holding promise in neurological applications.

MRI phantom for tissue simulation with respect to diffusion coefficient and kurtosis - Validation with injection of liposomal theranostics

This study presents gelatine-based and agar-based phantoms with an addition of glycerol, safflower oil, silicone oil and cellulose microcrystalline with a potential to cover the entire range of tissue diffusion coefficients and kurtosis values. Forty types of phantoms were prepared and examined for NMR relaxation times T₁ and T₂ and diffusional metrics D, K and ADC. Wide ranges of values of D (0.0003-0.0031 mm²s^-1), K (0.00-7.24) and ADC (0.0002-0.0031 mm²s^-1) were observed. Two of the phantoms closely mimic muscle and cortical gray matter with respect to water diffusion parameters. Although many of the presented phantoms display both D and K values within the range of human tissues, they match different tissues with respect to D and K. The imaging results for the gray matter simulating phantom injected with the liposomal solution, bear a resemblance to the particle size effect described in the literature. The phantoms presented in this work are simple in preparation and affordable tissue-simulating materials to be used primarily in development of diffusion kurtosis-based MRI methods and possibly in a preliminary assessment of MRI contrast agents. Further adjustments of the chemical compositions could potentially lead to development of new types of phantoms mimicking diffusional properties of more kinds of soft tissues.

Glioma-Specific Diffusion Signature in Diffusion Kurtosis Imaging

Purpose: This study aimed to assess the relationship between mean kurtosis (MK) and mean diffusivity (MD) values from whole-brain diffusion kurtosis imaging (DKI) parametric maps in preoperative magnetic resonance (MR) images from 2016 World Health Organization Classification of Tumors of the Central Nervous System integrated glioma groups. Methods: Seventy-seven patients with histopathologically confirmed treatment-naïve glioma were retrospectively assessed between 1 August 2013 and 30 October 2017. The area on scatter plots with a specific combination of MK and MD values, not occurring in the healthy brain, was labeled, and the corresponding voxels were visualized on the fluid-attenuated inversion recovery (FLAIR) images. Reversely, the labeled voxels were compared to those of the manually segmented tumor volume, and the Dice similarity coefficient was used to investigate their spatial overlap. Results: A specific combination of MK and MD values in whole-brain DKI maps, visualized on a two-dimensional scatter plot, exclusively occurs in glioma tissue including the perifocal infiltrative zone and is absent in tissue of the normal brain or from other intracranial compartments. Conclusions: A unique diffusion signature with a specific combination of MK and MD values from whole-brain DKI can identify diffuse glioma without any previous segmentation. This feature might influence artificial intelligence algorithms for automatic tumor segmentation and provide new aspects of tumor heterogeneity.

Gradient waveform design for tensor-valued encoding in diffusion MRI

Diffusion encoding along multiple spatial directions per signal acquisition can be described in terms of a b-tensor. The benefit of tensor-valued diffusion encoding is that it unlocks the 'shape of the b-tensor' as a new encoding dimension. By modulating the b-tensor shape, we can control the sensitivity to microscopic diffusion anisotropy which can be used as a contrast mechanism; a feature that is inaccessible by conventional diffusion encoding. Since imaging methods based on tensor-valued diffusion encoding are finding an increasing number of applications we are prompted to highlight the challenge of designing the optimal gradient waveforms for any given application. In this review, we first establish the basic design objectives in creating field gradient waveforms for tensor-valued diffusion MRI. We also survey additional design considerations related to limitations imposed by hardware and physiology, potential confounding effects that cannot be captured by the b-tensor, and artifacts related to the diffusion encoding waveform. Throughout, we discuss the expected compromises and tradeoffs with an aim to establish a more complete understanding of gradient waveform design and its impact on accurate measurements and interpretations of data.

The costs and benefits of estimating T1 of tissue alongside cerebral blood flow and arterial transit time in pseudo-continuous arterial spin labeling

Multi-post-labeling-delay pseudo-continuous arterial spin labeling (multi-PLD PCASL) allows for absolute quantification of the cerebral blood flow (CBF) as well as the arterial transit time (ATT). Estimating these perfusion parameters from multi-PLD PCASL data is a non-linear inverse problem, which is commonly tackled by fitting the single-compartment model (SCM) for PCASL, with CBF and ATT as free parameters. The longitudinal relaxation time of tissue T_1t is an important parameter in this model, as it governs the decay of the perfusion signal entirely upon entry in the imaging voxel. Conventionally, T_1t is fixed to a population average. This approach can cause CBF quantification errors, as T_1t can vary significantly inter- and intra-subject. This study compares the impact on CBF quantification, in terms of accuracy and precision, of either fixing T_1t , the conventional approach, or estimating it alongside CBF and ATT. It is shown that the conventional approach can cause a significant bias in CBF. Indeed, simulation experiments reveal that if T_1t is fixed to a value that is 10% off its true value, this may already result in a bias of 15% in CBF. On the other hand, as is shown by both simulation and real data experiments, estimating T_1t along with CBF and ATT results in a loss of CBF precision of the same order, even if the experiment design is optimized for the latter estimation problem. Simulation experiments suggest that an optimal balance between accuracy and precision of CBF estimation from multi-PLD PCASL data can be expected when using the two-parameter estimator with a fixed T_1t value between population averages of T_1t and the longitudinal relaxation time of blood T_1b .

An information theory model for optimizing quantitative magnetic resonance imaging acquisitions

Acquisition parameter selection is currently performed empirically for many quantitative MRI (qMRI) acquisitions. Tuning parameters for different scan times, tissues, and resolutions requires some amount of trial and error. There is an opportunity to quantitatively optimize these acquisition parameters in order to minimize variability of quantitative maps and post-processing techniques such as synthetic image generation. The objective of this work is to introduce and evaluate a quantitative method for selecting parameters that minimize image variability. An information theory framework was developed for this purpose and applied to a 3D-quantification using an interleaved Look-Locker acquisition sequence with T2 preparation pulse (3D-QALAS) signal model for qMRI. In this framework, mutual information is used to measure the information gained by a measurement as a function of acquisition parameters, quantifying the information content of potential acquisitions and allowing informed parameter selection. The information theory framework was tested on artificial data generated from a representative mathematical phantom, measurements acquired on a qMRI multiparametric imaging standard phantom, and in vivo measurements in a human brain. The phantom measurements showed that higher mutual information calculated by the model correlated with smaller coefficient of variation in the reconstructed parametric maps, and in vivo measurements demonstrated that information-based calibration of acquisition parameters resulted in a decrease in parametric map variability consistent with model predictions.

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.

A data-driven approach to optimising the encoding for multi-shell diffusion MRI with application to neonatal imaging

Diffusion MRI has the potential to provide important information about the connectivity and microstructure of the human brain during normal and abnormal development, noninvasively and in vivo. Recent developments in MRI hardware and reconstruction methods now permit the acquisition of large amounts of data within relatively short scan times. This makes it possible to acquire more informative multi-shell data, with diffusion sensitisation applied along many directions over multiple b-value shells. Such schemes are characterised by the number of shells acquired, and the specific b-value and number of directions sampled for each shell. However, there is currently no clear consensus as to how to optimise these parameters. In this work, we propose a means of optimising multi-shell acquisition schemes by estimating the information content of the diffusion MRI signal, and optimising the acquisition parameters for sensitivity to the observed effects, in a manner agnostic to any particular diffusion analysis method that might subsequently be applied to the data. This method was used to design the acquisition scheme for the neonatal diffusion MRI sequence used in the developing Human Connectome Project (dHCP), which aims to acquire high quality data and make it freely available to the research community. The final protocol selected by the algorithm, and currently in use within the dHCP, consists of 20 b=0 images and diffusion-weighted images at b = 400, 1000 and 2600 s/mm2 with 64, 88 and 128 directions per shell, respectively.

An optimal acquisition and post-processing pipeline for hybrid IVIM-DKI in head and neck

Purpose

To optimize the diffusion‐weighting b values and postprocessing pipeline for hybrid intravoxel incoherent motion diffusion kurtosis imaging in the head and neck region.

Methods

Optimized diffusion‐weighting b value sets ranging between 5 and 30 b values were constructed by optimizing the Cramér‐Rao lower bound of the hybrid intravoxel incoherent motion diffusion kurtosis imaging model. With this model, the perfusion fraction, pseudodiffusion coefficient, diffusion coefficient, and kurtosis were estimated. Sixteen volunteers were scanned with a reference b value set and 3 repeats of the optimized sets, of which 1 with volunteers swallowing on purpose.

The effects of (1) b value optimization and number of b values, (2) registration type (none vs. intervolume vs. intra‐ and intervolume registration), and (3) manual swallowing artifact rejection on the parameter precision were assessed.

Results

The SD was higher in the reference set for perfusion fraction, diffusion coefficient, and kurtosis by a factor of 1.7, 1.5, and 2.3 compared to the optimized set, respectively. A smaller SD (factor 0.7) was seen in pseudodiffusion coefficient. The sets containing 15, 20, and 30 b values had comparable repeatability in all parameters, except pseudodiffusion coefficient, for which set size 30 was worse. Equal repeatability for the registration approaches was seen in all parameters of interest. Swallowing artifact rejection removed the bias when present.

Conclusion

To achieve optimal hybrid intravoxel incoherent motion diffusion kurtosis imaging in the head and neck region, b value optimization and swallowing artifact image rejection are beneficial. The optimized set of 15 b values yielded the optimal protocol efficiency, with a precision comparable to larger b value sets and a 50% reduction in scan time.

Magnetic Resonance Diffusion Kurtosis Imaging versus Diffusion-Weighted Imaging in Evaluating the Pathological Grade of Hepatocellular Carcinoma

Purpose

To investigate the diagnostic efficacy of diffusion kurtosis imaging (DKI) and conventional diffusion-weighted imaging (DWI) for pathological grading.

Methods

From December 2015 to January 2017, consecutive patients suspected of having hepatocellular carcinoma (HCC) without prior treatment were prospectively enrolled in this study. MRI examinations were performed before surgical treatment. HCC patients confirmed by surgical pathology were included in the study. The mean diffusivity (MD) values, mean kurtosis (MK) values, and apparent diffusion coefficient (ADC) were calculated. The differences and correlations of these parameters among different pathological grades were analyzed. The diagnostic efficiency of DKI and DWI for predicting high-grade HCC was evaluated by receiver operating characteristic (ROC) curves. Logistic regression analyses were used to evaluate the predictive factors for pathological grade.

Results

A total of 128 patients (79 males and 49 females, age: 56.9±10.9 years, range, 32–80) with primary HCC were included: grade I: 22 (17.2%) patients, grade II: 37 (28.9%) patients, grade III: 43 (33.6%) patients, grade IV: 26 (20.3%) patients. The MK values of stage I, II, III, and IV were 0.86±0.13, 1.06±0.11, 1.27±0.17, and 1.57±0.13, respectively. The MK values were significantly higher in the high-grade group than in the low-grade group and were positively correlated with pathological grade (rho =0.7417, P<0.001). The MK value demonstrated a larger area under the curve (AUC), with a value of 0.93 than the MD value, which had an AUC of 0.815 (P<0.001), and ADC, which had an AUC of 0.662 (P=0.01). The MK value (>1.19), ADC (≤1.29×10^–3 mm²/s), and HBV (+) were independent predictors for the pathological grade of HCCs.

Conclusion

The MK values derived from DKI and the ADC values obtained from traditional DWI were more valuable than the MD values in predicting the histological grade of HCCs and could potentially guide clinical treatment before surgery.

Assessing White Matter Pathology in Early-Stage Parkinson Disease Using Diffusion MRI: A Systematic Review

Structural brain white matter (WM) changes such as axonal caliber, density, myelination, and orientation, along with WM-dependent structural connectivity, may be impacted early in Parkinson disease (PD). Diffusion magnetic resonance imaging (dMRI) has been used extensively to understand such pathological WM changes, and the focus of this systematic review is to understand both the methods utilized and their corresponding results in the context of early-stage PD. Diffusion tensor imaging (DTI) is the most commonly utilized method to probe WM pathological changes. Previous studies have suggested that DTI metrics are sensitive in capturing early disease-associated WM changes in preclinical symptomatic regions such as olfactory regions and the substantia nigra, which is considered to be a hallmark of PD pathology and progression. Postprocessing analytic approaches include region of interest-based analysis, voxel-based analysis, skeletonized approaches, and connectome analysis, each with unique advantages and challenges. While DTI has been used extensively to study WM disorganization in early-stage PD, it has several limitations, including an inability to resolve multiple fiber orientations within each voxel and sensitivity to partial volume effects. Given the subtle changes associated with early-stage PD, these limitations result in inaccuracies that severely impact the reliability of DTI-based metrics as potential biomarkers. To overcome these limitations, advanced dMRI acquisition and analysis methods have been employed, including diffusion kurtosis imaging and q-space diffeomorphic reconstruction. The combination of improved acquisition and analysis in DTI may yield novel and accurate information related to WM-associated changes in early-stage PD. In the current article, we present a systematic and critical review of dMRI studies in early-stage PD, with a focus on recent advances in DTI methodology. Yielding novel metrics, these advanced methods have been shown to detect diffuse WM changes in early-stage PD. These findings support the notion of early axonal damage in PD and suggest that WM pathology may go unrecognized until symptoms appear. Finally, the advantages and disadvantages of different dMRI techniques, analysis methods, and software employed are discussed in the context of PD-related pathology.

Determination of optimized set of b-values for Apparent Diffusion Coefficient mapping in liver Diffusion-Weighted MRI

In this manuscript we derive the Cramér-Rao Lower Bound (CRLB) of the monoexponential diffusion-weighted signal model under a realistic noise assumption, and propose a formulation to obtain optimized sets of b-values that maximize the noise performance of the Apparent Diffusion Coefficient (ADC) maps given a target ADC and a signal-to-noise ratio. Therefore, for various sets of parameters (S₀ and ADC), signal-to-noise ratios (SNR) and noise distribution, we computed optimized sets of b-values using CRLB-based analysis in two different ways: (i) through a greedy algorithm where b-values from a pool of candidates were added iteratively to the set, and (ii) through a two b-value search algorithm were all two b-value combinations of the pool of candidates were tested. Further, optimized sets of b-values were computed from synthetic data, phantoms, and in-vivo liver diffusion-weighted imaging (DWI) experiments to validate the CRLB-based analysis. The optimized sets of b-values obtained through the proposed CRLB-based analysis showed good agreement with the optimized sets obtained experimentally from synthetic, phantoms, and in-vivo liver data. The variance of the ADC maps decreased when employing the optimized set of b-values compared to various sets of b-values proposed in the literature for in-vivo liver DWI, although differences of notable magnitude between noise models and optimization strategies were not found. In addition, the higher b-values decreased for lower SNR under the Rician noise distribution. Optimization of the set b-values is critical to maximize the noise performance (i.e., maximize the precision and minimize the variance) of the estimated ADC maps in diffusion-weighted MRI. Hence, the proposed approach may help to optimize and standardize liver diffusion-weighted MRI acquisitions by computing optimized set of b-values for a given set of parameters.

Clinical MR Biomarkers

Non-invasive magnetic resonance imaging (MRI) techniques are increasingly applied in the clinic with a fast growing body of evidence regarding its value for clinical decision making. In contrast to biochemical or histological markers, the key advantages of imaging biomarkers are the non-invasive nature and the spatial and temporal resolution of these approaches. The following chapter focuses on clinical applications of novel MR biomarkers in humans with a strong focus on oncologic diseases. These include both clinically established biomarkers (part 1-4) and novel MRI techniques that recently demonstrated high potential for clinical utility (part 5-7).

A quick test of cognitive speed can predict development of dementia in Parkinson's disease

Parkinson’s disease (PD) patients frequently develop cognitive impairment. There is a need for brief clinical assessments identifying PD patients at high risk of progressing to dementia. In this study, we look into predicting dementia in PD and underlying structural and functional correlates to cognitive decline in PD. We included 175 patients with PD, 30 with PD dementia, 51 neurologically healthy controls and 121 patients with Alzheimer’s disease (AD) from Skane University Hospital, BIOFINDER cohorts. All underwent cognitive tests, including MMSE, 10-word list delayed recall (ADAS-cog), A Quick Test of cognitive speed (AQT), Letter S fluency, Clock Drawing Test (CDT) and pentagon copying. In non-demented patients with PD, abnormal AQT and CDT results predicted an increased risk of subsequent development of dementia (hazard ratio 2.2 for both). When comparing the cognitive profile between PD and AD, decreased performance on AQT, which measures attention and processing speed, was more typical in PD. Lastly, we investigated the underlying structural and functional correlates for the PD-specific test AQT with magnetic resonance imaging. In PD patients, decreased performance on AQT was associated with i) cortical thinning in temporoparietal regions, ii) changes in diffusion MRI, especially in the cingulum tract, and iii) decreased functional connectivity in posterior brain networks.

Diffusion Kurtosis Imaging as a Tool in Neurotoxicology

This commentary serves as an introduction to the magnetic resonance imaging (MRI) technique called diffusion kurtosis imaging (DKI) employed in the study by Arab et al. in the present issue of Neurotoxicology Research. In their study, DKI is employed for longitudinal investigation of a methamphetamine intoxication model of Parkinson's disease. The study employs an impressive number of animals and combines DKI with behavioral analysis at multiple time points. The commentary discusses some aspects of the study design especially the strength of combining behavioral analysis with MRI in an effort to provide as thorough a characterization and validity assessment of the animal model and cohort as possible. The potential clinical value of combining multiple MRI techniques (multimodal MRI) in PD is discussed as well as the benefit of multimodal MRI combined with behavioral analysis and subsequent histological analysis for in-depth characterization of animal models.

Brain Connectivity and Cognitive Flexibility in Nonirradiated Adult Survivors of Childhood Leukemia

Background

This study aimed to assess functional and structural brain connectivity in adult childhood leukemia survivors and the link with cognitive functioning and previously identified risk factors such as intrathecal methotrexate dose and age at start of therapy.

Methods

Thirty-one nonirradiated adult childhood leukemia survivors and 35 controls underwent cognitive testing and multimodal magnetic resonance imaging (resting state functional MRI, T1-weighted, diffusion-weighted, and myelin water imaging [MWI]). Analyses included dual regression, voxel-based morphometry, advanced diffusion, and MWI modeling techniques besides stepwise discriminant function analysis to identify the most affected executive cognitive domain. Correlations with discrete intrathecal MTX doses and (semi)continuous variables were calculated using Spearman's rank and Pearson's correlation, respectively. All correlation tests were two-sided. Positive and negative T-contrasts in functional and structural MRI analysis were one-sided.

Results

Survivors demonstrated lower functional connectivity between the default mode network (DMN) and inferior temporal gyrus (ITG; P < .008). Additionally, we observed higher fractional anisotropy (FA; P = .04) and lower orientation dispersion index (ODI; P = .008) at the left centrum semiovale, which could-given that several fiber bundles cross this region-suggest selective reduced integrity of the respective white matter tracts. Set shifting reaction time, a measure of cognitive flexibility, was mostly impaired and correlated with lower FA (r = -0.53, P = .003) and higher ODI (r = 0.40, P = .04) in survivors but not with DMN-ITG connectivity. There were no statistically significant differences between survivors and controls in WM or GM volume, nor was there a statistically significant correlation between imaging measurements and age at start of therapy or intrathecal methotrexate dose.

Conclusions

Adult, nonirradiated childhood leukemia survivors show altered brain connectivity, which is linked with cognitive flexibility.

Optimization of q-space sampling for mean apparent propagator MRI metrics using a genetic algorithm

Mean Apparent Propagator (MAP) MRI is a recently introduced technique to estimate the diffusion probability density function (PDF) robustly. Using the estimated PDF, MAP MRI then calculates zero-displacement and non-Gaussianity metrics, which might better characterize tissue microstructure compared to diffusion tensor imaging or diffusion kurtosis imaging. However, intensive q-space sampling required for MAP MRI limits its widespread adoption. A reduced q-space sampling scheme that maintains the accuracy of the derived metrics would make it more practical. A heuristic approach for acquiring MAP MRI with fewer q-space samples has been introduced earlier with scan duration of less than 10 minutes. However, the sampling scheme was not optimized systematically to preserve the accuracy of the model metrics. In this work, a genetic algorithm is implemented to determine optimal q-space subsampling schemes for MAP MRI that will keep total scan time under 10 min. Results show that the metrics derived from the optimized schemes more closely match those computed from the full set, especially in dense fiber tracts such as the corpus callosum.

Feasibility of Non-Gaussian Diffusion Metrics in Chronic Disorders of Consciousness

Diagnostic accuracy of different chronic disorders of consciousness (DOC) can be affected by the false negative errors in up to 40% cases. In the present study, we aimed to investigate the feasibility of a non-Gaussian diffusion approach in chronic DOC and to estimate a sensitivity of diffusion kurtosis imaging (DKI) metrics for the differentiation of vegetative state/unresponsive wakefulness syndrome (VS/UWS) and minimally conscious state (MCS) from a healthy brain state. We acquired diffusion MRI data from 18 patients in chronic DOC (11 VS/UWS, 7 MCS) and 14 healthy controls. A quantitative comparison of the diffusion metrics for grey (GM) and white (WM) matter between the controls and patient group showed a significant (p < 0.05) difference in supratentorial WM and GM for all evaluated diffusion metrics, as well as for brainstem, corpus callosum, and thalamus. An intra-subject VS/UWS and MCS group comparison showed only kurtosis metrics and fractional anisotropy differences using tract-based spatial statistics, owing mainly to macrostructural differences on most severely lesioned hemispheres. As a result, we demonstrated an ability of DKI metrics to localise and detect changes in both WM and GM and showed their capability in order to distinguish patients with a different level of consciousness.