Introduction

Patients with heart failure (HF) experience mortality and hospital readmission rates exceeding 15% and 30%, respectively, within 60 to 90 days post-discharge, placing a significant burden on medical systems [1]. HF can be classified into ischemic and nonischemic types according to the 2022 AHA/ACC/HF-SA Guidelines, both of which pose, complex challenges in clinical practice [2, 3]. Nonischemic HF (non-ICHF), characterized by myocardial hypertrophy, valve dysfunction, myocardial inflammation, and various structural and functional abnormalities, lacks effective prognostic biomarkers [4]. While biomarkers such as high-sensitivity C-reactive protein (hs-CRP), B-type natriuretic peptide (BNP), N-terminal pro-brain natriuretic peptide (NT-proBNP), and left ventricular ejection fraction (LVEF) have been broadly applied in assessing the prognosis of non-ICHFHF patients, their limitations in providing accurate prognostic information persist [5, 6].

The application of computed tomography coronary angiography (CCTA) in non-ICHF has emerged as a burgeoning focus within the field of cardiovascular medicine [7]. While initially primarily utilized for the assessment of coronary artery diseases, CCTA has expanded its scope to encompass comprehensive evaluations of cardiac structure and function, offering novel diagnostic and therapeutic opportunities for non-ICHF patients [8, 9]. First, CCTA employs high-resolution, three-dimensional imaging, which can be used to depict the anatomical structure of coronary arteries, facilitating the exclusion of coronary artery stenosis or obstruction [10]. Second, CCTA is employed for the assessment of cardiac structure, including myocardial thickness, ventricular wall motion abnormalities, valvular function, and the calculation of the cardiac lumen volume to muscle mass (V/M) ratio [6, 11], all of which are pivotal factors in the development and treatment of non-ICHF. Recent research has also indicated the potential of CCTA to measure epicardial adipose tissue (EAT) and pericoronary adipose tissue (PCTA), factors of emerging significance in non-ICHF investigations [12,13,14]. Therefore, CCTA provides non-ICHF patients with a multifaceted, noninvasive assessment option, integrating information from coronary arteries, cardiac structure and function, and cardiac adipose tissue. It holds potential as a potent tool for prognosis assessment and treatment decision-making, offering new hope for the management and prognosis of non-ICHF patients.

In this research, we utilize EAT and the V/M ratio to assess cardiac metabolic capacity, PCAT to assess cardiac inflammatory responses, and the CT-derived fractional flow reserve (CTFFR) to evaluate coronary hemodynamics. Our objective was to construct an integrated model for prognosticating the clinical outcomes of these patients.

Methods

Study Population

Between January 2019 and June 2022, a retrospective cohort study including 107 consecutive patients who underwent CCTA followed by invasive coronary angiography (ICA) within 3 months at The Second Affiliated Hospital, Zhejiang University School of Medicine, Hangzhou, Zhejiang, China, was performed. A total of 129 healthy controls were also included in the study to compare the CCTA-derived parameters of the patients with those in the normal population. A flowchart of this study is illustrated in Fig. 1.

Fig. 1
figure 1

Study flow chart

All participants were of Asian ethnicity. The inclusion criteria were as follows: (1) heart failure defined according to the European Society of Cardiology (ESC) guidelines [2]; (2) New York Heart Association [NYHA] class II-IV or signs consistent with the Framingham criteria; (3) no coronary artery stenosis on coronary angiography; and. (4) CT-FFR > 0.8. The exclusion criteria included (1) missing images or poor image quality; (2) incomplete laboratory tests; and (3)) incomplete ultrasound data and loss to follow-up. All patient information, including demographic characteristics, medical history, laboratory tests, and echocardiography results, was collected after admission.

All patients provided informed consent. This study was approved by the Human Research Ethics Committee of the Second Affiliated Hospital, Zhejiang University School of Medicine (I20221200).

Data collection

