Characterization of biomarker profiles in patients with coronary artery disease: A prospective coronary computed tomography angiography study

HIGHLIGHTS

  • Patients with different stages of coronary artery disease exhibit significant differences in novel biomarkers.

  • Lp(a) appears to play a key role in early plaque development in patients without coronary calcification.

  • Machine-learning-based high-risk signatures predict CAD severity.

ABSTRACT

Background and Aims

Conventional biomarkers such as low-density lipoprotein (LDL) and high-density lipoprotein may fail to identify patients’ risk for significant coronary artery disease (CAD). This study evaluates the associations between multiple biomarkers and different CAD phenotypes, exploring a machine-learning biomarker-based patient clustering.

Methods and Results

We included 787 patients on primary prevention from the prospective ACTION registry (January 2024 to June 2025). All patients underwent coronary CTA and simultaneous biomarker testing, including LDL, high-density lipoprotein, triglyceride, apolipoprotein A-1, apolipoprotein B, lipoprotein(a) [Lp(a)], glycated hemoglobin (HbA1c), and high-sensitivity C-reactive protein. Of 382 patients (48.5%) with coronary artery calcium = 0, 42 (11%) had coronary plaque. These patients showed higher Lp(a) vs those without plaque (16.5 vs 11.5, P =.030), despite comparable SCORE2 risk (3.5% vs 3.0%, P =.284). Three biomarker-defined groups were identified after a machine learning unsupervised clustering: Cluster 1 had a favorable lipid profile with the lowest prevalence of CAD-Reporting and Data System (RADS) ≥ 3 (9.9%). Cluster 2 and 3, despite their significant intercluster differences in terms of Lp(a), LDL, and HbA1c levels, both showed a significantly higher prevalence of CAD-RADS ≥ 3 compared to cluster 1 (respectively 21.8% and 17.9%; vs cluster 1, P =.001). High-risk biomarker signatures were significantly associated to the prevalence of CAD-RADS ≥ 3, independently from the baseline SCORE2 (adjusted odds ratio 2.25; 95% confidence interval 1.32-3.82).

Conclusions

Distinct biomarker signatures associate with distinct CAD prevalence and severity that conventional lipid markers fail to distinguish. Lp(a) appears relevant for early plaque detection in coronary artery calcium = 0 patients. A comprehensive biomarker evaluation may help identifying high-risk subgroups overlooked by a conventional assessment.

Graphical abstract

Abbreviations

ACS, Acute Coronary Syndrome; AI, Artificial Intelligence;ApoA-1, Apolipoprotein A-1; ApoB, Apolipoprotein B; CAC, Coronary Artery Calcium;CAD, Coronary Artery Disease; CAD-RADS, Coronary Artery Disease- Reporting and Data System; CI, Confidence Interval; CTA, Computed Tomography Angiography; eGFR, Estimated Glomerular Filtration Rate; HbA1c, Glycated Hemoglobin; HDL, High-Density Lipoprotein;hs-CRP, High-Sensitivity C-Reactive Protein; ICA, Invasive Coronary Angiography; LDL, Low-Density Lipoprotein; Lp(a), Lipoprotein(a); MI, Myocardial Infarction; OR, Odds Ratio;PCI, Percutaneous Coronary Intervention; SCORE2, SCORE2-OP, SCORE2-DM, Systematic Coronary Risk Evaluation.

Introduction

Coronary artery disease (CAD) is one of the world’s major causes of morbidity and mortality, according to the most recent World Health Organization reports. Necessary actions are being taken to identify coronary plaques early, especially those at high-risk of future cardiac events through noninvasive advanced coronary computed tomography angiography (CTA). Still, the selection of the initial noninvasive diagnostic test should be based on the pretest likelihood of obstructive CAD. Despite the potential of coronary CTA to investigate multiple morphological and functional features of CAD, coronary artery calcium (CAC) score still acts as one of the main measures for the total atherosclerotic burden. Numerous studies have demonstrated its prognostic importance and correlated its value with blood biomarkers linked to the risk of future cardiovascular (CV) events. ,,, Nonetheless, about one fourth of CV events involve patients with a zero CAC. These two opposing patient groups—those with and without coronary calcification—may differ substantially in the levels of biomarkers associated to the development and progression of CAD, especially if biomarkers are limited to low-density lipoprotein (LDL) or high-density lipoprotein (HDL) cholesterol. These conventional biomarkers are typically assessed before deciding whether to perform a coronary CTA, but they may potentially mislead clinical decisions, providing a false sense of reassurance.

