scholarly journals Identifying changes in immune cells and constructing prognostic models using immune-related genes in post-burn immunosuppression

PeerJ ◽  
2022 ◽  
Vol 10 ◽  
pp. e12680
Author(s):  
Peng Wang ◽  
Zexin Zhang ◽  
Bin Yin ◽  
Jiayuan Li ◽  
Cheng Xialin ◽  
...  

Background Burn patients are prone to infection as well as immunosuppression, which is a significant cause of death. Currently, there is a lack of prognostic biomarkers for immunosuppression in burn patients. This study was conducted to identify immune-related genes that are prognosis biomarkers in post-burn immunosuppression and potential targets for immunotherapy. Methods We downloaded the gene expression profiles and clinical data of 213 burn patients and 79 healthy samples from the Gene Expression Omnibus (GEO) database. Immune infiltration analysis was used to identify the proportion of circulating immune cells. Functional enrichment analyses were carried out to identify immune-related genes that were used to build miRNA-mRNA networks to screen key genes. Next, we carried out correlation analysis between immune cells and key genes that were then used to construct logistic regression models in GSE77791 and were validated in GSE19743. Finally, we determined the expression of key genes in burn patients using quantitative reverse transcription polymerase chain reaction (qRT-PCR). Results A total of 745 differently expressed genes were screened out: 299 were up-regulated and 446 were down-regulated. The number of Th-cells (CD4+) decreased while neutrophils increased in burn patients. The enrichment analysis showed that down-regulated genes were enriched in the T-cell activation pathway, while up-regulated genes were enriched in neutrophil activation response in burn patients. We screened out key genes (NFATC2, RORA, and CAMK4) that could be regulated by miRNA. The expression of key genes was related to the proportion of Th-cells (CD4+) and survival, and was an excellent predictor of prognosis in burns with an area under the curve (AUC) value of 0.945. Finally, we determined that NFATC2, RORA, and CAMK4 were down-regulated in burn patients. Conclusion We found that NFATC2, RORA, and CAMK4 were likely prognostic biomarkers in post-burn immunosuppression and potential immunotherapeutic targets to convert Th-cell dysfunction.

2020 ◽  
Author(s):  
An Zhi Zhang ◽  
Xin Yuan ◽  
Ya Li ◽  
Yu Fang Xie ◽  
Jiang Fen Li ◽  
...  

Abstract Background In recent years, immunotherapy has developed rapidly and has gradually become one of the important methods for treatment of gastric cancer. The research on immune cells and immune-related genes in the tumor microenvironment greatly encourages the development of immunotherapy. Methods : The devolution algorithm (CIBERSORT) was applied to infer the proportion of 22 immune infiltrating cells based on gene expression profiles of gastric cancer tissue, which were downloaded from TCGA and GEO databases. The TCGA database was utilized to analyze the differential expression of immune-related genes, and explore the potential molecular functions of these genes.ResultsWe have observed the enrichment of multiple immune cells in the microenvironment of gastric cancer. Some of these cells are closely related to Fuhrman grade and TNM staging. Survival analysis showed that the infiltration level of CD8 + T cells, activated CD4 + memory T cells and M2 macrophages was significantly related to the prognosis of gastric cancer patients. The functional enrichment analysis of immune-related genes revealed that these genes were mainly associated with cytokine activation and response. Four significant modules were screened by PPI network and 20 key genes were screened from the modules, and the expression levels of CALCR and PTH1R are strikingly related to the prognosis of gastric cancer patients.ConclusionsThe type and number of infiltrating immune cells in the microenvironment of gastric cancer, as well as immune-related genes are closely related to tumor progression, and can be used as important indicators for patient prognosis assessment.


BMC Cancer ◽  
2021 ◽  
Vol 21 (1) ◽  
Author(s):  
Wen Tan ◽  
Maomao Liu ◽  
Liangshan Wang ◽  
Yang Guo ◽  
Changsheng Wei ◽  
...  

