scholarly journals POS0363 IDENTIFICATION OF MOLECULAR PHENOTYPES AND IMMUNE CELL INFILTRATION IN PSORIATIC ARTHRITIS PATIENTS’ SKIN TISSUES BY INTEGRATED BIOINFORMATICS ANALYSIS

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 80 (Suppl 1) ◽  
pp. 619.1-619
Author(s):  
S. Song ◽  
S. X. Zhang ◽  
J. Qiao ◽  
R. Zhao ◽  
J. Shi ◽  
...  

Background:Systemic lupus erythematosus (SLE) is a chronic autoimmune disease with highly heterogeneous clinical presentation characterized by disease unpredictable flares and multi-systemic involvement1 2. This clinical heterogeneity calls for design a molecular stratification to improve clinical trial design and formulate personalization treatment therapies.Objectives:This research was conducted to develop a reliable method to stratify SLE patients combined gene expression information and disease status.Methods:The mRNA expression profile of GSE138458 (contained 307 patients and 23 controls) and GSE49454 (contained 111 patients and 16 controls) were downloaded from the publicly GEO databases. After background adjustment, batch correction, and other pre-procession, obtaining a big gene matrix to identify the differentially expressed genes (DEGs) in SLE compared with healthy controls, which were screened by P value < 0.01. SLE subtypes were identified by non-negative matrix factorization (NMF) based on DEGs. Acquired signature genes in different SLE subtypes were conducted to process pathway enrichment analysis in Metascape. SLEDAI score and immune cell infiltration was also performed between subtypes by software package R (version 4.0.3).Results:Total 1202 DEGs were imputed to NMF unsupervised machine learning method. Patients with SLE were stratified into two subsets based on 184 signature genes derived from obtained DEGs(Fig.1A, 1B). GO and KEGG enrichment analysis showed that signature genes were mainly involved in negative regulation of innate immune response, toll-like receptor signaling pathway, regulation of immune effector process and so on(Fig.1C). Patients in Sub1 group had severe disease activity measures compared with those in Sub2(Fig.1D). SLEDAI scores from GSE49454 dataset were also higher in Sub1 compare with Sub2(Fig.1E). Further, immune cell infiltration results revealed an insufficient of regulatory T cell, CD8 T cells and naive CD4 T cells in Sub1 and neutrophils cells in Sub2(P<0.05)(Fig.1F).Conclusion:Our findings indicate that patients with SLE could be stratified into 2 subtypes which had different lymphocyte status and closely related to disease activity. This phenotyping may help us understand the etiology of the disease, inform patient in the design of clinical trials and guide treatment decision.References:[1]Dorner T, Furie R. Novel paradigms in systemic lupus erythematosus. Lancet 2019;393(10188):2344-58. doi: 10.1016/S0140-6736(19)30546-X [published Online First: 2019/06/11].[2]Fanouriakis A, Tziolos N, Bertsias G, et al. Update οn the diagnosis and management of systemic lupus erythematosus. Annals of the rheumatic diseases 2021;80(1):14-25. doi: 10.1136/annrheumdis-2020-218272 [published Online First: 2020/10/15].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


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.


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 aimed to characterize infiltrating immune cells and genes associated with the 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 Gens and Genomes (KEGG) analysis were applied 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, which suggest novel molecular mechanisms of sarcoma metastasis.


2021 ◽  
Vol 12 (1) ◽  
Author(s):  
Li Zhang ◽  
Shangshang Hu ◽  
Jiasheng Chen ◽  
Shasha Ma ◽  
Fanghong Liu ◽  
...  

AbstractA growing number of studies have shown that competitive endogenous RNA (ceRNA) regulatory networks might play important roles during the process of hepatocellular carcinoma (HCC). This study assessed the role of the ceRNA network in immune cell infiltration in HCC. Immune-related gene sets were downloaded from Molecular Signatures Database, and differentially expressed genes were screened based on TCGA HCC transcriptome data. The corresponding miRNAs with low expression and good prognostic implications, and the corresponding lncRNAs with high expression and poor prognostic were identified to construct ceRNA networks. The networks were utilized for clinical correlation analysis and risk model construction, and the CIBERSORT algorithm was applied to assess immune cell infiltration. In this study, the mRNA-miRNA-lncRNA model was used to construct a ceRNA network in HCC using immune-related differentially expressed mRNAs. Assessment of the MIR4435-2HG/hsa-miR-1-3p/MMP9/hsa-miR-29-3p/DUXAP8 ceRNA network axis in HCC showed that a high risk/poor prognosis was significantly correlated with tumor stage and invasion depth. MMP9 was positively correlated with resting M0 macrophages and NK cells and negatively correlated with activated mast cells, resting mast cells, monocytes and activated NK cells. DUXAP8 was positively correlated with M2 macrophages and negatively correlated with MIR4435-2HG, which was positively correlated with M2 macrophages and negatively correlated with activated mast cells, CD8 T cells and follicular helper T cells. The correlation of the MIR4435-2HG/hsa-miR-1-3p/MMP9/hsa-miR-29-3p/DUXAP8 ceRNA network axis with immune cell infiltration provides further information on the mechanism of HCC development. The result might improve our understanding the interactions between immune related genes and non-coding RNAs in the occurrence and development of HCC, and the relevant RNAs might be used as diagnostic and prognostic biomarkers and molecular targets in HCC patients.


