Evaluation of Motion Compensation Methods for Noninvasive Cardiac Radioablation of Ventricular Tachycardia

Cardiorespiratory motion characteristics and their dosimetric impact on cardiac stereotactic body radiotherapy

BACKGROUND

Cardiac stereotactic body radiotherapy (CSBRT) is an emerging and promising noninvasive technique for treating refractory arrhythmias utilizing highly precise, single or limited-fraction high-dose irradiations. This method promises to revolutionize the treatment of cardiac conditions by delivering targeted therapy with minimal exposure to surrounding healthy tissues. However, the dynamic nature of cardiorespiratory motion poses significant challenges to the precise delivery of dose in CSBRT, introducing potential variabilities that can impact treatment efficacy. The complexities of the influence of cardiorespiratory motion on dose distribution are compounded by interplay and blurring effects, introducing additional layers of dose uncertainty. These effects, critical to the understanding and improvement of the accuracy of CSBRT, remain unexplored, presenting a gap in current clinical literature.

PURPOSE

To investigate the cardiorespiratory motion characteristics in arrhythmia patients and the dosimetric impact of interplay and blurring effects induced by cardiorespiratory motion on CSBRT plan quality.

METHODS

The position and volume variations in the substrate target and cardiac substructures were evaluated in 12 arrhythmia patients using displacement maximum (DMX) and volume metrics. Moreover, a four-dimensional (4D) dose reconstruction approach was employed to examine the dose uncertainty of the cardiorespiratory motion.

RESULTS

Cardiac pulsation induced lower DMX than respiratory motion but increased the coefficient of variation and relative range in cardiac substructure volumes. The mean DMX of the substrate target was 0.52 cm (range: 0.26-0.80 cm) for cardiac pulsation and 0.82 cm (range: 0.32-2.05 cm) for respiratory motion. The mean DMX of the cardiac structure ranged from 0.15 to 1.56 cm during cardiac pulsation and from 0.35 to 1.89 cm during respiratory motion. Cardiac pulsation resulted in an average deviation of -0.73% (range: -4.01%-4.47%) in V25 between the 3D and 4D doses. The mean deviations in the homogeneity index (HI) and gradient index (GI) were 1.70% (range: -3.10%-4.36%) and 0.03 (range: -0.14-0.11), respectively. For cardiac substructures, the deviations in D50 due to cardiac pulsation ranged from -1.88% to 1.44%, whereas the deviations in Dmax ranged from -2.96% to 0.88% of the prescription dose. By contrast, the respiratory motion led to a mean deviation of -1.50% (range: -10.73%-4.23%) in V25. The mean deviations in HI and GI due to respiratory motion were 4.43% (range: -3.89%-13.98%) and 0.18 (range: -0.01-0.47) (p < 0.05), respectively. Furthermore, the deviations in D50 and Dmax in cardiac substructures for the respiratory motion ranged from -0.28% to 4.24% and -4.12% to 1.16%, respectively.

CONCLUSIONS

Cardiorespiratory motion characteristics vary among patients, with the respiratory motion being more significant. The intricate cardiorespiratory motion characteristics and CSBRT plan complexity can induce substantial dose uncertainty. Therefore, assessing individual motion characteristics and 4D dose reconstruction techniques is critical for implementing CSBRT without compromising efficacy and safety.

Impact of abdominal compression on heart and stomach motion for stereotactic arrhythmia radioablation

PURPOSE

To evaluate the effectiveness of abdominal compression (AC) as a respiratory motion management method for the heart and stomach during stereotactic arrhythmia radioablation (STAR).

METHODS

