2066

Associations of Cardiac Rhythm with left atrial and left atrial appendage hemodynamics measured with 4D flow in Atrial Fibrillation
Amanda L DiCarlo1, Justin Baraboo1,2, Patrick McCarthy3, Rishi Arora4, Rod Passman4, Philip Greenland4, Daniel C Lee1,4, Daniel Kim1,2, and Michael Markl1,2
1Radiology, Northwestern University, Chicago, IL, United States, 2Biomedical Engineering, Northwestern University, Chicago, IL, United States, 3Cardiac Surgery, Northwestern University, Chicago, IL, United States, 4Cardiology, Northwestern University, Chicago, IL, United States

Synopsis

Atrial fibrillation (AF) is the most common cardiac arrhythmia, characterized by irregular heart rhythm which can lead to stroke due to thrombus formation in the left atrium (LA) or left atrial appendage (LAA). Reduced velocities and increased stasis measured with 4D flow MRI are risk metrics correlated with LA thrombus formation and potential predictors of stroke. This study analyzed the relationship between heart rate variability derived from a real time MRI acquisition and 4D flow hemodynamic risk metrics. LA peak velocity and stasis were significantly correlated with heart rate variability, suggesting a link between rhythm status and thrombus risk.

Introduction

Atrial fibrillation (AF) is the most common cardiac arrhythmia1 which can result in stroke due to thrombus formation in the left atrium (LA). Left atrial hemodynamics, including peak velocity and stasis, may correlate with risk for atrial thrombus formation and be useful for improving stroke prediction.2,3 4D flow MRI has been previously used to quantify LA flow, demonstrating increased stasis and decreased peak velocity in patients with AF.4,5 However, the impact of cardiac rhythm, or the level of heart rate variability, on 4D flow measured velocity and stasis has not been systematically investigated. The goal of this study was to analyze the relationship between cardiac rhythm and left atrial and left atrial appendage (LAA) peak velocity and stasis in patients with atrial fibrillation.

Methods

53 patients with AF underwent a pulmonary vein mapping MRI exam prior to an ablation procedure, including 4D flow (spatial resolution = 1.9-3.0 mm3, temporal resolution = 41-77 ms, venc = 50-100 cm/s, 1.5 and 3T MRI systems) and 3D contrast-enhanced (CE) MRA (spatial resolution = 1.3-1.8 mm3). A custom radial real time sequence (real time cine/real time phase contrast6 acquired in the LA) with continuous acquisition was also acquired and used for extracting RR interval information from the raw k-space data.

4D flow MRI processing: Following pre-processing (phase offset correction, velocity anti-aliasing), a 3D phase contrast MRA (PC-MRA) was generated from the 4D flow data and used to generate a manual segmentation of the LA (Mimics, Materialise). The LA and LAA were also segmented from the high-resolution CE-MRA. The CE-MRA LA segmentation was registered with the 4D flow LA segmentation7 and applied to both LA and LAA to be used for flow quantification. LA and LAA stasis (% time velocity < 0.1 m/s over cardiac cycle for a voxel) and peak velocity (calculated as the mean of the top 5% of velocities) were calculated. (Figure 1A)

RR-interval Derivation: The scanner patient monitoring unit (PMU) records timestamps and trigger information and embeds them into the raw data of the real-time sequence. Since radial k-space data is continuously acquired, the PMU information served as a continuous record of RR-intervals collected during the same MRI exam. Histograms of RR-interval length were generated for each patient (Figure 1B) and the mean and standard deviation were obtained. Patients were classified into low and high heart rate variability (HRV) groups based on their RR coefficient of variation (CV=RR standard deviation/mean RR interval) using the CV test method reported by Tateno and Glass.8 Rhythm status during the exam was classified as arrhythmia/high HRV for a CV = 0.156-0.3248 (n=27), while a CV below this range was classified as low HRV/sinus rhythm (n=21). For the group comparison, patients with CV>0.324 (n=5) were excluded.

Results

Figure 2 shows plots of LA stasis and peak velocity versus RR standard deviation. Both LA stasis (r=0.33, p<0.05) and LA peak velocity (r=-0.24, p<0.05) were correlated with RR standard deviation, a measure of RR variability. LAA stasis (r=0.12, p=0.38) and LAA peak velocity (r=-0.23, p=0.07) were not significantly correlated with RR variability. Between group comparisons of LA and LAA peak velocity and stasis are shown in Figure 3. LA stasis was significantly increased in the high HRV group compared to the low HRV group (t-test, p<0.05, high HRV: 52±14%, low HRV: 44±9%), and LA peak velocity was significantly decreased (t-test, p<0.05, high HRV: 32±6cm/s, low HRV: 35±5cm/s). Differences in LAA stasis and peak velocity between groups were non-significant.

