Up digital and personal: How heart digital twins can transform heart patient care

Precision medicine is the vision of health care where therapy is tailored to each patient. As part of this vision, digital twinning technology promises to deliver a digital representation of organs or even patients by using tools capable of simulating personal health conditions and predicting patient or disease trajectories on the basis of relationships learned both from data and from biophysics knowledge. Such virtual replicas would update themselves with data from monitoring devices and medical tests and assessments, reflecting dynamically the changes in our health conditions and the responses to treatment. In precision cardiology, the concepts and initial applications of heart digital twins have slowly been gaining popularity and the trust of the clinical community. In this article, we review the advancement in heart digital twinning and its initial translation to the management of heart rhythm disorders.

Predicting ventricular tachycardia circuits in patients with arrhythmogenic right ventricular cardiomyopathy using genotype-specific heart digital twins

Arrhythmogenic right ventricular cardiomyopathy (ARVC) is a genetic cardiac disease that leads to ventricular tachycardia (VT), a life-threatening heart rhythm disorder. Treating ARVC remains challenging due to the complex underlying arrhythmogenic mechanisms, which involve structural and electrophysiological (EP) remodeling. Here, we developed a novel genotype-specific heart digital twin (Geno-DT) approach to investigate the role of pathophysiological remodeling in sustaining VT reentrant circuits and to predict the VT circuits in ARVC patients of different genotypes. This approach integrates the patient's disease-induced structural remodeling reconstructed from contrast-enhanced magnetic-resonance imaging and genotype-specific cellular EP properties. In our retrospective study of 16 ARVC patients with two genotypes: plakophilin-2 (PKP2, n = 8) and gene-elusive (GE, n = 8), we found that Geno-DT accurately and non-invasively predicted the VT circuit locations for both genotypes (with 100%, 94%, 96% sensitivity, specificity, and accuracy for GE patient group, and 86%, 90%, 89% sensitivity, specificity, and accuracy for PKP2 patient group), when compared to VT circuit locations identified during clinical EP studies. Moreover, our results revealed that the underlying VT mechanisms differ among ARVC genotypes. We determined that in GE patients, fibrotic remodeling is the primary contributor to VT circuits, while in PKP2 patients, slowed conduction velocity and altered restitution properties of cardiac tissue, in addition to the structural substrate, are directly responsible for the formation of VT circuits. Our novel Geno-DT approach has the potential to augment therapeutic precision in the clinical setting and lead to more personalized treatment strategies in ARVC.

Predicting Ventricular Tachycardia Circuits in Patients with Arrhythmogenic Right Ventricular Cardiomyopathy using Genotype-specific Heart Digital Twins

Arrhythmogenic right ventricular cardiomyopathy (ARVC) is a genetic cardiac disease that leads to ventricular tachycardia (VT), a life-threatening heart rhythm disorder. Treating ARVC remains challenging due to the complex underlying arrhythmogenic mechanisms, which involve structural and electrophysiological (EP) remodeling. Here, we developed a novel genotype-specific heart digital twin (Geno-DT) approach to investigate the role of pathophysiological remodeling in sustaining VT reentrant circuits and to predict the VT circuits in ARVC patients of different genotypes. This approach integrates the patient's disease-induced structural remodeling reconstructed from contrast-enhanced magnetic-resonance imaging and genotype-specific cellular EP properties. In our retrospective study of 16 ARVC patients with two genotypes: plakophilin-2 (PKP2, n = 8) and gene-elusive (GE, n = 8), we found that Geno-DT accurately and non-invasively predicted the VT circuit locations for both genotypes (with 100%, 94%, 96% sensitivity, specificity, and accuracy for GE patient group, and 86%, 90%, 89% sensitivity, specificity, and accuracy for PKP2 patient group), when compared to VT circuit locations identified during clinical EP studies. Moreover, our results revealed that the underlying VT mechanisms differ among ARVC genotypes. We determined that in GE patients, fibrotic remodeling is the primary contributor to VT circuits, while in PKP2 patients, slowed conduction velocity and altered restitution properties of cardiac tissue, in addition to the structural substrate, are directly responsible for the formation of VT circuits. Our novel Geno-DT approach has the potential to augment therapeutic precision in the clinical setting and lead to more personalized treatment strategies in ARVC.

LASSNet: A Four Steps Deep Neural Network for Left Atrial Segmentation and Scar Quantification

Accurate quantification of left atrium (LA) scar in patients with atrial fibrillation is essential to guide successful ablation strategies. Prior to LA scar quantification, a proper LA cavity segmentation is required to ensure exact location of scar. Both tasks can be extremely time-consuming and are subject to inter-observer disagreements when done manually. We developed and validated a deep neural network to automatically segment the LA cavity and the LA scar. The global architecture uses a multi-network sequential approach in two stages which segment the LA cavity and the LA Scar. Each stage has two steps: a region of interest Neural Network and a refined segmentation network. We analysed the performances of our network according to different parameters and applied data triaging. 200+ late gadolinium enhancement magnetic resonance images were provided by the LAScarQS 2022 Challenge. Finally, we compared our performances for scar quantification to the literature and demonstrated improved performances.

Evaluation of a deep Learning-enabled automated computational heart modeling workflow for personalized assessment of ventricular arrhythmias

