scholarly journals Identification and validation of a novel zinc finger protein-related gene-based prognostic model for breast cancer

PeerJ ◽  
2021 ◽  
Vol 9 ◽  
pp. e12276
Author(s):  
Min Ye ◽  
Liang Li ◽  
Donghua Liu ◽  
Qiuming Wang ◽  
Yunuo Zhang ◽  
...  

Background Breast invasive carcinoma (BRCA) is a commonly occurring malignant tumor. Zinc finger proteins (ZNFs) constitute the largest transcription factor family in the human genome and play a mechanistic role in many cancers’ development. The prognostic value of ZNFs has yet to be approached systematically for BRCA. Methods We analyzed the data of a training set from The Cancer Genome Atlas (TCGA) database and two validation cohort from GSE20685 and METABRIC datasets, composed of 3,231 BRCA patients. After screening the differentially expressed ZNFs, univariate Cox regression, LASSO, and multiple Cox regression analysis were performed to construct a risk-based predictive model. ESTIMATE algorithm, single-sample gene set enrichment analysis (ssGSEA), and gene set enrichment analyses (GSEA) were utilized to assess the potential relations among the tumor immune microenvironment and ZNFs in BRCA. Results In this study, we profiled ZNF expression in TCGA based BRCA cohort and developed a novel prognostic model based on 14 genes with ZNF relations. This model was composed of high and low-score groups for BRCA classification. Based upon Kaplan-Meier survival curves, risk-status-based prognosis illustrated significant differences. We integrated the 14 ZNF-gene signature with patient clinicopathological data for nomogram construction with accurate 1-, 3-, and 5-overall survival predictive capabilities. We then accessed the Genomics of Drug Sensitivity in Cancer database for therapeutic drug response prediction of signature-defined BRCA patient groupings for our selected TCGA population. The signature also predicts sensitivity to chemotherapeutic and molecular-targeted agents in high- and low-risk patients afflicted with BRCA. Functional analysis suggested JAK STAT, VEGF, MAPK, NOTCH TOLL-like receptor, NOD-like receptor signaling pathways, apoptosis, and cancer-based pathways could be key for ZNF-related BRCA development. Interestingly, based on the results of ESTIMATE, ssGSEA, and GSEA analysis, we elucidated that our ZNF-gene signature had pivotal regulatory effects on the tumor immune microenvironment for BRCA. Conclusion Our findings shed light on the potential contribution of ZNFs to the pathogenesis of BRCA and may inform clinical practice to guide individualized treatment.

2021 ◽  
Author(s):  
Yuqin Qiu ◽  
Xuewei Qi ◽  
Zeyu Liu ◽  
Jingnan Xu ◽  
Xiaogang Wang ◽  
...  

Abstract Background: Hypoxia is widespread in solid tumors and is directly associated with colorectal cancer (CRC) aggressiveness, poor prognosis, and immunotherapy resistance. In this study, we aimed at developing a hypoxia-related marker to improve the prognosis prediction in CRC.Methods: We used gene expression data of CRC samples from the Cancer Genome Atlas Database and the hypoxia gene set to obtain a hypoxia gene expression matrice of 479 CRC patients. The prognostic model was constructed by screening hypoxia risk genes that were significantly associated with prognosis by univariate and multivariate Cox regression analysis. The predictive performance of the prognostic model was evaluated by Kaplan-Meier survival curve analyses and ROC curve analysis and validated in the GSE17536 dataset of Gene Expression Omnibus database. Finally, we analyzed the immune cell infiltration and expression of immunosuppressive genes in CRC patients at high and low risk of hypoxia.Results: We constructed a hypoxia risk prognostic model composed of two hypoxia-related genes (SLC2A3 and ENO3), which was proved to have better sensitivity and specificity after a series of validation. Independent prognostic analysis revealed that the risk score can serve as an independent prognostic factor for CRC. The infiltration of natural killer resting cells, activated master cells and T-cell regulatory cells were significantly increased in the hypoxia high-risk group, and Gene Set Enrichment Analysis showed that gene sets involved in tumor proliferation and differentiation, immune tolerance as well as hypoxia were also significantly enriched in this group. Negatively regulated genes in the Cancer Immunity Cycle, as well as immune checkpoints, were upregulated in the high hypoxia risk group, forming an immunosuppressive microenvironment, and mediating the immune escape. Conclusions: In summary, we constructed and validated a reliable hypoxia risk model that can independently predict the prognosis of CRC patients and reflect the status of the immune microenvironment, which is beneficial for screening CRC prognostic biomarkers and therapeutic targets.