Medical records including medical history, laboratory data, medication, and clinical course were reviewed. Laboratory blood analysis was performed within 48 h before CCTA. Based on previous literature, metabolic syndrome was diagnosed in patients with a body-mass index (BMI) greater than 25 kg/m2 and two or three of the following conditions [15]: (1) fasting plasma glucose ≥ 110 mg/dl or HbA1c (NSGP) ≥ 5.5%, (2) systolic blood pressure ≥ 130 mmHg or diastolic blood pressure ≥ 85 mmHg, and (3) triglycerides ≥ 150 mg/dl or high-density lipoprotein (HDL) cholesterol < 40 mg/dl. Throughout the entire echocardiographic assessment, a single-lead electrocardiogram was recorded. Following the guidelines of the American Society of Echocardiography, a phased-array transducer with a fusion frequency of 2–4 MHz was employed [16]. Imaging was conducted in M-mode and 2D modes, capturing images in parasternal long-axis, short-axis, apical four-chamber, and two-chamber views and at the aortic root for evaluation. Measurements included the thicknesses of the interventricular septum (IVS), the anterior wall of the left ventricle (LVID), and the posterior wall of the left ventricle (LVPW) during both systole and diastole of the heart. The diameter of the connection point between the aorta and the sinus of Valsalva (AO-stj) and anteroposterior diameter of the left atrium (LA-ap) were also measured. LVEF was calculated using the biplane Simpson's method, averaged over three consecutive heartbeats. BMI was calculated as weight (kilograms) divided by height (meters) squared.

Cardiac Computed Tomography Protocol and Image Acquisition

All participants underwent dual-source CT scans (Somatom Definition Flash or SOMATOM Force, Siemens Healthcare, Germany), which included coronary artery calcium scoring (CACS) and CCTA. The following parameters were used for the CACS protocol: tube voltage, 120 kV; tube current, 80 mAs with automated tube current modulation (CARE Dose 4D, Siemens Healthineers); and section thickness, 2 mm with a 1.5 mm increment. CCTA was conducted with a prospective ECG-gated sequence: CARE kV (reference tube voltage, 100 kV);); CARE Dose 4D (reference tube current, 288 mAs); acquisition phase, 30–75%; 65% R-R interval; and reconstructed slice thickness, 0.75 mm with a 0.5 mm increment. The images were reconstructed by advanced model iterative reconstruction (ADMIRE) level 3; kernel Bv40. Contrast medium was administered at a dose ranging from 30 to 60 mL, and the flow rate ranged from 4.0 to 6.0 mL/s according to the patients’ BMI, heart rate, and tube voltage. CCTA was performed after the injection of contrast medium (iodine 370 mg/mL [Ultravist, Bayer Schering Pharma, Berlin, Germany]) via an 18–20 G intravenous catheter, followed by a 30 mL saline flush using a dual-head power injector.

CT-derived and structural function index measurement

The CT-derived EAT (CTEAT), mass (CTMASS), V/M (CTV/M), FFR (CTFFR) and PCAT (CTPCAT) were assessed on a workstation (syngo.via VB40, Siemens Healthineers, Germany). Two radiologists who were blinded to the patient data (SX Lin [observer 1] and CJ Liu [observer 2], with over 5 years of experience in CT diagnostic imaging). CTFFR was calculated and analysed based on a software prototype (cFFR, syngo.via Frontier, version 3.0.1, Siemens Healthineers, Germany).

CTEAT, defined as the total amount of adipose tissue between the surface of the heart and the visceral layer of the pericardium, was measured by volumetry on short-axis slices with a thickness of 0.75 mm, ranging from the level of the pulmonary bifurcation to the apex and within a threshold range of -190 to -30 Hounsfield units (HU) to determine the total volume of tissue [17]. The adipose tissue around the coronary artery is part of the EAT and was defined as the attenuation coefficient of the fat tissue voxels within a distance from the coronary arterial wall equal to the corresponding vessel diameter. The attenuation coefficient was calculated separately for the right coronary artery (RCA), left anterior descending artery (LAD), and left circumflex artery (LCX) using the coronary artery analysis module [18].