Abstract Background Breast cancer is one of the most frequently diagnosed cancers among women worldwide. Alterations in the tumor microenvironment (TME) have been increasingly recognized as key in the development and progression of breast cancer in recent years. To deeply comprehend the gene expression profiling of the TME and identify immunological targets, as well as determine the relationship between gene expression and different prognoses is highly critical. Methods The stromal/immune scores of breast cancer patients from The Cancer Genome Atlas (TCGA) were employed to comprehensively evaluate the TME. Then, TME characteristics were assessed, overlapping genes of the top 3 Gene Ontology (GO) terms and upregulated differentially expressed genes (DEGs) were analyzed. Finally, through combined analyses of overall survival, time-dependent receiver operating characteristic (ROC), and protein-protein interaction (PPI) network, novel immune related genes with good prognosis were screened and validated in both TCGA and GEO database. Results Although the TME did not correlate with the stages of breast cancer, it was closely associated with the subtypes of breast cancer and gene mutations (CDH1, TP53 and PTEN), and had immunological characteristics. Based on GO functional enrichment analysis, the upregulated genes from the high vs low immune score groups were mainly involved in T cell activation, the external side of the plasma membrane, and receptor ligand activity. The top GO terms of the upregulated DEGs from the high vs low immune score groups exhibited better prognosis in breast cancer; 15 of them were related to good prognosis in breast cancer, especially CD226 and KLRC4-KLRK1. Conclusions High CD226 and KLRC4-KLRK1 expression levels were identified and validated to correlate with better overall survival in specific stages or subtypes of breast cancer. CD226, KLRC4-KLRK1 and other new targets seem to be promising avenues for promoting antitumor targeted immunotherapy in breast cancer.


2021 ◽  
Vol 21 (1) ◽  
Author(s):  
Ya-Nan Pi ◽  
Jun-Nan Guo ◽  
Ge Lou ◽  
Bin-Bin Cui

Abstract Background Cervical cancer (CC) is the leading cause of cancer-related death in women. A limited number of studies have investigated whether immune-prognostic features can be used to predict the prognosis of CC. This study aimed to develop an improved prognostic risk scoring model (PRSM) for CC based on immune-related genes (IRGs) to predict survival and determine the key prognostic IRGs. Methods We downloaded the gene expression profiles and clinical data of CC patients from the TCGA and GEO databases. The ESTIMATE algorithm was used to calculate the score for both immune and stromal cells. Differentially expressed genes (DEGs) in different subpopulations were analyzed by “Limma”. A weighted gene co-expression network analysis (WGCNA) was used to establish a DEG co-expression module related to the immune score. Immune-related gene pairs (IRGPs) were constructed, and univariate- and Lasso-Cox regression analyses were used to analyze prognosis and establish a PRSM. A log-rank test was used to verify the accuracy and consistency of the scoring model. Identification of the predicted key IRG was ensured by the application of functional enrichment, DisNor, protein–protein interactions (PPIs) and heatmap. Finally, we extracted the key prognostic immune-related genes from the gene expression data, validated the key genes by immunohistochemistry and analyzed the correlation between their expression and drug sensitivity. Results A new PRSM was developed based on 22 IRGPs. The prognosis of the low-risk group in the model group (P < 0.001) and validation group (P = 0.039) was significantly better than that in the high-risk group. Furthermore, M1 and M2 macrophages were highly expressed in the low-risk group. Retinoic acid-inducible gene-I (RIG-I)-like receptors (RLRs) and the Janus kinase-signal transducer and activator of transcription (JAK-STAT) signaling pathway were significantly enriched in the low-risk group. Three representative genes (CD80, CD28, and LCP2) were markers of CC prognosis. CD80 and CD28 may more prominent represent important indicators to improve patient prognosis. These key genes was positively correlated with drug sensitivity. Finally, we found that differences in the sensitivity to JNK inhibitors could be distinguished based on the use and risk grouping of this PRSM. Conclusions The prognostic model based on the IRGs and key genes have potential clinical significance for predicting the prognosis of CC patients, providing a foundation for clinical prognosis judgment and individualized treatment.


Nano LIFE ◽  
2019 ◽  
Vol 09 (01n02) ◽  
pp. 1940002
Author(s):  
Jichen Xu ◽  
Xianchun Zong ◽  
Qianshu Ren ◽  
Hongyu Wang ◽  
Lijuan Zhao ◽  
...  

