scholarly journals Identification of Disease-Specific Hub Biomarkers and Immune Infiltration in Osteoarthritis and Rheumatoid Arthritis Synovial Tissues by Bioinformatics Analysis

2021 ◽  
Vol 2021 ◽  
pp. 1-17
Author(s):  
Fei Sun ◽  
Jian lin Zhou ◽  
Pu ji Peng ◽  
Chen Qiu ◽  
Jia rui Cao ◽  
...  

Background. Osteoarthritis (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. Objective. This 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. Methods. Gene 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. Results. We 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. Conclusions. We 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.

2021 ◽  
Vol 12 ◽  
Author(s):  
Sheng Zhou ◽  
Hongcheng Lu ◽  
Min Xiong

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.


2021 ◽  
Vol 12 ◽  
Author(s):  
Rongguo Yu ◽  
Jiayu Zhang ◽  
Youguang Zhuo ◽  
Xu Hong ◽  
Jie Ye ◽  
...  

BackgroundRheumatoid arthritis (RA) refers to an autoimmune rheumatic disease that imposes a huge burden on patients and society. Early RA diagnosis is critical to preventing disease progression and selecting optimal therapeutic strategies more effectively. In the present study, the aim was at examining RA’s diagnostic signatures and the effect of immune cell infiltration in this pathology.MethodsGene Expression Omnibus (GEO) database provided three datasets of gene expressions. Firstly, this study adopted R software for identifying differentially expressed genes (DEGs) and conducting functional correlation analyses. Subsequently, we integrated bioinformatic analysis and machine-learning strategies for screening and determining RA’s diagnostic signatures and further verify by qRT-PCR. The diagnostic values were assessed through receiver operating characteristic (ROC) curves. Moreover, this study employed cell-type identification by estimating relative subsets of RNA transcript (CIBERSORT) website for assessing the inflammatory state of RA, and an investigation was conducted on the relationship of diagnostic signatures and infiltrating immune cells.ResultsOn the whole, 54 robust DEGs received the recognition. Lymphocyte-specific protein 1 (LSP1), Granulysin (GNLY), and Mesenchymal homobox 2 (MEOX2) (AUC = 0.955) were regarded as RA’s diagnostic markers and showed their statistically significant difference by qRT-PCR. As indicated from the immune cell infiltration analysis, resting NK cells, neutrophils, activated NK cells, T cells CD8, memory B cells, and M0 macrophages may be involved in the development of RA. Additionally, all diagnostic signatures might be different degrees of correlation with immune cells.ConclusionsIn conclusion, LSP1, GNLY, and MEOX2 are likely to be available in terms of diagnosing and treating RA, and the infiltration of immune cells mentioned above may critically impact RA development and occurrence.


2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Zitong Feng ◽  
Jingge Qu ◽  
Xiao Liu ◽  
Jinghui Liang ◽  
Yongmeng Li ◽  
...  

AbstractEsophageal squamous cell carcinoma (ESCC) is a life-threatening thoracic tumor with a poor prognosis. The role of molecular alterations and the immune microenvironment in ESCC development has not been fully elucidated. The present study aimed to elucidate key candidate genes and immune cell infiltration characteristics in ESCC by integrated bioinformatics analysis. Nine gene expression datasets from the Gene Expression Omnibus (GEO) database were analysed to identify robust differentially expressed genes (DEGs) using the robust rank aggregation (RRA) algorithm. Functional enrichment analyses showed that the 152 robust DEGs are involved in multiple processes in the tumor microenvironment (TME). Immune cell infiltration analysis based on the 9 normalized GEO microarray datasets was conducted with the CIBERSORT algorithm. The changes in macrophages between ESCC and normal tissues were particularly obvious. In ESCC tissues, M0 and M1 macrophages were increased dramatically, while M2 macrophages were decreased. A robust DEG-based protein–protein interaction (PPI) network was used for hub gene selection with the CytoHubba plugin in Cytoscape. Nine hub genes (CDA, CXCL1, IGFBP3, MMP3, MMP11, PLAU, SERPINE1, SPP1 and VCAN) had high diagnostic efficiency for ESCC according to receiver operating characteristic (ROC) curve analysis. The expression of all hub genes except MMP3 and PLAU was significantly related to macrophage infiltration. Univariate and multivariate regression analyses showed that a 7-gene signature constructed from the robust DEGs was useful for predicting ESCC prognosis. Our results might facilitate the exploration of potential targeted TME therapies and prognostic evaluation in ESCC.


2021 ◽  
Vol 2021 ◽  
pp. 1-9
Author(s):  
Ming Hu ◽  
Jianhua Wang

The purpose of this study was to identify hub genes closely correlated with Alzheimer's disease (AD) and their association with immune cell infiltration. In this work, 119 overlapping differentially expressed genes (DEGs) were obtained from GSE5281 and GSE122063 datasets through differential expression analysis. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed on the 119 DEGs, revealing some important biological functions and key pathways. AD immune cell infiltration analysis revealed a significant difference in the proportion of immune cells between the AD group and the control group. Finally, correlation analysis between target hub genes and immune cells indicated that GFAP had a positive or negative correlation with some specific immune cells. Our results provided useful clues, which will help to explain the molecular mechanism of AD and search for precise prognostic markers and potential therapeutic targets.


2021 ◽  
Author(s):  
shenglan li ◽  
Zhuang Kang ◽  
jinyi Chen ◽  
Can Wang ◽  
Zehao Cai ◽  
...  

Abstract Background Medulloblastoma is a common intracranial tumor among children. In recent years, research on cancer genome has established four distinct subtypes of medulloblastoma: WNT, SHH, Group3, and Group4. Each subtype has its own transcriptional profile, methylation changes, and different clinical outcomes. Treatment and prognosis also vary depending on the subtype. Methods Based on the methylation data of medulloblastoma samples, methylCIBERSORT was used to evaluate the level of immune cell infiltration in medulloblastoma samples and identified 10 kinds of immune cells with different subtypes. Combined with the immune database, 293 Imm-DEGs were screened. Imm-DEGs were used to construct the co-expression network, and the key modules related to the level of differential immune cell infiltration were identified. Three immune hub genes (GAB1, ABL1, CXCR4) were identified according to the gene connectivity and the correlation with phenotype in the key modules, as well as the PPI network involved in the genes in the modules. Results The subtype marker was recognized according to the immune hub, and the subtype marker was verified in the external data set, the methylation level of immune hub gene among different subtypes was compared and analyzed, at the same time, tissue microarray was used for immunohistochemical verification, and a multi-factor regulatory network of hub gene was constructed. Conclusions Identifying subtype marker is helpful to accurately identify the subtypes of medulloblastoma patients, and can accurately evaluate the treatment and prognosis, so as to improve the overall survival of patients.


2021 ◽  
Author(s):  
Di Cao ◽  
Jun Wang ◽  
Yan Lin ◽  
Guangwei Li

Abstract Background: The therapeutic efficacy of immune checkpoint inhibitor therapy is highly influenced by tumor mutation burden (TMB). The relationship between TMB and prognosis in lower-grade glioma is still unclear. We aimed to explore the relationships and mechanisms between them in lower-grade glioma.Methods: We leveraged somatic mutation data from The Cancer Genome Atlas (TCGA). Clinical cases were divided into high- and low-TMB groups based on the median of TMB. Infiltrating immune cells were analyzed using CIBERSORT and differential expression analysis between the prognostic groups performed. The key genes were identified as intersecting between immune-related genes. Cox regression and survival analysis were performed on the intersecting genes. A total of 7 hub genes were identified. The effect of somatic copy number alterations (SCNA) of the hub genes on immune cell infiltration was analyzed using TIMER, which was used to determine the risk factors and immune infiltration status in LGG. Subsequently, based on hub genes, a TMB Prognosis Index (TMBPI) model was constructed to predict the risk in LGG patients. Besides, this model was validated using data from TCGA and Chinese Glioma Genome Atlas (CGGA).Results: High-TMB favored worse prognosis (P<0.001) and macrophage infiltration was an independent risk factor (P<0.001). In the high-TMB group (P=0.033, P=0.009), the proportion of macrophages M0 and M2 increased and monocytes decreased (P=0.006). Besides, the SCNA of the hub genes affected the level of immune cell infiltration by varying degrees among which IGF2BP3, NPNT, and PLA2G2A had a significant impact. The AUC of the ROC curve at 1-, 3- and 5-years were all above 0.74.Conclusions: This study implies that high-TMB correlated with unfavorable prognosis in lower-grade glioma. And high-TMB may have an impact on prognosis by changing tumor microenvironment, caused by the SCNAs of genes. The TMBPI model accurately predicted prognosis in LGG patients.


2020 ◽  
Author(s):  
Jukun Song ◽  
Song He ◽  
Wei Wang ◽  
Jiaming Su ◽  
Dongbo Yuan ◽  
...  

Abstract Background Immune infiltration of Prostate cancer (PCa) was highly related to clinical outcomes. However, previous works failed to elucidate the diversity of different immune cell types that make up the function of the immune response system. The aim of the study was to uncover the composition of TIICs in PCa utilizing the CIBERSORT algorithm and further reveal the molecular characteristics of PCa subtypes. Method In the present work, we employed the CIBERSORT method to evaluate the relative proportions of immune cell profiling in PCa and adjacent samples, normal samples. We analyzed the correlation between immune cell infiltration and clinical information. The tumor-infiltrating immune cells of the TCGA PCa cohort were analyzed for the first time. The fractions of 22 immune cell types were imputed to determine the correlation between each immune cell subpopulation and clinical feature. Three types of molecular classification were identified via R-package of “CancerSubtypes”. The functional enrichment was analyzed in each subtype. The submap and TIDE algorithm were used to predict the clinical response to immune checkpoint blockade, and GDSC was employed to screen chemotherapeutic targets for the potential treatment of PCa. Results In current work, we utilized the CIBERSORT algorithm to assess the relative proportions of immune cell profiling in PCa and adjacent samples, normal samples. We investigated the correlation between immune cell infiltration and clinical data. The tumor-infiltrating immune cells in the TCGA PCa cohort were analyzed. The 22 immune cells were also calculated to determine the correlation between each immune cell subpopulation and survival and response to chemotherapy. Three types of molecular classification were identified. Each subtype has specific molecular and clinical characteristics. Meanwhile, Cluster I is defined as advanced PCa, and is more likely to respond to immunotherapy. Conclusions Our results demonstrated that differences in immune response may be important drivers of PCa progression and response to treatment. The deconvolution algorithm of gene expression microarray data by CIBERSOFT provides useful information about the immune cell composition of PCa patients. In addition, we have found a subtype of immunopositive PCa subtype and will help to explore the reasons for the poor effect of PCa on immunotherapy, and it is expected that immunotherapy will be used to guide the individualized management and treatment of PCa patients.


Sign in / Sign up

Export Citation Format

Share Document