The cardiac mass was determined using the syngo.via cardiac analysis module, which automatically traces the endocardial and epicardial borders to generate the volume of the heart. The myocardial volumes were converted to left ventricle mass (M) by assuming a constant tissue density of myocardium (1.05 g/cc). CTV/M was calculated by dividing the luminal volume of the coronary artery by the myocardial mass of the heart [9, 19]. The luminal volume of the coronary artery was extracted by a deep learning calculation method on the workstation. In this study, the right coronary artery was selected for coronary lumen volume analysis given its simplicity and thick diameter in the Chinese population.

Follow-up

Follow-up included follow-up data and results obtained from January 2019 and June 2022 in the HIS system as a follow-up of disease and care, clinic visits, and telephone interviews. The endpoints were patient readmission (defined as a hospital admission for which HF was the primary reason and requiring either diuretic, inotropic, or intravenous nitrate therapy). All statuses were reviewed by 2 independent investigators who used previously described criteria.

Statistical Analysis

Data were analysed using R software (Version 4.0). Continuous variables are reported as the mean ± SD or median (interquartile range, IQR) and were compared using independent samples t tests and one-way analysis of variance or Wilcoxon signed-rank and Kruskal‒Wallis tests. Categorical variables are presented as absolute values and percentages. The correlation between healthy controls and HF patients was analysed using the Pearson correlation. Inter- and intrao bserver agreement for the CTEAT, CTMASS, CTV/M, CTFFR and CTPCAT were performed using Bland–Altman analysis.

For survival analysis, survival and proportional hazards assumptions were estimated by the Kaplan‒Meier method, and any differences were evaluated with the stratified log-rank test. The optimal cut-off point was identified using the maximally selected rank statistic (maxstat) [20], which is a ranking statistic for maximum selection. Odds ratios (ORs) and 95% confidence intervals (CIs) were determined from univariate and multivariate logistic regression analyses to determine any factors associated with the endpoint. A favourable patient prognosis was taken as the dependent variable, and the statistically significant factors were included in the logistic regression model as independent variables for regression analysis. Three risk models were created based on the best traditional parameter variables (model 1), the CCTA-driven parameters (model 2), and the combination of the two sets of parameters (model 3). Receiver operating characteristic (ROC) curve analysis was performed to evaluate the prognostic accuracy of CTEAT, CTMASS, CTV/M, CTFFR and CTPCT for the endpoint.

Results

Baseline Patient Characteristics

The clinicopathological characteristics of the patients are summarized in Table 1. A total of 236 patients were enrolled in this clinical study. Compared to the healthy control group, the non-ICHF group was significantly different in age, height, diastolic blood pressure and metabolic syndrome. Among the laboratory data, CHO tended to be lower in the non-ICHF than in the control group (3.9 mmol/L vs. 3.6 mmol/L, P = 0.003), while TG (1.3 mmol/L vs. 1.1 mmol/L, P = 0.009), HDL-C (1.2 mmol/L vs. 1.1 mmol/L, P = 0.004), and LDL_C (1.9 mmol/L vs. 1.7 mmol/L, P = 0.0035) were significantly different (Table 1).

Table 1 Demographic characteristics of the healthy control and non-ischemic HF

CCTA-Derived Parameters

Except for CTV/M (26.2% [IQR, 24.4; 27.8] vs. 22.0% [IQR, 21.0; 24.0], P < 0.01), the CCTA-derived parameters in the healthy groups were significantly lower (P < 0.05) than those in the non-ICHF group (Table 1). In non-ICHF, the volume of adipose tissue around the proximal LAD (-79.8 HU [IQR, -78.2; -81.3]) was higher than that around the RCA (-68.9 HU [IQR, -65.3; -73.1]) and LCX (-73.3 HU [IQR, -67.1; -75.4]) (Table 1, Fig. 2).

