QSM reconstruction was performed by estimating the total field via a non-linear fit of the mGRE,41 (link) obtaining the local field by the Projection onto Dipole Fields (PDF) method,42 (link) and computing susceptibility with the Morphology Enabled Dipole Inversion with automatic uniform cerebrospinal fluid zero reference (MEDI+0) algorithm.43 (link)–46 (link) Cerebral blood flow (CBF) maps (mL/100g/min) were generated from the ASL data using FuncTool (GE Healthcare, Waukesha, WI, USA). All images were co-registered and interpolated to the resolution of the QSM maps using the FSL FLIRT algorithm.47 (link),48 (link)