Cancer Stemness Associated With Prognosis and the Efficacy of Immunotherapy in Adrenocortical Carcinoma

BackgroundCancer stem cells (CSCs) have been proven to influence drug resistance, recurrence, and metastasis in tumors. Our study aimed to identify stemness-related prognostic biomarkers for new therapeutic strategies in adrenocortical carcinoma.MethodsRNA-seq data and clinical characteristics were downloaded from The Cancer Genome Atlas (TCGA). The stemness indexes, mDNAsi and mRNAsi, were calculated to classify all samples into low-score and high-score groups. Two algorithms, based on the R language, ESTIMATE and single-sample Gene Set Enrichment Analysis (ssGSEA) were used to assess the immune cell infiltration states of adrenocortical carcinoma patients. Weighted Gene Co-expression Network Analysis (WGCNA) was used to find genes that were related to the stemness of cancer. By bioinformatics methods, the correlations between biomarkers capable of predicting immune checkpoint inhibitors (ICIs) responses and stemness of cancer were explored.ResultsHigh-mRNAsi predicted shorter overall survival (OS) and a higher metastatic trend in adrenocortical carcinoma (ACC) patients. Compared with the low-mRNAsi group, the high-mRNAsi group had a lower ImmuneScore and StromalScroe. Twenty-two stemness-related prognostic genes were obtained by WGCNA, which focused on the function of the cell cycle and cell mitosis. Immune cell infiltration, especially CD8+T cell, increased in the low-mRNAsi group compared with the high-mRNAsi group. Lower expression of PD-L1, CTLA-4, and TIGHT was evaluated in the high-mRNAsi group.ConclusionsACC patients with high-mRNAsi have poor prognosis and less immune cell infiltration. Combined with the finding of lower expression of CTLA-4, TIGHT, and PD-L1 in the high-mRNAsi group, we came to the conclusion that stemness index is a potential biomarker to predict the effectiveness of ICIs.

Download Full-text

Gene Alterations of N6-Methyladenosine (m6A) Regulators in Colorectal Cancer: A TCGA Database Study

BioMed Research International ◽

10.1155/2020/8826456 ◽

2020 ◽

Vol 2020 ◽

pp. 1-13

Author(s):

Qian Zhang ◽

Yuping Cai ◽

Vadim Kurbatov ◽

Sajid A. Khan ◽

Lingeng Lu ◽

...

Keyword(s):

Mrna Expression ◽

Messenger Rna ◽

Immune Cell ◽

Disease Free Survival ◽

Enrichment Analysis ◽

Gene Set Enrichment Analysis ◽

The Cancer Genome Atlas ◽

Cell Infiltration ◽

Immune Cell Infiltration ◽

The Relationship

N6-methyladenosine (m6A) plays an important role in many cancers. However, few studies have examined the role of m6A in colorectal CRC. To examine the effect of m6A on CRC, we studied the genome of 591 CRC cases from The Cancer Genome Atlas (TCGA). The relationship between the messenger RNA (mRNA) expression, copy number variation (CNVs), and mutations of m6A “Writers,” “Readers,” and “Erasers,” prognosis, immune cell infiltration, and genetic mutations in CRC cases were analyzed. CNVs and mutations were found in thirteen m6A regulators. As expected, gain and amplification of m6A regulators increased the mRNA expression of these regulators, while deletion led to reduction in the mRNA expression. Moreover, CNVs and mutation of these regulators were significantly associated with APC, TP53, and microsatellite instability (MSI) status ( p < 0.001 , p < 0.001 , and p = 0.029 , respectively). CNVs of m6A regulators also correlated with inferred immune cell infiltration in CRC tissues, especially in colon tissues. Additionally, alterations of RBM15, YTHDF2, YTHDC1, YTHDC2, and METTL14 genes were related to the worse overall survival and disease-free survival (DFS) of CRC patients. Specifically, the deletion status of “Writers” was also correlated to the DFS of CRC patients ( p = 0.02 ). Gene set enrichment analysis found that FTO was involved in mRNA 3 ′ end processing, polyubiquitin binding, and RNA polymerase promoter elongation, while YTHDC1 was related to interferon-alpha and gamma response. In conclusion, a novel relationship was identified between CNVs and mutations of m6A regulators with prognosis and inferred immune function of CRC. These findings will improve the understanding of the relationship of m6A in CRC.