Relying solely on serial measurements of conventional biomarkers, such as LDL-cholesterol, to achieve predefined targets may no longer be sufficient to fully monitor CV risk or even treatment efficacy. , Patients with low LDL still have a residual CV risk that can be identified by other biomarkers (eg, hs-CRP). Therefore, less conventional biomarkers, such as Lipoprotein(a) [Lp(a)], apolipoprotein B (ApoB), and the ApoB/apolipoprotein A-1 (ApoA-1) ratio, offer promise for improved risk prediction, although their routine clinical use remains limited. These markers are part of a complex network that also involves inflammation and glucose metabolism, and studies characterizing a comprehensive biomarker “signature” for detecting significant CAD at early stage of prevention using coronary CTA are currently lacking. This study is part of an ongoing single-center prospective registry using coronary CTA (ACTION registry) to improve the management of patients with suspected CAD, and examine the association between biomarkers of lipid metabolism, glucose metabolism, and inflammation with different phenotypes of CAD.

Methods

Study objectives

This study aims to characterize the association between multiple biomarkers related to CV event risk and the CAD phenotype, as defined by the CAC score and the CAD-Reporting and Data System (CAD-RADS) score. Furthermore, this study aims to define whether a comprehensive biomarker signature has the potential to add incremental predictive value to the presence of significant CAD at the stage of primary prevention, before undergoing coronary CTA.

Study design

This study is part of a prospective, ongoing, single-center registry on the use of advanced CTA-based imaging analysis for the improvement of patient’s care (ACTION registry).

The registry is currently enrolling an all-comers population with low-intermediate pretest probability referred by general cardiologists for suspected CAD. Patients with a past-medical history of CV events are not excluded from the ACTION registry, but this represented an exclusion criterion for the present analysis. Patients with low-quality CT scans were excluded from the ACTION registry.

All patients enrolled in the registry provided written informed consent and agreed to the research use of their data.

Each consenting patient completed a standardized baseline questionnaire and had blood drawn immediately before undergoing coronary CTA. Blood samples were stored for assessment of the following parameters: lipid profile [(total cholesterol, triglycerides, HDL, LDL, ApoA-1, ApoB, Lp(a)], glycated hemoglobin [HbA1c], high-sensitivity C-reactive protein [hs-CRP]), creatinine, electrolytes, and complete blood count. Estimated glomerular filtration rate (eGFR) was calculated using the chronic kidney disease-EPI 2021 equation. These biomarkers were selected for two main reasons: first, they are widely available, standardized, inexpensive, and are already integrated or have strong potential to be integrated into CV risk stratification in primary prevention. Second, they capture the three principal metabolic drivers of atherosclerosis: impaired glucose metabolism, disordered lipid metabolism, and systemic inflammation. eGFR was added as an additional adjustment factor, given the well-established impact of chronic kidney disease on atherosclerosis and vascular calcification. The ACTION registry received ethical approval from the local ethics committee (Identification code: CA2792).

Data collection and processing

Data were prospectively collected by the Cardiovascular Research Team in the Clinical Research Facility in a high-volume University Hospital. This study included consecutive patients enrolled in the ACTION registry from January 2024 to June 2025. Demographics and clinical data of each patient, including—among others—age, gender, ethnicity, past medical history, reason for referral by cardiologist, and CV risk factors were prospectively collected by the Cardiovascular Research Team members using a standard questionnaire, as well as vital parameters at the time of the CTA scan, including noninvasive blood pressure. One of the cardiologist researchers of the Cardiovascular Research Team had to be always present during the questionnaire administration to ensure that the clinical data were correctly entered. The SCORE2, SCORE2-OP, and SCORE2-DM were calculated for each patient, respectively for nondiabetic patients <70 years old, patient >70 years old, and diabetic patients. The resulting mean/median percentage of estimated 10-year risk was reported for all three subgroups under the generic term “SCORE2” in the tables. A SCORE2<2.5% ≥2.5% but <7.5%, and ≥7.5%, were defined respectively as low, intermediate, and high-risk scores.

Questionnaires were securely stored in the Clinical Research Facility according to the General Data Protection Regulation (GDPR), and data were transferred into an electronic case report form on RedCap. Blood samples were collected within 30 minutes before the CTA. Samples were then processed on site for most of the tests, except for Lp(a), ApoA-1, and ApoB, which required samples to be sent and processed by Tallaght University Hospital Laboratories. Both laboratories have ISL15189:2012 accreditation.

CTA protocol and analysis

