- Open Access
Cardiovascular magnetic resonance guided ablation and intra-procedural visualization of evolving radiofrequency lesions in the left ventricle
Journal of Cardiovascular Magnetic Resonance volume 20, Article number: 20 (2018)
Radiofrequency (RF) ablation has become a mainstay of treatment for ventricular tachycardia, yet adequate lesion formation remains challenging. This study aims to comprehensively describe the composition and evolution of acute left ventricular (LV) lesions using native-contrast cardiovascular magnetic resonance (CMR) during CMR-guided ablation procedures.
RF ablation was performed using an actively-tracked CMR-enabled catheter guided into the LV of 12 healthy swine to create 14 RF ablation lesions. T2 maps were acquired immediately post-ablation to visualize myocardial edema at the ablation sites and T1-weighted inversion recovery prepared balanced steady-state free precession (IR-SSFP) imaging was used to visualize the lesions. These sequences were repeated concurrently to assess the physiological response following ablation for up to approximately 3 h. Multi-contrast late enhancement (MCLE) imaging was performed to confirm the final pattern of ablation, which was then validated using gross pathology and histology.
Edema at the ablation site was detected in T2 maps acquired as early as 3 min post-ablation. Acute T2-derived edematous regions consistently encompassed the T1-derived lesions, and expanded significantly throughout the 3-h period post-ablation to 1.7 ± 0.2 times their baseline volumes (mean ± SE, estimated using a linear mixed model determined from n = 13 lesions). T1-derived lesions remained approximately stable in volume throughout the same time frame, decreasing to 0.9 ± 0.1 times the baseline volume (mean ± SE, estimated using a linear mixed model, n = 9 lesions).
Combining native T1- and T2-based imaging showed that distinctive regions of ablation injury are reflected by these contrast mechanisms, and these regions evolve separately throughout the time period of an intervention. An integrated description of the T1-derived lesion and T2-derived edema provides a detailed picture of acute lesion composition that would be most clinically useful during an ablation case.
Ventricular tachycardia (VT) ablation is now frequently performed, with the rate of use of these procedures increasing substantially in the last decade . Recurrence of VT after a single ablation procedure remains high, with 35% of patients who receive initially successful ablation treatment later experiencing VT recurrence during 6–23 month follow-up periods . VT recurrence after an ablation procedure is complex and may be related to the inability to localize the critical circuit, or, when the circuit is defined, inability to obliterate the critical channels . Once the critical isthmus of a VT circuit is isolated, radiofrequency (RF) ablation is performed with the hope of creating a necrotic lesion (permanent injury) at the putative isthmus, rendering VT non-inducible. The insult of ablation also leads to edema (reversible injury) surrounding the ablation site [4,5,6,7,8]. This reversible injury is thought to result in transient conduction block, with conduction recovering once the edema has resorbed [9,10,11], potentially leading to late arrhythmia recurrence. The ability to detect whether the arrhythmogenic substrate has been permanently destroyed, as evidenced by the presence of a lesion at a critical ablation site, may be invaluable and provide an additional intra-procedural endpoint to gauge long-term procedural success.
Cardiovascular magnetic resonance (CMR)-based identification of the critical isthmus sites may be used to guide ablation toward these targets , and the therapeutic lesions themselves can be directly visualized using native (non-contrast enhanced) CMR [4,5,6,7,8, 10, 13,14,15,16].
Ablation lesions visualized using native T1 reflect injury that is associated with lethal heating , that persists at least for 3 weeks , and that correlates to chronic lesions , therefore likely represents permanent injury. Conversely, myocardial edema is visualized with native T2  and is transient, largely resolved within 3 months as seen in follow-up T2-weighted imaging after atrial ablation . The T1-derived lesion and T2-derived edema are individually significant in clinical ablations as the distributions of these injured regions could determine whether electrical block will remain permanently or resolve after healing (leading to recurrent VT). To evaluate these regions, this study includes concurrent native T1- and T2-based imaging to compare their relative extents and to construct a comprehensive understanding of the acute lesion composition.
A major potential value of intra-procedural CMR imaging lies in the ability to directly visualize and interpret the functional effect of therapeutic ablation lesions during the period when intra-procedural modification is possible. It is therefore important to identify how different aspects of CMR information evolve in this time window and how this information might be related to conventional electrophysiology (EP) endpoints. This study aims to characterize the evolution of lesions, as previously suggested to take place , but using concurrent imaging of both the T2-derived edema and T1-derived lesion to provide volumetric measurements of each. We hypothesize that the acute ablation-induced T2-derived edema evolves during the acute time frame while the T1-derived lesion remains at a consistent size. Furthermore, adding T2 mapping to this imaging framework provides unambiguous quantification of edema severity. These techniques to study the acute lesion composition build upon prior work studying lesion characteristics via the surrogate measure of gadolinium kinetics . This study aims to show temporal characteristics of the native T1-derived lesion and T2-derived edema within a time frame of particular relevance to direct CMR guidance of ablative therapies, extending the current understanding of lesion evolution previously established using native-contrast CMR  to an earlier time frame. Insights from lesion evolution could be applied to both CMR-guided and traditional ablation cases.
CMR-guided RF ablation was performed in vivo in 12 healthy Yorkshire swine (58 ± 18 kg). The Animal Care Committee of Sunnybrook Research Institute approved all protocols. All animals received an intramuscular injection of ketamine (33 mg/kg) and atropine (0.05 mg/kg) pre-anaesthesia, followed by isoflurane gas (1–5%) continuously delivered via mechanical ventilation to maintain the surgical stage of anaesthesia. A sheath in one carotid artery acted as a port for catheter introduction. To mitigate arrhythmia, a bolus of amiodarone (75 mg) was given prior to the intervention, in addition to lidocaine (20 mg) as needed.
The entire intervention was performed within a 1.5 T wide-bore scanner (MR450w, General Electric Healthcare Waukesha, Wisconsin, USA). Figure 1 illustrates the interventional workflow and timing of data collection throughout. A 4-channel anterior cardiac array coil was used for all imaging, and was connected to the scanner bed separately from the catheter tip tracking coils. The CMR-enabled EP hardware configuration has been described previously .
Pre-procedure scans included short- and long-axis stacks of 2D balanced steady-state free-precession (bSSFP; CINE) images to serve as anatomical roadmaps (20 cardiac phases across 1 R-R interval, TR/TE = 5.2/1.9 ms, resolution = 1.25 × 1.5 mm, slice thickness = 6 mm). Both orientations were acquired at a rate of 1 breath-hold (approximately 14 s) per slice over 8 ± 3 min in total.
A 9 F CMR-enabled catheter (Imricor Medical Systems, Burnsville, Minnesota, USA) equipped with 2 micro-receive coils and 2 electrodes (3.7-mm tip electrode and 3.5-mm inter-electrode spacing) was used for the entire intervention. Active catheter tracking using a projection sequence  (FOV = 60 cm, flip angle = 5°, TR = 14.3 ms, tracking rate = 23 fps) was implemented in RTHawk software (HeartVista Inc., Menlo Park, California, USA). The catheter tip position was updated in real time and visualized overlaid on static anatomical roadmap images (Fig. 1b) simultaneously in VURTIGO, 4D visualization software for cardiovascular interventional guidance (Sunnybrook Research Institute, Toronto, Ontario, Canada) . Bipolar intracardiac electrogram (EGM; Fig. 1c) traces from the catheter tip’s electrodes were recorded by the Advantage-MR system (Imricor Medical Systems). EP traces and VURTIGO catheter navigation visualization were displayed on monitors in the CMR control room and on CMR-compatible monitors near the scanner bed to assist the interventionalist. Once the catheter was inserted into the heart to access the LV via the aortic valve, ablations were delivered from the catheter’s distal electrode (1500 T14 generator, St Jude Medical, St. Paul, Minnesota, USA) with 17 mL/min irrigation and with the dispersive electrode placed on the animal’s back. Ablations were delivered endocardially, in accordance with the more common clinical ablation approach (details given in Table 1). Conservative ablation parameters were chosen to avoid inducing arrhythmia during ablation. Ablation was generally performed in similar locations for consistency, but targeting accuracy was not directly evaluated as no specific targeting was used for lesion placement. Vital signs including electrocardiogram (ECG), end-tidal CO2, and peripheral capillary oxygen saturation were monitored using the Expression CMR Patient Monitor (Invivo Corp., Gainesville, Florida, USA). Peripheral pulse sensing from each animal’s foot was sufficient for robust cardiac gating, and respiratory motion was frozen during imaging by pausing mechanical ventilation.
CMR imaging protocol
Imaging began immediately after catheter withdrawal (as early as 3 min after the start of ablation) to visualize the lesion near the catheter contact point. The two primary native contrast (non-contrast enhanced) imaging sequences used were alternated regularly at consistent slice locations to evaluate RF lesion temporal evolution.
T2 mapping was performed using a previously validated T2-prepared spiral sequence to detect regional T2 changes near the ablation site reflecting inflamed edematous tissue (4 TEs between 3 and 184 ms, TR = 2 R-R intervals, 10 interleaves, 3072 points, FOV = 240 mm, readout bandwidth = 125 kHz, effective resolution = 1.3 × 1.3 mm, slice thickness = 6 mm). Images at one TE were acquired per breath-hold (approximately 16 s), and complete maps were acquired in 4 ± 2 min.
An IR-SSFP sequence generated images with various inversion times at 40 cardiac phases across 2 R-R intervals . Typical acquisition parameters were: views/segment = 16, flip angle = 45°, TR/TE = 5.6/2.0 ms, readout bandwidth = 62.5 kHz, FOV = 240 mm, matrix = 192 × 160, and slice thickness = 6 mm. Although SSFP sequences yield best quality with shorter TRs (e.g., TR = 2.7–4 ms at 1.5 T [7, 23]), this TR was longer than ideal in this wide-bore system with lower gradient specifications. This was the shortest achievable TR while maintaining other sequence parameters. Local shimming was performed to mitigate any off-resonance effects, and dark-band artifacts were not observed in the regions of interest.
RF ablation lesions were hyper-enhanced in IR-SSFP images due to shorter T1 compared to the healthy myocardium [7, 24]. TIs longer than approximately 700 ms (yielding optimal contrast for lesion core visualization [7, 13, 25]) were set to occur within diastole. Image acquisition was performed at 1 slice per breath-hold (each approximately 16 s), such that 5 slices with 2 rest breaths between each slice could be acquired within 2 min. Sequence parameters were adjusted slightly to accommodate animal size and heart rate (typically 75–95 bpm).
Native-contrast imaging was repeated up to approximately 3 h post-ablation. At the end of the CMR study a bolus of Gd-DTPA (0.2 mmol/kg, Magnevist, Bayer Healthcare Pharmaceuticals, Berlin, Germany) was injected for contrast-enhanced imaging to confirm the pattern of ablation. The IR-SSFP sequence was repeated after contrast agent injection (referred to as MCLE ), acquired at 1 slice per breath-hold, at 6 ± 5 min post-injection. Four of 12 animals succumbed to arrhythmia before the end of the imaging study when MCLE images were to be acquired. Animals were not moved during the CMR study; therefore, all images were well aligned post-ablation (confirmed via anatomical landmarks in corresponding images).
Ex vivo examination
After the imaging study (within hours of ablation), animals were euthanized and the hearts explanted immediately, then preserved in a 10% formalin solution. Hearts were embedded in dental alginate gel and sliced to 4-mm blocks aligned to the imaging plane for gross examination and measurement of lesion dimensions. Tissue blocks were sliced to 4-μm sections and stained with Hematoxylin and Eosin (H&E) and Masson’s trichrome (MT) to emphasize tissue morphological changes and viability, respectively. Slides were scanned under light microscopy (Leica SCN400 F, Leica Microsystems, Wetzlar, Germany) at 20× magnification (0.5 μm resolution) to observe features of healthy tissue, the lesion, and peripheral regions around the ablation site. Edema extent was not assessed ex vivo given the challenge of determining associated borders in tissue sections. CMR images were taken as a more accurate representation of the extent of edema in vivo. Observations from histopathology were qualitative in support of the image-based observations of ablation injury patterns.
Data analysis & statistics
All analysis was performed offline in MATLAB (MathWorks, Inc., Natick, Massachusetts, USA). Images were interpolated to a 256 × 256 matrix (DICOM format) before processing. T2 maps were generated using a previously-validated 3-parameter model  (approximately 15 s per map using non-optimized code). Regions of edematous tissue exhibiting increased T2 were segmented semi-automatically by applying a threshold of 3 standard deviations (SDs) above mean T2 in healthy tissue (similarly to the approach used in a prior ablation study ). Morphological operations were applied to this mask to remove noise and other erroneous pixels (MATLAB Image Processing Toolbox), identify edema as contiguous regions of segmented pixels, and fill these regions while preserving borders. This semi-automatic approach was implemented to yield reproducible volumes despite the often-diffuse quality of edema, and could be performed within approximately 2 min per slice.
IR-SSFP images acquired at TIs longer than approximately 700 ms were selected for lesion segmentation, with matching contrast in subsequent acquisitions. The T1-derived lesion volumes were segmented semi-automatically using a threshold of 2 SD above adjacent healthy myocardial signal intensity (SI), and correlated well with an expert viewer’s manual delineation of lesion volumes (Pearson’s r = 0.95, intercept = − 0.078, slope = 1.2, p < < 0.0001; 2 SD segmented volumes larger by 0.0017 ± 0.10 mL overall; 95% limits of agreement [− 0.030, 0.033]). Segmentation required approximately 3 min per slice.
All native T1-weighted images and T2 maps included in analysis were acquired in vivo. To assess temporal changes in the T1-derived lesions and T2-derived edema, CMR-based volume measurements were normalized to those from the initial imaging time point post-ablation. We report all post-ablation times with respect to the start of ablation, as lesion formation commences at this point. Lesion volume temporal evolution was evaluated using a linear mixed model (LMM) to account for clustering by animal and baseline differences (MATLAB Statistics Toolbox). The ratio of T2-derived edema and T1-derived lesion volumes was calculated to provide a direct comparison of these volumes across individual ablations at different time intervals, assessed using analysis of variance (ANOVA). T2 evolution was assessed in three key ROIs: healthy tissue; T1-derived lesions; and edematous regions (the largest edematous ROIs for each ablation, to examine changes in a consistent tissue region). LV wall thickness and lesion transmurality (lesion depth divided by wall thickness) were also compared at different time points.
Lesion diameters measured manually from the final in vivo IR-SSFP images acquired were compared to those from morphologically matched gross tissue slices (the gold standard for lesion formation). Lesions in MCLE images were segmented manually to compare lesion volumes derived from native-contrast and contrast-enhanced images. Both comparisons were performed using Bland-Altman analysis. All measurements are reported as the mean ± SD unless indicated otherwise, and p < 0.05 was considered statistically significant.
RF lesion temporal evolution
The native T1-derived lesion and T2-derived edema were clearly visualized and reflected the characteristic teardrop shape of RF lesions observed in corresponding contrast-enhanced images, gross pathology, and histopathology (Fig. 2). Elevated T2 indicating edema surrounding the ablation site was evident in the first T2 maps acquired as early as 3 min after the start of ablation. T2 maps were analyzed across 13 ablations visualized in 1–3 adjacent imaging slices in 11 animals. Each T2-derived edematous region was visualized at 5 ± 3 (median 4) time points following ablation. These regions tended to initially appear more focal and localized to the lesion border, assuming a more diffuse appearance by later time points (Fig. 3).
The volume of T2-derived edema was 0.77 ± 0.55 mL (n = 13) at the baseline measurements from the earliest T2 maps acquired (median 10 min post-ablation). The development of edema is shown in Fig. 4, illustrating the overall increase in volume compared to baseline. Normalized edema volume increased significantly beyond baseline (LMM slope = 0.003 ± 0.001, mean ± SE, p = 0.009, intercept = 1.20 ± 0.01, p < 0.0001). Using this model, T2-derived edema expanded to an estimated 1.7 ± 0.2 (mean ± SE) times the baseline volume by 180 min post-ablation. Although the overall trend indicated increasing volume, in 6/13 of T2-derived edematous regions the volume increased then later dropped (while still remaining 1.3 ± 0.6 times larger than baseline). The maximum T2-derived edema volume estimated across all lesions was 1.5 ± 1.0 mL. T2-derived edema did not consistently develop concentrically around the catheter contact point.
The TIs used for T1-derived lesion visualization using IR-SSFP were TI = 831 ± 150 ms. These lesion volumes were measured in 9 lesions in 9 animals, each visualized in 1–3 adjacent imaging slices at 5 ± 2 (median 4) time points post-ablation (summarized in Fig. 5). The mean volume from the baseline IR-SSFP images was 0.48 ± 0.23 mL (n = 9, median 20 min post-ablation). Cumulatively, a small trend towards decreasing normalized T1-derived lesion volume was detected up to 180 min post-ablation (LMM slope = − 0.0006 ± 0.0003, mean ± SE, p = 0.09, intercept = 1.01 ± 0.04, p < 0.0001). Using this model, T1-derived lesion volumes reached an estimated 0.9 ± 0.1 (mean ± SE) times the baseline volume by 180 min post-ablation.
The T2-derived edema consistently encompassed T1-derived lesions at individual ablation sites (Fig. 6), further supporting the distinctive tissue regions identified with each imaging contrast. The edema-lesion volume ratio was initially 2.1 ± 1.0 (0–25 min post-ablation, 95% CI [1.1, 3.0]). This ratio increased to 3.6 ± 1.1 (60–80 min; 95% CI [2.5, 4.8]), and reached 4.4 ± 3.2 (95% CI [1.1, 7.8]) by the final time interval at 80–185 min post-ablation. Although differences between intervals were not statistically significant (p = 0.2, ANOVA), volume ratios were greater than 1 within each interval (p < 0.05). The RF energy delivered did not correlate to the maximum measured T1-derived lesion or T2-derived edema volumes (r = − 0.2, 0.4; p = 0.7, 0.2, Spearman’s rho test); however, true energy deposition likely varied with catheter contact force, orientation, and blood flow.
T2 was evaluated in 3 critical regions: the T1-derived lesion; maximum T2-derived edematous region; and healthy tissue (Fig. 7). T2 within the T1-derived lesion ROIs increased between 0 and 25 min (53 ± 10 ms) and 60–80 min post-ablation (55 ± 8 ms, p = 0.04), then returned to the initial T2 at 80–185 min (53 ± 7 ms, p = 0.08). T2 within the largest edematous regions was significantly higher than in the T1-derived lesions overall (p < 0.001) and initially increased between 0 and 25 min (53 ± 15 ms) and 60–80 min post-ablation (58 ± 13 ms, p < 0.001), then decreased at 80–185 min, still remaining significantly higher than baseline (56 ± 15 ms, p < 0.001). T2 within both T1-derived lesions and largest edematous regions were significantly higher than in healthy tissue (39 ± 5 ms, p < 0.001), and T2 in healthy tissue did not change significantly during corresponding time intervals (p = 0.2, ANOVA).
The lesions were not generally transmural. T2-derived edema was 76 ± 18% transmural at baseline (15 ± 10 min post-ablation), reaching 79 ± 22% at final measurement (120 ± 67 min; p = 0.9). Similarly, transmurality of the T1-derived lesion did not change significantly from initial measurements of 52 ± 10% (21 ± 7 min post-ablation) to 50 ± 11% (162 ± 63 min; p = 0.3). LV wall swelling occurred at the ablation sites between pre-ablation (8 ± 2 mm) and baseline post-ablation measurements (11 ± 2 mm; p = 0.01), but further swelling through to the end of the studies was not substantial (12 ± 2 mm; p = 0.2). Wall thickness at remote sites remained consistent with pre-ablation measurements of 9 ± 2 mm, reaching 8 ± 2 mm initially post-ablation (p = 0.5) and 9 ± 2 mm (p = 0.9) at the end of the studies.
Final pattern of ablation
MCLE validation imaging reflected the lesion geometry consistently observed in corresponding intra-procedural images (Fig. 2). The MCLE-derived lesion volumes were 0.75 ± 0.39 mL (n = 9), significantly larger than those measured in the corresponding final IR-SSFP acquisitions (n = 6; bias = 0.41 ± 0.31 mL, p = 0.03; 95% limits of agreement [0.08, 0.74]), and smaller than the corresponding final T2-derived edema volumes (n = 9; bias = 0.45 ± 0.79 mL, p = 0.2; 95% limits of agreement [− 0.16, 1.1]).
Lesion diameters observed in gross pathology slices and in the final IR-SSFP acquisition were well correlated (n = 9 lesions, Pearson’s r = 0.87, slope = 0.9, p = 0.003), with a relatively small bias of 0.4 ± 1.0 mm (gross lesion diameter = 9.7 ± 2.0 mm, T1-derived lesion diameter = 10.1 ± 2.0 mm; 95% limits of agreement [− 0.41, 1.2]).
Histological characteristics of RF lesions
RF lesions exhibited the characteristic pale pink core of thermal injury bordered by a dark rim in gross pathology (Fig. 2d), also reflected in histological sections (Fig. 2e-f). Morphologic changes were emphasized using H&E (cytoplasm and extracellular matrix are stained pink, nuclei stained purple), and viability emphasized using MT (ischemic or necrotic tissue is stained purple, healthy viable myocytes red, and connective tissue blue).
The lesion core in H&E sections (Fig. 8a-b) exhibited disrupted cellular architecture consistent with thermal coagulation [28, 29] and the purple colour of the corresponding MT section (Fig. 8f) suggested non-viability of this tissue. The surrounding rim (Fig. 8c) was also distinguished by altered cellular architecture, and contained extravasated red blood cells, evidence of hemorrhage. The purple-to-red colour gradient outward from the lesion rim (Fig. 8g) suggested an increasing proportion of viable myocytes with distance from the catheter contact point. Interstitial space was increased conspicuously throughout the lesion (Fig. 8b-c) and extending beyond the lesion rim (compared to healthy tissue; Fig. 8d), consistent with observations of broad T2-derived edematous regions.
Native-contrast CMR was used to construct a comprehensive description of RF lesion composition and temporal evolution in the LV. Intra-procedural T1- and T2-based imaging performed concurrently each reflect a distinctive region of tissue injury, evolving separately throughout the time period of an intervention. Acute T1-derived lesions (thought to represent the permanent lesion) remained at a relatively stable size whereas the broader T2-derived edema (likely transient injury) was dynamic, tended to expand over time, and consistently extended beyond T1-derived lesions.
Lesion visualization immediately post-ablation
By performing ablation under CMR guidance, intra-procedural visualization of the T1-derived lesions and T2-derived edema was achieved within minutes after ablation. T1 contrast after ablation is believed to originate from oxidation of ferrous (Fe2+) to ferric (Fe3+) iron associated with the transformation of hemoglobin to paramagnetic methemoglobin occurring at 55–65 °C [7, 30], and a recent study using real-time CMR thermometry during ablation showed that native T1 contrast reflects tissue which received a lethal thermal dose .
The acute lesions appeared desiccated in the central core, from T2 maps showing a zone of shorter T2 than the surrounding tissue (Fig. 3). Ablation-induced injury may be explained by comparison to general descriptions of thermal injury and edema pathogenesis [28, 31]. A sub-lethal thermal dose in tissue surrounding the lesion core is likely the cause of edema formation there. Microvessel permeability may be initially increased around the lesion core due to the release of biochemical factors such as histamine into injured tissue, leading to accumulation of fluid (edema).
Acute lesion evolution minutes to hours post-ablation
Native-contrast CMR facilitated repeated measurement of both the T1-derived lesion and the T2-derived edema to assess their acute evolution throughout the interventional procedure. The evolution of native CMR contrast is of particular importance since lesion measurements can be repeated easily at different points in the procedure, while gadolinium-enhanced studies are limited to a single injection of contrast.
The acute T1-derived lesions appeared largely stable in volume, corroborating and extending existing studies that showed T1-derived lesions visualized acutely had a consistent extent which corresponded well to lesions visualized hours, weeks, and months post-ablation [5, 6, 13, 17]. By combining the results from this study with the existing data on T1-derived lesions, we suggest that the spatial extent of these lesions is established immediately after ablation, with damaged tissue later replaced by fibrotic scar. The subtle trend seen in this study toward decreasing T1-derived lesion size with time could be explained by passive edema diffusion toward the initially desiccated centre of the lesion, driven by reduced interstitial pressure there. Interstitial pressure changes in burned tissue are generally driven by the release of particles which drive osmotic pressure, drawing fluid into the interstitium . Increasing water content would elevate local T1, opposing the reduction of T1 due to iron transformation.
The interstitium in thermally injured tissue tends to become more compliant, perpetuating fluid accumulation, with damage to matrix molecules such as collagen and hyaluronic acid. Increased permeability of injured microvessels leads to continued fluid leakage into the interstitial space . Microvascular injury has previously been detected up to 6 mm beyond the pathological RF lesion core [29, 32], and the resulting ischemic effects likely contributed to the continuing increase in T2-derived edematous volumes and the increasing T2 within these regions up to approximately 60 min post-ablation (Fig. 7). The edema likely reached a stable extent once opposing interstitial and vascular pressures equalized, then later a more homogeneous fluid distribution with passive diffusion, consistent with the decreasing then stable T2 after the initial peak seen in edematous tissue towards the end of the CMR studies.
Non-concentric edema development of the T2-derived edema, and stable transmurality, could reflect preferential fluid spread occurring along cleavage planes between myocardial sheets, in line with LV muscle fibres that primarily run circumferentially. Beyond 120 min post-ablation, edematous regions appeared to have more diffuse borders (e.g., Fig. 3). Variable volume measurements at later time points could result from cases in which T2 elevation at the ablation site dropped below the 3 SD threshold; intensity-based segmentation is less robust to this diffuse pattern. Despite volume variability, the T2-derived edema volumes were still larger than immediate post-ablation dimensions (Fig. 4) and the corresponding T1-derived lesion volumes (Fig. 6).
Prior studies using T2-weighted CMR showed little change in edematous area beyond 3 min post-ablation  and that the maximum lesion SI occurred at approximately 12 min post-ablation in the right ventricle (RV) . In the current study, elevated T2 in the maximum T2-derived edema ROIs reached a maximum at approximately 60 min post-ablation, possibly arising from a greater capacity for edema development in the LV (compared to thin-walled RV and atria), or the approach used to delineate the edema.
Study limitations & future work
The earliest T2 map acquisition in this study was 3 min after the start of RF energy delivery, and delays typically included catheter withdrawal and waiting for any ablation-induced arrhythmia to settle. T1-weighted image acquisition was performed at 13 min at the earliest; therefore any earlier changes among T1-derived lesions were missed. Nevertheless, these lesions appeared consistently smaller than the T2-derived edema, and were stable or shrinking instead of expanding for time points beyond 13 min post-ablation.
The 2D imaging sequences used in this study were appropriate to maintain relatively short scan times for concurrent imaging of the T1-derived lesion and T2-derived edema. We anticipate the use of 3D MCLE, a novel sequence using compressed sensing to acquire images of isotropic resolution in a single breath-hold developed for infarct characterization, in future procedures for rapid lesion assessment . Moving towards a future clinical workflow, high-resolution imaging would likely be necessary for pre-ablation targeting--the importance of fine detail to describe complex re-entry circuits is well established in both interventional EP and imaging communities [34, 35]--as well as benefiting lesion assessment, as seen recently .
The pro-arrhythmic potential of tissue containing mixed viable and non-viable myocytes and potentially microvascular injury near the periphery of RF lesions was not investigated. However, prior work showed chronic lesions resembling well-circumscribed scar , a pattern associated with lower arrhythmogenic potential than heterogeneous scar. The potential for pro-arrhythmia is likely greater in healthy tissue, as was the case in the healthy swine model used for this study. Practically, in the presence of infarcted tissue (as is typical for VT in structural heart disease), lesion pro-arrhythmia may have a small effect since areas ablated are likely to be already ischemic and surrounded by chronic fibrotic tissue. Investigating the properties of ablations on or adjacent to infarcted tissue would build on promising recent results  and is an important direction for future studies. In ongoing CMR-based lesion assessment studies, we are also investigating the effect of ablation on the functional properties of local myocardial tissue , extending previous studies of infarcted tissue .
CMR-based intra-procedural lesion assessment
Native-contrast CMR imaging reflects ablation-induced tissue changes and facilitates consistent lesion assessment throughout MR-guided ablation. Considering the complexity of the arrhythmogenic substrate, lesion assessment should be repeatable when convenient for the interventionalist, supporting an iterative treatment approach.
Patients undergoing ablation for structural heart disease often exhibit multiple arrhythmia morphologies, suggestive of multiple re-entry circuits. CMR-based lesion assessment could be performed after delivering several ablations towards blocking a target re-entry circuit. Intermittent imaging exams and analysis would require 10–20 min to perform using the techniques employed in the current study (with segmentation not fully optimized for speed).
One proposed workflow for future clinical cases could involve identifying re-entry circuits in prior CMR-derived maps of scar and EP mapping, then delivering several ablations at a re-entry circuit while relying on EP data. T1-weighted CMR could be used to visualize a series of lesions created at the re-entry circuit. We suggest that ablation while incorporating this feedback should aim to produce lesions which appear continuous using T1-weighted imaging. Based on these T1-weighted images, remaining gaps identified at the re-entry circuit would clearly indicate that further ablation is needed, potentially constituting an additional procedural endpoint. T2 mapping could be used to interpret possibly discordant EP measures of procedural success, but should not be used alone for lesion assessment as the T2-derived edema likely reflects transient injury not contributing to long-term procedural success.
RF lesions could be visualized with the catheter held in place if lesions extend beyond the minimal artifact at the catheter tip, and others have proposed imaging and hardware based solutions for interventions [38, 39]. Catheter withdrawal before imaging (as in this study) may introduce challenges by causing the interventionalist to lose the catheter’s position at the ablation site such that it would need to be directed back precisely for further ablations.
Conscious patients under light sedation may be able to perform breath-holds during intra-procedural CMR imaging, but for those unable, existing respiratory-navigated or real-time sequences could be adapted for lesion imaging . Rapid intra-procedural registration of lesion images to prior scar maps would be needed to provide useful feedback to the interventionalist.
Finally, contrast-enhanced imaging should be reserved to confirm the pattern of ablation at the true completion of the procedure. An alternative sequence recently proposed visualizes both scar and RF lesions with native T1 contrast . For either approach, high-resolution coverage of an extended region of the heart, after ablating multiple re-entry circuits, would require a respiratory-navigated acquisition of about 10–20 min. The timing of this acquisition after contrast-agent injection should also be chosen carefully given known contrast agent kinetics.
Comprehensive understanding of acute lesion composition and evolution could be applied to both CMR-guided and non-CMR-guided ablation cases. In general, operators should be cognizant of the dynamic nature of edema, which could cause the acute appearance of a successful ablation procedure. The evolving, broad extent of T2-derived edema may mask the smaller T1-derived lesions (specifically indicated by the ratio of these volumes) which are thought to reflect true procedural success. T1-based imaging could potentially constitute an additional procedural endpoint. CMR-based lesion assessment could make a critical difference in scenarios where re-entry circuits appear to be blocked, for instance with more transmural T2-derived edema injuring deep epicardial re-entry circuits not reached by the T1-derived lesion. These findings may be used to help explain mechanisms of arrhythmia recurrence after acutely successful appearing procedures.
Native T1- and T2-based imaging performed concurrently throughout CMR-guided ablations demonstrated that acute T1-derived lesions remained at a stable size while the T2-derived edema was dynamic, expanded over time, and consistently extended beyond the T1-derived lesion. T2 quantification provides an unambiguous measure of edema development and severity near ablation sites. Lesion evolution is important for comparing CMR lesion visualization and EP-based endpoints, data which may be acquired at disparate time points. Integrated information on the T1-derived lesion and T2-derived edema provides a description of acute lesion composition that could be most useful during ablation procedures.
Analysis of variance
Beats per minute
balanced steady state free precession
cardiovascular magnetic resonance
dice similarity coefficient
field of view
frames per second
Hematoxylin and Eosin
Inversion-recovery balanced steady-state free precession
linear mixed model
left ventricle/left ventricular
multi-contrast late enhancement
region of interest
right ventricle/right ventricular
Palaniswamy C, Kolte D, Harikrishnan P, Khera S, Aronow WS, Mujib M, et al. Catheter ablation of postinfarction ventricular tachycardia: ten-year trends in utilization, in-hospital complications, and in-hospital mortality in the United States. Heart Rhythm. 2014;11:2056–63.
Mallidi J, Nadkarni GN, Berger RD, Calkins H, Nazarian S. Meta-analysis of catheter ablation as an adjunct to medical therapy for treatment of ventricular tachycardia in patients with structural heart disease. Heart Rhythm. 2011;8:503–10.
D’Avila A, Singh SM. Ventricular tachycardia ablation: are we winning the battle but losing the war? J Am Coll Cardiol. 2016;67:684–6.
Lardo AC, McVeigh ER, Jumrussirikul P, Berger RD, Calkins H, Lima J, et al. Visualization and temporal/spatial characterization of cardiac radiofrequency ablation lesions using magnetic resonance imaging. Circulation. 2000;102:698–705.
Dickfeld T, Kato R, Zviman M, Nazarian S, Dong J, Ashikaga H, et al. Characterization of acute and subacute radiofrequency ablation lesions with non-enhanced magnetic resonance imaging. Heart Rhythm. 2007;4:208–14.
Kholmovski EG, Ranjan R, Vijayakumar S, Silvernagel JM, Marrouche NF. Acute Assessment of Radiofrequency Ablation Cardiac Lesions by Non-Contrast MRI. In: Proc. ISMRM, 22nd Annual Meeting. Milan; 2014.
Celik H, Ramanan V, Barry J, Ghate S, Leber V, Oduneye S, et al. Intrinsic contrast for characterization of acute radiofrequency ablation lesions. Circ Arrhythmia Electrophysiol. 2014;7:718–27.
Arujuna A, Karim R, Caulfield D, Knowles B, Rhode K, Schaeffter T, et al. Acute pulmonary vein isolation is achieved by a combination of reversible and irreversible atrial injury after catheter ablation. Circ Arrhythmia Electrophysiol. 2012;5:691–700.
Yamada T, Murakami Y, Okada T, Okamoto M, Shimizu T, Toyama J, et al. Incidence, location, and cause of recovery of electrical connections between the pulmonary veins and the left atrium after pulmonary vein isolation. Europace. 2006;8:182–8.
Knowles BR, Caulfield D, Cooklin M, Rinaldi CA, Gill J, Bostock J, et al. 3-D visualization of acute RF ablation lesions using MRI for the simultaneous determination of the patterns of necrosis and edema. IEEE TBME. 2010;57:1467–75.
Ranjan R, Kato R, Zviman MM, Dickfeld TM, Roguin A, Berger RD, et al. Gaps in the ablation line as a potential cause of recovery from electrical isolation and their visualization using MRI. Circ Arrhythm Electrophysiol. 2011;4:279–86.
Piers SRD, Tao Q, De Riva Silva M, Siebelink HM, Schalij MJ, Van Der Geest RJ, et al. CMR-based identification of critical isthmus sites of ischemic and nonischemic ventricular tachycardia. JACC Cardiovasc Imaging. 2014;7:774–84.
Guttman MA, Tao S, Fink S, Kolandaivelu A, Halperin HR, Herzka DA. Non-contrast-enhanced T1 -weighted MRI of myocardial radiofrequency ablation lesions. Magn Reson Med. 2017;0:1–11.
Rassa AC, Kholmovski E, Suksaranjit P, Wilson BD, Akoum N, Marrouche N, et al. Dynamic T2 signal changes on MRI after radiofrequency ablation injury to the atrial myocardium. J Clin Trials Cardiol. 2015;2:1–5.
Vergara GR, Vijayakumar S, Kholmovski EG, Blauer JJE, Guttman MA, Gloschat C, et al. Real-time magnetic resonance imaging-guided radiofrequency atrial ablation and visualization of lesion formation at 3 tesla. Heart Rhythm. 2011;8:295–303.
Nordbeck P, Hiller KH, Fidler F, Warmuth M, Burkard N, Nahrendorf M, et al. Feasibility of contrast-enhanced and nonenhanced MRI for intraprocedural and postprocedural lesion visualization in interventional electrophysiology : animal studies and early delineation of isthmus ablation lesions in patients with typical atrial flutter. Circ Cardiovasc Imaging. 2011;4:282–94.
Toupin S, Bour P, Lepetit-Coiffé M, Ozenne V, Denis de Senneville B, Schneider R, et al. Feasibility of real-time MR thermal dose mapping for predicting radiofrequency ablation outcome in the myocardium in vivo. J Cardiovasc Magn Reson. 2017;19:14.
Higgins CB, Herfkens R, Lipton MJ, Sievers R, Sheldon P, Kaufman L, et al. Nuclear magnetic resonance imaging of acute myocardial infarction in dogs: alterations in magnetic relaxation times. Am J Cardiol. 1983;52:184–8.
Dickfeld T, Kato R, Zviman M, Lai S, Meininger G, Lardo AC, et al. Characterization of radiofrequency ablation lesions with gadolinium-enhanced cardiovascular magnetic resonance imaging. J Am Coll Cardiol. 2006;47:370–8.
Oduneye SO, Biswas L, Ghate S, Ramanan V, Barry J, Laish-FarKash A, et al. The feasibility of endocardial propagation mapping using magnetic resonance guidance in a swine model, and comparison with standard electroanatomic mapping. IEEE TMI. 2012;31:977–83.
Dumoulin CL, Souza SP, Darrow RD. Real-time position monitoring of invasive devices using magnetic resonance. Magn Reson Med. 1993;29:411–5.
Radau PE, Pintilie S, Flor R, Biswas L, Oduneye SO, Ramanan V, et al. VURTIGO : visualization platform for real-time, MRI-guided cardiac Electroanatomic mapping. Lect Notes Comput Sci. 2012;7085:244–53.
Detsky JS, Stainsby JA, Vijayaraghavan R, Graham JJ, Dick AJ, Wright GA. Inversion-recovery-prepared SSFP for cardiac-phase-resolved delayed-enhancement MRI. Magn Reson Med. 2007;58:365–72.
Herzka DA, Tao S, Fink S, Kolandaivelu A, Guttman MA, Halperin H. Assessment of RF ablation lesions with T1 mapping. In: 20th annual SCMR scientific sessions; 2017. p. 203.
Krahn P, Ramanan V, Biswas L, Yak N, Anderson K, Barry J, et al. Intrinsic MR visualization of RF lesions after RF ablation for ventricular arrhythmia. In: Proc. ISMRM, 24th Annual Meeting. Singapore; 2016.
Ghugre NR, Enriquez CM, Gonzalez I, Nelson MD, Coates TD, Wood JC. MRI detects myocardial iron in the human heart. Magn Reson Med. 2006;56:681–6.
Chubb H, Harrison JL, Weiss S, Krueger S, Koken P, Bloch LØ, et al. Development, pre-clinical validation, and clinical translation of a cardiac magnetic resonance-electrophysiology system with active catheter tracking for ablation of cardiac arrhythmia. JACC Clin Electrophysiol. 2017;3:89–103.
Thomsen S. Mapping of thermal injury in biologic tissues using quantitative pathologic techniques. In: proceedings of thermal treatment of tissue with image guidance. San Jose, California: SPIE; 1999. p. 82–95.
Nath S, Redick JA, Whayne JG, Haines DE. Ultrastructural observations in the myocardium beyond the region of acute coagulation necrosis following radiofrequency catheter ablation. J Cardiovasc Electrophysiol. 1994;5:838–45.
Farahani K, Saxton RE, Yoon H-C, De Salles AAF, Black KL, Lufkin RB. MRI of thermally denatured blood: methemoglobin formation and relaxation effects. Magn Reson Imaging. 1999;17:1489–94.
Demling RH. The burn edema process: current concepts. J Burn Care Rehabil. 2005;26:207–27.
Nath S, Whayne JG, Kaul S, Goodman NC, Jayaweera a R, Haines DE. Effects of radiofrequency catheter ablation on regional myocardial blood flow. Possible mechanism for late electrophysiological outcome. Circulation. 1994;89:2667–72.
Zhang L, Athavale P, Pop M, Wright GA. Multicontrast reconstruction using compressed sensing with low rank and spatially varying edge-preserving constraints for high-resolution MR characterization of myocardial infarction. Magn Reson Med. 2017;78:598–610.
Pop M, Ramanan V, Yang F, Zhang L, Newbigging S, Ghugre NR, et al. High-resolution 3-D T1*-mapping and quantitative image analysis of GRAY ZONE in chronic fibrosis. IEEE TBME. 2014;61:2930–8.
Ashikaga H, Sasano T, Dong J, Zviman MM, Evers R, Hopenfeld B, et al. Magnetic resonance-based anatomical analysis of scar-related ventricular tachycardia: implications for catheter ablation. Circ Res. 2007;101:939–47.
Krahn P, Ramanan V, Biswas L, Yak N, Anderson K, Barry J, et al. MRI-based myocardial ablation lesion extent relates to area of voltage reduction in MR-guided electroanatomical voltage maps. In: Proc. ISMRM, 25th Annual Meeting. Honolulu; 2017.
Oduneye SO, Pop M, Shurrab M, Biswas L, Ramanan V, Barry J, et al. Distribution of abnormal potentials in chronic myocardial infarction using a real time magnetic resonance guided electrophysiology system. J Cardiovasc Magn Reson. 2015;17:1–9.
Daniel BL, Butts K. The use of view angle tilting to reduce distortions in magnetic resonance imaging of cryosurgery. Magn Reson Imaging. 2000;18:281–6.
Dominguez-Viqueira W, Karimi H, Lam WW, Cunningham CH. A controllable susceptibility marker for passive device tracking. Magn Reson Med. 2014;72:269–75.
Detsky JS, Graham JJ, Vijayaraghavan R, Biswas L, Stainsby JA, Guttman MA, et al. Free-breathing, nongated real-time enhancement MRI of myocardial infarcts: a comparison with conventional delayed enhancement. J Magn Reson Imaging. 2008;28:621–5.
Lu Y, Radau P, Connelly K, Dick A, Wright GA. Segmentation of left ventricle in cardiac Cine MRI : an automatic image-driven method. Lect Notes Comput Sci. 2009;5528:339–47.
The authors sincerely thank Dr. Nilesh Ghugre for insights on image analysis and relaxometry, Sebastian Ferguson for assistance in CMR-guided ablation studies, Dr. Alex Kiss for helpful discussions on statistical methods, and Adebayo Adeeko and Taha Rashed for histopathology processing.
Funding provided by Ontario Graduate Scholarships, Imricor Medical Systems, the Federal Development Agency of Canada, and the Canadian Institute of Health Research (grant MOP-93531) are gratefully acknowledged.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Animal use was approved by the animal care committee at Sunnybrook Health Sciences Centre. All procedures followed institutional guidelines according to the approved Animal Use Protocol #16-235.
Consent for publication
This work was supported in part by Imricor Medical Systems and the authors report a pending patent related to RF lesion characterization using T1-weighted CMR imaging. The authors declare no other competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Krahn, P.R.P., Singh, S.M., Ramanan, V. et al. Cardiovascular magnetic resonance guided ablation and intra-procedural visualization of evolving radiofrequency lesions in the left ventricle. J Cardiovasc Magn Reson 20, 20 (2018). https://doi.org/10.1186/s12968-018-0437-z
- Image-guided intervention
- Catheter ablation
- Tissue characterization