Personalized, image-based computational heart modeling is a powerful technology that can be used to improve patient-specific arrhythmia risk stratification and ventricular tachycardia (VT) ablation targeting. However, most state-of-the-art methods still require manual interactions by expert users. The goal of this study is to evaluate the feasibility of an automated, deep learning-based workflow for reconstructing personalized computational electrophysiological heart models to guide patient-specific treatment of VT. Contrast-enhanced computed tomography (CE-CT) images with expert ventricular myocardium segmentations were acquired from 111 patients across 5 cohorts from 3 different institutions. A deep convolutional neural network (CNN) for segmenting left ventricular myocardium from CE-CT was developed, trained, and evaluated. From both CNN-based and expert segmentations in a subset of patients, personalized electrophysiological heart models were reconstructed, and rapid pacing was used to induce VTs. CNN-based and expert segmentations were more concordant in the middle myocardium than in the heart's base or apex. Wavefront propagation during pacing was similar between CNN-based and original heart models. Between most sets of heart models, VT inducibility was the same, the number of induced VTs was strongly correlated, and VT circuits co-localized. Our results demonstrate that personalized computational heart models reconstructed from deep learning-based segmentations even with a small training set size can predict similar VT inducibility and circuit locations as those from expertly-derived heart models. Hence, a user-independent, automated framework for simulating arrhythmias in personalized heart models could feasibly be used in clinical settings to aid VT risk stratification and guide VT ablation therapy. KEY POINTS: Personalized electrophysiological heart modeling can aid in patient-specific ventricular tachycardia (VT) risk stratification and VT ablation targeting. Current state-of-the-art, image-based heart models for VT prediction require expert-dependent, manual interactions that may not be accessible across clinical settings. In this study, we develop an automated, deep learning-based workflow for reconstructing personalized heart models capable of simulating arrhythmias and compare its predictions with that of expert-generated heart models. The number and location of VTs was similar between heart models generated from the deep learning-based workflow and expert-generated heart models. These results demonstrate the feasibility of using an automated computational heart modeling workflow to aid in VT therapeutics and has implications for generalizing personalized computational heart technology to a broad range of clinical centers. Abstract figure legend In this study, we evaluate whether an automated, deep learning-based computational electrophysiological heart models can predict similar arrhythmias as those of expert, manually-derived heart models. First, we build a deep neural network to automatically segment contrast-enhanced CT scans and demonstrate that predicted imaging metrics are comparable to that of manual segmentations. Second, electrophysiological heart models reconstructed from these automated segmentations predict similar wavefront propagation and VT circuits as those of expert-reconstructed heart models. This work represents an advancement towards construction of a user-independent, computational framework to aid in VT risk stratification and guide VT ablation. CT: computed tomography, VT: ventricular tachycardia. This article is protected by copyright. All rights reserved.

Wavefront directionality and decremental stimuli synergistically improve identification of ventricular tachycardia substrate: insights from personalized computational heart models

AIMS

Multiple wavefront pacing (MWP) and decremental pacing (DP) are two electroanatomic mapping (EAM) strategies that have emerged to better characterize the ventricular tachycardia (VT) substrate. The aim of this study was to assess how well MWP, DP, and their combination improve identification of electrophysiological abnormalities on EAM that reflect infarct remodelling and critical VT sites.

METHODS AND RESULTS

Forty-eight personalized computational heart models were reconstructed using images from post-infarct patients undergoing VT ablation. Paced rhythms were simulated by delivering an initial (S1) and an extra-stimulus (S2) from one of 100 locations throughout each heart model. For each pacing, unipolar signals were computed along the myocardial surface to simulate substrate EAM. Six EAM features were extracted and compared with the infarct remodelling and critical VT sites. Concordance of S1 EAM features between different maps was lower in hearts with smaller amounts of remodelling. Incorporating S1 EAM features from multiple maps greatly improved the detection of remodelling, especially in hearts with less remodelling. Adding S2 EAM features from multiple maps decreased the number of maps required to achieve the same detection accuracy. S1 EAM features from multiple maps poorly identified critical VT sites. However, combining S1 and S2 EAM features from multiple maps paced near VT circuits greatly improved identification of critical VT sites.

CONCLUSION

Electroanatomic mapping with MWP is more advantageous for characterization of substrate in hearts with less remodelling. During substrate EAM, MWP and DP should be combined and delivered from locations proximal to a suspected VT circuit to optimize identification of the critical VT site.

Fat infiltration in the infarcted heart as a paradigm for ventricular arrhythmias

Infiltrating adipose tissue (inFAT) has been recently found to co-localize with scar in infarcted hearts and may contribute to ventricular arrhythmias (VAs), a life-threatening heart rhythm disorder. However, the contribution of inFAT to VA has not been well-established. We investigated the role of inFAT versus scar in VA through a combined prospective clinical and mechanistic computational study. Using personalized computational heart models and comparing the results from simulations of VA dynamics with measured electrophysiological abnormalities during the clinical procedure, we demonstrate that inFAT, rather than scar, is a primary driver of arrhythmogenic propensity and is frequently present in critical regions of the VA circuit. We determined that, within the VA circuitry, inFAT, as opposed to scar, is primarily responsible for conduction slowing in critical sites, mechanistically promoting VA. Our findings implicate inFAT as a dominant player in infarct-related VA, challenging existing paradigms and opening the door for unexplored anti-arrhythmic strategies.

Biological substrate modification suppresses ventricular arrhythmias in a porcine model of chronic ischaemic cardiomyopathy

AIMS

Cardiomyopathy patients are prone to ventricular arrhythmias (VA) and sudden cardiac death. Current therapies to prevent VA include radiofrequency ablation to destroy slowly conducting pathways of viable myocardium which support re-entry. Here, we tested the reverse concept, namely that boosting local tissue viability in zones of slow conduction might eliminate slow conduction and suppress VA in ischaemic cardiomyopathy.

METHODS AND RESULTS

Exosomes are extracellular vesicles laden with bioactive cargo. Exosomes secreted by cardiosphere-derived cells (CDCEXO) reduce scar and improve heart function after intramyocardial delivery. In a VA-prone porcine model of ischaemic cardiomyopathy, we injected CDCEXO or vehicle into zones of delayed conduction defined by electroanatomic mapping. Up to 1-month post-injection, CDCEXO, but not the vehicle, decreased myocardial scar, suppressed slowly conducting electrical pathways, and inhibited VA induction by programmed electrical stimulation. In silico reconstruction of electrical activity based on magnetic resonance images accurately reproduced the suppression of VA inducibility by CDCEXO. Strong anti-fibrotic effects of CDCEXO, evident histologically and by proteomic analysis from pig hearts, were confirmed in a co-culture assay of cardiomyocytes and fibroblasts.

CONCLUSION

Biological substrate modification by exosome injection may be worth developing as a non-destructive alternative to conventional ablation for the prevention of recurrent ventricular tachyarrhythmias.

Personalized computational heart models with T1-mapped fibrotic remodeling predict sudden death risk in patients with hypertrophic cardiomyopathy

