Dataset Information

Cancer Stemness-Based Prognostic Immune-Related Gene Signatures in Lung Adenocarcinoma and Lung Squamous Cell Carcinoma

ABSTRACT:

Background

Cancer stem cells (CSCs) refer to cells with self-renewal capability in tumors. CSCs play important roles in proliferation, metastasis, recurrence, and tumor heterogeneity. This study aimed to identify immune-related gene-prognostic models based on stemness index (mRNAsi) in lung adenocarcinoma (LUAD) and lung squamous cell carcinoma (LUSC), respectively.

Methods

X-tile software was used to determine the best cutoff value of survival data in LUAD and LUSC based on mRNAsi. Tumor purity and the scores of infiltrating stromal and immune cells in lung cancer tissues were predicted with ESTIMATE R package. Differentially expressed immune-related genes (DEIRGs) between higher- and lower-mRNAsi subtypes were used to construct prognostic models.

Results

mRNAsi was negatively associated with StromalScore, ImmuneScore, and ESTIMATEScore, and was positively associated with tumor purity. LUAD and LUSC samples were divided into higher- and lower-mRNAsi groups with X-title software. The distribution of immune cells was significantly different between higher- and lower-mRNAsi groups in LUAD and LUSC. DEIRGs between those two groups in LUAD and LUSC were enriched in multiple cancer- or immune-related pathways. The network between transcriptional factors (TFs) and DEIRGs revealed potential mechanisms of DEIRGs in LUAD and LUSC. The eight-gene-signature prognostic model (ANGPTL5, CD1B, CD1E, CNTFR, CTSG, EDN3, IL12B, and IL2)-based high- and low-risk groups were significantly related to overall survival (OS), tumor microenvironment (TME) immune cells, and clinical characteristics in LUAD. The five-gene-signature prognostic model (CCL1, KLRC3, KLRC4, CCL23, and KLRC1)-based high- and low-risk groups were significantly related to OS, TME immune cells, and clinical characteristics in LUSC. These two prognostic models were tested as good ones with principal components analysis (PCA) and univariate and multivariate analyses. Tumor T stage, pathological stage, or metastasis status were significantly correlated with DEIRGs contained in prognostic models of LUAD and LUSC.

Conclusion

Cancer stemness was not only an important biological process in cancer progression but also might affect TME immune cell infiltration in LUAD and LUSC. The mRNAsi-related immune genes could be potential biomarkers of LUAD and LUSC. Evaluation of integrative characterization of multiple immune-related genes and pathways could help to understand the association between cancer stemness and tumor microenvironment in lung cancer.

SUBMITTER: Li N

PROVIDER: S-EPMC8567176 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:Accumulating evidence indicates that immunotherapy helped to improve the survival and quality-of-life of patients with lung adenocarcinoma (LUAD) or lung squamous cell carcinoma (LUSC) besides chemotherapy and gene targeting treatment. This study aimed to develop immune-related gene signatures in LUAD and LUSC subtypes, respectively. LUAD and LUSC samples were divided into high- and low-abundance groups of immune cell infiltration (Immunity_H and Immunity_L) based on the abundance of immune cell infiltrations. The distribution of immune cells was significantly different between the high- and low-immunity subtypes in LUAD and LUSC samples. The differentially expressed genes (DEGs) between those two groups in LUAD and LUSC contain some key immune-related genes, such as PDL1, PD1, CTLA-4, and HLA families. The DEGs were enriched in multiple immune-related pathways. Furthermore, the seven-immune-related-gene-signature (CD1B, CHRNA6, CLEC12B, CLEC17A, CLNK, INHA, and SLC14A2) prognostic model-based high- and low-risk groups were significantly associated with LUAD overall survival and clinical characteristics. The eight-immune-related-gene-signature (C4BPB, FCAMR, GRAPL, MAP1LC3C, MGC2889, TRIM55, UGT1A1, and VIPR2) prognostic model-based high- and low-risk groups were significantly associated with LUSC overall survival and clinical characteristics. The prognostic models were tested as good ones by receiver operating characteristic, principal component analysis, univariate and multivariate analysis, and nomogram. The verifications of these two immune-related-gene-signature prognostic models showed consistency in the train and test cohorts of LUAD and LUSC. In addition, patients with LUAD in the low-risk group responded better to immunotherapy than those in the high-risk group. This study revealed two reliable immune-related-gene-signature models that were significantly associated with prognosis and tumor microenvironment cell infiltration in LUAD and LUSC, respectively. Evaluation of the integrated characterization of multiple immune-related genes and pathways could help to predict the response to immunotherapy and monitor immunotherapy strategies.