2021 ◽  
Author(s):  
Congli Jia ◽  
Fu Yang ◽  
Ruining Li

Abstract Background: Breast cancer (BC) is the most common cancer among women, with high rates of metastasis and recurrence. Some studies have confirmed that pyroptosis is an immune-related programmed cell death. However, the correlation between the expression of pyroptosis-related genes in BC and its prognosis remains unclear. Methods: In this study, we identified 38 pyroptosis-related genes that were differentially expressed between BC and normal tissues. The prognostic value of each pyroptosis-related gene was evaluated using patient data from The Cancer Genome Atlas (TCGA). The Cox regression method was performed to establish a prognostic model for 16-gene signature, classifying all BC patients in the TCGA database into a low-or high-risk group. Results: The survival rate of BC patients in the high-risk group was significantly lower than that in the low-risk group (P<0.01). Prognostic model is independent prognostic factor for BC patients compared to clinical features. Single sample gene set enrichment analysis (ssGSEA) showed a decrease for immune cells and immune function in the high-risk group. Conclusions: Pyroptosis-related genes influence the tumor immune microenvironment and can predict the prognosis of BC.


2021 ◽  
Vol 27 ◽  
Author(s):  
Aoshuang Qi ◽  
Mingyi Ju ◽  
Yinfeng Liu ◽  
Jia Bi ◽  
Qian Wei ◽  
...  

Background: Complex antigen processing and presentation processes are involved in the development and progression of breast cancer (BC). A single biomarker is unlikely to adequately reflect the complex interplay between immune cells and cancer; however, there have been few attempts to find a robust antigen processing and presentation-related signature to predict the survival outcome of BC patients with respect to tumor immunology. Therefore, we aimed to develop an accurate gene signature based on immune-related genes for prognosis prediction of BC.Methods: Information on BC patients was obtained from The Cancer Genome Atlas. Gene set enrichment analysis was used to confirm the gene set related to antigen processing and presentation that contributed to BC. Cox proportional regression, multivariate Cox regression, and stratified analysis were used to identify the prognostic power of the gene signature. Differentially expressed mRNAs between high- and low-risk groups were determined by KEGG analysis.Results: A three-gene signature comprising HSPA5 (heat shock protein family A member 5), PSME2 (proteasome activator subunit 2), and HLA-F (major histocompatibility complex, class I, F) was significantly associated with OS. HSPA5 and PSME2 were protective (hazard ratio (HR) < 1), and HLA-F was risky (HR > 1). Risk score, estrogen receptor (ER), progesterone receptor (PR) and PD-L1 were independent prognostic indicators. KIT and ACACB may have important roles in the mechanism by which the gene signature regulates prognosis of BC.Conclusion: The proposed three-gene signature is a promising biomarker for estimating survival outcomes in BC patients.


2021 ◽  
Author(s):  
Chuan-Qi Xu ◽  
Kui-Sheng Yang ◽  
Shu-Xian Zhao ◽  
Jian Lv

Abstract Objective: Pancreatic cancer (PC) is one of the most malignant tumors. Cytosolic DNA sensing have been found to play an essential role in tumor. In this study, a cytosolic DNA sensing-related genes (CDSRGs) signature was constructed and the potential mechanisms also been discussed.Methods: The RNA expression and clinical data of PC were obtained from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO). Subsequently, univariate (UCR) and multivariate Cox regression (MCR) analyses were conducted to establish a prognostic model in the TCGA patients, which was verified by GEO patients. Cancer immune infiltrates were investigated via single sample gene set enrichment analysis (ssGSEA) and Tumor Immune Estimation Resource (TIMER). Finally, Gene Set Enrichment Analysis (GSEA) was used to investigate the related signaling pathways.Results: A prognostic model comprising four genes (POLR2E,IL18, MAVS, and FADD) was established. The survival rate of patients in the low-risk group was significantly higher than that of patients in the high-risk group. In addition, CDSRGs-risk score was proved as an independent prognostic factor in PC. Immune infiltrates and drug sensitivity are associated with POLR2E,IL18, MAVS, and FADD expression.Conclusions: In summary, we present and validated a CDSRGs risk model that is an independent prognostic factor and indicates the immune characteristics of PC. This prognostic model may facilitate the personalized treatment and monitoring.


