Parameter subset reduction for imaging-based digital twin generation of patients with left ventricular mechanical discoordination.

BACKGROUND: Integration of a patient's non-invasive imaging data in a digital twin (DT) of the heart can provide valuable insight into the myocardial disease substrates underlying left ventricular (LV) mechanical discoordination. However, when generating a DT, model parameters should be identifiable to obtain robust parameter estimations. In this study, we used the CircAdapt model of the human heart and circulation to find a subset of parameters which were identifiable from LV cavity volume and regional strain measurements of patients with different substrates of left bundle branch block (LBBB) and myocardial infarction (MI). To this end, we included seven patients with heart failure with reduced ejection fraction (HFrEF) and LBBB (study ID: 2018-0863, registration date: 2019-10-07), of which four were non-ischemic (LBBB-only) and three had previous MI (LBBB-MI), and six narrow QRS patients with MI (MI-only) (study ID: NL45241.041.13, registration date: 2013-11-12). Morris screening method (MSM) was applied first to find parameters which were important for LV volume, regional strain, and strain rate indices. Second, this parameter subset was iteratively reduced based on parameter identifiability and reproducibility. Parameter identifiability was based on the diaphony calculated from quasi-Monte Carlo simulations and reproducibility was based on the intraclass correlation coefficient ( ICC ) obtained from repeated parameter estimation using dynamic multi-swarm particle swarm optimization. Goodness-of-fit was defined as the mean squared error ( χ 2 ) of LV myocardial strain, strain rate, and cavity volume. RESULTS: A subset of 270 parameters remained after MSM which produced high-quality DTs of all patients ( χ 2 < 1.6), but minimum parameter reproducibility was poor ( ICC min = 0.01). Iterative reduction yielded a reproducible ( ICC min = 0.83) subset of 75 parameters, including cardiac output, global LV activation duration, regional mechanical activation delay, and regional LV myocardial constitutive properties. This reduced subset produced patient-resembling DTs ( χ 2 < 2.2), while septal-to-lateral wall workload imbalance was higher for the LBBB-only DTs than for the MI-only DTs (p < 0.05). CONCLUSIONS: By applying sensitivity and identifiability analysis, we successfully determined a parameter subset of the CircAdapt model which can be used to generate imaging-based DTs of patients with LV mechanical discoordination. Parameters were reproducibly estimated using particle swarm optimization, and derived LV myocardial work distribution was representative for the patient's underlying disease substrate. This DT technology enables patient-specific substrate characterization and can potentially be used to support clinical decision making.

A Lumped Two-Compartment Model for Simulation of Ventricular Pump and Tissue Mechanics in Ischemic Heart Disease.

Introduction: Computational modeling of cardiac mechanics and hemodynamics in ischemic heart disease (IHD) is important for a better understanding of the complex relations between ischemia-induced heterogeneity of myocardial tissue properties, regional tissue mechanics, and hemodynamic pump function. We validated and applied a lumped two-compartment modeling approach for IHD integrated into the CircAdapt model of the human heart and circulation. Methods: Ischemic contractile dysfunction was simulated by subdividing a left ventricular (LV) wall segment into a hypothetical contractile and noncontractile compartment, and dysfunction severity was determined by the noncontractile volume fraction ( N C V F ). Myocardial stiffness was determined by the zero-passive stress length ( L s 0 , p a s ) and nonlinearity ( k E C M ) of the passive stress-sarcomere length relation of the noncontractile compartment. Simulated end-systolic pressure volume relations (ESPVRs) for 20% acute ischemia were qualitatively compared between a two- and one-compartment simulation, and parameters of the two-compartment model were tuned to previously published canine data of regional myocardial deformation during acute and prolonged ischemia and reperfusion. In six patients with myocardial infarction (MI), the N C V F was automatically estimated using the echocardiographic LV strain and volume measurements obtained acutely and 6 months after MI. Estimated segmental N C V F values at the baseline and 6-month follow-up were compared with percentage late gadolinium enhancement (LGE) at 6-month follow-up. Results: Simulation of 20% of N C V F shifted the ESPVR rightward while moderately reducing the slope, while a one-compartment simulation caused a leftward shift with severe reduction in the slope. Through tuning of the N C V F , L s 0 , p a s , and k E C M , it was found that manipulation of the N C V F alone reproduced the deformation during acute ischemia and reperfusion, while additional manipulations of L s 0 , p a s and k E C M were required to reproduce deformation during prolonged ischemia and reperfusion. Out of all segments with LGE>25% at the follow-up, the majority (68%) had higher estimated N C V F at the baseline than at the follow-up. Furthermore, the baseline N C V F correlated better with percentage LGE than N C V F did at the follow-up. Conclusion: We successfully used a two-compartment model for simulation of the ventricular pump and tissue mechanics in IHD. Patient-specific optimizations using regional myocardial deformation estimated the N C V F in a small cohort of MI patients in the acute and chronic phase after MI, while estimated N C V F values closely approximated the extent of the myocardial scar at the follow-up. In future studies, this approach can facilitate deformation imaging-based estimation of myocardial tissue properties in patients with cardiovascular diseases.

In vivo analysis of the origin and characteristics of gaseous microemboli during catheter-mediated irreversible electroporation.

