Adityo Prakosa

Methodology for image-based reconstruction of ventricular geometry for patient-specific modeling of cardiac electrophysiology.

Patient-specific modeling of ventricular electrophysiology requires an interpolated reconstruction of the 3-dimensional (3D) geometry of the patient ventricles from the low-resolution (Lo-res) clinical images. The goal of this study was to implement a processing pipeline for obtaining the interpolated reconstruction, and thoroughly evaluate the efficacy of this pipeline in comparison with alternative methods. The pipeline implemented here involves contouring the epi- and endocardial boundaries in Lo-res images, interpolating the contours using the variational implicit functions method, and merging the interpolation results to obtain the ventricular reconstruction. Five alternative interpolation methods, namely linear, cubic spline, spherical harmonics, cylindrical harmonics, and shape-based interpolation were implemented for comparison. In the thorough evaluation of the processing pipeline, Hi-res magnetic resonance (MR), computed tomography (CT), and diffusion tensor (DT) MR images from numerous hearts were used. Reconstructions obtained from the Hi-res images were compared with the reconstructions computed by each of the interpolation methods from a sparse sample of the Hi-res contours, which mimicked Lo-res clinical images. Qualitative and quantitative comparison of these ventricular geometry reconstructions showed that the variational implicit functions approach performed better than others. Additionally, the outcomes of electrophysiological simulations (sinus rhythm activation maps and pseudo-ECGs) conducted using models based on the various reconstructions were compared. These electrophysiological simulations demonstrated that our implementation of the variational implicit functions-based method had the best accuracy.

Lack of regional association between atrial late gadolinium enhancement on cardiac magnetic resonance and atrial fibrillation rotors

BACKGROUND The extent of left atrial (LA) late gadolinium enhancement (LGE), as a surrogate for fibrosis, has been associated with atrial fibrillation (AF) recurrence after catheter ablation. Furthermore, there is ex vivo evidence that islands of fibrosis may anchor fibrillatory rotors. OBJECTIVE The purpose of this study was to examine the anatomical association of AF rotors with LA and right atrial (RA) LGE on cardiac magnetic resonance. METHODS The cohort included 9 patients with persistent AF (mean age 61.1 ± 9.7 years) who underwent LGE cardiac magnetic resonance before AF ablation using the focal impulse and rotor modulation system. The extent of LA and RA LGE was quantified globally and in each of the 7 sectors: LA posterior/inferior wall, anterior wall, roof, left and right pulmonary vein antra, and RA lateral and septal regions. The multivariable association of rotor incidence with global and per sector LGE extent was examined using multivariable Bernoulli logistic regression estimated by generalized estimating equations. RESULTS The mean RA and LA volumes were 113.2 ± 37.31 and 143.03 ± 58.25 mL, respectively. The mean RA and LA LGE burden was 17.2% ± 11.0% and 17.4% ± 14.4%, respectively. A total of 18 LA rotors and 9 RA rotors were identified in all patients. No univariable or multivariable association was observed between global or per sector LGE extent and focal impulse and rotor modulation rotor incidence. CONCLUSION In this cohort of patients, there was no association between AF rotor incidence and the global or regional extent of RA and LA LGE.

Accuracy of prediction of infarct-related arrhythmic circuits from image-based models reconstructed from low and high resolution MRI

Identification of optimal ablation sites in hearts with infarct-related ventricular tachycardia (VT) remains difficult to achieve with the current catheter-based mapping techniques. Limitations arise from the ambiguities in determining the reentrant pathways location(s). The goal of this study was to develop experimentally validated, individualized computer models of infarcted swine hearts, reconstructed from high-resolution ex-vivo MRI and to examine the accuracy of the reentrant circuit location prediction when models of the same hearts are instead reconstructed from low clinical-resolution MRI scans. To achieve this goal, we utilized retrospective data obtained from four pigs ~10 weeks post infarction that underwent VT induction via programmed stimulation and epicardial activation mapping via a multielectrode epicardial sock. After the experiment, high-resolution ex-vivo MRI with late gadolinium enhancement was acquired. The Hi-res images were downsampled into two lower resolutions (Med-res and Low-res) in order to replicate image quality obtainable in the clinic. The images were segmented and models were reconstructed from the three image stacks for each pig heart. VT induction similar to what was performed in the experiment was simulated. Results of the reconstructions showed that the geometry of the ventricles including the infarct could be accurately obtained from Med-res and Low-res images. Simulation results demonstrated that induced VTs in the Med-res and Low-res models were located close to those in Hi-res models. Importantly, all models, regardless of image resolution, accurately predicted the VT morphology and circuit location induced in the experiment. These results demonstrate that MRI-based computer models of hearts with ischemic cardiomyopathy could provide a unique opportunity to predict and analyze VT resulting for from specific infarct architecture, and thus may assist in clinical decisions to identify and ablate the reentrant circuit(s).

