Regional diastolic dysfunction in post-infarction heart failure: role of local mechanical load and SERCA expression.

AIMS: Regional heterogeneities in contraction contribute to heart failure with reduced ejection fraction (HFrEF). We aimed to determine whether regional changes in myocardial relaxation similarly contribute to diastolic dysfunction in post-infarction HFrEF, and to elucidate the underlying mechanisms. METHODS AND
RESULTS: Using the magnetic resonance imaging phase-contrast technique, we examined local diastolic function in a rat model of post-infarction HFrEF. In comparison with sham-operated animals, post-infarction HFrEF rats exhibited reduced diastolic strain rate adjacent to the scar, but not in remote regions of the myocardium. Removal of Ca2+ within cardiomyocytes governs relaxation, and we indeed found that Ca2+ transients declined more slowly in cells isolated from the adjacent region. Resting Ca2+ levels in adjacent zone myocytes were also markedly elevated at high pacing rates. Impaired Ca2+ removal was attributed to a reduced rate of Ca2+ sequestration into the sarcoplasmic reticulum (SR), due to decreased local expression of the SR Ca2+ ATPase (SERCA). Wall stress was elevated in the adjacent region. Using ex vivo experiments with loaded papillary muscles, we demonstrated that high mechanical stress is directly linked to SERCA down-regulation and slowing of relaxation. Finally, we confirmed that regional diastolic dysfunction is also present in human HFrEF patients. Using echocardiographic speckle-tracking of patients enrolled in the LEAF trial, we found that in comparison with controls, post-infarction HFrEF subjects exhibited reduced diastolic train rate adjacent to the scar, but not in remote regions of the myocardium.
CONCLUSION: Our data indicate that relaxation varies across the heart in post-infarction HFrEF. Regional diastolic dysfunction in this condition is linked to elevated wall stress adjacent to the infarction, resulting in down-regulation of SERCA, disrupted diastolic Ca2+ handling, and local slowing of relaxation.

1. Introduction

Heart failure may result from either inadequate ejection of blood during systole (systolic dysfunction) or impaired filling during diastole (diastolic dysfunction). Many heart failure patients exhibit significant systolic dysfunction and are classified as having heart failure with reduced ejection fraction (HFrEF). An important cause of systolic dysfunction in HFrEF is regional heterogeneity in contraction across the left ventricle, and in patients with a myocardial infarction contractile dysfunction is particularly prominent adjacent to the scar. While several mechanisms may promote remodelling near the infarct, including neurohumoral and paracrine signalling, we recently demonstrated that elevated wall stress in this region directly triggers disruption of cardiomyocyte t-tubules, Ca2+ homoeostasis, and local contractile function.

In addition to systolic impairments, HFrEF patients also frequently exhibit left ventricular diastolic dysfunction., Importantly, abnormal ventricular filling has been shown to be an independent predictor of mortality in HFrEF, suggesting that diastolic dysfunction is a key pathophysiological mediator of disease progression. As with systolic dysfunction, diastolic impairment in HFrEF has been traced to altered cardiomyocyte Ca2+ homoeostasis. Specifically, reduced expression of the sarcoplasmic reticulum (SR) Ca2+ ATPase 2 (SERCA) and, in some cases, decreased activity of the Na+/Ca2+ exchanger (NCX), lead to slowed removal of cytosolic Ca2+ and slowed, incomplete relaxation. However, while local differences in systolic function are established to be a central feature of HFrEF, it is unknown whether global diastolic dysfunction results from regional changes in myocardial relaxation. We presently investigated diastolic function across the hearts of post-infarction HFrEF rats and patients, and hypothesized that relaxation is particularly impaired adjacent to the scar due to wall stress-induced disruption of diastolic Ca2+ handling.

2. Methods

An expanded methods section is available in the Supplementary material online.

2.1 Animal model

All animal experiments were approved by the Norwegian Animal Research Authority and performed in accordance with the Norwegian Animal Welfare Act and NIH Guidelines (NIH publication No. 85-23, revised 2011). Myocardial infarction was induced by coronary artery ligation in anaesthetized male Wistar rats (endotracheal ventilation with a mixture of 2.5% isoflurane, 97.5% O2). Sham-operated rats served as controls. Post-infarction rats that had progressed into heart failure 6 weeks after surgery were selected based on established echocardiographic criteria (large infarct size and left atrial diameter >5 mm). Mean infarct size was 40.5% of the left ventricular wall (standard deviation 6.6%, range 27.3–51.5%). Other animal characteristics are presented in Table . The current work employed an expanded cohort of sham and HFrEF animals compared with the work previously published by Frisk et al.

2.2 In vivo characterization of HFrEF and sham rats

Rats were characterized in vivo by echocardiography and magnetic resonance imaging (MRI). For echocardiography, anaesthesia was induced by O2 and isoflurane (4.0%), and maintained in freely breathing animals using O2 and isoflurane (1.75%). Doppler inflow recordings were used to calculate left ventricular peak trans-mitral flow (E) and peak trans-mitral deceleration rate.