Fig. 2 
figure 2

The Boxplot for CCTA-derived parameters in healthy control and non-ischemic Heart failure; (A) CT derived epicardial adipose tissue; (B) CT derived cardiac muscle mass; (C) CT derived cardiac muscle mass to lumen volume; (D) peri-coronary adipose in the right coronary artery; (E) peri-coronary adipose in the left anterior descending artery; (F) peri-coronary adipose in the circumflex branch of the left coronary artery

There was a moderate to strong correlation between healthy patients and non-ICHFHF patients in the attenuation of the PCAT of the LAD (r = 0.55; P < 0.0001), the right circumflex artery and LCX (r = 0.422; P < 0.0001), the coronary lumen volume (r = -0.57; P < 0.0001) and CT-derived cardiac V/M ratio (r = -0.72; P < 0.0001) (Supplementary Fig. 1).

Clinical Endpoints

The last follow-up was in October 2022. During a mean follow-up of 6 ± 0.9 months (range 2–28 months), 40 patients were rehospitalized for heart failure. These patients demonstrated a significant difference in terms of HDL-C (1.1 mmol/L [0.9;1.2] vs. 1.2 mmol/L [1.1;1.4], P = 0.005) and CRP (4.8 mg/L [2.5;14.1] vs. 3.5 mg/L [1.5;7.0], P = 0.033) (Table 2). Univariate logistic regression analysis showed that HDL-C (OR: 7.68, 95% CI: 1.57–3.76, P = 0.08) and CRP (OR: 0.92, 95% CI: 0.85–0.99, P = 0.08) were associated with readmission. The multivariable regression model showed that HDL-C (OR: 2.35, 95% CI: 1.31–4.29, P = 0.08) and TG (OR: 2.16,: 1.14–3.12, P < 0.01) were significantly associated with readmission in the study population (Table 3).

Table 2 Demographic characteristics of the readmission non-ischemic HF
Table 3 Univariable and Multivariable logistic Analysis for the Prediction of non-ischemic HF readmission

Non-ICHF Patients sustaining readmission also had higher EAT (63.9 cm3 [-55.5; -69.5] vs. -51.7 cm3 [-47.1; -57.2], P < 0.01), PCATRCA (-66.5 HU [-64.1; -71.1] vs. -71.2 HU [-67.8; -77.6], P < 0.001), PCATLAD (-79.0 HU ± 2.0 vs -81.1 ± 2.0, P < 0.001), PCATLCX and CTM (-79.4 HU ± 6.4 vs. -88.6 HU ± 7.5, P < 0.001). However, there was no difference between the groups in terms of the CTFFR (0.9 [0.8;0.9] vs. 0.9 [0.8;0.9], P = 0.987) (Table 2, Fig. 3).

Fig. 3
figure 3

The Boxplot for CCTA-derived parameters in hospital readmission non-ischemic Heart failure. (A) CT derived epicardial adipose tissue; (B) CT derived cardiac muscle mass; (C) CT derived cardiac muscle mass to lumen volume; (D) peri-coronary adipose in the right coronary artery; (E) peri-coronary adipose in the left anterior descending artery; (F) peri-coronary adipose in the circumflex branch of the left coronary artery

In univariate logistic regression analysis, the CCTA-derived parameters were associated with readmission among non-ICHF patients (Table 3). The results of multivariate logistic regression indicated that EAT (OR: 1.15, 95% CI: 0.97–1.52, P = 0.012), PCATRCA (OR: 1.44, 95% CI: 1.02–2.52, P = 0.035), PCATLAD (OR: 3.1, 95% CI: 1.19–2.22, P = 0.014), PCATLCX (OR: 2.68, 95% CI: 1.52–8.91, P < 0.01) and CTV/M (OR: 0.25, 95% CI: 0.03–0.74, P = 0.008).