The aim of this paper is to identify key genes in lung adenocarcinoma (LUAD) through weighted gene co-expression network analysis (WGCNA), and to further understand the molecular mechanism of LUAD. 107 gene expression profiles were downloaded from GSE10072 in the GEO database. We performed rigorous processing of the initial gene expression profile data. Subsequently, we used WGCNA to identify disease-driven modules and enforced functional enrichment analysis. The key genes were defined as the most connected genes in the driver module and were validated using the GSE75037 and TCGA database. GSE10072 removed 41 unpaired lung samples and 4 outliers. By analyzing the 62 samples using WGCNA, we obtained 26 modules and identified the brown and magenta modules as the driving modules for the LUAD. We found that the “Cell cycle”, “Oocyte meiosis” and “Progesterone-mediated oocyte maturation” pathways may be related to the occurrence of LUAD. GSE75037 removed 8 outlier and obtained 2909 differentially expressed genes (DEGs), 26 genes (9 genes in the brown module, 17 genes in the magenta module) overlap with key genes in the driver module. The results of the survival analysis suggest that 19 genes were significantly correlated with the patient’s survival time, including KPNA2, FEN1, RRM2, TOP2A, CENPF, MCM4, BIRC5, MELK, MAD2L1, CCNB1, CCNA2, KIF11, CDKN3, NUSAP1, CEP55, AURKA, NEK2, KIF14 and CDCA8, which may be potential biomarkers or therapeutic targets for LUAD. In this study, we provide a theoretical basis for further understanding the biological mechanism of LUAD through bioinformatics analysis of LUAD.


2021 ◽  
Author(s):  
Xiong-Wen Wang ◽  
Qian Yan ◽  
Bao-Qian Ye ◽  
Bo-Qing Wang ◽  
Wen-Jiang Zheng

Abstract Background: The combination of epigenetic drugs and immunotherapy should be able to develop an optimal treatment plan for hepatocellular carcinoma (HCC), yet its mechanism is still in the preliminary exploration stage. The purpose of this study is to analyze the DNA methylation and gene expression profiles of immune-related CpG sites to identify the molecular subtypes and CpG sites related to the prognosis of HCC. Methods: In this study, the DNA methylation and gene expression datasets were downloaded from The Cancer Genome Atlas database, together with immune-related genes downloaded from the immunology database and analysis portal database to explore the prognostic molecular subtypes of HCC. Univariate and multivariate survival analysis was used for selecting the significant methylation sites, and the consensus clustering was performed to find the best molecular subtype associated with the survival of HCC. Next, we used the least absolute shrinkage and selection operator (LASSO) algorithm to construct a prognostic-related model and performed internal verification. Finally, we explored the levels of 16 immune-related genes expression correlate with the infiltration levels of immune cells in HCC. Results: By performing consistent clustering analysis on 830 immune-related CpG sites in 231 samples of a training set, we identified seven subgroups with significant differences in overall survival. Finally, 16 classifiers of immune-related CpG sites were constructed and used in the testing set to verify the prognosis of DNA methylation subgroups, and the results were consistent with the training set. Using the TIMER database, we analyzed 16 immune-related CpG sites expression with the abundance of six types of immune infiltrating cells and found that most are positively correlated with the level of infiltration of multiple immune cells in HCC. Conclusions: This study screened potential immune-related prognostic methylation sites and established a new prognosis model of HCC based on DNA methylation molecular subtype, which may help in the early diagnosis of HCC and developing more effective personalized treatments.


2021 ◽  
Vol 11 (11) ◽  
pp. 1177
Author(s):  
Shao-Hua Yu ◽  
Jia-Hua Cai ◽  
De-Lun Chen ◽  
Szu-Han Liao ◽  
Yi-Zhen Lin ◽  
...  