MRI was performed on a 9.4T preclinical MRI system with hardware dedicated for rat cardiac imaging. Anaesthesia was induced in a chamber with O2/4.0% isoflurane, and adjusted during examination using O2/1.5–2.0% isoflurane to maintain stable anaesthesia. Respiration, electrocardiogram, and body temperature were continuously monitored, and the latter maintained at 37°C using heated air. As described in detail in the Supplementary material online, three short-axis slices (basal, mid-ventricular, and apical) were acquired in each animal using phase-contrast MRI. The myocardium was segmented semiautomatically and the viable myocardium divided into three equally sized regions adjacent, medial, and remote to the infarct (Figure ). In each region, circumferential strain and strain rate were calculated as previously described. MRI-based estimates of local curvature and wall thickness and catheter-based measurements of intraventricular pressure enabled calculation of wall stress (expressed as an integration of diastolic and systolic values). Hearts were then excised from sedated animals, immediately immersed in ice-cold saline, weighed, and used for further experiments.

Slowing of relaxation and elevated wall stress in the adjacent region of HFrEF rats. MRI images illustrate the division of HFrEF rat hearts into three regions (adjacent, medial, and remote) according to their proximity to the infarction (A). Comparison was made to equivalent regions in sham-operated controls. Representative recordings (B) and mean data (C) revealed lower peak circumferential systolic strain adjacent to the infarct in HFrEF. The rate of relaxation was dyssynchronous across the failing heart, as indicated by measurements of circumferential strain rate (D). Specifically, peak diastolic strain rate was only reduced in the adjacent region of HFrEF hearts, and was increased in the remote region (E). Wall stress was calculated across sham and HFrEF hearts based on intraventricular pressure measurements and local ventricular geometry (see Methods section). Measurements are presented as bullseye plots (F) and tracings of wall stress values during the cardiac cycle (G). Mean data demonstrate that while integrated wall stress was high in all regions of HFrEF hearts, values were particularly elevated in the adjacent region (H). (nhearts = 14, 29 in sham, HFrEF). *P < 0.05 calculated by two-way ANOVA with a post hoc Bonferroni t-test. AVc, aortic valve closure; AVo, aortic valve opening; MVo, mitral valve opening.

2.3 Cardiomyocyte Ca2+ imaging

For single cell experiments, left ventricular cardiomyocytes were isolated as described previously. Briefly, excised hearts were perfused in a Langendorff setup with an isolation buffer containing collagenase. The remote, medial, and adjacent regions were then separated and individual cardiomyocytes isolated by a secondary digestion and agitation.

Whole-cell, wide-field Ca2+ transients were recorded in cardiomyocytes pre-incubated with 10 µmol/L fluo 4-AM. Cells were plated on the stage of an inverted microscope and superfused with a Hepes Tyrode’s solution containing 1.8 mmol/L CaCl2 at 37°C, and field-stimulated at a range of frequencies (0.5, 1, 2, 4, and 6 Hz). By rapidly applying 10 mmol/L caffeine, we used mono-exponential fits of the decay of the elicited transient to define the rate constant of Ca2+ extrusion (1/τcaff). The NCX contribution to Ca2+ extrusion was defined as the change in rate constant resulting from application of 5 mmol/L Ni2+, an NCX antagonist, while remaining Ca2+ flux was attributed to slow extrusion pathways. The rate of SERCA-dependent Ca2+ reuptake was calculated by subtracting the extrusion rate constant from the rate constant of Ca2+ decay measured during 1 Hz Ca2+ transients (1/τ1Hz − 1/τcaff).

Diastolic [Ca2+]i was determined in cardiomyocytes loaded with 5 µM fura-2 AM, and these values were employed to calibrate fluo-4 recordings as described previously.

2.4 Modelling

Rat ventricular cardiomyocyte electrophysiology and Ca2+ dynamics were simulated using the model of Gattoni et al., to investigate the consequences of alterations in Ca2+ fluxes observed experimentally. The isolated effects of a 45% increase in NCX conductance or 25% reduction in SERCA conductance on the overall Ca2+ transient were determined by allowing the model to run to steady-state during 6 Hz stimulation.

2.5 Papillary muscle experiments

We examined the effect of mechanical stress on myocardial diastolic function in papillary muscles excised from the left ventricle and mounted in a myobath system as previously described. Muscles were stretched to achieve either normal (≈4 kN/m2) or high diastolic stress (≈15–20 kN/m2), and electrically field-stimulated to develop isometric force during 48 h of culture (Figure ). Only muscles that developed ≥1.5 mN force were included. Muscle core ischaemia was avoided by employing a relatively slow pacing rate (0.5 Hz) and low temperature (22°C); indeed, lactate values (in mmol/kg wet weight) in normal (10.8 ± 2.7) and high stress groups (11.2 ± 3.5) were similar to values in freshly isolated muscles (7.5 ± 0.53, P = NS). After 48 h, muscles subjected to high diastolic stress were returned to a normal stress level (≈4 kN/m2) and allowed to stabilise for 10 min while force was recorded. After completion of the experimental protocol, papillary muscles were snap frozen in liquid nitrogen and stored for molecular analysis.

