Heterogeneous myocyte remodelling and spatial heterogeneity of repolarization within the myocardial infarction border zone

Funding Acknowledgements

Type of funding sources: Public grant(s) – National budget only. Main funding source(s): Fund for Scientific Research-Flanders (FWO)

Background

Sudden cardiac death due to ventricular arrhythmias is a major cause of mortality after myocardial infarction (MI). The border zone (BZ) surrounding the infarct is the dominant source of arrhythmias. Here a substrate of heterogeneous repolarization is implicated, which could be due to heterogeneous myocyte remodelling.

Objective

To examine myocyte remodelling within the BZ, in comparison to the remote myocardium, and evaluate the local profile of repolarization of these regions in vivo.

Methods

MI was induced by 120-minute occlusion of the left anterior descending coronary artery followed by reperfusion in 6 domestic pigs. After 4 weeks, magnetic resonance imaging was performed to assess infarct remodelling and local wall thickness. Within 3 days, electro-anatomical mapping was performed. A non-contact recording of a 64-electrode array was translated to 2048 electrograms distributed over the LV and local activation-recovery-interval (ARI) determined by custom software. After recovery (2-4 days), the pigs were sacrificed, and samples collected from the BZ and remote region for RNA analysis and single cardiomyocyte isolation. Cell dimensions were measured and cellular AP duration (APD) was optically recorded using a fluorescent voltage dye, Di-8-Annepps (stimulation at 1Hz, 37°C). Expression and variability of cardiomyocyte hypertrophy biomarkers were extracted from single nuclear RNA sequencing data (10x Genomics).

Results

Cardiomyocyte APD in large population samples (> 100 cells per region in each pig) revealed higher heterogeneity in the BZ than the remote region, quantified as the standard deviation (SD) (BZ: 105.9 ± 17.0ms vs remote: 73.9 ± 8.6ms, P = 0.001). Cellular APD heterogeneity correlated strongly with in vivo local ARI heterogeneity, which demonstrated increased heterogeneity in the BZ (R2 = 0.67, P = 0.002). BZ myocytes were hypertrophied with greater increase in cell width than length, and cellular hypertrophy was more heterogeneous by SD in the BZ (BZ: 12.9 ± 2.4μm vs remote: 8.3 ± 1.1μm, P < 0.001). NPPB transcripts reporting on hypertrophic remodelling were higher in BZ than remote (mean lognorm gene expression, BZ: 0.431 ± 0.014 vs remote: 0.107 ± 0.004, P < 0.001), and showed greater heterogeneity in expression between cells by proportion of hypertrophic (NPPB +ve) cells (BZ: 30.86% vs remote: 8.37%, P < 0.001). Wall thickness variance was higher in the BZ compared to the remote region (anterior BZ: 0.15 ± 0.02mm, septal BZ: 0.16 ± 0.04mm vs remote: 0.04 ± 0.02mm, P < 0.001), contributing to increased heterogeneity of local wall stress in BZ.

Conclusion

Cardiomyocyte remodelling in the BZ is heterogeneous, possibly related to differences in local wall stress, which may contribute to heterogeneous repolarization in vivo and underlie arrhythmia vulnerability within the BZ.

Adrenergic stimulation amplifies the difference in beat-to-beat variability between the scar border zone and remote region

Funding Acknowledgements

Type of funding sources: Public Institution(s). Main funding source(s): KU Leuven BOF-C1 “Blood pressure induced premature ventricular beats as triggers for ventricular arrhythmia in ischemic cardiomyopathy”

Background

Myocardial infarction (MI) results in a dense scar region surrounded by a heterogeneous region of fibrosis and remodeled myocytes called the border zone (BZ). Beta-adrenergic stimulation results in increased beat-to-beat variability of repolarization (BVR) which could increase spatial heterogeneity and arrhythmia vulnerability.

Objective

To examine the effect of adrenergic stimulation on the beat-to-beat variability in the BZ, compared to the remote region, using novel methodology for determining spatially dense activation-repolarization intervals.

Methods

Anterior-septal myocardial infarction (MI) was induced in 10 domestic pigs by 120-minute occlusion of the left anterior descending artery followed by reperfusion. Electro-anatomical mapping was performed after one month. The BZ was defined using contact mapping as the region with bipolar voltage between 0.5 and1.5mV. A non-contact recording of a 64-electrode array was translated to 2048 non-contact electrograms distributed over the LV (EnSite PrecisionTM, St. Jude/Abbott Medical). Electrophysiological recordings were made during baseline and during an isoproterenol (ISO) infusion (incremental doses of 0.01µg/kg until 0.04µg/kg). In each of the 2048 points non-contact electrograms over 25 consecutive beats were processed to determine the BVR using a custom-made algorithm, validated against monophasic action potential recordings.

