- Open Access
Non-invasive assessment of cardiac function and pulmonary vascular resistance in an canine model of acute thromboembolic pulmonary hypertension using 4D flow cardiovascular magnetic resonance
Journal of Cardiovascular Magnetic Resonance volume 16, Article number: 23 (2014)
The purpose of this study was to quantify right (RV) and left (LV) ventricular function, pulmonary artery flow (QP), tricuspid valve regurgitation velocity (TRV), and aorta flow (QS) from a single 4D flow cardiovascular magnetic resonance (CMR) (time-resolved three-directionally motion encoded CMR) sequence in a canine model of acute thromboembolic pulmonary hypertension (PH).
Acute PH was induced in six female beagles by microbead injection into the right atrium. Pulmonary arterial (PAP) and pulmonary capillary wedge (PCWP) pressures and cardiac output (CO) were measured by right heart catheterization (RHC) at baseline and following induction of acute PH. Pulmonary vascular resistance (PVRRHC) was calculated from RHC values of PAP, PCWP and CO (PVRRHC = (PAP-PCWP)/CO). Cardiac magnetic resonance (CMR) was performed on a 3 T scanner at baseline and following induction of acute PH. RV and LV end-diastolic (EDV) and end-systolic (ESV) volumes were determined from both CINE balanced steady-state free precession (bSSFP) and 4D flow CMR magnitude images. QP, TRV, and QS were determined from manually placed cutplanes in the 4D flow CMR flow-sensitive images in the main (MPA), right (RPA), and left (LPA) pulmonary arteries, the tricuspid valve (TRV), and aorta respectively. MPA, RPA, and LPA flow was also measured using two-dimensional flow-sensitive (2D flow) CMR.
Biases between 4D flow CMR and bSSFP were 0.8 mL and 1.6 mL for RV EDV and RV ESV, respectively, and 0.8 mL and 4 mL for LV EDV and LV ESV, respectively. Flow in the MPA, RPA, and LPA did not change after induction of acute PAH (p = 0.42-0.81). MPA, RPA, and LPA flow determined with 4D flow CMR was significantly lower than with 2D flow (p < 0.05). The correlation between QP/TRV and PVRRHC was 0.95. The average QP/QS was 0.96 ± 0.11.
Using both magnitude and flow-sensitive data from a single 4D flow CMR acquisition permits simultaneous quantification of cardiac function and cardiopulmonary hemodynamic parameters important in the assessment of PH.
Pulmonary hypertension (PH) is an ultimately fatal disease characterized by an abnormal increase in the mean pulmonary artery pressure (mPAP) that ultimately leads to right ventricular (RV) failure. The initial evaluation of patients with suspected PH includes non-invasive imaging with echocardiography because of its ability to non-invasively estimate pulmonary artery pressures, pulmonary vascular resistance (PVR) and cardiac function. Depending upon the results of the initial assessment with echocardiography, further imaging work-up to identify the underlying cause and determine appropriate management of PH may include ventilation/perfusion scintigraphy in patients with known or suspected chronic thromboembolic PH, computed tomography to assess for diffuse lung disease, and invasive right heart catheterization (RHC) to confirm the diagnosis. Currently, RHC is considered the reference standard for assessing PH patients and providing definite diagnosis, because the categorization of the type of PH is based on the observed pulmonary artery pressures (mPAP > 25 mmHg), pulmonary capillary wedge pressure (PCWP < than 15 mmHg) and pulmonary vascular resistance (PVR > 3WU).
Cardiovascular magnetic resonance (CMR) is increasingly being used to monitor pulmonary hemodynamics and cardiac function in patients with PH. Flow-sensitive CMR sequences can be used to measure pulmonary artery flow (QP) flow and TRV, and thereby estimate PVR (PVRRHC ∝ TRV/Qp) using a method analogous to that used by Abbas et al. for echocardiography. Using this model, TRV is used as a surrogate of peak systolic pulmonary arterial pressure and reflects the trans-tricuspid gradient. With respect to cardiac size and function, CMR is considered the gold standard for quantification of left[6–8] and right[9–11] ventricular size and function with time-resolved “CINE” balanced steady-state free precession (bSSFP) imaging due to its high reproducibility.
In this study we have investigated an alternative CMR approach to assess flow and ventricular function using a single 4D flow MR sequence that uses a radially undersampled, time-resolved, 3-dimensional, 3-directionally velocity-encoded imaging scheme. Ventricular size and function were measured using the time-resolved magnitude images and compared with values obtained using standard CINE bSSFP and two-dimensional flow-sensitive (2D flow) imaging. Tricuspid valve, pulmonary artery and aorta flow was quantified using the time-resolved phase images and compared with ventricular stroke volumes obtained using volumetric methods. In addition, the ratio of TRV/QP was correlated to PVRRHC to determine if this method could also be used to estimate PVR non-invasively. A benefit of using the same sequence for measuring cardiac chamber volumes and flow includes an overall shortened examination acquisition time (particularly in cases where numerous flow measurements are required). In addition, a free-breathing acquisition is beneficial in patients who have difficulty with the numerous breath holds required using standard cardiac MR sequences. Although not explored in this study, 4D flow MR has been shown to be of use in assessing additional hemodynamic characteristics of normal and abnormal flow in a variety of cardiovascular diseases.
Acute canine thromboembolic pulmonary hypertension model
All studies were approved by the institutional local animal care and use committee (RARC). Acute PAH was induced by injection of micro-beads (150–500 μm) in the right atrium and ventricle in six adult female beagles (8.3 ± 2.4 kg). Dogs were anesthetized with propofol (10 mg/kg body weight), intubated, and maintained under anesthesia with isoflurane (1 to 3%) with 100% oxygen. During anesthesia, ventilation was adjusted to keep end-tidal CO2 within normal limits (30–50 mmHg). A femoral arterial catheter was inserted and systemic arterial pressure (SAP) and arterial blood gases were monitored. The femoral and external jugular veins were catheterized for RHC access, delivery of emboli, delivery of contrast for angiography and blood sampling.
While maintaining the dog under anesthesia, CMR and RHC were repeated twice, the first prior to induction of PH to obtain baseline data, and the second repeating all acquisitions after successful induction of thromboembolic PH (confirmed by > twofold increase in baseline mPAP). The dogs were euthanized according to the RARC protocol after the post-embolization CMR study.
Right heart catheterization
After induction of anesthesia, dogs were granted a resting period of approximately 30 minutes to ensure adaptation of the circulation to anesthesia. Baseline pulmonary arterial pressure (PAP), pulmonary capillary wedge pressure (PCWP) and right ventricular cardiac output (QP) were measured. Measurements were made using a 7.5-French, fluid-filled catheter connected to an analog pressure recorder. PVR [in Woods units (WU)] was calculated using the formula, analogous to Ohm’s law for electrical circuits: PVRRHC = ∆P/QP, where ∆P is the trans-pulmonary pressure gradient (∆P = mPAP - PCWP) and QP is the flow in the pulmonary artery measured by thermodilution.
Magnetic resonance imaging
CMR measurements were performed prior to and following induction of PH resulting in a total of 12 measurements on a 3.0 T clinical systems (MR750, GE Healthcare, Waukesha, WI). Contiguous axial CINE bSSFP slices were obtained covering the entire heart. Parameters for CINE bSSFP imaging included: 310 × 190 mm field of view, 224 × 192 acquisition matrix, 7 mm slice thickness, 0 mm gap, ±125 kHz bandwidth, 45° flip angle, TR/TE = 3.1/1.1 ms (fractional echo readout), and a prospectively gated, k-space segmented acquisition (12 views per segment), for an acquired temporal resolution of 37 ms. 20 temporally interpolated time frames were reconstructed at each slice location. Between 13 and 16 slices were acquired depending on subject anatomy. Each slice was acquired within an 8 to 10 second time interval of suspended ventilation.
2D and 4D flow CMR were performed following the administration of 0.1 mmol/kg of Gd-based intravenous contrast (gadobenate dimeglumine, Bracco Diagnostics, Inc., Princeton, NJ). A standard 2D flow CMR sequence was used to assess flow through the main, right, and left pulmonary arteries (MPA, RPA and LPA, respectively). Parameters for 2D flow imaging were: 160-220 × 160-176 mm field of view, 256 × 128 acquisition matrix, 5 mm slice thickness, TR/TE = 5.53-5.93/3.17-3.44 ms, and a prospectively gated, k-space segmented acquisition. Temporal resolution was 22.11-35.44 ms and data were interpolated to 20 reconstructed time frames for each location. 2D flow CMR acquisitions were performed during a breath-hold at end-expiration.
4D flow CMR was performed using a previously described 3D radial sequence. Parameters were: imaging volume = 320 × 320 × 220 mm, readout length = 256 samples, TR/TE = 6.7/2.4 ms, flip angle 10-20°, spatial resolution = 1.3 mm isotropic. Retrospective ECG-gating was used. Respiratory triggering with bellows allowed for breathing compensation with an acceptance rate of 50%, resulting in a scan duration of 10–12 min during free breathing. Data were retrospectively sorted into 20 time frames according to their position in the cardiac cycle. Subsequently, image reconstruction was performed utilizing a compressed sensing reconstruction and a temporal filter for view sharing. To minimize the number of slices needing manual segmentation for measuring ventricular volumes, three contiguous slices were averaged. In addition, time-resolved 4D flow CMR magnitude images were reformatted into the left-ventricular short-axis (SA) orientation using home-built MatLab software, for a second, more convenient LV segmentation.
RV and LV function
RV and LV end-diastolic (EDV) and end-systolic (ESV) volumes were determined from manually segmented contours of end-diastolic and end-systolic bSSFP and 4D flow CMR images (Figure 1), respectively. Segmentation of bSSFP images was accomplished using ReportCard (GE Healthcare, Waukesha, WI). Segmentation of 4D flow CMR images was conducted using Osirix (Pixmeo, Geneva, Switzerland). Stroke volume (SV = EDV-ESV) and ejection fraction (EF = SV/EDV) were also determined for both bSSFP and 4D flow CMR.
Analysis of the 2D flow CMR data was performed with CV flow (Medis, Leiden, the Netherlands). The MPA, RPA, and LPA were manually segmented at each phase of the cardiac cycle. Net flow per heart beat was recorded for each location. Average velocity was recorded for the MPA only. PVR was estimated from 2D flow and 2D CINE bSSFP data using the following equation PVRCMR (WU) = 19.38 – [4.62 · lnPA average velocity] – [0.08 · RVEF].
Quantitative flow analysis of the 4D flow CMR datasets was performed with a previously described Matlab-based home built software (The MathWorks, Natick, MA, USA). Two-dimensional cutplanes were manually placed perpendicular to the direction of the flow in the tricuspid valve, MPA, RPA, LPA (Figure 2) and aorta (Figure 3) using Ensight (CEI, Apex, NC, USA). Specifically, the tricuspid valve plane was generated using flow streamlines, which allowed for localization of the peak tricuspid regurgitation velocity with no need for annular motion compensation. From these analyses, we recorded (a) the peak tricuspid valve regurgitation velocity (TRV) for the generated two-dimensional tricuspid valve cutplane; (b) net flow, time-to-peak flow, and time-to-peak acceleration in the MPA, RPA, and LPA; and (c) net flow in the aorta. To assess the internal consistency of the flow measurements, we determined the relative flow through the RPA and LPA, as a percentage of MPA flow, and the ratio of pulmonary (QP) to aorta (QS) flow.
Values are reported as mean ± standard deviation. Differences between flow (2D and 4D) and CINE bSSFP and between the two flow techniques were assessed using Bland-Altman analysis. Linear regression analysis was used to determine the Pearson correlation coefficients between TRV/QP and PVRRHC and mPAP. A linear regression equation was derived to calculate PVR4Dflow from the imaging parameters. Bland-Altman analysis was used to assess the differences between the PVRRHC and PVR4Dflow and PVRCMR. Values obtained prior to and following induction of PH were compared used a paired Student’s t-test. P-values less than 0.05 were considered statistically significant.
RV and LV function
Mean values for RV and LV volumes using CINE bSSFP and 4D flow CMR are summarized in Table 1. Differences in RV EDV, ESV, and SV were not significant (P > 0.05). LV EDV determined with both axial and short-axis reconstructions did not significantly differ from that measured using axial CINE bSSFP (P > 0.05). LV ESV and SV were significantly larger with 4D flow CMR than CINE bSSFP using both axial and short-axis reconstructions (P < 0.05). LV EDV and SV were smaller (P = 0.02 and 0.001, respectively) and LV ESV was larger (P = 0.10) when calculated from the short-axis reconstructions than the axial reconstructions.
Results of the flow analysis in the MPA, RPA and LPA are summarized in Tables 2 and3. One pre- and one post-embolization 4D flow datasets were excluded due to inadequate imaging quality. Flows were significantly higher with 2D flow than with 4D flow. Relative flow through the RPA and LPA was 61% and 39%, respectively, with 4D flow and 59% and 41%, respectively, with 2D flow. The average flow through the RPA and LPA combined (12.6 ± 3.5 mL/cycle with 4D flow and 18.4 ± 6.4 mL/cycle with 2D flow) was not significantly different from the flow measured through the MPA (12.0 ± 3.9 mL/cycle with 4D flow, P = 0.96, and 18.9 ± 3.1 mL/cycle with 2D flow, P = 0.70).
Comparing flow data pre- and post-embolization, the differences in MPA, RPA, and LPA flows were not significant. Although there was a trend toward faster time-to-peak flow and time-to-peak acceleration values post-embolization in the pulmonary arteries, these changes were not significant (Table 3). Time-to-peak acceleration was shorter post-embolization than pre-embolization, although this difference was not significant (Table 3).
The average aortic flow (QS) using 4D flow CMR was 11.6 ± 4.0 mL/cycle (P = 0.06 with bSSFP and 0.08 with 4D flow CMR volumetry). The mean difference between QS with 4D flow CMR and LV SV calculated from 4D flow CMR was -0.7 mL. The difference between QS and MPA flow (QP) calculated from 4D flow CMR phase-sensitive data was not statistically significant (P = 0.30). Furthermore, the average QP/QS was 0.96 ± 0.11 (range: 0.87-1.21). Differences in aortic flow pre- and post-embolization, 12.0 ± 3.8 mL/cycle and 11.1 ± 4.4 mL/cycle, respectively, did not reach statistical significance (P = 0.81).
Tricuspid regurgitation was detected in all canines pre- and post-embolization, enabling calculation of PVRCMR in all cases. PVRRHC values pre- and post-embolization were 2.4 ± 0.9 WU and 9.8 ± 5.7 WU, respectively (Figure 4). The Pearson correlation coefficient between TRV/QP and PVRRHC was 0.95 for combined pre- and post-embolization data (Figure 5). When analysing the data pre- and post-embolization separately, the correlation between TRV/QP and PVRRHC was higher post-embolization (R = 0.99) than pre-embolization (R = 0.26). The Pearson correlation coefficient between TRV/QP and mPAP was 0.65. The equation for calculating PVR4Dflow using the combined data was PVR4Dflow (WU) = 0.16(TRV/QP) – 7.25. Using this equation, PVR4Dflow values pre- and post-embolization were 3.3 ± 1.8 WU (P = 0.35 for PVR4Dflow vs. PVRRHC) and 9.0 ± 6.2 WU (p = 0.15 for PVR4Dflow vs. PVRRHC), respectively (Figure 4). In one case, PVR4Dflow was lower post embolization. The mean difference between PVR4Dflow and PVRRHC was 0 (WU) with positive and negative levels of agreement of 3.52 and -3.52 (WU), respectively (Figure 6).
Using 2D flow and 2D CINE bSSFP, PVRCMR was 5.0 ± 1.6 (WU) pre-embolization (P = 0.051 with PVRRHC and 0.21 with PVR4Dflow) and 7.5 ± 1.3 (WU) post-embolization (P = 0.38 with PVRRHC and 0.66 with PVR4Dflow) (Figure 4). PVRCMR increased post-embolization in all cases. The Pearson correlation coefficients between PVRCMR and PVRRHC was 0.78 for combined pre- and post-embolization data, 0.03 for pre-embolization data and 0.91 for post-embolization data.. The mean difference between PVRCMR and PVRRHC was 0.18 (WU) with positive and negative levels of agreement of 8.4 and -8.0 (WU), respectively.
In this study we demonstrated the feasibility of measuring RV and LV function and pulmonary vascular resistance from a single, free breathing 4D flow CMR sequence. RV and LV volumes were measured using the time-resolved magnitude data from the same 4D flow CMR acquisition. Using the three-directional velocity information, we found a strong correlation between the ratio of the peak tricuspid regurgitation velocity to the flow through the pulmonary arteries (TRV/QP) and PVR determined at RHC, using an approach that is analogous to that used for echocardiography.
RV and LV function
Cardiac magnetic resonance using breath-hold CINE balanced steady-state free precession imaging (bSSFP) is considered the gold standard for quantification of ventricular size and function[6–11]. However, in patients with PH, dyspnea can be severe, and the multiple breath-holds required to scan the entire heart can be difficult, or impossible, for many patients. An alternative approach is to use free-breathing acquisition methods with respiratory navigation or triggering to coordinate data acquisition with the respiration, as was done here with the 4D flow MR approach. Using a volumetric (3D) approach to assessing cardiac function would also be beneficial in that it would enable quantification of RV and left ventricular (LV) volumes simultaneously, which is important because the most common cause of pulmonary hypertension is left heart disease. Quantification of LV and RV function using time-resolved 3D acquisitions has previously been demonstrated using both free-breathing and breath-hold bSSFP techniques[20, 21]. Although free-breathing approaches to cardiac function analysis are longer than single breath-hold acquisitions, in our clinical experience, patients with PH frequently have difficulty holding their breath for even more than 10 seconds. Future studies in PH patients comparing these two approaches to whole-heart functional analysis are warranted. LV ESV was larger, and therefore LV SV was smaller, with time-resolved magnitude reconstructions from the 4D flow CMR dataset than with the standard 2D bSSFP acquisition. This may be partially related to the higher spatial resolution with the 4D flow CMR data, resulting in improved delineation of the blood-myocardial interface, which is frequently a challenge when using standard 2D bSSFP. The results from our study are encouraging in that the bias we observed using a 3D radially undersampled approach was of the same magnitude as previously published studies.
Estimation of PVR
Currently, Doppler echocardiography is used clinically in patients with PH to estimate PVR from the ratio of the TRV to the velocity time integral in the right ventricular outflow tract (VTIRVOT), which is a surrogate of the flow through the pulmonary artery (QP). Abbas et al. reported a strong linear correlation between TRV/VTIRVOT and PVRRHC and proposed calculating PVR from the following equation: PVR = 10 × TRV/VTIRVOT in Woods units (WU). However, a potential important limitation of this approach is that Doppler ultrasound does not measure flow directly, but rather the velocity time integral which does not take into account the area of the vessel of interest. Furthermore, echocardiography-based methods of assessing PVR can be limited by body habitus or other anatomical factors obstructing the acoustic window. In contrast, CMR-based flow sequences can be used to measure flow, in addition to velocity, without the limitations of echocardiography. As a result, other investigators have used CMR to estimate PVR, generating similar models based on cardiac function and pulmonary artery flow. Garcia-Alvarez et al. derived an equation for calculating PVR using the average MPA velocity and the RVEF. This model was derived in a cohort of human PH subjects and subsequently validated in a second cohort and in a more recent canine model of acute and chronic thromboembolic PH. Swift et al. created a model to estimate PVR based on left atrial size (as an estimate of PCWP) and a linear regression equation for calculating mPAP from the interventricular septal angle and ventricular mass index. As with our study, limits of agreement for determining PVR using these models were also fairly wide, -6.0 to 4.9 WU (Garcia-Alvarez) and -5.1 to 4.6 WU (Swift). Although the bias we observed for PVRCMR method using 2D flow and 2D CINE bSSFP was minimal (0.18 WU), the limits of agreement were fairly wide, wider than what we observed using the TRV/VTIRVOT method and wider than what was reported by Garcia-Alvarez et al.
Flow quantification with 4D flow
4D flow CMR has been used previously to evaluate PH. Qualitatively assessing flow patterns in the MPA, Reiter et al. reported that in patients with pulmonary arterial hypertension abnormal separation between the boundary layer next to the wall of the MPA resulted in abnormal vortex development in the MPA. In addition, they reported delayed acceleration times in the MPA. These hemodynamic phenomena should influence PA-RV interactions by changing both PVR and RV function. Interestingly, we did not detect a significant difference in time-to-peak acceleration. In fact, we found an earlier time to peak flow after pulmonary artery embolization and induction of PH. The different observations between this study and that of the Reiter group are presumably related to differences in the etiology of PH, acute in the current study and chronic in the study by Reiter et al..
A benefit of using 4D flow CMR to quantify blood flow is that quantification can be performed of any vessel within the imaging volume a posteriori. In cases where flow quantification is performed in numerous vessels, the overall scan time is typically less using accelerated 4D flow CMR methods than using standard 2D flow CMR techniques. The stroke volumes measured with the 4D flow CMR flow data were slightly lower than the stroke volumes determined using Cine bSSFP. MPA flow volumes with 4D flow CMR were substantially lower than values observed with 2D flow CMR. Frydrychowicz et al. also reported underestimation of Qp and Qs using the same 4D flow CMR study as that used in the current study. Some of these differences may be related to differences in temporal resolution and differences in flow during free breathing (4D flow) and end-expiration (2D flow).
Limitations of the current study include the fact that the reference standard used for LV volumes were axial CINE bSSFP images. Although the LV short-axis orientation is typically used to calculate LV volumes, Fratz et al. recently reported that axial images were more reproducible than short-axis images for calculating LV volumes in patients with repaired tetralogy of Fallot. Future studies comparing short-axis reconstructed 4D flow CMR data to short-axis 2D CINE bSSFP will have to be conducted to confirm the small bias observed in this study.
Other limitations of this study include (a) a small sample size without a cohort to validate the model, (b) the estimation of the trans-pulmonary pressure gradient assumes that the modified Bernoulli equation is valid and requires the presence of tricuspid valve regurgitant jet, and (c) the fact that the equation for calculating PVR is based on RHC data obtained during acute embolic PH, which will most likely not be applicable to chronic PAH. In acute pulmonary embolism, right atrial pressures are elevated, reducing the RA-RV pressure gradient and, therefore, lowering the peak TRV. Although acute pulmonary embolism is the most common cause of acute PH, it is not clear from this study if these results would be transferable to acute pulmonary thromboembolism in humans or to other etiologies of PH. Subsequent studies will have to be conducted to assess the relationship between TRV/QP in other causes of PH, including in humans.
In conclusion, this study demonstrates the feasibility of comprehensively estimating pulmonary vascular resistance, assessing ventricular function and pulmonary artery hemodynamics in a canine model of acute embolic PH using 4D flow CMR. Although not explored in this study, another potential benefit of using 4D flow CMR to assess PH is the derivation of additional hemodynamic parameters such as wall shear stress[16, 18, 27, 28] and pulse wave velocities[29, 30], which could provide further insights into pulmonary artery remodeling and interactions between pulmonary arterial stiffening and RV dysfunction.
D’Alonzo GE, Barst RJ, Ayres SM, Bergofsky EH, Brundage BH, Detre KM, Fishman AP, Goldring RM, Groves BM, Kernis JT, Levy PS, Pietra GG, Reid LM, Reeves JT, Rich S, Vreim CE, Williams GW, Wu M: Survival in patients with primary pulmonary hypertension. Results from a national prospective registry. Ann Intern Med. 1991, 115: 343-9. 10.7326/0003-4819-115-5-343.
Rudski LG, Lai WW, Afilalo J, Hua L, Handschumacher MD, Chandrasekaran K, Solomon SD, Louie EK, Schiller NB: Guidelines for the echocardiographic assessment of the right heart in adults: a report from the American society of echocardiography endorsed by the European association of echocardiography, a registered branch of the European society of cardiology, and the Canadian society of echocardiography. J Am Soc Echocardiogr. 2010, 23: 685-713. 10.1016/j.echo.2010.05.010. quiz 786–688
Shah SJ: Pulmonary hypertension. JAMA. 2012, 308: 1366-74. 10.1001/jama.2012.12347.
Simonneau G, Robbins IM, Beghetti M, Channick RN, Delcroix M, Denton CP, Elliott CG, Gaine SP, Gladwin MT, Jing ZC, Krowka MJ, Langleben D, Nakanishi N, Souza R: Updated clinical classification of pulmonary hypertension. J Am Coll Cardiol. 2009, 54: S43-54. 10.1016/j.jacc.2009.04.012.
Abbas AE, Fortuin FD, Schiller NB, Appleton CP, Moreno CA, Lester SJ: A simple method for noninvasive estimation of pulmonary vascular resistance. J Am Coll Cardiol. 2003, 41: 1021-7. 10.1016/S0735-1097(02)02973-X.
Clay S, Alfakih K, Messroghli DR, Jones T, Ridgway JP, Sivananthan MU: The reproducibility of left ventricular volume and mass measurements: a comparison between dual-inversion-recovery black-blood sequence and SSFP. Eur Radiol. 2006, 16: 32-7. 10.1007/s00330-005-2853-3.
Hudsmith LE, Petersen SE, Francis JM, Robson MD, Neubauer S: Normal human left and right ventricular and left atrial dimensions using steady state free precession magnetic resonance imaging. J Cardiovasc Magn Reson. 2005, 7: 775-82. 10.1080/10976640500295516.
Maceira AM, Prasad SK, Khan M, Pennell DJ: Normalized left ventricular systolic and diastolic function by steady state free precession cardiovascular magnetic resonance. J Cardiovasc Magn Reson. 2006, 8: 417-26. 10.1080/10976640600572889.
Grothues F, Moon JC, Bellenger NG, Smith GS, Klein HU, Pennell DJ: Interstudy reproducibility of right ventricular volumes, function, and mass with cardiovascular magnetic resonance. Am Heart J. 2004, 147: 218-23. 10.1016/j.ahj.2003.10.005.
Catalano O, Antonaci S, Opasich C, Moro G, Mussida M, Perotti M, Calsamiglia G, Frascaroli M, Baldi M, Cobelli F: Intra-observer and interobserver reproducibility of right ventricle volumes, function and mass by cardiac magnetic resonance. J Cardiovasc Med (Hagerstown). 2007, 8: 807-14. 10.2459/JCM.0b013e32801105ef.
Maceira AM, Prasad SK, Khan M, Pennell DJ: Reference right ventricular systolic and diastolic function normalized to age, gender and body surface area from steady-state free precession cardiovascular magnetic resonance. Eur Heart J. 2006, 27: 2879-88. 10.1093/eurheartj/ehl336.
Gu T, Korosec FR, Block WF, Fain SB, Turk Q, Lum D, Zhou Y, Grist TM, Haughton V, Mistretta CA: PC VIPR: a high-speed 3D phase-contrast method for flow quantification and high-resolution angiography. AJNR Am J Neuroradiol. 2005, 26: 743-9.
Markl M, Frydrychowicz A, Kozerke S, Hope M, Wieben O: 4D flow MRI. J Magn Reson Imaging. 2012, 36: 1015-36. 10.1002/jmri.23632.
Wood P: Pulmonary hypertension. Br Med Bull. 1952, 8: 348-53.
Frydrychowicz A, Wieben O, Niespodzany E, Reeder SB, Johnson KM, Francois CJ: Quantification of thoracic blood flow using volumetric magnetic resonance imaging with radial velocity encoding: in vivo validation. Invest Radiol. 2013, 48: 819-25. 10.1097/RLI.0b013e31829a4f2f.
Liu J, Wieben O, Jung Y, Samsonov AA, Reeder SB, Block WF: Single breathhold cardiac CINE imaging with multi-echo three-dimensional hybrid radial SSFP acquisition. J Magn Reson Imaging. 2010, 32: 434-40. 10.1002/jmri.22269.
Garcia-Alvarez A, Fernandez-Friera L, Mirelis JG, Sawit S, Nair A, Kallman J, Fuster V, Sanz J: Non-invasive estimation of pulmonary vascular resistance with cardiac magnetic resonance. Eur Heart J. 2011, 32: 2438-45. 10.1093/eurheartj/ehr173.
Stalder AF, Russe MF, Frydrychowicz A, Bock J, Hennig J, Markl M: Quantitative 2D and 3D phase contrast MRI: optimized analysis of blood flow and vessel wall parameters. Magn Reson Med. 2008, 60: 1218-31. 10.1002/mrm.21778.
Oudiz RJ: Pulmonary hypertension associated with left-sided heart disease. Clin Chest Med. 2007, 28: 233-41. 10.1016/j.ccm.2006.12.001.
Greil GF, Boettger T, Germann S, Klumpp B, Baltes C, Kozerke S, Bialkowski A, Urschitz MS, Miller S, Wolf I, Meinzer HP, Sieverding L: Quantitative assessment of ventricular function using three-dimensional SSFP magnetic resonance angiography. J Magn Reson Imaging. 2007, 26: 288-95. 10.1002/jmri.20967.
Peters DC, Korosec FR, Grist TM, Block WF, Holden JE, Vigen KK, Mistretta CA: Undersampled projection reconstruction applied to MR angiography. Magn Reson Med. 2000, 43: 91-101. 10.1002/(SICI)1522-2594(200001)43:1<91::AID-MRM11>3.0.CO;2-4.
Garcia-Alvarez A, Fernandez-Friera L, Garcia-Ruiz JM, Nuno-Ayala M, Pereda D, Fernandez-Jimenez R, Guzman G, Sanchez-Quintana D, Alberich-Bayarri A, Pastor-Escuredo D, Sanz-Rosa D, Garcia-Preito J, Gonzalez-Mirelis JG, Pizarro G, Jimenez-Borreguero LJ, Fuster V, Sanz J, Ibanez B: Noninvasive monitoring of serial changes in pulmonary vascular resistance and acute vasodilator testing using cardiac magnetic resonance. J Am Coll Cardiol. 2013, 62: 1621-31. 10.1016/j.jacc.2013.07.037.
Swift AJ, Rajaram S, Hurdman J, Hill C, Davies C, Sproson TW, Morton AC, Capener D, Elliot C, Condliffe R, Wild JM, Kiely DG: Noninvasive estimation of PA pressure, flow, and resistance with CMR imaging: derivation and prospective validation study from the ASPIRE registry. JACC Cardiovasc Imaging. 2013, 6: 1036-47. 10.1016/j.jcmg.2013.01.013.
Reiter G, Reiter U, Kovacs G, Kainz B, Schmidt K, Maier R, Olschewski H, Rienmueller R: Magnetic resonance-derived 3-dimensional blood flow patterns in the main pulmonary artery as a marker of pulmonary hypertension and a measure of elevated mean pulmonary arterial pressure. Circ Cardiovasc Imaging. 2008, 1: 23-30. 10.1161/CIRCIMAGING.108.780247.
Fratz S, Schuhbaeck A, Buchner C, Busch R, Meierhofer C, Martinoff S, Hess J, Stern H: Comparison of accuracy of axial slices versus short-axis slices for measuring ventricular volumes by cardiac magnetic resonance in patients with corrected tetralogy of fallot. Am J Cardiol. 2009, 103: 1764-9. 10.1016/j.amjcard.2009.02.030.
Hui-li G: The management of acute pulmonary arterial hypertension. Cardiovasc Ther. 2011, 29: 153-75. 10.1111/j.1755-5922.2009.00095.x.
Bieging ET, Frydrychowicz A, Wentland A, Landgraf BR, Johnson KM, Wieben O, Francois CJ: In vivo three-dimensional MR wall shear stress estimation in ascending aortic dilatation. J Magn Reson Imaging. 2011, 33: 589-97. 10.1002/jmri.22485.
Frydrychowicz A, Berger A, Russe MF, Stalder AF, Harloff A, Dittrich S, Hennig J, Langer M, Markl M: Time-resolved magnetic resonance angiography and flow-sensitive 4-dimensional magnetic resonance imaging at 3 Tesla for blood flow and wall shear stress analysis. J Thorac Cardiovasc Surg. 2008, 136: 400-7. 10.1016/j.jtcvs.2008.02.062.
Markl M, Wallis W, Brendecke S, Simon J, Frydrychowicz A, Harloff A: Estimation of global aortic pulse wave velocity by flow-sensitive 4D MRI. Magn Reson Med. 2010, 63: 1575-82. 10.1002/mrm.22353.
Wentland AL, Wieben O, Francois CJ, Boncyk C, Munoz Del Rio A, Johnson KM, Grist TM, Frydrychowicz A: Aortic pulse wave velocity measurements with undersampled 4D flow-sensitive MRI: comparison with 2D and algorithm determination. J Magn Reson Imaging. 2013, 37: 853-9. 10.1002/jmri.23877.
This study was supported by funding from the following funding sources: National Institutes of Health: R01HL072260, R01HL105598, R01HL086939.
Department of Radiology Research and Development Fund.
The project was supported by the Clinical and Translational Science Award (CTSA) program, previously through the National Center for Research Resources (NCRR) grant 1UL1RR025011, and now by the National Center for Advancing Translational Sciences (NCATS), grant 9U54TR000021.
A.F. received an educational stipend from Bracco Diagnostics.
The Departments of Radiology and Medical Physics receive support by GE Healthcare.
AR – Study design, data acquisition and analysis, manuscript preparation and editing. AF – Study design, data acquisition, manuscript editing. KMJ – Study design, manuscript editing. HK – Study design, data acquisition, manuscript editing. NCC – Study design, data acquisition, manuscript editing. OW – Study design, manuscript editing. CJF – Study design, data acquisition and analysis, manuscript preparation and editing. All authors read and approved the final manuscript.
Authors’ original submitted files for images
About this article
Cite this article
Roldán-Alzate, A., Frydrychowicz, A., Johnson, K.M. et al. Non-invasive assessment of cardiac function and pulmonary vascular resistance in an canine model of acute thromboembolic pulmonary hypertension using 4D flow cardiovascular magnetic resonance. J Cardiovasc Magn Reson 16, 23 (2014). https://doi.org/10.1186/1532-429X-16-23
- Pulmonary hypertension
- Pulmonary vascular resistance
- Heart function
- 4D flow cardiovascular magnetic resonance
- Thromboembolic pulmonary hypertension
- Canine model