This study was approved by the Institutional Review Board of the Affiliated Hospital of Nantong University (approval No. 2025-K179-01), and the requirement for informed consent was waived. The software and the equipment used are listed in the Table of Materials.
1. Study subjects and grouping
Between January 2016 and May 2023, data were collected for 752 patients with T2DM who underwent CCTA at the Affiliated Hospital of Nantong University. T2DM was diagnosed according to the Guidelines for the Prevention and Treatment of Type 2 Diabetes in China (2020 Edition). The exclusion criteria were as follows: (1) poor CCTA image quality; (2) congenital coronary anomalies, coronary stent implantation, or post-coronary artery bypass surgery; (3) pregnancy or hemoglobinopathy; (4) cardiogenic diseases other than coronary heart disease; (5) incomplete baseline data. A total of 185 subjects were included in the final analysis. A total of 360 diseased coronary arteries were analyzed, among which the major branches with coronary plaques included the left main (LM), left anterior descending (LAD), left circumflex artery (LCX), and right coronary artery (RCA). Clinical baseline data, including fasting plasma glucose (FPG), HbA1c, total cholesterol (TC), triglycerides (TG), high-density lipoprotein cholesterol (HDL-C), low-density lipoprotein cholesterol (LDL-C), and uric acid (UA), were collected within 48 h before and after the CCTA examination.
Since HGI is a linear regression residual, the predicted HbA1c can be generated by inputting FPG into the univariate linear regression equation between HbA1c and FPG in the study population. HGI is defined as the difference between the actual HbA1c and the predicted HbA1c. Patients and their affected coronary vessels were stratified into three groups according to HGI tertiles: H1 (low group) ≤-0.91, -0.91< H2 (medium group) <0.49, H3 (high group) ≥0.49. The study design and patient selection criteria are shown in Figure 1.
2. CCTA examination
CCTA examinations were performed using a 256-slice computed tomography (CT) scanner with the following protocol: detector collimation of 160 mm × 0.5 mm, gantry rotation speed of 0.28 s/rotation, slice thickness of 0.625 mm, tube current-time product of 350–635 mAs, and tube voltage of 100–120 kVp. All patients were scanned using retrospective electrocardiogram-gated technology. The scanning was performed from 1 cm inferior to the tracheal bifurcation to the diaphragmatic surface of the heart.
Contrast-enhanced imaging was performed following the intravenous administration of a non-ionic contrast agent into the brachial vein at a flow rate of 5 mL/s, immediately followed by a 20 mL saline flush delivered at the same rate. The dosage of the contrast agent was 0.8 mL per kilogram of body weight. For enhanced scans, contrast agent tracking triggering technology was employed, automatically detecting the CT value of the aortic root, with a triggering threshold of 120 Hounsfield units (HU).
3. Analysis of FAI
The original CCTA images were uploaded to the Deepwise Vascular Analysis Software (DW-CACTAS, Deepwise Corporation, China) by a radiologist with over three years of experience, who was blind to the patient’s clinical and CT information. The software automatically delineates the PCAT around the proximal 40 mm of LAD and LCX and 10–50 mm of RCA and allows for manual adjustments to the automatic depiction of the coronary artery walls, subsequently calculating the FAI for these three coronary arteries. PCAT was characterized as the adipose tissue situated radially at a distance equal to the coronary vessel’s diameter from its outer wall, setting its threshold value between -190 and -30 HU14. FAI was defined as the mean CT attenuation value of PCAT (Figure 2).
4. Quantitative analysis of CCTA plaque
Coronary atherosclerotic plaques refer to a type of tissue structure existing within or adjacent to the coronary artery lumen, with an area exceeding 1 mm2, distinct from pericardial tissue, epicardial fat tissue, and the vascular lumen15. CCTA images were transferred to a commercial workstation for plaque quantitative analysis. In CCTA images, the CT threshold values for different plaque components were as follows: lipid component CT value was set at -100–30 HU, fibrous component CT value at 31–130 HU, and calcified component CT value at 131–300 HU16, and the total plaque volume of each diseased coronary artery was measured.
The four CT vulnerability characteristics of high-risk coronary artery plaques were positive remodeling, low attenuation plaque, spotty calcification, and the napkin-ring sign. Coronary artery plaques that exhibited two or more of the above characteristics were defined as CT high-risk plaques17.
5. Statistical analysis
Categorical data were expressed as frequency and percentage [%], analyzed using the chi-square test, and further compared between groups using the row-by-column chi-square partition (corrected P = 0.0167). The Kolmogorov-Smirnov test was used to examine quantitative data normality. Normally distributed quantitative data were presented as mean ± standard deviation (SD), and statistical differences were assessed using analysis of variance (ANOVA). Quantitative data that were not normally distributed are presented as the median and interquartile range [M (Q1, Q3)], with the Kruskal-Wallis test applied for comparisons across multiple groups. Univariate and multivariate binary logistic regression analyses were conducted to identify independent risk factors associated with CT high-risk plaques. Statistical analysis was conducted using SPSS 26.0 software. A P-value of <0.05 was considered statistically significant.
6. Linear regression analysis of HGI with perivascular FAI
Multivariable linear regression was performed to assess the association between HGI (continuous) and perivascular FAI (LAD-FAI and LCX-FAI), adjusting for the same covariates as in the primary analysis.
7. Sensitivity analysis using HGI quartiles
To assess the robustness of the findings, HGI was categorized into quartiles, and its association with CT high-risk plaques was re-evaluated using multivariable logistic regression, adjusting for the same covariates as in the primary analysis.
8. Subgroup and interaction analyses
Subgroup analyses were performed according to age, gender, and glycemic control level using multivariable logistic regression, adjusting for the same covariates as in the primary analysis except the stratification variable itself. Interaction terms between HGI and each stratification variable were tested to assess effect modification.
9. Restricted cubic spline analysis
To explore potential dose-response and nonlinear relationships between HGI and the outcomes, restricted cubic spline (RCS) analyses were performed with HGI modeled as a continuous variable. For CT high-risk plaque, RCS logistic regression models were fitted, whereas for LAD-FAI, RCS linear regression models were used. All models were adjusted for the same covariates as in the primary analysis. The overall association and the nonlinear component were assessed using the Wald test.