Prevalence of variables associated with outcome

Based on the maxstat, the optimal cut-off values of the EAT and V/M were 54.49 cm3 (HR: 1.05; 95% CI: 1.03–1.07; P < 0.001) and 20%, respectively, (HR: 0.59; 95% CI: 0.48–0.72; P < 0.001) for the endpoint of readmission (Supplementary Fig. 2). Patients with a PCATRCA lower than -64.68 HU were more likely to be readmitted (HR: 1.1; 95% CI: 1.03–1.16; P = 0.002). Patients with PCATLAD < -81.07 HU (HR: 1.3; 95% CI: 1.1–1.53; P = 0.002) or PCATLCX < -73.42 HU (HR: 1.33; 95% CI: 1.18–1.51; P < 0.001) were at the greatest risk of readmission (Table 4, Fig. 4).

Table 4 The K-M Analysis for the Prediction of non-ischemic HF readmission
Fig. 4
figure 4

The K-M survival curves of CCTA-derived parameters. (A) CTEAT; (B) CTM; (C) CTV/M; (D) PCATLAD; (E) PCATRCA; (F) PCATLCX

To test whether the CT-derived parameters had predictive value, we compared the AUC curves for 1) the model based on sex and Ao Asc (AUC = 0.71 95% CI: 0.601–0.812); 2) the model based on CTEAT, PCATRCA, PCATLAD, PCATLCX, CTM and CTV/M (AUC = 0.768, 95% CI: 0.67–0.86); and 3) the combined model (AUC = 0.819, 95% CI: 0.735–0.904) (Table 5, Fig. 5). The CCTA-derived volume was not associated with an increased risk of readmission.

Table 5 Results of ROC calculations for the three model
Fig. 5 
figure 5

The ROC analysis for the three models

Discussion

Several novel and clinically important findings were identified in our study. First, higher epicardial fat and pericoronary adipose tissue were observed in non-ICHF patients than in healthy people and are related to readmission. Second, the coronary lumen volume to cardiac muscle mass ratio (V/M) and left ventricular myocardial mass (M) were independently associated with poor prognosis. To our knowledge, this is the first reporting of the associations between readmission and risk factors, including CTEAT, CTV/M, CTFFR and CTPCAT, for non-ICHF.

Obesity is prevalent among HF patients, but whether it can serve as an independent prognostic factor for HF outcomes remains a topic of debate [21]. Some relevant studies have suggested that overweight and mildly to moderately obese HF patients may have longer survival rates than normal-weight or underweight HF patients [22, 23]. The potential explanation for this difference may be related to whether fat leads to excessive accumulation of WAT and myocardial fat infiltration [24, 25]. We hypothesized that, compared to BMI, EAT might be a more precise indicator of obesity, and thus, measuring EAT could better assess disease progression. In fact, previous research has confirmed that EAT participates in the occurrence and progression of cardiovascular diseases by synthesizing and secreting proinflammatory mediators and neurohormones such as IL-1β, IL-6, TNF-α, MCP-1, resistin, visfatin, and others [23]. Currently, research on EAT in HF patients is primarily focused on those classified based on ejection fraction (EF) [26]. Studies by Van Woerden and Obokata, utilizing cardiac magnetic resonance (CMR) and echocardiography, found that EAT mass and thickness were higher in HF with preserved EF (HFpEF) patients than in normal individuals [27, 28]. Liu et al. discovered that measuring EAT volume and density via CT scans could better predict the prognosis of patients with HFpEF [29]. Therefore, our research findings, based on CCTA, reveal that EAT levels in patients with non-ICHF often exceed those in individuals with normal cardiac function, consistent with the aforementioned research results.

