Identification of Biomarkers Related to CD8+ T Cell Infiltration With Gene Co-expression Network in Lung Squamous Cell Carcinoma.
Ontology highlight
ABSTRACT:
SUBMITTER: Tang M
PROVIDER: S-EPMC8012856 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Clear cell renal cell carcinoma (ccRCC) is an extremely common kind of kidney cancer in adults. Immunotherapy and targeted therapy are particularly effective at treating ccRCC. In this study, weighted gene co-expression network analysis and a deconvolution algorithm that quantifies the cellular composition of immune cells were used to analyze ccRCC expression data from the Gene Expression Omnibus database, and identify modules related to CD8+ T cells. Ten hub genes (LCK, CD2, CD3D, CD3G, IRF1, IFNG, CCR5, CD8A, CCL5, and CXCL9) were identified by co-expression network and protein-protein interactions network analysis. Datasets obtained from The Cancer Genome Atlas were analyzed and the data revealed that the hub genes were meaningfully up-regulated in tumor tissues and correlated with promotion of tumor progression. After Kaplan-Meier analysis and Oncomine meta-analysis, CCL5 was selected as a prognostic biomarker. Finally, the experimental results show that reduced expression of CCL5 decreased cell proliferation and invasion in the ccRCC cell line. Various analyses were performed and verified, CCL5 is a potential biomarker and therapeutic target which related to CD8+ T cell infiltration in ccRCC.
| S-EPMC7066925 | biostudies-literature
Project description:The five-year survival rate of lung squamous cell carcinoma is significantly lower than that of other cancer types. It is therefore urgent to discover novel prognosis biomarkers and therapeutic targets and understand their correction with infiltrating immune cells to improve the prognosis of patients with lung squamous cell carcinoma. In this study, we employed robust rank aggregation algorithms to overcome the shortcomings of small sizes and potential bias in each Gene Expression Omnibus dataset of lung squamous cell carcinoma and identified 513 robust differentially expressed genes including 220 upregulated and 293 downregulated genes from six microarray datasets. Functional enrichment analysis showed that these robust differentially expressed genes were obviously involved in the extracellular matrix and structure organization, epidermis development, cell adhesion molecule binding, p53 signaling pathway, and interleukin-17 signaling pathway to affect the progress of lung squamous cell carcinoma. We further identified six hub genes from 513 robust differentially expressed genes by protein-protein interaction network and 10 topological analyses. Moreover, the results of immune cell infiltration analysis from six integrated Gene Expression Omnibus datasets and our sequencing transcriptome data demonstrated that the abundance of monocytes was significantly lower in lung squamous cell carcinoma compared to controls. Immune correlation analysis and survival analysis of hub genes suggested that three hub genes, collagen alpha-1(VII) chain, mesothelin, and chordin-like protein 1, significantly correlated with tumor-infiltrating monocytes as well as may be potential prognostic biomarkers and therapy targets in lung squamous cell carcinoma. The investigation of the correlation of hub gene markers and infiltrating monocytes can also improve to well understand the molecular mechanisms of lung squamous cell carcinoma development.
| S-EPMC9388304 | biostudies-literature
Project description:BackgroundAlthough immune checkpoint blockade (ICB) therapy has brought survival benefits to patients with specific cancer types, most of cancer patients remain refractory to the ICB therapy, which is largely attributed to the immunosuppressive tumor microenvironment. Thereby, it is urgent to profile key molecules and signal pathways responsible for modification of tumor microenvironment.MethodsMultiple databases of esophageal squamous cell carcinoma (ESCC) were integratively analyzed to screen candidate genes responsible for infiltration of CD8+ T cells. Expression of pescadillo ribosomal biogenesis factor 1 (PES1) in clinical ESCC samples was examined by qRT-PCR, western blotting, and immunohistochemistry. The mechanisms of PES1 were investigated via RNA sequencing and mass spectrometry followed by immunoprecipitation and proximity ligation assay. The clinical and therapeutic significance of PES1 in ESCC was comprehensively investigated using ESCC cells and mouse model.ResultsPES1 was significantly upregulated and correlated with poor prognosis in ESCC patients. PES1 knockdown decreased ESCC cell growth in vitro and in vivo and enhanced the efficacy of ICB therapy in mouse model, which was established through subcutaneous inoculation with ESCC cells. Analyses on RNA sequencing and mass spectrometry suggested that PES1 expression was negatively correlated with IL15 and ILF3 was one of the PES1-associated proteins. It has been known that ILF3 interacts with and stabilizes IL15 mRNA to increase IL15 protein level. Our data further indicated that PES1 interfered with the interaction between ILF3 and IL15 mRNA and impaired ILF3-mediated stabilization of IL15 mRNA, which eventually reduced the protein level of IL15. Interestingly, the inhibitory effect of ICB therapy boosted by PES1 knockdown dramatically antagonized by knockdown of IL15, which suppressed the tumor-infiltrated CD8+ T cells in ESCC. Finally, we confirmed the relationships among PES1, IL15, and CD8+ T cell infiltration in 10 locally advanced ESCC patients receiving ICB neoadjuvant therapy and demonstrated that ICB therapy would be more effective in those with low expression of PES1.ConclusionsAltogether, our findings herein provided novel insights on biological function and clinical significance of PES1 and suggested that high expression of PES1 could suppress ILF3-IL15 axis-mediated immunosurveillance and promote resistance to ICB through restraining tumor-infiltrated CD8+ T cells.
| S-EPMC10037800 | biostudies-literature
Project description:Lung squamous cell carcinoma (LUSC) is associated with a worse prognosis than other histological subtypes of non-small cell lung cancer. Due to the vital role of CD8+ T cells in anti-tumor immunity, the characterization of CD8+ T cell infiltration-related (CTLIR) gene signature in LUSC is worthy of in-depth exploration. In our study, tumor tissues of LUSC patients from Renmin Hospital of Wuhan University were stained by multiplex immunohistochemistry to evaluate the density of infiltrated CD8+ T cells and explore the correlation with immunotherapy response. We found that the proportion of LUSC patients who responded to immunotherapy was higher in the high density of CD8+ T cell infiltration group than in the low density of CD8+ T cell infiltration group. Subsequently, we collected bulk RNA-sequencing data from The Cancer Genome Atlas (TCGA) database. The abundance of infiltrating immune cells in LUSC patients was analyzed by using CIBERSORT algorithm, and weighted correlation network analysis was performed to identify the co-expressed gene modules related to CD8+ T cells. We then developed a prognostic gene signature based on CD8+ T cell co-expressed genes and calculated the CTLIR risk score, which stratified LUSC patients into high-risk and low-risk groups. With univariate and multivariate analyses, the gene signature was identified as an independent prognostic factor in LUSC patients. The overall survival of LUSC patients in the high-risk group was significantly shorter than that of the low-risk group in the TCGA cohort, which was validated in Gene Expression Omnibus datasets. We analyzed immune cell infiltration in the tumor microenviroment and found fewer CD8+ T cells and more regulatory T cell infiltration in the high-risk group, which is characterized as an immunosuppressive phenotype. Furthermore, the LUSC patients in the high-risk group were predicted to have a better response to immunotherapy than those in the low-risk group when treated with PD-1 and CTLA4 inhibitors. In conclusion, we performed a comprehensive molecular analysis of the CTLIR gene signature in LUSC and constructed a risk model for LUSC patients to predict prognosis and immunotherapy response.
| S-EPMC10246359 | biostudies-literature
Project description:Intrinsic cancer cells and the tumor-infiltrating immune cells (TIICs) recruited to the immune microenvironment define the malignant phenotype of lung squamous cell carcinoma (LUSC). Understanding more about the immune microenvironment of LUSC enables the selection of high-risk patients who would derive benefit from immunotherapy. Based on large public LUSC cohorts obtained from TCGA and GEO datasets, 22 types of infiltrating immune cell subgroups were evaluated by CIBERSORT. Meta-analysis, principal component analysis (PCA), single-sample gene set enrichment analysis (ssGSEA), and hierarchical clustering analysis were used to evaluate specific immune responses of LUSC. The distribution of TIICs of LUSC was entirely different from normal. TIIC subpopulations were also found to be closely associated with clinical features and molecular subtypes. Unsupervised clustering analysis revealed that three distinct TIIC subgroups existed with different survival patterns. TIICs are extensively implicated in the pathogenesis and development of LUSC. Characterizing the composition of TIICs influences the metabolism, pathological stage, and survival of tumor patients. It is hoped that this immune landscape could provide a more accurate understanding of the development and immunotherapy of LUSC.
| S-EPMC7569448 | biostudies-literature
Project description:Hepatocellular carcinoma (HCC) is a common cancer with poor prognosis. Due to the lack of effective biomarkers and its complex immune microenvironment, the effects of current HCC therapies are not ideal. In this study, we used the GSE57957 microarray data from Gene Expression Omnibus database to construct a co-expression network. The weighted gene co-expression network analysis and CIBERSORT algorithm, which quantifies cellular composition of immune cells, were used to identify modules related to immune cells. Four hub genes (EFTUD2, GAPDH, NOP56, PA2G4) were identified by co-expression network and protein-protein interactions network analysis. We examined these genes in TCGA database, and found that the four hub genes were highly expressed in tumor tissues in multiple HCC groups, and the expression levels were significantly correlated with patient survival time, pathological stage and tumor progression. On the other hand, methylation analysis showed that the up-regulation of EFTUD2, GAPDH, NOP56 might be due to the hypomethylation status of their promoters. Next, we investigated the correlations between the expression levels of four hub genes and tumor immune infiltration using Tumor Immune Estimation Resource (TIMER). Gene set variation analysis suggested that the four hub genes were associated with numerous pathways that affect tumor progression or immune microenvironment. Overall, our results showed that the four hub genes were closely related to tumor prognosis, and may serve as targets for treatment and diagnosis of HCC. In addition, the associations between these genes and immune infiltration enhanced our understanding of tumor immune environment and provided new directions for the development of drugs and the monitoring of tumor immune status.
| S-EPMC8262220 | biostudies-literature
Project description:BackgroundHead and neck squamous cell carcinoma (HNSCC) is a frequently lethal malignancy, and the mortality is considerably high. The tumor microenvironment (TME) has been identified as a critical participation in cancer development, treatment, and prognosis. However, competing endogenous RNA (ceRNA) networks grouping with immune/stromal scores of HNSCC patients need to be further illustrated. Therefore, our study aimed to provide clues for searching promising prognostic markers of TME in HNSCC.Materials and methodsESTIMATE algorithm was used to calculate immune scores and stromal scores of the enrolled HNSCC patients. Differentially expressed genes (DEGs), lncRNAs (DELs), and miRNAs (DEMs) were identified by comparing the expression difference between high and low immune/stromal scores. Then, a ceRNA network and protein-protein interaction (PPI) network were constructed for selecting hub regulators. In addition, survival analysis was performed to access the association between immune scores, stromal scores, and differentially expressed RNAs in the ceRNA network and the overall survival (OS) of HNSCC patients. Then, the GSE65858 datasets from Gene Expression Omnibus (GEO) database was used for verification. At last, the difference between the clinical characteristics and immune cell infiltration in different expression groups of IL10RA, PRF1, and IL2RA was analyzed.ResultsSurvival analysis showed a better OS in the high immune score group, and then we constructed a ceRNA network composed of 97 DEGs, 79 DELs and 22 DEMs. Within the ceRNA network, FOXP3, IL10RA, STAT5A, PRF1, IL2RA, miR-148a-3p, miR-3065-3p, and lncRNAs, including CXCR2P1, HNRNPA1P21, CTA-384D8.36, and IGHV1OR15-2, were closely correlated with the OS of HNSCC patients. Especially, using the data from GSE65858, we successfully verified that IL10RA, PRF1, and IL2RA were not only significantly upregulated in patients high immune scores, but also their high expressions were associated with longer survival time. In addition, stratified analysis showed that PRF1 and IL2RA might be involved in the mechanism of tumor progress.ConclusionIn conclusion, we constructed a ceRNA network related to the TME of HNSCC, which provides candidates for therapeutic intervention and prognosis evaluation.
| S-EPMC9463596 | biostudies-literature
Project description:Cytotoxic T cells expressing cell surface CD8 played a key role in anti-cancer immunotherapy, including kidney renal clear cell carcinoma (KIRC). Here we set out to comprehensively analyze and evaluate the significance of CD8+ T cell-related markers for patients with KIRC. We checked immune cell response in KIRC and identified cell type-specific markers and related pathways in the tumor-infiltrating CD8+ T (TIL-CD8T) cells. We used these markers to explore their prognostic signatures in TIL-CD8+ T by evaluating their prognostic efficacy and group differences at various levels. Through pan-cancer analysis, 12 of 63 up-regulated and 162 of 396 down-regulated genes in CD8+ T cells were found to be significantly correlated with the survival prognosis. Based on our highly integrated multi-platform analyses across multiple datasets, we constructed a 6-gene risk scoring model specific to TIL-CD8T. In this model, high TIL-CD8 sig score was corresponding to a higher incidence frequency of copy number variation and drug sensitivity to sorafenib. Moreover, the prognosis of patients with the same or similar immune checkpoint gene levels could be distinguished from each other by TIL-CD8 sig score.
| S-EPMC8544304 | biostudies-literature
Project description:Background: CD8+ T cells, a critical component of the tumor immune microenvironment, have become a key target of cancer immunotherapy. Considering the deficiency of robust biomarkers for head and neck squamous cell carcinoma (HNSCC), this study aimed at establishing a molecular signature associated with CD8+T cells infiltration. Methods: Single-cell RNA sequencing data retrieved from the Gene Expression Omnibus (GEO) database was analyzed to obtain the different cell types. Next, the cell proportions were investigated through deconvolution of RNA sequencing in the Cancer Genome Atlas (TCGA) database, and then the immune-related genes (IRGs) were identified by weighted gene co-expression network analysis (WGCNA). LASSO-Cox analysis was employed to establish a gene signature, followed by validation using a GEO dataset. Finally, the molecular and immunological properties, and drug responses between two subgroups were explored by applying "CIBERSORT", "ESTIMATE", and single sample gene set enrichment analysis (ssGSEA) methods. Results: A total of 215 differentially expressed IRGs were identified, of which 45 were associated with the overall survival of HNSCC. A risk model was then established based on eight genes, including DEFB1, AICDA, TYK2, CCR7, SCARB1, ULBP2, STC2, and LGR5. The low-risk group presented higher infiltration of memory activated CD4+ T cells, CD8+ T cells, and plasma cells, as well as a higher immune score, suggesting that they could benefit more from immunotherapy. On the other hand, the high-risk group showed higher abundance of activated mast cells and M2 macrophages, as well as a lower immune score. Conclusion: It was evident that the 8-gene signature could accurately predict HNSCC prognosis and thus it may serve as an index for clinical treatment.
| S-EPMC9355512 | biostudies-literature
Project description:Rheumatoid arthritis (RA) is an autoimmune rheumatic disease, which do not respond well to current treatment partially. Therefore, further in-depth elucidation of the molecular mechanism and pathogenesis of RA is urgently needed for the diagnosis, personalized therapy and drug development. Herein, we collected 111 RA samples from Gene Expression Omnibus (GEO) database, and conducted differentially expressed genes and GESA analysis. Abnormal activation and imbalance of immune cells in RA were observed. WGCNA was utilized to explore the gene modules and CD8+ T cell-related genes (CRGs) were chosen for KEGG and GO analysis. Besides, to explore biomarkers of RA in depth, machine learning algorithms and bioinformatics analysis were used, and we identified GDF15, IGLC1, and IGHM as diagnostic markers of RA, which was confirmed by clinical samples. Next, ssGSEA algorithms were adopted to investigate the differences in immune infiltration of 23 immune cell subsets between RA and healthy control group. Finally, optimal classification analysis based on consensus clustering combined with ssGSEA algorithms were conducted. GDF15 was revealed that to be positively correlated with mast cells and type 2 T helper cells, but negatively correlated with most other immune cells. On the other hand, IGHM and IGLC1 were negatively correlated with CD56dim natural killer cells, while positively associated with other immune cells. Finally, RA samples in subtype A exhibited a higher immune infiltration status. This study could provide guidance for individualized treatment of RA patients and provide new targets for drug design.
| S-EPMC10866417 | biostudies-literature