4D computed tomography (4DCT) scans of patients imaged with AC or without AC (free-breathing: FB) were obtained from ventricular-tachycardia (VT) (n = 3), lung cancer (n = 18), and liver cancer (n = 18) patients. Patients treated for VT were imaged both FB and with AC. Lung and liver patients were imaged once with FB or with AC, respectively. The heart, left ventricle (LV), LV components (LVCs), and stomach were contoured on each phase of the 4DCTs. Centre of mass (COM) translations in the left/right (LR), ant/post (AP), and sup/inf (SI) directions were measured for each structure. Minimum distances between LVCs and the stomach over the respiratory cycle were also measured on each 4DCT phase. Mann-Whitney U-tests were performed between AC and FB datasets with a significance of α = 0.05.

RESULTS

No statistical difference (all p values were >0.05) was found in COM translations between FB and AC patient datasets for all contoured cardiac structures. A reduction in COM translation with AC relative to FB was patient, direction, and structure specific for the three VT patients. A significant decrease in the AP range of motion of the stomach was observed under AC compared to FB. No statistical difference was found between minimum distances to the stomach and LVCs between FB and AC.

CONCLUSIONS

AC was not a consistent motion management method for STAR, nor does not uniformly affect the separation distance between LVCs and the stomach. If AC is employed in future STAR protocols, the motion of the target volume and its relative distance to the stomach should be compared on two 4DCTs: one while the patient is FB and one under AC.

Radiation therapy for ventricular arrhythmias

Ventricular arrhythmias (VA) can be life-threatening arrhythmias that result in significant morbidity and mortality. Catheter ablation (CA) is an invasive treatment modality that can be effective in the treatment of VA where medications fail. Recurrence occurs commonly following CA due to an inability to deliver lesions of adequate depth to cauterise the electrical circuits that drive VA or reach areas of scar responsible for VA. Stereotactic body radiotherapy is a non-invasive treatment modality that allows volumetric delivery of energy to treat circuits that cannot be reached by CA. It overcomes the weaknesses of CA and has been successfully utilised in small clinical trials to treat refractory VA. This article summarises the current evidence for this novel treatment modality and the steps that will be required to bring it to the forefront of VA treatment.

Analysis of left ventricle regional myocardial motion for cardiac radioablation: Left ventricular motion analysis

PURPOSE

Left ventricle (LV) regional myocardial displacement due to cardiac motion was assessed using cardiovascular magnetic resonance (CMR) cine images to establish region-specific margins for cardiac radioablation treatments.

METHODS

CMR breath-hold cine images and LV myocardial tissue contour points were analyzed for 200 subjects, including controls (n = 50) and heart failure (HF) patients with preserved ejection fraction (HFpEF, n = 50), mid-range ejection fraction (HFmrEF, n = 50), and reduced ejection fraction (HFrEF, n = 50). Contour points were divided into segments according to the 17-segment model. For each patient, contour point displacements were determined for the long-axis (all 17 segments) and short-axis (segments 1-12) directions. Mean overall, tangential (longitudinal or circumferential), and normal (radial) displacements were calculated for the 17 segments and for each segment level.

RESULTS

The greatest overall motion was observed in the control group-long axis: 4.5 ± 1.2 mm (segment 13 [apical anterior] epicardium) to 13.8 ± 3.0 mm (segment 6 [basal anterolateral] endocardium), short axis: 4.3 ± 0.8 mm (segment 9 [mid inferoseptal] epicardium) to 11.5 ± 2.3 mm (segment 1 [basal anterior] endocardium). HF patients exhibited lesser motion, with the smallest overall displacements observed in the HFrEF group-long axis: 4.3 ± 1.7 mm (segment 13 [apical anterior] epicardium) to 10.6 ± 3.4 mm (segment 6 [basal anterolateral] endocardium), short axis: 3.9 ± 1.3 mm (segment 8 [mid anteroseptal] epicardium) to 7.4 ± 2.8 mm (segment 1 [basal anterior] endocardium).

CONCLUSIONS

This analysis provides an estimate of epicardial and endocardial displacement for the 17 segments of the LV for patients with normal and impaired LV function. This reference data can be used to establish treatment planning margin guidelines for cardiac radioablation. Smaller margins may be used for patients with higher degree of impaired heart function, depending on the LV segment.