2.6 Western blotting

Frozen tissue from rat left ventricles was homogenized and protein concentrations quantified as previously described. Primary antibodies for immunoblotting were SERCA2, NCX1, phospholamban, phospholamban Ser16, phospholamban Thr17, Cav1.2, and vinculin. We also used western blotting to assess MEK and ERK signalling pathways, as described in the Supplementary material online. Secondary antibodies were anti-goat IgG HRP-conjugated antibody, anti-rabbit, or anti-mouse IgG HRP-linked whole antibody.

2.7 PCR

Papillary muscle tissue was gently defrosted and mRNA extracted using an RNeasy Mini Kit (Qiagen, Hilden, Germany). qPCR was then performed to quantify SERCA2 and NCX1 gene expression, normalized to the expression of housekeeping gene Ribosomal Protein L4.

2.8 Human echocardiography

We examined local diastolic function in a selection of patients enrolled in the LEAF (LEvosimendan in Acute heart Failure following myocardial infarction) trial (NCT00324766). Briefly, this study randomized patients with percutaneous coronary intervention-treated ST-elevation myocardial infarction complicated with symptomatic acute heart failure to 24 h infusion of levosimendan or placebo. The Regional Ethics Committee approved the study, and it was conducted in accordance with the principles of the Declaration of Helsinki. All patients provided written informed consent. Levosimendan treatment was observed to only transiently augment cardiac function in the initial days following treatment., We presently examined echocardiographic four-chamber apical recordings from study participants at 6 weeks of follow-up, and included those with high acoustic quality of recordings, apical infarctions, and reduced ejection fraction (<40%, Figure ). The viable myocardium of the lateral and septal walls was divided into three equally sized regions (Figure ), and speckle-tracking was performed on kernels in the middle of the remote and adjacent regions (B-mode recordings; frame rate = 68 ± 2.7 frames/s). Peak early diastolic strain rate was calculated as the average of values from the septal and lateral walls. Measurements were compared with regionally matched values from patients in the same study cohort with an ejection fraction ≥50% at 6 weeks of follow-up, and who were without a visible infarction in four-chamber view (Figures ). Trans-mitral and pulmonary vein velocities, left ventricular ejection fraction, and end-diastolic volume were calculated as previously described. Blood pressure was measured with a brachial cuff at the same time point. Patient characteristics are presented in Table .

Figure 6

Selection of LEAF trial participants. HFrEF patients were compared with non-failing ‘control’ individuals (EF ≥50%) who were without a detectable myocardial infarction at 6 weeks of follow-up.

2.9 Statistical analysis

All results are expressed as mean values ± standard error of the mean. Statistical significance was calculated by the Student’s t-test for unpaired data. Nested ANOVA or two-way ANOVA analysis with a post hoc Bonferroni t-test was used where appropriate, as indicated in the figure legends. Statistical significance was defined as P < 0.05.

3. Results

3.1 Regional diastolic dysfunction and elevation of wall stress in HFrEF rats

Regional diastolic function was investigated in an experimental rat model of HFrEF. Animals that developed heart failure 6 weeks following myocardial infarction exhibited characteristic increases in end-diastolic left ventricular pressure and diameter, and marked reduction in global systolic and diastolic function (Table ). MRI-based segmentation of the viable myocardium into regions adjacent, medial, and remote to the infarction (Figure ) revealed marked differences in regional function. Peak circumferential systolic strain was most markedly reduced in the adjacent region, but equivalent to sham values in the remote region (Figure ). Furthermore, we observed that local diastolic dysfunction, as assessed by diastolic strain rate, also occurred only in regions adjacent to the infarction (Figure ). As high wall stress may be a trigger for remodelling, we examined local left ventricular wall stress across HFrEF and sham hearts (wall stress = pressure × local radius/2 × wall thickness). Wall stress was increased in HFrEF hearts, and was particularly elevated adjacent to the myocardial scar tissue due to local flattening of curvature and thinning of the ventricular wall (Figure ).

3.2 Impaired diastolic Ca2+ removal adjacent to the infarction