Hypertrophic cardiomyopathy (HCM) is associated with risk of sudden cardiac death (SCD) due to ventricular arrhythmias (VAs) arising from the proliferation of fibrosis in the heart. Current clinical risk stratification criteria inadequately identify at-risk patients in need of primary prevention of VA. Here, we use mechanistic computational modeling of the heart to analyze how HCM-specific remodeling promotes arrhythmogenesis and to develop a personalized strategy to forecast risk of VAs in these patients. We combine contrast-enhanced cardiac magnetic resonance imaging and T1 mapping data to construct digital replicas of HCM patient hearts that represent the patient-specific distribution of focal and diffuse fibrosis and evaluate the substrate propensity to VA. Our analysis indicates that the presence of diffuse fibrosis, which is rarely assessed in these patients, increases arrhythmogenic propensity. In forecasting future VA events in HCM patients, the imaging-based computational heart approach achieved 84.6%, 76.9%, and 80.1% sensitivity, specificity, and accuracy, respectively, and significantly outperformed current clinical risk predictors. This novel VA risk assessment may have the potential to prevent SCD and help deploy primary prevention appropriately in HCM patients.

Arrhythmia in hypertrophic cardiomyopathy: Risk prediction using contrast enhanced MRI, T1 mapping, and personalized virtual heart technology

BACKGROUND

Hypertrophic cardiomyopathy (HCM), a disease with myocardial fibrosis manifestation, is a common cause of sudden cardiac death (SCD) due to ventricular arrhythmias (VA). Current clinical risk stratification criteria are inadequate in identifying patients who are at risk for VA and in need of an implantable cardioverter defibrillator (ICD) for primary prevention.

OBJECTIVE

We aimed to develop a risk prediction approach based on imaging biomarkers from the combination of late gadolinium contrast-enhanced (LGE) MRI and T1 mapping. We then aimed to compare the prediction to a virtual heart computational risk assessment approach based on LGE-T1 virtual heart models.

METHODS

The methodology involved combining short-axis LGE-MRI with post-contrast T1 maps to define personalized thresholds for diffuse and dense fibrosis. The combined LGE-T1 maps were used to evaluate imaging biomarkers for VA risk prediction. The risk prediction capability of the biomarkers was compared with that of the LGE-T1 virtual heart arrhythmia inducibility simulation. VA risk prediction performance from both approaches was compared to clinical outcome (presence of clinical VA).

RESULTS

Image-based biomarkers, including hypertrophy, signal intensity heterogeneity, and fibrotic border complexity, could not discriminate high vs low VA risk. LGE-T1 virtual heart technology outperformed all the image-based biomarker metrics and was statistically significant in predicting VA risk in HCM.

CONCLUSIONS

We combined two MR imaging techniques to analyze imaging biomarkers in HCM. Raw and processed image-based biomarkers cannot discriminate patients with VA from those without VA. Hybrid LGE-T1 virtual heart models could correctly predict VA risk for this cohort and may improve SCD risk stratification to better identify HCM patients for primary preventative ICD implantation.

Spatial dispersion analysis of LGE-CMR for prediction of ventricular arrhythmias in patients with cardiac sarcoidosis

BACKGROUND

Patients with cardiac sarcoidosis (CS) are at increased risk of life-threatening ventricular arrhythmias (VA). Current approaches to risk stratification have limited predictive value.

OBJECTIVES

To assess the utility of spatial dispersion analysis of late gadolinium enhancement cardiac magnetic resonance (LGE-CMR), as a quantitative measure of myocardial tissue heterogeneity, in risk stratifying patients with CS for VA and death.

METHODS

Sixty two patients with CS underwent LGE-CMR. LGE images were segmented and dispersion maps of the left and right ventricles were generated as follows. Based on signal intensity (SI), each pixel was categorized as abnormal (SI ≥3SD above the mean), intermediate (SI 1-3 SD above the mean) or normal (SI <1SD above the mean); and each pixel was then assigned a value of 0 to 8 based on the number of adjacent pixels of a different category. Average dispersion score was calculated for each patient. The primary endpoint was VA during follow up. The composite of VA or death was assessed as a secondary endpoint.

RESULTS

During 4.7 ± 3.5 years of follow up, six patients had VA, and five without documented VA died. Average dispersion score was significantly higher in patients with VA versus those without (0.87 ± 0.08 vs. 0.71 ± 0.16; p = .002) and in patients with events versus those without (0.83 ± 0.08 vs. 0.70 ± 0.16; p = .003). Patients at higher tertiles of dispersion score had a higher incidence of VA (p = .03) and the composite of VA or death (p = .01).

CONCLUSIONS

Increased substrate heterogeneity, quantified by spatial dispersion analysis of LGE-CMR, may be helpful in risk-stratifying patients with CS for adverse events, including life-threatening arrhythmias.

Analyzing the Role of Repolarization Gradients in Post-infarct Ventricular Tachycardia Dynamics Using Patient-Specific Computational Heart Models

Aims: Disease-induced repolarization heterogeneity in infarcted myocardium contributes to VT arrhythmogenesis but how apicobasal and transmural (AB-TM) repolarization gradients additionally affect post-infarct VT dynamics is unknown. The goal of this study is to assess how AB-TM repolarization gradients impact post-infarct VT dynamics using patient-specific heart models.

Method: 3D late gadolinium-enhanced cardiac magnetic resonance images were acquired from seven post-infarct patients. Models representing the patient-specific scar and infarct border zone distributions were reconstructed without (baseline) and with repolarization gradients along both the AB-TM axes. AB only and TM only models were created to assess the effects of each ventricular gradient on VT dynamics. VTs were induced in all models via rapid pacing.

Results: Ten baseline VTs were induced. VT inducibility in AB-TM models was not significantly different from baseline (p>0.05). Reentry pathways in AB-TM models were different than baseline pathways due to alterations in the location of conduction block (p<0.05). VT exit sites in AB-TM models were different than baseline VT exit sites (p<0.05). VT inducibility of AB only and TM only models were not significantly different than that of baseline (p>0.05) or AB-TM models (p>0.05). Reentry pathways and VT exit sites in AB only and TM only models were different than in baseline (p<0.05). Lastly, repolarization gradients uncovered multiple VT morphologies with different reentrant pathways and exit sites within the same structural, conducting channels.

Conclusion: VT inducibility was not impacted by the addition of AB-TM repolarization gradients, but the VT reentrant pathway and exit sites were greatly affected due to modulation of conduction block. Thus, during ablation procedures, physiological and pharmacological factors that impact the ventricular repolarization gradient might need to be considered when targeting the VTs.

Electromechanical modeling of human ventricles with ischemic cardiomyopathy: numerical simulations in sinus rhythm and under arrhythmia

