The shared biomarkers and immune landscape in psoriatic arthritis and rheumatoid arthritis: Findings based on bioinformatics, machine learning and single-cell analysis

by Kaiyi Zhou, Siyu Luo, Qinxiao Wang, Sheng Fang

Psoriatic arthritis (PsA) and rheumatoid arthritis (RA) are the most common types of inflammatory musculoskeletal disorders that share overlapping clinical features and complications. The aim of this study was to identify shared marker genes and mechanistic similarities between PsA and RA.

We utilized datasets from the Gene Expression Omnibus (GEO) database to identify differentially expressed genes (DEGs) and perform functional enrichment analyses. To identify the marker genes, we applied two machine learning algorithms: the least absolute shrinkage and selection operator (LASSO) and the support vector machine recursive feature elimination (SVM-RFE). Subsequently, we assessed the diagnostic capacity of the identified marker genes using the receiver operating characteristic (ROC) curve and decision curve analysis (DCA). A transcription factor (TF) network was constructed using data from JASPAR, HumanTFDB, and GTRD. We then employed CIBERSORT to analyze the abundance of immune infiltrates in PsA and RA, assessing the relationship between marker genes and immune cells. Additionally, cellular subpopulations were identified by analyzing single-cell sequencing data from RA, with T cells examined for trajectory and cellular communication using Monocle and CellChat, thereby exploring their linkage to marker genes.

A total of seven overlapping DEGs were identified between PsA and RA. Gene enrichment analysis revealed that these genes were associated with mitochondrial respiratory chain complex IV, Toll-like receptors, and NF-κB signaling pathways. Both machine learning algorithms identified Ribosomal Protein L22-like 1 (RPL22L1) and Lymphocyte Antigen 96 (LY96) as potential diagnostic markers for PsA and RA. These markers were validated using test sets and experimental approaches. Furthermore, GSEA analysis indicated that gap junctions may play a crucial role in the pathogenesis of both conditions. The TF network suggested a potential association between marker genes and core enrichment genes related to gap junctions. The application of CIBERSORT and single-cell RNA sequencing provided a comprehensive understanding of the role of marker genes in immune cell function. Our results indicated that RPL22L1 and LY96 are involved in T cell development and are associated with T cell communication with NK cells and monocytes. Notably, high expression of both RPL22L1 and LY96 was linked to enhanced VEGF signaling in T cells.

Our study identified RPL22L1 and LY96 as key biomarkers for PsA and RA. Further investigations demonstrated that these two marker genes are closely associated with gap junction function, T cell infiltration, differentiation, and VEGF signaling. Collectively, these findings provide new insights into the diagnosis and treatment of PsA and RA.