Download Full-text

Hypoxic Characteristic in the Immunosuppressive Microenvironment of Hepatocellular Carcinoma

Frontiers in Immunology ◽

10.3389/fimmu.2021.611058 ◽

2021 ◽

Vol 12 ◽

Author(s):

Zhuomao Mo ◽

Daiyuan Liu ◽

Dade Rong ◽

Shijun Zhang

Keyword(s):

Hepatocellular Carcinoma ◽

Immune Cell ◽

Cox Regression ◽

Cancer Genome ◽

Gene Set Enrichment Analysis ◽

The Cancer Genome Atlas ◽

Cell Infiltration ◽

Immune Checkpoints ◽

Immune Cell Infiltration ◽

Immunosuppressive Microenvironment

Background: Generally, hepatocellular carcinoma (HCC) exists in an immunosuppressive microenvironment that promotes tumor evasion. Hypoxia can impact intercellular crosstalk in the tumor microenvironment. This study aimed to explore and elucidate the underlying relationship between hypoxia and immunotherapy in patients with HCC.Methods: HCC genomic and clinicopathological datasets were obtained from The Cancer Genome Atlas (TCGA-LIHC), Gene Expression Omnibus databases (GSE14520) and International Cancer Genome Consortium (ICGC-LIRI). The TCGA-LIHC cases were divided into clusters based on single sample gene set enrichment analysis and hierarchical clustering. After identifying patients with immunosuppressive microenvironment with different hypoxic conditions, correlations between immunological characteristics and hypoxia clusters were investigated. Subsequently, a hypoxia-associated score was established by differential expression, univariable Cox regression, and lasso regression analyses. The score was verified by survival and receiver operating characteristic curve analyses. The GSE14520 cohort was used to validate the findings of immune cell infiltration and immune checkpoints expression, while the ICGC-LIRI cohort was employed to verify the hypoxia-associated score.Results: We identified hypoxic patients with immunosuppressive HCC. This cluster exhibited higher immune cell infiltration and immune checkpoint expression in the TCGA cohort, while similar significant differences were observed in the GEO cohort. The hypoxia-associated score was composed of five genes (ephrin A3, dihydropyrimidinase like 4, solute carrier family 2 member 5, stanniocalcin 2, and lysyl oxidase). In both two cohorts, survival analysis revealed significant differences between the high-risk and low-risk groups. In addition, compared to other clinical parameters, the established score had the highest predictive performance at both 3 and 5 years in two cohorts.Conclusion: This study provides further evidence of the link between hypoxic signals in patients and immunosuppression in HCC. Defining hypoxia-associated HCC subtypes may help reveal potential regulatory mechanisms between hypoxia and the immunosuppressive microenvironment, and our hypoxia-associated score could exhibit potential implications for future predictive models.

Download Full-text

SERPINE1 associated with remodeling of the tumor microenvironment in colon cancer progression: a novel therapeutic target

BMC Cancer ◽

10.1186/s12885-021-08536-7 ◽

2021 ◽

Vol 21 (1) ◽

Author(s):

Shaokun Wang ◽

Li Pang ◽

Zuolong Liu ◽

Xiangwei Meng

Keyword(s):

Colorectal Cancer ◽

Colon Cancer ◽

Immune Cells ◽

Immune Cell ◽

Cox Regression ◽

Gene Set Enrichment Analysis ◽

The Cancer Genome Atlas ◽

Cell Infiltration ◽

Immune Cell Infiltration ◽

The Relationship

