Dataset Information

Identifying Immune Cell Infiltration and Effective Diagnostic Biomarkers in Rheumatoid Arthritis by Bioinformatics Analysis.

ABSTRACT:

Background

Rheumatoid arthritis (RA) is a chronic systemic autoimmune disorder characterized by inflammatory cell infiltration, leading to persistent synovitis and joint destruction. The pathogenesis of RA remains unclear. This study aims to explore the immune molecular mechanism of RA through bioinformatics analysis.

Methods

Five microarray datasets and a high throughput sequencing dataset were downloaded. CIBERSORT algorithm was performed to evaluate immune cell infiltration in synovial tissues between RA and healthy control (HC). Wilcoxon test and Least Absolute Shrinkage and Selection Operator (LASSO) regression were conducted to identify the significantly different infiltrates of immune cells. Differentially expressed genes (DEGs) were screened by "Batch correction" and "RobustRankAggreg" methods. Functional correlation of DEGs were analyzed by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG). Candidate biomarkers were identified by cytoHubba of Cytoscape, and their diagnostic effectiveness was predicted by Receiver Operator Characteristic Curve (ROC) analysis. The association of the identified biomarkers with infiltrating immune cells was explored using Spearman's rank correlation analysis in R software.

Results

Ten significantly different types of immune cells between RA and HC were identified. A total of 202 DEGs were obtained by intersection of DEGs screened by two methods. The function of DEGs were significantly associated with immune cells. Five hub genes (CXCR4, CCL5, CD8A, CD247, and GZMA) were screened by R package "UpSet". CCL5+CXCR4 and GZMA+CD8A were verified to have the capability to diagnose RA and early RA with the most excellent specificity and sensitivity, respectively. The correlation between immune cells and biomarkers showed that CCL5 was positively correlated with M1 macrophages, CXCR4 was positively correlated with memory activated CD4+ T cells and follicular helper T (Tfh) cells, and GZMA was positively correlated with Tfh cells.

Conclusions

CCL5, CXCR4, GZMA, and CD8A can be used as diagnostic biomarker for RA. GZMA-Tfh cells, CCL5-M1 macrophages, and CXCR4- memory activated CD4+ T cells/Tfh cells may participate in the occurrence and development of RA, especially GZMA-Tfh cells for the early pathogenesis of RA.

SUBMITTER: Zhou S

PROVIDER: S-EPMC8411707 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:Heart failure (HF) and osteoarthritis (OA) are medical conditions that can significantly impact daily activities. Evidence has shown that HF and OA may share some pathogenic mechanisms. However, the underlying genomic mechanisms remain unclear. This study aimed to explore the underlying molecular mechanism and identify diagnostic biomarkers for HF and OA. With the cutoff criteria of fold change (FC) > 1.3 and P < .05, 920, 1500, 2195, and 2164 differentially expressed genes (DEGs) were identified in GSE57338, GSE116250, GSE114007, and GSE169077, respectively. After making the intersection of DEGs, we obtained 90 upregulated DEGs and 51 downregulated DEGs in HF datasets and 115 upregulated DEGs and 75 downregulated DEGs in OA datasets. Afterward, we conducted genome ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses, protein-protein interaction (PPI) networks, and hub genes screening based on DEGs. Then, 4 common DEGs (fibroblast activation protein alpha [FAP], secreted frizzled-related protein 4 (SFRP4), Thy-1 cell surface antigen (THY1), matrix remodeling associated 5 [MXRA5]) between HF and OA were screened and validated in GSE5406 and GSE113825 datasets, based on which we established the support vector machine (SVM) models. The combined area under the receiver operating characteristic curve (AUC) of THY1, FAP, SFRP4, and MXRA5 in the HF training and test sets reached 0.949 and 0.928. While in the OA training set and test set, the combined AUC of THY1, FAP, SFRP4, and MXRA5 reached 1 and 1, respectively. The analysis of immune cells in HF revealed high levels of dendritic cell (DC), B cells, natural killer T cell (NKT), Type 1 regulatory T cell (Tr1), cytotoxic T cell (Tc), exhausted T cell (Tex), and mucosal-associated invariant T cell (MAIT), while displaying lower levels of monocytes, macrophages, NK, CD4 + T, gamma delta T (γδ T), T helper type 1 (Th1), T helper type 2 (Th2), and effector memory T cell (Tem). Moreover, the 4 common DEGs were positively correlated with DCs and B cells and negatively correlated with γδ T. In OA patients, the abundance of monocyte, macrophage, CD4 + naïve, and natural T regulatory cell (nTreg) was higher, while the infiltration of CD8 + T, γδ T, CD8 + naïve, and MAIT was lower. The expression of THY1 and FAP was significantly correlated with macrophage, CD8 + T, nTreg, and CD8 + naïve. SFRP4 was correlated with monocyte, CD8 + T, γδ T, CD4 + naïve, nTreg, CD8 + naïve and MAIT. MXRA5 was correlated with macrophage, CD8 + T, nTreg and CD8 + naïve. FAP, THY1, MXRA5, and SFRP4 may be diagnostic biomarkers for both HF and OA, and their correlation with immune cell infiltrations suggests shared immune pathogenesis.