Multimodal fusion workflow for target delineation in cardiac radioablation of ventricular tachycardia

BACKGROUND

Cardiac radioablation (CR) is an innovative treatment to ablate cardiac arrythmia sources by radiation therapy. CR target delineation is a challenging task requiring the exploitation of highly different imaging modalities, including cardiac electro-anatomical mapping (EAM).

PURPOSE

In this work, a data integration process is proposed to alleviate the tediousness of CR target delineation by generating a fused representation of the heart, including all the information of interest resulting from the analysis and registration of electro-anatomical data, PET scan and planning computed tomography (CT) scan. The proposed process was evaluated by cardiologists during delineation trials.

METHODS

The data processing pipeline was composed of the following steps. The cardiac structures of interest were segmented from cardiac CT scans using a deep learning method. The EAM data was registered to the cardiac CT scan using a point cloud based registration method. The PET scan was registered using rigid image registration. The EAM and PET information, as well as the myocardium thickness, were projected on the surface of the 3D mesh of the left ventricle. The target was identified by delineating a path on this surface that was further projected to the thickness of the myocardium to create the target volume. This process was evaluated by comparison with a standard slice-by-slice delineation with mental EAM registration. Four cardiologists delineated targets for three patients using both methods. The variability of target volumes, and the ease of use of the proposed method, were evaluated.

RESULTS

All cardiologists reported being more confident and efficient using the proposed method. The inter-clinician variability in delineated target volume was systematically lower with the proposed method (average dice score of 0.62 vs. 0.32 with a classical method). Delineation times were also improved.

CONCLUSIONS

A data integration process was proposed and evaluated to fuse images of interest for CR target delineation. It effectively reduces the tediousness of CR target delineation, while improving inter-clinician agreement on target volumes. This study is still to be confirmed by including more clinicians and patient data to the experiments.

Stereotactic Arrhythmia Radioablation (STAR): Assessment of cardiac and respiratory heart motion in ventricular tachycardia patients - A STOPSTORM.eu consortium review

AIM

To identify the optimal STereotactic Arrhythmia Radioablation (STAR) strategy for individual patients, cardiorespiratory motion of the target volume in combination with different treatment methodologies needs to be evaluated. However, an authoritative overview of the amount of cardiorespiratory motion in ventricular tachycardia (VT) patients is missing.

METHODS

In this STOPSTORM consortium study, we performed a literature review to gain insight into cardiorespiratory motion of target volumes for STAR. Motion data and target volumes were extracted and summarized.

RESULTS

Out of the 232 studies screened, 56 provided data on cardiorespiratory motion, of which 8 provided motion amplitudes in VT patients (n = 94) and 10 described (cardiac/cardiorespiratory) internal target volumes (ITVs) obtained in VT patients (n = 59). Average cardiac motion of target volumes was < 5 mm in all directions, with maximum values of 8.0, 5.2 and 6.5 mm in Superior-Inferior (SI), Left-Right (LR), Anterior-Posterior (AP) direction, respectively. Cardiorespiratory motion of cardiac (sub)structures showed average motion between 5-8 mm in the SI direction, whereas, LR and AP motions were comparable to the cardiac motion of the target volumes. Cardiorespiratory ITVs were on average 120-284% of the gross target volume. Healthy subjects showed average cardiorespiratory motion of 10-17 mm in SI and 2.4-7 mm in the AP direction.

CONCLUSION

This review suggests that despite growing numbers of patients being treated, detailed data on cardiorespiratory motion for STAR is still limited. Moreover, data comparison between studies is difficult due to inconsistency in parameters reported. Cardiorespiratory motion is highly patient-specific even under motion-compensation techniques. Therefore, individual motion management strategies during imaging, planning, and treatment for STAR are highly recommended.

