The Use of Magnetic Resonance Imaging in Radiation Therapy Treatment Simulation and Planning.

Ever since its introduction as a diagnostic imaging tool the potential of magnetic resonance imaging (MRI) in radiation therapy (RT) treatment simulation and planning has been recognized. Recent technical advances have addressed many of the impediments to use of this technology and as a result have resulted in rapid and growing adoption of MRI in RT. The purpose of this article is to provide a broad review of the multiple uses of MR in the RT treatment simulation and planning process, identify several of the most used clinical scenarios in which MR is integral to the simulation and planning process, highlight existing limitations and provide multiple unmet needs thereby highlighting opportunities for the diagnostic MR imaging community to contribute and collaborate with our oncology colleagues. EVIDENCE LEVEL: 5 TECHNICAL EFFICACY: Stage 5.

Magnetic Resonance Imaging Artifact Associated With Transcutaneous Bone Conduction Implants: Cholesteatoma and Vestibular Schwannoma Surveillance.

OBJECTIVE: To evaluate the magnetic resonance (MR) image artifact and image distortion associated with the two transcutaneous bone conduction implants currently available in the United States. STUDY DESIGN: Cadaveric study. METHODS: Two cadaveric head specimens (1 male, 1 female) were unilaterally implanted according to manufacturer guidelines and underwent MR imaging (General Electric and Siemens 1.5 T scanners) under the following device conditions: (1) no device, (2) Cochlear Osia with magnet and headwrap, (3) Cochlear Osia without magnet, and (4) MED-EL Bonebridge with magnet. Maximum metal mitigation techniques were employed in all conditions, and identical sequences were obtained. Blinded image scoring (diagnostic vs nondiagnostic image) was performed by experienced neuroradiologists according to anatomical subsites. RESULTS: All device conditions produced artifact and image distortion. The Osia with magnet produced diagnostic T1- and T2-weighted images of the ipsilateral temporal bone, however, non-echo planar imaging diffusion-weighted imaging (DWI) was nondiagnostic. The Osia without magnet scanned on the Siemens MR imaging demonstrated the least amount of artifact and was the only condition that allowed for diagnostic imaging of the ipsilateral temporal bone on DWI. The Bonebridge produced a large area of artifact and distortion with the involvement of the ipsilateral and contralateral temporal bones. CONCLUSION: In summary, of the three device conditions (Osia with magnet, Osia without magnet, and Bonebridge), Osia without magnet offered the least amount of artifact and distortion and was the only condition in which diagnostic DWI was available for the middle ear and mastoid regions on the Siemens MR imaging scanner.

Evaluation of a New, Highly Flexible Radiofrequency Coil for MR Simulation of Patients Undergoing External Beam Radiation Therapy.

PURPOSE: To evaluate the performance of a new, highly flexible radiofrequency (RF) coil system for imaging patients undergoing MR simulation. METHODS: Volumetric phantom and in vivo images were acquired with a commercially available and prototype RF coil set. Phantom evaluation was performed using a silicone-filled humanoid phantom of the head and shoulders. In vivo assessment was performed in five healthy and six patient subjects. Phantom data included T1-weighted volumetric imaging, while in vivo acquisitions included both T1- and T2-weighted volumetric imaging. Signal to noise ratio (SNR) and uniformity metrics were calculated in the phantom data, while SNR values were calculated in vivo. Statistical significance was tested by means of a non-parametric analysis of variance test. RESULTS: At a threshold of p = 0.05, differences in measured SNR distributions within the entire phantom volume were statistically different in two of the three paired coil set comparisons. Differences in per slice average SNR between the two coil sets were all statistically significant, as well as differences in per slice image uniformity. For patients, SNRs within the entire imaging volume were statistically significantly different in four of the nine comparisons and seven of the nine comparisons performed on the per slice average SNR values. For healthy subjects, SNRs within the entire imaging volume were statistically significantly different in seven of the nine comparisons and eight of the nine comparisons when per slice average SNR was tested. CONCLUSIONS: Phantom and in vivo results demonstrate that image quality obtained from the novel flexible RF coil set was similar or improved over the conventional coil system. The results also demonstrate that image quality is impacted by the specific coil configurations used for imaging and should be matched appropriately to the anatomic site imaged to ensure optimal and reproducible image quality.

