Free breathing three-dimensional cardiac quantitative susceptibility mapping for differential cardiac chamber blood oxygenation – initial validation in patients with cardiovascular disease inclusive of direct comparison to invasive catheterization

Background Differential blood oxygenation between left (LV) and right ventricles (RV; ΔSaO2) is a key index of cardiac performance; LV dysfunction yields increased RV blood pool deoxygenation. Deoxyhemoglobin increases blood magnetic susceptibility, which can be measured using an emerging cardiovascular magnetic resonance (CMR) technique, Quantitative Susceptibility Mapping (QSM) – a concept previously demonstrated in healthy subjects using a breath-hold 2D imaging approach (2DBHQSM). This study tested utility of a novel 3D free-breathing QSM approach (3DNAVQSM) in normative controls, and validated 3DNAVQSM for non-invasive ΔSaO2 quantification in patients undergoing invasive cardiac catheterization (cath). Methods Initial control (n = 10) testing compared 2DBHQSM (ECG-triggered 2D gradient echo acquired at end-expiration) and 3DNAVQSM (ECG-triggered navigator gated gradient echo acquired in free breathing using a phase-ordered automatic window selection algorithm to partition data based on diaphragm position). Clinical testing was subsequently performed in patients being considered for cath, including 3DNAVQSM comparison to cine-CMR quantified LV function (n = 39), and invasive-cath quantified ΔSaO2 (n = 15). QSM was acquired using 3 T scanners; analysis was blinded to comparator tests (cine-CMR, cath). Results 3DNAVQSM generated interpretable QSM in all controls; 2DBHQSM was successful in 6/10. Among controls in whom both pulse sequences were successful, RV/LV susceptibility difference (and ΔSaO2) were not significantly different between 3DNAVQSM and 2DBHQSM (252 ± 39 ppb [17.5 ± 3.1%] vs. 211 ± 29 ppb [14.7 ± 2.0%]; p = 0.39). Acquisition times were 30% lower with 3DNAVQSM (4.7 ± 0.9 vs. 6.7 ± 0.5 min, p = 0.002), paralleling a trend towards lower LV mis-registration on 3DNAVQSM (p = 0.14). Among cardiac patients (63 ± 10y, 56% CAD) 3DNAVQSM was successful in 87% (34/39) and yielded higher ΔSaO2 (24.9 ± 6.1%) than in controls (p < 0.001). QSM-calculated ΔSaO2 was higher among patients with LV dysfunction as measured on cine-CMR based on left ventricular ejection fraction (29.4 ± 5.9% vs. 20.9 ± 5.7%, p < 0.001) or stroke volume (27.9 ± 7.5% vs. 22.4 ± 5.5%, p = 0.013). Cath measurements (n = 15) obtained within a mean interval of 4 ± 3 days from CMR demonstrated 3DNAVQSM to yield high correlation (r = 0.87, p < 0.001), small bias (− 0.1%), and good limits of agreement (±8.6%) with invasively measured ΔSaO2. Conclusion 3DNAVQSM provides a novel means of assessing cardiac performance. Differential susceptibility between the LV and RV is increased in patients with cine-CMR evidence of LV systolic dysfunction; QSM-quantified ΔSaO2 yields high correlation and good agreement with the reference of invasively-quantified ΔSaO2.


