Explore the vast range of titles available.

Tumour-infiltrating immune cells are a source of important prognostic information for patients with resectable colon cancer. We developed a novel immune model based on systematic assessments of the immune landscape inferred from bulk tumor transcriptomes of stage I–III colon cancer patients. The “Cell type Identification By Estimating Relative Subsets Of RNA Transcripts (CIBERSORT)” algorithm was used to estimate the fraction of 22 immune cell types from six microarray public datasets. The random forest method and least absolute shrinkage and selection operator model were then used to establish immunoscores for diagnosis and prognosis. By comparing immune cell compositions in samples of 870 colon cancer patients and 70 normal controls, we constructed a diagnostic model, designated the diagnostic immune risk score (dIRS), that showed high specificity and sensitivity in both the training [area under the curve (AUC) = 0.98, p < 0.001] and validation (AUC 0.96, p < 0.001) sets. We also established a prognostic immune risk score (pIRS) that was found to be an independent prognostic factor for relapse-free survival in every series (training: HR 2.23; validation: HR 1.65; entire: HR 2.01; p < 0.001 for all), which showed better prognostic value than TNM stage. In addition, integration of the pIRS with clinical characteristics in a composite nomogram showed improved accuracy of relapse risk prediction, providing a higher net benefit than TNM stage, with well-fitted calibration curves. The proposed dIRS and pIRS models represent promising novel signatures for the diagnosis and prognosis prediction of colon cancer.

Immune cells have essential auxiliary functions and influence clinical outcomes in cancer, with high immune infiltration being associated with improved clinical outcomes and better response to treatment in breast cancer (BC). However, studies to date have not fully considered the tumor-infiltrating immune cell (TIIC) landscape in tumors. This study investigated potential biomarkers based on TIICs to improve prognosis and treatment effect in BC. We enrolled 5112 patients for analysis and used cell type identification by estimating relative subsets of RNA transcripts (CIBERSORT), a new computational algorithm, to quantify 22 TIICs in primary BC. From the results of univariate Cox regression, 12 immune cells were determined to be significantly related to the overall survival (OS) of BC patients. Furthermore, least absolute shrinkage and selection operator (LASSO) and multivariate Cox regression analyses were applied to construct an immune prognostic model based on six potential biomarkers. By dividing patients into low- and high-risk groups, a significant distinction in OS was found in the training cohort, with 20-year survival rates of 42.6% and 26.3%, respectively. Applying a similar protocol to validation and test cohorts, we found that OS was significantly shorter in the high-risk group than in the low-risk group, regardless of the molecular subtype of BC. Using the immune score model to predict the effect of BC patients to chemotherapy, the survival advantage for the low-risk group was evident among those who received chemotherapy, regardless of the chemotherapy regimen. In evaluating the predictive value of the nomogram, a decision curve showed better predictive accuracy than the standard tumor-node-metastasis (TNM) staging system. The immune cell infiltration-based immune score model can be effectively and efficiently used to predict the prognosis of BC patients as well as the effect of chemotherapy.

Tumor‐infiltrating immune cells play a crucial role in tumor progression and response to treatment. However, the limited studies on infiltrating immune cells have shown inconsistent and even controversial results for osteosarcoma (OS). In addition, the dynamic changes of infiltrating immune cells after neoadjuvant chemotherapy are largely unknown. We downloaded the RNA expression matrix and clinical information of 80 OS patients from the TARGET database. CIBERSORT was used to evaluate the proportion of 22 immune cell types in patients based on gene expression data. M2 macrophages were found to be the most abundant immune cell type and were associated with improved survival in OS. Another cohort of pretreated OS samples was evaluated by immunohistochemistry to validate the results from CIBERSORT analysis. Matched biopsy and surgical samples from 27 patients were collected to investigate the dynamic change of immune cells and factors before and after neoadjuvant chemotherapy. Neoadjuvant chemotherapy was associated with increased densities of CD3+ T cells, CD8+ T cells, Ki67 + CD8+ T cells and PD‐L1+ immune cells. Moreover, HLA‐DR‐CD33+ myeloid‐derived suppressive cells (MDSC) were decreased after treatment. We determined that the application of chemotherapy may activate the local immune status and convert OS into an immune “hot” tumor. These findings provide rationale for investigating the schedule of immunotherapy treatment in OS patients in future clinical trials. Host anti–tumor immune response boosted by neoadjuvant chemotherapy. Following neoadjuvant chemotherapy, CD3+ T cells increased significantly and there was a trend of increased cytotoxic T cells. CD8+ T cells in both tumor center and stroma also increased remarkably. Importantly, activated CD8+ T cells, defined as Ki67 + CD8+ T cells, were more abundant in post–chemotherapy samples, and were negatively correlated with the proliferation ability of tumor cells.

