Introduction

Follow up of patients after acute myocarditis is important in order to detect persisting or developing heart failure (HF) as well as enabling a safe return to physical activity [1]. In order to identify patients at risk and optimize treatment, current ESC guidelines recommend follow-up by cardiovascular magnetic resonance (CMR) imaging if the initial CMR showed acute inflammation [2]. Besides inconspicuous troponin and ECG, current guidelines demand normalized left ventricular ejection fraction (LVEF) as prerequisite for return to physical activity/sport [2]. Beyond LVEF, current guidelines advocate global longitudinal strain (GLS) as an additional parameter with incremental value to quantify LV-dysfunction [3]. Recent studies showed the capability of left ventricular global longitudinal strain (LV-GLS) to predict outcome in patients with chronic heart failure (HF) due to ischemic and non-ischemic dilated cardiomyopathy as well as in patients with myocarditis [4,5,6]. CMR strain analysis was recently applied in acute myocarditis and revealed impaired LV, left atrial (LA) and right ventricular (RV) systolic function [7,8,9,10,11]. In particular, Luetkens et al. recently suggested CMR derived longitudinal LV-strain as new parameter to predict functional recovery after acute myocarditis and further demonstrated improved LV- and RV-strain at the time of follow-up [10].

Besides ventricular function, the relevance of atrial function is increasingly recognized as an important factor in patients with HF: Left atrial (LA) dilatation represents a strong predictive marker in patients with dilatative cardiomyopathy (DCM) [12]. A recent meta-analysis of studies investigating heart failure with preserved ejection fraction (HFpEF) revealed that reduced LA-reservoir strain precedes changes in LVEF [13]. Furthermore, decreased LA reservoir strain was found to be of prognostic value in HFpEF, but also in heart failure with reduced ejection fraction (HFrEF) [13, 14]. The diagnostic value of LA-strain parameters has been evaluated in acute myocarditis before, but there are currently no data on atrial strain after acute myocarditis [8, 9]. Right atrial (RA) function has been neglected so far. However, CMR feature tracking (FT) recently revealed a prognostic role of RA function in HF patients [15]. In patients with acute myocarditis, Dick et al. recently showed a trend towards reduced RA reservoir strain values [8]. In summary, atrial function is of increasingly recognized value in several clinical settings, but little is known on development of atrial function after acute myocarditis. Therefore, this study evaluated CMR-FT derived myocardial strain of all four cardiac chambers in order to assess atrial and ventricular function within the first 3 months after acute myocarditis.

Methods

Study population

The local ethics committee approved the study (PV3987) and written informed consent was obtained from all participants. This study is based on additional analyses in a subgroup of 35 consecutive patients of a recently published study population with strictly defined acute/active myocarditis, who underwent CMR at baseline (BL) and at 3 months follow-up (FU). We have published a detailed description of the entire study population before [11]. Briefly, 78 patients with clinically suspected myocarditis were prospectively included. Besides a structured interview, cardiac evaluation included the collection of blood samples as well as the performance of CMR, which was repeated after 3 months at FU. Of the 78 included patients with clinically suspected myocarditis, 48 fulfilled a strict definition of acute myocarditis: Acute myocarditis was defined by endomyocardial biopsy (EMB) in patients with “cardiomyopathy-like” (CLM) presentation and by the combination of recent onset of chest pain, dynamically elevated troponin T values and typical, focal, non-ischemic LGE patterns by CMR in patients who presented with “infarct-like” myocarditis (ILM) [11, 16]. From 39 patients, who underwent 3-months FU, 4 patients had to be excluded from this analysis due to insufficient image quality to perform atrial strain measurements of the left and right atrium thus leaving 35 patients in the final study population for this study [17]. Medical treatment of patients was determined independent from participation in this study in agreement with current recommendations [1]. Thirty patients (86%) were treated with beta-blockers, 28 patients (80%) with angiotensin converting enzyme inhibitors/angiotensin II antagonists, 11 patients (31%) with aldosterone antagonists and 11 patients (31%) with loop diuretics. Two patients (6%) received immunosuppressive therapy guided by endomyocardial biopsy.

CMR protocol and data analyses

The CMR protocol was conducted as previously described [11, 18]. Briefly, CMR was performed at 1.5-T (Achieva, Philips Medical Systems, Best, The Netherlands). Besides conventional cine SSFP sequences, the protocol included edema-sensitive, early myocardial enhancement and late gadolinium enhancement (LGE) sequences as well as T1 and T2 mapping [11]. Two patients of our study population had a history of atrial fibrillation, but none of the patients had prevalent atrial fibrillation at the time of CMR acquisition. In addition to our recently published analyses, CMR-FT measurements were performed using dedicated software Medis Suite MR (Medis Medical Imaging, Leiden, The Netherlands) as described before [19]. Myocardial LV-Strain parameters were measured in short- as well as long-axes orientations, employing endocardial as well as epicardial contours according to current standard [20, 21]. RV and RA strain were measured in the 4-chamber view, LA strain in the 2- and 4-chamber view respectively. Biplane left atrial GLS (LABiGLS) and volumes were calculated as the mean value of two and four-chamber measurements. In the following paragraphs, data referred to as “LA” represent biplane values, if not further specified. LV global myocardial strain was assessed as longitudinal (GLS), circumferential (GCS) as well as radial strain (GRS) employing endo- as well as epicardial contours. GLS for the thin walls of LA, RA and RV was obtained from endocardial contours, since epicardial contours are only available for LV strain by the Medis software (Medis Medical Imaging, Leiden, The Netherlands) (Fig. 1). Presented strain parameters represent peak systolic strain values. The maximum of LA/RA was defined at LV end-systole whereas the minimum was defined at LV end-diastole. Due to the superior reproducibility, reported atrial strain values represent total atrial strain (εS) and therefore describe atrial reservoir function [22].

Fig. 1
figure 1

Exemplary strain curves for peak systolic atrial strain. CMR-FT of peak-systolic atrial strain at baseline (A) and follow-up (B). Strain contours are displayed at end-diastole (D) and end-systole (S). Myocardial deformation is visualized through the red dots following the green lines. In this patient RA-GLS improved from 25.8 at BL to 41.6% at FU. LA-GLS (combined from two-chamber and four-chamber LA-GLS) increased from 26.4 at BL to 56.6% at FU. Please note the different scale of displayed graphs, automatically generated by the software

Echocardiography

In order to classify diastolic dysfunction for this sub-cohort-study, echocardiographic data, derived during the patients’ hospitalization, was utilized where available. The classification of diastolic dysfunction was performed as recommended by Nagueh et al. [23]. Unfortunately tissue-tracking was not included in the echo protocol.

Statistical analysis

Statistical analyses were performed using MedCalc for Windows, version 13.3.3.0 (MedCalc Software, Ostend, Belgium). Continuous data are presented as median with first and third quartiles (Q1–Q3). BL and FU parameters were compared using Wilcoxon Test (paired samples) test. Diastolic dysfunction of ILM and CLM was compared by Fisher’s exact test. Twenty randomly selected data sets were analyzed independently by two experienced observers and inter-observer agreement was assessed by calculating intra-class-correlation coefficients (ICC). Two-sided P-values were calculated and P < 0.05 defined as statistically significant. We did not apply a correction for multiple testing. P < 0.05 was defined as statistically significant.

Results

Inter-observer agreement of strain measurements

LA-GLS provided an excellent intraclass correlation coefficient (ICC) with 0.943 (95% CI 0.863; 0.977). RA-GLS also offered a very good ICC with 0.878 (95% CI 0.717; 0.950). ICCs were excellent as well for LV-strain analysis with LVSAXGCS 0.958 (95% CI 0.564; 0.989) and LVLAXGLS 0.922 (95%-CI 0.434; 0.979). Modest ICCs were found for LV-GRS with LVSAXGRS 0.575 (0.077; 0.824) and LVLAXGRS 0.686 (0.366; 0.862).

Atrial recovery

Median left atrial ejection fraction (LAEF) significantly increased from BL to FU in the entire study population (49.5% (31.1; 59.0) vs. 59.5% (40.8; 66.5); P = 0.0019) as well as in the subgroups (Tables 1 and 2). Median LA volume index (LAVi), showed a significant decrease from BL to FU in both subgroups with the most profound reduction in CLM patients from 62.5 (51.4; 77.4) at BL to 20.1 ml/m2 (18.3; 27.2) at FU (Table 2b). Median LA-GLS increased from 33.2 (14.5; 39.2) at BL to 37.0% (25.2; 44.1, P = 0.0018) at FU in the entire study population. The increase in median LA-GLS remained significant after stratification into both subgroups: Median LA-GLS increased from 36.7 (26.5; 42.3) at BL to 41.3% (34.5; 44.8, P = 0.0262) at FU in ILM patients and from 11.3 (6.4; 21.1) at BL to 21.4% (14.2; 30.7, P = 0.0186) at FU in CLM patients.

Table 1 CMR parameters at Baseline (BL) and at 3 months Follow-up (FU)
Table 2 CMR parameters at BL and at 3 months Follow-up (FU) in patients with infarct-like (a) and cardiomyopathy-like myocarditis (b)

There was no significant change in median right atrial ejection fraction (RAEF) and median right atrial volume index (RAVi) from BL to FU, neither in the entire study population, nor in the subgroup analyses (Tables 1 and 2). However, median RA-GLS significantly increased from BL with 30.8 (22.5; 37.0) to FU with 33.7% (26.8; 45.4, P = 0.0027) in the entire study population. This finding remained consistent in both subgroups: RA-GLS increased from 32.7 (25.8; 41.0) to 35.8% (27.7; 48.0, P = 0.0495) in ILM patients and from 22.8 (13.1; 33.9) to 31.0% (26.0; 40.8) in CLM patients (P = 0.0266), (Table 2, Fig. 2).

Fig. 2
figure 2

Atrial strain at baseline (BL) and follow-up (FU). Left atrial (LA) and right atrial (RA) strain at baseline (BL) and at follow up (FU). Box-Whisker plots of median left atrial global longitudinal strain (LA-GLS) and right atrial global longitudinal strain (RA-GLS) at BL and FU for infarct-like myocarditis (ILM) and cardiomyopathy-like myocarditis (CLM)

There was a significant, inverse correlation of changes in heart rate between BL and FU with changes in LA-GLS (− 0.515 (− 0.0803; 0.630) P = 0.0022) and RA-GLS (− 0.351 (− 0.613; − 0.0201) P = 0.0387).

Ventricular recovery

Median LVEF significantly improved in the entire study population (Table 1). In CLM patients, median LVEF significantly improved from 33 (16.8; 43.3) at BL to 49% (34.8; 65.5, P = 0.0007) at FU. In ILM patients, median LVEF increased from 61 (55; 66) at BL to 64% (62; 72) at FU (P = 0.0015). Median left ventricular end-diastolic volume index (LVEDVi) significantly decreased from BL to FU in CLM and ILM (Table 2). LV-strain showed a significant increase from BL to FU in the entire study population for all parameters (Table 1): Median LVLAXGLS increased from BL to FU in CLM with − 10.4% (− 13.3; − 5.8) vs. − 13.0% (− 19.8; − 9.8, P = 0.0093), but also in ILM with − 17.6% (− 19.3; − 16.4) vs. − 20.8% (− 21.17; − 18.9, P = 0.0017) (Table 2). Median LVSAXGCS also significantly increased from BL to FU in both subgroups (Table 2). Sufficient data to classify LV diastolic dysfunction (DD) by echocardiography were available in 25 patients. All 25 patients had DD at BL. Eighteen (72%) patients were classified as DD I° and 7 (28%) as DD ≥ II°. At baseline, severe DD ≥ II° was more frequent in CLM patients (n = 6, 60%) compared to ILM patients (n = 1, 7%, P = 0.0068). There was a significant decrease in E/e’ from 8.3 (5.4; 12.1) at BL to 5.9 (5.1; 8.1, P = 0.0046) at FU.