Results

During baseline conditions the maximal BVR was increased in the BZ compared to the remote region (BZ: 3.28±0.90 ms vs remote: 2.61±0.67 ms, P=0.002). During ISO infusion the maximal BVR was also increased in the BZ (BZ: 3.55±0.74 ms vs remote: 2.21±0.60 ms, P<0.001). During baseline the BZ exhibited a larger spatial variance of BVR than the remote region (BZ: 0.20±0.11 ms2 vs remote: 0.087±0.055 ms2, P=0.002). During ISO infusion the spatial variance of BVR was larger in the BZ (BZ: 0.23±0.12 ms2 vs remote: 0.083±0.056 ms2, P=0.001). The maximal BVR was not significantly different during baseline and ISO in the BZ, nor the remote region (P>0.05). However, the difference of the maximal BVR between BZ and remote regions was significantly increased during ISO (baseline: 0.67±0.48 ms vs ISO: 1.34±0.49ms, P=0.001).

Conclusion

The MI BZ showed increased temporal heterogeneity in repolarization that could serve as functional substrate for re-entry. Adrenergic stimulation amplified this vulnerability by increasing the difference in maximal BVR between BZ and remote regions.

Myocardial stiffness assessed by shear wave elastography relates to pressure-volume loop derived measurements of chamber stiffness

Funding Acknowledgements

Type of funding sources: Public grant(s) – National budget only. Main funding source(s): Research Foundations Flanders

Background

Increased myocardial stiffness is an important cause of diastolic dysfunction. Currently, invasive pressure-volume loop analysis is the gold standard method for the assessment of the left ventricular (LV) chamber stiffness. Its non-invasive assessment in the clinic is cumbersome, requires the combination of several parameters and has limited reliability. Shear wave elastography (SWE) is a novel method that evaluates the propagation of shear waves travelling along the myocardium using high frame rate echocardiography. The propagation speed is directly related to myocardial stiffness. Shear waves can be induced naturally by mitral valve closure (MVC). So far, the in vivo validation of SWE against an invasive gold standard reference method is still lacking.

Purpose

To compare myocardial stiffness assessed by shear wave propagation speed after MVC to invasive pressure-volume loop derived measurements of chamber stiffness.

Methods

Fifteen pigs (31.2 ± 4.1 kg) were included in the study. The instantaneous stiffness of the myocardium was altered by performing the following interventions: 1) preload reduction, 2) afterload increase, 3) preload increase and 4) induction of ischemia/reperfusion (I/R) injury by balloon occlusion of the proximal LAD for 90 min. with subsequent reperfusion of 40 min. To obtain the end-diastolic pressure-volume loop relation (EDPVR), a set of pressure-volume loops was acquired under preload reduction. From the EDPVR, the chamber stiffness constant β and operating chamber stiffness dP/dV were derived. SWE datasets in a parasternal long-axis view were acquired with an experimental ultrasound scanner at an average frame rate of 1304 ± 115 Hz. Shear waves after MVC were visualized on tissue acceleration maps by drawing an M-mode line along the interventricular septum (Figure 1A). The propagation speed was calculated by semi-automatically measuring the spatiotemporal slope.

Results

The chamber stiffness constant β significantly increased after the induction of the I/R injury (0.05 ± 0.01 1/ml vs. 0.09 ± 0.03 1/ml; p < 0.001). The operating chamber stiffness dP/dV decreased by reducing preload and increased by increasing afterload, increasing preload or by inducing an I/R injury (0.50 ± 0.18 mmHg/ml vs. 0.09 ± 0.05 mmHg/ml, 0.67 ± 0.19 mmHg/ml, 0.78 ± 0.35 mmHg/ml and 1.09 ± 0.38 mmHg/ml, respectively; p < 0.01). Likewise, shear wave propagation speed after MVC increased by increasing pre- and afterload (p = 0.001) and by inducing I/R injury (p < 0.001) (Figure 1B). Preload reduction had no significant influence (p = 0.118). Shear wave speed had a strong positive correlation with β (r = 0.63; p < 0.001) (Figure 1C) and dP/dV (r = 0.81; p < 0.001) (Figure 1D).