We developed a novel patient-specific computational model for the numerical simulation of ventricular electromechanics in patients with ischemic cardiomyopathy (ICM). This model reproduces the activity both in sinus rhythm (SR) and in ventricular tachycardia (VT). The presence of scars, grey zones and non-remodeled regions of the myocardium is accounted for by the introduction of a spatially heterogeneous coefficient in the 3D electromechanics model. This 3D electromechanics model is firstly coupled with a 2-element Windkessel afterload model to fit the pressure-volume (PV) loop of a patient-specific left ventricle (LV) with ICM in SR. Then, we employ the coupling with a 0D closed-loop circulation model to analyze a VT circuit over multiple heartbeats on the same LV. We highlight similarities and differences on the solutions obtained by the electrophysiology model and those of the electromechanics model, while considering different scenarios for the circulatory system. We observe that very different parametrizations of the circulation model induce the same hemodynamical considerations for the patient at hand. Specifically, we classify this VT as unstable. We conclude by stressing the importance of combining electrophysiological, mechanical and hemodynamical models to provide relevant clinical indicators in how arrhythmias evolve and can potentially lead to sudden cardiac death.

Characterization of the Electrophysiologic Remodeling of Patients With Ischemic Cardiomyopathy by Clinical Measurements and Computer Simulations Coupled With Machine Learning

Rationale

Patients with ischemic cardiomyopathy (ICMP) are at high risk for malignant arrhythmias, largely due to electrophysiological remodeling of the non-infarcted myocardium. The electrophysiological properties of the non-infarcted myocardium of patients with ICMP remain largely unknown.

Objectives

To assess the pro-arrhythmic behavior of non-infarcted myocardium in ICMP patients and couple computational simulations with machine learning to establish a methodology for the development of disease-specific action potential models based on clinically measured action potential duration restitution (APDR) data.

Methods and Results

We enrolled 22 patients undergoing left-sided ablation (10 ICMP) and compared APDRs between ICMP and structurally normal left ventricles (SNLVs). APDRs were clinically assessed with a decremental pacing protocol. Using genetic algorithms (GAs), we constructed populations of action potential models that incorporate the cohort-specific APDRs. The variability in the populations of ICMP and SNLV models was captured by clustering models based on their similarity using unsupervised machine learning. The pro-arrhythmic potential of ICMP and SNLV models was assessed in cell- and tissue-level simulations. Clinical measurements established that ICMP patients have a steeper APDR slope compared to SNLV (by 38%, p < 0.01). In cell-level simulations, APD alternans were induced in ICMP models at a longer cycle length compared to SNLV models (385–400 vs 355 ms). In tissue-level simulations, ICMP models were more susceptible for sustained functional re-entry compared to SNLV models.

Conclusion

Myocardial remodeling in ICMP patients is manifested as a steeper APDR compared to SNLV, which underlies the greater arrhythmogenic propensity in these patients, as demonstrated by cell- and tissue-level simulations using action potential models developed by GAs from clinical measurements. The methodology presented here captures the uncertainty inherent to GAs model development and provides a blueprint for use in future studies aimed at evaluating electrophysiological remodeling resulting from other cardiac diseases.

Predicting risk of sudden cardiac death in patients with cardiac sarcoidosis using multimodality imaging and personalized heart modeling in a multivariable classifier

Cardiac sarcoidosis (CS), an inflammatory disease characterized by formation of granulomas in the heart, is associated with high risk of sudden cardiac death (SCD) from ventricular arrhythmias. Current "one-size-fits-all" guidelines for SCD risk assessment in CS result in insufficient appropriate primary prevention. Here, we present a two-step precision risk prediction technology for patients with CS. First, a patient's arrhythmogenic propensity arising from heterogeneous CS-induced ventricular remodeling is assessed using a novel personalized magnetic-resonance imaging and positron-emission tomography fusion mechanistic model. The resulting simulations of arrhythmogenesis are fed, together with a set of imaging and clinical biomarkers, into a supervised classifier. In a retrospective study of 45 patients, the technology achieved testing results of 60% sensitivity [95% confidence interval (CI): 57-63%], 72% specificity [95% CI: 70-74%], and 0.754 area under the receiver operating characteristic curve [95% CI: 0.710-0.797]. It outperformed clinical metrics, highlighting its potential to transform CS risk stratification.

Impact of augmented-reality improvement in ablation catheter navigation as assessed by virtual-heart simulations of ventricular tachycardia ablation

BACKGROUND

Recently, an augmented reality (AR) solution allows the physician to place the ablation catheter at the designated lesion site more accurately during cardiac electrophysiology studies. The improvement in navigation accuracy may positively affect ventricular tachycardia (VT) ablation termination, however assessment of this in the clinic would be difficult. Novel personalized virtual heart technology enables non-invasive identification of optimal lesion targets for infarct-related VT. This study aims to evaluate the potential impact of such catheter navigation accuracy improvement in virtual VT ablations.

METHODS

2 MRI-based virtual hearts with 2 in silico induced VTs (VT 1, VT 2) were included. VTs were terminated with virtual "ground truth" endocardial ablation lesions. 106 navigation error values that were previously assessed in a clinical study evaluating the improvement of ablation catheter navigation accuracy guided with AR (53 with, 53 without) were used to displace the "ground truth" ablation targets. The corresponding ablations were simulated based on these errors and VT termination for each simulation was assessed.

RESULTS

In 54 VT 1 ablation simulations, smaller error with AR significantly resulted in more VT termination (25) compared to the error without AR (16) (P < 0.01). In 52 VT 2 ablation simulations, no significant difference was observed from error with (11) and without AR (13) (P = 0.58). The substrate characteristic may impact the effect of improved accuracy to an improved VT termination.

CONCLUSION

Virtual heart shows that the increased catheter navigation accuracy provided by AR guidance can affect the VT termination.

Feasibility study shows concordance between image-based virtual-heart ablation targets and predicted ECG-based arrhythmia exit-sites

INTRODUCTION

We recently developed two noninvasive methodologies to help guide VT ablation: population-derived automated VT exit localization (PAVEL) and virtual-heart arrhythmia ablation targeting (VAAT). We hypothesized that while very different in their nature, limitations, and type of ablation targets (substrate-based vs. clinical VT), the image-based VAAT and the ECG-based PAVEL technologies would be spatially concordant in their predictions.

OBJECTIVE

The objective is to test this hypothesis in ischemic cardiomyopathy patients in a retrospective feasibility study.

METHODS

Four post-infarct patients who underwent LV VT ablation and had pre-procedural LGE-CMRs were enrolled. Virtual hearts with patient-specific scar and border zone identified potential VTs and ablation targets. Patient-specific PAVEL based on a population-derived statistical method localized VT exit sites onto a patient-specific 238-triangle LV endocardial surface.