This study aimed to find potential diagnostic markers for osteoarthritis (OA) and analyze the role of immune cells infiltration in this pathology. We used OA datasets from the Gene Expression Omnibus database. First, R software was used to identify differentially expressed genes (DEGs) and perform functional correlation analysis. Then least absolute shrinkage and selection operator (LASSO) logistic regression and support vector machine-recursive feature elimination algorithms were used to screen and verify the diagnostic markers of OA. Finally, CIBERSORT was used to evaluate the infiltration of immune cells in OA tissues, and the correlation between diagnostic markers and infiltrating immune cells was analyzed. A total of 458 DEGs were screened in this study. GRB10 and E2F3 (AUC = 0.962) were identified as diagnostic markers of OA. Immune cell infiltration analysis found that resting mast cells, T regulatory cells, CD4 memory resting T cells, activated NK cells, and eosinophils may be involved in the OA process. In addition, GRB10 was correlated with NK resting cells, naive CD4 + T cells, and M1 macrophages, while E2F3 was correlated with resting mast cells. In conclusion, GRB10 and E2F3 can be used as diagnostic markers of osteoarthritis, and immune cell infiltration plays an important role in the occurrence and progression of OA.

Background Tongue squamous cell carcinoma (TSCC) is characterized by aggressive invasion and poor prognosis. Currently, immune checkpoint inhibitors may prolong overall survival compared with conventional treatments. However, PD1/PDL1 remain inapplicable in predicting the prognosis of TSCC; thus, it is urgent to explore the genetic characteristics of TSCC. Materials and methods We utilized single‐sample gene set enrichment analysis (ssGSEA) to classify TSCC patients from the TCGA database into clusters with different immune cell infiltrations. ESTIMATE (immune‐related scores) and CIBERSORT (immune cell distribution) analyses were used to evaluate the immune landscape among clusters. GO, KEGG, and GSEA analyses were performed to analyze the different underlying molecular mechanisms in the clusters. Based on the immune characteristics, we applied the LASSO Cox regression to select hub genes and construct a prognostic risk model. Finally, we established an interactive network among these hub genes by using Cytoscape, and a pan‐cancer analysis to further verify and decipher the innate function of these genes. Results Using ssGSEA, we constructed three functional clusters with different overall survival and immune‐cell infiltration. ESTIMATE and CIBERSORT analyses revealed the different distributions of immune cells (T cells, B cells, and macrophages) with diverse immune‐related scores (ESTIMATE, immune, stromal, and tumor purity scores). Moreover, pathways including those of the interferon‐gamma response, hypoxia, and glycolysis of the different subtypes were investigated to elucidate their involvement in mediating the heterogeneous immune characteristics. Subsequently, after LASSO Cox regression, a signature of 15 immune‐related genes was established that is more prognostically effective than the TNM stage. Furthermore, three hub genes—PGK1, GPI, and RPE—were selected using Cytoscape evaluation and verified by immunohistochemistry. PGK1, the foremost regulator, was a comprehensively profiled pan‐cancer, and a PGK1‐based interactive network was established. Conclusion Our results suggest that immune‐related genes and clusters in TSCC have the potential to guide individualized treatments. Our findings suggested that these immune‐related genes may act as promising indicators for mediating the immune microenvironment characterizations and providing novel insight into immunotherapy.