Conclusions

Shear wave speed after MVC is strongly related to invasive pressure-volume loop derived measures of chamber stiffness. The results of this study indicate the potential of SWE as a novel non-invasive method for the assessment of the instantaneous stiffness of the myocardium. Abstract Figure.

Myocardial stiffness assessed by natural shear wave elastography is related to pressure-volume loop derived parameters

Background

Several cardiovascular disorders are accompanied by a stiffening of the myocardium and may result in diastolic heart failure. The non-invasive assessment of myocardial stiffness could therefore improve the understanding of the pathophysiology and guide treatment. Shear wave elastography (SWE) is a recent technique with tremendous potential for evaluating myocardial stiffness in a non-invasive way. Using high frame rate echocardiography, the propagation speed of shear waves is evaluated, which is directly related to the stiffness of the myocardium. These waves are induced by for instance mitral valve closure (MVC) and propagate throughout the cardiac muscle. However, validation of SWE against an invasive gold standard method is lacking.

Purpose

The aim of this study was to compare echocardiographic shear wave elastography against invasive pressure-volume loops, a gold standard reference method for assessing chamber stiffness.

Methods

In 15 pigs (31.2±4.1 kg) stiffness of the myocardium was acutely changed by inducing ischemia/reperfusion (I/R) injury. For this, the proximal LAD was balloon occluded for 90 minutes with subsequent reperfusion for 40 minutes. Conventional and high frame rate echocardiographic images were acquired simultaneously with pressure-volume loops during baseline conditions and after the induction of the I/R injury. Preload was reduced in order to acquire a set of pressure-volume loops to derive the end-diastolic pressure volume relation (EDPVR). From the EDPVR, the stiffness coefficient β and the operating chamber stiffness dP/dV were obtained. High frame rate echocardiographic datasets of the parasternal long axis view were acquired with an experimental ultrasound scanner (HD-PULSE) at an average frame rate of 1304±115 Hz. Tissue acceleration maps were obtained by drawing an M-mode line along the interventricular septum in order to visualize shear waves after MVC (at end-diastole). The propagation speed was assessed by semi-automatically measuring the slope (Figure A).

Results

I/R injury led to an elevated chamber stiffness constant β (0.09±0.03 1/ml vs. 0.05±0.01 1/ml; p<0.001) and operating chamber stiffness dP/dV (1.09±0.38 mmHg/ml vs. 0.50±0.18 mmHg/ml; p<0.01). Likewise, shear wave speed after MVC increased after the induction of the I/R injury in comparison to baseline (6.1±1.2 m/s vs. 3.2±0.8 m/s; p<0.001). Shear wave speed had a moderate positive correlation with β (r=0.63; p<0.001) (Figure B) and a strong positive correlation with dP/dV (r=0.81; p<0.001) (Figure C).

Conclusion

End-diastolic shear wave speed is strongly related to chamber stiffness, assessed invasively by pressure-volume loops. These results indicate that shear wave propagation speed could be used as a novel non-invasive measurement of the mechanical properties of the ventricle.

Funding Acknowledgement

Type of funding sources: Public grant(s) – National budget only. Main funding source(s): FWO - Research Foundation Flanders

Repolarization heterogeneity within the myocardial infarction border zone correlates with variability of myocyte remodeling

Background

Myocardial infarction (MI) results in a regional scar, with a border zone (BZ) of surviving myocytes interspersed with fibrosis providing an anatomical substrate for re-entry. Heterogeneous repolarization within the BZ may add a functional component aggravating re-entrant arrhythmias.

Purpose

We studied BZ heterogeneity and developed novel methodology for high resolution mapping of local in vivo activation-repolarization intervals (ARI) within the BZ and for studying the relation to cellular action potential (AP) profiles of cells isolated from the BZ.

Methods

Anterior-septal myocardial infarction was induced in 5 domestic pigs by 120-minute occlusion of the left anterior descending artery followed by reperfusion (18.9±4.7% of the left ventricle). After 1-month, electro-anatomical mapping was performed. Contact mapping was used to define the BZ (bipolar voltage 0.5–1.5mV). A non-contact recording of a 64-electrode array was translated to 2048 non-contact electrograms distributed over the LV. The non-contact electrograms were processed to determine the ARIs using a custom-made algorithm, validated against monophasic action potential recordings. After 2–4 days recovery, single cardiomyocytes were enzymatically isolated from the anterior-septal BZs and remote regions. Cardiomyocytes were field stimulated at 1Hz at 37°C and cellular AP duration (APD) was optically recorded (fluorescent voltage-sensitive dye Di-8-Annepps).