2022 ◽  
Vol 12 ◽  
Author(s):  
Lan-Xin Mu ◽  
You-Cheng Shao ◽  
Lei Wei ◽  
Fang-Fang Chen ◽  
Jing-Wei Zhang

Purpose: This study aims to reveal the relationship between RNA N6-methyladenosine (m6A) regulators and tumor immune microenvironment (TME) in breast cancer, and to establish a risk model for predicting the occurrence and development of tumors.Patients and methods: In the present study, we respectively downloaded the transcriptome dataset of breast cancer from Gene Expression Omnibus (GEO) database and The Cancer Genome Atlas (TCGA) database to analyze the mutation characteristics of m6A regulators and their expression profile in different clinicopathological groups. Then we used the weighted correlation network analysis (WGCNA), the least absolute shrinkage and selection operator (LASSO), and cox regression to construct a risk prediction model based on m6A-associated hub genes. In addition, Immune infiltration analysis and gene set enrichment analysis (GSEA) was used to evaluate the immune cell context and the enriched gene sets among the subgroups.Results: Compared with adjacent normal tissue, differentially expressed 24 m6A regulators were identified in breast cancer. According to the expression features of m6A regulators above, we established two subgroups of breast cancer, which were also surprisingly distinguished by the feature of the immune microenvironment. The Model based on modification patterns of m6A regulators could predict the patient’s T stage and evaluate their prognosis. Besides, the low m6aRiskscore group presents an immune-activated phenotype as well as a lower tumor mutation load, and its 5-years survival rate was 90.5%, while that of the high m6ariskscore group was only 74.1%. Finally, the cohort confirmed that age (p < 0.001) and m6aRiskscore (p < 0.001) are both risk factors for breast cancer in the multivariate regression.Conclusion: The m6A regulators play an important role in the regulation of breast tumor immune microenvironment and is helpful to provide guidance for clinical immunotherapy.


2021 ◽  
Author(s):  
Jing Liu ◽  
Ting Ye ◽  
Xue fang Zhang ◽  
Yong jian Dong ◽  
Wen feng Zhang ◽  
...  

Abstract Most of the malignant melanomas are already in the middle and advanced stages when they are diagnosed, which is often accompanied by the metastasis and spread of other organs.Besides, the prognosis of patients is bleak. The characteristics of the local immune microenvironment in metastatic melanoma have important implications for both tumor progression and tumor treatment. In this study, data on patients with metastatic melanoma from the TCGA and GEO datasets were selected for immune, stromal, and estimate scores, and overlapping differentially expressed genes (DEGs) were screened. A nine-IRGs prognostic model (ALOX5AP, ARHGAP15, CCL8, FCER1G, GBP4, HCK, MMP9, RARRES2 and TRIM22) was established by univariate COX regression, LASSO and multivariate COX regression. Receiver operating characteristic (ROC) curves were used to test the predictive accuracy of the model. Immune infiltration was analyzed by using CIBERSORT, Xcell and ssGSEA in high-risk and low-risk groups. The immune infiltration of the high-risk group was significantly lower than that of the low-risk group. Immune checkpoint analysis revealed that the expression of PDCD1, CTLA4, TIGIT, CD274, HAVR2 and LAG3 were significantly different in groups with different levels of risk scores. WGCNA analysis found that the yellow-green module contained seven genes (ALOX5AP, FCER1G, GBP4, HCK, MMP9, RARRES2 and TRIM22) from the nine-IRG prognostic model, of which the yellow-green module had the highest correlation with risk scores. The results of GO and KEGG suggested that the genes in the yellow-green module were mainly enriched in immune-related biological processes. Finally, we analyzed the prognostic ability and expression characteristics of ALOX5AP, ARHGAP15, CCL8, FCER1G, GBP4, HCK, MMP9, RARRES2 and TRIM22 in metastatic melanoma. Overall, a prognostic model for metastatic melanoma based on the characteristics of the tumor immune microenvironment was established, which was helpful for further studies.It could function well in helping people to understand the characteristics of the immune microenvironment in metastatic melanoma and to find possible therapeutic targets.