AIMS: Irreversible electroporation (IRE) ablation is a non-thermal ablation method based on the application of direct current between a multi-electrode catheter and skin electrode. The delivery of current through blood leads to electrolysis. Some studies suggest that gaseous (micro)emboli might be associated with myocardial damage and/or (a)symptomatic cerebral ischaemic events. The aim of this study was to compare the amount of gas generated during IRE ablation and during radiofrequency (RF) ablation. METHODS AND RESULTS: In six 60-75 kg pigs, an extracorporeal femoral shunt was outfitted with a bubble-counter to detect the size and total volume of gas bubbles. Anodal and cathodal 200 J IRE applications were delivered in the left atrium (LA) using a 14-electrode circular catheter. The 30 and 60 s 40 W RF point-by-point ablations were performed. Using transoesophageal echocardiography (TOE), gas formation was visualized. Average gas volumes were 0.6 ± 0.6 and 56.9 ± 19.1 µL (P < 0.01) for each anodal and cathodal IRE application, respectively. Also, qualitative TOE imaging showed significantly less LA bubble contrast with anodal than with cathodal applications. Radiofrequency ablations produced 1.7 ± 2.9 and 6.7 ± 7.4 µL of gas, for 30 and 60 s ablation time, respectively. CONCLUSION: Anodal IRE applications result in significantly less gas formation than both cathodal IRE applications and RF applications. This finding is supported by TOE observations.

Deformation imaging to assess global and regional effects of cardiac regenerative therapy in ischaemic heart disease: A systematic review.

Currently, left ventricular ejection fraction (LVEF) is the most common endpoint in cardiovascular stem cell therapy research. However, this global measure of cardiac function might not be suitable to detect the regional effects sorted by this therapy and is hampered by high operator variability and loading dependency. Deformation imaging might be more accurate in detecting potential regional functional improvements by cardiac regenerative therapy. The aim of this systematic review is to provide a comprehensive overview of current literature on the value of deformation imaging in cardiac regenerative therapy. A systematic review of current literature available on PubMed, Embase, and Cochrane databases was performed regarding both animal and patient studies in which deformation imaging was used to study cardiac cell therapy. After critical appraisal, outcomes regarding study design, type of cell therapy, procedural characteristics, outcome measure, method for measuring strain, and efficacy on both LVEF and deformation parameters were depicted. A total of 30 studies, 15 preclinical and 15 clinical, were included for analysis. Deformation outcomes improved significantly in 14 out of 15 preclinical studies and in 10 out of 15 clinical studies, whereas LVEF improved in 12 and 4 articles, respectively. Study designs and used deformation outcomes varied significantly among the included papers. Six studies found a positive effect on deformation outcomes without LVEF improvement. Hence, deformation imaging seems at least equal, and perhaps superior, to LVEF measurement in the assessment of cardiac regenerative therapy. However, strategies varied substantially and call for a standardized approach.

Lower retention after retrograde coronary venous infusion compared with intracoronary infusion of mesenchymal stromal cells in the infarcted porcine myocardium.

BACKGROUND: Commonly used strategies for cell delivery to the heart are intramyocardial injection and intracoronary (IC) infusion, both having their advantages and disadvantages. Therefore, alternative strategies, such as retrograde coronary venous infusion (RCVI), are explored. The aim of this confirmatory study was to compare cardiac cell retention between RCVI and IC infusion. As a secondary end point, the procedural safety of RCVI is assessed. METHODS: Four weeks after myocardial infarction, 12 pigs were randomised to receive mesenchymal stromal cells, labelled with Indium-111, via RCVI (n=6) or IC infusion (n=6). Four hours after cell administration, nuclear imaging was performed to determine the number of cells retained in the heart both in vivo and ex vivo. Procedure-related safety measures were reported. RESULTS: Cardiac cell retention is significantly lower after RCVI compared with IC infusion (in vivo: RCVI: median 2.89% vs IC: median 13.74%, p=0.002, ex vivo: RCVI: median 2.55% vs IC: median 39.40%, p=0.002). RCVI led to development of pericardial fluid and haematomas on the frontal wall of the heart in three cases. Coronary venous dissection after RCVI was seen in three pigs, of which one also developed pericardial fluid and a haematoma. IC infusion led to no flow in one pig. CONCLUSION: RCVI is significantly less efficient in delivering cells to the heart compared with IC infusion. RCVI led to more procedure-related safety issues than IC infusion, with multiple cases of venous dissection and development of haematomas and pericardial fluid collections.

Leukocyte-Associated Immunoglobulin-like Receptor-1 is regulated in human myocardial infarction but its absence does not affect infarct size in mice.

Heart failure after myocardial infarction (MI) depends on infarct size and adverse left ventricular (LV) remodelling, both influenced by the inflammatory response. Leukocyte-associated immunoglobulin-like receptor 1 (LAIR-1) is an inhibitory receptor of ITAM-dependent cell activation, present on almost all immune cells. We investigated regulation of LAIR-1 leukocyte expression after MI in patients and hypothesized that its absence in a mouse model of MI would increase infarct size and adverse remodelling. In patients, LAIR-1 expression was increased 3 days compared to 6 weeks after MI on circulating monocytes (24.8 ± 5.3 vs. 21.2 ± 5.1 MFI, p = 0.008) and neutrophils (12.9 ± 4.7 vs. 10.6 ± 3.1 MFI, p = 0.046). In WT and LAIR-1-/- mice, infarct size after ischemia-reperfusion injury was comparable (37.0 ± 14.5 in WT vs. 39.4 ± 12.2% of the area at risk in LAIR-1-/-, p = 0.63). Remodelling after permanent left coronary artery ligation did not differ between WT and LAIR-1-/- mice (end-diastolic volume 133.3 ± 19.3 vs. 132.1 ± 27.9 µL, p = 0.91 and end-systolic volume 112.1 ± 22.2 vs. 106.9 ± 33.5 µL, p = 0.68). Similarly, no differences were observed in inflammatory cell influx or fibrosis. In conclusion, LAIR-1 expression on monocytes and neutrophils is increased in the acute phase after MI in patients, but the absence of LAIR-1 in mice does not influence infarct size, inflammation, fibrosis or adverse cardiac remodelling.