Project description:BackgroundOral squamous cell carcinoma (OSCC) is the most common type of oral cancer, which remains a major cause of morbidity and mortality in patients with head and neck cancers. However, the critical immune-related signatures and their prognostic values have rarely been investigated.Materials and methodsGene differential analysis was used to measure the differences of gene expression between the groups. Correlation analysis was used to assess the association between the gene expression levels and immune-related risk score/DNA methylation levels. The gene set enrichment analysis (GSEA) was used to identify the pathways or cell types enriched by those identified differentially expressed genes (DEGs).ResultsIn this study, we identified four immune-related gene signatures, including CTSG, TNFRSF4, LCORL, and PLAU, that were significantly associated with the overall survival in OSCC patients from the Cancer Genome Atlas (TCGA) OSCC cohort. Moreover, these four immune-related signatures were differentially expressed between the OSCC and nontumor tissues. The two groups (high and low risk) stratified by the immune-related risk scores had significantly different OS and mortality rates. The gene expression patterns and prognostic values of these immune-related signatures were also verified in two independent validation cohorts. Furthermore, the downregulated genes in the high-risk group (which were also upregulated in the low-risk group) were significantly enriched in the cell type-specific signatures of type 2 T helper cell (Th2), plasmacytoid dendritic cell (pDC), and memory B cell. In contrast, the upregulated genes in the high-score group were enriched in growth factor receptor-related signaling pathways, such as the VEGFA-VEGFR2 signaling pathway, PI3K-Akt signaling pathway, focal adhesion-PI3K-Akt-mTOR signaling pathway, and PDGF pathway, suggesting that those pathways were inversely correlated with immune cell infiltration.ConclusionIn summary, the immune-related signatures had the potential for predicting the risk of OSCC patients. Moreover, the present study also improved our understanding of the association between the growth factor receptor pathways and immune cell infiltration in OSCC.

Project description:BackgroundCuproptosis is a novel form of programmed cell death that differs from other types such as pyroptosis, ferroptosis, and autophagy. It is a promising new target for cancer therapy. Additionally, immune-related genes play a crucial role in cancer progression and patient prognosis. Therefore, our study aimed to create a survival prediction model for lung adenocarcinoma patients based on cuproptosis and immune-related genes. This model can be utilized to enhance personalized treatment for patients.MethodsRNA sequencing (RNA-seq) data of lung adenocarcinoma (LUAD) patients were collected from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. The levels of immune cell infiltration in the GSE68465 cohort were determined using gene set variation analysis (GSVA), and immune-related genes (IRGs) were identified using weighted gene coexpression network analysis (WGCNA). Additionally, cuproptosis-related genes (CRGs) were identified using unsupervised clustering. Univariate COX regression analysis and least absolute shrinkage selection operator (LASSO) regression analysis were performed to develop a risk prognostic model for cuproptosis and immune-related genes (CIRGs), which was subsequently validated. Various algorithms were utilized to explore the relationship between risk scores and immune infiltration levels, and model genes were analyzed based on single-cell sequencing. Finally, the expression of signature genes was confirmed through quantitative real-time PCR (qRT-PCR), immunohistochemistry (IHC), and Western blotting (WB).ResultsWe have identified 5 Oncogenic Driver Genes namely CD79B, PEBP1, PTK2B, STXBP1, and ZNF671, and developed proportional hazards regression models. The results of the study indicate significantly reduced survival rates in both the training and validation sets among the high-risk group. Additionally, the high-risk group displayed lower levels of immune cell infiltration and expression of immune checkpoint compared to the low-risk group.