We next investigated whether variable wall stress and local slowing of relaxation in HFrEF were associated with regional differences in cardiomyocyte Ca2+ handling. Myocytes isolated from the three regions of HFrEF hearts and equivalent regions in sham were field-stimulated to elicit Ca2+ transients across a range of frequencies (Figure ). The magnitude of Ca2+ transients was reduced in the medial and remote region in HFrEF, but maintained in the adjacent region (Figure ). However, the declining phase of the Ca2+ transient was markedly slowed in the adjacent region in HFrEF compared with sham (Figure ). Furthermore, adjacent region myocytes exhibited a significant accumulation of diastolic Ca2+ at high pacing frequencies when time available for Ca2+ removal was limited (Figure ). Elevated end-diastolic Ca2+ limits the extent of cardiomyocyte relaxation, whereas reduced speed of Ca2+ removal slows the rate of relaxation; both alterations are consistent with the diastolic dysfunction observed in vivo adjacent to the infarct. By comparison, maintained diastolic function in the remote region of HFrEF hearts was associated with preserved Ca2+ transient decline and significantly lower diastolic Ca2+ levels (Figure ). Diastolic Ca2+ handling in the medial zone was intermediate between values in the adjacent and remotes zones, mirroring observations of local in vivo diastolic function. For clarity of presentation, statistical analysis of the frequency response of Ca2+ homoeostasis across HFrEF hearts is presented in Supplementary material online, .

Impaired diastolic Ca2+ handling in cardiomyocytes from the adjacent region. Representative recordings of whole-cell, wide-field Ca2+ transients (A, 1 Hz and 6 Hz stimulation, fluo-4), and mean data (B) revealed that, relative to sham, Ca2+ transient magnitude was maintained in cardiomyocytes isolated from the adjacent zone of HFrEF hearts, but reduced in the medial and remote zones. However, HFrEF cells from the adjacent region exhibited significantly slower Ca2+ transient decay (C), and at high-pacing frequencies, a significant elevation of diastolic Ca2+ levels (D). (nhearts = 4, 5 in sham, HFrEF). Comparisons between HFrEF regions are presented in Supplementary material online, . *P < 0.05 vs. sham calculated with nested ANOVA.

3.3 Reduced SERCA expression and function in the adjacent region

Function of the main mediators of cytosolic Ca2+ removal, SERCA, and NCX, was assessed by rapid caffeine application in isolated cells (Figure ). Comparison of the decay kinetics of action potential- and caffeine-elicited Ca2+ transients revealed a reduced rate of SR Ca2+ reuptake in the adjacent region (Figure ). Protein levels of SERCA were also significantly down-regulated in this region (Figure ). Expression of phospholamban, the endogenous SERCA inhibitor, was unaltered (Figure ), and no change in phosphorylation status was observed at either serine-16 or threonine-17 (data not shown). Expression of other key mediators of Ca2+ homoeostasis was similarly unchanged (Figure ). Lowered expression and activity of SERCA are consistent with the slowed Ca2+ removal in the adjacent region, and with elevated resting Ca2+ levels observed at higher pacing frequencies (Figure ), where SERCA activity is known to play an augmenting role in setting diastolic [Ca2+]i.

Reduced SERCA activity in the adjacent region. Representative recordings of Ca2+ transients are illustrated for cardiomyocytes during 1 Hz pacing, followed by rapid application of 10 mmol/L caffeine (A). The magnitude of caffeine-elicited Ca2+ release, an indicator of SR Ca2+ content, tended to be lower in the medial and distal regions (B). Fits of the declining phases of 1 Hz and caffeine transients revealed slowed SR Ca2+ reuptake in HFrEF cardiomyocytes from the adjacent region (C), while the rate of NCX Ca2+ extrusion tended to be higher in the remote zone (D) (nhearts = 4, 5 in sham, HFrEF). (E) Caffeine transients recorded in the presence of 5 mmol/L Ni2+ revealed no alterations in Ca2+ fluxes via slow extrusion pathways (nhearts = 2, 3 in sham, HFrEF). (F) Simulations of Ca2+ transients at physiological frequency (6 Hz) demonstrate the consequences of enhancing NCX activity (blue), as observed in myocytes from the distal zone, or reducing SERCA activity (red), as observed in the adjacent zone. *P < 0.05, calculated by nested ANOVA.

Lower SERCA2 expression in the adjacent region. Consistent with impaired SERCA function in the adjacent region, representative immunoblots (A, vinculin as loading control), and mean data (B) show that SERCA2 expression was significantly reduced in the adjacent region. By contrast, expression of phospholamban (PLB, C), NCX (D), the L-type Ca2+ channel (Cav1.2, E), and the ryanodine receptor (RyR, F) were not markedly altered across HFrEF and sham hearts. Complete blots are presented in Supplementary material online, (nhearts = 6, 6 in sham, HFrEF). *P < 0.05, calculated by two-way ANOVA with a post hoc Bonferroni t-test.