Cardiac radioablation of incessant ventricular tachycardia in patients with terminal heart failure under permanent left ventricular assist device therapy-description of two cases

PURPOSE

Cardiac radioablation (cRA) using a stereotactic single-session radioablative approach has recently been described as a possible treatment option for patients with otherwise untreatable recurrent ventricular tachycardia (VT). There is very limited experience in cRA for patients undergoing left ventricular assist device (LVAD) therapy. We present clinical experiences of two patients treated with cRA for incessant VT under long-term LVAD therapy.

METHODS

Two male patients (54 and 61 years old) with terminal heart failure under LVAD therapy (both patients for 8 years) showed incessant VT despite extensive antiarrhythmic drug therapy and repeated catheter ablation. cRA with a single dose of 25 Gy was applied as a last resort strategy under compassionate use in both patients following an electroanatomical mapping procedure.

RESULTS

Both patients displayed ongoing VT during and after the cRA procedure. Repeated attempts at post-procedural rhythm conversion failed in both patients; however, one patient was hemodynamically stabilized and could be discharged home for several months before falling prey to a fatal bleeding complication. The second patient initially stabilized for a few days following cRA before renewed acceleration of running VT required bilateral ablation of the stellate ganglion; the patient died 50 days later. No immediate side effects of cRA were detected in either patient.

CONCLUSION

cRA might serve as a last resort strategy for patients with terminal heart failure undergoing LVAD therapy and displaying incessant VT. Intermediate- and long-term outcomes of these seriously ill patients often remain poor; therefore, best supportive care strategies should also be evaluated as long as no clear beneficial effects of cRA procedures can be shown. For patients treated with cRA under running ventricular rhythm abnormality, strategies for post-procedural generation of stabilized rhythm have to be established.

Cardiac motion and its dosimetric impact during radioablation for refractory ventricular tachycardia

INTRODUCTION

Cardiac radioablation (CR) is a noninvasive treatment option for patients with refractory ventricular tachycardia (VT) during which high doses of radiation, typically 25 Gy, are delivered to myocardial scar. In this study, we investigate motion from cardiac cycle and evaluate the dosimetric impact in a cohort of patients treated with CR.

METHODS

This retrospective study included eight patients treated at our institution who had respiratory-correlated and ECG-gated 4DCT scans acquired within 2 weeks of CR. Deformable image registration was applied between maximum systole (SYS) and diastole (DIAS) CTs to assess cardiac motion. The average respiratory-correlated CT (AVG_resp ) was deformably registered to the average cardiac (AVG_cardiac ), SYS, and DIAS CTs, and contours were propagated using the deformation vector fields (DVFs). Finally, the original treatment plan was recalculated on the deformed AVG_resp CT for dosimetric assessment.

RESULTS

Motion magnitudes were measured as the mean (SD) value over the DVFs within each structure. Displacement during the cardiac cycle for all chambers was 1.4 (0.9) mm medially/laterally (ML), 1.6 (1.0) mm anteriorly/posteriorly (AP), and 3.0 (2.8) mm superiorly/inferiorly (SI). Displacement for the 12 distinct clinical target volumes (CTVs) was 1.7 (1.5) mm ML, 2.4 (1.1) mm AP, and 2.1 (1.5) SI. Displacements between the AVG_resp and AVG_cardiac scans were 4.2 (2.0) mm SI and 5.8 (1.4) mm total. Dose recalculations showed that cardiac motion may impact dosimetry, with dose to 95% of the CTV dropping from 27.0 (1.3) Gy on the AVG_resp to 20.5 (7.1) Gy as estimated on the AVG_cardiac .

CONCLUSIONS

Cardiac CTV motion in this patient cohort is on average below 3 mm, location-dependent, and when not accounted for in treatment planning may impact target coverage. Further study is needed to assess the impact of cardiac motion on clinical outcomes.

Evaluation of the feasibility of cardiac gating for SBRT of ventricular tachycardia based on real-time ECG signal acquisition