All coronary CTA scans were performed using a 256-slice multidetector Revolution Apex CT scanner (GE Healthcare, Waukesha, WI) with prospective ECG gating. Collimation of 256 × 0.625 mm with a z-coverage of 16 cm was used with a display field of view of 25 cm. Gantry rotation time was between 350 and 280 ms, dependent on patient size. Tube voltage was between 80 and 120 kV based on patient size (Auto Prescription, GE Healthcare, Waukesha, WI), and an adaptive tube current (SmartmA, GE Healthcare, Waukesha, WI) was used. Patient heart rates were reduced with oral and/or intravenous beta-blockers (Metoprolol) if average heart rate was over 55BPM on arrival to the imaging department, and provided no contraindications were present. Sublingual nitroglycerine (800 mcg) was given routinely immediately before the CT scan, provided no contraindications were present. An amount of 40 to 70 mL of nonionic iodinated contrast media, Iohexol 350 mg/mL (Omnipaque) was administered through an 18 g cannula at a rate of between 4.0 and 7 mL/sec, dependent on scan kV. A 40 mL saline bolus was administered following the contrast at the same infusion rate.

Bolus tracking (Smart Prep, GE Healthcare, Waukesha, WI) was used to monitor the contrast at the level of the carina, and the scan was started manually by the radiographer at peak opacification in the ascending aorta.

For patients with regular heart rates below 60 BPM, the scan was triggered at 75% of one R-R interval during a breath hold. For faster or irregular heart rates, a systolic acquisition at 45% of one R-R interval was also acquired. Both phases utilized padding to allow for digital motion artefact correction (Snap Shot Freeze, GE Healthcare).

All scans were reconstructed at 0.625 mm using a standard soft tissue kernel with DLIR (TrueFidelity, GE Healthcare) set to high.

Scans were analyzed by two highly experienced cardiologists. The scans were reported according to the CAD-RADS recommendations.

Statistical analysis

The study population was divided into two main subgroups—patients with CAC score = 0 and patients with CAC score >0. Differences within these two subgroups were tested between patients with or without plaque and patients with CAC score ≥ or <101, respectively. Continuous variables were tested for identifying normal distributions according to the Shapiro–Wilk test and Q-Q plot visual examination. They were compared with an unpaired Student t test or with the Mann–Whitney U nonparametric test, when appropriate. Data were expressed as means ± standard deviation or median and interquartile range. Categorical variables were compared using χ 2 or Fisher’s exact test as appropriate and expressed as frequencies and percentages. Continuous variables were compared across more than two groups using the Kruskal–Wallis test. When the overall test was statistically significant, pairwise comparisons were performed with Bonferroni-adjusted post hoc tests to account for multiple testing. An unsupervised k-means clustering algorithm was applied to standardized values of laboratory variables (HbA1c, Lp(a), apoA-1, apoB, ratio, LDL, HDL, triglycerides, hs-CRP, eGFR) to identify natural groupings within the study population. k-means is a way to group data into a number of clusters by finding “centers” in the average values of the data provided. The optimal number of clusters was set to three to ensure interpretability of the subsequent statistical comparisons and to avoid creating clusters that were too small. Quantitative clustering diagnostics, including elbow and silhouette analyses, supported K = 3 as a reasonable compromise between cluster compactness and interpretability. Inspection of the elbow plot revealed a clear inflection around K = 3, beyond which further increases in cluster number resulted in less relevant reductions in within-cluster variance (Supplementary Figure S1). The silhouette graph (Supplementary Figure S2) showed that an increased number of clusters would not have reduced clusters overlap, and, instead, would have just increased it for K = 5 or 6. Clustering stability was assessed using a resampling approach. The k-means clustering procedure ( K = 3) was repeated on multiple randomly selected subsets of the study population using the same biomarker variables. Agreement between original and resampled cluster assignments was quantified using the adjusted Rand index (ARI). After a repeated random subsampling generating 15 independent resamples (each resample included 250 randomly selected patients), the following results showed a moderate-to-good stability (mean ARI: 0.51; median ARI: 0.47; ARI min-max range: 0.30-0.83). k-means with K = 3 splits patients into 3 groups (clusters) based on overall similarity across the values of biomarkers. Each group has a centroid (the average profile). In other words, patients in one cluster are more alike than those in other groups, in terms of biomarker values. Each patient was assigned to one of the three resulting clusters for subsequent analyses. To simplify data interpretation, the binary values of CAD-RADS score ≥3 and CAC score ≥101 were chosen to explore the presence and severity of CAD in the three different clusters. In order to test the consistency of the results, a sensitivity analysis was performed, repeating the clustering and the subsequent analyses after excluding patients on lipid-lowering therapy at the time of the coronary CTA.