Fits of the decays of caffeine-elicited transients were employed to assess rates of non-SERCA mediated Ca2+ removal, with sensitivity to 5 mmol/L Ni2+ used to differentiate between NCX vs. slow extrusion pathways, including the plasma membrane Ca2+ exchanger and mitochondrial Ca2+ uniporter (see Methods section). Removal of Ca2+ by slow pathways was similar across sham and HFrEF hearts (Figure ). However, while Ca2+ extrusion via NCX was also maintained in the adjacent region of HFrEF hearts, NCX activity tended to be increased in the remote zone (Figure ). As NCX is a key regulator of resting [Ca2+]i, especially at lower stimulation frequencies, the graded NCX activity observed across HFrEF is consistent with observed differences in diastolic Ca2+ levels (Figure ).

Diastolic Ca2+ handling directly impacts systolic [Ca2+]i, as competing NCX and SERCA activities and resting Ca2+ levels importantly determine the extent of SR Ca2+ reuptake. Tendencies for higher NCX activity in the remote zone and lowered diastolic [Ca2+]i in the remote and medial zones were associated with decreased SR Ca2+ content (Figure ) and smaller Ca2+ transients (Figure ) in these regions. This finding was confirmed by mathematical modelling, which reproduced a reduction in both resting [Ca2+]i and Ca2+ transient amplitude when NCX activity was matched to values in the remote region (Figure ). Conversely, no reduction in SR Ca2+ content or release (Figure ) was observed in the adjacent region, despite down-regulation of SERCA in this region. However, mathematical simulations verified that reducing SERCA activity by 25% while maintaining NCX activity only modestly affected Ca2+ transient amplitude, but slowed Ca2+ transient decline and elevated diastolic [Ca2+]i (Figure ), paralleling findings in adjacent zone myocytes.

3.4 High mechanical stress down-regulates SERCA and slows relaxation ex vivo

In HFrEF, high in vivo wall stress adjacent to the infarction (Figure ) coincided with down-regulation of SERCA, disrupted cardiomyocyte Ca2+ handling, and impaired myocardial relaxation. To investigate the causality of these associations, we stretched explanted rat left ventricular papillary muscles to reproduce the elevated wall stress values observed in vivo in the adjacent region, and compared with muscles exposed to normal load (Figure ). Forty-eight hours of exposure to high load conditions triggered a marked reduction in SERCA gene expression (Figure ), and slowing of force decline during isometric electrical stimulation (Figure ). These observations suggest that high wall stress adjacent to the myocardial infarction directly promotes local diastolic dysfunction by signalling reduced SERCA expression and impaired removal of cytosolic Ca2+. Interestingly, Figure  also demonstrates a slowing of tension development in the high stress group, consistent with stress-induced disruption of t-tubule structure and systolic function described in our previous work.

High mechanical stress down-regulates SERCA2 and slows relaxation in ex vivo myocardial tissue. Excised rat papillary muscles were stretched to reproduce high wall stress values comparable to those in the adjacent region of HFrEF hearts, and maintained during 48 h of culture. Comparison was made with muscles subjected to low wall stress conditions approximating those present in the normal heart. Average tension during the protocol is illustrated in A for representative muscles. High stress triggered a significant down-regulation of SERCA2 mRNA (B) (nmuscles = 4 in low group, 3 in high). Returning high stress muscles to normal levels at the completion of the protocol revealed that SERCA2 down-regulation was associated with marked slowing of force decline (representative recordings in C; mean data in D) (nmuscles = 11, 9 in low, high). *P < 0.05 vs. low stress calculated with Student’s t-test.

3.5 Regional diastolic dysfunction in patients with post-infarction HFrEF

To investigate whether regional diastolic dysfunction is also a feature of human post-infarction HFrEF, we examined cardiac function in patients enrolled in the LEAF trial that presented with reduced ejection fraction (HFrEF) 6 weeks following a myocardial infarction (Figure , Table ). Compared with controls, global systolic function was reduced in HFrEF, and speckle-tracking revealed marked differences in local function (Figure ). Specifically, in agreement with previous work,, we observed reduced contraction magnitude (peak systolic longitudinal strain) adjacent to the infarction (Figure ). Although less markedly altered, strain values were also somewhat reduced in the remote region of HFrEF patients, which may reflect an expected rise in filling pressures and a previously established load-dependence of strain measurements. In line with our observations made in the rat model, we observed that local diastolic function varied across the failing human heart. Peak strain rate measurements revealed slowed relaxation adjacent to the infarct, but not at remote sites (Figure ). This spatially dyssynchronous pattern of relaxation in HFrEF patients is expected to impair early ventricular filling, and indicates that measurements of local diastolic myocardial function detect diastolic abnormalities not evident from global parameters (Table ).

HFrEF patients exhibited regional diastolic dysfunction adjacent to the infarction. Four-chamber, apical echocardiographic images of HFrEF patients were examined to divide the viable myocardium into regions adjacent and remote to the myocardial infarction (schematically illustrated in A, right panel). Using speckle-tracking to examine local strain at indicated regions (arrows), comparison was made with regionally matched values in control hearts (A, left panel). Peak strain was reduced in HFrEF hearts, particularly in the adjacent region (B). Representative recordings of strain rate (C) and mean measurements (D) reveal that diastolic strain rate was also reduced adjacent to the myocardial infarction in HFrEF, but similar to control values in the remote region (nhearts = 12, 9 in control, HFrEF). *P < 0.05 with two-way ANOVA with a post hoc Bonferroni t-test.