Discussion and Conclusions

The results presented show a trend towards slower LA velocities and elevated LA stasis with higher heart rate variability, suggesting a link between rhythm status and thrombosis risk. 4D flow is an averaged technique and does not capture beat-to-beat variations in flow metrics, however an overall increase in slow flow is observed as the RR-interval distribution disperses and variability increases. The present study did not consider potentially confounding factors associated with disease progression, such as LA volume or morphology. LAA velocity and stasis were not significantly correlated with HRV. This could be attributed to increased measurement noise and difficulty detecting the small LAA during segmentation (mitigated by using CE-MRA in the present study), or reflective of physiological differences between the LA and LAA hemodynamics with arrhythmia. Future work including an expanded patient cohort and follow-up (post ablation) MRI data will further investigate the relationship between rhythm status and hemodynamics in AF.

Acknowledgements

No acknowledgement found.

References

1. January CT, Wann LS, Alpert JS, Calkins H, Cigarroa JE, Cleveland JC Jr, et al. 2014 AHA/ACC/HRS guideline for the management of patients with atrial fibrillation. J Am Coll Cardiol. 2014; 64(21):2246–80.

2. Goldman ME, Pearce LA, Hart RG, Zabalgoitia M, Asinger RW, Safford R, Halperin JL. Pathophysiologic correlates of thromboembolism in nonvalvular atrial fibrillation: I. Reduced flow velocity in the left atrial appendage (The Stroke Prevention in Atrial Fibrillation [SPAF-III] study). J Am Soc Echocardiogr. 1999; 12:1080–1087.

3. Fatkin D, Kelly RP, Feneley MP. Relations between left atrial appendage blood flow velocity, spontaneous echocardiographic contrast and thromboembolic risk in vivo. J Am Coll Cardiol. 1994; 23:961–969.

4. Markl M, Lee DC, Ng J, Carr M, Carr J and Goldberger JJ. Left Atrial 4D Flow MRI: Stasis and Velocity Mapping in Patients with Atrial Fibrillation. Invest Radiol. 2016; 51:147.

5. Markl M, Lee DC, Furiasse N, et al. Left Atrial and Left Atrial Appendage 4D Blood Flow Dynamics in Atrial Fibrillation. Circ Cardiovasc Imaging. 2016; 9(9):e004984.

6. Haji-Valizadeh H, Feng L, Ma LE, Shen D, Block KT, Robinson JD, Markl M, Rigsby CK, Kim D. Highly accelerated, real-time phase-contrast MRI using radial k-space sampling and GROG-GRASP reconstruction: a feasibility study in pediatric patients with congenital heart disease. NMR Biomed 2020; 33(5):e4240.

7. Jenkinson M, Smith S. A global optimisation method for robust affine registration of brain images. Med Image Anal. 2001; 5(2):143-156.

8. Tateno and Glass. Automatic detection of atrial fibrillation using the coefficient of variation and density histograms of RR and ∆RR Intervals. Med. Biol. Eng. Comput. 2001; 39:664-671

Figures

Figure 1: A) Workflow for segmentation and registration of 4D flow and CE-MRA data. From left to right: the LA was segmented from the 4D flow PC-MRA and the LA and LAA were segmented from the CE-MRA. The left atria are registered together, and the resulting transform is applied to the LA and LAA of the CE-MRA. Flow quantifications (stasis, peak velocities) are then quantified. B) Sample of the RR signal extracted from a continuously acquired real time phase contrast sequence and RR histogram

Figure 2: Correlation plots of 4D flow derived LA (left) and LAA (right) mean stasis (top row) and peak velocity (bottom row) vs RR interval standard deviation obtained from a separate real time acquisition.

Figure 3: Box plots illustrating differences in mean stasis (top row) and peak velocity (bottom row) for the LA (left) and LAA (right) between low (n=21) and high (n=27) heart rate variability groups. *p<0.05

Proc. Intl. Soc. Mag. Reson. Med. 29 (2021)
2066