2021 ◽  
Vol 8 ◽  
Author(s):  
Yi Liu ◽  
Juan Xiang ◽  
Gang Peng ◽  
Chenfu Shen

PDZ-binding kinase (PBK) is known to regulate tumor progression in some cancer types. However, its relationship to immune cell infiltration and prognosis in different cancers is unclear. This was investigated in the present study by analyzing data from TCGA, GEO, GETx, TIMER, CPTAC, GEPIA2, cBioPortal, GSCALite, PROGNOSCAN, PharmacoDB, STRING, and ENCORI databases. PBK was overexpressed in most tumors including adenocortical carcinoma (hazard ratio [HR] = 2.178, p &lt; 0.001), kidney renal clear cell carcinoma (KIRC; HR = 1.907, p &lt; 0.001), kidney renal papillary cell carcinoma (HR = 3.024, p &lt; 0.001), and lung adenocarcinoma (HR = 1.255, p &lt; 0.001), in which it was associated with poor overall survival and advanced pathologic stage. PBK methylation level was a prognostic marker in thyroid carcinoma (THCA). PBK expression was positively correlated with the levels of BIRC5, CCNB1, CDC20, CDK1, DLGAP5, MAD2L1, MELK, PLK1, TOP2A, and TTK in 32 tumor types; and with the levels of the transcription factors E2F1 and MYC, which regulate apoptosis, the cell cycle, cell proliferation and invasion, tumorigenesis, and metastasis. It was also negatively regulated by the microRNAs hsa-miR-101-5p, hsa-miR-145-5p, and hsa-miR-5694. PBK expression in KIRC, liver hepatocellular carcinoma, THCA, and thymoma was positively correlated with the infiltration of immune cells including B cells, CD4+T cells, CD8+ T cells, macrophages, monocytes, and neutrophils. The results of the functional enrichment analysis suggested that PBK and related genes contribute to tumor development via cell cycle regulation. We also identified 20 drugs that potentially inhibit PBK expression. Thus, PBK is associated with survival outcome in a variety of cancers and may promote tumor development and progression by increasing immune cell infiltration into the tumor microenvironment. These findings indicate that PBK is a potential therapeutic target and has prognostic value in cancer treatment.


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 ◽  
Vol 20 (1) ◽  
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 ◽  
Vol 2020 ◽  
pp. 1-17
Author(s):  
Pan Liao ◽  
Shuzhen Han ◽  
Honglan Qu

Objective. We investigated the expression patterns, potential functions, unique prognostic value, and potential therapeutic targets of E2Fs in brain and CNS cancer and tumor-infiltrating immune cell microenvironments. Methods. We analyzed E2F mRNA expression levels in diverse cancer types via Oncomine and GEPIA databases, respectively. Moreover, we evaluated the prognostic values using GEPIA database and TCGAportal database and the correlation of E2F expression with immune infiltration and the correlation between immune cell infiltration and GBM and LGG prognosis via TIMER database. Then, cBioPortal, GeneMANIA, and DAVID databases were used for mutation analysis, PPI network analysis of coexpressed gene, and functional enrichment analysis. Results. E2F1-8 expression increased in most cancers, including brain and CNS cancer. Higher expression in E2F1, 2, 4, 6, 7, and 8 indicated poor OS of LGG. Higher E2F3–6 and E2F1–8 expressions correlated with poor prognosis and increased immune infiltration levels in CD8+ T cells, macrophages, neutrophils, and DCs in GBM and CD8+ T cells, B cells, CD4+ T cells, neutrophils, macrophages, and DCs in LGG, respectively. Conclusion. E2F1–8 and E2F2–8 could be hopeful prognostic biomarkers of GBM and LGG, respectively. E2F3–6 and E2F1–8 could be likely therapeutic targets in patients with immune cell infiltration of GBM and LGG, respectively.


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.


Sign in / Sign up

Export Citation Format

Share Document