RESULTS

Ten induced VTs were analyzed and 9-exit sites were localized by PAVEL onto the patient-specific LV endocardial surface. All nine predicted VT exit sites were in the scar border zone defined by voltage mapping and spatially correlated with successful clinical lesions. There were 2.3 ± 1.9 VTs per patient in the models. All five VAAT lesions fell within regions ablated clinically. VAAT targets correlated well with 6 PAVEL-predicted VT exit sites. The distance between the center of the predicted VT-exit-site triangle and nearest corresponding VAAT ablation lesion was 10.7 ± 7.3 mm.

CONCLUSIONS

VAAT targets are concordant with the patient-specific PAVEL-predicted VT exit sites. These findings support investigation into combining these two complementary technologies as a noninvasive, clinical tool for targeting clinically induced VTs and regions likely to harbor potential VTs.

Personalized Digital-Heart Technology for Ventricular Tachycardia Ablation Targeting in Hearts With Infiltrating Adiposity

BACKGROUND

Infiltrating adipose tissue (inFAT) is a newly recognized proarrhythmic substrate for postinfarct ventricular tachycardias (VT) identifiable on contrast-enhanced computed tomography. This study presents novel digital-heart technology that incorporates inFAT from contrast-enhanced computed tomography to noninvasively predict VT ablation targets and assesses the capability of the technology by comparing its predictions with VT ablation procedure data from patients with ischemic cardiomyopathy.

METHODS

Digital-heart models reflecting patient-specific inFAT distributions were reconstructed from contrast-enhanced computed tomography. The digital-heart identification of fat-based ablation targeting (DIFAT) technology evaluated the rapid-pacing-induced VTs in each personalized inFAT-based substrate. DIFAT targets that render the inFAT substrate noninducible to VT, including VTs that arise postablation, were determined. DIFAT predictions were compared with corresponding clinical ablations to assess the capabilities of the technology.

RESULTS

DIFAT was developed and applied retrospectively to 29 ischemic cardiomyopathy patients with contrast-enhanced computed tomography. DIFAT ablation volumes were significantly less than the estimated clinical ablation volumes (1.87±0.35 versus 7.05±0.88 cm³, P<0.0005). DIFAT targets overlapped with clinical ablations in 79% of patients, mostly in the apex (72%) and inferior/inferolateral (74%). In 3 patients, DIFAT targets colocalized with redo ablations delivered years after the index procedure.

CONCLUSIONS

DIFAT is a novel digital-heart technology for individualized VT ablation guidance designed to eliminate VT inducibility following initial ablation. DIFAT predictions colocalized well with clinical ablation locations but provided significantly smaller lesions. DIFAT also predicted VTs targeted in redo procedures years later. As DIFAT uses widely accessible computed tomography, its integration into clinical workflows may augment therapeutic precision and reduce redo procedures.

Accurate Conduction Velocity Maps and Their Association With Scar Distribution on Magnetic Resonance Imaging in Patients With Postinfarction Ventricular Tachycardias

BACKGROUND

Characterizing myocardial conduction velocity (CV) in patients with ischemic cardiomyopathy (ICM) and ventricular tachycardia (VT) is important for understanding the patient-specific proarrhythmic substrate of VTs and therapeutic planning. The objective of this study is to accurately assess the relation between CV and myocardial fibrosis density on late gadolinium-enhanced cardiac magnetic resonance imaging (LGE-CMR) in patients with ICM.

METHODS

We enrolled 6 patients with ICM undergoing VT ablation and 5 with structurally normal left ventricles (controls) undergoing premature ventricular contraction or VT ablation. All patients underwent LGE-CMR and electroanatomic mapping (EAM) in sinus rhythm (2960 electroanatomic mapping points analyzed). We estimated CV from electroanatomic mapping local activation time using the triangulation method that provides an accurate estimate of CV as it accounts for the direction of wavefront propagation. We evaluated the association between LGE-CMR intensity and CV with multilevel linear mixed models.

RESULTS

Median CV in patients with ICM and controls was 0.41 m/s and 0.65 m/s, respectively. In patients with ICM, CV in areas with no visible fibrosis was 0.81 m/s (95% CI, 0.59-1.12 m/s). For each 25% increase in normalized LGE intensity, CV decreased by 1.34-fold (95% CI, 1.25-1.43). Dense scar areas have, on average, 1.97- to 2.66-fold slower CV compared with areas without dense scar. Ablation lesions that terminated VTs were localized in areas of slow conduction on CV maps.

CONCLUSIONS

CV is inversely associated with LGE-CMR fibrosis density in patients with ICM. Noninvasive derivation of CV maps from LGE-CMR is feasible. Integration of noninvasive CV maps with electroanatomic mapping during substrate mapping has the potential to improve procedural planning and outcomes. Visual Overview: A visual overview is available for this article.

Substrate Spatial Complexity Analysis for the Prediction of Ventricular Arrhythmias in Patients With Ischemic Cardiomyopathy

BACKGROUND

Transition zones between healthy myocardium and scar form a spatially complex substrate that may give rise to reentrant ventricular arrhythmias (VAs). We sought to assess the utility of a novel machine learning approach for quantifying 3-dimensional spatial complexity of grayscale patterns on late gadolinium enhanced cardiac magnetic resonance images to predict VAs in patients with ischemic cardiomyopathy.

METHODS

One hundred twenty-two consecutive ischemic cardiomyopathy patients with left ventricular ejection fraction ≤35% without prior history of VAs underwent late gadolinium enhanced cardiac magnetic resonance images. From raw grayscale data, we generated graphs encoding the 3-dimensional geometry of the left ventricle. A novel technique, adapted to these graphs, assessed global regularity of signal intensity patterns using Fourier-like analysis and generated a substrate spatial complexity profile for each patient. A machine learning statistical algorithm was employed to discern which substrate spatial complexity profiles correlated with VA events (appropriate implantable cardioverter-defibrillator firings and arrhythmic sudden cardiac death) at 5 years of follow-up. From the statistical machine learning results, a complexity score ranging from 0 to 1 was calculated for each patient and tested using multivariable Cox regression models.

RESULTS

At 5 years of follow-up, 40 patients had VA events. The machine learning algorithm classified with 81% overall accuracy and correctly classified 86% of those without VAs. Overall negative predictive value was 91%. Average complexity score was significantly higher in patients with VA events versus those without (0.5±0.5 versus 0.1±0.2; P<0.0001) and was independently associated with VA events in a multivariable model (hazard ratio, 1.5 [1.2-2.0]; P=0.002).

