Multi-echo quantitative susceptibility mapping: how to combine echoes for accuracy and precision at 3 Tesla

Abstract

MRI; Multi-echo QSM; Quantitative susceptibility mapping

Imágen por resonancia magnética; QSM de ecos múltiples; Mapeo cuantitativo de susceptibilidad

Imatge per ressonància magnètica; QSM de ressò múltiple; Mapeig quantitatiu de susceptibilitat

Purpose To compare different multi-echo combination methods for MRI QSM. Given the current lack of consensus, we aimed to elucidate how to optimally combine multi-echo gradient-recalled echo signal phase information, either before or after applying Laplacian-base methods (LBMs) for phase unwrapping or background field removal. Methods Multi-echo gradient-recalled echo data were simulated in a numerical head phantom, and multi-echo gradient-recalled echo images were acquired at 3 Tesla in 10 healthy volunteers. To enable image-based estimation of gradient-recalled echo signal noise, 5 volunteers were scanned twice in the same session without repositioning. Five QSM processing pipelines were designed: 1 applied nonlinear phase fitting over TEs before LBMs; 2 applied LBMs to the TE-dependent phase and then combined multiple TEs via either TE-weighted or SNR-weighted averaging; and 2 calculated TE-dependent susceptibility maps via either multi-step or single-step QSM and then combined multiple TEs via magnitude-weighted averaging. Results from different pipelines were compared using visual inspection; summary statistics of susceptibility in deep gray matter, white matter, and venous regions; phase noise maps (error propagation theory); and, in the healthy volunteers, regional fixed bias analysis (Bland–Altman) and regional differences between the means (nonparametric tests). Results Nonlinearly fitting the multi-echo phase over TEs before applying LBMs provided the highest regional accuracy of and the lowest phase noise propagation compared to averaging the LBM-processed TE-dependent phase. This result was especially pertinent in high-susceptibility venous regions. Conclusion For multi-echo QSM, we recommend combining the signal phase by nonlinear fitting before applying LBMs.

Supported by the UK Engineering and Physical Sciences Research Council (EPSRC), award number: 1489882 (e.b.); by the EPSRC-funded UCL Centre for Doctoral Training in Medical Imaging, grant EP/L016478/1 (a.k.), and the Department of Health's National Institute for Health Research funded Biomedical Research Centre at University College London Hospitals (a.k.); by the UCL Leonard Wolfson Experimental Neurology Centre, grant PR/ylr/18575 (d.l.t) The Queen Square MS Centre, where part of the MRI scans for this work were performed, is supported by grants from the UK MS Society and by the National Institute for Health Research University College London Hospitals Biomedical Research Centre (UCLH/BRC). F. Grussu was supported by PREdICT, a study at the Vall d'Hebron Institute of Oncology in Barcelona funded by AstraZeneca (f.g.), and funding from the postdoctoral fellowships program Beatriu de Pinós (2020 BP 00117), funded by the Secretary of Universities and Research, Government of Catalonia (f.g.)