Abstract Background The change of immune cell infiltration essentially influences the process of colorectal cancer development. The infiltration of immune cells can be regulated by a variety of genes. Thus, modeling the immune microenvironment of colorectal cancer by analyzing the genes involved can be more conducive to the in-depth understanding of carcinogenesis and the progression thereof. Methods In this study, the number of stromal and immune cells in malignant tumor tissues were first estimated by using expression data (ESTIMATE) and cell-type identification with relative subsets of known RNA transcripts (CIBERSORT) to calculate the proportion of infiltrating immune cell and stromal components of colon cancer samples from the Cancer Genome Atlas database. Then the relationship between the TMN Classification and prognosis of malignant tumors was evaluated. Results By investigating differentially expressed genes using COX regression and protein-protein interaction network (PPI), the candidate hub gene serine protease inhibitor family E member 1 (SERPINE1) was found to be associated with immune cell infiltration. Gene Set Enrichment Analysis (GSEA) further projected the potential pathways with elevated SERPINE1 expression to carcinogenesis and immunity. CIBERSORT was subsequently utilized to investigate the relationship between the expression differences of SERPINE1 and immune cell infiltration and to identify eight immune cells associated with SERPINE1 expression. Conclusion We found that SERPINE1 plays a role in the remodeling of the colon cancer microenvironment and the infiltration of immune cells.

Download Full-text

ARID1A Is a Prognostic Biomarker and Associated with Immune Infiltrates in Hepatocellular Carcinoma

Canadian Journal of Gastroenterology and Hepatology ◽

10.1155/2022/3163955 ◽

2022 ◽

Vol 2022 ◽

pp. 1-14

Author(s):

Yuanyuan Feng ◽

Xinfang Tang ◽

Changcheng Li ◽

Ying Su ◽

Xiaoyu Wang ◽

...

Keyword(s):

Hepatocellular Carcinoma ◽

Overall Survival ◽

Immune Cell ◽

Enrichment Analysis ◽

Gene Set Enrichment Analysis ◽

Cell Infiltration ◽

Migration And Invasion ◽

Risk Indicator ◽

Immune Cell Infiltration

Objective. ARID1A has been discovered as a potential cancer biomarker. But its role in hepatocellular carcinoma (HCC) is subject to considerable dispute. Methods. The relationship between ARID1A and clinical factors was investigated. Clinicopathological variables related to overall survival in HCC subjects were identified using Cox and Kaplan–Meier studies. The connection between immune infiltrating cells and ARID1A expression was investigated using the tumor Genome Atlas (TCGA) dataset for gene set enrichment analysis (GSEA). Finally, a cell experiment was used to confirm it. Results. The gender and cancer topography (T) categorization of HCC were linked to increased ARID1A expression. Participants with advanced levels of ARID1A expression had a worse prognosis than someone with lower levels. ARID1A was shown to be a risk indicator of overall survival on its own. ARID1A expression is inversely proportional to immune cell infiltration. In vitro, decreasing ARID1A expression substantially slowed the cell cycle and decreased HCC cell proliferation, migration, and invasion. Conclusion. The expression of ARID1A could be used to predict the outcome of HCC. It is closely related to tumor immune cell infiltration.

Download Full-text

RNA N6-Methyladenosine Regulator-Mediated Methylation Modifications Pattern and Immune Infiltration Features in Glioblastoma

Frontiers in Oncology ◽

10.3389/fonc.2021.632934 ◽

2021 ◽

Vol 11 ◽

Author(s):

Yimin Pan ◽

Kai Xiao ◽

Yue Li ◽

Yuzhe Li ◽

Qing Liu

Keyword(s):

Immune Cell ◽

Enrichment Analysis ◽

Gene Set Enrichment Analysis ◽

Cell Infiltration ◽

Immune Cell Infiltration ◽

Lasso Regression ◽

Gene Set Enrichment ◽

Clinical Prognosis ◽

Prognostic Assessment ◽

The Relationship