The aim of this study is to identify potential biomarkers for early diagnosis of gynecologic cancer in order to improve survival. Cervical cancer (CC) and endometrial cancer (EC) are the most common malignant tumors of gynecologic cancer among women in the world. As the underlying molecular mechanisms in both cervical and endometrial cancer remain unclear, a comprehensive and systematic bioinformatics analysis is required. In our study, gene expression profiles of GSE9750, GES7803, GES63514, GES17025, GES115810, and GES36389 downloaded from Gene Expression Omnibus (GEO) were utilized to analyze differential gene expression between cancer and normal tissues. A total of 78 differentially expressed genes (DEGs) common to CC and EC were identified to perform the functional enrichment analyses, including gene ontology and pathway analysis. KEGG pathway analysis of 78 DEGs indicated that three main types of pathway participate in the mechanism of gynecologic cancer such as drug metabolism, signal transduction, and tumorigenesis and development. Furthermore, 20 diagnostic signatures were confirmed using the least absolute shrink and selection operator (LASSO) regression with 10-fold cross validation. Finally, we used the GEPIA2 online tool to verify the expression of 20 genes selected by the LASSO regression model. Among them, the expression of PAMR1 and SLC24A3 in tumor tissues was downregulated significantly compared to the normal tissue, and found to be statistically significant in survival rates between the CC and EC of patients (p < 0.05). The two genes have their function: (1.) PAMR1 is a tumor suppressor gene, and many studies have proven that overexpression of the gene markedly suppresses cell growth, especially in breast cancer and polycystic ovary syndrome; (2.) SLC24A3 is a sodium–calcium regulator of cells, and high SLC24A3 levels are associated with poor prognosis. In our study, the gene signatures can be used to predict CC and EC prognosis, which could provide novel clinical evidence to serve as a potential biomarker for future diagnosis and treatment.


2021 ◽  
Vol 9 (Suppl 1) ◽  
pp. A12.1-A12
Author(s):  
Y Arjmand Abbassi ◽  
N Fang ◽  
W Zhu ◽  
Y Zhou ◽  
Y Chen ◽  
...  

Recent advances of high-throughput single cell sequencing technologies have greatly improved our understanding of the complex biological systems. Heterogeneous samples such as tumor tissues commonly harbor cancer cell-specific genetic variants and gene expression profiles, both of which have been shown to be related to the mechanisms of disease development, progression, and responses to treatment. Furthermore, stromal and immune cells within tumor microenvironment interact with cancer cells to play important roles in tumor responses to systematic therapy such as immunotherapy or cell therapy. However, most current high-throughput single cell sequencing methods detect only gene expression levels or epigenetics events such as chromatin conformation. The information on important genetic variants including mutation or fusion is not captured. To better understand the mechanisms of tumor responses to systematic therapy, it is essential to decipher the connection between genotype and gene expression patterns of both tumor cells and cells in the tumor microenvironment. We developed FocuSCOPE, a high-throughput multi-omics sequencing solution that can detect both genetic variants and transcriptome from same single cells. FocuSCOPE has been used to successfully perform single cell analysis of both gene expression profiles and point mutations, fusion genes, or intracellular viral sequences from thousands of cells simultaneously, delivering comprehensive insights of tumor and immune cells in tumor microenvironment at single cell resolution.Disclosure InformationY. Arjmand Abbassi: None. N. Fang: None. W. Zhu: None. Y. Zhou: None. Y. Chen: None. U. Deutsch: None.


PeerJ ◽  
2019 ◽  
Vol 7 ◽  
pp. e7821 ◽  
Author(s):  
Xiaoming Zhang ◽  
Jing Zhuang ◽  
Lijuan Liu ◽  
Zhengguo He ◽  
Cun Liu ◽  
...  

Background Cumulative evidence suggests that long non-coding RNAs (lncRNAs) play an important role in tumorigenesis. This study aims to identify lncRNAs that can serve as new biomarkers for breast cancer diagnosis or screening. Methods First, the linear fitting method was used to identify differentially expressed genes from the breast cancer RNA expression profiles in The Cancer Genome Atlas (TCGA). Next, the diagnostic value of all differentially expressed lncRNAs was evaluated using a receiver operating characteristic (ROC) curve. Then, the top ten lncRNAs with the highest diagnostic value were selected as core genes for clinical characteristics and prognosis analysis. Furthermore, core lncRNA-mRNA co-expression networks based on weighted gene co-expression network analysis (WGCNA) were constructed, and functional enrichment analysis was performed using the Database for Annotation, Visualization and Integrated Discovery (DAVID). The differential expression level and diagnostic value of core lncRNAs were further evaluated by using independent data set from Gene Expression Omnibus (GEO). Finally, the expression status and prognostic value of core lncRNAs in various tumors were analyzed based on Gene Expression Profiling Interactive Analysis (GEPIA). Results Seven core lncRNAs (LINC00478, PGM5-AS1, AL035610.1, MIR143HG, RP11-175K6.1, AC005550.4, and MIR497HG) have good single-factor diagnostic value for breast cancer. AC093850.2 has a prognostic value for breast cancer. AC005550.4 and MIR497HG can better distinguish breast cancer patients in early-stage from the advanced-stage. Low expression of MAGI2-AS3, LINC00478, AL035610.1, MIR143HG, and MIR145 may be associated with lymph node metastasis in breast cancer. Conclusion Our study provides candidate biomarkers for the diagnosis and prognosis of breast cancer, as well as a bioinformatics basis for the further elucidation of the molecular pathological mechanism of breast cancer.