CONCLUSIONS

Substrate spatial complexity analysis of late gadolinium enhanced cardiac magnetic resonance images may be helpful in refining VA risk in patients with ischemic cardiomyopathy, particularly to identify low-risk patients who may not benefit from prophylactic implantable cardioverter-defibrillator therapy. Visual Overview: A visual overview is available for this article.

Ventricular arrhythmia risk prediction in repaired Tetralogy of Fallot using personalized computational cardiac models

BACKGROUND

Adults with repaired tetralogy of Fallot (rTOF) are at increased risk for ventricular tachycardia (VT) due to fibrotic remodeling of the myocardium. However, the current clinical guidelines for VT risk stratification and subsequent implantable cardioverter-defibrillator deployment for primary prevention of sudden cardiac death in rTOF remain inadequate.

OBJECTIVE

The purpose of this study was to determine the feasibility of using an rTOF-specific virtual-heart approach to identify patients stratified incorrectly as being at low VT risk by current clinical criteria.

METHODS

This multicenter retrospective pilot study included 7 adult rTOF patients who were considered low risk for VT based on clinical criteria. Patient-specific computational heart models were generated from late gadolinium enhanced magnetic resonance imaging (LGE-MRI), incorporating the individual distribution of rTOF fibrotic remodeling in both ventricles. Simulations of rapid pacing determined VT inducibility. Model creation and simulations were performed by operators blinded to clinical outcome.

RESULTS

Two patients in the study experienced clinical VT. The virtual hearts constructed from LGE-MRI scans of 7 rTOF patients correctly predicted reentrant VT in the models from VT-positive patients and no arrhythmia in those from VT-negative patients. There were no statistically significant differences in clinical criteria commonly used to assess VT risk, including QRS duration and age, between patients who did and those who did not experience clinical VT.

CONCLUSION

This study demonstrates the feasibility of image-based virtual-heart modeling in patients with congenital heart disease and structurally abnormal hearts. It highlights the potential of the methodology to improve VT risk stratification in patients with rTOF.

How personalized heart modeling can help treatment of lethal arrhythmias: A focus on ventricular tachycardia ablation strategies in post-infarction patients

Precision Cardiology is a targeted strategy for cardiovascular disease prevention and treatment that accounts for individual variability. Computational heart modeling is one of the novel approaches that have been developed under the umbrella of Precision Cardiology. Personalized computational modeling of patient hearts has made strides in the development of models that incorporate the individual geometry and structure of the heart as well as other patient-specific information. Of these developments, one of the potentially most impactful is the research aimed at noninvasively predicting the targets of ablation of lethal arrhythmia, ventricular tachycardia (VT), using patient-specific models. The approach has been successfully applied to patients with ischemic cardiomyopathy in proof-of-concept studies. The goal of this paper is to review the strategies for computational VT ablation guidance in ischemic cardiomyopathy patients, from model developments to the intricacies of the actual clinical application. To provide context in describing the road these computational modeling applications have undertaken, we first review the state of the art in VT ablation in the clinic, emphasizing the benefits that personalized computational prediction of ablation targets could bring to the clinical electrophysiology practice. This article is characterized under: Analytical and Computational Methods > Computational Methods Models of Systems Properties and Processes > Organ, Tissue, and Physiological Models Translational, Genomic, and Systems Medicine > Translational Medicine.

Characterizing Conduction Channels in Postinfarction Patients Using a Personalized Virtual Heart

Patients with myocardial infarction have an abundance of conduction channels (CC); however, only a small subset of these CCs sustain ventricular tachycardia (VT). Identifying these critical CCs (CCCs) in the clinic so that they can be targeted by ablation remains a significant challenge. The objective of this study is to use a personalized virtual-heart approach to conduct a three-dimensional (3D) assessment of CCCs sustaining VTs of different morphologies in these patients, to investigate their 3D structural features, and to determine the optimal ablation strategy for each VT. To achieve these goals, ventricular models were constructed from contrast enhanced magnetic resonance imagings of six postinfarction patients. Rapid pacing induced VTs in each model. CCCs that sustained different VT morphologies were identified. CCCs' 3D structure and type and the resulting rotational electrical activity were examined. Ablation was performed at the optimal part of each CCC, aiming to terminate each VT with a minimal lesion size. Predicted ablation locations were compared to clinical. Analyzing the simulation results, we found that the observed VTs in each patient model were sustained by a limited number (2.7 ± 1.2) of CCCs. Further, we identified three types of CCCs sustaining VTs: I-type and T-type channels, with all channel branches bounded by scar, and functional reentry channels, which were fully or partially bounded by conduction block surfaces. The different types of CCCs accounted for 43.8, 18.8, and 37.4% of all CCCs, respectively. The mean narrowest width of CCCs or a branch of CCC was 9.7 ± 3.6 mm. Ablation of the narrowest part of each CCC was sufficient to terminate VT. Our results demonstrate that a personalized virtual-heart approach can determine the possible VT morphologies in each patient and identify the CCCs that sustain reentry. The approach can aid clinicians in identifying accurately the optimal VT ablation targets in postinfarction patients.

Sensitivity of Ablation Targets Prediction to Electrophysiological Parameter Variability in Image-Based Computational Models of Ventricular Tachycardia in Post-infarction Patients

Ventricular tachycardia (VT), which could lead to sudden cardiac death, occurs frequently in patients with myocardial infarction. Computational modeling has emerged as a powerful platform for the non-invasive investigation of lethal heart rhythm disorders in post-infarction patients and for guiding patient VT ablation. However, it remains unclear how VT dynamics and predicted ablation targets are influenced by inter-patient variability in action potential duration (APD) and conduction velocity (CV). The goal of this study was to systematically assess the effect of changes in the electrophysiological parameters on the induced VTs and predicted ablation targets in personalized models of post-infarction hearts. Simulations were conducted in 5 patient-specific left ventricular models reconstructed from late gadolinium-enhanced magnetic resonance imaging scans. We comprehensively characterized all possible pre-ablation and post-ablation VTs in simulations conducted with either an “average human VT”-based electrophysiological representation (i.e., EP_avg) or with ±10% APD or CV (i.e., EP_var); additional simulations were also executed in some models for an extended range of these parameters. The results showed that: (1) a subset of reentries (76.2–100%, depending on EP parameter set) conducted with ±10% APD/CV was observed in approximately the same locations as reentries observed in EP_avg cases; (2) emergent VTs could be induced sometimes after ablation in EP_avg models, and these emergent VTs often corresponded to the pre-ablation reentries in simulations with EP_var parameter sets. These findings demonstrate that the VT ablation target uncertainty in patient-specific ventricular models with an average representation of VT-remodeled electrophysiology is relatively low and the ablation targets stable, as the localization of the induced VTs was primarily driven by the remodeled structural substrate. Thus, personalized ventricular modeling with an average representation of infarct-remodeled electrophysiology may uncover most targets for VT ablation.

