Evaluation of deep learning-based reconstruction late gadolinium enhancement images for identifying patients with clinically unrecognized myocardial infarction

Background The presence of infarction in patients with unrecognized myocardial infarction (UMI) is a critical feature in predicting adverse cardiac events. This study aimed to compare the detection rate of UMI using conventional and deep learning reconstruction (DLR)-based late gadolinium enhancement (LGEO and LGEDL, respectively) and evaluate optimal quantification parameters to enhance diagnosis and management of suspected patients with UMI. Methods This prospective study included 98 patients (68 men; mean age: 55.8 ± 8.1 years) with suspected UMI treated at our hospital from April 2022 to August 2023. LGEO and LGEDL images were obtained using conventional and commercially available inline DLR algorithms. The myocardial signal-to-noise ratio (SNR), contrast-to-noise ratio (CNR), and percentage of enhanced area (Parea) employing the signal threshold versus reference mean (STRM) approach, which correlates the signal intensity (SI) within areas of interest with the average SI of normal regions, were analyzed. Analysis was performed using the standard deviation (SD) threshold approach (2SD–5SD) and full width at half maximum (FWHM) method. The diagnostic efficacies based on LGEDL and LGEO images were calculated. Results The SNRDL and CNRDL were two times better than the SNRO and CNRO, respectively (P < 0.05). Parea−DL was elevated compared to Parea−O using the threshold methods (P < 0.05); however, no intergroup difference was found based on the FWHM method (P > 0.05). The Parea−DL and Parea−O also differed except between the 2SD and 3SD and the 4SD/5SD and FWHM methods (P < 0.05). The receiver operating characteristic curve analysis revealed that each SD method exhibited good diagnostic efficacy for detecting UMI, with the Parea−DL having the best diagnostic efficacy based on the 5SD method (P < 0.05). Overall, the LGEDL images had better image quality. Strong diagnostic efficacy for UMI identification was achieved when the STRM was ≥ 4SD and ≥ 3SD for the LGEDL and LGEO, respectively. Conclusions STRM selection for LGEDL magnetic resonance images helps improve clinical decision-making in patients with UMI. This study underscored the importance of STRM selection for analyzing LGEDL images to enhance diagnostic accuracy and clinical decision-making for patients with UMI, further providing better cardiovascular care. Supplementary Information The online version contains supplementary material available at 10.1186/s12880-024-01308-2.


Background
Myocardial infarction (MI) is diagnosed based on the detection of acute myocardial injury according to cardiac biomarker abnormalities in the context of acute myocardial ischemia [1].Unrecognized MI (UMI) is a type of MI that has yet to be clinically diagnosed, with the prevalence increasing by 10.0% every decade [2].Delayed detection due to atypical symptoms can delay treatment, leading to poor prognosis [3].Failure to achieve reperfusion within a few hours after blood flow cessation may cause myocardial apoptosis in vessel-supplied regions.Therefore, determining the presence or absence of MI and quantifying related variables are crucial in improving the diagnosis, treatment, and prognosis [4,5].
Cardiac magnetic resonance (CMR) imaging is a promising tool for MI detection because of good tissue contrast and spatial resolution.However, patient compliance is challenging for several reasons, such as the requirement to acquire each high-resolution slice and the need for stable respiration; furthermore, certain conditions, including unstable heartbeat and arrhythmia, can cause motion artifacts on free-breathing scans.As relatively shorter breath-holds are required to acquire more slices, higher-spatial resolution late gadolinium enhancement (LGE) is most frequently utilized in magnetic resonance imaging (MRI) to observe and quantify the degree of myocardial necrosis and microvascular occlusion.Although the enhancement is achieved semi-automatically using post-processing software, the initial sketch of the endocardium, epicardium, enhanced myocardium, and remote normal myocardium relies on the reader's experience to some extent [6].Additionally, a previous study reported that LGE could identify only 23 of the 872 participants (2.6%) with UMI [7].The clinical significance of UMI has been reported using different imaging techniques in diagnosing, refining risk stratification, and guiding clinical decisions for treatments.All underscored the role of CMR in improving the detection accuracy of UMIs, which may affect adverse cardiac outcomes and optimize cardiovascular disease management [8][9][10].Therefore, timely and accurate UMI identification and assessment are fundamental for patient stratification and therapeutic planning [4,5,11].In practice, despite many applications of standard deviation (SD) and full width at half maximum (FWHM) techniques, no consensus exists for quantifying scars on LGE images; this challenge persists across different cardiac diseases [12][13][14].Obviously, a gap exists in current diagnostic frameworks for analyzing myocardium delayed enhancement.
Deep learning (DL) methods can improve image quality and eliminate intra-and inter-observer variability, enabling more accurate diagnosis and treatment strategies [15,16] and segmentation for precisely sketched lesions [17][18][19][20][21], among others.However, no DL reconstruction (DLR)-based magnetic resonance (MR) studies have evaluated patients with suspected UMI.Therefore, this study aimed to explore the feasibility and diagnostic performance of DLR-based LGE imaging (LGE DL ) for patients with UMI compared with that of conventional imaging (LGE O ) and propose an appropriate signal threshold versus reference mean (STRM) for analyzing LGE DL .