Left-Right Intensity Asymmetries Vary Depending on Scanner Model for FLAIR and T1 Weighted MRI Images.

BACKGROUND: Localized regions of left-right image intensity asymmetry (LRIA) were incidentally observed on T2 -weighted (T2 -w) and T1 -weighted (T1 -w) diagnostic magnetic resonance imaging (MRI) images. Suspicion of herpes encephalitis resulted in unnecessary follow-up imaging. A nonbiological imaging artifact that can lead to diagnostic uncertainty was identified. PURPOSE: To investigate whether systematic LRIA exist for a range of scanner models and to determine if LRIA can introduce diagnostic uncertainty. STUDY TYPE: A retrospective study using the Alzheimer's Disease Neuroimaging Initiative (ADNI) data base. SUBJECTS: One thousand seven hundred fifty-three (median age: 72, males/females: 878/875) unique participants with longitudinal data were included. FIELD STRENGTH: 3T. SEQUENCES: T1 -w three-dimensional inversion-recovery spoiled gradient-echo (IR-SPGR) or magnetization-prepared rapid gradient-echo (MP-RAGE) and T2 -w fluid-attenuated inversion recovery (FLAIR) long tau fast spin echo inversion recovery (LT-FSE-IR). Only General Electric, Philips, and Siemens' product sequences were used. ASSESSMENT: LRIA was calculated as the left-right percent difference with respect to the mean intensity from automated anatomical atlas segmented regions. Three neuroradiologists with 37 (**), 32 (**), and 3 (**) years of experience rated the clinical impact of 30 T2 -w three-dimensional FLAIR exams with LRIA to determine the diagnostic uncertainty. Statistical comparisons between retrospective intensity normalized T1 m and original T1 -w images were made. STATISTICAL TESTS: For each image type, a linear mixed effects model was fit using LRIA scores from all scanners, regions, and participants as the outcome and age and sex as predictors. Statistical significance was defined as having a P-value <0.05. RESULTS: LRIA scores were significantly different from zero on most scanners. All clinicians were uncertain or recommended definite diagnostic follow-up in 62.5% of cases with LRIA >10%. Individuals with acute brain pathology or focal neurologic deficits are not enrolled in ADNI; therefore, focal signal abnormalities were considered false positives. DATA CONCLUSION: LRIA is system specific, systematic, creates diagnostic uncertainty, and impacts IR-SPGR, MP-RAGE, and LT-FSE-IR product sequences. LEVEL OF EVIDENCE: 2 Technical Efficacy Stage: 3.

Findings of the AAPM Ad Hoc committee on magnetic resonance imaging in radiation therapy: Unmet needs, opportunities, and recommendations.

The past decade has seen the increasing integration of magnetic resonance (MR) imaging into radiation therapy (RT). This growth can be contributed to multiple factors, including hardware and software advances that have allowed the acquisition of high-resolution volumetric data of RT patients in their treatment position (also known as MR simulation) and the development of methods to image and quantify tissue function and response to therapy. More recently, the advent of MR-guided radiation therapy (MRgRT) - achieved through the integration of MR imaging systems and linear accelerators - has further accelerated this trend. As MR imaging in RT techniques and technologies, such as MRgRT, gain regulatory approval worldwide, these systems will begin to propagate beyond tertiary care academic medical centers and into more community-based health systems and hospitals, creating new opportunities to provide advanced treatment options to a broader patient population. Accompanying these opportunities are unique challenges related to their adaptation, adoption, and use including modification of hardware and software to meet the unique and distinct demands of MR imaging in RT, the need for standardization of imaging techniques and protocols, education of the broader RT community (particularly in regards to MR safety) as well as the need to continue and support research, and development in this space. In response to this, an ad hoc committee of the American Association of Physicists in Medicine (AAPM) was formed to identify the unmet needs, roadblocks, and opportunities within this space. The purpose of this document is to report on the major findings and recommendations identified. Importantly, the provided recommendations represent the consensus opinions of the committee's membership, which were submitted in the committee's report to the AAPM Board of Directors. In addition, AAPM ad hoc committee reports differ from AAPM task group reports in that ad hoc committee reports are neither reviewed nor ultimately approved by the committee's parent groups, including at the council and executive committee level. Thus, the recommendations given in this summary should not be construed as being endorsed by or official recommendations from the AAPM.