In our study, we observed that patients with elevated amounts of EAT often exhibited an increased risk of readmission. The potential explanation for this observation may be linked to the detrimental impact of EAT on cardiac structure and function through endocrine and vascular pathways. Research indicates that the accumulation of EAT can lead to increased myocardial wall thickness, affecting cardiac diastolic function [30]. Excessive fat accumulation in EAT may also promote inflammation by releasing proinflammatory cytokines (such as adiponectin) into adjacent myocardium, thereby fostering microvascular inflammation and fibrosis, ultimately leading to atrial and ventricular functional impairment and, consequently, myocardial dysfunction [31].

The results of the NeXtsTeps (NXT) trial support the hypothesis of vascular involvement and microcirculatory dysfunction in HF patients [32]. However, in clinical practice, the assessment of microcirculatory dysfunction following myocardial infarction is often conducted using CMR or CT myocardial perfusion imaging, which requires advanced equipment and significant patient cooperation, making it challenging for widespread implementation [6]. In this study, we employed the ratio of coronary artery luminal volume to left ventricular mass, which reflects cardiac metabolic capacity without the need for additional pharmacological agents. Among patients with microvascular angina (MVA), V/M, as an indicator of myocardial oxygen demand, was significantly lower than in healthy individuals, consistent with the findings of this study [33]. We hypothesize that in readmitted non-ICHF patients, the microvessels undergo endothelium-dependent vasodilatory injury, thereby impacting myocardial contractility and relaxation. Calculating CTV/M can further enhance our understanding of the pathophysiology of cardiovascular diseases.

Inflammation has been shown to be linked to the development, progression, and worsening prognosis of HF [33, 34]. Assessment of pericoronary adipose tissue is a validated technique that allows the direct detection of coronary inflammation on CCTA. Proinflammatory cytokines released by the coronary arteries in the presence of inflammation inhibit adipocyte maturation, leading to smaller adipocytes and more water content [18]. In the present study, we found for the first time a reduction in the pericoronary fat attenuation coefficient, especially in the RCA, in readmitted patients with non-ICHF. Interestingly, the CRISP-CT study also demonstrated that the fat attenuation index (FAI) around the RCA can effectively predict myocardial ischaemia and the risk of high-risk plaques in coronary artery disease patients [33, 35].

Hence, our research outcomes are in concurrence with the predominant body of studies, demonstrating that elevated amounts of PCATs are linked to elevated risk, particularly in the context of the RCA. This phenomenon may be due to the accumulation of adipose tissue surrounding the coronary arteries, resulting in microvascular circulatory blockage, which influences the prognostic outlook of individuals suffering from non-ICHF.

In conclusion, CTEAT, CTPCAT, and CTV/M are associated with the prognosis of non-ICHF patients. Therefore, based on the above parameters, CCTA provides a multilevel and noninvasive evaluation option for non-ICHF patients that combines information on coronary arteries, cardiac structure and function, cardiac fat and other aspects and is expected to become a powerful tool for prognostic assessment and treatment decision-making, providing new hope for improving the management and prognosis of non-ICHF patients.

Limitations

First, this was a retrospective study conducted at a single centre, with relatively small sample sizes in each cohort, constraining the data we could analyse, such as waist circumference, waist-to-hip ratio, and educational level. Second, we focused solely on the presence of non-ICHF and did not perform subgroup analyses based on EF among HF patients. Third, due to radiation dose considerations in this study, we opted for prospectively collected CCTA data, limiting the assessment to EAT, PCAT, and V/M measurements. Comprehensive cardiac functional analysis based on the entire cardiac cycle was not feasible. Fourth, regarding the determination of optimal cut-off values, we chose to utilize maxstat for the calculations [20]. Compared to ROC curves, the maxstat method is typically less influenced by subjective judgements since it relies on the inherent statistical properties of the data, obviating the need for researcher subjectivity.

Future directions

In the future, a well-designed prospective study incorporating a larger number of cases and patient data will be essential to corroborate our findings. Additionally, long-term follow-up data can help enhance and refine model construction.