A multivariable logistic regression model was built to test the independent association of the cluster classification with the presence of CAD-RADS score ≥3.

Three other separate multivariable logistic regression models were built for the overall population and the subgroups of patients with a CAC score of 0 and those with a CAC score >0, to explore significant predictors of the presence of pure lipid plaque and the presence of a CAC score ≥101, respectively, using the single biomarkers.

Variable selection for multivariable model building was performed using the least absolute shrinkage and selection operator (LASSO) regression. All candidate laboratory variables were standardized prior to analysis. The penalty parameter ( λ ) was chosen via 10-fold cross-validation to minimize prediction error, and only variables with nonzero coefficients in the optimal model were retained for inclusion in the final multivariable logistic regression. Odds ratio with 95% confidence interval (CI) was used as the summary effect size. Missing values did not exceed 6% for any variable and occurred exclusively in laboratory results, due to rare errors in blood sample collection. The percentage of missing values for each laboratory variable is reported in Supplementary Figure S3.

For the k-means clustering analysis, listwise deletion was applied, and patients with any missing laboratory value were excluded. Patients excluded due to missing of complete laboratory values ( n = 91, 11.5%) showed no clinically meaningful differences in age, sex, BMI, SCORE2 risk, or prevalence of CAD-RADS ≥3 compared with included patients (Supplementary Table S1). Interobserver variability between the two expert analysts for the adjudication of CAD-RADS ≥ 3 was assessed in a sample of 35 patients, yielding a kappa coefficient of 0.83 (95% CI 0.69-0.98) and a percentage of agreement of 95%.

All data analyses were performed using STATANoW version 18.0 (Stata Corporation, LLC, College Station, TX).

Results

Between January 2024 and June 2025, a total of 856 patients were enrolled in the ACTION registry. A total of 69 patients were excluded since they had a past-medical history of acute CV events or prior coronary revascularization procedures, resulting into a final population of 787 patients. Among them, 382 patients (48.5%) had a CAC score equal to 0. Out of 382 patients with 0 CAC score, 42 (11%) showed the presence of a plaque. Only 5 patients out of 42 (11%) were categorized according to the CAD-RADS classification as having a CAD-RADS ≥3. Out of 405 patients with CAC score >0 (51.5%), 171 (42.2%) had a CAC score ≥101 and 123 patients (30.4%) had a CAD-RADS ≥3, with a significant difference between patients with CAC score between 1 and 100 and patients with a CAC score ≥101 (CAD-RADS ≥ 3: 16.0% vs 50.3% P =.000).

Difference in biomarkers between the four subgroups

Patients’ clinical characteristics across the four subgroups are presented in Table 1 , while biomarker values for these subgroups are shown in Table 2 . Patients with plaque despite a zero CAC score were comparable with patients without plaque regarding almost all the baseline clinical characteristics, including the median 10-year risk percentage according to the SCORE2. In terms of biomarkers concentration, patients with plaque showed higher HbA1c levels (39 vs 37 mmol/mol, P =.009), despite a similar percentage of diabetic patients, a higher Lp(a) concentrations (16.5 vs 11.5 mg/dL, P =.030), despite a comparable LDL and HDL level, and lower eGFR (88 vs 90 mL/min/1.73 m², P =.022) compared with those without plaque. In the subgroup with CAC score >0, patients with a CAC score≥101 ( n = 171) were older (64 vs 60 years, P <.001), more frequently male (70% vs 56%, P =.003), and had a higher SCORE2 compared with those with CAC score between 1 and 100. In terms of biomarkers, patients with CAC score ≥101 had lower HDL (1.3 vs 1.5 mmol/L, P =.009) and lower ApoA-1 (1.6 vs 1.7 g/L, P =.043) compared with those with CAC score between 1 and 100, with similar values of all the other biomarkers.

Table 1

General characteristics of the population stratified for the four subgroups of interest.

Patients with zero calcium score P -value Patients with coronary calcium P -value
Zero calcium and no plaque ( n = 340) Zero calcium and plaque ( n = 42) CAC score ≥1 and ≤100 ( n = 234) CAC score ≥101 ( n = 171)
Age (y) 53 (45-60) 56.0 (46.7-62.5) .888 60 (52-66) 64.0 (57.5-69.0) <.001
Height (cm) 169 (163-177) 169.0 (164.0-174.5) .596 172.0 (164.0-180.0) 173.0 (165.0-180.4) .403
Weight (kg) 79 (68-91) 80.2 (69.2-95.1) .481 83.4 (71.2-98.3) 87.8 (75.7-101.5) .873
BMI (kg/m 2) 27.3 (24.2-30.8) 28.3 (25.5-31.8) .122 28.0 (25.6-31.6) 29.0 (26.0-33.1) .057
Gender
Male 131 (38.5) 15 (35.7) .723 130 (55.6) 120 (70.2) .003
Female 209 (61.5) 27 (64.3) 104 (44.4) 51 (29.8)
Smoking (%) .555 .174
Current smoker 45 (13.3) 5 (12.2) 37 (15.8) 26 (15.2)
Nonsmoker 212 (62.7) 23 (56.2) 136 (58.1) 86 (50.3)
Ex-smoker 81 (24.0) 13 (31.7) 61 (26.1) 59 (34.5)
Hypertension 102 (30.2) 22 (52.4) .014 109 (46.6) 92 (53.8) .079
Diabetes 17 (5.0) 5 (12.2) .136 27 (11.6) 24 (14.0) .527
Family history of cardiovascular disease 261 (77.9) 31 (77.5) .417 166 (74.4) 123 (73.7) .194
Statins use 107 (31.4) 26 (61.9) .000 113 (48.3) 90 (52.6) .388
SCORE2 10-y risk (%) 3.0 (1.8-4.7) 3.5 (2.0-5.4) .284 4.9 (2.9-7.2) 6.1 (3.9-9.5) .001
CAD-RADS ≥ 3 5 (11) 37 (16.0) 86 (50.3) .000
Reasons for referral to coronary CTA
Chest pain 105 (31.0) 16 (38.0) .343 84 (35.9) 68 (39.7) .455
Primary prevention without chest pain symptoms 235 (69.0) 26 (62.0) .465 150 (64.1) 103 (60.2) .536

Bold values are statistically significant p <0.05.

Values are n (%), mean ± SD, or median (interquartile range).

BMI , body mass index; CAC , coronary artery calcium; CAD-RADS , Coronary Artery Disease-Reporting and Data System.

Table 2

Laboratory biomarker values in each of the four subgroups.

Patients with zero calcium score P -value Patients with coronary calcification P -value
Zero CAC score without plaque ( n = 340) Zero CAC score with plaque ( n = 42) CAC score ≥1 and ≤100 ( n = 234) CAC score ≥101 ( n = 171)
Hb1Ac (mmol/mol) 37 (35-39) 39 (36-41) .009 39 (37-41) 40.0 (37.0-43.5) .119
Hs-CRP (mg/L) 1.1 (0.6-2.6) 1.6 (1.1-2.5) .101 1.4 (0.7-2.7) 1.6 (0.9-3.7) .140
Total cholesterol (mmol/L) 4.9 (4.2-5.6) 5.0 (3.8-5.9) .255 4.8 (4.0-5.5) 4.4 (3.7-5.3) .069
Triglycerides (mmol/L) 1.1 (0.8-1.8) 1.2 (0.9-1.7) .736 1.4 (0.9-1.8) 1.5 (1.0-2.0) .092
HDL (mmol/L) 1.5 (1.2-1.9) 1.5 (1.2-1.9) .754 1.5 (1.2-1.7) 1.3 (1.1-1.6) .009
LDL (mmol/L) 2.7 (2.0-3.2) 2.7 (1.9-3.6) .426 2.6 (2.0-3.3) 2.3 (1.9-3.2) .076
Apo A-I (g/L) 1.6 (1.4-1.9) 1.5 (1.4-1.9) .713 1.7 (1.4-1.9) 1.6 (1.4-1.8) .043
Apo B (g/L) 0.9 (0.7-1.0) 0.9 (0.7-1.1) .868 0.9 (0.7-1.1) 0.9 (0.7-1.1) .746
Apo B/AboA-1 ratio 0.5 (0.4-0.6) 0.5 (0.4-0.7) .886 0.6 (0.4-0.7) 0.5 (0.4-0.7) .569
Lp (a) (mg/dL) 11.5 (8.8-38.0) 16.5 (8.8-86.6) .030 10.9 (8.8-51.9) 16.4 (8.8-65.0) .089
eGFR (mL/min/1.73 m²) 90 (81-90) 88 (74-90) .022 90.0 (77.5-90.0) 87 (76-90) .346
Only gold members can continue reading. Log In or Register to continue

Stay updated, free articles. Join our Telegram channel

Jun 27, 2026 | Posted by in CARDIOLOGY | Comments Off on Characterization of biomarker profiles in patients with coronary artery disease: A prospective coronary computed tomography angiography study

Full access? Get Clinical Tree

Get Clinical Tree app for offline access