Exploration of effective biomarkers for venous thrombosis embolism in Behçet’s disease based on comprehensive bioinformatics analysis

Abstract

Behçet’s disease (BD) is a multifaceted autoimmune disorder affecting multiple organ systems. Vascular complications, such as venous thromboembolism (VTE), are highly prevalent, affecting around 50% of individuals diagnosed with BD. This study aimed to identify potential biomarkers for VTE in BD patients. Three microarray datasets (GSE209567, GSE48000, GSE19151) were retrieved for analysis. Differentially expressed genes (DEGs) associated with VTE in BD were identified using the Limma package and weighted gene co-expression network analysis (WGCNA). Subsequently, potential diagnostic genes were explored through protein–protein interaction (PPI) network analysis and machine learning algorithms. A receiver operating characteristic (ROC) curve and a nomogram were constructed to evaluate the diagnostic performance for VTE in BD patients. Furthermore, immune cell infiltration analyses and single-sample gene set enrichment analysis (ssGSEA) were performed to investigate potential underlying mechanisms. Finally, the efficacy of listed drugs was assessed based on the identified signature genes. The limma package and WGCNA identified 117 DEGs related to VTE in BD. A PPI network analysis then selected 23 candidate hub genes. Four DEGs ( E2F1, GATA3, HDAC5, and MSH2) were identified by intersecting gene sets from three machine learning algorithms. ROC analysis and nomogram construction demonstrated high diagnostic accuracy for these four genes (AUC: 0.816, 95% CI: 0.723–0.909). Immune cell infiltration analysis revealed a positive correlation between dysregulated immune cells and the four hub genes. ssGSEA provided insights into potential mechanisms underlying VTE development and progression in BD patients. Additionally, therapeutic agent screening identified potential drugs targeting the four hub genes. This study employed a systematic approach to identify four potential hub genes ( E2F1, GATA3, HDAC5, and MSH2) and construct a nomogram for VTE diagnosis in BD. Immune cell infiltration analysis revealed dysregulation, suggesting potential macrophage involvement in VTE development. ssGSEA provided insights into potential mechanisms underlying BD-induced VTE, and potential therapeutic agents were identified.