Integration of bulk RNA-seq and scRNA-seq reveals transcriptomic signatures associated with deep vein thrombosis

Bao-Ze PANMing-Jun JiangLiming DengJie ChenXian-Peng DaiZi-Xuan WuZhihe DengDong-Yang LuoYang-Yi-Jing WangDan NingGuo-Zuo XiongGuoshan Bi^*

• Second Affiliated Hospital of University of South China, Hengyang, China

Background: Deep vein thrombosis (DVT) is a prevalent peripheral vascular disease. The intricate and multifaceted nature of the associated mechanisms hinders a comprehensive understanding of disease-relevant targets. This study aimed to identify and examine the most distinctive genes linked to DVT.In this study, the bulk RNA sequencing (bulk RNA-seq) analysis was conducted on whole blood samples from 11 DVT patients and 6 control groups.Topology analysis was performed using 7 protein-protein interaction (PPI) network algorithms. The combination of weighted correlation network analysis (WGCNA) and clinical prediction models was employed to validate hub DEGs. Furthermore, single-cell RNA sequencing (scRNA-seq) was performed on peripheral blood samples from 3 DVT patients and 3 control groups to probe the cellular localization of target genes. Based on the same methodology as the internal test set, 12 DVT patients and 6 control groups were collected to construct an external test set and validated using machine learning (ML) algorithms and immunofluorescence (IF).Concurrently, the examination of the pathways in disparate cell populations was conducted on the basis of the CellChat pathway.Results: A total of 193 DEGs were identified in the internal test set. Additionally, a total of 8 highly characteristic genes (including TLR1, TLR7, TLR8, CXCR4, DDX58, TNFSF10, FCGR1A and CD36) were identified by the PPI network algorithm. In accordance with the WGCNA model, the aforementioned genes were all situated within the blue core module, exhibiting a correlation coefficient of 0.84. The model demonstrated notable disparities in TLR8 (P = 0.018, AUC = 0.847), CXCR4 (P = 0.00088, AUC = 1.000), TNFSF10 (P = 0.00075, AUC = 0.958), and FCGR1A (P = 0.00022, AUC = 0.986). Furthermore, scRNA-seq demonstrated that B cells, T cells and monocytes play an active role in DVT. In the external validation set, CXCR4 was validated as a potential target by the ML algorithm and IF. In the context of the CellChat pathway, it indicated that MIF -(CD74 + CXCR4) plays a potential role.The findings of this study indicate that CXCR4 may serve as a potential genetic marker for DVT, with MIF -(CD74 + CXCR4) potentially implicated in the regulatory mechanisms underlying DVT.