Whole-brain 3D MR fingerprinting brain imaging: clinical validation and feasibility to patients with meningioma.

PURPOSE: MR fingerprinting (MRF) is a MR technique that allows assessment of tissue relaxation times. The purpose of this study is to evaluate the clinical application of this technique in patients with meningioma. MATERIALS AND METHODS: A whole-brain 3D isotropic 1mm3 acquisition under a 3.0T field strength was used to obtain MRF T1 and T2-based relaxometry values in 4:38 s. The accuracy of values was quantified by scanning a quantitative MR relaxometry phantom. In vivo evaluation was performed by applying the sequence to 20 subjects with 25 meningiomas. Regions of interest included the meningioma, caudate head, centrum semiovale, contralateral white matter and thalamus. For both phantom and subjects, mean values of both T1 and T2 estimates were obtained. Statistical significance of differences in mean values between the meningioma and other brain structures was tested using a Friedman's ANOVA test. RESULTS: MR fingerprinting phantom data demonstrated a linear relationship between measured and reference relaxometry estimates for both T1 (r2 = 0.99) and T2 (r2 = 0.97). MRF T1 relaxation times were longer in meningioma (mean ± SD 1429 ± 202 ms) compared to thalamus (mean ± SD 1054 ± 58 ms; p = 0.004), centrum semiovale (mean ± SD 825 ± 42 ms; p < 0.001) and contralateral white matter (mean ± SD 799 ± 40 ms; p < 0.001). MRF T2 relaxation times were longer for meningioma (mean ± SD 69 ± 27 ms) as compared to thalamus (mean ± SD 27 ± 3 ms; p < 0.001), caudate head (mean ± SD 39 ± 5 ms; p < 0.001) and contralateral white matter (mean ± SD 35 ± 4 ms; p < 0.001) CONCLUSIONS: Phantom measurements indicate that the proposed 3D-MRF sequence relaxometry estimations are valid and reproducible. For in vivo, entire brain coverage was obtained in clinically feasible time and allows quantitative assessment of meningioma in clinical practice.

Whole brain 3D MR fingerprinting in multiple sclerosis: a pilot study.

BACKGROUND: MR fingerprinting (MRF) is a novel imaging method proposed for the diagnosis of Multiple Sclerosis (MS). This study aims to determine if MR Fingerprinting (MRF) relaxometry can differentiate frontal normal appearing white matter (F-NAWM) and splenium in patients diagnosed with MS as compared to controls and to characterize the relaxometry of demyelinating plaques relative to the time of diagnosis. METHODS: Three-dimensional (3D) MRF data were acquired on a 3.0T MRI system resulting in isotropic voxels (1 × 1 × 1 mm3) and a total acquisition time of 4 min 38 s. Data were collected on 18 subjects paired with 18 controls. Regions of interest were drawn over MRF-derived T1 relaxometry maps encompassing selected MS lesions, F-NAWM and splenium. T1 and T2 relaxometry features from those segmented areas were used to classify MS lesions from F-NAWM and splenium with T-distributed stochastic neighbor embedding algorithms. Partial least squares discriminant analysis was performed to discriminate NAWM and Splenium in MS compared with controls. RESULTS: Mean out-of-fold machine learning prediction accuracy for discriminant results between MS patients and controls for F-NAWM was 65 % (p = 0.21) and approached 90 % (p < 0.01) for the splenium. There was significant positive correlation between time since diagnosis and MS lesions mean T2 (p = 0.015), minimum T1 (p = 0.03) and negative correlation with splenium uniformity (p = 0.04). Perfect discrimination (AUC = 1) was achieved between selected features from MS lesions and F-NAWM. CONCLUSIONS: 3D-MRF has the ability to differentiate between MS and controls based on relaxometry properties from the F-NAWM and splenium. Whole brain coverage allows the assessment of quantitative properties within lesions that provide chronological assessment of the time from MS diagnosis.

Magnetic resonance biomarkers in radiation oncology: The report of AAPM Task Group 294.

