Dataset Information

Assessing immune infiltration and the tumor microenvironment for the diagnosis and prognosis of sarcoma.

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.

SUBMITTER: Zhu N

PROVIDER: S-EPMC7709254 | biostudies-literature | 2020 Dec

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Assessing immune infiltration and the tumor microenvironment for the diagnosis and prognosis of sarcoma.

Zhu Naiqiang N Hou Jingyi J

Cancer cell international 20201202 1

<h4>Background</h4>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.<h4>Method</h4>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 ...[more]

PMID: 33292275

Similar Datasets

Project description:Objective: Because of the modest immunotherapeutic response among ovarian carcinoma (OC) patients, it is significant to evaluate antitumor immune response and develop more effective precision immunotherapeutic regimens. Here, this study aimed to determine diverse immune subtypes of OC. Methods: This study curated the expression profiles of prognostic immunologically relevant genes and conducted consensus clustering analyses for determining immune subtypes among OC patients in TCGA cohort. With Boruta algorithm, characteristic genes were screened for conducting an immune scoring system through principal component analysis algorithm. The single-sample gene set enrichment analysis and ESTIAMTE methods were adopted for quantifying the immune infiltrates and responses to chemotherapeutic agents were estimated with pRRophetic algorithm. Two immunotherapeutic cohorts were used for investigating the efficacy of immune score in predicting therapeutic benefits. Results: Two immune subtypes were conducted among 377 OC patients. Immune subtype 2 was characterized by worse clinical prognosis, more frequent genetic variations and mutations, enhanced immune infiltrates, and increased expression of MHC molecules and programmed cell death protein 1/programmed death ligand 1 (PD-1/PD-L1). In total, 30 prognosis-relevant characteristic immune subtype-derived genes were identified for constructing the immune score of OC patients. High immune score was linked with more dismal prognosis, decreased immune infiltrations, and expression of MHC molecules. High immune score presented favorable sensitivity to doxorubicin and vinorelbine and reduced sensitivity to cisplatin. In addition, immune score possessed the potential in predicting benefits from anti-PD-1/PD-L1 therapy. Conclusion: Collectively, our findings propose two complex and diverse immune subtypes of OC. Quantitative assessment of immune subtypes in individual patients strengthens the understanding of tumor microenvironment features and promotes more effective immunotherapeutic regimens.

Project description:The neutrophils exhibit both anti-tumor and pro-tumor effects in cancers. The correlation between neutrophils and tumor development in lung adenocarcinoma (LUAD) is still uncertain, possibly due to a lack of specific neutrophil infiltration evaluation methods. In this study, we identified 30 hub genes that were significantly associated with neutrophil infiltration in LUAD through data mining, survival analysis, and multiple tumor-infiltrating immune cells (TICs) analysis, including TIMER, CIBERSORT, QUANTISEQ, XCELL, and MCPCOUNTER. Consensus clustering analysis showed that these 30 hub genes were correlated with clinical features in LUAD. We further developed a neutrophil scoring system based on these hub genes. The neutrophil score was significantly correlated with prognosis and tumor immune microenvironment (TIME) in LUAD. It was also positively associated with PD-L1 expression and negatively associated with tumor mutational burden (TMB). When combined with the neutrophil score, the predictive capacity of PD-L1 and TMB for prognosis was significantly improved. Thus, the 30 hub genes might play an essential role in the interaction of neutrophils and LUAD, and the neutrophil scoring system might effectually assess the infiltration of neutrophils. Furthermore, we verified the expression of these 30 genes in the LUAD tumor tissues collected from our department. We further found that overexpressed TNFAIP6 and TLR6 and downregulated P2RY13, SCARF1, DPEP2, PRAM1, CYP27A1, CFP, GPX3, and NCF1 in LUAD tissue might be potentially associated with neutrophils pro-tumor effects. The following in vitro experiments demonstrated that TNFAIP6 and TLR6 were significantly overexpressed, and P2RY13 and CYP27A1 were significantly downregulated in LUAD cell lines, compared to BEAS-2B cells. Knocking down TNFAIP6 in A549 and PC9 resulted in the upregulation of FAS, CCL3, and ICAM-1, and the downregulation of CCL2, CXCR4, and VEGF-A in neutrophils when co-culturing with the conditioned medium (CM) from LUAD cells. Knocking down TNFAIP6 in LUAD also led to an elevated early apoptosis rate of neutrophils. Therefore, overexpressed TNFAIP6 in LUAD cancer cells might lead to neutrophils "N2" polarization, which exhibited pro-tumor effects. Further research based on the genes identified in this pilot study might shed light on neutrophils' effects on LUAD in the future.