Glioblastoma (GBM) is a group of intracranial neoplasms with intra-tumoral heterogeneity. RNA N6-methyladenosine (m6A) methylation modification reportedly plays roles in immune response. The relationship between the m6A modification pattern and immune cell infiltration in GBM remains unknown. Utilizing expression data of GBM patients, we thoroughly explored the potential m6A modification pattern and m6A-related signatures based on 21 regulators. Thereafter, the m6A methylation modification-based prognostic assessment pipeline (MPAP) was constructed to quantitatively assess GBM patients’ clinical prognosis combining the Robustness and LASSO regression. Single-sample gene-set enrichment analysis (ssGSEA) was used to estimate the specific immune cell infiltration level. We identified two diverse clusters with diverse m6A modification characteristics. Based on differentially expressed genes (DEGs) within two clusters, m6A-related signatures were identified to establish the MPAP, which can be used to quantitatively forecast the prognosis of GBM patients. In addition, the relationship between 21 m6A regulators and specific immune cell infiltration was demonstrated in our study and the m6A regulator ELAVL1 was determined to play an important role in the anticancer response to PD-L1 therapy. Our findings indicated the relationship between m6A methylation modification patterns and tumor microenvironment immune cell infiltration, through which we could comprehensively understand resistance to multiple therapies in GBM, as well as accomplish precise risk stratification according to m6A-related signatures.

Download Full-text

Prediction of Hepatocellular Carcinoma Prognosis and Immune Cell Infiltration Using Gene Signature Associated with Inflammatory Response

Computational and Mathematical Methods in Medicine ◽

10.1155/2022/2415129 ◽

2022 ◽

Vol 2022 ◽

pp. 1-24

Author(s):

Bin-Bin Da ◽

Shuai Luo ◽

Ming Huang ◽

Fei Song ◽

Rong Ding ◽

...

Keyword(s):

Hepatocellular Carcinoma ◽

Inflammatory Response ◽

Immune Cell ◽

Cox Regression ◽

Inflammatory Responses ◽

Gene Set Enrichment Analysis ◽

The Cancer Genome Atlas ◽

Cell Infiltration ◽

Immune Cell Infiltration ◽

Differential Analysis

It has been demonstrated that the inflammatory response influences cancer development and can be used as a prognostic biomarker in various tumors. However, the relevance of genes associated with inflammatory responses in hepatocellular carcinoma (HCC) remains unknown. The Cancer Genome Atlas (TCGA) database was analyzed using weighted gene coexpression network analysis (WGCNA) and differential analysis to discover essential inflammatory response-related genes (IFRGs). Cox regression studies, both univariate and multivariate, were employed to develop a prognostic IFRGs signature. Additionally, Gene Set Enrichment Analysis (GSEA) was used to deduce the biological function of the IFRGs signature. Finally, we estimated immune cell infiltration using a single sample GSEA (ssGSEA) and x-cell. Our results revealed that, among the major HCC IFRGs, two (DNASE1L3 and KLKB1) were employed to create a predictive IFRG signature. The IFRG signature could correctly predict overall survival (O.S) as per Kaplan-Meier time-dependent roc curves analysis. It was also linked to pathological tumor stage and T stage and might be used as a prognostic predictor in HCC. GSEA analysis concluded that the IFRG signature might influence the immune response in HCC. Immunological cell infiltration and immune checkpoint molecule expression differed in the high-risk and low-risk groups. As a result of our findings, DNASILE may play a role in the tumor microenvironment. However, more research is necessary to confirm the role of DNASE1L3 and KLKB1.

Download Full-text

Landscape of Immune Cell Infiltration in Cervical Cancer

10.21203/rs.3.rs-584271/v1 ◽

2021 ◽

Author(s):

Shasha Shi ◽

Fu Peng ◽

Chenghao Yu

Keyword(s):

Cervical Cancer ◽

Immune Cell ◽

Checkpoint Inhibitors ◽

Expression Profiles ◽

The United States ◽

The Cancer Genome Atlas ◽

Cell Infiltration ◽

Immune Checkpoints ◽

Immune Cell Infiltration ◽

Cancer Tissues

