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

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.

PeerJ ◽  
2020 ◽  
Vol 8 ◽  
pp. e9996
Author(s):  
Yongyong Wang ◽  
Jianji Guo

Background Squamous cell lung carcinoma (LUSC) was closely associated with smoking which was known to have a distant immunosuppression effect. In this study, we aimed to explore the relationship between immune cells and clinical outcomes of LUSC patients with smoking history. Methods The immune cell infiltration and RNA expression profiles of LUSC patients were collected from The Cancer Genome Atlas (TCGA). Then, the correlation between immune cell infiltration and clinical characteristics was explored. According to the level of immune cell infiltration, LUSC patients with smoking history were divided into high or low group to screen the differentially expressed lncRNAs and mRNAs. The prediction of target genes was performed by miRanda. Finally, the prognostic value of a certain signature was confirmed in an independent dataset. Results Higher abundance of tumor-infiltrating T follicular helper (Tfh) cells together with a lower abundance of resting memory CD4 T cells had been found in LUSC current reformed smokers for ≤15 years and current smoking patients. Moreover, Tfh cell infiltration was not only associated with better overall survival (OS) but also varied from different degrees of TNM stage. Low expression of lncRNA PWRN1 and its potential regulating genes DMRTB1, PIRT, APOBEC1, and ZPBP2 were associated with better OS. Combining PWRN1 and four regulating genes as a signature, patients with higher-level expression of the signature had shorter survival time in not only the TCGA but also in the GEO dataset. Conclusions It was found that Tfh cells presented higher infiltration in LUSC current reformed smokers for ≤15 years and current smokers, while resting memory CD4 T cells had lower infiltration. The signature consisting of PWRN1 as well as its predicted targeted mRNAs was dysregulated in different levels of Tfh cell infiltration and might indicate patients’ OS.


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.


2020 ◽  
Author(s):  
Zeyu Yang ◽  
Tianjing Du ◽  
Qiao Xiong ◽  
Weiwei Zhang ◽  
Chao Wang ◽  
...  

Abstract Background The role of immune cell infiltration in tumor biology and the potential of immunotherapy for the treatment of several cancers have been proven. However, the immunogenomic landscape and immune cell infiltration need to be comprehensively analyzed in bladder cancer (BC). This study aimed to explore the immune-related genes (IRGs) in BC to create a prognostic risk assessment model and gain some insights into the molecular underpinnings of BC. Methods Based on the datasets retrieved from The Cancer Genome Atlas (TCGA) database, we identified survival-associated IRGs via univariate Cox analysis. Then, we created an immune-related gene-based prognostic factor (IRGPF) and validated it by multivariable Cox analysis. We displayed the profiles of 22 types of immune cells by using CIBERSORT and explored the relationship between IRGs and immune cell infiltration. Results Altogether, 58 differentially expressed IRGs were found to be associated with the prognosis of patients with BC. We constructed a prognostic assessment model as an IRGPF with IRGs (THBS1, CXCL9, CXCL11, FABP6, BIRC5, and PPY). Profiles of the infiltrating immune cells confirmed their significance based on clinical factors and individual differences. The IRGPF was related to immune cell infiltration, and the key gene was identified as THBS1. Conclusions Our findings confirmed that IRGs could act as independent prognostic factors and immune-driver factors. Patients with high levels of activated memory CD4 T cells but low levels of resting memory CD4 T cells had a better prognosis. This study indicates the possibility of developing new immunotherapeutic strategies and individualized treatment based on this approach.


2021 ◽  
Author(s):  
Meng Wang ◽  
Ruijie Zhang ◽  
Qiongfeng Guan ◽  
Yindan Yao ◽  
Liyuan Han