Project description:BackgroundChromobox family genes (CBXs) are known to play roles in numerous modifications of the chromatin in order to inhibit the transcription of target genes. CBXs have been shown to be expressed at high levels in many types of cancer and can also serve as a target gene for therapeutic purposes. However, little is known about the expression and prognostic value of CBXs in human sarcomas.MethodsThe transcription level of CBXs was analyzed using the Oncomine dataset, and the differential expression of CBXs in sarcoma was reported by the Gene Expression Profiling Interactive Analysis (GEPIA) dataset. We also used the CCLE dataset to evaluate the expression of CBXs in a sarcoma cell line. The prognostic value of CBXs was analyzed using GEPIA and Kaplan-Meier analysis. In addition, the corrections between CBXs and their co-expressed genes were reported using Oncomine and GEPIA datasets. DAVID was used to perform GO function enrichment analysis for the CBXs and their co-expression genes. Finally, TIMER was used to analyze the immune cell infiltration of CBXs in patients with sarcoma.ResultsHP1-α/β/γ (CBX1/3/5) and CBX4/6/8 were found to be overexpressed in human sarcoma, and CBXs were upregulated in almost all the sarcoma cell line. The expression levels of HP1-α/β/γ (CBX1/3/5) and CBX7 were associated with overall survival (OS) in patients with sarcoma, while high expression levels of CBX7 were related to disease-free survival (DFS). In addition, the expression levels of CBX2/6/7 were related to recurrence-free survival (RFS). We also found that the CBX family was positively correlated with the infiltration of immune cells, including CD8+ T cells, CD4+ T cells, B cells, macrophages, neutrophils, and dendritic cells, in sarcoma.ConclusionsThe results from the present study indicated that CBXs were significantly associated with prognosis and immunological status in sarcoma. These data suggest that CBXs could serve as potential biomarkers for prognosis and immune infiltration in human sarcoma.

Project description:Recent studies have highlighted that N6-methyladenosine (m6A) plays a significant role in tumorigenicity and progression. However, the mechanism of m6A modifications in the tumor microenvironment (TME) immune cell infiltration in cervical cancer (CC) remains unclear. Clinical and RNA sequencing data of 25 m6A RNA methylation regulators were acquired from the Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) database. LASSO Cox regression analysis was used to generate a prognostic risk signature. m6A modification patterns were identified based on the expression of 25 m6A regulators, and their correlation with TME immune cell-infiltrating characterization was analyzed. Principal component analysis was used to construct an m6A-scoring signature (m6A score) to evaluate the m6A modification patterns of individual CC samples and guide the selection of more effective immunotherapeutic strategies. Genetic and expression alterations of 25 m6A regulators were highly heterogeneous between CC and normal tissues. METTL14 and IGF2BP1 were selected to conduct the prognostic risk signature. Three m6A modification patterns were identified in 659 CC samples, which were associated with distinct clinical outcomes and biological pathways. The TME immune cell-infiltrating characterization of the three m6A modification patterns was highly consistent with 3 tumor immune phenotypes, including immune-excluded, immune-inflamed, and immune-desert phenotypes. Due to the heterogeneity of m6A modification patterns, an m6A scoring signature was established to evaluate the m6A modification patterns of individual CC samples. Univariate and multivariate Cox regression analyses revealed that the m6A score is a robust and independent prognostic biomarker for assessing the prognosis of CC patients. A low m6A score, characterized by higher somatic mutation and higher expression of proliferation-related and DNA repair-related genes, indicated poor overall survival. Activation of immune infiltration was exhibited by the high m6A score, which was likely to have a good response and clinical benefits to antiPD-1/L1 immunotherapy. This study highlights the prognostic value of 25 m6A regulators in CC. The m6A modification is related to immune regulation and the formation of TME heterogeneity and complexity. An m6A scoring signature to clarify the individual m6A modification pattern could enhance our understanding of TME immune cell-infiltrating characterization and guide immunotherapy strategies.