2021 ◽  
Author(s):  
Zhiyuan Zheng ◽  
Wei Wu ◽  
Zehang Lin ◽  
Shuhan Liu ◽  
Qiaoqian Chen ◽  
...  

Abstract Background: Ferroptosis is a newly discovered type of programmed cell death that participates in the biological processes of various cancers. However, the mechanism by which ferroptosis modulates acute myeloid leukemia (AML) remains unclear. This study aimed to investigate the role of ferroptosis-related long non-coding RNAs (lncRNAs) in AML and establish a corresponding prognostic model.Methods: RNA-sequencing data and clinicopathological characteristics were obtained from The Cancer Genome Atlas database, and ferroptosis-related genes were obtained from the FerrDb database. The “limma” R package, Cox regression, and the least absolute shrinkage and selection operator were used to determine the ferroptosis-related lncRNA signature with the lowest Akaike information criteria (AIC). The risk score of ferroptosis-related lncRNAs was calculated and patients with AML were divided into high- and low-risk groups based on the median risk score. The Kaplan-Meier curve and Cox regression were used to evaluate the prognostic value of the risk score. Finally, gene set enrichment analysis (GSEA) and single-sample gene set enrichment analysis (ssGSEA) were performed to explore the biological functions of the ferroptosis-related lncRNAs.Results: Seven ferroptosis-related lncRNA signatures were identified in the training group, and Kaplan-Meier and Cox regression analyses confirmed that risk scores were independent prognostic predictors of AML in both the training and validation groups (All P < 0.05). In addition, the area under the curve (AUC) analysis confirmed that the signatures had a good predictive ability for the prognosis of AML. GSEA and ssGSEA showed that the seven ferroptosis-related lncRNAs were related to glutathione metabolism and tumor immunity.Conclusions: In this study, seven novel ferroptosis-related lncRNA signatures (AP001266.2, AC133961.1, AF064858.3, AC007383.2, AC008906.1, AC026771.1, and KIF26B-AS1) were established. These signatures were shown to accurately predict the prognosis of AML, which would provide new insights into strategies for the development of new AML therapies.


2021 ◽  
Vol 19 (1) ◽  
pp. 169-190
Author(s):  
Peiyuan Li ◽  
◽  
Gangjie Qiao ◽  
Jian Lu ◽  
Wenbin Ji ◽  
...  

<abstract> <p>Plasmacytoma variant translocation 1 (PVT1) is involved in multiple signaling pathways and plays an important regulatory role in a variety of malignant tumors. However, its role in the prognosis and immune invasion of bladder urothelial carcinoma (BLCA) remains unclear. This study investigated the expression of PVT1 in tumor tissue and its relationship with immune invasion, and determined its prognostic role in patients with BLCA. Patients were identified from the cancer genome atlas (TCGA). The enrichment pathway and function of PVT1 were explained by gene ontology (GO) term analysis, gene set enrichment analysis (GSEA) and single-sample gene set enrichment analysis (ssGSEA), and the degree of immune cell infiltration was quantified. Kaplan–Meier analysis and Cox regression were used to analyze the correlation between PVT1 and survival rate. PVT1-high BLCA patients had a lower 10-year disease-specific survival (DSS P &lt; 0.05) and overall survival (OS P &lt; 0.05). Multivariate Cox regression analysis showed that PVT1 (high vs. low) (P = 0.004) was an independent prognostic factor. A nomogram was used to predict the effect of PVT1 on the prognosis. PVT1 plays an important role in the progression and prognosis of BLCA and can be used as a medium biomarker to predict survival after cystectomy.</p> </abstract>


2021 ◽  
Author(s):  
HUA HUANG ◽  
Shanshan Xu ◽  
Youran Li ◽  
Yunfei Gu ◽  
Lijiang Ji