Introduction
Differential blood oxygenation between the left and right heart (ΔSaO 2 ) is an established index of cardiac performance; left ventricular (LV) dysfunction results in stagnant blood flowresulting in increased time for organ extraction of oxygen from blood and delivery of a greater fraction of deoxygenated blood to the right heart. Increased ΔSaO 2 has been shown to predict adverse prognosis in patients with heart failure with and without pulmonary hypertension [1][2][3] for whom it is commonly used to guide management [4,5]. However, in current clinical practice, oxygen saturation is measured by invasive catheterization (cath). Non-invasive imaging methods to measure oxygenation in the heart are limited, prohibiting non-invasive quantification of cardiac blood pool oxygenation as part of routine clinical evaluation [6][7][8][9][10][11][12][13]. Given the fact that invasive catheterization entails procedural risks and can be challenging in critically ill patients [14][15][16], a non-invasive imaging method to accurately measure cardiac oxygenation would be of substantial clinical utility.
Quantitative susceptibility mapping (QSM) is an emerging cardiovascular magnetic resonance (CMR) technique that enables quantification of diamagnetic and paramagnetic materials [17][18][19][20][21][22][23]. Iron is a magnetically active element contained in hemoglobin that is central to oxygen transport -it is weakly diamagnetic when bound to oxygen, and paramagnetic when deoxygenated [24]. This change in magnetic susceptibility by deoxyheme [25,26], provides a metric by which QSM can measure blood oxygen saturation. Prior work has validated QSM tissue characterization, including liver and brain iron content [22,[27][28][29]. Regarding blood oxygenation, a pilot study by our group showed cardiac QSM to be feasible in healthy subjects [30]. However, a 2D acquisition strategy was employed, which is suboptimal for imaging patients in whom breath-holding is often compromised. To address this, a free-breathing 3D QSM approach was developed that uses diaphragmatic navigator gating to track respiratory position. This study compared 3D free-breathing QSM (3D NAV QSM) to 2D breath held QSM (2D BH QSM) in controls, as well as to the reference of ΔSaO 2 measured in patients undergoing invasive cardiac catheterization.

Study population
3D NAV QSM was first compared to 2D BH QSM among healthy subjects without clinically reported cardiovascular conditions or associated risk factors to test the relative performance in a cohort able to undergo prolonged imaging inclusive of both pulse sequences. Next, a second group of healthy subjects were scanned and rescanned to test the reproducibility of 3D NAV QSM ΔSaO 2 measurement. The second group of healthy subjects was asked to get off and then get back on the table between the two scans. All healthy subjects were without self-reported cardiovascular disease or atherosclerosis risk factors.
After initial testing, 3D NAV QSM was then performed among clinical patients who were being considered for or had undergone invasive catheterization to quantify blood oxygen saturation. Catheterization was performed by experienced physicians using standard techniques; intracardiac blood samples were obtained under baseline conditions (without supplemental O 2 ) and used to calculate ΔSaO 2 between the left and right heart. To test the effect of gadolinium in ΔSaO 2 measurement, QSM were obtained in 3 healthy subjects and in 4 patients both pre-contrast and~30 min post-contrast administration.
This study was performed at Weill Cornell Medicine (WCM; New York, New York, USA). All participants (controls and patients) provided written informed consent for research participation. This protocol was performed with the approval from the WCM Institutional Review Board.

Data acquisition
CMR was performed using commercial 3 T scanners (750/SIGNA, General Electric Healthcare, Waukesha Wisconsin, USA). 3D NAV QSM and 2D BH QSM imaging parameters were identical between healthy subjects and patients: 1 st TE/ΔTE/#TE/TR/BW = 2.3 ms/3.6 ms/5/20 ms/111.1 kHz, acquisition matrix = 192 × 144, slice thickness = 5 mm, views per heartbeat = 10, parallel imaging factor = 2. Full 3D flow compensation was implemented for both 2D BH QSM and 3D NAV QSM to minimize the phase generated by intra-chamber blood flow [26]. To shorten 3D NAV QSM scan time, 75% partial Fourier acquisition was applied in the phase and slice encoding direction. Typical resolution is 1.5 × 1.5x5mm 3 , 40cm 2 FOV, and 20 slices per scan. Furthermore, 2D BH QSM used electrocardiographic (ECG) gating and breath-holding, and 3D NAV QSM used ECG and respiratory gating to ensure the acquisition of data at a consistent cardiac and respiratory phase. QSM was performed using non-contrast CMR in healthy subjects, and at the end of clinical exams in patients, which was approximately~30 min postgadolinium (Dotarem [gadoterate meglumine]; 0.2 mmol/ kg) infusion. As stated above, in 3 healthy subjects and 4 patients, QSM was acquired both pre-and post-contrast.
2D BH QSM employed a conventional ECG-triggered multi-echo gradient echo sequence, for which data was acquired during end-expiration (~12 s per breathhold). 3D NAV QSM employed a tailored ECG-triggered navigator gated multi-echo gradient echo sequence, for which data was acquired during free breathing: a cross pair diaphragmatic navigator was used to track respiratory motion [31,32]. A 2-bin phase-ordered automatic window selection (PAWS) gating algorithm (4 mm effective gating window) was used to tailor data acquisition according to diaphragm position in realtime [22]. In PAWS, each diaphragm position falls within a 2 mm bin for which k-space is acquired from alternating directions from bin to bin. The scan is complete when two adjacent bins have acquired all necessary k-space lines. Two navigator echoes were used in each heartbeat. The first navigator was acquired immediately before acquisition and used for PAWS gating. The second navigator was acquired immediately after data acquisition and used to provide additional motion suppression: if the difference between diaphragm positions detected by the two navigator echoes was > 4 mm for a given heartbeat, then the data acquired in that heartbeat were discarded and scheduled to be reacquired in a later heartbeat.
Ancillary imaging was performed to test QSM in relation to conventional cardiac functional/ remodeling indices. Cine-CMR was performed using a conventional balanced steady state free precession (bSSFP) pulse sequence with typical parameters: TE = 1.4 ms, TR = 3.8 ms, FA = 60 o , bandwidth = 781 Hz/Px, resolution: 1.5 × 1.5x6mm 3 , and SENSE acceleration R = 2, which was acquired in contiguous long and short axis images, the latter of which were segmented to quantify LV end-diastolic and end-systolic chamber volumes for calculation of LV and right ventricular (RV) ejection fraction (EF) as well as stroke volume.
QSM post processing QSM maps were reconstructed by first obtaining a total field map (containing both the local and background field) with the contributions of fat chemical shift removed. This was done using both graph-cut based phase unwrapping [33] and IDEAL water/fat separation [34] with iterative chemical shift update [35]. Next, a susceptibility map was obtained using the preconditioned total field inversion method [36]. In this work, two regularization terms, similar to the regularization terms described in MEDI+ 0 [37], were added to the inversion to restrict the susceptibility variations within the RV and LV as follows: The first two terms are the data fidelity term and structure consistency regularization term, respectively, where w is the signal to noise (SNR) weighting, f is the total field, d is the dipole kernel, * is the convolution operator, P is the preconditioner, λ is the regularization parameter, M G is a binary edge mask constructed by retaining the highest 70% of gradients of the T2*w image obtained by taking the square root of the sum of squares gradient recalled echo (GRE) images across echoes), and ∇ is the gradient operator [20], P is a binary mask that is 1 inside the region of interest (ROI) and a larger value, P outside , outside of the ROI (see below). The final QSM map, χ, is then χ = Py * . The last two terms constrain the susceptibility variation within the RV and the LV blood pools, where λ RV and λ LV are the regularization parameters, M RV and M LV are the mask for RV and LV obtained through manual segmentation on the GRE images. y RV and y LV are the average susceptibility over the RV and LV blood pools, respectively. In this study, the values of P outside = 20, λ = 1/1000, and λ RV = λ LV = 1/20 were empirically determined in an initial study in healthy subjects via visual inspection of the corresponding QSM, and then fixed for subsequent subjects.
To account for potential field errors from water/fat separation, an iterative reweighted least squares fitting method, MERIT [38], was implemented to modify noise weighting of fat voxels (fat fraction> 30%) in each Gauss-Newton iteration to account for residual fat chemical shift not removed from the total field in the water/fat separation step. After iteration i in the Gauss-Newton solver, the noise weighting, w, in fat voxels for the next iteration i + 1 was recalculated as the voxel-by-voxel data term residual for the i th iteration, and σ i is the standard deviation of ρ i over all voxels.
The differential susceptibility between RV and LV blood pools (Δχ) was converted to blood oxygenation difference (ΔSaO 2 ) using an established formula [30]: Where χ deoxyheme is the molar susceptibility of deoxyheme such that 4χ deoxyheme ¼ 151:054 ppb ml μmol is the molar susceptibility of a fully deoxygenated deoxyhemoglobin [39]. H ¼ 4Hct M Hb is the heme concentration in blood (in μmol/ml), where Hct is the hematocrit, ρ RBC;Hb ¼ 0:34 g ml is the mass concentration of hemoglobin in a red blood cell, and M Hb ¼ 64450 Â 10 −6 g μmol is the molar mass of deoxyhemoglobin. For controls, Hct was assumed to be 47% in men and 42% in women. These values were obtained by taking the average of the range in men and in women observed in a prior study [40]. For patients, Hct data was obtained from peripheral blood samples. Note that the current approach only measures the susceptibility difference between RV and LV blood pools (therefore only measures the oxygen saturation difference between the Fig. 1 Representative Examples of Cardiac quantitative susceptibility mapping (QSM) in Healthy Subjects. a Unsuccessful 2D BH QSM due to slice mis-registration (white arrows) attributable to inconsistent breath-hold positions, resulting in non-diagnostic QSM map. Corresponding 3D NAV QSM was successful, yielding physiologic differential oxygen saturation between the left ventricle (LV) and right ventricle (RV). b Successful 2D BH QSM and 3D NAV QSM, resulting in equivalent QSM maps RV and LV blood pools), bypassing the need to reference the blood susceptibility to a susceptibility reference (typically chosen to be water), which may or may not replicate in-vivo conditions. Post-processing was performed using MATLAB (MathWorks, Natick, Massachusetts, USA).