PURPOSE: A magnetic resonance (MR) biologic marker (biomarker) is a measurable quantitative characteristic that is an indicator of normal biological and pathogenetic processes or a response to therapeutic intervention derived from the MR imaging process. There is significant potential for MR biomarkers to facilitate personalized approaches to cancer care through more precise disease targeting by quantifying normal versus pathologic tissue function as well as toxicity to both radiation and chemotherapy. Both of which have the potential to increase the therapeutic ratio and provide earlier, more accurate monitoring of treatment response. The ongoing integration of MR into routine clinical radiation therapy (RT) planning and the development of MR guided radiation therapy systems is providing new opportunities for MR biomarkers to personalize and improve clinical outcomes. Their appropriate use, however, must be based on knowledge of the physical origin of the biomarker signal, the relationship to the underlying biological processes, and their strengths and limitations. The purpose of this report is to provide an educational resource describing MR biomarkers, the techniques used to quantify them, their strengths and weakness within the context of their application to radiation oncology so as to ensure their appropriate use and application within this field.

TURBINE-MRE: A 3D hybrid radial-Cartesian EPI acquisition for MR elastography.

PURPOSE: To develop a novel magnetic resonance elastography (MRE) acquisition using a hybrid radial EPI readout scheme (TURBINE), and to demonstrate its feasibility to obtain wave images and stiffness maps in a phantom and in vivo brain. METHOD: The proposed 3D TURBINE-MRE is based on a spoiled gradient-echo MRE sequence with the EPI readout radially rotating about the phase-encoding axis to sample a full 3D k-space. A polyvinyl chloride phantom and 6 volunteers were scanned on a compact 3T GE scanner with a 32-channel head coil at 80 Hz and 60 Hz external vibration, respectively. For comparison, a standard 2D, multislice, spin-echo (SE) EPI-MRE acquisition was also performed with the same motion encoding and resolution. The TURBINE-MRE images were off-line reconstructed with iterative SENSE algorithm. The regional ROI analysis was performed on the 6 volunteers, and the median stiffness values were compared between SE-EPI-MRE and TURBINE-MRE. RESULTS: The 3D wave-field images and the generated stiffness maps were comparable between TURBINE-MRE and standard SE-EPI-MRE for the phantom and the volunteers. The Bland-Altman plot showed no significant difference in the median regional stiffness values between the two methods. The stiffness measured with the 2 methods had a strong linear relationship with a Pearson correlation coefficient of 0.943. CONCLUSION: We demonstrated the feasibility of the new TURBINE-MRE sequence for acquiring the desired 3D wave-field data and stiffness maps in a phantom and in-vivo brains. This pilot study encourages further exploration of TURBINE-MRE for functional MRE, free-breathing abdominal MRE, and cardiac MRE applications.

Application of Adaptive Image Receive Coil Technology for Whole-Brain Imaging.

OBJECTIVE. The Adaptive Image Receive (AIR) radiofrequency coil is an emergent technology that is lightweight and flexible and exhibits electrical characteristics that overcome many of the limitations of traditional rigid coil designs. The purpose of this study was to apply the AIR coil for whole-brain imaging and compare the performance of a prototype AIR coil array with the performance of conventional head coils. SUBJECTS AND METHODS. A phantom and 15 healthy adult participants were imaged. A prototype 16-channel head AIR coil was compared with conventional 8-and 32-channel head coils using clinically available MRI sequences. During consensus review, two board-certified neuroradiologists graded the AIR coil compared with an 8-channel coil and a 32-channel coil on a 5-point ordinal scale in multiple categories. One- and two-sided Wilcoxon signed rank tests were performed. Noise covariance matrices and geometry factor (g-factor) maps were calculated. RESULTS. The signal-to-noise ratio, structural sharpness, and overall image quality scores of the prototype 16-channel AIR coil were better than those of the 8-channel coil but were not as good as those of the 32-channel coil. Noise covariance matrices showed stable performance of the AIR coil across participants. The median g-factors for the 16-channel AIR coil were, overall, less than those of the 8-channel coil but were greater than those of the 32-channel coil. CONCLUSION. On average, the prototype 16-channel head AIR coil outperformed a conventional 8-channel head coil but did not perform as well as a conventional 32-channel head coil. This study shows the feasibility of the novel AIR coil technology for imaging the brain and provides insight for future coil design improvements.

Artificial neural networks for magnetic resonance elastography stiffness estimation in inhomogeneous materials.