Abstract Background: Colorectal cancer (CRC), the commonly seen malignancy, ranks the 3rd place among the causes of cancer-associated mortality. As suggested by more and more studies, long coding RNAs (lncRNAs) have been considered as prognostic biomarkers for CRC. But the significance of hypoxic lncRNAs in predicting CRC prognosis remains unclear.Methods: The gene expressed profiles for CRC cases were obtained based on the Cancer Genome Atlas (TCGA) and applied to estimate the hypoxia score using a single-sample gene set enrichment analysis (ssGSEA) algorithm. Overall survival (OS) of high- and low-hypoxia score group was analyzed by the Kaplan–Meier (KM) plot. To identify differentially expressed lncRNAs (DELs) between two hypoxia score groups, this study carried out differential expression analysis, and then further integrated with the DELs between controls and CRC patients to generate the hypoxia-related lncRNAs for CRC. Besides, prognostic lncRNAs were screened by the univariate Cox regression, which were later utilized for constructing the prognosis nomogram for CRC by adopting the least absolute shrinkage and selection operator (LASSO) algorithm. In addition, both accuracy and specificity of the constructed prognostic signature were detected through the receiver operating characteristic (ROC) analysis. Moreover, our constructed prognosis signature also was validated in the internal testing test. This study operated gene set enrichment analysis (GSEA) for exploring potential biological functions associated with the prognostic signature. Finally, the ceRNA network of the prognostic lncRNAs was constructed.Results: Among 2299 hypoxia-related lncRNAs of CRC in total, LINC00327, LINC00163, LINC00174, SYNPR-AS1, and MIR31HG were identified as prognostic lncRNAs by the univariate Cox regression, and adopted for constructing the prognosis signature for CRC. ROC analysis showed the predictive power and accuracy of the prognostic signature. Additionally, the GSEA revealed that ECM-receptor interaction, PI3K-Akt pathway, phagosome, and Hippo pathway were mostly associated with the high-risk group. 352 miRNAs-mRNAs pairs and 177 lncRNAs-miRNAs were predicted.Conclusion: To conclude , we identified 5 hypoxia-related lncRNAs to establish an accurate prognostic signature for CRC, providing important prognostic markers and therapeutic target.


2021 ◽  
Vol 19 (2) ◽  
pp. 1448-1470
Author(s):  
Jing Liu ◽  
◽  
Xuefang Zhang ◽  
Ting Ye ◽  
Yongjian Dong ◽  
...  

<abstract> <p>Most of the malignant melanomas are already in the middle and advanced stages when they are diagnosed, which is often accompanied by the metastasis and spread of other organs. Besides, the prognosis of patients is bleak. The characteristics of the local immune microenvironment in metastatic melanoma have important implications for both tumor progression and tumor treatment. In this study, data on patients with metastatic melanoma from the TCGA and GEO datasets were selected for immune, stromal, and estimate scores, and overlapping differentially expressed genes were screened. A nine-IRGs prognostic model (ALOX5AP, ARHGAP15, CCL8, FCER1G, GBP4, HCK, MMP9, RARRES2 and TRIM22) was established by univariate COX regression, LASSO and multivariate COX regression. Receiver operating characteristic curves were used to test the predictive accuracy of the model. Immune infiltration was analyzed by using CIBERSORT and Xcell in high-risk and low-risk groups. The immune infiltration of the high-risk group was significantly lower than that of the low-risk group. Immune checkpoint analysis revealed that the expression of PDCD1, CTLA4, TIGIT, CD274, HAVR2 and LAG3 demonstrated the visible difference in groups with different levels of risk scores. WGCNA analysis found that the yellow-green module contained seven genes from the nine-IRG prognostic model, and the yellow-green module had the highest correlation with risk scores. The results of GO and KEGG suggested that the genes in the yellow-green module were mainly enriched in immune-related biological processes. Finally, the expression characteristics of ALOX5AP, ARHGAP15, CCL8, FCER1G, GBP4, HCK, MMP9, RARRES2 and TRIM22 were analyzed between metastatic melanoma and normal samples. Overall, a prognostic model for metastatic melanoma based on the tumor immune microenvironment characteristics was established, which left plenty of space for further studies. It could function well in helping people to understand characteristics of the immune microenvironment in metastatic melanoma.</p> </abstract>


Sign in / Sign up

Export Citation Format

Share Document