Dataset Information

Stemness Related Genes Revealed by Network Analysis Associated With Tumor Immune Microenvironment and the Clinical Outcome in Lung Adenocarcinoma.

ABSTRACT: Lung adenocarcinoma (LUAD) is one of the leading fatal malignancy with high morbidity and mortality worldwide. However, due to its complicated mechanism and lack of effective clinical therapeutics, early diagnosis and prognosis are still unsatisfactory. Most of the previous studies focused on cancer stem cells (CSCs), the relationship between cancer stemness (stem-like characteristics) and anti-tumor immunity has not been clearly revealed. Therefore, this study aimed to comprehensively analyze the role of cancer stemness and tumor microenvironment (TME) in LUAD using weighted gene co-expression network analysis (WGCNA). We constructed a gene co-expression network, identified key modules, and hub genes, and further explored the relationship between hub gene expression and cancer immunological characteristics through a variety of algorithms, including Estimation of STromal and Immune cells in MAlignant Tumor tissues using Expression data (ESTIMATE) and Gene Set Enrichment Analysis (GSEA). The hub genes were renamed stemness related genes (SRGs), whose functions were examined at the transcription and protein levels through survival analysis with additional samples, Oncomine database, immunohistochemistry, single cell RNA sequencing (scRNA-seq) and single-sample Gene Set Enrichment Analysis (ssGSEA). Subsequently, Tumor Immune Dysfunction and Exclusion (TIDE) and Connectivity Map (CMap) were implemented for treatment and prognosis analyses. As a result, 15 co-expressed SRGs (CCNA2, CCNB1, CDC20, CDCA5, CDCA8, FEN1, KIF2C, KPNA2, MCM6, NUSAP1, RACGAP1, RRM2, SPAG5, TOP2A, and TPX2) were identified. The overexpression of which was discovered to be associated with reduced immune infiltration in LUAD. It was discovered that there was a general negative correlation between cancer stemness and immunity. The expression of SRGs could probably affect our tumor occurrence, progression, the efficacy of chemotherapy and immunotherapy, and clinical outcomes. In conclusion, the 15 SRGs reported in our study may be used as potential candidate biomarkers for prognostic indicators and therapeutic targets after further validation.

SUBMITTER: Zeng H

PROVIDER: S-EPMC7525184 | biostudies-literature | 2020

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Stemness Related Genes Revealed by Network Analysis Associated With Tumor Immune Microenvironment and the Clinical Outcome in Lung Adenocarcinoma.

Zeng Hao H Ji Jianrui J Song Xindi X Huang Yeqian Y Li Hui H Huang Juan J Ma Xuelei X

Frontiers in genetics 20200916

Lung adenocarcinoma (LUAD) is one of the leading fatal malignancy with high morbidity and mortality worldwide. However, due to its complicated mechanism and lack of effective clinical therapeutics, early diagnosis and prognosis are still unsatisfactory. Most of the previous studies focused on cancer stem cells (CSCs), the relationship between cancer stemness (stem-like characteristics) and anti-tumor immunity has not been clearly revealed. Therefore, this study aimed to comprehensively analyze t ...[more]

PMID: 33193623

Similar Datasets