Median right ventricular ejection fraction (RVEF) significantly improved from 48.0 (33.8; 55.0) at BL to 57.0% (47.8; 59.5, P = 0.0398) at FU in CML, but not in ILM (Table 2). Median right ventricular end-diastolic volume index (RVEDVi) significantly decreased from BL to FU in both subgroups. There were no significant differences in median RV-global longitudinal strain (RV-GLS) between BL with − 21.5% (− 24.2; − 17.8) and FU with − 22.7% (− 24.6; − 19.9, P = 0.1792) in the entire study population, as well as in the subgroup analyses (Tables 1 and 2).

Discussion

This study evaluated atrial and ventricular function after acute myocarditis by CMR-FT. To the best of our knowledge, development atrial function by FT after acute myocarditis has not been investigated so far. Because of its superior reproducibility, we report total atrial strain (εS) as parameter for atrial reservoir function [22]. Our major findings were: First, there was a significant increase in median LA und RA strain from BL to FU. Second, we found a significant increase in LV, but not in RV strain from BL to FU. Third, atrial and LV functional improvement was found in patients with “cardiomyopathy-like” as well as in patients with “infarct-like” clinical presentation.

Left atrial recovery

In our analyses, median LA-GLS significantly increased from BL to FU (Table 1, Fig. 2). This finding could be related to two different aspects of myocardial recovery after acute myocarditis: First, receding inflammation of atrial myocardium resulted in recovery of myocardial function in analogy to LV myocardial recovery. However, only few reports have actually provided histological data on atrial involvement in myocardial inflammation [24]. Moreover, atrial tissue characterization by CMR is difficult and requires high resolution sequences [25]. Interestingly, there was a modest to moderate inverse correlation of decrease in heart rate with the increase in LA-strain from BL to FU, which certainly reflects myocardial healing and recovery of cardiac function, but could also constitute a potential confounder.[26]. Second, receding myocardial inflammation on the ventricular level [11] resulted in improved systolic and diastolic function with subsequent decline in LA load [27]. LA-strain was recently found to be associated with the extent of DD in patients with preserved LVEF and to be a sensitive marker for severe DD irrespective of LV systolic function [28, 29]. Furthermore, LA-strain appears to be superior compared to LAVi and E/e’-ratio in classifying DD [29]. In our analysis, CLM patients had a significantly lower median BL LA-GLS with 11.3% (6.4; 21.1) compared to ILM patients with 36.7% (26.5; 42.3). Correspondingly, LV diastolic dysfunction at BL was more pronounced in CLM patients compared to ILM patients. At FU, LAVi and E/E’ showed a significant reduction in the overall study population, indicating an improved diastolic function with cessation of myocardial inflammation. Accordingly, we found a significant increase in overall LA-strain at FU (Table 1). In summary, CMR derived LA-strain seems to reflect improvement of systolic and diastolic LV dysfunction after acute myocarditis [30].

However, recent data indicate an independent, incremental value of LA-strain beyond LV function, e.g. as an early indicator of risk for the development of HF, occurrence of atrial fibrillation, but also stroke risk [31,32,33]. Of note, LA-strain offers additional information in HFpEF patients [28, 29], but also HFrEF patients, indicating a role for LA-strain measurements in HF independent from LVEF [13, 14]. In particular, Deferm et al. recently demonstrated a prognostic role of LA reservoir strain in patients with acute HFrEF, independent from changes in LA volumes and LV-function [34]. Accordingly, there could be a role for LA-strain measurements to guide therapy in acute HF beyond conventional parameters such as LVEF [35]. In analogy to these findings in other settings, monitoring LA-strain could improve the individual assessment of healing in patients after acute myocarditis, which is one of the major challenges in clinical routine [2, 16]. Nevertheless, larger studies with long-term follow-up are required before clinical implementation of LA-strain measurements.