Project description:Lung adenocarcinoma (LUAD), a malignancy with high incidence and mortality rates worldwide, contains multiple genomic and epigenomic abnormalities. And the useful tumor markers associated with these abnormalities need further investigation. Whereas apoptosis is a form of programmed cell death, the expression of apoptosis-related genes in LUAD and its relationship with prognosis is unclear. In the present study, we identified 64 differentially expressed apoptosis-related genes (DEARGs) that were differentially expressed between LUAD tissue and normal lung tissue. Based on these DEARGs, all LUAD cases were classified into two subtypes using The Cancer Genome Atlas (TCGA) cohort to assess the prognostic value of apoptosis-related genes for survival. An 11-gene signature was established by applying the Least Absolute Shrinkage and Selection Operator (LASSO) Cox regression method to construct a multigene prediction model and classify all LUAD patients in the TCGA cohort into high or low AS-score groups. Patients in the low AS-score group had significantly higher survival and prognosis than those in the high AS-score group. Taking the median risk score of the AS-score, LUAD patients in the GSE68465 cohort were divided into two risk groups, low and high. The overall survival (OS) time was longer in the low AS-score group. Combined with clinical characteristics, the AS-score was an independent predictor of LUAD patients. Gene ontology (GO) and Kyoto Encylopedia of Genes and Genomes (KEGG) analyses showed that the differential genes between the two groups were mainly enriched in cellular immunity. Further analysis revealed higher immune checkpoint protein expression and higher tumor mutational burden (TMB) in the high AS-score group, suggesting better efficacy of immunotherapy in the high AS-score group than the low AS-score group. And the high AS-score group was better in chemotherapy and targeted therapy efficiency. In conclusion, the AS-score constructed based on apoptosis-related genes can predict the prognosis of LUAD patients and provide some guidance for the antitumor treatment of LUAD patients.

Project description:Background: Accumulating evidence shows that pyroptosis plays a crucial role in hepatocellular carcinoma (HCC). However, the relationship between pyroptosis-related long non-coding RNAs (lncRNAs) and HCC tumor characteristics remains enigmatic. We aimed to explore the predictive effect of pyroptosis-related lncRNAs (PRLs) in the prognosis of HCC. Methods: We comprehensively analyzed the role of the PRLs in the tumor microenvironment and HCC prognosis by integrating genomic data from patients of HCC. Consensus clustering analysis of PRLs was applied to identify HCC subtypes. A prognostic model was then established with a training cohort from The Cancer Genome Atlas (TCGA) using univariate and least absolute shrinkage and selection operator (LASSO) Cox regression analysis. Further, we evaluated the accuracy of this predictive model using a validation set. We predicted IC50s of commonly used chemotherapeutic and targeted drugs through the R package pRRophetic. Results: Based on pyroptosis-related lncRNAs, a prognostic risk signature composed of seven PRLs (MKLN1AS, AL031985.3, SNHG4, GHRLOS, AC005479.2, AC099850.4, and AC026412.3) was established. For long-term prognosis of HCC patients, our model shows excellent accuracy to forecast overall survival of HCC individuals both in training set and testing set. We found a significant correlation between clinical features and the risk score. Patients in the high-risk group had tumor characteristics associated with progression such as aggressive pathological grade and stage. Besides that, gene set enrichment analysis (GSEA) showed that cell cycle and focal adhesion were significantly enriched in the high-risk group. Conclusion: The association of the risk model constituted by these seven pyroptosis-related lncRNAs with clinical prognosis, tumor microenvironment, chemotherapy and small molecule drugs was evaluated. Our study provides strong evidence for individualized prediction of prognosis, shedding light on immunotherapy in HCC patients.

Dataset Information

Assessing immune infiltration and the tumor microenvironment for the diagnosis and prognosis of sarcoma.

Publications

Assessing immune infiltration and the tumor microenvironment for the diagnosis and prognosis of sarcoma.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets