
==== Front
Proteomes
Proteomes
proteomes
Proteomes
2227-7382
MDPI

10.3390/proteomes12030024
proteomes-12-00024
Article
Investigating the Prognostic Potential of Plasma ST2 in Patients with Peripheral Artery Disease: Identification and Evaluation
https://orcid.org/0000-0002-7191-1034
Li Ben Methodology Validation Formal analysis Investigation Writing – original draft Writing – review & editing 1234†
Shaikh Farah Methodology Validation Investigation Data curation Writing – review & editing 2†
Zamzam Abdelrahman Methodology Software Validation Formal analysis Investigation Writing – review & editing Visualization 2
Abdin Rawand Conceptualization Methodology Validation Investigation Resources Writing – review & editing Supervision Project administration 5
Qadura Mohammad Conceptualization Methodology Validation Investigation Resources Writing – review & editing Supervision Project administration Funding acquisition 1236*
1 Department of Surgery, University of Toronto, Toronto, ON M5S 1A1, Canada; benx.li@mail.utoronto.ca
2 Division of Vascular Surgery, St. Michael’s Hospital, Unity Health Toronto, University of Toronto, Toronto, ON M5B 1W8, Canada; farah.shaikh@unityhealth.to (F.S.); abdelrahman.zamzam@gmail.com (A.Z.)
3 Institute of Medical Science, University of Toronto, Toronto, ON M5S 1A1, Canada
4 Temerty Centre for Artificial Intelligence Research and Education in Medicine (T-CAIREM), University of Toronto, Toronto, ON M5S 1A1, Canada
5 Department of Medicine, McMaster University, Hamilton, ON L8S 4L8, Canada; rawand.abdin@medportal.ca
6 Li Ka Shing Knowledge Institute, St. Michael’s Hospital, Unity Health Toronto, University of Toronto, Toronto, ON M5B 1W8, Canada
* Correspondence: mohammad.qadura@utoronto.ca
† These authors contributed equally to this work.

29 8 2024
9 2024
12 3 2416 5 2024
27 8 2024
27 8 2024
© 2024 by the authors.
2024
https://creativecommons.org/licenses/by/4.0/ Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).
Soluble interleukin 1 receptor-like 1 (ST2) is a circulating protein demonstrated to be associated with cardiovascular diseases; however, it has not been studied as a biomarker for peripheral artery disease (PAD). Using a prospectively recruited cohort of 476 patients (312 with PAD and 164 without PAD), we conducted a prognostic study of PAD using clinical/biomarker data. Plasma concentrations of three circulating proteins [ST2, cytokine-responsive gene-2 (CRG-2), vascular endothelial growth factor (VEGF)] were measured at baseline and the cohort was followed for 2 years. The outcome of interest was a 2-year major adverse limb event (MALE; composite of major amputation, vascular intervention, or acute limb ischemia). Using 10-fold cross-validation, a random forest model was trained using clinical characteristics and plasma ST2 levels. The primary model evaluation metric was the F1 score. Out of the three circulating proteins analyzed, ST2 was the only one that was statistically significantly higher in individuals with PAD compared to patients without PAD (mean concentration in plasma of 9.57 [SD 5.86] vs. 11.39 [SD 6.43] pg/mL, p < 0.001). Over a 2-year period, 28 (9%) patients with PAD experienced MALE. Our predictive model, incorporating clinical features and plasma ST2 levels, achieved an F1 score of 0.713 for forecasting 2-year MALE outcomes. Patients identified as high-risk by this model showed a significantly increased likelihood of developing MALE (HR 1.06, 95% CI 1.02–1.13, p = 0.003). By combining clinical characteristics and plasma ST2 levels, our proposed predictive model offers accurate risk assessment for 2-year MALE in PAD patients. This algorithm supports risk stratification in PAD, guiding clinical decisions regarding further vascular evaluation, specialist referrals, and appropriate medical or surgical interventions, thereby potentially enhancing patient outcomes.

soluble interleukin 1 receptor-like 1
biomarkers
prognosis
peripheral artery disease
Blair FoundationThis research was funded by the Blair Foundation, funding number not applicable.
==== Body
pmc1. Introduction

Over 200 million people worldwide suffer from peripheral artery disease (PAD) due to atherosclerosis of the lower extremity arteries [1,2]. Although PAD is associated with significant limb loss and mortality, it remains undertreated [3]. An important contributor to this problem is the lack of validated and standardized prognostic biomarkers that can help identify high-risk patients and guide subsequent management [3]. Therefore, there is an important need to identify and validate novel biomarkers for PAD [4].

Soluble interleukin 1 receptor-like 1 (ST2), cytokine-responsive gene-2 (CRG-2), and vascular endothelial growth factor (VEGF) are circulating protein species that have been shown to be correlated with diseases of the cardiovascular system [5,6,7]. We selected the 3 aforementioned proteins for analysis because previous studies have demonstrated strong associations between these proteoforms and cardiovascular diseases, suggesting their potential relevance to PAD [5,6,7]. Specifically, Stojkovic et al. (2020) demonstrated associations between ST2 and platelet activation in PAD patients undergoing endovascular revascularization [5], while Horrevoets and colleagues showed that CRG fragments from endothelial cells were involved in the inflammatory reaction of atherosclerotic lesions [6]. Elsewhere, Ganta et al. (2021) demonstrated that inhibition of VEGF isoforms in ischemic muscle may promote perfusion recovery in preclinical PAD [7]. While previous research has demonstrated correlations between these proteins and cardiovascular conditions, few studies have assessed their prognostic implications for PAD [5,6,7]. Furthermore, prior studies have primarily focused on individual proteins, with no previous exploration of protein biomarkers in combination with clinical characteristics for prognostic purposes.

PAD is a complex disease involving multiple metabolic pathways that are related to different proteoforms; therefore, we hypothesized that biomarkers integrated with clinical features can enhance prognostic accuracy compared to the analysis of individual proteins alone [8]. By amalgamating biomarker data with clinical features linked to PAD outcomes, there is potential to develop highly accurate and specific predictive algorithms for adverse limb events associated with PAD [9,10,11]. The goal of this study was to identify a PAD-specific biomarker and integrate clinical and circulating protein species biomarker data through predictive modeling techniques to support PAD prognosis.

2. Materials and Methods

2.1. Ethics Approval

Approval for this study was granted by the research ethics board at Unity Health Toronto, University of Toronto, Canada on 2 August 2017 (approval code REB # 16-365). Prior to participation, informed consent was obtained from all patients, and all procedures adhered to the principles outlined in the Declaration of Helsinki [12].

2.2. Design

This prognostic study adhered to the Transparent Reporting of a Multivariable Prediction Model for Individual Prognosis or Diagnosis + Artificial Intelligence (TRIPOD + AI) statement [13], ensuring the transparent and comprehensive reporting of its findings.

2.3. Patient Recruitment

This study prospectively enrolled patients from ambulatory clinics at our institution between September 2020 and February 2022, including individuals with and without PAD. PAD was diagnosed based on an Ankle–Brachial Index (ABI) below 0.9 or a Toe Brachial Index (TBI) below 0.67, along with absent or diminished pedal pulses [14]. Patients with acute coronary syndrome, elevated troponin levels, or acute limb ischemia within the prior three months were excluded from the study.

2.4. Baseline Characteristics

Baseline characteristics included gender, age, dyslipidemia, hypertension, diabetes, smoking status, presence of congestive heart failure (CHF), coronary artery disease (CAD), history of stroke, and medication use including statins, acetylsalicylic acid (ASA), angiotensin II receptor blocker (ARB) or angiotensin-converting enzyme inhibitor (ACE-I), calcium channel blockers, beta-blockers, furosemide or hydrochlorothiazide, insulin, and oral antihyperglycemic agents. The definitions for cardiovascular risk factors and medications followed guidelines from the American College of Cardiology [15,16].

2.5. Quantification of Plasma Protein Levels

A phlebotomist collected blood samples from participants through the median cubital vein and centrifugation was used to isolate the plasma, which was then aliquoted and stored at −80 °C with no freeze–thaw cycles. On the day of analysis, the plasma was thawed to room temperature. Specifically, citrate plasma samples were used to prevent coagulation of the samples and to preserve the proteins of interest (ST2, CRG-2, and VEGF) for quantification. Using the LUMINEX assay (Bio-Techne, Minneapolis, MO, USA), in ref. [17], concentrations of three circulating proteins—ST2, CRG-2, and VEGF—were measured due to their involvement in metabolic pathways linked to atherosclerosis and cardiovascular diseases. Before analysis, the MagPix analyzer [18] was calibrated using bead kits from Fluidics Verification and Calibration (Luminex Corp; Austin, TX, USA) [19]. To ensure consistency, all tests were performed on the same day to minimize variability between assays. The inter- and intra-assay coefficients of variability for the samples were less than 10%. Luminex xPonent software version 4.3 [20] was used to evaluate a minimum of 50 beads per protein. The lower limit of detection (LLOD) and lower limit of quantification (LLOQ) of the Luminex multiplex assays were in the 0.1–1 pg/mL range for the proteins of interest [21]. Given that the mean plasma concentrations of the proteins of interest (ST2, CRG-2, and VEGF) ranged from 9.57 to 42.5 pg/mL in our cohort, all samples were measurable.

2.6. Follow-Up and Outcomes

Follow-up clinic visits were conducted at 12- and 24-months post-recruitment. The primary outcome of interest was major adverse limb events (MALE) occurring within 2 years of recruitment. MALE was defined as the need for major amputation of the lower extremity above the ankle, vascular intervention (either open or endovascular revascularization of the lower extremity arteries), or acute limb ischemia (sudden decrease in limb perfusion occurring within 14 days due to arterial thrombosis or embolism). Initial analysis revealed that all adverse limb events occurred exclusively in PAD patients. Therefore, prognostic models were developed solely within the PAD cohort.

2.7. Model Development and Evaluation

Random forest was chosen as the predictive algorithm, which is an ensemble model that works through several decision trees [22]. Decision trees can classify samples into branch-like segments and leverage data from multiple covariates to develop prediction algorithms for an outcome of interest [23]. Random forest can leverage complex datasets effectively because of its non-parametric nature [23]. This algorithm was chosen because it is well described in literature and performs well for predicting human health outcomes [24,25,26].

The dataset was partitioned randomly into 70% for training and 30% for testing purposes. Utilizing the random forest algorithm, which underwent training with 10-fold cross-validation, our aim was to predict 2-year MALE. The input features comprised clinical characteristics such as gender, age, diabetes, dyslipidemia, hypertension, smoking status, CAD, CHF, previous stroke, and medication use (including statins, ASA, ARB or ACE-I, beta-blockers, calcium channel blockers, furosemide or hydrochlorothiazide, insulin, and oral antihyperglycemic agents), alongside plasma ST2 levels. Post-training, the models were evaluated on unseen test set data. Variable importance scores (gain) were calculated to identify the most influential predictive factors, gauging their relative impact on prediction outcomes [27].

2.8. Statistical Analysis

We summarized the baseline characteristics using means and standard deviations (SDs) for continuous variables or numbers and proportions for categorical variables. Group differences were assessed using independent t-tests for continuous variables and chi-square tests for categorical variables. Protein levels were compared between PAD and non-PAD patients using independent t-tests. Proteins that showed differential expression in patients with versus without PAD were further utilized in model development. Two-year event rates were compared between patients with vs. without PAD using chi-square tests. Hazard ratios were calculated to determine the association between ST2 and adverse limb events controlling for all baseline characteristics. The F1 score was the primary model evaluation metric, which measures the harmonic mean of the precision and recall values by which ST2 levels predicted adverse limb events [28]. The F1 score, defined by the formula 2 × (precision × recall)/(precision + recall), ranges from 0 to 1, where 0 indicates minimal precision and/or recall, and 1 signifies optimal precision and recall [28]. Using our prognostic model, patients were stratified into low- or high-risk categories for developing 2-year MALE based on the Youden Index, which maximizes the model’s performance (sensitivity and specificity) through receiver operating characteristic (ROC) curve analysis [29]. To assess freedom from MALE over 2 years between low- and high-risk groups, Kaplan–Meier curves were plotted and compared using Cox proportional hazards analysis, adjusting for all baseline characteristics. This stratified analysis aimed to elucidate the clinical relevance of risk predictions derived from the prognostic model, offering insights into how the trajectories of low-risk versus high-risk patients differ in terms of MALE risk over a 2-year period. Statistical significance was set at a two-tailed p-value < 0.05. All statistical analyses were conducted using SPSS software version 23 (SPSS Inc., Chicago, IL, USA) [30].

3. Results

3.1. Patients

In total, 476 patients participated in this study, with 312 having PAD and 164 without PAD. Those with PAD were older (mean age 71 [SD 10] vs. 65 [SD 12] years; p < 0.001) and a higher percentage had CAD (38% vs. 21%; p < 0.001), diabetes (42% vs. 21%; p < 0.001), dyslipidemia (84% vs. 61%; p < 0.001), hypertension (82% vs. 59%; p < 0.001), and a history of stroke (16% vs. 8%; p = 0.011). They were also more likely to be current/past smokers (80% vs. 64%; p = 0.002). Additionally, individuals with PAD were more frequently prescribed risk-reduction medications such as ACE-I/ARB (66% vs. 45%, p = 0.001), statins (73% vs. 57%, p < 0.001), ASA (80% vs. 60%, p < 0.001), and beta blockers (41% vs. 30%, p = 0.001) (Table 1).

3.2. Plasma Protein Concentrations

Of the three proteins tested, ST2 was the only protein that was statistically significantly higher in individuals with PAD compared to patients without PAD: mean concentration in plasma of 11.39 [SD 6.43] vs. 9.57 [SD 5.86] pg/mL; p < 0.001 (Table 2). ST2 was therefore included in further analyses.

3.3. Adverse Limb Events

All adverse limb events observed during the 2-year follow-up period were exclusively found in patients with PAD: MALE (n = 28, 9%), major amputation (n = 17, 5%), and vascular intervention (n = 19, 6%). No patients developed acute limb ischemia (Table 3). There were significant associations between ST2 and adverse events over 2 years of follow-up, including MALE (HR 1.06 [95% CI 1.02–1.13], p = 0.005) and need for vascular intervention (HR 1.07 [95% CI 1.01–1.12], p = 0.003) (Table 4). Of note, Table 1 and Table 3 are similar to our previous publication because we used a similar patient cohort to investigate PAD biomarkers [31]. While our previous paper identified IL-7 as a myokine biomarker for PAD, the current study is innovative in that we discovered a novel PAD biomarker, ST2, which has a unique biological relationship with PAD through its involvement in various mechanistic pathways of atherosclerosis, thrombosis, and plaque vulnerability. Therefore, the current study adds to our previous work by unveiling a new PAD biomarker that can further guide our understanding of PAD pathophysiology and prognosis.

3.4. Model Performance

Using a combination of clinical features and plasma ST2 levels, the random forest model achieved an F1 score of 0.713 for predicting 2-year MALE. The most important predictive features for the model were (1) ST2, (2) age, (3) current/past smoking, (4) diabetes, (5) hypertension, (6) gender, (7) CAD, (8) dyslipidemia, and (9) CHF (Figure 1). Of note, ST2 was the most important predictive feature for 2-year MALE when evaluated alongside patient age, gender, and comorbidities as input features.

3.5. Risk Stratification Using Model

Youden’s Index calculation identified 0.60 as the optimal threshold for predicting 2-year MALE using our model. This threshold was employed to categorize our cohort into high and low-risk groups for adverse limb events. Over the 2-year follow-up period, patients classified as high-risk exhibited reduced freedom from MALE compared to those classified as low-risk (HR 1.06, 95% CI 1.02–1.13, p = 0.005) (Figure 2).

4. Discussion

4.1. Summary of Findings

In this study, we identified ST2 as a specific proteoform biomarker for PAD and developed a robust prognostic model combining clinical characteristics and plasma ST2 levels to accurately predict PAD prognosis. Several key findings emerged from our analysis. Firstly, among the three circulating protein species studied, ST2 was the sole proteoform found to be elevated in individuals with PAD compared to those without PAD. Although this study investigated multiple proteins as potential PAD biomarkers, only ST2 emerged as the singular protein that was significantly correlated with PAD-related adverse limb events. Therefore, ST2 was identified as the most robust PAD-specific prognostic biomarker. Secondly, our predictive model, integrating clinical features and plasma ST2 levels, demonstrated strong performance in forecasting PAD prognosis. Feature importance analysis underscored ST2 levels as the most significant predictor, highlighting the pivotal role of proteoform biomarkers in predicting outcomes in PAD. Although additional research is needed to confirm our findings, our study demonstrates that ST2 can be considered a PAD-specific biomarker, as it was the only protein that was significantly elevated in patients with PAD and predicted the development of adverse limb events over a 2-year follow-up period. Given the prognostic relevance of ST2, further research is warranted to explore its biological relationship with PAD development and progression, aiming to enhance our understanding of proteome complexity and inform targeted therapeutic strategies. Thirdly, utilizing our prognostic model, we stratified patients into low- and high-risk categories for adverse events. Kaplan–Meier analysis illustrated that patients identified as high-risk by our model were more susceptible to developing 2-year MALE compared to those classified as low-risk. This underscores the clinical utility of our model in aiding clinicians to predict the future risk trajectory of PAD patients concerning adverse limb events.

4.2. Comparison to Existing Literature

Ross and colleagues (2019) developed a prognostic model for predicting major adverse cardiac and cerebrovascular events in PAD patients using electronic health records data [32]. However, their model did not incorporate biomarker-based input features, despite the demonstrated impact of proteoform biomarkers on PAD prognosis in our study and others [5,6,7]. Our research addressed this limitation by integrating protein species biomarker data into our prediction models. We achieved strong performance in predicting 2-year MALE among PAD patients, with an F1 score of 0.713. Thus, our study underscores the importance of utilizing proteoform biomarker data to enhance prognostic models, potentially improving their predictive accuracy compared to models based solely on clinical features.

4.3. Explanation of Findings

There are several potential explanations for our findings. First, ST2 was found to be significantly elevated in patients with PAD in this study and was an important predictor of PAD prognosis. Discovered in 1989, the ST2 gene is found on chromosome 2q12 as part of the wider interleukin 1 (IL-1) gene cluster [33]. ST2 is the IL-33 receptor and facilitates downregulation of hypertrophy and fibrosis in tissues that are mechanically strained [34]. ST2 has been demonstrated to be associated with the pathogenesis of atherosclerosis and thrombosis in previous studies [35,36]. Specifically, in monocytes, ST2 proteoforms induce tissue factor expression and the release of prothrombotic vesicles [37]. Consequently, increased levels of IL-33/ST2 are associated with an increased risk of in-stent restenosis and poor outcomes in patients with CAD [38,39]. Similarly, in patients with carotid artery stenosis, ST2 levels have been found to be elevated and correlated with atherosclerotic plaque vulnerability [40]. ST2 has been shown to be a potential biomarker in individuals with heart failure and myocardial infarction as it is correlated with adverse cardiovascular events and death [41,42]. Taken together, these findings explain the potential mechanism by which ST2 is involved in PAD development. In comparison to our previously published work [43,44,45,46,47], we identified ST2 as a novel prognostic biomarker for PAD that has excellent predictive performance for adverse limb events in combination with clinical features. Importantly, ST2 has a strong biological relationship to PAD through its involvement in various mechanistic pathways of atherosclerosis, thrombosis, and plaque vulnerability. Therefore, there is potential to combine ST2 with our previously identified PAD biomarkers to form a comprehensive panel of proteins to support PAD management. This may form the basis for a larger future validation study. Secondly, our analysis demonstrated a high incidence of adverse limb events in PAD patients, with close to 10% of the cohort experiencing MALE over 2 years. This highlights the need for new strategies to reduce complications in this patient population, emphasizing the value of developing more effective prognostic tools. Third, our prognostic model achieved excellent performance for several potential reasons. Contrary to traditional statistics such as logistic regression, which assumes that there is a linear correlation between the covariates and the logit of the dependent variable, advanced modeling methods are not limited by the linearity assumption and can therefore better model non-linear relationships [48,49]. This is valuable in proteoform biomarker-based prognostic models, where different protein species may be involved in various biological pathways that may interact in complex ways because of proteome complexity to contribute to a disease condition [50]. Our random forest algorithm likely achieved good performance because it is an ensemble model that combines several decision trees [51]. This approach minimizes overfitting, reduces variance, and effectively handles complex datasets [51]. Our study highlights the advantages of developing a prognostic model that includes proteoform biomarker information, which may enhance model performance in comparison to solely using clinical information. Because PAD is a complex condition that involves multiple biological pathways, previous papers have shown the value of an integrated approach to support the prognosis of PAD [52]. This investigation affirms that by applying advanced predictive modeling methods to assess clinical information in conjunction with biomarker data, accurate PAD prognostic tools can be generated.

4.4. Implications

Our predictive models offer practical implications for guiding clinical decision-making across various scenarios. First, our algorithm can help screen individuals for asymptomatic PAD, which is especially valuable in family practice settings. General practitioners use the predictive model during their clinical assessments to understand a patient’s PAD risk [53]. Patients who screen positive for being at elevated risk of PAD-associated adverse outcomes may be sent for additional vascular evaluation, including arterial duplex ultrasound, to confirm a PAD diagnosis [54]. Patients identified as low-risk may receive ongoing treatment from their family physician with a focus on risk factor optimization such as lifestyle modifications and prescription of statins and ASA [55]. Patients predicted to be at elevated risk for MALE should be referred to a vascular surgeon for evaluation and treatment [56]. Vascular surgeons can integrate our model with their clinical judgment to identify individuals at heightened risk of adverse limb outcomes, potentially benefiting from the following: (1) additional imaging studies to evaluate vascular anatomy and disease severity [57], (2) low-dose rivaroxaban therapy [58], or (3) surgical interventions for limb salvage [59,60]. Our predictive model has the potential to enhance care for individuals with PAD across both generalist and specialist healthcare settings. It can facilitate PAD screening, risk stratification, and early identification of patients at increased risk for adverse limb outcomes. This approach may reduce unnecessary specialist referrals, improve PAD management outcomes, and decrease healthcare expenditures [61].

4.5. Limitations

Our study has several limitations. Firstly, recruitment occurred at a single center, and validation in other settings may help determine model generalizability. Notably, the non-PAD group had a smaller sample size compared to the PAD group because the focus of this study was to identify a prognostic biomarker for PAD patients. However, it may be prudent for future studies to include a larger non-PAD cohort as a control group. Secondly, the outcomes in our study were recorded based on a 2-year follow-up period. Additional follow-up may improve our understanding of the predictive ability of our model given that PAD is a chronic condition. Given that the goal of this study was to identify prognostic biomarkers for PAD, the correlation between ST2 and other cardiovascular diseases was not investigated. Thirdly, laboratory results such as lipid levels, creatinine, liver function tests, glucose, and coagulation parameters were not measured in this study. Future models incorporating these variables may further improve predictive performance. Fourthly, the protein species biomarkers evaluated in this paper are used predominantly in research environments. Additional work is required to demonstrate the clinical feasibility of including these proteoform biomarkers in the routine care of patients with PAD.

5. Conclusions

We identified ST2 as a protein species biomarker for PAD and used plasma concentrations of ST2 in combination with clinical characteristics to build an algorithm that accurately predicts PAD prognosis. Our model holds promise for PAD risk stratification, supporting prompt identification and targeted treatment of PAD. Specifically, the model can identify high-risk patients who may be referred for additional vascular work-up and may benefit from more aggressive medical and/or surgical management. Additionally, our paper highlights the need for translational research assessing the biological relationship between ST2 and PAD progression, offering potential advancements in our understanding of the underlying pathogenesis of PAD and informing targeted therapeutic strategies.

Author Contributions

Conceptualization, R.A. and M.Q.; methodology, B.L., F.S., A.Z., R.A. and M.Q.; software, A.Z.; validation, B.L., F.S., A.Z., R.A. and M.Q.; formal analysis, B.L. and A.Z.; investigation, B.L., F.S., A.Z., R.A. and M.Q.; resources, R.A. and M.Q.; data curation, F.S.; writing—original draft preparation, B.L.; writing—review and editing, B.L., F.S., A.Z., R.A. and M.Q.; visualization, A.Z.; supervision, R.A. and M.Q.; project administration, R.A. and M.Q.; funding acquisition, M.Q. All authors have read and agreed to the published version of the manuscript.

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki, and approved by the Institutional Review Board of Unity Health Toronto, University of Toronto, Canada on 2 August 2017 (approval code REB # 16-365).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding author.

Conflicts of Interest

The authors declare no conflicts of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Figure 1 Variable importance scores (gain) for the random forest model input features for prognosis of 2-year major adverse limb events in individuals with peripheral artery disease. Abbreviations: soluble interleukin 1 receptor-like 1 (ST2), congestive heart failure (CHF), coronary artery disease (CAD), diabetes mellitus (DM), hypercholesterolemia (HC), hypertension (HTN).

Figure 2 Kaplan–Meier analysis of the freedom from major adverse limb events among patients categorized as low versus high risk by the random forest model. Classification into these risk groups was based on a threshold of 0.60 derived from the Youden Index, optimizing the prediction model’s sensitivity and specificity through receiver operating characteristic curve analysis. Cox proportional hazards analysis adjusted for sex, age, dyslipidemia, hypertension, diabetes, smoking history, coronary artery disease, congestive heart failure, previous stroke, use of statins, acetylsalicylic acid, angiotensin II receptor blocker or angiotensin-converting enzyme inhibitor, calcium channel blocker, beta blocker, furosemide or hydrochlorothiazide, oral antihyperglycemic agent, and insulin. Abbreviations: CI (confidence interval), HR (hazard ratio).

proteomes-12-00024-t001_Table 1 Table 1 Baseline characteristics.

	PAD
(n = 312)	Non-PAD
(n = 164)	p	
Age, mean (SD)	71 (10)	65 (12)	<0.001	
Female sex	109 (35)	67 (41)	0.204	
Hypertension	257 (82)	96 (59)	<0.001	
Diabetes	131 (42)	34 (21)	<0.001	
Dyslipidemia	263 (84)	100 (61)	<0.001	
Current smoking	78 (25)	35 (21)	0.002	
Past smoking	171 (55)	71 (43)	0.001	
Coronary artery disease	118 (38)	34 (21)	<0.001	
Congestive heart failure	11 (4)	4 (2)	0.519	
Previous stroke	51 (16)	13 (8)	0.011	
Statin	229 (73)	93 (57)	<0.001	
Acetylsalicylic acid	251 (80)	99 (60)	<0.001	
Beta blocker	134 (41)	50 (30)	0.001	
ACE-I/ARB	216 (66)	74 (45)	0.001	
Calcium channel blocker	82 (25)	34 (21)	0.079	
Hydrochlorothiazide or furosemide	41 (13)	17 (10)	0.190	
Insulin	22 (7)	6 (4)	0.255	
Oral antihyperglycemic agent	24 (8)	8 (5)	0.201	
Values are presented as number (percentage) unless specified otherwise. ARB (angiotensin II receptor blocker), ACE-I (angiotensin-converting enzyme inhibitor), PAD (peripheral artery disease), SD (standard deviation).

proteomes-12-00024-t002_Table 2 Table 2 Protein plasma concentrations in individuals with vs. without peripheral artery disease.

	Non-PAD
(n = 164)	PAD
(n = 312)		
	Mean	Standard
Deviation	Mean	Standard
Deviation	p	
ST2	9.57	5.86	11.39	6.43	<0.001	
CRG-2	38.56	28.56	42.5	30.46	0.085	
VEGF	16.61	12.24	20.81	26.96	0.096	
Protein concentrations reported in pg/mL. Abbreviations: soluble interleukin 1 receptor-like 1 (ST2), cytokine-responsive gene-2 (CRG-2), vascular endothelial growth factor (VEGF).

proteomes-12-00024-t003_Table 3 Table 3 Two-year adverse limb events.

	PAD
(n = 312)	Non-PAD
(n = 164)	p	
Major adverse limb event	28 (9)	0 (0)	0.001	
Major amputation	17 (5)	0 (0)	0.002	
Vascular intervention	19 (6)	0 (0)	0.001	
Acute limb ischemia	0 (0)	0 (0)	N/A	
Values are presented as number (percentage) unless specified otherwise. PAD (peripheral artery disease).

proteomes-12-00024-t004_Table 4 Table 4 Adjusted hazard ratios for associations between soluble interleukin 1 receptor-like 1 (ST2) and 2-year major adverse limb events in patients with peripheral artery disease.

	Hazard Ratio [95% CI] *	p-Value	
Major adverse limb event	1.06 [1.02–1.13]	0.005	
Vascular intervention	1.07 [1.01–1.12]	0.003	
Major amputation	1.00 [0.99–1.11]	0.084	
* Adjusted for age, sex, hypertension, dyslipidemia, diabetes, past/current smoking, congestive heart failure, coronary artery disease, previous stroke, acetylsalicylic acid, statin, angiotensin-converting enzyme inhibitor or angiotensin II receptor blocker, beta blocker, calcium channel blocker, hydrochlorothiazide or furosemide, oral antihyperglycemic agent, and insulin.

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
==== Refs
References

1. Zemaitis M.R. Boll J.M. Dreyer M.A. Peripheral Arterial Disease StatPearls StatPearls Publishing Treasure Island, FL, USA 2021
2. Olin J.W. Sealove B.A. Peripheral Artery Disease: Current Insight into the Disease and Its Diagnosis and Management Mayo Clin. Proc. 2010 85 678 692 10.4065/mcp.2010.0133 20592174
3. Mehta A. Dhindsa D.S. Hooda A. Nayak A. Massad C.S. Rao B. Makue L.F. Rajani R.R. Alabi O. Quyyumi A.A. Premature Atherosclerotic Peripheral Artery Disease: An Underrecognized and Undertreated Disorder with a Rising Global Prevalence Trends Cardiovasc. Med. 2021 31 351 358 10.1016/j.tcm.2020.06.005 32565142
4. Zamzam A. Syed M.H. Rotstein O.D. Eikelboom J. Klein D.J. Singh K.K. Abdin R. Qadura M. Validating Fatty Acid Binding Protein 3 as a Diagnostic and Prognostic Biomarker for Peripheral Arterial Disease: A Three-Year Prospective Follow-up Study eClinicalMedicine 2023 55 101766 10.1016/j.eclinm.2022.101766 36531981
5. Stojkovic S. Demyanets S. Kopp C.W. Hengstenberg C. Wojta J. Eichelberger B. Panzer S. Gremmel T. Association of Soluble Suppression of Tumorigenesis 2 (sST2) with Platelet Activation, Monocyte Tissue Factor and Ischemic Outcomes Following Angioplasty and Stenting Front. Cardiovasc. Med. 2020 7 605669 10.3389/fcvm.2020.605669 33415128
6. Horrevoets A.J. Fontijn R.D. van Zonneveld A.J. de Vries C.J. ten Cate J.W. Pannekoek H. Vascular Endothelial Genes That Are Responsive to Tumor Necrosis Factor-Alpha in Vitro Are Expressed in Atherosclerotic Lesions, Including Inhibitor of Apoptosis Protein-1, Stannin, and Two Novel Genes Blood 1999 93 3418 3431 10.1182/blood.V93.10.3418.410k23_3418_3431 10233894
7. Ganta V.C. Annex B.H. Peripheral Vascular Disease: Preclinical Models and Emerging Therapeutic Targeting of the Vascular Endothelial Growth Factor Ligand-Receptor System Expert Opin. Ther. Targets 2021 25 381 391 10.1080/14728222.2021.1940139 34098826
8. Parvar S.L. Fitridge R. Dawson J. Nicholls S.J. Medical and Lifestyle Management of Peripheral Arterial Disease J. Vasc. Surg. 2018 68 1595 1606 10.1016/j.jvs.2018.07.027 30360849
9. Cornejo Del Río V. Mostaza J. Lahoz C. Sánchez-Arroyo V. Sabín C. López S. Patrón P. Fernández-García P. Fernández-Puntero B. Vicent D. Prevalence of Peripheral Artery Disease (PAD) and Factors Associated: An Epidemiological Analysis from the Population-Based Screening PRE-Diabetes and Type 2 DIAbetes (SPREDIA-2) Study PLoS ONE 2017 12 e0186220 10.1371/journal.pone.0186220 29073236
10. Duval S. Massaro J.M. Jaff M.R. Boden W.E. Alberts M.J. Califf R.M. Eagle K.A. D’Agostino R.B. Pedley A. Fonarow G.C. An Evidence-Based Score to Detect Prevalent Peripheral Artery Disease (PAD) Vasc. Med. 2012 17 342 351 10.1177/1358863X12445102 22711750
11. Gouda P. Ramasundarahettige C. Anand S. Muhlhoffer E. Berkowitz S. Fox K.A. Eikelboom J. Welsh R. Clinical Factors Associated with Peripheral Artery Disease in Patients with Documented Coronary Artery Disease: A Post Hoc Analysis of the COMPASS Trial Atherosclerosis 2021 331 38 44 10.1016/j.atherosclerosis.2021.07.003 34340829
12. World Medical Association World Medical Association Declaration of Helsinki: Ethical Principles for Medical Research Involving Human Subjects JAMA 2013 310 2191 2194 10.1001/jama.2013.281053 24141714
13. Collins G.S. Moons K.G.M. Dhiman P. Riley R.D. Beam A.L. Calster B.V. Ghassemi M. Liu X. Reitsma J.B. van Smeden M. TRIPOD + AI Statement: Updated Guidance for Reporting Clinical Prediction Models That Use Regression or Machine Learning Methods BMJ 2024 385 e078378 10.1136/bmj-2023-078378 38626948
14. Gul F. Janzer S.F. Peripheral Vascular Disease StatPearls StatPearls Publishing Treasure Island, FL, USA 2021
15. Grundy S.M. Stone N.J. Bailey A.L. Beam C. Birtcher K.K. Blumenthal R.S. Braun L.T. de Ferranti S. Faiella-Tommasino J. Forman D.E. 2018 AHA/ACC/AACVPR/AAPA/ABC/ACPM/ADA/AGS/APhA/ASPC/NLA/PCNA Guideline on the Management of Blood Cholesterol J. Am. Coll. Cardiol. 2019 73 e285 e350 10.1016/j.jacc.2018.11.003 30423393
16. Whelton P.K. Carey R.M. Aronow W.S. Casey D.E. Collins K.J. Dennison H.C. DePalma S.M. Gidding S. Jamerson K.A. Jones D.W. 2017 ACC/AHA/AAPA/ABC/ACPM/AGS/APhA/ASH/ASPC/NMA/PCNA Guideline for the Prevention, Detection, Evaluation, and Management of High Blood Pressure in Adults J. Am. Coll. Cardiol. 2018 71 e127 e248 10.1016/j.jacc.2017.11.006 29146535
17. Luminex Assays, Multiplex Immunoassays Available online: https://www.bio-techne.com/ (accessed on 6 May 2023)
18. MAGPIX® System|xMAP Instrument|Luminex Corporation Available online: https://www.luminexcorp.com/magpix-system/ (accessed on 18 December 2021)
19. Luminex Assays—CA Available online: https://www.thermofisher.com/ca/en/home/life-science/antibodies/immunoassays/procartaplex-assays-luminex.html (accessed on 18 December 2021)
20. xPONENT® Software for xMAP® Instruments Luminex Corporation Austin, TX, USA 2024
21. Understanding LLD, LLOQ, and ULOQ in Luminex Multiplex Cytokine Assays—Creative Proteomics Available online: https://cytokine.creative-proteomics.com/llq-lloq-uloq-luminex-cytokine-assays.htm (accessed on 15 July 2024)
22. Rigatti S.J. Random Forest J. Insur. Med. 2017 47 31 39 10.17849/insm-47-01-31-39.1 28836909
23. Song Y.-Y. Lu Y. Decision Tree Methods: Applications for Classification and Prediction Shanghai Arch. Psychiatry 2015 27 130 135 10.11919/j.issn.1002-0829.215044 26120265
24. Elfanagely O. Toyoda Y. Othman S. Mellia J.A. Basta M. Liu T. Kording K. Ungar L. Fischer J.P. Machine Learning and Surgical Outcomes Prediction: A Systematic Review J. Surg. Res. 2021 264 346 361 10.1016/j.jss.2021.02.045 33848833
25. Bektaş M. Tuynman J.B. Costa Pereira J. Burchell G.L. van der Peet D.L. Machine Learning Algorithms for Predicting Surgical Outcomes after Colorectal Surgery: A Systematic Review World J. Surg. 2022 46 3100 3110 10.1007/s00268-022-06728-1 36109367
26. Senders J.T. Staples P.C. Karhade A.V. Zaki M.M. Gormley W.B. Broekman M.L.D. Smith T.R. Arnaout O. Machine Learning and Neurosurgical Outcome Prediction: A Systematic Review World Neurosurg. 2018 109 476 486.e1 10.1016/j.wneu.2017.09.149 28986230
27. Loh W.-Y. Zhou P. Variable Importance Scores J. Data Sci. 2021 19 569 592 10.6339/21-JDS1023
28. Hicks S.A. Strümke I. Thambawita V. Hammou M. Riegler M.A. Halvorsen P. Parasa S. On Evaluation Metrics for Medical Applications of Artificial Intelligence Sci. Rep. 2022 12 5979 10.1038/s41598-022-09954-8 35395867
29. Ruopp M.D. Perkins N.J. Whitcomb B.W. Schisterman E.F. Youden Index and Optimal Cut-Point Estimated from Observations Affected by a Lower Limit of Detection Biom. J. 2008 50 419 430 10.1002/bimj.200710415 18435502
30. SPSS Software Available online: https://www.ibm.com/analytics/spss-statistics-software (accessed on 18 December 2021)
31. Li B. Shaikh F. Zamzam A. Syed M.H. Abdin R. Qadura M. The Identification and Evaluation of Interleukin-7 as a Myokine Biomarker for Peripheral Artery Disease Prognosis J. Clin. Med. 2024 13 3583 10.3390/jcm13123583 38930112
32. Ross E.G. Jung K. Dudley J.T. Li L. Leeper N.J. Shah N.H. Predicting Future Cardiovascular Events in Patients with Peripheral Artery Disease Using Electronic Health Record Data Circ. Cardiovasc. Qual. Outcomes 2019 12 e004741 10.1161/CIRCOUTCOMES.118.004741 30857412
33. Kakkar R. Lee R.T. The IL-33/ST2 Pathway: Therapeutic Target and Novel Biomarker Nat. Rev. Drug Discov. 2008 7 827 840 10.1038/nrd2660 18827826
34. Januzzi J.L. ST2 as a Cardiovascular Risk Biomarker: From the Bench to the Bedside J. Cardiovasc. Transl. Res. 2013 6 493 500 10.1007/s12265-013-9459-y 23558647
35. Stojkovic S. Kaun C. Basilio J. Rauscher S. Hell L. Krychtiuk K.A. Bonstingl C. de Martin R. Gröger M. Ay C. Tissue Factor Is Induced by Interleukin-33 in Human Endothelial Cells: A New Link between Coagulation and Inflammation Sci. Rep. 2016 6 25171 10.1038/srep25171 27142573
36. Demyanets S. Konya V. Kastl S.P. Kaun C. Rauscher S. Niessner A. Pentz R. Pfaffenberger S. Rychli K. Lemberger C.E. Interleukin-33 Induces Expression of Adhesion Molecules and Inflammatory Activation in Human Endothelial Cells and in Human Atherosclerotic Plaques Arterioscler. Thromb. Vasc. Biol. 2011 31 2080 2089 10.1161/ATVBAHA.111.231431 21737781
37. Stojkovic S. Thulin Å. Hell L. Thaler B. Rauscher S. Baumgartner J. Gröger M. Ay C. Demyanets S. Neumayer C. IL-33 Stimulates the Release of Procoagulant Microvesicles from Human Monocytes and Differentially Increases Tissue Factor in Human Monocyte Subsets Thromb. Haemost. 2017 117 1379 1390 10.1160/TH16-10-0784 28492698
38. Demyanets S. Tentzeris I. Jarai R. Katsaros K.M. Farhan S. Wonnerth A. Weiss T.W. Wojta J. Speidl W.S. Huber K. An Increase of Interleukin-33 Serum Levels after Coronary Stent Implantation Is Associated with Coronary in-Stent Restenosis Cytokine 2014 67 65 70 10.1016/j.cyto.2014.02.014 24725541
39. Demyanets S. Speidl W.S. Tentzeris I. Jarai R. Katsaros K.M. Farhan S. Krychtiuk K.A. Wonnerth A. Weiss T.W. Huber K. Soluble ST2 and Interleukin-33 Levels in Coronary Artery Disease: Relation to Disease Activity and Adverse Outcome PLoS ONE 2014 9 e95055 10.1371/journal.pone.0095055 24751794
40. Stankovic M. Ljujic B. Babic S. Maravic-Stojkovic V. Mitrovic S. Arsenijevic N. Radak D. Pejnovic N. Lukic M.L. IL-33/IL-33R in Various Types of Carotid Artery Atherosclerotic Lesions Cytokine 2019 120 242 250 10.1016/j.cyto.2019.05.010 31132589
41. Dhillon O.S. Narayan H.K. Quinn P.A. Squire I.B. Davies J.E. Ng L.L. Interleukin 33 and ST2 in Non-ST-Elevation Myocardial Infarction: Comparison with Global Registry of Acute Coronary Events Risk Scoring and NT-proBNP Am. Heart J. 2011 161 1163 1170 10.1016/j.ahj.2011.03.025 21641364
42. Aimo A. Januzzi J.L. Bayes-Genis A. Vergaro G. Sciarrone P. Passino C. Emdin M. Clinical and Prognostic Significance of sST2 in Heart Failure: JACC Review Topic of the Week J. Am. Coll. Cardiol. 2019 74 2193 2203 10.1016/j.jacc.2019.08.1039 31648713
43. Li B. Zamzam A. Syed M.H. Jahanpour N. Jain S. Abdin R. Qadura M. Urinary Fatty Acid Binding Protein 3 (uFABP3) Has Diagnostic and Prognostic Value in Peripheral Artery Disease J. Vasc. Surg. 2022 75 e325 10.1016/j.jvs.2022.03.761
44. Li B. Djahanpour N. Zamzam A. Syed M.H. Jain S. Arfan S. Abdin R. Qadura M. The Prognostic Capability of Inflammatory Proteins in Predicting Peripheral Artery Disease Related Adverse Events Front. Cardiovasc. Med. 2022 9 1073751 10.3389/fcvm.2022.1073751 36582735
45. Li B. Shaikh F. Zamzam A. Syed M.H. Abdin R. Qadura M. A Machine Learning Algorithm for Peripheral Artery Disease Prognosis Using Biomarker Data iScience 2024 27 109081 10.1016/j.isci.2024.109081 38361633
46. Li B. Nassereldine R. Zamzam A. Syed M.H. Mamdani M. Al-Omran M. Abdin R. Qadura M. Development and Evaluation of a Prediction Model for Peripheral Artery Disease-Related Major Adverse Limb Events Using Novel Biomarker Data J. Vasc. Surg. 2024 80 490 497.e1 10.1016/j.jvs.2024.03.450 38599293
47. Li B. Zamzam A. Syed M.H. Jahanpour N. Jain S. Abdin R. Qadura M. Urinary Cystatin C Has Prognostic Value in Peripheral Artery Disease Biomolecules 2022 12 860 10.3390/biom12070860 35883416
48. Stoltzfus J.C. Logistic Regression: A Brief Primer Acad. Emerg. Med. 2011 18 1099 1104 10.1111/j.1553-2712.2011.01185.x 21996075
49. Kia B. Mendes A. Parnami A. George R. Mobley K. Ditto W.L. Nonlinear Dynamics Based Machine Learning: Utilizing Dynamics-Based Flexibility of Nonlinear Circuits to Implement Different Functions PLoS ONE 2020 15 e0228534 10.1371/journal.pone.0228534 32126089
50. Smith B.J. Silva-Costa L.C. Martins-de-Souza D. Human Disease Biomarker Panels through Systems Biology Biophys. Rev. 2021 13 1179 1190 10.1007/s12551-021-00849-y 35059036
51. Couronné R. Probst P. Boulesteix A.-L. Random Forest versus Logistic Regression: A Large-Scale Benchmark Experiment BMC Bioinform. 2018 19 270 10.1186/s12859-018-2264-5 30016950
52. Flores A.M. Demsas F. Leeper N.J. Ross E.G. Leveraging Machine Learning and Artificial Intelligence to Improve Peripheral Artery Disease Detection, Treatment, and Outcomes Circ. Res. 2021 128 1833 1850 10.1161/CIRCRESAHA.121.318224 34110911
53. Nakar S. Vinker S. Neuman S. Kitai E. Yaphe J. Baseline Tests or Screening: What Tests Do Family Physicians Order Routinely on Their Healthy Patients? J. Med. Screen. 2002 9 133 134 10.1136/jms.9.3.133 12370326
54. Akbari C.M. Stone L. Accreditation and Credentialing in the Vascular Laboratory Semin. Vasc. Surg. 2002 15 178 181 10.1016/S0895-7967(02)70029-X 12228866
55. Burns P. Gough S. Bradbury A.W. Management of Peripheral Arterial Disease in Primary Care BMJ 2003 326 584 588 10.1136/bmj.326.7389.584 12637405
56. Koh C.E. Walker S.R. Vascular Surgery Consults: A Significant Workload ANZ J. Surg. 2007 77 352 354 10.1111/j.1445-2197.2007.04058.x 17497974
57. Nishimiya K. Matsumoto Y. Shimokawa H. Recent Advances in Vascular Imaging Arterioscler. Thromb. Vasc. Biol. 2020 40 e313 e321 10.1161/ATVBAHA.120.313609 33054393
58. Eikelboom J.W. Connolly S.J. Bosch J. Dagenais G.R. Hart R.G. Shestakovska O. Diaz R. Alings M. Lonn E.M. Anand S.S. Rivaroxaban with or without Aspirin in Stable Cardiovascular Disease N. Engl. J. Med. 2017 377 1319 1330 10.1056/NEJMoa1709118 28844192
59. Farber A. Menard M.T. Conte M.S. Kaufman J.A. Powell R.J. Choudhry N.K. Hamza T.H. Assmann S.F. Creager M.A. Cziraky M.J. Surgery or Endovascular Therapy for Chronic Limb-Threatening Ischemia N. Engl. J. Med. 2022 387 2305 2316 10.1056/NEJMoa2207899 36342173
60. Caetano A.P. Conde Vasco I. Veloso Gomes F. Costa N.V. Luz J.H. Spaepen E. Formiga A. Coimbra É. Neves J. Bilhim T. Successful Revascularization Has a Significant Impact on Limb Salvage Rate and Wound Healing for Patients with Diabetic Foot Ulcers: Single-Centre Retrospective Analysis with a Multidisciplinary Approach Cardiovasc. Interv. Radiol. 2020 43 1449 1459 10.1007/s00270-020-02604-4 32743744
61. Margolis J. Barron J.J. Grochulski W.D. Health Care Resources and Costs for Treating Peripheral Artery Disease in a Managed Care Population: Results from Analysis of Administrative Claims Data J. Manag. Care Pharm. 2005 11 727 734 10.18553/jmcp.2005.11.9.727 16300416