Image-based reconstruction of three-dimensional myocardial infarct geometry for patient-specific modeling of cardiac electrophysiology.

PURPOSE Accurate three-dimensional (3D) reconstruction of myocardial infarct geometry is crucial to patient-specific modeling of the heart aimed at providing therapeutic guidance in ischemic cardiomyopathy. However, myocardial infarct imaging is clinically performed using two-dimensional (2D) late-gadolinium enhanced cardiac magnetic resonance (LGE-CMR) techniques, and a method to build accurate 3D infarct reconstructions from the 2D LGE-CMR images has been lacking. The purpose of this study was to address this need. METHODS The authors developed a novel methodology to reconstruct 3D infarct geometry from segmented low-resolution (Lo-res) clinical LGE-CMR images. Their methodology employed the so-called logarithm of odds (LogOdds) function to implicitly represent the shape of the infarct in segmented image slices as LogOdds maps. These 2D maps were then interpolated into a 3D image, and the result transformed via the inverse of LogOdds to a binary image representing the 3D infarct geometry. To assess the efficacy of this method, the authors utilized 39 high-resolution (Hi-res) LGE-CMR images, including 36 in vivo acquisitions of human subjects with prior myocardial infarction and 3 ex vivo scans of canine hearts following coronary ligation to induce infarction. The infarct was manually segmented by trained experts in each slice of the Hi-res images, and the segmented data were downsampled to typical clinical resolution. The proposed method was then used to reconstruct 3D infarct geometry from the downsampled images, and the resulting reconstructions were compared with the manually segmented data. The method was extensively evaluated using metrics based on geometry as well as results of electrophysiological simulations of cardiac sinus rhythm and ventricular tachycardia in individual hearts. Several alternative reconstruction techniques were also implemented and compared with the proposed method. RESULTS The accuracy of the LogOdds method in reconstructing 3D infarct geometry, as measured by the Dice similarity coefficient, was 82.10% ± 6.58%, a significantly higher value than those of the alternative reconstruction methods. Among outcomes of electrophysiological simulations with infarct reconstructions generated by various methods, the simulation results corresponding to the LogOdds method showed the smallest deviation from those corresponding to the manual reconstructions, as measured by metrics based on both activation maps and pseudo-ECGs. CONCLUSIONS The authors have developed a novel method for reconstructing 3D infarct geometry from segmented slices of Lo-res clinical 2D LGE-CMR images. This method outperformed alternative approaches in reproducing expert manual 3D reconstructions and in electrophysiological simulations.

A feasibility study of arrhythmia risk prediction in patients with myocardial infarction and preserved ejection fraction

AIM To predict arrhythmia susceptibility in myocardial infarction (MI) patients with left ventricular ejection fraction (LVEF)  >35% using a personalized virtual heart simulation approach. METHODS AND RESULTS A total of four contrast enhanced magnetic resonance imaging (MRI) datasets of patient hearts with MI and average LVEF of 44.0 ± 2.6% were used in this study. Because of the preserved LVEF, the patients were not indicated for implantable cardioverter defibrillator (ICD) insertion. One patient had spontaneous ventricular tachycardia (VT) prior to the MRI scan; the others had no arrhythmic events. Simulations of arrhythmia susceptibility were blind to clinical outcome. Models were constructed from patient MRI images segmented to identify myocardium, grey zone, and scar based on pixel intensity. Grey zone was modelled as having altered electrophysiology. Programmed electrical stimulation (PES) was performed to assess VT inducibility from 19 bi-ventricular sites in each heart model. Simulations successfully predicted arrhythmia risk in all four patients. For the patient with arrhythmic event, in-silico PES resulted in VT induction. Simulations correctly predicted that VT was non-inducible for the three patients with no recorded VT events. CONCLUSIONS Results demonstrate that the personalized virtual heart simulation approach may provide a novel risk stratification modality to non-invasively and effectively identify patients with LVEF  >35% who could benefit from ICD implantation.