Background Ovarian cancer (OC) is the most malignant tumor in the female reproductive system. About 75% of OC in complete remission of clinical symptoms still develop a recurrence. Therefore, searching for new treatment methods plays an important role in improving the prognosis of OC. Methods We downloaded the MAF files, RNA-seq data and clinical information from the TCGA database. The “maftools” package in R software was used to visualize the OC mutation data. We calculated the tumor mutation burden (TMB) of OC and analyzed its correlation with clinicopathological parameters and prognostic value. Tumor mutation burden related signature model was constructed to predict the overall survival (OS) of OC. Results The results revealed that there was a statistical correlation between TMB and FIGO stage, grade and tumor residual size of ovarian cancer patients. The Kaplan–Meier curve indicated that a high TMB is associated with better clinical outcomes of OC. The difference analysis indicated 24 upregulated genes and 619 downregulated genes in the high-TMB group compared with the low-TMB group. Besides, the TMBRS model based on five hub genes (RBMS3, PLA2G5, CDH2, AMHR2 and ADAMTS8) was constructed to predict the OS of OC. The ROC curve and validation data sets all revealed that the TMBRS model was reliable in predicting recurrence risk. Immune microenvironment analysis indicated the correlations between TMB and infiltrating immune cells. Conclusions Our results suggest that TMB plays an important role in the prognosis and guiding immunotherapy of OC. By detecting the TMB of OC, clinicians can more accurately treat patients with immunotherapy, thereby improving their survival rate.

Annexin A1 (ANXA1) is a phospholipid-linked protein involved in inflammation, immune response, and mast cell reactivity. Recently, we reported that ANXA1 is associated with aggressive features of triple-negative breast cancer (TNBC); however, its clinical relevance remains controversial. We hypothesized that human TNBC with high expression of ANXA1 mRNA is associated with pro-cancerous immune cell infiltration, including mast cells, and with an aggressive phenotype. Clinical and RNA-seq data were obtained from The Cancer Genome Atlas (TCGA, n = 1079) and Molecular Taxonomy of Breast Cancer International Consortium (METABRIC) (n = 1904). TNBC patients had significantly higher levels of ANXA1 expression compared to the other subtypes in both TCGA and METABRIC cohorts (p < 0.001). ANXA1 protein expression was assessed by immunohistochemistry in Japanese TNBC patient cohort (n = 48), where 17 cases (35.4%) had positive ANXA1 staining, and their overall survival was significantly shorter compared with negative staining group (p = 0.008). The CIBERSORT algorithm was used to calculate immune cell infiltrations. ANXA1 high tumors were associated with activated mast cells and M2 macrophages (p > 0.01), but did not show any association with tumor heterogeneity nor cytolytic activity. High expression of ANXA1 group enriched inflammation, epithelial-to-mesenchymal transition (EMT), and angiogenesis-related genes in a gene set enrichment assay in both cohorts. To our knowledge, this is the first study to demonstrate that ANXA1 is associated with infiltration of mast cells and inflammation that is associated with the aggressive phenotype of TNBC, such as EMT and angiogenesis.

Background Kawasaki disease (KD) is a febrile systemic vasculitis involvingchildren younger than five years old. However, the specific biomarkers and precise mechanisms of this disease are not fully understood, which can delay the best treatment time, hence, this study aimed to detect the potential biomarkers and pathophysiological process of KD through bioinformatic analysis. Methods The Gene Expression Omnibus database (GEO) was the source of the RNA sequencing data from KD patients. Differential expressed genes (DEGs) were screened between KD patients and healthy controls (HCs) with the “ limma ” R package. Weighted gene correlation network analysis (WGCNA) was performed to discover the most corresponding module and hub genes of KD. The node genes were obtained by the combination of the least absolute shrinkage and selection operator (LASSO) regression model with the top 5 genes from five algorithms in CytoHubba, which were further validated with the receiver operating characteristic curve (ROC curve). CIBERSORTx was employed to discover the constitution of immune cells in KDs and HCs. Functional enrichment analysis was performed to understand the biological implications of the modular genes. Finally, competing endogenous RNAs (ceRNA) networks of node genes were predicted using online databases. Results A total of 267 DEGs were analyzed between 153 KD patients and 92 HCs in the training set, spanning two modules according to WGCNA. The turquoise module was identified as the hub module, which was mainly enriched in cell activation involved in immune response, myeloid leukocyte activation, myeloid leukocyte mediated immunity, secretion and leukocyte mediated immunity biological processes; included type II diabetes mellitus, nicotinate and nicotinamide metabolism, O-glycan biosynthesis, glycerolipid and glutathione metabolism pathways. The node genes included ADM, ALPL, HK3, MMP9 and S100A12, and there was good performance in the validation studies. Immune cell infiltration analysis revealed that gamma delta T cells, monocytes, M0 macrophage, activated dendritic cells, activated mast cells and neutrophils were elevated in KD patients. Regarding the ceRNA networks, three intact networks were constructed: NEAT1/NORAD/XIST-hsa-miR-524-5p-ADM, NEAT1/NORAD/XIST-hsa-miR-204-5p-ALPL, NEAT1/NORAD/XIST-hsa-miR-524-5p/hsa-miR-204-5p-MMP9. Conclusion To conclude, the five-gene signature and three ceRNA networks constructed in our study are of great value in the early diagnosis of KD and might help to elucidate our understanding of KD at the RNA regulatory level.