PURPOSE: To test the hypothesis that removing the assumption of material homogeneity will improve the spatial accuracy of stiffness estimates made by Magnetic Resonance Elastography (MRE). METHODS: An artificial neural network was trained using synthetic wave data computed using a coupled harmonic oscillator model. Material properties were allowed to vary in a piecewise smooth pattern. This neural network inversion (Inhomogeneous Learned Inversion (ILI)) was compared against a previous homogeneous neural network inversion (Homogeneous Learned Inversion (HLI)) and conventional direct inversion (DI) in simulation, phantom, and in-vivo experiments. RESULTS: In simulation experiments, ILI was more accurate than HLI and DI in predicting the stiffness of an inclusion in noise-free, low-noise, and high-noise data. In the phantom experiment, ILI delineated inclusions ≤ 2.25 cm in diameter more clearly than HLI and DI, and provided a higher contrast-to-noise ratio for all inclusions. In a series of stiff brain tumors, ILI shows sharper stiffness transitions at the edges of tumors than the other inversions evaluated. CONCLUSION: ILI is an artificial neural network based framework for MRE inversion that does not assume homogeneity in material stiffness. Preliminary results suggest that it provides more accurate stiffness estimates and better contrast in small inclusions and at large stiffness gradients than existing algorithms that assume local homogeneity. These results support the need for continued exploration of learning-based approaches to MRE inversion, particularly for applications where high resolution is required.

Role and future of MRI in radiation oncology.

Technical innovations and developments in areas such as disease localization, dose calculation algorithms, motion management and dose delivery technologies have revolutionized radiation therapy resulting in improved patient care with superior outcomes. A consequence of the ability to design and accurately deliver complex radiation fields is the need for improved target visualization through imaging. While CT imaging has been the standard of care for more than three decades, the superior soft tissue contrast afforded by MR has resulted in the adoption of this technology in radiation therapy. With the development of real time MR imaging techniques, the problem of real time motion management is enticing. Currently, the integration of an MR imaging and megavoltage radiation therapy treatment delivery system (MR-linac or MRL) is a reality that has the potential to provide improved target localization and real time motion management during treatment. Higher magnetic field strengths provide improved image quality potentially providing the backbone for future work related to image texture analysis-a field known as Radiomics-thereby providing meaningful information on the selection of future patients for radiation dose escalation, motion-managed treatment techniques and ultimately better patient care. On-going advances in MRL technologies promise improved real time soft tissue visualization, treatment margin reductions, beam optimization, inhomogeneity corrected dose calculation, fast multileaf collimators and volumetric arc radiation therapy. This review article provides rationale, advantages and disadvantages as well as ideas for future research in MRI related to radiation therapy mainly in adoption of MRL.

Quantifying Tumor Stiffness With Magnetic Resonance Elastography: The Role of Mechanical Properties for Detection, Characterization, and Treatment Stratification in Oncology.

The viscoelastic properties of tissue are significantly altered with the development of tumors and these alterations can be assessed with magnetic resonance elastography (MRE). Accurate detection and characterization of malignant and benign lesions can be obtained by quantifying tumor stiffness, improving the specificity and diagnostic accuracy of conventional magnetic resonance imaging. Furthermore, MRE can be used to stratify patients for treatment based on risk of normal tissue toxicity and surgical considerations including consistency and adherence of the tumor to surrounding structures. MRE is a reliable reproducible technique demonstrated in studies that include both patients with cancer and normal volunteers, and an average technical failure rate of <1%. The addition of MRE into a multiparametric magnetic resonance imaging assessment may improve the diagnosis and treatment of cancer.

Characterization and evaluation of a flexible MRI receive coil array for radiation therapy MR treatment planning using highly decoupled RF circuits.