QSM performance scores
Image quality (of GRE images) was assessed using two different approaches. Image quality for individual short axis slices (in-plane data) was scored semi-quantitatively based on visually assessed motion artifact/endocardial blurring (0 = severe, 1 = moderate, 2 = negligible). Scoring was performed by consensus of three experienced physicians (JWW, JK, JK). Through plane image quality was measured quantitatively based on the magnitude of slice misregistration (from both in-plane and throughplane motion), which was measured as the standard deviation of the second derivative along heart surface curves that were obtained from two perpendicular reformatted long-axis images that depict the heart-lung interface [30]: a higher standard deviation corresponded to lower through plane image quality (less smooth heartlung interface).

Statistical methods
Continuous variables were compared between groups using Student's t-tests (expressed as mean ± standard deviation). 2D and 3D image quality scores were compared using a two-tailed Wilcoxon paired-sample signed rank test. Pearson correlation coefficients, Deming linear regression [41], and the Bland Altman plots were used to test the reproducibility in the two scans from the second group of healthy subjects, to test the associations between the QSM based ΔSaO 2 and the invasively quantified ΔSaO 2 , and to test the reproducibility in scans with and without contrast. Two-sided p < 0.05 was deemed indicative of statistical significance. Statistical analyses were performed using MATLAB and Prism 7 (GraphPad Software, La Jolla, California, USA).
3D NAV QSM successfully generated interpretable QSM in all 10 healthy subjects, whereas 2D BH QSM was successful in only 6/10 (60%) of cases. Figure 1 provides a representative example of as an interpretable QSM dataset acquired by 3D NAV QSM despite non-interpretable 2D BH QSM (1A), In all cases for which 2D BH QSM was uninterpretable, failure was due to slice misregistration between sequential LV short axis datasets. Consistent with this, quantitative slice misregistration was over 2-fold higher in cases for which 2D BH QSM failed (n = 4) compared to cases (n = 6) in which 2D BH QSM yielded diagnostic results (2.3 ± 0.4 vs. 1.0 ± 0.2, p < 0.001). There was no significant difference in slice misregistration between 3D NAV QSM and 2D BH QSM in exams for which both sequences were successful (1.0 ± 0.1 vs 1.0 ± 0.2, p = 1.0).
Regarding data acquisition time, results demonstrated 3D NAV QSM to yield a 30% reduction compared to 2D BH QSM (4.7 ± 0.9 vs. 6.7 ± 0.5 min, p = 0.002) attributable to the interval time between each breath-hold. Navigator efficiency for 3D NAV QSM was 54 ± 12%. Reduced scan time yielded by 3D NAV QSM remained significant even among controls in whom both pulse sequences produced diagnostic results and acquisitions were reduced by an average of 37% (4.3 ± 1.4 vs. 6.8 ± 0.4 min, p = 0.002), with navigator efficiency for 3D NAV QSM 56 ± 16%.
Reproducibility of 3D NAV QSM (as tested in 5 healthy subjects) was high, as evidenced by small mean differences (− 0.4%) and reasonable limits of agreement (± 2.2%) between data acquired during two separate scans (Fig. 2).