Study population
This study prospectively recruited 98 patients (68 men and 30 women, mean age: 55.8 ± 8.1 years) who presented at our hospital between April 2022 and August 2023 without typical MI symptoms, such as angina pectoris of cardiogenic origin but with suspected UMI after a physical examination.Based on the guidelines of European and American associations and previous reports [1,7,22], the inclusion criteria were as follows: (i) the absence of typical angina symptoms; (ii) the presence of elevated or decreased serum cardiac troponin (cTn) levels, with at least one instance of elevation above the upper limit of the normal value (the 99th percentile of the reference value's upper limit); (iii) prior evidence of MI on electrocardiography in the absence of left ventricular hypertrophy and left bundle branch block; and (iv) no prior history of oncological disease or surgery for cardiovascular diseases.The exclusion criteria were as follows: (i) clinically unstable condition, decompensated heart failure, contraindication to CMR, an estimated glomerular filtration rate ≤ 30 mL/min, and contraindication to the use of gadolinium contrast; and (ii) LGE images that could not be used for clinical diagnosis and objective assessments (Fig. 1).

CMR examination and image construction
All patients underwent a routine cardiac MRI examination, including a short-axis LGE imaging sequence, on a 3.0-T MRI scanner (Signa Architect, GE Healthcare, Waukesha, WI, USA) at our hospital.A new commercial inline deep-learning-based reconstruction (DLR, brand name: AIR™ Recon DL, DV29.1_R04,GE Healthcare, USA) employs no bias terms and rectified linear unit activations to identify 4.4 million features on directly received image data immediately after scanning on an MR console computer to reduce noise and Gibbs artifacts, and further eliminate intra-and interobserver differences [13,16].The parameters for the LGE sequence were as follows: echo time = 2.7 ms; repetition time = 5.6 ms; flip angle = 25°; field of view = 34 cm × 34 cm; matrix = 260 × 174; slice thickness = 8 mm; slice spacing = 2 mm; receiver bandwidth = 83.33 kHz; views per segment = 24; number of excitations = 1; and theoretical acquisition time = 8 s×nine heart beats.The LGE O and LGE DL were simultaneously generated using conventional inline reconstruction and AIR™ Recon DL algorithms.Fifteen minutes before LGE sequence scanning, a single bolus of 0.1 mmol/kg (0.2 ml/kg) Gadobenate Dimeglumine (Bracco Imaging S.P.A., Milano, Italy) was administered, followed by 20-mL saline flush at a flow rate of 2 ml/s [23].This dosage was selected based on its efficacy of myocardial enhancement for visualization under the condition of patient safety.

Assessment of myocardial enhancement area and diagnostic efficacy
Ultimately, data from 61 patients with myocardial enhancement were included in the analysis (43 men [70.5%] and 18 women [29.5%]), with a mean age of 55.9 ± 8.7 years (Fig. 1).The percentage of whole-heart myocardial enhancement area (P area ) in segments S1-S16 was assessed semi-quantitatively to diagnose cardiovascular disease using Circle Cardiovascular Imaging Inc. (cvi 42 , Circle Cardiovascular Imaging Inc., Calgary, AB, Canada).The delayed enhancement area (i.e., scar size) was subsequently quantified based on threshold methods, which involve adding 2-5 times SD to the mean signal intensity (SI) of the reference myocardium, and the FWHM method, which identifies the half maximum SI at the full width of SI distribution within one region of interest (ROI) in the myocardial tissue.The P area was calculated as the scar size divided by the myocardial volume.Furthermore, the diagnostic efficacy of the P area of LGE DL and LGE O images (P area−DL and P area−O , respectively) in differentiating patients with UMI was assessed, with the clinical diagnosis of UMI as the gold standard.

Theory/calculation CMR image assessment
Qualitative and quantitative imaging evaluations were performed double-blindedly by two radiologists with > 5 years of experience in CMR diagnosis.Moreover, one of the radiologists repeated the assessment 1 month later.

Image quality
For the objective evaluation of image quality, ROIs were on LGE O and LGE DL images to determine the SI of the normal myocardium (SI Myo−O and SI Myo−DL , respectively) and myocardial delayed enhancement area (SI MDEA−O and SI MDEA−DL , respectively), as well as the SD of the background noise at the corner of the images (SD BG−O and SD BG−DL , respectively) and the myocardial delayed enhancement area (SD MDEA−O and SD MDEA−DL , respectively) (Fig. 2).Additionally, for LGE O and LGE DL images, the myocardial signal-to-noise ratios (SNRs) (SNR O and SNR DL , respectively) and contrast-to-noise ratios (CNRs) (CNR O and CNR DL , respectively) were calculated [9,10,24,25] using the following formulae: The short-axis LGE O and LGE DL images were divided into 16 segments based on the American Heart Association criteria, and the SNR and CNR of each segment were calculated.

Statistical analysis
All data were statistically analyzed using R-project software (version 4.0.4,http://www.r-project.org).Quantitative data are expressed as either the‾x ± SD or median (interquartile range).All quantitative data were analyzed using either a paired t-test or a Wilcoxon signed-rank test depending on the results of the Shapiro-Wilk and Levene's tests, which were used to assess variance homogeneity and data normality, respectively.To control the false discovery rate, we applied the Benjamini--Hochberg method for multiple comparison corrections.The intraclass correlation coefficients (ICCs) of the objective quantitative indicators, including the SNR, CNR, SD, and P area for LGE O and LGE DL images (SNR O , SNR DL , CNR O , CNR DL , SD O , SD DL , P area−O and P area−DL , respectively) were quantified to assess the degree of intra-and interobserver agreement.Receiver operating characteristic (ROC) curves for P area−DL and P area−O were constructed using the different threshold methods to determine and compare their diagnostic efficacies for the UMI or non-UMI groups based on the area under the curve (AUC).All statistical significance was set at P < 0.05.

Objective evaluation of image quality
The SDs of the normal myocardium, delayed myocardial enhancement areas, and background of the images are presented in Table 1 2a-b-.I, Supplementary Fig. 1a-d.I).The SI DL values were higher than the respective SI O values for all segments (P < 0.05) except for S7-S9 and S11.In comparing the SI Myo−DL and SI Myo−O values, the SI DL values were higher than the corresponding SI O values for S1-S6, S10, and S12-S16 (P < 0.05).The SI DL values were slightly higher than the corresponding SI O values for S7-S9 or S11; however, the difference was not significant (P > 0.05) (Fig. 3a).

P area assessment
The myocardial enhancement area was semi-quantitatively analyzed using various SD thresholds and the FWHM method.For the 2SD (Fig. 3b.I, Supplementary Figs.1), 3SD (Figs. 2 and 3b.II, Supplementary Fig. 1), and 5D methods (Figs. 2 and 3b.IV, Supplementary Fig. 1), the P area−DL values for the overall myocardium were higher than the corresponding P area−O values for all 16 segments.For the 4SD method, the P area−DL values of the overall myocardium were higher than the corresponding P area−O values only in S1-S12 (Figs. 2 and 3b.III, Supplementary Fig. 1).For the FWHM method (Figs. 2  and 3.c, Supplementary Fig. 1), the P area−DL values were slightly higher than the corresponding P area−O values for all segments.
Regarding the DLR-based P area , the overall different threshold and FWHM-based P area−DL values were higher than those based on any other approach (all P < 0.05).Regarding the P area−O , the values for the 2SD threshold were significantly higher than those based on other approaches (all P < 0.05) (Table 2).

Assessment of the consistency of the quantitative measurements
The degree of intra-and inter-observer agreement for the objective measurements (SD Myo , SD MDEA , SD BG , SNR, CNR, and SI Myo ) and P area between LGE DL and LGE O images was good based on the various SD and FWHM methods (for objective measurements: all ICCs > 0.60, all P < 0.05; for P area : all ICCs > 0.70, P < 0.05).These measurements were better for LGE DL images than for LGE O images (Figs. 4 and 5).

Analysis and comparison of diagnostic efficacy
All SD methods exhibited good diagnostic efficacy for UMI, with AUC values of the ROC curves ≥ 0.78.The   P area−DL based on the 5SD threshold method exhibited the optimal diagnostic efficacy of 0.891 (sensitivity = 0.688 and specificity = 1).For the conventional imaging enhancement, the P area−O based on the 3SD method exhibited the optimal diagnostic efficacy of 0.840.The diagnostic efficacy was better for LGE DL images than for LGE O images for UMI detection for every SD threshold method, whereas it was not different between LGE DL and LGE O parameters based on the FWHM method (Table 3; Fig. 6).

Discussion
This study compared LGE O and LGE DL images based on different SD thresholds and the FWHM method.The significant differences in P area values between LGE O and LGE DL images for the SD threshold methods but not for the FWHM method suggested that the STRM should be ≥ 3, regardless of whether conventional or DLRbased LGE images are used, as previously reported.An STRM ≥ 4 and P area−DL values based on the 5SD threshold exhibited the highest diagnostic efficacy for detecting UMI.Additionally, the LGE DL images generated in this study could display the delayed enhancement area in patients with UMI for the first time, with significantly better image quality than was previously achievable with LGE O images, such as artifacts in the myocardium, intensified foci and lower background noise, lower SD, and higher SNR and CNR values in all patients with UMI.Thus, LGE DL imaging can improve diagnostic confidence without impacting diagnostic efficacy.
The presence of an infarction in patients with UMI is a critical feature for predicting adverse cardiac events [26][27][28].The P area on LGE images is the most frequently used direct indicator of irreversible damage at the pathological tissue level and can predict the treatment response to cardioprotective interventions [29,30].However, the clinical approach for quantifying the myocardial enhancement area is not uniform, with SD thresholds used in some instances and the FWHM method employed in others.Additionally, the generation of LGE images using conventional reconstruction and DLR-based methods is inconsistent.Generally, an STRM ≥ 3SD is the optimal reference threshold for clinical use.Quantifying the SD thresholds depends predominantly on the SI and SD of the ROIs drawn in the distal normal myocardium; however, the image quality of the remote normal myocardium may affect the visual sketching of the area to avoid the delayed lesion intensification on LGE DL images [31].For example, using a lower SD threshold of the distal myocardium leads to a significantly lower threshold for encompassing the extent of delayed enhancement, resulting in underestimation [13].The SD values, including SD Myo , SD MDEA , and SD BG of the LGE DL images, showed  [13,16,26,32,33].This was the first study to evaluate and directly compare LGE DL and LGE O images of delayed intensification foci in patients with UMI.The diagnostic performance of the P area−DL was higher than that of the P area−O for the threshold approaches, especially for the P area−DL based on the 5SD threshold, which exhibited the best AUC (0.891).For LGE O images, the P area−O based on the 3SD threshold exhibited the optimal AUC of 0.840, consistent with data from previous studies recommending using an STRM ≥ 3SD for infarct size.This study recruited patients with UMI without clinically significant cardiogenic chest pain and with a relatively small range of reinforcing foci; these results confirm that the 3SD threshold is sufficient for conventional LGE images.In contrast, a threshold ≥ 4SD should be used for DLR LGE images to optimize the intra-and inter-reader agreement and diagnostic efficacy.The diagnosis of the extent of infarction in UMI-related cases using the 4SD threshold was possibly a more reliable parameter for LGE O and LGE DL images despite the better diagnostic efficacy of the 5SD threshold for LGE DL imaging.Furthermore, the detection rate of UMI was 67% (63/91); this rate was similar for LGE O and LGE DL images despite the better image quality and more reliable assessment of pathological features on LGE DL imaging.
This study has some limitations.First, all participants were recruited using a single-center design, and only those who underwent an MR examination were included for analysis, limiting the generalizability of our results.Despite LGE images with high diagnostic accuracy of MI detection, the final diagnosis relies on experienced radiologists due to the lack of pathological validation for delayed enhancement areas on LGE images.Therefore, to enhance the robustness of result generalization, multicenter and large data, including comparison of P area−DL and P area−O using various SD and FWHM methods and validation of the accuracy and reliability for UMI diagnosis should be considered for future LGE O or LGE DL .

Conclusions
The selection of SD thresholds for LGE DL (≥ 4SD) and LGEO (≥ 3SD) images was recommended for future research, as the difference between P area−DL and P area−O affected diagnostic efficacy and clinical decision-making in patients with UMI.Moreover, P area−DL and P area−O were similar when the FWHM method was used, implying LGE DL images retained informational integrity.Despite the same UMI detection rates between LGE O and LGE DL images, the LGE DL images showed superior image quality and reliable features for diagnosis with more confidence.Therefore, STRM selection and diagnostic outcomes should be carefully utilized and interpreted, particularly for DLR-based CMR images.

Fig. 2
Fig. 2 Schematic diagram of P area using accordingly (a.II) and (b.III) 4SD, (b.II) 3SD, (a.III) 5SD, (a.IV) and (b.IV) FWHM methods for (a) LGE DL images, (b) LGE O images, and (a.V) electrocardiogram of a patient with UMI. Figure 2(a) shows clearer, less noisy, more uniform normal myocardial signal and better contrast between the enhancement area and normal myocardium than Fig. 2(b).The patient with UMI underwent stress perfusion myocardium and received an intravenous injection of 20 mCi 99mTc-MIBI.The stress perfusion maps as Fig. 2(b.V) supported our P area maps with clearer myocardium enhancement in the enlarged left ventricle, with the morphological anomaly, relatively light sparsity of 20 mCi 99mTc-MIBI (a radiation tracker, RT) in the middle and basal segments of the anterior wall and the middle segment of the anteroseptal wall, relatively strong sparsity of RTs in the apex, the apical segment of the septal wall, the middle and basal segments of the posteroseptal wall, the apical, middle, and basal segments of the inferior wall, the apical segment of the lateral wall, and the middle and basal segments of posterolateral, and normal perfusion in the remaining myocardium.Note: SD: standard deviation; 2, 3, 4, and 5SD threshold methods: mean P area respectively adding 2, 3, 4, and 5 times of standard deviation of P area as the threshold for myocardial enhancement area; FWHM: full width at half maximum; LGE DL : deep learning-based reconstruction late gadolinium enhancement; LGE O : conventionally constructed late gadolinium enhancement; UMI: unrecognized myocardial infarction Fig. 2 Schematic diagram of P area using accordingly (a.II) and (b.III) 4SD, (b.II) 3SD, (a.III) 5SD, (a.IV) and (b.IV) FWHM methods for (a) LGE DL images, (b) LGE O images, and (a.V) electrocardiogram of a patient with UMI. Figure 2(a) shows clearer, less noisy, more uniform normal myocardial signal and better contrast between the enhancement area and normal myocardium than Fig. 2(b).The patient with UMI underwent stress perfusion myocardium and received an intravenous injection of 20 mCi 99mTc-MIBI.The stress perfusion maps as Fig. 2(b.V) supported our P area maps with clearer myocardium enhancement in the enlarged left ventricle, with the morphological anomaly, relatively light sparsity of 20 mCi 99mTc-MIBI (a radiation tracker, RT) in the middle and basal segments of the anterior wall and the middle segment of the anteroseptal wall, relatively strong sparsity of RTs in the apex, the apical segment of the septal wall, the middle and basal segments of the posteroseptal wall, the apical, middle, and basal segments of the inferior wall, the apical segment of the lateral wall, and the middle and basal segments of posterolateral, and normal perfusion in the remaining myocardium.Note: SD: standard deviation; 2, 3, 4, and 5SD threshold methods: mean P area respectively adding 2, 3, 4, and 5 times of standard deviation of P area as the threshold for myocardial enhancement area; FWHM: full width at half maximum; LGE DL : deep learning-based reconstruction late gadolinium enhancement; LGE O : conventionally constructed late gadolinium enhancement; UMI: unrecognized myocardial infarction

Fig. 3 Fig. 4
Fig. 3 (a) Signal intensity of the left ventricular myocardial on LGE DL and LGE O images.Percentage areas of left ventricular myocardial enhancement in LGE DL and LGE O images using (b.I) 2SD, (b.II) 3SD, (b.III) 4SD, (b.IV) 5SD, and (c) FWHM methods for quantification.Note: SI: signal intensity; WM: whole myocardium; P area : percentage of myocardial enhancement area; LGE DL : deep learning-based reconstruction late gadolinium enhancement; LGE O : conventionally constructed late gadolinium enhancement;; SD: standard deviation; 2, 3, 4, and 5SD threshold methods: mean P area respectively adding 2, 3, 4, and 5 times of standard deviation of P area as the threshold for myocardial enhancement area; FWHM: full width at half maximum; DL, deep learning late gadolinium enhancement; O: original late gadolinium enhancement

Table 1
Objective evaluation of image quality for LGE DL and LGE O Note: SD, standard deviation; SNR, signal-to-noise ratio; WM, whole myocardium; MDEA, myocardium delayed enhancement area; BG, background

Table 2
Differences between different-threshold and FWHM methods

Table 3
Area under the curve (AUC) for differentiation of UMI or non-UMI groups 6ig.6Diagnostic efficacy for UMI.Note: UMI: unrecognised myocardial infarction; SD: standard deviation; P area-DL , percentage of myocardial enhancement area with deep learning late gadolinium enhancement; P area-O , percentage of myocardial enhancement area with original late gadolinium enhancement; 2, 3, 4, and 5SD threshold methods, mean P area respectively adding 2, 3, 4, and 5 times of standard deviation of P area as the threshold for myocardial enhancement area; FWHM: full width at half maximum