Computational Identification of Ventricular Arrhythmia Risk in Pediatric Myocarditis

Children with myocarditis have increased risk of ventricular tachycardia (VT) due to myocardial inflammation and remodeling. There is currently no accepted method for VT risk stratification in this population. We hypothesized that personalized models developed from cardiac late gadolinium enhancement magnetic resonance imaging (LGE-MRI) could determine VT risk in patients with myocarditis using a previously-validated protocol. Personalized three-dimensional computational cardiac models were reconstructed from LGE-MRI scans of 12 patients diagnosed with myocarditis. Four patients with clinical VT and eight patients without VT were included in this retrospective analysis. In each model, we incorporated a personalized spatial distribution of fibrosis and myocardial fiber orientations. Then, VT inducibility was assessed in each model by pacing rapidly from 26 sites distributed throughout both ventricles. Sustained reentrant VT was induced from multiple pacing sites in all patients with clinical VT. In the eight patients without clinical VT, we were unable to induce sustained reentry in our simulations using rapid ventricular pacing. Application of our non-invasive approach in children with myocarditis has the potential to correctly identify those at risk for developing VT.

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 radiofrequency 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 (5 swine) and human studies (21 patients) and in a prospective feasibility study (5 patients). We first assessed in retrospective studies (one of which included a proportion of clinical images with artifacts) 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 prior to the clinical procedure.

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 (RD_sim and RD_FIRM) 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 RD_sim sites. Regions harboring RD_sim and RD_FIRM were compared. RD_sim were found in 38 atrial regions (median [inter-quartile range (IQR)] = 4 [3; 4] per model). RD_FIRM were identified and subsequently ablated in 24 atrial regions (2 [1; 3] per patient), which was significantly fewer than the number of RD_sim-harboring regions in corresponding models (p < 0.05). Computational modeling predicted RD_sim in 20 of 24 (83%) atrial regions identified as RD_FIRM-harboring during clinical mapping. In a large number of cases, we uncovered RD_sim-harboring regions in which RD_FIRM were never observed (18/22 regions that differed between the two modalities; 82%); we termed such cases “latent” RD_sim sites. During follow-up (230 [180; 326] days), AF recurrence occurred in 7/11 (64%) individuals. Interestingly, latent RD_sim sites were observed in all seven computational models corresponding to patients who experienced recurrent AF (2 [2; 2] per patient); in contrast, latent RD_sim 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.

Field of view of mapping catheters quantified by electrogram associations with radius of myocardial attenuation on contrast-enhanced cardiac computed tomography

BACKGROUND

Contrast-enhanced cardiac computed tomography (CE-CT) provides useful substrate characterization in patients with ventricular tachycardia (VT).

OBJECTIVE

The purpose of this study was to describe the association between endocardial electrogram measurements and myocardial characteristics on CE-CT, in particular the field of view of electrogram features.

METHODS

Fifteen patients with postinfarct VT who underwent catheter ablation with preprocedural CE-CT were included. Electroanatomic maps were registered to CE-CT, and myocardial attenuation surrounding each endocardial point was measured at a radius of 5, 10, and 15 mm. The association between endocardial voltage and attenuation was assessed using a multilevel random effects linear regression model, clustered by patient, with best model fit defined by highest log likelihood.

RESULTS

A total of 4698 points were included. There was a significant association of bipolar and unipolar voltage with myocardial attenuation at all radii. For unipolar voltage, the best model fit was at an analysis radius of 15 mm regardless of the mapping catheter used. For bipolar voltage, the best model fit was at an analysis radius of 15 mm for points acquired with a conventional ablation catheter. In contrast, the best model fit for points acquired with a multipolar mapping catheter was at an analysis radius of 5 mm.

CONCLUSION

Myocardial attenuation on CE-CT indicates a smaller myocardial field of view of bipolar electrograms using multipolar catheters with smaller electrodes in comparison to standard ablation catheters despite similar interelectrode spacing. Smaller electrodes may provide improved spatial resolution for the definition of myocardial substrate for VT ablation.

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.

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.

Feasibility of using patient-specific models and the "minimum cut" algorithm to predict optimal ablation targets for left atrial flutter

BACKGROUND

Left atrial flutter (LAFL) occurs in patients after atrial fibrillation ablation. Identification of optimal ablation targets to terminate LAFL remains challenging.

OBJECTIVE

The purpose of this study was to use patient-specific models to simulate LAFL and predict optimal ablation targets using a novel approach based on flow network theory.

METHODS

Late gadolinium-enhanced cardiac magnetic resonance scans from 10 patients with LAFL were used to construct atrial models incorporating fibrosis by investigators blinded to procedural findings. Rapid pacing was applied in silico to induce LAFL. In each LAFL, we represented reentrant wave propagation as an electric flow network and identified the "minimum cut" (MC), which was the smallest amount of tissue that separated the flow into 2 discontinuous components. In silico ablation was applied at MCs, and targets were compared to those that terminated LAFL during catheter ablation.

RESULTS

Patient-specific atrial models were successfully generated from patient scans. LAFL was induced in 7 of 10 models. Ablation of MCs terminated LAFL in 4 models and produced new, slower LAFL morphologies in the other 3. For the latter cases, flow analysis was repeated to identify MCs of emergent LAFLs. Ablation of these MCs terminated emergent LAFLs. The MC-based ablation lesions in simulations were similar in length and location to ablation targets that terminated LAFL during catheter ablation for these 7 patients.

CONCLUSION

Personalized atrial simulations can predict ablation targets for LAFL. These simulations provide a powerful tool for planning ablation procedures and may reduce procedural times and complications.

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).

Cardiac electrophysiological activation pattern estimation from images using a patient-specific database of synthetic image sequences