4. Discussion

While the pathophysiology of HFrEF has traditionally been attributed to weakening of contraction, impaired left ventricular filling is highly prevalent. The present study has provided new insight into the underlying mechanisms. We observed that in a rat model of post-infarction HFrEF, regional diastolic function varied significantly across the left ventricle with slowing of relaxation adjacent to the infarction scar. Local impairment of relaxation was associated with reduction in SERCA expression and activity, and impaired diastolic Ca2+ removal in cardiomyocytes. Ex vivo experiments demonstrated that elevated wall stress, as present in the adjacent region, directly triggers reduction in SERCA expression and diastolic dysfunction. Finally, we confirmed that human patients with post-infarction HFrEF also exhibit a similar pattern of regional diastolic dysfunction. These findings reveal a novel mechanistic link between altered mechanical load, active (Ca2+-dependent) myocardial stiffness, and regional variation in ventricular relaxation in HFrEF.

4.1 Regional diastolic dysfunction in post-infarction HFrEF

It is well established that global systolic and diastolic ventricular performance are depressed in HFrEF, as indicated by functional parameters such as ejection fraction and transmitral filling patterns., However, mechanical function at the global level is in essence dependent on regional myocardial function. Indeed, global systolic dysfunction in HFrEF has been previously linked with differences in the regional timing and magnitude of contraction across the ventricle. In contrast, regional diastolic function has not been widely examined and, to the best of our knowledge, the present study represents the first description of local diastolic impairment in post-infarction HFrEF. While homogenous relaxation across the left ventricle normally ensures a rapid fall in intraventricular pressure during early diastole, we propose that regional disparities in relaxation critically contribute to global diastolic dysfunction in HFrEF. This notion is supported by previous work examining acute ischaemia, where slowing of relaxation in ischaemic regions was observed to decrease the mitral-to-apical pressure gradient, and disrupt the normal pattern of ventricular filling. Similarly, our present observations in HFrEF rats showed depression of global diastole with elevated end-diastolic pressures (Table ), but impairment of local diastolic function only adjacent to the infarction (Figure ). We similarly observed slowed relaxation only adjacent to the infarction of HFrEF patients. Of note, these measurements were compared with patients recovering from acute heart failure, i.e. not healthy individuals, which may explain why global diastolic function was not significantly different between the cohorts. We have presented only longitudinal strain and strain rate measurements in these patients since the adjacent region was often not included in short axis sections, making circumferential strain difficult to calculate. Radial strain, on the other hand, less validly reflects myocardial function since the myofibres are not arranged in this orientation.

4.2 Role of SERCA down-regulation

Our investigations in the HFrEF rat model suggest a key role of reduced SERCA expression in promoting diastolic dysfunction in the adjacent region. While other factors such as the sensitivity of the myofilaments to Ca2+ also affect cellular relaxation, our findings support previous work linking SERCA down-regulation and global diastolic dysfunction in human HFrEF. Furthermore, knock-out of SERCA in mice significantly slows Ca2+ uptake and impairs diastolic function, while restoration of SERCA expression in heart failure improves relaxation. Thus, regional slowing of relaxation in post-infarction HFrEF may be a marker of reduced SERCA expression.

Our data indicate that there is a complex balance of SERCA and NCX activity across HFrEF hearts which critically determines local diastolic function. We observed that SERCA loss adjacent to the infarction was not accompanied by altered NCX expression (Figure ), which appears to parallel previous whole-heart observations in human HFrEF patients with diastolic dysfunction. PCR data from papillary muscles subjected to high wall stress also did not reveal significant changes in NCX at the transcriptional level, although a tendency towards down-regulation was observed (relative mRNA expression = 1.0/0.52 in low-stress/high stress). As we have previously observed marked t-tubule reorganization in the adjacent zone, it may be speculated that resulting displacement of NCX from sites of Ca2+ release contributes to the tendency for somewhat slower NCX-mediated Ca2+ release and diastolic dysfunction in this region (Figure ). By contrast, HFrEF patients with preserved diastolic function are reported to exhibit maintained or even increased NCX levels, with the latter proposed to compensate for reduced SERCA-mediated Ca2+ removal. Interestingly, we did presently observe a strong tendency towards increased NCX activity in the remote region, which was linked to reduced resting Ca2+ levels both experimentally (Figure ) and via mathematical modelling (Figure ). We believe that this represents a compensatory response of the remote myocardium to lower active resting tension. Indeed, local acceleration of re-lengthening in the remote zone was observed in vivo. The mechanisms underlying NCX augmentation are presently unclear, as no change in NCX protein levels was observed. However, other factors not presently investigated such as intracellular sodium homoeostasis, transmembrane potential, and the regulatory protein phospholemman are also important determinants of NCX activity. Of note, while SERCA and NCX activities were observed to vary considerably across HFrEF hearts, we did not observe alterations in slow Ca2+ extrusion pathways, which remain present when SERCA and NCX function are blocked (caffeine and Ni2+ treatment, respectively, Figure ). Thus, although for instance mitochondrial structure and function are suggested to be disrupted in HFrEF, we do not expect that changes in function of the mitochondrial Ca2+ uniporter, or plasmalemmal Ca2+ ATPase, are important contributors to differences in local diastolic function in HFrEF.