Results

In vivo, regional ARIs tended to be longer in the BZs than remote. ARI heterogeneity, quantified as the standard deviation of ARIs in a neighborhood of 1cm radius, was increased in the BZ (anterior BZ: 3.4±1.0 ms, P=0.052, septal BZ: 3.6±1.7 ms, P=0.027 vs remote: 2.0±0.5 ms). Cellular APD was measured in large population samples (>100 cells per region in each pig) and was longer in BZ myocytes compared to the remote region. Cellular APD heterogeneity, measured as the standard deviation within cell population samples pooled by region per animal, was increased in the BZ (anterior BZ: 105.9±17.0 ms, P=0.0010; septal BZ: 98.1±20.8 ms, P=0.0127 vs remote: 73.9±8.6 ms). Cell APD correlated to in vivo ARI (R2=0.34, P=0.021) and cellular heterogeneity correlated strongly with in vivo heterogeneity (R2=0.67, P=0.002).

Conclusion

In the BZ of MI, in vivo regional heterogeneity adds a functional substrate for re-entry that may result from heterogeneous cellular remodeling and increased cell-cell APD variability.

Funding Acknowledgement

Type of funding sources: Public Institution(s). Main funding source(s): KU Leuven BOF-C1 “Blood pressure induced premature ventricular beats as triggers for ventricular arrhythmia in ischemic cardiomyopathy”

Natural shear wave propagation speed is influenced by both changes in myocardial structural properties as well as loading conditions

Funding Acknowledgements

Type of funding sources: Public grant(s) – National budget only. Main funding source(s): Research Foundation - Flanders (FWO)

Background

Shear wave elastography (SWE) is a promising tool for the non-invasive assessment of myocardial stiffness. It is based on the evaluation of the propagation speed of shear waves by high frame rate echocardiography. These waves can be induced by for instance mitral valve closure (MVC) and the speed at which they travel is related to the instantaneous stiffness of the myocardium. Myocardial stiffness is defined by the local slope of the stress-strain relation and can therefore be altered by both changes in structural properties of the myocardium as well as loading conditions.

Purpose

The aim of this study was to investigate how changes in myocardial structural properties as well as loading conditions affect shear wave speed after MVC.

Methods

Until now, 8 pigs (weight: 33.6 ± 5.4 kg) were included. The following interventions were performed: 1) preload was reduced by balloon occlusion of the vena cava inferior, 2) afterload was increased by balloon occlusion of descending aorta, 3) preload was increased by intravenous administration of 500 ml of saline and 4) ischemia/reperfusion injury (I/R injury) was induced in the septal wall by balloon occlusion of the LAD for 90 min. with subsequent reperfusion for 40 min. Echocardiographic and left ventricular pressure recordings were simultaneously obtained during each intervention. Left ventricular parasternal long-axis views were acquired with an experimental high frame rate ultrasound scanner (average frame rate: 1279 ± 148 Hz). Shear waves were visualized on tissue acceleration maps by drawing an M-mode line along the interventricular septum. Shear wave propagation speed after MVC was calculated by assessing the slope of the wave pattern on the tissue acceleration map (Figure A).

Results

The change in left ventricular end-diastolic pressure (LVEDP) and shear wave speed after MVC between baseline and each intervention are shown in Figure B and C, respectively. Preload reduction resulted in significant lower LVEDP compared to baseline (p < 0.01), while the other loading changes did not have a significant effect. Shear wave speed after MVC significantly increased by afterload and preload increase (p < 0.01). I/R injury resulted in increased shear wave speed (p < 0.01) without significantly altering LVEDP. There was a good positive correlation between the change in LVEDP and the change in shear wave speed induced by loading changes (r = 0.76; p < 0.001) (Figure D). However, the correlation became less strong if data of I/R injury was taken into account as well (r = 0.63; p < 0.001).

Conclusion

Our results suggest that SWE is capable to characterize myocardial tissue properties and besides has the potential as a novel method for the estimation of left ventricular filling pressures. However, in the presence of structural changes of the myocardium, care should be taken when estimating filling pressures based on shear wave propagation speed.

Abstract Figure.