Abstract Background: This study aimed to identify potential diagnostic markers of ischemic stroke (IS) and discuss the function of immune cell infiltration during the pathological process. Methods: We used IS datasets from the Gene Expression Omnibus database. Differentially expressed genes (DEGs) were identified, and functional correlation analysis was performed. We then screened and verified the diagnostic markers of IS. We evaluated the infiltration of immune cells in infarcts using CIBERSORT and analyzed the correlation between diagnostic markers and infiltrating immune cells. Results: A total of 366 DEGs were screened in this study. Genes encoding CTSG, F13A1, PABPC1, ECHDC2, BIRC2 and infiltrating monocytes, M0 macrophages, activated dendritic cells, and neutrophils (area under the curve [AUC] = 0.945) were identified as diagnostic markers of IS. Immune cell infiltration analysis suggested that memory B cells, regulatory T cells, M0 macrophages, CD8 + T cells, γδT cells, activated natural killer cells, monocytes, activated mast cells, and neutrophils were involved in the IS process. Analysis of correlations between expressed genes and infiltrating immune cells found that CTSG was positively associated with M0 macrophages, F13A1 was positively associated with monocytes, PABPC1 was positively associated with activated dendritic cells, eosinophils were negatively associated with neutrophils, ECHDC2 was negatively associated with monocytes, and BIRC2 was positively associated with eosinophils. Conclusion: five genes and four types of immune cells were identified as diagnostic markers of IS, and immune cell infiltration may play an important role in the progression of IS.


2021 ◽  
Author(s):  
Zhihao Chen ◽  
Liubing Li ◽  
Ziyuan Li ◽  
Xi Wang ◽  
Mingxiao Han ◽  
...  

Abstract Background: The potential functions of circular RNAs (circRNAs) and micro RNAs (miRNAs) in osteosarcoma (OS) have not been fully elucidated. Especially, the behavior and mechanism of immune responses in OS development and progression have not been fully demonstrated. It was reported that circRNAs and miRNAs can serve as biomarkers for the diagnosis, prognosis, and therapy of many cancers. This study aimed to identify novel key serum biomarkers to diagnose and predict metastasis of OS based on the analysis of immune cell infiltration characteristics.Methods: The differentially-expressed circRNAs (DEcircRNAs), differentially-expressed miRNAs (DEmiRNAs),and differentially-expressed mRNAs (DEmRNAs) of human OS were investigated based on the microarray data downloaded from Gene Expression Omnibus (GEO) datasets. Then, we analyzed immune characteristics pattern of tumor-infiltrating immune cells in OS. On this basis, we identified statistically-significant transcription factors and performed pathway enrichment analysis. Subsequently, we constructed protein-protein interaction (PPI) and competitive endogenous RNA (ceRNA) networks. Moreover, the biological characteristic of targets in ceRNA networks was proposed. Finally, the expression and diagnostic capability of these potential biomarkers from ceRNA network were confirmed by RT-qPCR in patients’ serum.Results: Seven differentially-expressed circRNAs (DEcircRNAs), 166 differentially-expressed miRNAs (DEmiRNAs) and 175 differentially-expressed mRNAs (DEmRNAs) were identified in total. The highest level of infiltration in OS patients were M0 macrophages, M2 macrophages and CD8+ T cells. Further, M0 macrophages and CD8+ T cells were showed the largest negative correlation coefficients. These significant immune characteristics pattern of tumor-infiltrating immune cells were revealed by the principal component analysis in OS. Moreover, we found 185 statistically-significant transcription factors in which the main significant molecules show the potential in immunotherapy of OS. Hsa-circ-0010220, hsa-miR-326, hsa-miR-338-3p, and FAM98A from ceRNA networks associated with immune cell infiltration were confirmed as the potential novel biomarkers for OS diagnosis, of which FAM98A could distinguish and predict metastasis. Most importantly, a novel diagnostic model consisting of the four promising biomarkers (hsa-circ-0010220, hsa-miR-326, hsa-miR-338-3p, and FAM98A) was highlighted with 0.928 AUC value.Conclusions: In summary, the potenial serum biomarkers to diagnose and predict metastasis of OS based on the analysis of immune cell infiltration characteristics were found, and a novel diagnostic model consisting of four promising serum biomarkers was proposed firstly. These results provided a new perspective for the immunotherapy of OS.


2021 ◽  
Vol 80 (Suppl 1) ◽  
pp. 411.1-411
Author(s):  
T. Cheng ◽  
S. X. Zhang ◽  
J. Qiao ◽  
R. Zhao ◽  
S. Song ◽  
...  