The growth in the use of magnetic resonance imaging (MRI) data for radiation therapy (RT) treatment planning has been facilitated by scanner hardware and software advances that have enabled RT patients to be imaged in treatment position while providing morphologic and functional assessment of tumor volumes and surrounding normal tissues. Despite these advances, manufacturers have been slow to develop radiofrequency (RF) coils that closely follow the contour of a RT patient undergoing MR imaging. Instead, relatively large form surface coil arrays have been adapted from diagnostic imaging. These arrays can be challenging to place on, and in general do not conform to the patient's body habitus, resulting in sub optimal image quality. The purpose of this study is to report on the characterization of a new flexible and highly decoupled RF coil for use in MR imaging of RT patients. Coil performance was evaluated by performing signal-to-noise ratio (SNR) and noise correlation measurements using two coil (SNR) and four coil (noise correlation) element combinations as a function of coil overlap distance and comparing these values to those obtained using conventional coil elements. In vivo testing was performed in both normal volunteers and patients using a four and 16 element RF coil. Phantom experiments demonstrate the highly decoupled nature of the new coil elements when compared to conventional RF coils, while in vivo testing demonstrate that these coils can be integrated into extremely flexible and form fitting substrates that follow the exact contour of the patient. The new coil design addresses limitations imposed by traditional surface coil arrays and have the potential to significantly impact MR imaging for both diagnostic and RT applications.

Regional assessment of in vivo myocardial stiffness using 3D magnetic resonance elastography in a porcine model of myocardial infarction.

PURPOSE: The stiffness of a myocardial infarct affects the left ventricular pump function and remodeling. Magnetic resonance elastography (MRE) is a noninvasive imaging technique for measuring soft-tissue stiffness in vivo. The purpose of this study was to investigate the feasibility of assessing in vivo regional myocardial stiffness with high-frequency 3D cardiac MRE in a porcine model of myocardial infarction, and compare the results with ex vivo uniaxial tensile testing. METHODS: Myocardial infarct was induced in a porcine model by embolizing the left circumflex artery. Fourteen days postinfarction, MRE imaging was performed in diastole using an echocardiogram-gated spin-echo echo-planar-imaging sequence with 140-Hz vibrations and 3D MRE processing. The MRE stiffness and tensile modulus from uniaxial testing were compared between the remote and infarcted myocardium. RESULTS: Myocardial infarcts showed increased in vivo MRE stiffness compared with remote myocardium (4.6 ± 0.7 kPa versus 3.0 ± 0.6 kPa, P = 0.02) within the same pig. Ex vivo uniaxial mechanical testing confirmed the in vivo MRE results, showing that myocardial infarcts were stiffer than remote myocardium (650 ± 80 kPa versus 110 ± 20 kPa, P = 0.01). CONCLUSIONS: These results demonstrate the feasibility of assessing in vivo regional myocardial stiffness with high-frequency 3D cardiac MRE. Magn Reson Med 79:361-369, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Cardiac MR elastography using reduced-FOV, single-shot, spin-echo EPI.

PURPOSE: To implement a reduced field of view (rFOV) technique for cardiac MR elastography (MRE) and to demonstrate the improvement in image quality of both magnitude images and post-processed MRE stiffness maps compared to the conventional full field of view (full-FOV) acquisition. METHODS: With Institutional Review Board approval, 17 healthy volunteers underwent both full-FOV and rFOV cardiac MRE scans using 140-Hz vibrations. Two cardiac radiologists blindly compared the magnitude images and stiffness maps and graded the images based on several image quality attributes using a 5-point ordinal scale. Fisher's combined probability test was performed to assess the overall evaluation. The octahedral shear strain-based signal-to-noise ratio (OSS-SNR) and median stiffness over the left ventricular myocardium were also compared. RESULTS: One volunteer was excluded because of an inconsistent imaging resolution during the exam. In the remaining 16 volunteers (9 males, 7 females), the rFOV scans outperformed the full-FOV scans in terms of subjective image quality and ghosting artifacts in the magnitude images and stiffness maps, as well as the overall preference. The quantitative measurements showed that rFOV had significantly higher OSS-SNR (median: 1.4 [95% confidence interval (CI): 1.2-1.5] vs. 2.1 [95% CI: 1.8-2.4]), P < 0.05) compared to full-FOV. Although no significant change was found in the median myocardial stiffness between the 2 scans, we observed a decrease in the stiffness variation within the myocardium from 2.1 kPa (95% CI: [1.9, 2.3]) to 1.9 kPa (95% CI: [1.7, 2.0]) for full-FOV and rFOV, respectively (P < 0.05) in a subgroup of 7 subjects with ghosting present in the myocardium. CONCLUSION: This pilot volunteer study demonstrated that rFOV cardiac MRE has the capability to reduce ghosting and to improve image quality in both MRE magnitude images and stiffness maps. Magn Reson Med 80:231-238, 2018. © 2017 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution NonCommercial License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited and is not used for commercial purposes.