It is important to note that strain rate reflects not only active, Ca2+-dependent myocardial relaxation, but also early diastolic load (i.e. filling pressure) and restoring forces. Lengthening load was increased in the adjacent region, as evidenced by elevated diastolic wall stress, which is expected to promote faster re-lengthening in opposition to the effects of slowed Ca2+ removal. Restoring force on the other hand, is difficult to determine due to its opposing and simultaneous action with active cross-bridge tension. In general, however, restoring forces are set by the viscoelastic properties of the myocardium and the degree of systolic strain. Indeed, reductions in systolic strain as observed adjacent to the infarct may contribute to decreases in diastolic re-lengthening in this region. Likewise, diastolic dysfunction may be promoted by collagen accumulation which is reported to be particularly prominent adjacent to an infarction., Ventricular sections from HFrEF did indeed show marked fibrosis in the adjacent region (Supplementary material online, ), in agreement with previous work., Tethering to the immobile infarct may also act to reduce strain and strain rate in the adjacent region, although strain measurements based on speckle-tracking are less influenced by tethering effects than Doppler tissue imaging. Furthermore, no correlation was observed between infarct size and diastolic strain rate (Supplementary material online, ). Nevertheless, while we believe that the diastolic deficit in the adjacent region is consistent with greater active myocardial stiffness following SERCA down-regulation, complex alterations in passive stiffness and ventricular geometry may additionally contribute to local slowing of re-lengthening.

Beyond diastolic dysfunction, decreased SERCA expression has often been linked to declining systolic performance in HFrEF, particularly in humans and large animal models, where SERCA loss is associated with reduced SR Ca2+ content and release. Presently, however, we observed that the reduction in SERCA activity was not accompanied by any reduction in SR Ca2+ load (Figure ), Ca2+ transient magnitude (Figure ), or fractional SR Ca2+ release (Supplementary material online, ). This is in line with previous results in small rodents, showing that Ca2+ transient magnitude is less markedly affected by SERCA reduction during heart failure, since action potential prolongation in these species augments triggering of Ca2+ release by L-type Ca2+ current. In addition, the pronounced rise in diastolic Ca2+ observed at high pacing frequencies in the adjacent HFrEF region makes more Ca2+ available for SR reuptake, and thus release. Indeed, mathematical modelling revealed that transients simulated at 6 Hz with a 25% reduction in SERCA activity, only exhibited a modest reduction Ca2+ transient amplitude (Figure ). Conversely, in the remote region, the tendency towards higher NCX activity makes less Ca2+ available for SR reuptake, resulting in reduced SR Ca2+ content and smaller Ca2+ transients. Thus, while augmented NCX activity has compensatory effects on diastolic Ca2+ handling, such actions come at the expense of systolic Ca2+ release. The interplay between systolic and diastolic Ca2+ handling is also complicated by the fact that Ca2+ release and removal fluxes in fact overlap near the peak of the Ca2+ transient. We have previously shown that Ca2+ release is slowed and dyssynchronous in the adjacent region, which is expected to delay and desynchronize the early phase of Ca2+ transient decline. However, due to the short time span of Ca2+ release compared with removal, we expect that such effects would be minor at late stages of the Ca2+ transient where we have estimated Ca2+ fluxes.

Other detrimental consequences of SERCA loss to be considered include increased susceptibility to pro-arrhythmic Ca2+ waves, impaired mechanoenergetics, and triggering of hypertrophy, SR/ER stress, and apoptotic signalling pathways. As both remodelling and arrhythmogenesis are particularly pronounced in regions adjacent to an infarction, a local reduction in SERCA expression can be speculated to locally contribute to these processes.

4.3 Role of elevated wall stress