Background:Psoriatic arthritis (PsA) is an inflammatory musculoskeletal disease associated with cutaneous psoriasis1. Heterogeneity of clinical manifestation often makes differential diagnosis difficult 2. Thus, the underlying molecular pathogenesis of PsA need to be further studied to diagnose early and ensure optimal management of arthritis and key comorbidities.Objectives:This research was conducted to identify molecular phenotypes and immune infiltration in the skin tissues of psoriatic arthritis patients according to bioinformatics analysis.Methods:The mRNA expression profiles of GSE13355 (116 samples), GSE14905 (56 samples) and GSE30999 (162 samples) were obtained from the publicly GEO databases. Non-negative matrix factorization (NMF), functional enrichment and cibersort algorithm were applied to illustrate the conditions of PsA patients’ skin tissues for classification after screening the differentially expressed genes (DEGs) between lesion biopsy and non-lesion biopsy.Results:Two subsets (Sub1 and Sub2) were identified and validated by NMF typing of 612 detected DEGs (Figure 1a). A total of 54 signature genes (18 in Sub1 and 36 in Sub2) were obtained (Figure 1b). GO and KEGG enrichment analysis showed the signature genes in Sub1 were mainly involved in proliferation and differentiation of immune cells, whereas genes in Sub2 were related to humoral immune response mediated by antimicrobial peptide (Figure 1c.1d). Further, immune cell infiltration results revealed Sub2 had higher levels of resting NK cells (P<0.001), macrophages M1(P<0.001), resting mast cells (P<0.001) and regulatory T cells (P<0.001) but lower concentrations of activated CD4+ memory T cells (P<0.001), activated NK cells (P<0.05), activated dendritric cells(P<0.001), eosinophils (P<0.05) and neutrophil (P<0.001) (Figure 1e).Conclusion:The pathogenesis of psoriatic arthritis is related to both cellular immunity and humoral immunity. It is indispensable to adjust the treatment strategies according to patient’s immune status.References:[1]Ritchlin CT, Colbert RA, Gladman DD. Psoriatic Arthritis. The New England journal of medicine 2017;376(10):957-70. doi: 10.1056/NEJMra1505557 [published Online First: 2017/03/09].[2]Veale DJ, Fearon U. The pathogenesis of psoriatic arthritis. Lancet (London, England) 2018;391(10136):2273-84. doi: 10.1016/s0140-6736(18)30830-4 [published Online First: 2018/06/13].Acknowledgements:This project was supported by National Science Foundation of China (82001740), Open Fund from the Key Laboratory of Cellular Physiology (Shanxi Medical University) (KLCP2019) and Innovation Plan for Postgraduate Education in Shanxi Province (2020BY078).Disclosure of Interests:None declared


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.


2020 ◽  
Author(s):  
Naiqiang Zhu ◽  
Jingyi Hou

Abstract Background: Sarcomas, cancers originating from mesenchymal cells, are comprehensive tumors with poor prognoses, yet their tumorigenic mechanisms are largely unknown. In this study, we characterize infiltrating immune cells and analyze immune scores to identify the molecular mechanism of immunologic response to sarcomas.Method: The “CIBERSORT” algorithm was used to calculate the amount of L22 immune cell infiltration in sarcomas. Then, the “ESTIMATE” algorithm was used to assess the “Estimate,” “Immune,” and “Stromal” scores. Weighted gene co-expression network analysis (WGCNA) was utilized to identify the significant module related to the immune therapeutic target. Gene ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed using the “clusterProfiler” package in R for annotation and visualization. Results: Macrophages were the most common immune cells infiltrating sarcomas. The number of CD8 T cells was negatively associated with that of M0 and M2 macrophages, and positively associated with M macrophages in sarcomas samples. The clinical parameters (disease type, gender) significantly increased with higher Estimate, Immune, and Stromal scores, and with a better prognosis. The blue module was significantly associated with CD8 T cells. Functional enrichment analysis showed that the blue module was mainly involved in chemokine signaling and the PI3K-Akt signaling pathway. CD48, P2RY10 and RASAL3 were identified and validated at the protein level.Conclusion: Based on the immune cell infiltration and immune microenvironment, three key genes were identified, thus presenting novel molecular mechanisms of sarcoma metastasis.


2020 ◽  
Author(s):  
Naiqiang Zhu ◽  
Jingyi Hou