Personalized virtual-heart technology for guiding the ablation of infarct-related ventricular tachycardia

Ventricular tachycardia (VT), which can lead to sudden cardiac death, occurs frequently in patients with myocardial infarction. Catheter-based radio-frequency ablation of cardiac tissue has achieved only modest efficacy, owing to the inaccurate identification of ablation targets by current electrical mapping techniques, which can lead to extensive lesions and to a prolonged, poorly tolerated procedure. Here, we show that personalized virtual-heart technology based on cardiac imaging and computational modelling can identify optimal infarct-related VT ablation targets in retrospective animal (five swine) and human studies (21 patients), as well as in a prospective feasibility study (five patients). We first assessed, using retrospective studies (one of which included a proportion of clinical images with artefacts), the capability of the technology to determine the minimum-size ablation targets for eradicating all VTs. In the prospective study, VT sites predicted by the technology were targeted directly, without relying on prior electrical mapping. The approach could improve infarct-related VT ablation guidance, where accurate identification of patient-specific optimal targets could be achieved on a personalized virtual heart before the clinical procedure.A personalized virtual-heart model that determines optimal radio-frequency ablation targets for infarct-related tachycardia is validated in retrospective large-animal and patient studies, and in a prospective study in patients.

The Fibrotic Substrate in Persistent Atrial Fibrillation Patients: Comparison Between Predictions From Computational Modeling and Measurements From Focal Impulse and Rotor Mapping

Focal impulse and rotor mapping (FIRM) involves intracardiac detection and catheter ablation of re-entrant drivers (RDs), some of which may contribute to arrhythmia perpetuation in persistent atrial fibrillation (PsAF). Patient-specific computational models derived from late gadolinium-enhanced magnetic resonance imaging (LGE-MRI) has the potential to non-invasively identify all areas of the fibrotic substrate where RDs could potentially be sustained, including locations where RDs may not manifest during mapped AF episodes. The objective of this study was to carry out multi-modal assessment of the arrhythmogenic propensity of the fibrotic substrate in PsAF patients by comparing locations of RD-harboring regions found in simulations and detected by FIRM (RDsim and RDFIRM) and analyze implications for ablation strategies predicated on targeting RDs. For 11 PsAF patients who underwent pre-procedure LGE-MRI and FIRM-guided ablation, we retrospectively simulated AF in individualized atrial models, with geometry and fibrosis distribution reconstructed from pre-ablation LGE-MRI scans, and identified RDsim sites. Regions harboring RDsim and RDFIRM were compared. RDsim were found in 38 atrial regions (median [inter-quartile range (IQR)] = 4 [3; 4] per model). RDFIRM were identified and subsequently ablated in 24 atrial regions (2 [1; 3] per patient), which was significantly fewer than the number of RDsim-harboring regions in corresponding models (p < 0.05). Computational modeling predicted RDsim in 20 of 24 (83%) atrial regions identified as RDFIRM-harboring during clinical mapping. In a large number of cases, we uncovered RDsim-harboring regions in which RDFIRM were never observed (18/22 regions that differed between the two modalities; 82%); we termed such cases “latent” RDsim sites. During follow-up (230 [180; 326] days), AF recurrence occurred in 7/11 (64%) individuals. Interestingly, latent RDsim sites were observed in all seven computational models corresponding to patients who experienced recurrent AF (2 [2; 2] per patient); in contrast, latent RDsim sites were only discovered in two of four patients who were free from AF during follow-up (0.5 [0; 1.5] per patient; p < 0.05 vs. patients with AF recurrence). We conclude that substrate-based ablation based on computational modeling could improve outcomes.