PURPOSE

To investigate the feasibility of cardiac synchronized gating in stereotactic body radiation therapy (SBRT) of ventricular tachycardia (VT) using a real-time electrocardiogram (ECG) signal acquisition.

METHODS AND MATERIALS

Stability of beam characteristics during simulated ECG gating was examined by developing a microcontroller interface to a Varian Clinac iX linear accelerator allowing gating at frequencies and duty cycles relevant to cardiac rhythm. Delivery accuracy was evaluated by measuring dose linearity with an ionization chamber, and flatness and symmetry with a two-dimensional detector array, for different gating windows within typical human cardiac cycle periods. To establish a practical method of gating based on actual ECG signals, an AD8232 Heart Monitor board was used to acquire the ECG signal and synchronize the beam delivery. Real-time cardiac gated delivery measurements were performed for a single 10 × 10 cm² field and for a VT-SBRT plan using intensity-modulated radiation therapy (IMRT).

RESULTS AND DISCUSSION

Dose per monitor unit (MU) values were found to be linear within most gating windows investigated with maximum differences relative to non-gated delivery of <2% for gating windows ≥200 ms and for >10 MUs. Beam profiles for both gated and non-gated modes were also found to agree with maximum differences of 0.5% relative to central axis dose for all sets of beam-on/beam-off combinations. Comparison of dose distributions for intensity-modulated SBRT plans between non-gating and cardiac gating modes provided a gamma passing rate of 97.2% for a 2% 2 mm tolerance.

CONCLUSIONS

Beam output is stable with respect to linearity, flatness, and symmetry for gating windows within cardiac cycle periods. Agreement between dose distributions for VT-SBRT using IMRT in non-gated and cardiac cycle gated delivery modes shows that the proposed methodology is feasible. Technically, gating for delivery of SBRT for VT is possible with regard to beam stability.

Stereotactic Radioablation for Ventricular Tachycardia in the Setting of Electrical Storm

BACKGROUND

Stereotactic body radiotherapy (SBRT) has been reported as a safe and efficient therapy for treating refractory ventricular tachycardia (VT) despite optimal medical treatment and catheter ablation. However, data on the use of SBRT in patients with electrical storm (ES) is lacking. The aim of this study was to assess the clinical outcomes associated with SBRT in the context of ES.

METHODS

This retrospective study included patients who underwent SBRT in the context of ES from March 2020 to March 2021 in one tertiary center (CHU Lille). The target volume was delineated according to a predefined workflow. The efficacy was assessed with the following end points: sustained VT recurrence, VT reduced with antitachycardia pacing, and implantable cardioverter defibrillator shock.

RESULTS

Seventeen patients underwent SBRT to treat refractory VT in the context of ES (mean 67±12.8 age, 59% presenting ischemic heart disease, mean left ventricular ejection fraction: 33.7± 9.7%). Five patients presented with ES related to incessant VT. Among these 5 patients, the time to effectiveness ranged from 1 to 7 weeks after SBRT. In the 12 remaining patients, VT recurrences occurred in 7 patients during the first 6 weeks following SBRT. After a median 12.5 (10.5-17.8) months follow-up, a significant reduction of the VT burden was observed beyond 6 weeks (-91% [95% CI, 78-103]), P<0.0001). The incidence of implantable cardioverter defibrillator shock and antitachycardia pacing was 36% at 1 year.

CONCLUSIONS

SBRT is associated with a significant reduction of the VT burden in the event of an ES; however, prospective randomized control trials are needed. In patients without incessant VT, recurrences are observed in half of patients during the first 6 weeks. VT tolerance and implantable cardioverter defibrillator programming adjustments should be integrated as part of an action plan defined before SBRT for each patient.

Case Report: Treatment Planning Study to Demonstrate Feasibility of Transthoracic Ultrasound Guidance to Facilitate Ventricular Tachycardia Ablation With Protons