Clinical patients
3D NAV QSM was acquired in 39 patients whose population characteristics are shown in Table 1. In this group, QSM data was successfully obtained in 87% (34/39) of cases. In 5 cases, 3D NAV QSM yielded non-diagnostic results (nonphysiological equivalence between LV and RV blood oxygenation)all of which had substantial motion artifact. In the remainder of patients (n = 34), RV/LV susceptibility difference calculated using 3D NAV QSM was substantial (298 ± 72 ppb), corresponding to a ΔSaO 2 of 24.9 ± 6.1% (p < 0.001 vs. controls). Image acquisition time in the overall clinical cohort was 6.9 ± 2.8 min; increased acquisition time tended to be longer among patients compared to controls (4.7 ± 0.9 min; p = 0.04) due to greater respiratory variability and lower navigator efficiency in patients (36 ± 12% vs. 54 ± 12%).
In a subgroup of 15 patients who had successful QSM, invasive cardiac catheterization (cath) was available as a reference standard for heart chamber oxygenation: all patients underwent catheterization for evaluation of known/ suspected heart failure -47% had LV systolic dysfunction and 53% had primary pulmonary hypertension. Mean interval between tests (cath, CMR) was 4 ± 3 days (range 0-12 days). Figure 3b provides representative patient examples, including close agreement with invasively quantified ΔSaO 2 and increased magnitude of difference in context of LV dysfunction. As shown in Fig. 4a, QSM yielded good correlation with invasively quantified ΔSaO 2 (r = 0.87, p < 0.001); corresponding to small bias (− 0.1%) and reasonable limits of agreement (±8.6%) between the two tests (Fig. 4b). Cardiac QSM before and after contrast administration 3D NAV QSM was obtained successfully in all 7 subjects who were scanned both before and after contrast administration. As shown in Fig. 5, the ΔSaO 2 measured from pre-and post-contrast QSM matched very well (slope =

Discussion
This is the first study to test the free breathing 3D cardiac QSM for non-invasive measurement of blood oxygenation, inclusive of validation data provided by cine-CMR quantified cardiac remodeling and invasively quantified oxygen saturation in cardiac patients.  While our data validate cardiac QSM for differential chamber oxygenation, it is important to recognize that prior research has applied different CMR approaches for this purpose. Most previous approaches are based on measurement of blood CMR relaxation times (T2, T2*, and T1) [6-10, 12, 13, 42]. However, conventional methods based on longitudinal (T1) or transverse (T2) relaxation properties can be challenging to apply clinically. For example, the dependence of spin echo T2 on oxygenation is well understood, but in practice requires measuring several model parameters in addition to oxygenation, thereby potentially limiting accuracy or complicating clinical implementation. Recently, an oxygen saturation measurement based on acquiring multiple T2 maps using a 2D T2 prepared bSSFP sequence designed to overcome these limitations was shown to provide good agreement with invasive catheterization based measurement in an animal study. A comparison between this promising approach and our proposed QSM approach is warranted in a future study [11]. An alternative approach consists of quantifying the magnetic susceptibility of blood. The physical model relating blood susceptibility to oxygen saturation is simpler than that for T2 as it is linear with the slope a known physical constant. The magnetic susceptibility of venous blood can be computed from the CMR image phase by geometric modeling [43][44][45][46][47]. Mapping of magnetic susceptibility throughout the 3D field of view, as is done in QSM, enables measurement of oxygen saturation of any vascular structure (including the heart) by simple ROI analysis [25,26,48]. Our current data extends on prior work by our group that has shown QSM to provide an index of hemorrhage [26,28], as well as an index of metabolism and oxygen utilization in the brain [49,50].
Our QSM results regarding differential LV and RV blood oxygen saturation are consistent with values reported in prior literature as well as expected differences between subjects with and without cardiovascular disease. Regarding control data, ΔSaO 2 measured from 3D NAV QSM (17.5 ± 3.1%) was in agreement with a prior study that reported ΔSaO 2 in healthy subjects undergoing invasive cardiac catheterization [51], in which a mean difference of 18.8% was reported (arterial: 97.3%, venous: 78.5%). ΔSaO 2 as measured in controls were also lower than that in patients with cardiovascular disease (17.5 ± 3.1% vs. 24.9 ± 6.1%, p < 0.001), consistent with expected physiological differences between the two groups. Among our subgroup of patients with cath validation (n = 15), QSM derived ΔSaO 2 demonstrated a linear relationship with invasive measurements, and bias between the two oxygenation measurement approaches was small. While our data cannot be interpreted in context of an exact partition value with which to guide clinical decision making, it should be noted that an array of studies have demonstrated prognosis to vary in relation to magnitude of differential oxygen saturations between the left and right heart [1][2][3][4][5]. In this context, our findings suggest promise for development of QSM derived ΔSaO 2 as a non-invasive risk stratification tool, such that clinically stable patients with normal values would be screened out as low risk, whereas those with elevated QSM derived ΔSaO 2 are referred for confirmatory invasive testing. Note that, in this work, QSM was performed during non-contrast imaging in healthy subjects, and at the end of contrast-enhanced exams (~30 min post gadolinium infusion) in patients. Whereas gadolinium changes blood susceptibility, prior data has shown gadolinium to be near completely cleared from myocardium, and to be well mixed in the intravascular space at our imaging time point [52,53]. Accordingly, QSM calculations were based on the premise that contrast blood pool concentration was constant across cardiac chambers, such that gadolinium contributions to susceptibility in the LV and RV cancel out when computing differential oxygenation. Indeed, our data demonstrate that cardiac QSM derived ΔSaO 2 measurements agree well between pre-and post-contrast.
One key technical innovation in our current study concerns use of navigator technology for free breathing 3D QSM. The conventional cardiac 2D BH QSM approach generates data by imaging LV (short axis) slices individually via sequential breath-holds in order to reconstruct a 3D field map. When one or more of these breath-holds is acquired at a different respiratory position than the others, the resulting field map will not be a true 3D volume, and the QSM map will contain artifacts as the reconstruction model assumes a continuous 3D dataset as input. The susceptibility inversion problem is inherently 3D, given that susceptibility changes within one region of a given structure (e.g. RV blood pool) affects the field in the surrounding areas in all three (x, y, z) spatial directions: this is because the field is a 3D convolution of the underlying 3D susceptibility distribution with the dipole kernel. The 3D NAV QSM approach has no inherent mis-registration limitation, and the success of 3D NAV QSM acquisition depends on navigator accuracy in motion tracking. These aspects of QSM reconstruction may explain our seemingly discordant finding regarding GRE image quality and diagnostic performance of the two QSM pulse sequences tested in our study: even though GRE images from the 2D BH QSM sequence were assigned higher image scores than images from the 3D NAV QSM sequence, 2D BH QSM failed to generate interpretable QSM more often than did 3D NAV QSM and this failure was primarily attributable to slice mis-registration.
There are several limitations of our study. First, scan time of 3D NAV QSM sequence is long (~4-7 min) compared to other routine clinical sequences. One potential approach to shorten scan time is to increase parallel acceleration factor, use compressed sensing [54,55], and/ or apply data acquisition strategies such as echo planar readout [56] or echo sharing methods [57]. Non-Cartesian acquisitions that allow for self-gating [58] and multi-phase reconstruction [59] are also alternatives to shorten cardiac QSM. A second limitation to the current 3D NAV QSM approach is that the respiratory motion is tracked by a diaphragmatic navigator, which is known to occasionally fail [60]. More advanced navigator techniques such as a fat navigator can be used to improve success rates [31,32,[61][62][63]. A third issue to consider is that the QSM validation in this study was derived from cardiac remodeling indices in 39 cardiac patients of whom 15 had invasive catheterization for evaluation of known or suspected heart failure; further validation in a larger clinical cohort including patients undergoing catheterization for different indications is warranted. Finally, whereas the interval between QSM and invasive catheterization was relatively short (mean 4 ± 3 days, median 2 days [IQR 1-5 days]), cardiac chamber oxygenation could have varied during this time, thus resulting in discordance between tests. Future research with simultaneous invasive and CMR oxygenation measurements in animal models could be of utility for further validation of QSM, as well as comparisons with alternative pulse sequences for oxygenation quantification (i.e. T2). Nevertheless, our current findings that 3D navigator QSM was feasible among a clinical cohort in whom it generally agreed with differential oxygenation saturation as measured invasively is of substantial importance with respect to translational application of this technique, and provides initial proof of concept with respect to clinical implementation.
Heart rate could have impacted our QSM results. This has been shown to be the case for oxygenation methods such as BOLD, which relies on the magnitude of the CMR signal. On the other hand, QSM, which relies on the phase of the CMR signal, is expected to be relatively insensitive to heart rate because the contributions to the phase (field) induced by the difference in LV and RV magnetic susceptibility (oxygenation) are not affected by heart rate. Further research is warranted to specifically assess physiologic factors impacting QSM, as well as to validate this pulse sequence in a larger clinical cohort including patients undergoing catheterization for different indications.
In conclusion, we provide validation of free breathing 3D cardiac QSM as an index of cine-CMR evidenced LV dysfunction and differential LV/RV oxygen saturation as measured by invasive catheterization. Future research is necessary to test accelerated free-breathing QSM strategies, refine cardiac QSM for myocardial tissue characterization, and validate QSM-derived blood oxygenation for non-invasive stratification of heart failure symptoms and prognostic outcomes.