Abstract
The pathogenesis of Crohn's disease (CD) involves abnormal immune cell infiltration and dysregulated immune response. Therefore, thorough research on immune cell abnormalities in CD is crucial for improved treatment of this disease. Single-cell RNA sequencing (scRNA-seq) and bulk RNA-seq data of CD were obtained from the Gene Expression Omnibus (GEO) database. Cell-type identification by estimating relative subsets of RNA transcripts (CIBERSORT), weighted gene co-expression network analysis (WGCNA), protein-protein interaction (PPI) networks evaluated the proportion of immune infiltrating cells, constructed co-expression network and identified key genes, respectively. Based on the dataset (GSE134809), 15 cell clusters were defined and labeled as different cell types. Among the 11 modules, the yellow module had the closest relationship with plasma cells (cluster 5). Confirmed using RNA sequencing and IHC assay, the expression of COL5A2 in CD samples was higher than that in control samples. Furthermore, the COL5A2 protein expression remarkably decreased in the group of patients who responded to anti-tumor necrosis factor (TNF) treatments, compared to the non-response group. The comprehensive analyses described here provided novel insight into the landscape of CD-associated immune environment. In addition, COL5A2 were identified as potential diagnostic indicators for CD, as well as promising predictive markers for CD patients.