Background

Cardiac arrhythmias, such as ventricular tachycardia, are disruptions in the normal cardiac function that originate from problems in the electrical conduction of signals inside the heart. Recently, a non-invasive treatment option based on external photon or proton beam irradiation has been used to ablate the arrhythmogenic structures. Especially in proton therapy, based on its steep dose gradient, it is crucial to monitor the motion of the heart in order to ensure that the radiation dose is delivered to the correct location. Transthoracic ultrasound imaging has the potential to provide guidance during this treatment delivery. However, it has to be noted that the presence of an ultrasound probe on the chest of the patient introduces constraints on usable beam angles for both protons and photon treatments. This case report investigates the possibility to generate a clinically acceptable proton treatment plan while the ultrasound probe is present on the chest of the patient.

Case

A treatment plan study was performed based on a 4D cardiac-gated computed tomography scan of a 55 year-old male patient suffering from refractory ventricular tachycardia who underwent cardiac radioablation. A proton therapy treatment plan was generated for the actual treatment target in presence of an ultrasound probe on the chest of this patient. The clinical acceptability of the generated plan was confirmed by evaluating standard target dose-volume metrics, dose to organs-at-risk and target dose conformity and homogeneity.

Conclusion

The generation of a clinically acceptable proton therapy treatment plan for cardiac radioablation of ventricular tachycardia could be performed in the presence of an ultrasound probe on the chest of the patient. These results establish a basis and justification for continued research and product development for ultrasound-guided cardiac radioablation.

First experimental exploration of real-time cardiorespiratory motion management for future stereotactic arrhythmia radioablation treatments on the MR-linac

Objective.Stereotactic arrhythmia radioablation (STAR) is a novel, non-invasive treatment for refractory ventricular tachycardia (VT). The VT isthmus is subject to both respiratory and cardiac motion. Rapid cardiac motion presents a unique challenge. In this study, we provide first experimental evidence for real-time cardiorespiratory motion-mitigated MRI-guided STAR on the 1.5 T Unity MR-linac (Elekta AB, Stockholm, Sweden) aimed at simultaneously compensating cardiac and respiratory motions.Approach.A real-time cardiorespiratory motion-mitigated radiotherapy workflow was developed on the Unity MR-linac in research mode. A 15-beam intensity-modulated radiation therapy treatment plan (1 × 25 Gy) was created in Monaco v.5.40.01 (Elekta AB) for the Quasar MRI^4Dphantom (ModusQA, London, ON). A film dosimetry insert was moved by combining either artificial (cos⁴, 70 bpm, 10 mm peak-to-peak) or subject-derived (59 average bpm, 15.3 mm peak-to-peak) cardiac motion with respiratory (sin, 12 bpm, 20 mm peak-to-peak) motion. A balanced 2D cine MRI sequence (13 Hz, field-of-view = 400 × 207 mm², resolution = 3 × 3 × 15 mm³) was developed to estimate cardiorespiratory motion. Cardiorespiratory motion was estimated by rigid registration and then deconvoluted into cardiac and respiratory components. For beam gating, the cardiac component was used, whereas the respiratory component was used for MLC-tracking. In-silico dose accumulation experiments were performed on three patient data sets to simulate the dosimetric effect of cardiac motion on VT targets.Main results.Experimentally, a duty cycle of 57% was achieved when simultaneously applying respiratory MLC-tracking and cardiac gating. Using film, excellent agreement was observed compared to a static reference delivery, resulting in a 1%/1 mm gamma pass rate of 99%. The end-to-end gating latency was 126 ms on the Unity MR-linac. Simulations showed that cardiac motion decreased the target's D98% dose between 0.1 and 1.3 Gy, with gating providing effective mitigation.Significance.Real-time MRI-guided cardiorespiratory motion management greatly reduces motion-induced dosimetric uncertainty and warrants further research and development for potential future use in STAR.