PeerJ ◽  
2020 ◽  
Vol 8 ◽  
pp. e8831 ◽  
Author(s):  
Xiaojiao Guan ◽  
Yao Yao ◽  
Guangyao Bao ◽  
Yue Wang ◽  
Aimeng Zhang ◽  
...  

Esophageal cancer is a common malignant tumor in the world, and the aim of this study was to screen key genes related to the development of esophageal cancer using a variety of bioinformatics analysis tools and analyze their biological functions. The data of esophageal squamous cell carcinoma from the Gene Expression Omnibus (GEO) were selected as the research object, processed and analyzed to screen differentially expressed microRNAs (miRNAs) and differential methylation genes. The competing endogenous RNAs (ceRNAs) interaction network of differentially expressed genes was constructed by bioinformatics tools DAVID, String, and Cytoscape. Biofunctional enrichment analysis was performed using Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG). The expression of the screened genes and the survival of the patients were verified. By analyzing GSE59973 and GSE114110, we found three down-regulated and nine up-regulated miRNAs. The gene expression matrix of GSE120356 was calculated by Pearson correlation coefficient, and the 11696 pairs of ceRNA relation were determined. In the ceRNA network, 643 lncRNAs and 147 mRNAs showed methylation difference. Functional enrichment analysis showed that these differentially expressed genes were mainly concentrated in the FoxO signaling pathway and were involved in the corresponding cascade of calcineurin. By analyzing the clinical data in The Cancer Genome Atlas (TCGA) database, it was found that four lncRNAs had an important impact on the survival and prognosis of esophageal carcinoma patients. QRT-PCR was also conducted to identify the expression of the key lncRNAs (RNF217-AS1, HCP5, ZFPM2-AS1 and HCG22) in ESCC samples. The selected key genes can provide theoretical guidance for further research on the molecular mechanism of esophageal carcinoma and the screening of molecular markers.


2021 ◽  
Author(s):  
Jian Lin ◽  
Yuanhua Lu ◽  
Bizhou Wang ◽  
Ping Jiao ◽  
Jie Ma

Abstract Background Type 1 diabetes mellitus (T1DM) is a chronic autoimmune disease caused by severe loss of pancreatic β cells. Immune cells are key mediators of β cell destruction. This study attempted to investigate the role of immune cells and immune-related genes in the occurrence and development of T1DM. Methods The raw gene expression profile of the samples from 12 T1DM patients and 10 normal controls was obtained from Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) were identified by Limma package in R. The least absolute shrinkage and selection operator (LASSO) - support vector machines (SVM) were used to screen the hub genes. CIBERSORT algorithm was used to identify the different immune cells in distribution between T1DM and normal samples. Correlation of the hub genes and immune cells was analyzed by Spearman, and gene-GO-BP and gene-pathway interaction networks were constructed by Cytoscape plug-in ClueGO. Receiver operating characteristic (ROC) curves were used to assess diagnostic value of genes in T1DM. Results The 50 immune-related DEGs were obtained between the T1DM and normal samples. Then, the 50 immune-related DEGs were further screened to obtain the 5 hub genes. CIBERSORT analysis revealed that the distribution of plasma cells, resting mast cells, resting NK cells and neutrophils had significant difference between T1DM and normal samples. Natural cytotoxicity triggering receptor 3 (NCR3) was significantly related to the activated NK cells, M0 macrophages, monocytes, resting NK cells, and resting memory CD4+ T cells. Moreover, tumor necrosis factor (TNF) was significantly associated with naive B cell and naive CD4+ T cell. NCR3 [Area under curve (AUC) = 0.918] possessed a higher accuracy than TNF (AUC = 0.763) in diagnosis of T1DM. Conclusions The immune-related genes (NCR3 and TNF) and immune cells (NK cells) may play a vital regulatory role in the occurrence and development of T1DM, which possibly provide new ideas and potential targets for the immunotherapy of diabetes mellitus (DM).


Sign in / Sign up

Export Citation Format

Share Document