BackgroundSystemic lupus erythematosus (SLE) is an autoimmune illness caused by a malfunctioning immunomodulatory system. China has the second highest prevalence of SLE in the world, from 0.03% to 0.07%. SLE is diagnosed using a combination of immunological markers, clinical symptoms, and even invasive biopsy. As a result, genetic diagnostic biomarkers for SLE diagnosis are desperately needed.MethodFrom the Gene Expression Omnibus (GEO) database, we downloaded three array data sets of SLE patients’ and healthy people’s peripheral blood mononuclear cells (PBMC) (GSE65391, GSE121239 and GSE61635) as the discovery metadata (nSLE = 1315, nnormal = 122), and pooled four data sets (GSE4588, GSE50772, GSE99967, and GSE24706) as the validate data set (nSLE = 146, nnormal = 76). We screened the differentially expressed genes (DEGs) between the SLE and control samples, and employed the least absolute shrinkage and selection operator (LASSO) regression, and support vector machine recursive feature elimination (SVM-RFE) analyze to discover possible diagnostic biomarkers. The candidate markers’ diagnostic efficacy was assessed using the receiver operating characteristic (ROC) curve. The reverse transcription quantitative polymerase chain reaction (RT-qPCR) was utilized to confirm the expression of the putative biomarkers using our own Chinese cohort (nSLE = 13, nnormal = 10). Finally, the proportion of 22 immune cells in SLE patients was determined using the CIBERSORT algorithm, and the correlations between the biomarkers’ expression and immune cell ratios were also investigated.ResultsWe obtained a total of 284 DEGs and uncovered that they were largely involved in several immune relevant pathways, such as type І interferon signaling pathway, defense response to virus, and inflammatory response. Following that, six candidate diagnostic biomarkers for SLE were selected, namely ABCB1, EIF2AK2, HERC6, ID3, IFI27, and PLSCR1, whose expression levels were validated by the discovery and validation cohort data sets. As a signature, the area under curve (AUC) values of these six genes reached to 0.96 and 0.913, respectively, in the discovery and validation data sets. After that, we checked to see if the expression of ABCB1, IFI27, and PLSCR1 in our own Chinese cohort matched that of the discovery and validation sets. Subsequently, we revealed the potentially disturbed immune cell types in SLE patients using the CIBERSORT analysis, and uncovered the most relevant immune cells with the expression of ABCB1, IFI27, and PLSCR1.ConclusionOur study identified ABCB1, IFI27, and PLSCR1 as potential diagnostic genes for Chinese SLE patients, and uncovered their most relevant immune cells. The findings in this paper provide possible biomarkers for diagnosing Chinese SLE patients.

Sepsis is the major cause of mortality in the intensive care unit. The aim of this study was to identify the key prognostic biomarkers of abnormal expression and immune infiltration in sepsis. In this study, a total of 36 differentially expressed genes were identified to be mainly involved in a number of immune-related Gene Ontology terms and Kyoto Encyclopedia of Genes and Genomes pathways. The hub genes (MMP9 and C3AR1) were significantly related to the prognosis of sepsis patients. The immune infiltration analysis indicated a significant difference in the relative cell content of naive B cells, follicular Th cells, activated NK cells, eosinophils, neutrophils and monocytes between sepsis and normal controls. Weighted gene co-expression network analysis and a de-convolution algorithm that quantifies the cellular composition of immune cells were used to analyse the sepsis expression data from the Gene Expression Omnibus database and to identify modules related to differential immune cells. CEBPB is the key immune-related gene that may be involved in sepsis. Gene set enrichment analysis revealed that CEBPB is involved in the processes of T cell selection, B cell–mediated immunity, NK cell activation and pathways of T cells, B cells and NK cells. Therefore, CEBPB may play a key role in the biological and immunological processes of sepsis.