While abnormal patterns of cardiac electrophysiological activation are at the origin of important cardiovascular diseases (e.g., arrhythmia, asynchrony), the only clinically available method to observe detailed left ventricular endocardial surface activation pattern is through invasive catheter mapping. However, this electrophysiological activation controls the onset of the mechanical contraction; therefore, important information about the electrophysiology could be deduced from the detailed observation of the resulting motion patterns. In this paper, we present the study of this inverse cardiac electrokinematic relationship. The objective is to predict the activation pattern knowing the cardiac motion from the analysis of cardiac image sequences. To achieve this, we propose to create a rich patient-specific database of synthetic time series of the cardiac images using simulations of a personalized cardiac electromechanical model, in order to study this complex relationship between electrical activity and kinematic patterns in the context of this specific patient. We use this database to train a machine-learning algorithm which estimates the depolarization times of each cardiac segment from global and regional kinematic descriptors based on displacements or strains and their derivatives. Finally, we use this learning to estimate the patient’s electrical activation times using the acquired clinical images. Experiments on the inverse electrokinematic learning are demonstrated on synthetic sequences and are evaluated on clinical data with promising results. The error calculated between our prediction and the invasive intracardiac mapping ground truth is relatively small (around 10 ms for ischemic patients and 20 ms for nonischemic patient). This approach suggests the possibility of noninvasive electrophysiological pattern estimation using cardiac motion imaging.

3D strain assessment in ultrasound (Straus): a synthetic comparison of five tracking methodologies

This paper evaluates five 3D ultrasound tracking algorithms regarding their ability to quantify abnormal deformation in timing or amplitude. A synthetic database of B-mode image sequences modeling healthy, ischemic and dyssynchrony cases was generated for that purpose. This database is made publicly available to the community. It combines recent advances in electromechanical and ultrasound modeling. For modeling heart mechanics, the Bestel-Clement-Sorine electromechanical model was applied to a realistic geometry. For ultrasound modeling, we applied a fast simulation technique to produce realistic images on a set of scatterers moving according to the electromechanical simulation result. Tracking and strain accuracies were computed and compared for all evaluated algorithms. For tracking, all methods were estimating myocardial displacements with an error below 1 mm on the ischemic sequences. The introduction of a dilated geometry was found to have a significant impact on accuracy. Regarding strain, all methods were able to recover timing differences between segments, as well as low strain values. On all cases, radial strain was found to have a low accuracy in comparison to longitudinal and circumferential components.

Benchmarking framework for myocardial tracking and deformation algorithms: an open access database

In this paper we present a benchmarking framework for the validation of cardiac motion analysis algorithms. The reported methods are the response to an open challenge that was issued to the medical imaging community through a MICCAI workshop. The database included magnetic resonance (MR) and 3D ultrasound (3DUS) datasets from a dynamic phantom and 15 healthy volunteers. Participants processed 3D tagged MR datasets (3DTAG), cine steady state free precession MR datasets (SSFP) and 3DUS datasets, amounting to 1158 image volumes. Ground-truth for motion tracking was based on 12 landmarks (4 walls at 3 ventricular levels). They were manually tracked by two observers in the 3DTAG data over the whole cardiac cycle, using an in-house application with 4D visualization capabilities. The median of the inter-observer variability was computed for the phantom dataset (0.77 mm) and for the volunteer datasets (0.84 mm). The ground-truth was registered to 3DUS coordinates using a point based similarity transform. Four institutions responded to the challenge by providing motion estimates for the data: Fraunhofer MEVIS (MEVIS), Bremen, Germany; Imperial College London - University College London (IUCL), UK; Universitat Pompeu Fabra (UPF), Barcelona, Spain; Inria-Asclepios project (INRIA), France. Details on the implementation and evaluation of the four methodologies are presented in this manuscript. The manually tracked landmarks were used to evaluate tracking accuracy of all methodologies. For 3DTAG, median values were computed over all time frames for the phantom dataset (MEVIS=1.20mm, IUCL=0.73 mm, UPF=1.10mm, INRIA=1.09 mm) and for the volunteer datasets (MEVIS=1.33 mm, IUCL=1.52 mm, UPF=1.09 mm, INRIA=1.32 mm). For 3DUS, median values were computed at end diastole and end systole for the phantom dataset (MEVIS=4.40 mm, UPF=3.48 mm, INRIA=4.78 mm) and for the volunteer datasets (MEVIS=3.51 mm, UPF=3.71 mm, INRIA=4.07 mm). For SSFP, median values were computed at end diastole and end systole for the phantom dataset(UPF=6.18 mm, INRIA=3.93 mm) and for the volunteer datasets (UPF=3.09 mm, INRIA=4.78 mm). Finally, strain curves were generated and qualitatively compared. Good agreement was found between the different modalities and methodologies, except for radial strain that showed a high variability in cases of lower image quality.

Generation of synthetic but visually realistic time series of cardiac images combining a biophysical model and clinical images

We propose a new approach for the generation of synthetic but visually realistic time series of cardiac images based on an electromechanical model of the heart and real clinical 4-D image sequences. This is achieved by combining three steps. The first step is the simulation of a cardiac motion using an electromechanical model of the heart and the segmentation of the end diastolic image of a cardiac sequence. We use biophysical parameters related to the desired condition of the simulated subject. The second step extracts the cardiac motion from the real sequence using nonrigid image registration. Finally, a synthetic time series of cardiac images corresponding to the simulated motion is generated in the third step by combining the motion estimated by image registration and the simulated one. With this approach, image processing algorithms can be evaluated as we know the ground-truth motion underlying the image sequence. Moreover, databases of visually realistic images of controls and patients can be generated for which the underlying cardiac motion and some biophysical parameters are known. Such databases can open new avenues for machine learning approaches.

Synthetic echocardiographic image sequences for cardiac inverse electro-kinematic learning

In this paper, we propose to create a rich database of synthetic time series of 3D echocardiography (US) images using simulations of a cardiac electromechanical model, in order to study the relationship between electrical disorders and kinematic patterns visible in medical images. From a real 4D sequence, a software pipeline is applied to create several synthetic sequences by combining various steps including motion tracking and segmentation. We use here this synthetic database to train a machine learning algorithm which estimates the depolarization times of each cardiac segment from invariant kinematic descriptors such as local displacements or strains. First experiments on the inverse electrokinematic learning are demonstrated on the synthetic 3D US database and are evaluated on clinical 3D US sequences from two patients with Left Bundle Branch Block.