Project description:Background: Atopic dermatitis (AD) is a dermatological disorder characterized by symptoms such as chronically inflamed skin and frequently intolerable itching. The mechanism underlying AD development is still unclear. Our study aims to identify the diagnostic and therapeutic biomarkers for AD and provide insight into immune mechanisms at the molecular level through bioinformatics analysis. Methods: The GSE6012, GSE32924, and GSE36842 gene expression profiles were obtained for analysis from the Gene Expression Omnibus database. Differentially expressed genes (DEGs) were segregated using the “Batch correction” and “RobustRankAggreg” methods. Weighted gene co-expression network analysis (WGCNA) was performed to screen for module genes with AD traits. Then, common DEGs (co-DEGs) were screened out via combined differential expression analysis and WGCNA. Functional enrichment analysis was performed for these co-DEGs using Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG), followed by protein-protein interaction network analysis. Candidate hub genes were identified using the “cytoHubba” plugin in Cytoscape, and their value for AD diagnosis was validated using receiver operating characteristic curve analysis in the external database GSE120721. Immunohistochemical staining was performed for further validation. The CIBERSORT algorithm was used to evaluate skin samples obtained from healthy controls (HCs) and lesions of AD patients, to determine the extent of immune cell infiltration. The association between the identified hub genes and significant differential immune cells was analyzed using Pearson correlation analysis. Results: A total of 259 DEGs were acquired from the intersection of DEGs obtained by the two independent procedures, and 331 AD-trait module genes were separated out from the blue module via WGCNA analysis. Then, 169 co-DEGs arising from the intersection of the 259 DEGs and the 331 AD-trait module genes were obtained. We found that co-DEGs were significantly enhanced in the type I interferon and IL-17 signal transduction pathways. Thirteen potential hub genes were identified using Cytoscape. Five hub genes (CCR7, CXCL10, IRF7, MMP1, and RRM2) were identified after screening via external dataset validation and immunohistochemical analysis. We also identified four significant differential immune cells, i.e., activated dendritic cells, plasma cells, resting mast cells, and CD4+ naïve T cells, between AD patients and HCs. Moreover, the relationship between the identified hub genes and significant differential immune cells was analyzed. The results showed that the CCR7 expression level was positively correlated with the number of CD4+ naïve T cells (R = 0.42, p = 0.011). Conclusion: CCR7, CXCL10, IRF7, MMP1, and RRM2 could be potential diagnostic and therapeutic biomarkers for AD. CCR7 expression level was positively correlated with the number of CD4+ naïve T cells in AD. These findings need to be corroborated in future studies.

Project description:BackgroundOsteoarthritis (OA) and rheumatoid arthritis (RA) are well-known cause of joint disability. Although they have shown the analogous clinical features involving chronic synovitis that progresses to cartilage and bone destruction, the pathogenesis that initiates and perpetuates synovial lesions between RA and OA remains elusive.ObjectiveThis study is aimed at identifying disease-specific hub genes, exploring immune cell infiltration, and elucidating the underlying mechanisms associated with RA and OA synovial lesion.MethodsGene expression profiles (GSE55235, GSE55457, GSE55584, and GSE12021) were selected from Gene Expression Omnibus for analysis. Differentially expressed genes (DEGs) were identified by the "LIMMA" package in Bioconductor. The DEGs were identified by Gene Ontology (GO) and KEGG pathway analysis. A protein-protein interaction network was constructed to identify candidate hub genes by using STRING and Cytoscape. Hub genes were identified by validating from GSE12021. Furthermore, we employed the CIBERSORT website to assess immune cell infiltration between OA and RA. Finally, we explored the correlation between the levels of hub genes and relative proportion of immune cells in OA and RA.ResultsWe identified 68 DEGs which were mainly enriched in immune response and chemokine signaling pathway. Six hub genes with a cutoff of AUC > 0.80 by ROC analysis and relative expression of P < 0.05 were identified successfully. Compared with OA, the RA synovial tissues consisted of a higher proportion of 7 immune cells, whereas 4 immune cells were found in relatively lower proportion (P < 0.05). In addition, the levels of 6 hub genes were closely associated with relative proportion of 11 immune cells in OA and RA.ConclusionsWe used bioinformatics analysis to identify hub genes and explored immune cell infiltration of immune microenvironment in synovial tissues. Our results should offer insights into the underlying molecular mechanisms of synovial lesion and provide potential target for immune-based therapies of OA and RA.