Abstract BackgroundCervical cancer is a life-threatening cancer among women. It is the second most prevalent malignant tumor in women. It ranks high in cancer deaths among women worldwide, including in the United States. Immune checkpoint inhibitors have emerged as an important therapeutic approach to treat several cancers, including cervical cancer. Notably, the development and progress of cervical cancer may be related to sustained immune response. This underlines the need to clarify immune cell infiltration (ICI) in cervical cancer tissues. MethodsIn this study, disease-related information of 964 cervical cancer patients was first retrieved from Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA) database. We utilized bioinformatics data to analyze the expression profiles of immune genes in cervical cancer tissues. ResultsPatients were divided into high and low groups according to ICI score. High ICI scores corresponded with activation of immune signaling pathways and high tumor mutation burden (TMB), which was related to better prognosis of G1-2 cervical cancer. In addition, most immune checkpoints and immuno-related genes such as CD274, CD8A, CXCL10, etc. were over-expressed in the high ICI group. ConclusionsThis study demonstrated that ICI score can accurately predict the prognosis of cervical cancer. Understanding ICI patterns will deepen our understanding of tumor microenvironment (TME) of cervical cancer, which may create the foundation for the development of efficient immunotherapeutic strategies against the cancer.

Download Full-text

A hypoxia-related signature for clinically predicting diagnosis, prognosis and immune microenvironment of hepatocellular carcinoma patients.

10.21203/rs.3.rs-17783/v1 ◽

2020 ◽

Author(s):

Baohui Zhang ◽

Bufu Tang ◽

Jianyao Gao ◽

Jiatong Li ◽

Lingming Kong ◽

...

Keyword(s):

Hepatocellular Carcinoma ◽

Immune Cell ◽

Cox Regression ◽

Enrichment Analysis ◽

High Risk Group ◽

Cancer Genome ◽

Gene Set Enrichment Analysis ◽

The Cancer Genome Atlas ◽

Diagnostic Model ◽

Potential Biomarker

Abstract Background Hypoxia plays an indispensable role in the development of hepatocellular carcinoma (HCC). However, there are few studies on the application of hypoxia molecules in the prognosis predicting of HCC. We aimed to identify the hypoxia-related genes in HCC and construct reliable models for diagnosis, prognosis and recurrence of HCC patients as well as exploring the potential mechanism.Methods Differentially expressed genes (DEGs) analysis was performed using The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) database and four clusters were determined by a consistent clustering analysis. Three DEGs closely related to overall survival（OS）were identified using Cox regression and LASSO analysis and the hypoxia-related signature was developed and validated in TCGA and International Cancer Genome Consortium (ICGC) database. Then the Gene Set Enrichment Analysis (GSEA) was performed to explore signaling pathways regulated by the signature and the CIBERSORT was used for estimating the fractions of immune cell types.Results A total of 397 hypoxia-related DEGs were detected and three genes (PDSS1, CDCA8 and SLC7A11) were selected to construct a prognosis, recurrence and diagnosis model. Then patients were divided into high- and low-risk groups. Our hypoxia-related signature was significantly associated with worse prognosis and higher recurrence rate. The diagnostic model also accurately distinguished HCC from normal samples and nodules. Furthermore, the hypoxia-related signature could positively regulate immune response and the high-risk group had higher fractions of macrophages, B memory cells and follicle-helper T cells, and exhibited higher expression of immunocheckpoints such as PD1and PDL1.Conclusions Altogether, our study showed that hypoxia-related signature is a potential biomarker for diagnosis, prognosis and recurrence of HCC, and it provided an immunological perspective for developing personalized therapies.

Download Full-text

P4HA1, a Prognostic Biomarker that Correlates With Immune Infiltrates in Lung Adenocarcinoma and Pan-Cancer

Frontiers in Cell and Developmental Biology ◽

10.3389/fcell.2021.754580 ◽

2021 ◽

Vol 9 ◽

Author(s):

Qi Zhao ◽

Junfeng Liu

Keyword(s):

Prognostic Value ◽

Immune Cell ◽

Prognostic Biomarker ◽

Enrichment Analysis ◽

Tissue Expression ◽

The Cancer Genome Atlas ◽

Cell Infiltration ◽

Pathway Enrichment Analysis ◽

Immune Cell Infiltration ◽

Pan Cancer