Project description:PurposeLung cancer tissue includes tumor tissue, stromal cells, immune cells, and epithelial cells. These nontumor cells dilute the tumor purity in lung cancer tissues. Tumor purity plays an essential role in the immune response to lung cancer. At present, the biological processes related to the purity of lung cancer tumors remains unclear.MethodsWe measured tumor purity in 486 lung carcinoma tissues from TCGA-LUAD FPKM by using the "estimate" R package. Lung carcinoma tumor mutation burden was calculated by analyzing TCGA single nucleotide polymorphism data. The immune cell proportion was also evacuated via the CIBERSORT method. Lung carcinoma samples with P < 0.05 were considered significant. Based on the tumor purity and lung carcinoma gene matrix, we performed weighted gene coexpression network analysis (WGCNA), and the tumor purity-related module was identified. Then, we analyzed the functions of the factors involved in the module. We screened the coexpressed factors related to clinical outcome and immunophenotype. Finally, expression levels of these factors were measured at tissue and single-cell levels.ResultsA lung cancer tumor purity correlated coexpression network was determined. Five coexpressed genes (CD4, CD53, EVI2B, PLEK, and SASH3) were identified as tumor purity coexpressed genes that negatively correlated with tumor purity. Because the factors in the coexpression network often participate in similar biological processes, we found that CD4, CD53, EVI2B, PLEK, and SASH3 were most related to positive regulation of cytokine production and interleukin-2 production through functional enrichment. In a clinical phenotype analysis, we found that these five factors can be used as independent prognostic risk factors. We found that these factors were significantly negatively correlated with tumor purity and positively correlated with the immune score in the immunophenotyping analysis. Using GSEA analysis, we found that the antigen processing and presentation pathway were related to the five tumor coexpressed genes mentioned above. SASH3 and CD53 were used to conduct a prognostic model based on the interaction analysis of the Support Vector Machine and the Least Absolute Shrinkage and Selection Operator. SASH3 was verified to be related to CD8A using a single-cell analysis.ConclusionTumor purity-related coexpression factors in the tumor microenvironment have essential clinical, genomic, and biological significance in lung cancer. These coexpression factors (SASH3 and CD53) can be used to classify tumor purity phenotypes and to predict clinical outcomes.

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:BackgroundLung adenocarcinoma (LUAD) is the most common subtype of lung cancer which typically exhibits a diverse progression trajectory. Our study sought to explore the cell differentiation trajectory of LUAD and its clinical relevance.MethodsUtilizing a single-cell RNA-sequencing dataset (GSE117570), we identified LUAD cells of distinct differential status along with differentiation-related genes (DRGs). DRGs were applied to the analysis of bulk-tissue RNA-sequencing dataset (GSE72094) to classify tumors into different subtypes, whose clinical relevance was further analyzed. DRGs were also applied to gene co-expression network analysis (WGCNA) using another bulk-tissue RNA-sequencing dataset (TCGA-LUAD). Genes from modules that demonstrated a significant correlation with clinical traits and were differentially expressed between normal tissue and tumors were identified. Among these, genes with significant prognostic relevance were used for the development of a prognostic nomogram, which was tested on TCGA-LUAD dataset and validated in GSE72094. Finally, CCK-8, EdU, cell apoptosis, cell colony formation, and Transwell assays were used to verify the functions of the identified genes.ResultsFour clusters of cells with distinct differentiation status were characterized, whose DRGs were predominantly correlated with pathways of immune regulation. Based on DRGs, tumors could be clustered into four subtypes associated with distinct immune microenvironment and clinical outcomes. DRGs were categorized into four modules. A total of nine DRGs (SFTPB, WFDC2, HLA-DPA1, TIMP1, MS4A7, HLA-DQA1, VCAN, KRT8, and FABP5) with most significant survival-predicting power were integrated to develop a prognostic model, which outperformed the traditional parameters in predicting clinical outcomes. Finally, we verified that knockdown of WFDC2 inhibited proliferation, migration, and invasion but promoted the apoptosis of A549 cells in vitro.ConclusionThe cellular composition and cellular differentiation status of tumor mass can predict the clinical outcomes of LUAD patients. It also plays an important role in shaping the tumor immune microenvironment.

Project description:Background:Prognostic genes in the tumor microenvironment play an important role in immune biological processes and the response of cancer to immunotherapy. Thus, we aimed to assess new biomarkers that are associated with immune/stromal cells in lung adenocarcinomas (LUAD) using the ESTIMATE algorithm, which also significantly affects the prognosis of cancer. Methods:The RNA sequencing (RNA-Seq) and clinical data of LUAD were downloaded from the the Cancer Genome Atlas (TCGA ). The immune and stromal scores were calculated for each sample using the ESTIMATE algorithm. The LUAD gene chip expression profile data and the clinical data (GSE37745, GSE11969, and GSE50081) were downloaded from the Gene Expression Omnibus (GEO) for subsequent validation analysis. Differentially expressed genes were calculated between high and low score groups. Univariate Cox regression analysis was performed on differentially expressed genes (DEGs) between the two groups to obtain initial prognosis genes. These were verified by three independent LUAD cohorts from the GEO database. Multivariate Cox regression was used to identify overall survival-related DEGs. UALCAN and the Human Protein Atlas were used to analyze the mRNA /protein expression levels of the target genes. Immune cell infiltration was evaluated using the Tumor Immune Estimation Resource (TIMER) and CIBERSORT methods, and stromal cell infiltration was assessed using xCell. Results:In this study, immune scores and stromal scores are significantly associated with the clinical characteristics of LUAD, including T stage, M stage, pathological stage, and overall survival time. 530 DEGs (18 upregulated and 512 downregulated) were found to coexist in the difference analysis with the immune scores and stromal scores subgroup. Univariate Cox regression analysis showed that 286 of the 530 DEGs were survival-related genes (p < 0.05). Of the 286 genes initially identified, nine prognosis-related genes (CSF2RB, ITK, FLT3, CD79A, CCR4, CCR6, DOK2, AMPD1, and IGJ) were validated from three separate LUAD cohorts. In addition, functional analysis of DEGs also showed that various immunoregulatory molecular pathways, including regulation of immune response and the chemokine signaling pathways, were involved. Five genes (CCR6, ITK, CCR4, DOK2, and AMPD1) were identified as independent prognostic indicators of LUAD in specific data sets. The relationship between the expression levels of these genes and immune genes was assessed. We found that CCR6 mRNA and protein expression levels of LUAD were greater than in normal tissues. We evaluated the infiltration of immune cells and stromal cells in groups with high and low levels of expression of CCR6 in the TCGA LUAD cohort. In summary, we found a series of prognosis-related genes that were associated with the LUAD tumor microenvironment.

Dataset Information

Stemness Related Genes Revealed by Network Analysis Associated With Tumor Immune Microenvironment and the Clinical Outcome in Lung Adenocarcinoma.

Publications

Stemness Related Genes Revealed by Network Analysis Associated With Tumor Immune Microenvironment and the Clinical Outcome in Lung Adenocarcinoma.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets