Comprehensive bioinformatics analysis reveals the crosstalk genes and immune relationship between the systemic lupus erythematosus and venous thromboembolism.

Abstract

It is well known that patients with systemic lupus erythematosus (SLE) had a high risk of venous thromboembolism (VTE). This study aimed to identify the crosstalk genes between SLE and VTE and explored their clinical value and molecular mechanism initially. We downloaded microarray datasets of SLE and VTE from the Gene Expression Omnibus (GEO) dataset. Differential expression analysis was applied to identify the crosstalk genes (CGs). Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed on the shared genes. The shared diagnostic biomarkers of the two diseases were further screened from CGs using least absolute shrinkage and selection operator (Lasso) regression. Two risk scores for SLE and VTE were constructed separately to predict the likelihood of illness according to the diagnostic biomarkers using a logical regression algorithm. The immune infiltration levels of SEL and VTE were estimated via the CIBERSORT algorithm and the relationship of CGs with immune cell infiltration was investigated. Finally, we explored potential phenotype subgroups in SLE and VTE based on the expression level of CGs through the consensus clustering method and studied immune cell infiltration in different subtypes. A total of 171 CGs were obtained by the intersection of differentially expressed genes (DEGs) between SLE and VTE cohorts. The functional enrichment shown these CGs were mainly related to immune pathways. After screening by lasso regression, we found that three hub CGs (RSAD2, HSP90AB1, and FPR2) were the optimal shared diagnostic biomarkers for SLE and VTE. Based on the expression level of RSAD2 and HSP90AB1, two risk prediction models for SLE and VTE were built by multifactor logistic regression and systemically validated in internal and external validation datasets. The immune infiltration results revealed that CGs were highly correlated with multiple infiltrated immunocytes. Consensus clustering was used to respectively regroup SLE and VTE patients into C1 and C2 clusters based on the CGs expression profile. The levels of immune cell infiltration and immune activation were higher in C1 than in C2 subtypes. In our study, we further screen out diagnostic biomarkers from crosstalk genes SLE and VTE and built two risk scores. Our findings reveal a close relationship between CGs and the immune microenvironment of diseases. This provides clues for further exploring the common mechanism and interaction between the two diseases.