Abstract Background: Sarcomas, cancers originating from mesenchymal cells, are comprehensive tumors with poor prognoses, yet their tumorigenic mechanisms are largely unknown. In this study, we characterize infiltrating immune cells and analyze immune scores to identify the molecular mechanism of immunologic response to sarcomas.Method: The “CIBERSORT” algorithm was used to calculate the amount of L22 immune cell infiltration in sarcomas. Then, the “ESTIMATE” algorithm was used to assess the “Estimate,” “Immune,” and “Stromal” scores. Weighted gene co-expression network analysis (WGCNA) was utilized to identify the significant module related to the immune therapeutic target. Gene ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed using the “clusterProfiler” package in R for annotation and visualization. Results: Macrophages were the most common immune cells infiltrating sarcomas. The number of CD8 T cells was negatively associated with that of M0 and M2 macrophages, and positively associated with M macrophages in sarcomas samples. The clinical parameters (disease type, gender) significantly increased with higher Estimate, Immune, and Stromal scores, and with a better prognosis. The blue module was significantly associated with CD8 T cells. Functional enrichment analysis showed that the blue module was mainly involved in chemokine signaling and the PI3K-Akt signaling pathway. CD48, P2RY10 and RASAL3 were identified and validated at the protein level.Conclusion: Based on the immune cell infiltration and immune microenvironment, three key genes were identified, thus presenting novel molecular mechanisms of sarcoma metastasis.


2021 ◽  
Vol 23 (Supplement_4) ◽  
pp. iv8-iv8
Author(s):  
Evyn Woodhouse ◽  
Liyam Laraba ◽  
Charlotte Lespade ◽  
Marie Srotyr ◽  
Alison C Lloyd ◽  
...  

Abstract Aims Previous work has shown that increased numbers of macrophages are associated with more rapid schwannoma tumour growth and we are interested in signals that control entry of macrophages and other immune cells into these tumours. Activation of the Raf-kinase domain and the Raf/MEK/ERK pathway within Schwann cells has been observed to induce an inflammatory response in peripheral nerves in the absence of injury. Activation of an inducible Raf-kinase transgene in Schwann cells allows modelling of acute demyelination of peripheral nerves without nerve injury. This Raf-oestrogen receptor fusion protein (Raf-TR) is activated by the oestrogen analogue Tamoxifen and so allows targeted, controlled activation of the Raf/MEK/ERK pathway within the Schwann cells. Here, in order to understand drivers of tumour formation, we assess the effect of MAPK activation in Merlin-null Schwann cells upon immune cell infiltration within the PNS. Method RafTR-P0CRE-NF2fl/fl mice of 4-6 weeks age were injected daily (IP) with 2mg of 4-hydroxy-tamoxifen or vehicle (corn oil) control for 5 consecutive days. RafTR was activated on either a Merlin (NF2) wild-type (NF2 fl/fl, P0-CRE-) or NF2 null (NF2 fl/fl, P0-CRE+) background and effects on immune cell infiltration studied in each condition. Immunofluorescence was performed in the dorsal root ganglia (DRGs) and sciatic nerves of mice to identify various immune cell infiltrates at various timepoints. These will include neutrophils, mast cells, T-Cells and macrophages using the cell markers Csf3r, C-kit, CD3 and IBA1 respectively. Results At 21 days post treatment, a significantly increased infiltration of macrophages within the sciatic nerve and dorsal root ganglia was observed in mice treated with Tamoxifen when compared to vehicle controls. Loss of NF2 led to a massive increase in the number of macrophages recruited to peripheral nerves in tamoxifen-treated mice compared to Cre- mice and Cre+ treated with vehicle alone. Further assessment of other immune cell infiltration including neutrophils, mast cells and T cells are ongoing. Conclusion Raf/MEK/ERK signalling, in the absence of tumour suppressor Merlin, significantly increases the infiltration of inflammatory cells such as macrophages into peripheral nerves even in the absence of injury. As this effect is enhanced in NF2 null mice, this suggests that Merlin plays an important role in inhibiting the inflammatory response in peripheral nerves. It also suggests that Merlin could be involved in maintaining the blood nerve barrier (BNB), as in its absence the greater influx of immune cells into the nerves and DRGs suggests a more complete loss of BNB function than just activation of the Raf/MEK/ERK cascade alone.


Sign in / Sign up

Export Citation Format

Share Document