Identifying the hub genes and immune infiltration related to pyroptosis in rheumatoid arthritis.
Ontology highlight
ABSTRACT: Abstract
Rheumatoid arthritis (RA) is one of the most common autoimmune joint disorders globally, but its pathophysiological mechanisms have not been thoroughly investigated. Pyroptosis significantly correlates with programmed cell death. However, targeting pyroptosis has not been considered as a therapeutic strategy in RA due to a lack of systematic studies on validated biomarkers. The present study aimed to identify hub pyroptosis biomarkers and immune infiltration in RA. The gene expression profiles of synovial tissues were obtained from the Gene Expression Omnibus (GEO) database to identify differentially expressed pyroptosis genes (DEPGs). Meanwhile, the CIBERSORT algorithm was used to explore the association between immune infiltration and RA. Consequently, two hub DEPGs (EGFR and JUN) were identified as critical genes in RA. Through gene ontology and pathway enrichment analysis. EGFR and JUN were found to be primarily involved in the ErbB signaling pathway, PD-1 checkpoint pathway, GnRH signaling pathway, etc. Furthermore, for immune infiltration analysis, the pyroptosis genes EGFR and JUN were closely connected with four and one immune cell types, respectively. Overall, this study presents a novel method to identify hub DEPGs and their correlation with immune infiltration, which may provide novel perspectives into the diagnosis and treatment of patients with RA.
SUBMITTER: Xie W
PROVIDER: S-EPMC8677948 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:BackgroundRheumatoid 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.MethodsFive 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.ResultsTen 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.ConclusionsCCL5, 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.
| S-EPMC8411707 | biostudies-literature
Project description:IntroductionPeriodontitis is an inflammatory disease and its molecular mechanisms is not clear. A recently discovered cell death pathway called cuproptosis, may related to the disease.MethodsThe datasets GSE10334 of human periodontitis and control were retrieved from the Gene Expression Omnibus database (GEO) for analysis.Following the use of two machine learning algorithms, least absolute shrinkage and selection operator (LASSO) and support vector machine-recursive feature removal (SVM-RFE) were used to find CRG-based signature. Then the Receiver operating characteristic (ROC) curves was used to evaluate the gene signature's discriminatory ability. The CIBERSORT deconvolution algorithm was used to study the link between hub genes and distinct types of immune cells. Next, the association of the CRGs with immune cells in periodontitis and relevant clusters of cuproptosis were found. The link between various clusters was ascertained by the GSVA and CIBERSORT deconvolution algorithm. Finally, An external dataset (GSE16134) was used to confirm the diagnosis capacity of the identified biomarkers. In addition, clinical samples were examined using qRT-PCR and immunohistochemistry to verifiy the expression of genes related to cuprotosis in periodontitis and the signature may better predict the periodontitis.Results15 periodontitis-related DE-CRGs were found,then 11-CRG-based signature was found by using of LASSO and SVM-RFE. ROC curves also supported the value of signature. CIBERSORT results of immune cell signature in periodontitis showed that signature genes is a crucial component of the immune response.The relevant clusters of cuproptosis found that the NFE2L2, SLC31A1, FDX1,LIAS, DLD, DLAT, and DBT showed a highest expression levels in Cluster2 ,while the NLRP3, MTF1, and DLST displayed the lowest level in Cluster 2 but the highest level in Cluster1. The GSVA results also showed that the 11 cuproptosis diagnostic gene may regulate the periodontitis by affecting immune cells. The external dataset (GSE16134) confirm the diagnosis capacity of the identified biomarkers, and clinical samples examined by qRT-PCR and immunohistochemistry also verified that these cuprotosis related signiture genes in periodontitis may better predict the periodontitis.ConclusionThese findings have important implications for the cuproptosis and periodontitis, and highlight further research is needed to better understand the mechanisms underlying this relationship between the cuproptosis and periodontitis.
| S-EPMC10196202 | biostudies-literature
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.
| S-EPMC8152926 | biostudies-literature
Project description:Breast cancer (BC) is the second leading cause of death among women and is highly heterogeneous. Three pyroptosis (PR) subtypes were identified in patients with BC from the Cancer Genome Atlas Database (TCGA) cohorts using 20 PR-related regulators, which illustrate a strong association between BC prognosis and PR. Lung metastasis commonly occurs in the advanced stages of BC, resulting in a poor quality of life. Eight differentially expressed (DE) lncRNAs were identified using LASSO-Cox analysis between PR-related and BC lung metastasis. Moreover, a BRCA risk-score (RS) model was established using multivariate Cox regression, which correlated with prognosis in TCGA-BRCA. Clinical characteristics, tumor mutation burden, and tumor immune cell infiltration were extensively investigated between high- and low-RS groups. Similarly, a lower RS implied longer overall survival, greater inflammatory cell infiltration, and better immunotherapeutic response to PD-1 blockers. Our findings provide a foundation for future studies targeting PR and confirme that RS could predict the prognosis of patients with BC.
| S-EPMC8931340 | biostudies-literature
Project description:BackgroundThe current study is aimed at identifying the cross-talk genes between periodontitis (PD) and rheumatoid arthritis (RA), as well as the potential relationship between cross-talk genes and pyroptosis-related genes.MethodsDatasets for the PD (GSE106090, GSE10334, GSE16134) and RA (GSE55235, GSE55457, GSE77298, and GSE1919) were downloaded from the GEO database. After batch correction and normalization of datasets, differential expression analysis was performed to identify the differentially expressed genes (DEGs). The cross-talk genes linking PD and RA were obtained by overlapping the DEGs dysregulated in PD and DEGs dysregulated in RA. Genes involved in pyroptosis were summarized by reviewing literatures, and the correlation between pyroptosis genes and cross-talk genes was investigated by Pearson correlation coefficient. Furthermore, the weighted gene coexpression network analysis (WGCNA) was carried out to identify the significant modules which contained both cross-talk genes and pyroptosis genes in both PD data and RA data. Thus, the core cross-talk genes were identified from the significant modules. Receiver-operating characteristic (ROC) curve analysis was performed to identify the predictive accuracy of these core cross-talk genes in diagnosing PD and RA. Based on the core cross-talk genes, the experimentally validated protein-protein interaction (PPI) and gene-pathway network were constructed.ResultsA total of 40 cross-talk genes were obtained. Most of the pyroptosis genes were not differentially expressed in disease and normal samples. By selecting the modules containing both cross-talk genes or pyroptosis genes, the blue module was identified to be significant module. Three genes, i.e., cross-talk genes (TIMP1, LGALS1) and pyroptosis gene-GPX4, existed in the blue module of PD network, while two genes (i.e., cross-talk gene-VOPP1 and pyroptosis gene-AIM2) existed in the blue module of RA network. ROC curve analysis showed that three genes (TIMP1, VOPP1, and AIM2) had better predictive accuracy in diagnosing disease compared with the other two genes (LGALS1 and GPX4).ConclusionsThis study revealed shared mechanisms between RA and PD based on cross-talk and pyroptosis genes, supporting the relationship between the two diseases. Thereby, five modular genes (TIMP1, LGALS1, GPX4, VOPP1, and AIM2) could be of relevance and might serve as potential biomarkers. These findings are a basis for future research in the field.
| S-EPMC8731299 | biostudies-literature
Project description:Melanoma is an aggressive skin cancer that has gained attention worldwide. Growing evidence has highlighted that the tumor microenvironment (TME) is an important feature of carcinogenesis and contributes to therapeutic efficacy in melanoma. However, additional advances in melanoma immuno-oncology are necessary to achieve a comprehensive knowledge of the immune infiltrate population and to identify accurate and readily measurable biomarkers. In this study, we analyzed gene expression of 468 melanoma cases from the TCGA database, which led to the identification of three melanoma clusters (representedby low, median and high infiltration) that display unique immune features. We found that the microenvironment clusters had substantial prognostic efficacy. The median cluster was characterized by an inability to draw immune cells, highlighting possible immune escape mechanisms, and lower CXCL9 and CXCL10 expression, which was correlated to poor prognosis. Deep molecular characterization of immune cells, cytolytic-activity and tumor-inflammatory status revealed diversity of the local immune infiltration landscape in the melanoma clusters. Differentially expressed genes related to TME were extracted from each infiltration cluster. Functional annotations revealed that these genes were mainly related to immune system activation and the processes of immunoreaction. The top ten hub genes in immune infiltration-related protein-protein interaction (PPI) networks were selected for further prognostic investigation. Further validation showed that five of ten hub genes were good prognostic biomarkers for melanoma in two independent groups from the Gene Expression Omnibus database. In brief, these data highlight that systemic characterization of melanoma could uncover tumor infiltrate characteristics, which can help select the most adequate treatment and identify consistent and important indicators of the local immune tumor microenvironment in melanoma patients.
| S-EPMC7778557 | biostudies-literature
Project description:Although several genes (including a strong effect in the human leukocyte antigen (HLA) region) and some environmental factors have been implicated to cause susceptibility to rheumatoid arthritis (RA), the etiology of the disease is not completely understood. The ability to screen the entire genome for association to complex diseases has great potential for identifying gene effects. However, the efficiency of gene detection in this situation may be improved by methods specifically designed for high-dimensional data. The aim of this study was to compare how three different statistical approaches, multifactor dimensionality reduction (MDR), random forests (RF), and an omnibus approach, worked in identifying gene effects (including gene-gene interaction) associated with RA. We developed a test set of genes based on previous linkage and association findings and tested all three methods. In the presence of the HLA shared-epitope factor, other genes showed weaker effects. All three methods detected SNPs in PTPN22 and TRAF1-C5 as being important. But we did not detect any new genes in this study. We conclude that the three high-dimensional methods are useful as an initial screening for gene associations to identify promising genes for further modeling and additional replication studies.
| S-EPMC2795981 | biostudies-literature
Project description:BackgroundThe disability rate associated with rheumatoid arthritis (RA) ranks high among inflammatory joint diseases. However, the cause and potential molecular events are as yet not clear. Here, we aimed to identify differentially expressed genes (DEGs), pathways and immune infiltration involved in RA utilizing integrated bioinformatics analysis and investigating potential molecular mechanisms.Materials and methodsThe expression profiles of GSE55235, GSE55457, GSE55584 and GSE77298 were downloaded from the Gene Expression Omnibus database, which contained 76 synovial membrane samples, including 49 RA samples and 27 normal controls. The microarray datasets were consolidated and DEGs were acquired and further analyzed by bioinformatics techniques. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses of DEGs were performed using R (version 3.6.1) software, respectively. The protein-protein interaction (PPI) network of DEGs were developed utilizing the STRING database. Finally, the CIBERSORT was used to evaluate the infiltration of immune cells in RA.ResultsA total of 828 DEGs were recognized, with 758 up-regulated and 70 down-regulated. GO and KEGG pathway analyses demonstrated that these DEGs focused primarily on cytokine receptor activity and relevant signaling pathways. The 30 most firmly related genes among DEGs were identified from the PPI network. The principal component analysis showed that there was a significant difference between the two tissues in infiltration immune.ConclusionThis study shows that screening for DEGs, pathways and immune infiltration utilizing integrated bioinformatics analyses could aid in the comprehension of the molecular mechanisms involved in RA development. Besides, our study provides valuable data related to DEGs, pathways and immune infiltration of RA and may provide new insight into the understanding of molecular mechanisms.
| S-EPMC7784358 | biostudies-literature
Project description:Rheumatoid arthritis (RA) is an autoimmune disease characterized by erosive arthritis, which has not been thoroughly cured yet, and standardized treatment is helpful for alleviating clinical symptoms. Here, various bioinformatics analysis tools were comprehensively utilized, aiming to identify critical biomarkers and possible pathogenesis of RA. Three gene expression datasets profiled by microarray were obtained from GEO database. Dataset GSE55235 and GSE55457 were merged for subsequent analyses. We identified differentially expressed genes (DEGs) in RStudio with limma package, performing functional enrichment analysis based on GSEA software and clusterProfiler package. Next, protein-protein interaction (PPI) network was set up through STRING database and Cytoscape. Moreover, CIBERSORT website was used to assess the inflammatory state of RA. Finally, we validated the candidate hub genes with dataset GSE77298. As a result, we identified 106 DEGs (72 upregulated and 34 downregulated genes). Through GO, KEGG, and GSEA analysis, we found that DEGs were mainly involved in immune response and inflammatory signaling pathway. With the help of Cytoscape software and MCODE plug-in, the most prominent subnetwork was screened out, containing 14 genes and 45 edges. For ROC curve analysis, eight genes with AUC >0.80 were considered as hub genes of RA. In conclusion, compared with healthy controls, the DEGs and their closely related biological functions were analyzed, and we held that chemokines and immune cells infiltration promote the progression of rheumatoid arthritis. Targeting the eight biomarkers we identified may be useful for the diagnosis and treatment of rheumatoid arthritis.
| S-EPMC7424395 | biostudies-literature
Project description:BackgroundRheumatoid arthritis (RA), a chronic autoimmune inflammatory disease, is often characterized by persistent morning stiffness, joint pain, and swelling. Early diagnosis and timely treatment of RA can effectively delay the progression of the condition and significantly reduce the incidence of disability. In the study, we explored the function of pyroptosis-related genes (PRGs) in the diagnosis and classification of rheumatoid arthritis based on Gene Expression Omnibus (GEO) datasets.MethodWe downloaded the GSE93272 dataset from the GEO database, which contains 35 healthy controls and 67 RA patients. Firstly, the GSE93272 was normalized by the R software "limma" package. Then, we screened PRGs by SVM-RFE, LASSO, and RF algorithms. To further investigate the prevalence of RA, we established a nomogram model. Besides, we grouped gene expression profiles into two clusters and explored their relationship with infiltrating immune cells. Finally, we analyzed the relationship between the two clusters and the cytokines.ResultCHMP3, TP53, AIM2, NLRP1, and PLCG1 were identified as PRGs. The nomogram model revealed that decision-making based on established model might be beneficial for RA patients, and the predictive power of the nomogram model was significant. In addition, we identified two different pyroptosis patterns (pyroptosis clusters A and B) based on the 5 PRGs. We found that eosinophil, gamma delta T cell, macrophage, natural killer cell, regulatory T cell, type 17 T helper cell, and type 2 T helper cell were significant high expressed in cluster B. And, we identified gene clusters A and B based on 56 differentially expressed genes (DEGs) between pyroptosis cluster A and B. And we calculated the pyroptosis score for each sample to quantify the different patterns. The patients in pyroptosis cluster B or gene cluster B had higher pyroptosis scores than those in pyroptosis cluster A or gene cluster A.ConclusionIn summary, PRGs play vital roles in the development and occurrence of RA. Our findings might provide novel views for the immunotherapy strategies with RA.
| S-EPMC10057543 | biostudies-literature