PurposeQSM provides insight into healthy brain aging and neuropathologies such as multiple sclerosis (MS), traumatic brain injuries, brain tumors, and neurodegenerative diseases. Phase data for QSM are usually acquired from 3D gradient‐echo (3D GRE) scans with long acquisition times that are detrimental to patient comfort and susceptible to patient motion. This is particularly true for scans requiring whole‐brain coverage and submillimeter resolutions. In this work, we use a multishot 3D echo plannar imaging (3D EPI) sequence with shot‐selective 2D CAIPIRIHANA to acquire high‐resolution, whole‐brain data for QSM with minimal distortion and blurring.MethodsTo test clinical viability, the 3D EPI sequence was used to image a cohort of MS patients at 1‐mm isotropic resolution at 3 T. Additionally, 3D EPI data of healthy subjects were acquired at 1‐mm, 0.78‐mm, and 0.65‐mm isotropic resolution with varying echo train lengths (ETLs) and compared with a reference 3D GRE acquisition.ResultsThe appearance of the susceptibility maps and the susceptibility values for segmented regions of interest were comparable between 3D EPI and 3D GRE acquisitions for both healthy and MS participants. Additionally, all lesions visible in the MS patients on the 3D GRE susceptibility maps were also visible on the 3D EPI susceptibility maps. The interplay among acquisition time, resolution, echo train length, and the effect of distortion on the calculated susceptibility maps was investigated.ConclusionWe demonstrate that the 3D EPI sequence is capable of rapidly acquiring submillimeter resolutions and providing high‐quality, clinically relevant susceptibility maps.