Objective: Prolyl 4-hydroxylase, alpha polypeptide I (P4HA1), a key enzyme in collagen synthesis, comprises two identical alpha subunits and two beta subunits. However, the immunomodulatory role of P4HA1 in tumor immune microenvironment (TIME) remains unclear. This study aimed to evaluate the prognostic value of P4HA1 in pan-cancer and explore the relationship between P4HA1 expression and TIME.Methods: P4HA1 expression, clinical features, mutations, DNA methylation, copy number alteration, and prognostic value in pan-cancer were investigated using the Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression data. Pathway enrichment analysis of P4HA1 was performed using R package “clusterProfiler.” The correlation between immune cell infiltration level and P4HA1 expression was analyzed using three sources of immune cell infiltration data, including ImmuCellAI database, TIMER2 database, and a published work.Results: P4HA1 was substantially overexpressed in most cancer types. P4HA1 overexpression was associated with poor survival in patients. Additionally, we discovered that P4HA1 expression was positively associated with infiltration levels of immunosuppressive cells, such as tumor-associated macrophages, cancer-associated fibroblasts, nTregs, and iTregs, and negatively correlated with CD8+ T and NK cells in pan-cancer.Conclusions: Our results highlighted that P4HA1 might serve as a potential prognostic biomarker in pan-cancer. P4HA1 overexpression is indicative of an immunosuppressive microenvironment. P4HA1 may be a potential target of immunotherapy.

Download Full-text

Comprehensive Analysis of the Immune and Prognostic Implication of COL6A6 in Lung Adenocarcinoma

Frontiers in Oncology ◽

10.3389/fonc.2021.633420 ◽

2021 ◽

Vol 11 ◽

Author(s):

Yi Ma ◽

Mantang Qiu ◽

Haifa Guo ◽

Haiming Chen ◽

Jiawei Li ◽

...

Keyword(s):

T Cell ◽

Lung Adenocarcinoma ◽

Immune Cell ◽

Enrichment Analysis ◽

Gene Signature ◽

Gene Set Enrichment Analysis ◽

Cell Infiltration ◽

Gene Copy ◽

Immune Cell Infiltration ◽

Term Survival

Collagen type VI alpha 6 chain (COL6A6), a novel collagen, has been considered as a tumor suppressor and therapeutic target in several tumors. However, the functional role of COL6A6 in immune cell infiltration and prognostic value in lung adenocarcinoma (LUAD) remains unknown. Here, we evaluated COL6A6 expression and its impact on survival among LUAD patients from The Cancer Genome Atlas (TCGA) and several other databases. COL6A6 was downregulated in LUAD tissues compared to normal tissues at both mRNA and protein levels. COL6A6 expression was negatively associated with pathological stage, tumor stage, and lymph node metastasis. High COL6A6 expression was a favorable prognostic factor in LUAD. Next, we explored the associations between COL6A6 expression and immune cell infiltration. COL6A6 expression was positively associated with the infiltration of B cells, T cells, neutrophils and dendritic cells. Additionally, the immune cell infiltration levels were associated with COL6A6 gene copy number in LUAD. Consistently, gene set enrichment analysis showed that various immune pathways were enriched in the LUAD samples with high COL6A6 expression, including pathways related to T cell activation and T cell receptor signaling. The impacts of COL6A6 on immune activity were further assessed by enrichment analysis of 50 COL6A6-associated immunomodulators. Thereafter, using Cox regression, we identified a seven-gene risk prediction signature based on the COL6A6-associated immunomodulators. The resulting risk score was an independent prognostic predictor in LUAD. Receiver operating characteristic curve analysis confirmed that the seven-gene signature had good prognostic accuracy in the TCGA-LUAD cohort and a Gene Expression Omnibus dataset. Finally, we constructed a clinical nomogram to predict long-term survival probabilities, and calibration curves verified its accuracy. Our findings highlight that COL6A6 is involved in tumor immunity, suggesting COL6A6 may be a potential immunotherapeutic target in LUAD. The proposed seven-gene signature is a promising prognostic biomarker in LUAD.

Download Full-text