Right atrial recovery

We were able to demonstrate a small but significant increase in RA-peak-systolic-GLS from BL to FU in patients with ILM and CLM (Table 2). Analogous to LA-strain, this finding could be related to recovery of atrial myocardial function, but also to improved RV diastolic and systolic function. In patients with acute myocarditis, Dick et al. were able to show a trend towards reduced RA reservoir strain values before [8]. Recently, right heart dysfunction has been found to be of particular relevance in patients with HFpEF, contributing to poor outcome [36]. Jain et al. described significantly reduced right atrial conduit and reservoir strain in patients with HFpEF as well as HFrEF. RA conduit and reservoir strain were found to be independent predictors of mortality irrespective of LVEF or HF status [15]. Tough these studies did not investigate patients with myocarditis in particular, patients with ILM and CLM experience acute HF and can be attributed to HFpEF or HFrEF categories. These findings underline the potential value of monitoring right atrial function in patients after acute myocarditis. Whether RA function in myocarditis provides incremental prognostic information and whether RA strain is able to assess RV diastolic function, needs to be investigated in future studies.

Left ventricular recovery

Recovery of RV- and LV- strain in patients with acute myocarditis has been described before [10]. Luetkens et al. suggested LV-GLS as a predictor of functional recovery in myocarditis [10]. We were able to confirm an increase of global LV-strain parameters from BL to FU in our data (Table 1). Furthermore, LVLAXGLS and LVSAXGCS improved from BL to FU independent of clinical presentation (Table 2). The absolute increase of LVSAXGCS however was smaller compared with LVLAXGLS in the overall study population and just barely significant in ILM patients (Table 2). In CLM, absolute LVSAXGCS at baseline was much lower than in ILM patients, mirroring reduced LVEF and increased LV dilatation (Table 2). The results agree with previous findings that showed a more pronounced deterioration of LVLAXGLS per reduction of LVEF in percent compared to LVSAXGCS [37]. Especially in patients with preserved LVEF, LVLAXGLS therefore probably offers a more sensitive description of myocardial contractility compared to LVSAXGCS and LVEF, making LVLAXGLS the optimal parameter to describe subclinical changes in systolic function [37]. This explains the discrepancy between absolute change of LVLAXGLS and LVSAXGCS in ILM, where LVEF at baseline was preserved (Table 2a). Furthermore, CMR feature-tracking derived myocardial LV-GLS has been independently associated with mortality in patients with dilated cardiomyopathy [5]. A recent meta-analysis indicated a superior prognostic value of LV-GLS over LVEF in predicting major adverse cardiac events [38]. In particular in patients with acute myocarditis, LV-GLS was superior compared to clinical features, LVEF and LGE in this context [6]. It should therefore be considered to adopt LV strain analyses into routine follow-up of patients after acute myocarditis.

Right ventricular recovery

Reduced RV function has been described as independent predictor of adverse outcomes [39]. In our cohort, global RV-GLS did not change significantly between BL and FU. Up to now, conflicting data on RV strain have been published: In contrast to our findings, Luetkens et al. reported an increase in RV-GLS from BL to FU in a similar study population [10]. Baessler et al. reported a paradoxically increased, “supernormal” basal RV-strain-rate in patients with acute myocarditis, whereas no change was found in global RV-strain-rate [7]. These divergent findings could be related to different study populations with different degrees of RV involvement, but also to different extent of LV impairment with subsequently increased “secondary” RV load. Further studies are necessary to better understand the role of RV-GLS changes after acute myocarditis.

Clinical presentation: ILM vs. CLM

