- 1Department of Cardiology, Cardiovascular Research Institute Maastricht (CARIM), Maastricht University Medical Center+, Maastricht, Netherlands
- 2Department of Advanced Computing Sciences, Maastricht University, Maastricht, Netherlands
- 3Biomedical Research Institute, Hasselt University, Diepenbeek, Belgium
- 4Department of Physiology, Cardiovascular Research Institute Maastricht (CARIM), Maastricht University, Maastricht, Netherlands
- 5Department of Medical Biology, Amsterdam University Medical Center (UMC), Amsterdam Medical Center (AMC), Amsterdam, Netherlands
- 6Department of Radiology and Nuclear Medicine, Cardiovascular Research Institute Maastricht (CARIM), Maastricht University Medical Center+, Maastricht, Netherlands
- 7Philips Research, Eindhoven, Netherlands
Introduction: Patients with ventricular tachyarrhythmias (VT) are at high risk of sudden cardiac death. When appropriate, catheter ablation is modestly effective, with relatively high VT recurrence and complication rates. Personalized models that incorporate imaging and computational approaches have advanced VT management. However, 3D patient-specific functional electrical information is typically not considered. We hypothesize that incorporating non-invasive 3D electrical and structural characterization in a patient-specific model improves VT-substrate recognition and ablation targeting.
Materials and methods: In a 53-year-old male with ischemic cardiomyopathy and recurrent monomorphic VT, we built a structural-functional model based on high-resolution 3D late-gadolinium enhancement (LGE) cardiac magnetic resonance imaging (3D-LGE CMR), multi-detector computed tomography (CT), and electrocardiographic imaging (ECGI). Invasive data from high-density contact and pace mapping obtained during endocardial VT-substrate modification were also incorporated. The integrated 3D electro-anatomic model was analyzed off-line.
Results: Merging the invasive voltage maps and 3D-LGE CMR endocardial geometry led to a mean Euclidean node-to-node distance of 5 ± 2 mm. Inferolateral and apical areas of low bipolar voltage (<1.5 mV) were associated with high 3D-LGE CMR signal intensity (>0.4) and with higher transmurality of fibrosis. Areas of functional conduction delay or block (evoked delayed potentials, EDPs) were in close proximity to 3D-LGE CMR-derived heterogeneous tissue corridors. ECGI pinpointed the epicardial VT exit at ∼10 mm from the endocardial site of origin, both juxtaposed to the distal ends of two heterogeneous tissue corridors in the inferobasal left ventricle. Radiofrequency ablation at the entrances of these corridors, eliminating all EDPs, and at the VT site of origin rendered the patient non-inducible and arrhythmia-free until the present day (20 months follow-up). Off-line analysis in our model uncovered dynamic electrical instability of the LV inferolateral heterogeneous scar region which set the stage for an evolving VT circuit.
Discussion and conclusion: We developed a personalized 3D model that integrates high-resolution structural and electrical information and allows the investigation of their dynamic interaction during arrhythmia formation. This model enhances our mechanistic understanding of scar-related VT and provides an advanced, non-invasive roadmap for catheter ablation.
1. Introduction
Catheter ablation of arrhythmogenic substrate and triggers of ventricular tachycardia (VT) reduces arrhythmia burden and implantable cardioverter-defibrillator (ICD) shocks (1, 2) but recurrences [30–40% within 1–2 years (3, 4)], costs and complication rates remain high (5, 6). This limited procedural efficacy is attributed to the difficulty of identifying heterogeneous scar corridors related to the critical VT isthmus based on electrogram characteristics and pace mapping, especially during substrate-based ablation.
Recent advances in the field of catheter ablation for scar-related VT have primarily focused on defining target locations for ablation through pre-procedural characterization of the ventricular scar architecture by combining imaging techniques [such as computed tomography (CT) or cardiac magnetic resonance imaging (CMR)] with electrocardiogram (ECG)-VT algorithms. Image-guided or -aided substrate ablation may increase long-term success, decrease recurrence rates and reduce procedural duration (7, 8). Imaging-based approaches mostly focus on structural myocardial targets, ignoring 3D functional information such as activation/propagation maps, repolarization maps, wavefront direction or arrhythmia triggers. Electrocardiographic imaging (ECGI) combines body-surface potential mapping and anatomical imaging to provide a non-invasive reconstruction of electrical activity. It enables a beat-to-beat non-invasive 3D assessment of epicardial activation, repolarization and propagation, besides the pinpointing of the VT exit site with relatively high spatial and temporal resolution (9, 10).
In this study, we combined ECGI with high-resolution 3D dark-blood late gadolinium enhancement (3D-LGE) CMR and cardiac CT in a case with ischemic cardiomyopathy and recurrent monomorphic VT, providing a patient-specific structural-functional model to improve mechanistic understanding and to provide a road map for future personalized catheter ablation for VT. Additionally, we examined the accuracy of non-invasively determined parameters compared to standard invasive VT-substrate characteristics.
2. Materials and methods
We constructed a personalized 3D cardiac model that integrates the individualized structural-functional characteristics of this case to identify mechanisms of arrhythmia formation and to develop an individualized framework for future VT ablation guidance.
The model incorporates (1) 3D dark-blood LGE CMR and cardiac CT delineating the myocardial anatomy and scar (structural aspects); and (2) functional-electrical aspects by non-invasive ECGI. Model outcomes were compared to data from invasive electro-anatomical mapping (EAM), and arrhythmia mechanisms were investigated.
2.1. Patient characteristics
A 53-year-old male was admitted to our hospital for electrical storm due to recurrent hemodynamically tolerated, sustained monomorphic VTs with an inferobasal to mid, septal left-ventricular (LV) origin (Figure 1A) (11). He was previously known with anterior and inferolateral myocardial infarctions (LV ejection fraction 31%), also evident from the 12-lead electrocardiogram (Figure 1B), and coronary artery bypass grafting. The VT was managed with intravenous infusion of amiodarone. A 3D-LGE CMR and cardiac CT scan were performed to aid catheter ablation for VT. Additional electrical information was obtained by body-surface potential mapping before and during the VT ablation procedure.
Figure 1. Ensure that all the figure captions are correct, and that all figures are of the highest quality/resolution. If necessary, you may upload improved figures to the Production Forum. Please note that figures must be cited sequentially, per the author guidelines. The patient’s electrocardiogram during ventricular tachycardia (VT) (A) and sinus rhythm (B). (A) Monomorphic VT demonstrating a wide QRS complex of 151 bpm with a right bundle branch block pattern (transition in V3) and a superior axis. The positive polarities in leads I and aVL suggest an inferomedial to inferobasal, septal origin (47). (B) Sinus rhythm, 72 bpm, left axis deviation, a QRS width of 112 ms, pathological Q waves in inferior leads and V3–V6 and negative T waves in inferior leads and in V4–V6. Q waves are indicative of previous inferior/apical infarctions.
2.2. A 3D dark-blood late gadolinium-enhancement CMR
The patient underwent a free-breathing 3D-LGE CMR with high isotropic resolution (acquired resolution 1.6 mm × 1.6 mm × 1.6 mm, reconstructed 0.8 mm × 0.8 mm × 0.8 mm) and dark-blood nulling (12) on a 1.5 T MR system (Ingenia, Philips Healthcare, Best, Netherlands). This allowed for more accurate scar demarcation and image quality with respect to conventional bright-blood nulled LGE CMR (13). Details of the dark-blood LGE mechanism without using additional magnetization preparation has were described earlier (14).
Cardiac magnetic resonance-image processing was performed using a research software tool (ADAS3D Medical SL, Barcelona, Spain) (15). Pixel-signal intensity (PSI) maps were created from the 3D delayed enhancement acquisitions, and normalized for maximal PSI. A normalized PSI-based cut-off above 0.6 was defined as dense scar; border zone tissues were allocated in case of values between 0.4 and 0.6% (16). The LV myocardium was divided in nine three-dimensional transmural layers color-coded for scar architecture. Heterogeneous tissue corridors, defined as continuous corridors of border zone tissue connecting two regions of healthy tissue, bordered by dense scar and/or an anatomical barrier, were automatically calculated and visualized on the shells by the research tool (15). Subsequently, 3D-LGE CMR scar transmurality was defined as the percentage of transmural layers having either dense scar (weight 1) or border zone (weight 0.5). Scar transmurality was projected onto the endocardial shell.
2.3. Computed tomography and electrocardiographic imaging
Prior to the VT ablation procedure, an ECG-gated helical dual-source CT scan (Somatom Force, Siemens Healthineers, Forchheim, Germany) using intravenous iodine contrast medium was performed during breath-hold, to visualize cardiac geometry. Additionally, the patient underwent a low-dose thoracic CT scan to locate the body-surface electrodes, required for ECGI.
The ECGI was performed as previously described and validated (17, 18). Briefly, 224 Ag-AgCl electrodes (BioSemi, Amsterdam, The Netherlands) were attached to the patient’s torso to record body-surface potentials at a 2048 Hz sampling rate starting 20 min prior to, and during, the electrophysiological study. From the thoracic CT scan, a 224-electrode torso was manually digitized. A ∼2000-node diastolic cardiac geometry was semi-automatically segmented (19). Body-surface potential measurements were performed during sinus rhythm, atrial pacing, during programed ventricular stimulation for VT induction and for evoked delayed potential (EDP) identification (20), and during VT.
Epicardial ventricular unipolar electrograms (UEGs) were reconstructed for selected beats through previously detailed methods (17). For each UEG, activation times (ATs) were determined from the steepest downslope (of the UEG QRS complex); this was done with a spatiotemporal approach that considers the spatial flow of current (21). ATs of each beat were calculated with respect to the average of the first 25 epicardial activations according to ECGI. UEGs were disregarded at the valvular base of the ventricles. VT dynamics were assessed by calculating the absolute activation time differences between consecutive VT beats.
2.4. Catheter ablation procedure for VT
The patient underwent continued body-surface potential mapping during the VT ablation procedure (Supplementary Figure 1). The procedure was performed under general anesthesia without preemptive hemodynamic support. The CARTO navigation system (Biosense Webster Inc., Diamond Bar, CA, USA) was used for guidance of the procedure. Both retrograde and anterograde access were obtained. The ADAS3D-derived anatomy and scar architecture were merged with the aortic root, main stem and pulmonary veins. After positioning of a quadripolar diagnostic catheter (6 F, CRD-2™, Abbott, Chicago, IL, USA) in the right-ventricular (RV) apex, programed RV stimulation was performed at three driving cycle lengths (600, 500, and 400 ms) with up to three ventricular extrastimuli (≥200 ms) to induce sustained VT (>30 s or necessitating termination because of hemodynamic instability) and to determine the RV effective refractory period (ERP). During right-atrial (RA) pacing, a detailed electroanatomical voltage map was obtained using a high-density mapping (PentaRay, Biosense Webster Inc., Diamond Bar, CA, USA) and an irrigated contact-force sensing 3.5 mm-tip catheter (SmartTouch, Biosense Webster Inc., Diamond Bar, CA, USA). Near-field peak-to-peak bipolar voltage amplitude below 1.5 mV indicated scar, with dense scar <0.5 mV and border zone 0.5–1.5 mV. EDPs, characterized by low-amplitude near-field potentials with 10-ms delay or block after RV extrastimulus (S2 50 ms above the ERP), were tagged in the 3D system, as previously described (22). The VT site-of-origin was identified using pace mapping aimed at a mean correlation coefficient between the VT-QRS morphology and the paced-QRS of >90% (23). Substrate modification was performed aiming at complete EDP elimination and targeting the VT site-of-origin using radiofrequency energy delivery (power 50 W, temperature limit 43°C, flow rate 20–30 mL/min). After the last lesion set, the substrate was reassessed for the abolition of EDPs, and programed ventricular stimulation was repeated to test for VT inducibility.
2.5. Postprocedural image integration and analysis
The geometries from 3D-LGE CMR, CT, ECGI and invasive EAM were manually aligned and registered digitally in 3D Slicer (24). The most important anatomic landmarks used for integration were: the LV and RV epicardium, LV apex and aorta from CT; LV apex, epicardium, endocardium and aorta from 3D-LGE CMR; and LV endocardium, apex and aorta from EAM. The accuracy of the image integration was based on Euclidean node-to-node distance calculation and visual evaluation. The EAM geometry was converted isotropic (i.e., equal node-to-node distances throughout the geometry) to ensure that all regions of the heart had equal weight in downstream analyses. Further analyses were performed in MATLAB R2020b (MathWorks, Natick, MA, USA).
Information from 3D-LGE CMR and ECGI was combined to examine triggers, substrate, and their interactions. 3D-LGE CMR scar delineation was compared to areas of low invasive bipolar voltage, by finding the corresponding nearest neighbors of the CMR anatomy and the invasive electroanatomic anatomy. Any neighbors with a corresponding distance of >10 mm from each other were disregarded. The colocalization of endocardial EDPs with MRI-derived heterogeneous tissue corridors was also evaluated, disregarding EDPs that were tagged >10 mm from the endocardial 3D-LGE CMR shell.
3. Results
3.1. Ablation procedure and outcome
A sustained monomorphic VT arising from the inferobasal LV, closely resembling the documented arrhythmia, was repeatedly induced by programed electrical stimulation (by minimally two ventricular premature beats) and led to hemodynamic compromise. The high-density LV-endocardial bipolar voltage map (5,375 points) revealed two remote areas of scar in the lateral and apical regions of the LV, colocalizing with areas of increased endocardial PSI based on 3D-LGE CMR. 21 EDPs were tagged. Pace mapping at the inferobasal LV resulted in a 92% match with the induced VT. Radiofrequency delivery targeted at the inferolateral EDPs and the site of VT origin, eliminated the EDPs and rendered the patient non-inducible at the end of the procedure. Total procedural time was 330 min; fluoroscopy time 10 min. The procedure was without any complications. A single-chamber ICD was implanted afterward. No VT had recurred at 20 months follow-up.
3.2. Contrast-enhanced CT and 3D dark-blood late gadolinium-enhancement CMR
Contrast-enhanced multidetector CT showed wall thinning in the inferolateral LV and apex (Figure 2A and Supplementary Figure 2). The CMR revealed a partly transmural infarcted area in the inferolateral and apical segments with aneurysm formation (Figure 2B and Supplementary Videos 1, 2). Dense scar and border zone, as detected through CMR, consisted of non-uniform transmural patterns (Figure 2C). Heterogeneous tissue corridors, identified from CMR (Figure 2D) were localized in the basal to mid (infero) lateral, and apical aspects of the LV.
Figure 2. Cardiac imaging of the patient prior to the ablation procedure. (A) Contrast-enhanced dual-source computed tomography (CT), performed for anatomical reference showing wall thinning in the apical aneurysm and laterobasal aspect of the left ventricle (LV). (B) 3D dark-blood late gadolinium enhancement cardiac magnetic resonance imaging (LGE CMR) LGE CMR showing areas with hyperenhancement in the laterobasal and apical aspects of the LV. (C) Two cross-sectional views of LV scar, detected by 3D dark-blood LGE CMR. (D) CMR-derived healthy and scar tissue, including corridors and their centerlines, depicted in nine transmural layers (interlayer distance approximately 0.8 mm) of the ventricular myocardium. Colors correspond to colors in panel (C).
Manual image integration of invasive EAM and 3D-LGE CMR endocardial geometry led to a mean Euclidean node-to-node distance 5 ± 2 mm. Increased PSI (>0.4) of the innermost endocardial layer of 3D dark-blood LGE CMR closely resembled reduced endocardial bipolar voltage (<1.5 mV) using previously defined cutoffs (15), while areas of non-agreement were near the scar-to-healthy border (Figure 3A). Both modalities agreed (both “scar” or both “healthy”) for 68% of tissue, and disagreed for the remaining 32%, mostly concerning the border zone. According to both modalities, the majority of scar consisted of border zone (76% according to EAM, 85% according to 3D-LGE CMR). Predicting bipolar scar through the CMR endocardial signal intensity rendered a receiver operating characteristic (ROC) area under the curve of 0.72. Further comparison of 3D-LGE CMR scar transmurality to EAM bipolar voltage showed that low voltages were associated with higher scar transmurality, and vice versa (Figure 3B).
Figure 3. Comparison of 3D-late gadolinium enhancement cardiac magnetic resonance imaging (LGE CMR) to electroanatomical mapping (EAM). (A) Invasive EAM bipolar voltage map (left) and 3D dark-blood LGE CMR endocardial signal intensity (middle), and the comparison of both (right). Areas of agreement regarding healthy tissue are shown in green, agreement regarding scar tissue are shown in red, and differences between both maps are shown in gray. (B) CMR scar transmurality compared to EAM bipolar voltage. Higher transmurality of scar (detected by 3D-LGE CMR) was associated with lower bipolar voltage in EAM.
A total of 17/21 EDPs were within 10 mm of the endocardial CMR-wall. 15/17 (88%) of these EDPs were located in an area of endocardial scar detected by 3D dark-blood LGE CMR (Figures 4A, B). EDPs were often identified within or close to (transmural) CMR-derived corridors. When projecting all EDPs and corridors onto the endocardial CMR layer, the average Euclidean distance from EDP to the nearest corridor was 3 (0–6) mm [median (first and third quartiles)].
Figure 4. A 3D-late gadolinium enhancement cardiac magnetic resonance imaging (LGE CMR) heterogeneous tissue corridors and evoked delayed potentials (EDPs). (A) Annotation of EDPs during the electroanatomical mapping (EAM) procedure, compared to a superimposed anatomical view of transmural heterogeneous tissue corridors detected through 3D LGE CMR. All EDPs except one coincided with ventricular scar or border zone. EDPs qualitatively coincided well with corridors close to the endocardium. (B) EDPs on locations 1–4 in panel A. (C) Investigating of ventricular tachycardia (VT) working mechanisms through structural-functional image overlay. Transmural layers of scar and border zone with centerlines of slow-conducting heterogeneous tissue corridors are superimposed and colored. VT exits were determined by pace mapping during the procedure. Invasively determined endocardial (brown) and non-invasively determined epicardial (white) exits coincide at the end of epicardial and endocardial corridors. Corridors of interest leading to the VT exit have an increased width. Ablation of the substrate (purple), including the identified corridors, rendered the patient arrhythmia-free.
3.3. Structural-functional substrate characterization
Structural-functional overlay of the 3D-LGE CMR and ECGI data (body-surface potential mapping during the ablation procedure) allowed the investigation of the VT trajectory and mechanism. In Figure 4C, the endocardial (invasive EAM) and epicardial (ECGI) locations of the VT exit according to pace mapping are shown, together with transmural fibrosis and PSI-derived heterogeneous tissue corridors. ECGI pinpointed the VT exit in close proximity (spatial accuracy 10/14/13 mm) to the invasive endocardial site with the highest pace-mapping agreement (0.90/0.92/0.92, respectively). When projecting the ECGI-exit onto the endocardial geometry, i.e., omitting the endocardial-to-epicardial distance, the distance between invasive and non-invasive VT exit decreased to 5/12/12 mm, respectively. Both epicardial and endocardial heterogeneous tissue corridors were adjacent to the VT exit, possibly serving as slow-conducting channels (isthmuses) for the VT. Radiofrequency ablation to, amongst other sites, the entrance of these corridors, rendered the patient arrhythmia-free (20 months follow-up).
Although the monomorphic VT (with exception of its first beat) appeared stable on the ECG (Figure 5A), ECGI revealed significant beat-to-beat differences in activation patterns during the onset of VT (Figure 5B). The VT stabilized over time, causing beat-to-beat differences in its activation pattern to decline. The anterolateral area of early activation migrated inferolaterally as the VT evolved. This area of shifting early activation correlated with CMR-detected corridors leading to the VT exit (Figure 4B). Beat-to-beat-differences in activation pattern were present during all inductions (episodes) of VT, higher than during sinus rhythm and ventricular pacing (Figure 5C).
Figure 5. Dynamics of epicardial activation patterns through electrocardiographic imaging (ECGI). (A) ventricular tachycardia (VT) during procedure in which beats 1 through 7 are depicted. (B) Activation patterns of VT beats 1 through 7, with absolute difference in activation time (AT) of consecutive beats shown below. The pentagram in VT-1 illustrates the origin of the premature ventricular complex (PVC) (Supplementary Figure 3). For reference, dynamics of activation patterns in sinus rhythm (SR) and pacing are also shown. (C) Average absolute difference in activation times of consecutive beats, over the entire epicardium. Dynamics in activation patterns decrease over time for several episodes (inductions) of VT, while SR and ventricular pacing show stable activation patterns.
Prior to the electrophysiological study, a premature ventricular complex (PVC) was recorded with ECGI (Supplementary Figure 3A). The PVC origin matched significantly, although not completely, with the VT exit. ECGI pinpointed the PVC to the lateral LV wall (Supplementary Figure 3B). The crowded activation isochrones superior to the first-activated area suggest a conduction block in the basolateral LV (blue arrow). CMR-detected heterogeneous tissue corridors (numbers 1–3, Supplementary Figure 3C) connected the PVC and VT origins and coincided with the area of shifting early activation during VT (Supplementary Figure 3C). Hence, we hypothesize the PVC could be involved in VT initiation.
4. Discussion
In this proof-of-concept study, we developed a 3D structural-functional model for scar-related VT by combining 3D dark-blood LGE CMR and ECGI (of PVC and VT) with high spatial and temporal resolution. This model was developed in an exemplary case with recurrent monomorphic VT in the setting of ischemic cardiomyopathy. Patient-specific VT substrates, triggers, and their interactions were detailed, in a non-invasive manner, serving as a roadmap for future individualized ablation strategies. Integrating high-density invasive EAM advanced the understanding of possible VT mechanisms.
In our integrated model, the (invasive and non-invasive) VT exit colocalized at the end of CMR-detected heterogeneous tissue corridors. Furthermore, the patient’s spontaneous PVC originated in close proximity to the VT exit and CMR-detected corridors. This suggests that the PVC may be involved in patient’s VT formation. Moreover, EDPs were found close to heterogeneous tissue corridors, suggesting a strong link between functional endocardial conduction slowing/block and structural border zone channels.
4.1. 3D dark-blood LGE CMR
Acquiring 3D-LGE CMR with high isotropic resolution visualizes myocardial fibrosis from every desired direction compared to a fixed viewing angle when using (non-isotropic) 3D or standard 2D imaging (25). Subtle scar architectures, such as submillimeter tissue corridors, may be missed using these latter acquisition techniques. Additionally, the dark-blood nulling CMR, recently validated against histopathology in an experimental large-animal model (13), allows for improved detection and quantification of (sub) endocardial fibrosis, often encountered in ischemic etiologies, compared to conventional bright-blood LGE (14). We have shown that 3D dark-blood LGE CMR has a reasonable correlation, at least qualitatively, with the endocardial bipolar voltage. We found a 68% agreement between reduced (<1.5 mV) endocardial bipolar voltage and increased (>0.4) 3D-LGE CMR-signal intensity when applying a point-by-point comparison between both modalities. Areas of non-agreement were predominantly found at the healthy-to-scar-border zone. The ROC area under the curve to define areas of low voltage through 3D-LGE CMR was 0.72, which compared with the low range of values found in a recent study, although different methods of analysis may have affected the results (26). The observed electrostructural mismatch stems from the fundamentally different assessment of scar characteristics, with varying spatial resolution and dissimilar gating for cardiac phases and respiratory motion, besides the variable impact of modulators (e.g., contact force, fiber orientation, wavefront direction, and sympathetic hyperinnervation). The scar border zone may be especially susceptible to these factors.
Additionally, we have shown a good correlation between the 3D dark-blood LGE CMR-derived myocardial heterogeneous tissue corridors and the invasively measured areas of functional electrical conduction slowing or block (15). These corridors may comprise slow-conducting pathways, or protected isthmuses, that may facilitate reentrant excitation (15) upon unidirectional conduction delay or block (22).
4.2. ECGI
ECGI has been validated in multiple studies. For a recent overview we refer to Stoks and Cluitmans (27). Its strength lies in the reconstruction of (paced) rhythm of single origin (28). In a structurally normal heart, ECGI was more accurate than the 12-lead ECG to pinpoint the origin of VT or PVC (29). The 12-lead ECG was even less precise in the presence of myocardial scarring (30). We have shown that ECGI-based reconstructions localized the site of ventricular activation (adjacent to scar border zone) by pace mapping at ∼10 mm from the endocardial location (if endocardial-to-epicardial distance was omitted). This was more accurate than in a recent study using a commercial system (30) but comparable to other reports (27).
Moreover, ECGI can generate a complete 3D-epicardial activation map in a single beat, while several VT cycle lengths are needed during invasive contact mapping EAM to construct an activation or propagation map (which is often not possible since VT is often polymorphic or hemodynamically not tolerated). This may be relevant for investigating the dynamic behavior of VT, especially in polymorphic VT but also during the onset of monomorphic VTs. In our case, ECGI uncovered marked initial activation dynamicity, only stabilizing after several beats. An area of electrical instability (shifting early activation) during VT was present close to the VT exit and structural-functional corridors (Figures 4C, 5B).
4.2. Structural-functional VT substrate characterization
Non-invasive modalities identifying critical components of the VT circuit to guide catheter ablation have evolved considerably over the last decades. Historical milestones in the human heart are shown in Figure 6.
Figure 6. Timeline of invasive and non-invasive techniques to identify critical components of the ventricular tachycardia (VT) circuit and guide catheter ablation (7, 15, 31–33, 35–38, 43, 44). Note that the time axis is not linear, which is reflected by the color gradients. ECG, electrocardiogram; PVC, premature ventricular complex; EAM, electroanatomical mapping; ECGI, electrocardiographic imaging; CMR, cardiac magnetic resonance imaging; CT, computed tomography.
In 1972, Wellens et al. found that, in (mostly) post-infarction patients, VTs were based on reentrant excitation, that could be terminated through electrical stimulation (31). In 1981, Josephson et al. compared 12-lead ECG QRS morphologies of VTs to endocardial VT exit sites in 34 patients (32), demonstrating that the presence of myocardial scar diminishes the ECG accuracy to pinpoint the endocardial VT exit.
In 1999, Wittkampf et al. first validated an impedance-based 3D catheter navigation system in humans for EAM (33). In 2000, De Groot et al. introduced a 3D-catheter localization system using ultrasound transducers to guide VT ablation (34), that allowed real-time display of the catheter tip and repositioning to previously marked sites. In the same year, Marchlinski et al. introduced the magnetic mapping-based CARTO system to perform linear substrate ablation which proved effective to significantly reduce VT recurrence (35). These 3D-navigation systems could be considered the foundation of the modern-day EAM.
In 2005, pioneering work on image integration was done by Dickfeld et al. who coregistered MRI acquisitions with the electromagnetic catheter positioning system to guide catheter navigation in the human right atrium and ventricle (36). In the same year, Ghanem et al. showed that near-simultaneous ECGI-derived (non-invasive) electrograms had a moderate correlation with invasive potentials recorded in sinus rhythm and during endocardial and epicardial pacing in patients (37). Five years later, Tian et al. (38) coregistered contrast-enhanced CT relying on distinct anatomic, dynamic and perfusion characteristics indicative of myocardial scar and border zone with electrogram voltage obtained during endocardial point-by-point mapping. A high segmental accuracy was reported. Though contrast-enhanced CMR still has superior contrast resolution compared to CT, the latter is more easily applicable for scar demarcation in patients with an ICD and allows for increased temporal and spatial resolution (38). Furthermore, wall thickness analysis with CT accurately predicted isthmuses in postinfarction VT (39).
In 2008, Codreanu et al. performed image integration of EAM and LGE CMR to conclude that transmural scar (by CMR) could be detected by spiky electrograms, reduced voltage amplitudes or prolonged bipolar electrogram durations (40). A year later, Desjardins et al. investigated the relationship between reduced bipolar voltage in EAM and automatically segmented scar in LGE CMR, to find that sites critical to VT in humans are located within areas of delayed enhancement in CMR. In 2011, Wijnmaalen et al. performed a head-to-head comparison of LGE CMR with electroanatomical voltage mapping through image integration. They found that increased scar transmurality led to decreased endocardial bipolar voltage and that non-transmural and border zone scar was undetected by EAM (41). The combination of LGE CMR with multidetector CT is superior for structure-function relationship determination in scar-related VT (42). With the availability of 3D-LGE CMR acquisitions with high isotropic resolution, Andreu et al. (15) identified heterogeneous tissue corridors that colocalize in 79% with invasively recorded conduction channels (2015). VT ablation procedures guided by CMR led to increased procedure success, and decreased VT recurrence, inducibility and procedure time (7). In recent years, promising computational approaches have arisen. For example, Prakosa et al. performed CMR-based computational electrophysiological modeling, which allowed for prospective ablation guidance in a small cohort (43). Moreover, Lozova et al. used machine-learning algorithms to pinpoint ablation targets (44).
The aforementioned advances in the field of catheter-based ablation of scar-related VT have primarily focused on defining target locations for ablation through analysis of structural imaging modalities, typically CMR (43, 44) and/or CT (42). Though such structural detailing holds promise and allows to non-invasively investigate the substrate and perform virtual stress tests, these strategies usually do not consider patient-specific 3D functional electrical parameters. Our current proof-of-concept patient-tailored model encompasses 3D-functional and structural information and their interactions. As an example, we have shown that invasively- and non-invasively determined VT exits align colocalize with endings of heterogeneous tissue corridors, and are connected with the PVC site of origin. Moreover, we have demonstrated dynamic electrical instability of the heterogeneous scar region where the scar heterogeneity set the stage for an evolving VT circuit.
Combined, we argue that the integration (during spontaneous or induced PVCs and VT, including their spatiotemporal dynamics) of ECGI with high resolution 3D-LGE dark-blood CMR may be synergistic for the development of (non-invasive) individualized VT models to improve our understanding of VT mechanisms and treatment.
4.3. Future perspectives
Prospective studies are needed to investigate the procedural efficacy and cost-effectiveness of such personalized structural-functional models. This would require optimization of the modeling pipeline, including automated segmentation and merging of different modalities and automation of analyses, which is currently being developed.
Ideally, in future cases, ECGI recordings should be performed pre-procedurally during spontaneous VT or non-invasive programed stimulation (in the absence of deep sedation or general anesthesia) to allow preprocedural structural-functional analysis of the VT trajectory and its mechanism. This would guide targeted RF delivery to critical components of the VT circuit without the need for extensive intraprocedural voltage and pace mapping. For instance, CMR-derived heterogeneous tissue corridors (aligning to EDPs) connecting to ECGI-derived VT exit sites could be ablated. ECGI activation maps should be compared to invasively acquired (endo- and epicardial) VT activation/propagation maps to investigate excitation accuracy, VT dynamicity and trajectories.
In absence of an ICD, CMR may replace CT for localization of ECGI electrodes and the heart geometry, thereby omitting radiation. Beyond ischemic cardiomyopathies, this model could be extended other cardiac diseases with increased susceptibility for VT/VF, including non-ischemic and idiopathic etiologies. In the latter, ECGI uncovered previously undetected functional repolarization abnormalities, crucial for arrhythmia formation (18). Next, the model could be improved by the incorporation of static and dynamic ventricular repolarization. Furthermore, integration of our model with existing computational frameworks (digital twinning) may aid in identifying ablation targets and predict VT recurrences (risk stratification). Finally, structural-functional image modeling may be of value in novel emerging non-invasive ablation strategies, such as stereotactic cardiac radioablation (45, 46).
4.4. Study limitations
This was a proof-of-concept study based on one exemplary case. We aimed to first develop the model before recruiting more patients, which is currently ongoing. Firstly, although the ECGI VT activation map and exit site were acquired during the procedure, the study data were merged and analyzed retrospectively. The VT ablation procedure was not guided by our model. Secondly, anatomies were merged based on visual inspection. Advanced merging methods may further improve results of invasive EAM vs. ECGI/3D-LGE CMR comparisons (26). Thirdly, we used an epicardial-only formulation of ECGI (17) which does not actively consider transmural fibrosis, although this method is most widely applied and the most extensively validated (27). Lastly, we used bright-blood validated LGE CMR thresholds for scar delineation in the ADAS3D program. However, upon further analysis, in this patient, altering this threshold affected EAM vs. CMR scar comparisons by <1%.
5. Conclusion
In this proof-of-concept study, we combined 3D-LGE CMR, CT and ECGI imaging modalities to develop a patient-tailored high-resolution non-invasive structural-functional model to investigate VT substrate and trigger in a case of recurrent VT in the setting of ischemic cardiomyopathy. The model accurately pinpointed the VT exit site, collocated areas of structural and functional corridors and provided insight into scar transmurality and endocardial voltage characteristics. Beyond existing methods, our model encompasses beat-by-beat functional electrical information, such as dynamic 3D activation patterns of spontaneous ventricular ectopics and the VT circuit.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The studies involving human participants were reviewed and approved by the METC azM/UM. The patients/participants provided their written informed consent to participate in this study.
Author contributions
JS, MC, PV, and RTB: study setup. JS: ECGI data acquisition and processing, data integration and analysis, and manuscript draft. YK: CT wall thickness analysis. RH and SG: CMR data acquisition. RTB: EAM and clinical data acquisition. JS, MC, BH, BB, PV, and RTB: conceptualization. MC, BH, RH, SG, YK, KV, BB, PV, and RTB: critical review of manuscript. All authors have read and approved the final manuscript and data interpretation.
Funding
This study was supported by the Special Research Fund (BOF) of Hasselt University (BOF17DOCMA15) and the Maastricht University Medical Center (MUMC+) to JS, the Hein Wellens Foundation, Health Foundation Limburg (Maastricht, The Netherlands), and a Veni grant from the Netherlands Organization for Scientific Research (TTW16772) to MC, the Netherlands CardioVascular Research Initiative (CVON2017-13 VIGILANCE and CVON2018B030 PREDICT2), Den Haag, The Netherlands to PV, and a Veni grant from the Netherlands Organization for Scientific Research (NWO/ZonMw 0915016181013) to RTB.
Conflict of interest
MC is part-time employed by Philips Research.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcvm.2023.1112980/full#supplementary-material
References
1. Reddy V, Reynolds M, Neuzil P, Richardson A, Taborsky M, Jongnarangsin K, et al. Prophylactic catheter ablation for the prevention of defibrillator therapy. N Engl J Med. (2007) 357:2657–65. doi: 10.1056/NEJMoa065457
2. Kuck K, Schaumann A, Eckardt L, Willems S, Ventura R, Delacrétaz E, et al. Catheter ablation of stable ventricular tachycardia before defibrillator implantation in patients with coronary heart disease (VTACH): a multicentre randomised controlled trial. Lancet. (2010) 375:31–40. doi: 10.1016/S0140-6736(09)61755-4
3. Santangeli P, Muser D, Maeda S, Filtz A, Zado E, Frankel D, et al. Comparative effectiveness of antiarrhythmic drugs and catheter ablation for the prevention of recurrent ventricular tachycardia in patients with implantable cardioverter-defibrillators: a systematic review and meta-analysis of randomized controlled trials. Heart Rhythm. (2016) 13:1552–9. doi: 10.1016/j.hrthm.2016.03.004
4. Briceño D, Romero J, Villablanca P, Londoño A, Diaz J, Maraj I, et al. Long-term outcomes of different ablation strategies for ventricular tachycardia in patients with structural heart disease: systematic review and meta-analysis. EP Eur. (2018) 20:104–15. doi: 10.1093/europace/eux109
5. Coyle K, Coyle D, Nault I, Parkash R, Healey J, Gray C, et al. Cost effectiveness of ventricular tachycardia ablation versus escalation of antiarrhythmic drug therapy. JACC Clin Electrophysiol. (2018) 4:660–8. doi: 10.1016/j.jacep.2018.01.007
6. Calkins H, Epstein A, Packer D, Arria A, Hummel J, Gilligan D, et al. Catheter ablation of ventricular tachycardia in patients with structural heart disease using cooled radiofrequency energy. J Am Coll Cardiol. (2000) 35:1905–14. doi: 10.1016/S0735-1097(00)00615-X
7. Soto-Iglesias D, Penela D, Jáuregui B, Acosta J, Fernández-Armenta J, Linhart M, et al. Cardiac magnetic resonance-guided ventricular tachycardia substrate ablation. JACC Clin Electrophysiol. (2020) 6:436–47. doi: 10.1016/j.jacep.2019.11.004
8. Hendriks A, Kis Z, Glisic M, Bramer W, Szili-Torok T. Pre-procedural image-guided versus non-image-guided ventricular tachycardia ablation—a review. Neth Heart J. (2020) 28:573–83. doi: 10.1007/s12471-020-01485-z
9. Cronin E, Bogun F, Maury P, Peichl P, Chen M, Namboodiri N, et al. 2019 HRS/EHRA/APHRS/LAHRS expert consensus statement on catheter ablation of ventricular arrhythmias. EP Eur. (2019) 21:1143–4. doi: 10.1093/europace/euz132
10. Wang Y, Cuculich P, Zhang J, Desouza K, Vijayakumar R, Chen J, et al. Noninvasive electroanatomic mapping of human ventricular arrhythmias with electrocardiographic imaging. Sci Transl Med. (2011) 3:98ra84. doi: 10.1126/scitranslmed.3002152
11. Andreu D, Fernández-Armenta J, Acosta J, Penela D, Jáuregui B, Soto-Iglesias D, et al. axis–based algorithm to identify the origin of scar-related ventricular tachycardia in the 17-segment American heart association model. Heart Rhythm. (2018) 15:1491–7. doi: 10.1016/j.hrthm.2018.06.013
12. Holtackers R, Gommers S, Van De Heyning C, Mihl C, Smink J, Higgins D, et al. Steadily increasing inversion time improves blood suppression for free-breathing 3D late gadolinium enhancement MRI with optimized dark-blood contrast. Invest Radiol. (2021) 56:335–40. doi: 10.1097/RLI.0000000000000747
13. Holtackers R, Gommers S, Heckman L, Van De Heyning C, Chiribiri A, Prinzen F. Histopathological validation of dark-blood late gadolinium enhancement MRI without additional magnetization preparation. J Magn Reson Imaging. (2022) 55:190–7. doi: 10.1002/jmri.27805
14. Holtackers R, Chiribiri A, Schneider T, Higgins D, Botnar R. Dark-blood late gadolinium enhancement without additional magnetization preparation. J Cardiovasc Magn Reson. (2017) 19:64. doi: 10.1186/s12968-017-0372-4
15. Andreu D, Ortiz-Pérez J, Fernández-Armenta J, Guiu E, Acosta J, Prat-González S, et al. 3D delayed-enhanced magnetic resonance sequences improve conducting channel delineation prior to ventricular tachycardia ablation. EP Eur. (2015) 17:938–45. doi: 10.1093/europace/euu310
16. Andreu D, Berruezo A, Ortiz-Pérez J, Silva E, Mont L, Borràs R, et al. Integration of 3D electroanatomic maps and magnetic resonance scar characterization into the navigation system to guide ventricular tachycardia ablation. Circ Arrhythm Electrophysiol. (2011) 4:674–83. doi: 10.1161/CIRCEP.111.961946
17. Cluitmans M, Bonizzi P, Karel J, Das M, Kietselaer B, de Jong M, et al. In vivo validation of electrocardiographic imaging. JACC Clin Electrophysiol. (2017) 3:232–42. doi: 10.1016/j.jacep.2016.11.012
18. Cluitmans M, Bear L, Nguyên U, van Rees B, Stoks J, ter Bekke RMA, et al. Noninvasive detection of spatiotemporal activation-repolarization interactions that prime idiopathic ventricular fibrillation. Sci Transl Med. (2021) 13:eabi9317. doi: 10.1126/scitranslmed.abi9317
19. Ecabert O, Peters J, Walker M, Ivanc T, Lorenz C, von Berg J, et al. Segmentation of the heart and great vessels in CT images using a model-based adaptation framework. Med Image Anal. (2011) 15:863–76. doi: 10.1016/j.media.2011.06.004
20. Jackson N, Gizurarson S, Viswanathan K, King B, Massé S, Kusha M, et al. Decrement evoked potential mapping. Circ Arrhythm Electrophysiol. (2015) 8:1433–42. doi: 10.1161/CIRCEP.115.003083
21. Cluitmans M, Coll-Font J, Erem B, Bear L, Nguyên U, ter Bekke R, et al. Spatiotemporal approximation of cardiac activation and recovery isochrones. J Electrocardiol. (2022) 71:1–9. doi: 10.1016/j.jelectrocard.2021.12.007
22. de Riva M, Naruse Y, Ebert M, Androulakis A, Tao Q, Watanabe M, et al. Targeting the hidden substrate unmasked by right ventricular extrastimulation improves ventricular tachycardia ablation outcome after myocardial infarction. JACC Clin Electrophysiol. (2018) 4:316–27. doi: 10.1016/j.jacep.2018.01.013
23. de Chillou C, Groben L, Magnin-Poull I, Andronache M, Abbas M, Zhang N, et al. Localizing the critical isthmus of postinfarct ventricular tachycardia: the value of pace-mapping during sinus rhythm. Heart Rhythm. (2014) 11:175–81. doi: 10.1016/j.hrthm.2013.10.042
24. Fedorov A, Beichel R, Kalpathy-Cramer J, Finet J, Fillion-Robin J, Pujol S, et al. 3D slicer as an image computing platform for the quantitative imaging network. Magn Reson Imaging. (2012) 30:1323–41. doi: 10.1016/j.mri.2012.05.001
25. Holtackers R, Emrich T, Botnar R, Kooi M, Wildberger J, Kreitner K. Late gadolinium enhancement cardiac magnetic resonance imaging: from basic concepts to emerging methods. RoFo Fortschritte Gebiet Rontgenstrahlen Bildgeb Verfahren. (2022) 194:491–504. doi: 10.1055/a-1718-4355
26. Orini M, Seraphim A, Graham A, Bhuva A, Zacur E, Kellman P, et al. Detailed assessment of low-voltage zones localization by cardiac mri in patients with implantable devices. JACC Clin Electrophysiol. (2022) 8:225–35. doi: 10.1016/j.jacep.2021.10.002
27. Stoks J, Cluitmans M. Electrocardiographic imaging: history, applications, and future perspectives. In: W Van Etten, C Swenne, E Van Der Velde editors. IEEE Milestone Award ‘String Galvanometer’: The Heritage and the Promise of Electrocardiography and Electrophysiology. Leiden: Boerhaave Continuing Medical Education (2022). p. 65–84. ISBN: 978 90 67677 82 0.
28. Duchateau J, Sacher F, Pambrun T, Derval N, Chamorro-Servent J, Denis A, et al. Performance and limitations of noninvasive cardiac activation mapping. Heart Rhythm. (2019) 16:435–42. doi: 10.1016/j.hrthm.2018.10.010
29. Erkapic D, Neumann T. Ablation of premature ventricular complexes exclusively guided by three-dimensional noninvasive mapping. Card Electrophysiol Clin. (2015) 7:109–15. doi: 10.1016/j.ccep.2014.11.010
30. Graham A, Orini M, Zacur E, Dhillon G, Daw H, Srinivasan N, et al. Evaluation of ECG imaging to map hemodynamically stable and unstable ventricular arrhythmias. Circ Arrhythm Electrophysiol. (2020) 13:e007377. doi: 10.1161/CIRCEP.119.007377
31. Wellens H, Schuilenburg R, Durrer D. Electrical stimulation of the heart in patients with ventricular tachycardia. Circulation. (1972) 46:216–26. doi: 10.1161/01.CIR.46.2.216
32. Josephson M, Horowitz L, Waxman H, Cain M, Spielman S, Greenspan A, et al. Sustained ventricular tachycardia: role of the 12-lead electrocardiogram in localizing site of origin. Circulation. (1981) 64:257–72. doi: 10.1161/01.CIR.64.2.257
33. Wittkampf F, Wever E, Derksen R, Wilde A, Ramanna H, Hauer R, Robles de medina EO. LocaLisa: new technique for real-time 3-dimensional localization of regular intracardiac electrodes. Circulation. (1999) 99:1312–7. doi: 10.1161/01.CIR.99.10.1312
34. De Groot N, Bootsma M, Van Der Velde E, Schalij M. Three-dimensional catheter positioning during radiofrequency ablation in patients: first application of a real-time position management system. J Cardiovasc Electrophysiol. (2000) 11:1183–92. doi: 10.1046/j.1540-8167.2000.01183.x
35. Marchlinski F, Callans D, Gottlieb C, Zado E. Linear ablation lesions for control of unmappable ventricular tachycardia in patients with ischemic and nonischemic cardiomyopathy. Circulation. (2000) 101:1288–96. doi: 10.1161/01.CIR.101.11.1288
36. Dickfeld T, Calkins H, Bradley D, Solomon S. Stereotactic catheter navigation using magnetic resonance image integration in the human heart. Heart Rhythm. (2005) 2:413–5. doi: 10.1016/j.hrthm.2004.11.023
37. Ghanem R, Jia P, Ramanathan C, Ryu K, Markowitz A, Rudy Y. Noninvasive electrocardiographic imaging (ECGI): comparison to intraoperative mapping in patients. Heart Rhythm. (2005) 2:339–54. doi: 10.1016/j.hrthm.2004.12.022
38. Tian J, Jeudy J, Smith M, Jimenez A, Yin X, Bruce P, et al. Three-dimensional contrast-enhanced multidetector CT for anatomic, dynamic, and perfusion characterization of abnormal myocardium to guide ventricular tachycardia ablations. Circ Arrhythm Electrophysiol. (2010) 3:496–504. doi: 10.1161/CIRCEP.109.889311
39. Takigawa M, Duchateau J, Sacher F, Martin R, Vlachos K, Kitamura T, et al. Are wall thickness channels defined by computed tomography predictive of isthmuses of postinfarction ventricular tachycardia? Heart Rhythm. (2019) 16:1661–8. doi: 10.1016/j.hrthm.2019.06.012
40. Codreanu A, Odille F, Aliot E, Marie P, Magnin-Poull I, Andronache M, et al. Electroanatomic characterization of post-infarct scars. J Am Coll Cardiol. (2008) 52:839–42. doi: 10.1016/j.jacc.2008.05.038
41. Wijnmaalen A, Van Der Geest R, Van Huls Van Taxis CFB, Siebelink HMJ, Kroft LJM, Bax JJ, et al. Head-to-head comparison of contrast-enhanced magnetic resonance imaging and electroanatomical voltage mapping to assess post-infarct scar characteristics in patients with ventricular tachycardias: real-time image integration and reversed registration. Eur Heart J. (2011) 32:104–14. doi: 10.1093/eurheartj/ehq345
42. Cochet H, Komatsu Y, Sacher F, Jadidi A, Scherr D, Riffaud M, et al. Integration of merged delayed-enhanced magnetic resonance imaging and multidetector computed tomography for the guidance of ventricular tachycardia ablation: a pilot study. J Cardiovasc Electrophysiol. (2013) 24:419–26. doi: 10.1111/jce.12052
43. Prakosa A, Arevalo H, Deng D, Boyle P, Nikolov P, Ashikaga H, et al. Personalized virtual-heart technology for guiding the ablation of infarct-related ventricular tachycardia. Nat Biomed Eng. (2018) 2:732–40. doi: 10.1038/s41551-018-0282-2
44. Lozoya R, Berte B, Cochet H, Jaïs P, Ayache N, Sermesant M. Model-Based feature augmentation for cardiac ablation target learning from images. IEEE Trans Biomed Eng. (2019) 66:30–40. doi: 10.1109/TBME.2018.2818300
45. Cuculich P, Schill M, Kashani R, Mutic S, Lang A, Cooper D, et al. Noninvasive cardiac radiation for ablation of ventricular tachycardia. N Engl J Med. (2017) 377:2325–36. doi: 10.1056/nejmoa1613773
46. Grehn M, Mandija S, Andratschke N, Zeppenfeld K, Blamek S, Fast M, et al. Survey results of the STOPSTORM consortium about stereotactic arrhythmia radioablation in Europe. EP Eur. (2022) 24:40849. doi: 10.1093/europace/euac053.376
Keywords: VT, image integration, multi-modality, electroanatomical mapping, ECGI, CMR, personalized medicine, CT
Citation: Stoks J, Hermans BJM, Boukens BJD, Holtackers RJ, Gommers S, Kaya YS, Vernooy K, Cluitmans MJM, Volders PGA and ter Bekke RMA (2023) High-resolution structural-functional substrate-trigger characterization: Future roadmap for catheter ablation of ventricular tachycardia. Front. Cardiovasc. Med. 10:1112980. doi: 10.3389/fcvm.2023.1112980
Received: 30 November 2022; Accepted: 03 February 2023;
Published: 16 February 2023.
Edited by:
Bart Mulder, University Medical Center Groningen, NetherlandsReviewed by:
Moniek Cox, University Medical Center Groningen, NetherlandsAnn-Kathrin Rahm, Heidelberg University Hospital, Germany
Georgios Kollias, Krankenhaus der Elisabethinen, Austria
Copyright © 2023 Stoks, Hermans, Boukens, Holtackers, Gommers, Kaya, Vernooy, Cluitmans, Volders and ter Bekke. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Rachel M. A. ter Bekke, cmFjaGVsLnRlci5iZWtrZUBtdW1jLm5s