What drives local remodelling adjacent to a myocardial infarction? Paracrine factors and neurohumoral activation likely play important roles. However, another potent trigger of remodelling—ventricular load—is also strikingly altered in the adjacent region, as tapering of the myocardium and curvature flattening toward the infarction result in locally high wall stress., We recently demonstrated that elevated wall stress triggers disorganization of cardiomyocyte t-tubule structure, leading to slowed, de-synchronized Ca2+ release; changes which are linked to reduced contractility in the adjacent region (Figure ). We presently demonstrate that heightened wall stress also down-regulates SERCA, further supporting the paradigm that excessive mechanical load plays a central role in disrupting myocardial structure and function during HFrEF. Intriguingly, in concentric hypertrophy, which acts to relieve wall stress, we recently observed that SERCA activity was enhanced and diastolic Ca2+ homoeostasis compensated. Diastolic dysfunction in these hearts instead resulted from enhanced passive stiffness. In apparent agreement with wall-stress regulated expression of SERCA, it was previously observed that the expression of SERCA in dyssynchronous heart failure was down-regulated in the late-activated lateral wall which is subjected to excessive strain and stress, whereas unloading of the heart by cardiac resynchronization therapy has been shown to restore SERCA expression. Unloading of the failing heart by left ventricular assist devices also induces reverse remodelling, and in some cases up-regulation of SERCA, providing further support for this paradigm. Thus, while it was previously recognized that high wall stress, particularly in the form of late systolic afterload, acutely slows relaxation, our data indicate that heightened stress also has chronic effects on myocardial relaxation through remodelling of diastolic Ca2+ homoeostasis. In an effort to identify mechanosensitive signalling pathways which may trigger SERCA down-regulation, we investigated previously proposed regulators of SERCA expression, namely the ERK/MEK-pathway, mir25, and TNFα signalling. While components of all three pathways tended to be up-regulated compared with sham, there was no evidence for marked specific up-regulation in the adjacent region compared with the more remote regions (Supplementary material online, ). Thus, although the mechanosensitive mechanism leading to local SERCA loss remains unclear, an intriguing possibility is that locally distinct Ca2+ regulation adjacent to the infarction may contribute to such signalling.

5. Conclusion

In conclusion, we have demonstrated that regional variation in relaxation is associated with diastolic dysfunction in post-infarction HFrEF. Our data indicate that slowing of relaxation adjacent to the infarction is caused by elevated local wall stress, which triggers down-regulation of SERCA expression and slower diastolic Ca2+ removal. These findings establish a new paradigm implicating mechanical stress and cardiomyocyte Ca2+ homoeostasis as local determinants of diastolic function in health and disease.

Click here for additional data file.

Table 1

Animal characteristics

	Sham	HFrEF
Post-mortem (n = 14, 21)
Heart to body weight ratio (mg/g)	3.10 ± 0.13	5.51 ± 1.21*
Lung to body weight ratio (mg/g)	3.37 ± 0.17	10.55 ± 0.53*
Echocardiography (n = 13, 21)
Heart rate (min−1)	409.8 ± 9.7	362.1 ± 7.5*
Peak trans-mitral flow (E) (mm/s)	894.7 ± 16.1	1085.7 ± 54.8*
Peak trans-mitral deceleration rate (cm/s2)	2968 ± 277	5438 ± 600*
MRI (n = 14, 22)
Ejection fraction (%)	69.6 ± 1.53	28.3 ± 8.7*
End-diastolic volume (µL)	467.0 ± 13.9	1052 ± 33*
End-systolic volume (µL)	142.6 ± 8.6	755.6 ± 26.3*
Left ventricular mass (g)	0.70 ± 0.02	0.77 ± 0.02*
Left ventricular catheterization (n= 14, 22)
End-diastolic pressure (mmHg)	2.14 ± 0.27	26.5 ± 0.8*
dP/dtmin (mmHg/ms)	30 792 ± 1524	18 595 ± 551*

Rats which developed HFrEF 6 weeks following myocardial infarction were functionally compared with sham-operated controls.

P < 0.05 vs. sham measured by Student’s t-test.

Table 2

Patient characteristics

	Non-failing (n = 12)	HFrEF (n = 9)
Age (years)	64 ± 3.7	62 ± 3.7
Heart rate (min−1)	62 ± 2.5	77 ± 4.2*
End-diastolic volume (mL)	90 ± 12.1	143 ± 11.6*
End-systolic volume (mL)	37 ± 5.0	95 ± 9.1*
Left ventricular ejection fraction (%)	59 ± 1.4	34 ± 1.6*
E/A	1.5 ± 0.20	1.1 ± 0.17
E/eʹ	10.4 ± 1.0	11.9 ± 1.2
S/D	1.14 ± 0.07	1.07 ± 0.13
Diastolic blood pressure (mmHg)	78 ± 2.4	78 ± 3.1
Systolic blood pressure (mmHg)	123 ± 6.5	116 ± 4.8
Infarct size (%)	32 ± 6.0	54 ± 2.6

Participants in the LEAF trial exhibiting HFrEF were compared with non-failing individuals at 6 weeks of follow-up (see Fig. 6 for inclusion criteria).

E/A, mitral valve E-wave velocity/A-wave velocity, E/eʹ, mitral valve E velocity/mitral annular e′ velocity; S/D, pulmonary vein S-wave/D-wave velocity.

P < 0.05 vs. control measured by Student’s t-test.