Depending on clinical presentation, patients with acute myocarditis can be divided into two major groups: Briefly, patients with chest pain, elevated Troponin levels and ST-alterations on ECG can be classified as having “infarct-like myocarditis” (ILM), whereas patients presenting with symptoms of new-onset heart-failure can be categorized as having “cardiomyopathy-like myocarditis” (CLM) [1, 40, 41]. Despite preserved LVEF, major adverse cardiac events such as recurrent myocarditis, sustained ventricular tachycardia and sudden cardiac death have been attributed to an infarct-like pattern in patients with acute myocarditis [42]. Assessment of cardiac function at baseline and during follow up is routinely performed by echocardiography which offers LVEF as established parameter in characterizing systolic function [1]. Since LVEF can be near normal or preserved in the acute stage of ILM (Table 2a, [41]) monitoring of ejection fraction alone often does not allow for sufficient estimation of cardiac function in all patients. CMR is able to provide complementary tools in this context: While T1 and T2 mapping CMR are able to monitor myocardial inflammation, strain analysis allows the monitoring of myocardial function beyond LVEF [37]. In our cohort, the change in T1 and T2 mapping parameters did not correlate significantly with the observed changes in strain or LVEF. We could observe a significant increase in median LV-GLS from BL to FU in both subgroups, indicating functional improvement independent from LVEF. As mentioned above, LVLAXGLS allows for a better description of subclinical changes in contractile dysfunction than LVEF and LVSAXGCS, especially in patients with ILM [37]. LVLAXGLS therefore constitutes an important parameter depicting LV-systolic dysfunction independent of clinical phenotype in myocarditis. Moreover, the increase in LA-GLS from BL to FU in both subgroups suggests improvement of diastolic function with subsequent relief of LA load beyond recovery of atrial myocardium itself [28, 29]. Upon this premise, CMR-FT seems to offer a phenotype-independent option to assess myocardial function in healing myocarditis [2]. Our group recently published normal values for atrial strain values using the same method and software (Medis Medical Imaging, Leiden, The Netherlands) [19]. In this publication, median two/four chamber LA peak-systolic global longitudinal strain (GLS) were 38.2 (33.0; 43.7)/33.4% (28.4–37.3) and 29.8% (24.1–35.1) for RA-GLS in the control group. Briefly, our findings indicate normalization of median LA- (37.0% (25.2; 44.1)) and RA-GLS (33.7% (26.8; 45.4)) at FU for the overall study population. (Table 1). However, this observation was primarily driven by ILM patients: In ILM patients, we observed an increase of median LA-GLS from (36.7% (26.5; 42.3)) at BL to (41.3% (34.5; 44.8)) at FU (Table 2a). In contrast, CLM patients had lower median LA- (11.3% (6.4; 21.1)) and RA-GLS (22.8% (13.1; 33.9)) values at BL, which improved, but did not fully normalize at FU (LA-GLS 21.4% (14.2; 30.7); RA-GLS 31.0% (26.0; 40.8)) (Table 2b). In line with these findings, strain values of all four cardiac chambers and LVEF remained reduced in CLM patients at FU. Though healing improves cardiac strain also in patients with CML compared to the acute phase, patients need close clinical follow up in order to offer early medical therapy according to current guidelines. [2]

Limitations

First, our study population was relatively small. Furthermore, the study population is affected by an inherent selection bias, since patients with clinical worsening and subsequent ICD-implantation or heart transplantation were not able to participate at FU. However, our study population is representative for the majority of patients after acute myocarditis, who undergo regular follow-up visits in outpatient settings. Nevertheless, larger studies with long-term follow-up are warranted to better understand the relevance of functional recovery by atrial CMR-FT compared to myocardial tissue characterization, but also conventional biomarkers after acute myocarditis. We observed a significant correlation between changes in LA- and RA-strain with changes in heart rate from BL to FU, which could be explained by myocardial healing. However, changes in heart rate could constitute a potential technical confounder for the observed changes in LA-strain [26]. In addition, though we could demonstrate a very good inter-observer-agreement for the measured strain values, test–retest (inter-scan) variability was not addressed in this study [43].

Conclusion

Our findings demonstrate recovery of LA and RA function by CMR-FT strain analyses in patients after acute myocarditis independent from clinical presentation. Monitoring atrial strain could be an important tool for an individual assessment of healing after acute myocarditis.