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:BackgroundOsteoarthritis (OA) is the most common chronic degenerative joint disorder globally that is characterized by synovitis, cartilage degeneration, joint space stenosis, and sub-cartilage bone hyperplasia. However, the pathophysiologic mechanisms of OA have not been thoroughly investigated.MethodsIn this study, we conducted various bioinformatics analyses to identify hub biomarkers and immune infiltration in OA. The gene expression profiles of synovial tissues from 29 healthy controls and 36 OA samples were obtained from the gene expression omnibus database to identify differentially expressed genes (DEGs). The CIBERSORT algorithm was used to explore the association between immune infiltration and arthritis.ResultsEighteen hub DEGs were identified as critical biomarkers for OA. Through gene ontology and pathway enrichment analyses, it was found that these DEGs were primarily involved in PI3K-Akt signaling pathway and Rap1 signaling pathway. Furthermore, immune infiltration analysis revealed differences in immune infiltration between patients with OA and healthy controls. The hub gene ZNF160 was closely related to immune cells, especially mast cell activation in OA.ConclusionOverall, this study presented a novel method to identify hub DEGs and their correlation with immune infiltration, which may provide novel insights into the diagnosis and treatment of patients with OA.
| S-EPMC8527722 | biostudies-literature
Project description:BackgroundRheumatoid arthritis (RA) is an autoimmune disease that affects individuals of all ages. The basic pathological manifestations are synovial inflammation, pannus formation, and erosion of articular cartilage, bone destruction will eventually lead to joint deformities and loss of function. However, the specific molecular mechanisms of synovitis tissue in RA are still unclear. Therefore, this study aimed to screen and explore the potential hub genes and immune cell infiltration in RA.MethodsThree microarray datasets (GSE12021, GSE55457, and GSE55235), from the Gene Expression Omnibus (GEO) database, have been analyzed to explore the potential hub genes and immune cell infiltration in RA. First, the LIMMA package was used to screen the differentially expression genes (DEGs) after removing the batch effect. Then the clusterProfiler package was used to perform functional enrichment analyses. Second, through weighted coexpression network analysis (WGCNA), the key module was identified in the coexpression network of the gene set. Third, the protein-protein interaction (PPI) network was constructed through STRING website and the module analysis was performed using Cytoscape software. Fourth, the CIBERSORT and ssGSEA algorithm were used to analyze the immune status of RA and healthy synovial tissue, and the associations between immune cell infiltration and RA-related diagnostic biomarkers were evaluated. Fifth, we used the quantitative reverse transcription-polymerase chain reaction (qRT-PCR) to validate the expression levels of the hub genes, and ROC curve analysis of hub genes for discriminating between RA and healthy tissue. Finally, the gene-drug interaction network was constructed using DrugCentral database, and identification of drug molecules based on hub genes using the Drug Signature Database (DSigDB) by Enrichr.ResultsA total of 679 DEGs were identified, containing 270 downregulated genes and 409 upregulated genes. DEGs were primarily enriched in immune response and chemokine signaling pathways, according to functional enrichment analysis of DEGs. WGCNA explored the co-expression network of the gene set and identified key modules, the blue module was selected as the key module associated with RA. Seven hub genes are identified when PPI network and WGCNA core modules are intersected. Immune infiltration analysis using CIBERSORT and ssGSEA algorithms revealed that multiple types of immune infiltration were found to be upregulated in RA tissue compared to normal tissue. Furthermore, the levels of 7 hub genes were closely related to the relative proportions of multiple immune cells in RA. The results of the qRT-PCR demonstrated that the relative expression levels of 6 hub genes (CD27, LCK, CD2, GZMB, IL7R, and IL2RG) were up-regulated in RA synovial tissue, compared with normal tissue. Simultaneously, ROC curves indicated that the above 6 hub genes had strong biomarker potential for RA (AUC >0.8).ConclusionsThrough bioinformatics analysis and qRT-PCR experiment, our study ultimately discovered 6 hub genes (CD27, LCK, CD2, GZMB, IL7R, and IL2RG) that closely related to RA. These findings may provide valuable direction for future RA clinical diagnosis, treatment, and associated research.
| S-EPMC9868484 | biostudies-literature
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:Rheumatoid arthritis (RA) is a common autoimmune disease that can lead to severe joint damage and disability. And early diagnosis and treatment of RA can avert or substantially slow the progression of joint damage in up to 90% of patients, thereby preventing irreversible disability. Previous research indicated that 50% of the risk for the development of RA is attributable to genetic factors, but the pathogenesis is not well understood. Thus, it is urgent to identify biomarkers to arrest RA before joints are irreversibly damaged. Here, we first use the Robust Rank Aggregation method (RRA) to identify the differentially expressed genes (DEGs) between RA and normal samples by integrating four public RA patients' mRNA expression data. Subsequently, these DEGs were used as the input for the weighted gene co-expression network analysis (WGCNA) approach to identify RA-related modules. The function enrichment analysis suggested that the RA-related modules were significantly enriched in immune-related actions. Then the hub genes were defined as the candidate genes. Our analysis showed that the expression levels of candidate genes were significantly associated with the RA immune microenvironment. And the results indicated that the expression of the candidate genes can use as predictors for RA. We hope that our method can provide a more convenient approach for the early diagnosis of RA.
| S-EPMC9899220 | 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:BackgroundRheumatoid arthritis (RA) is recognized as a chronic inflammatory disease featured by pathological synovial inflammation. Currently, the underlying pathophysiological mechanisms of RA remain unclear. In the study, we attempted to explore the underlying mechanisms of RA and provide potential targets for the therapy of RA via bioinformatics analysis.MethodsWe downloaded four microarray datasets (GSE77298, GSE55235, GSE12021, and GSE55457) from the GEO database. Firstly, GSE77298 and GSE55457 were identified DEGs by the "limma" and "sva" packages of R software. Then, we performed GO, KEGG, and GSEA enrichment analyses to further analyze the function of DEGs. Hub genes were screened using LASSO analysis and SVM-RFE analysis. To further explore the differences of the expression of hub genes in healthy control and RA patient synovial tissues, we calculated the ROC curves and AUC. The expression levels of hub genes were verified in synovial tissues of normal and RA rats by qRT-PCR and western blot. Furthermore, the CIBERSORTx was implemented to assess the differences of infiltration in 22 immune cells between normal and RA synovial tissues. We explored the association between hub genes and infiltrating immune cells.ResultsCRTAM, CXCL13, and LRRC15 were identified as RA's potential hub genes by machine learning and LASSO algorithms. In addition, we verified the expression levels of three hub genes in the synovial tissue of normal and RA rats by PCR and western blot. Moreover, immune cell infiltration analysis showed that plasma cells, T follicular helper cells, M0 macrophages, M1 macrophages, and gamma delta T cells may be engaged in the development and progression of RA.ConclusionsIn brief, our study identified and validated that three hub genes CRTAM, CXCL13, and LRRC15 might involve in the pathological development of RA, which could provide novel perspectives for the diagnosis and treatment with RA.
| S-EPMC9617710 | 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:The effect of pyroptosis-related genes (PRGs) on the tumor microenvironment (TME) in lung adenocarcinoma (LUAD) remains unclear. Thus, this study is aimed at evaluating the prognostic value of PRGs in patients with LUAD and to elucidate their role in the TME and their effect on immunotherapy. Transcriptomic and clinical data were obtained from the Cancer Genome Atlas and the Gene Expression Omnibus databases (GSE3141, GSE31210). Patients with LUAD were classified using consistent clustering, and the differences in the TME for each type were determined using the ESTIMATE and CIBERSORT algorithms. PRGs were screened using univariate regression analysis, and a prognostic risk model was constructed using LASSO regression analysis. The tumor mutational burden and the tumor immune dysfunction and exclusion algorithms were used to predict therapeutic sensitivity in LUAD patients. Then, we evaluated the potential therapeutic interventions using the GDSC database. LUAD patients in cluster 2 had significantly shorter overall survival and progression-free survival rates, lower immune scores, and higher infiltration of T follicular helper cells than those in cluster 1. We used five PRGs to classify patients with LUAD into different risks groups and found that the high-risk group is sensitive to immunotherapy; however, its immune-related pathways were inhibited, which may be related to tumor metabolic reprogramming. Lastly, we identified several potential therapeutic drugs for application in low-risk patients who were less sensitive to immunotherapy. Overall, our findings showed that PRGs can be used to predict prognosis and may aid in the development of personalized therapeutic strategies in LUAD patients.
| S-EPMC9581708 | biostudies-literature
Project description:Little is known on the relationship between the expression of pyroptosis related genes (PRGs) and prognosis of hepatocellular carcinoma (HCC). In this study, a specific PRGs prognostic model was developed with an aim to improve therapeutic efficiency among HCC patients. In total, 42 PRGs that were differentially expressed between HCC tissues and adjacent tissues and we exhibited the mutation frequency, classification, the location of copy number variation (CNV) alteration and the CNV variation frequency of PRGs. Two clusters were distinguished by the consensus clustering analysis based on the 42 differentially expressed genes (DEGs). There were significant differences in clinical features including T stage, grade, gender, and stage among different clusters. Kaplan-Meier curve analysis showed that cluster 1 had a better prognosis than cluster 2. The prognostic value of PRGs for survival was evaluated to construct a multigene signature using The Cancer Genome Atlas (TCGA) cohort. Based on the univariate analysis and multivariate analysis, a 10-gene signature was built and all HCC patients in the TCGA cohort were divided into low-risk group and high-risk group. HCC patients in the high-risk group showed significantly lower survival possibilities than those in the low-risk group (P < 0.001). Utilizing the median risk score from the TCGA cohort, HCC patients from International Cancer Genome Consortium (ICGC)-LIRI-JP cohort and Gene Expression Omnibus (GEO) cohort (GSE14520) were divided into two risk subgroups. The result showed that overall survival (OS) time was decreased in the high-risk group. Combined with the clinical characteristics, the risk score was an independent factor for predicting the OS of HCC patients. Then, ROC curve and survival analysis were performed to evaluate the prognostic prediction value of the model. Finally, we constructed a PRGs clinical characteristics nomogram to further predict HCC patient survival probability. There were significant differences in immune cell infiltration, GSEA enrichment pathway, IC50 of chemotherapeutics, PRGs mutation frequency between high-risk group and low-risk group. This work suggests PRGs signature played a crucial role in predicting the prognosis, infiltration of cancer microenvironment, and sensitivity of chemotherapeutic agents.
| S-EPMC9491002 | biostudies-literature