Dataset Information

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

ABSTRACT:

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.

SUBMITTER: Ye M

PROVIDER: S-EPMC8530103 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:Previous literatures have suggested the importance of inflammatory response during lung adenocarcinoma (LUAD) development. This study aimed at exploring the inflammation-related genes and developing a prognostic signature for predicting the prognosis of LUAD. Survival‑associated inflammation-related genes were identified by univariate Cox regression analysis in the dataset of The Cancer Genome Atlas (TCGA). The least absolute shrinkage and selection operator (LASSO) penalized Cox regression model was used to derive a risk signature which is significantly negatively correlated with OS and divide samples into high-, medium- and low-risk group. Univariate and multivariate Cox analyses suggested that the level of risk group was an independent prognostic factor of the overall survival (OS). Time-dependent receiver operating characteristic (ROC) curve indicated the AUC of 1-, 3- and 5-years of the risk signature was 0.715, 0.719, 0.699 respectively. A prognostic nomogram was constructed by integrating risk group and clinical features. The independent dataset GSE30219 of Gene Expression Omnibus (GEO) was used for verification. We further explored the differences among risk groups in Gene set enrichment analysis (GSEA), tumor mutation and tumor microenvironment. Furthermore, Single Sample Gene Set Enrichment Analysis (ssGSEA) and the results of Cell-type Identification By Estimating Relative Subsets Of RNA Transcripts (CIBERSORT) suggested the status of immune cell infiltration was highly associated with risk groups. We demonstrated the prediction effect of CTLA-4 and PD-1/PD-L1 inhibitors in the low-risk group was better than that in the high-risk group using two methods of immune score include immunophenoscore from The Cancer Immunome Atlas (TCIA) and TIDE score from Tumor Immune Dysfunction and Exclusion (TIDE). In addition, partial targeted drugs and chemotherapy drugs for lung cancer had higher drug sensitivity in the high-risk group. Our findings provide a foundation for future research targeting inflammation-related genes to predictive prognosis and some reference significance for the selection of immunotherapy and drug regimen for lung adenocarcinoma.

Project description:BackgroundMultiple myeloma (MM) is a malignant hematopoietic disease that is usually incurable. RNA-binding proteins (RBPs) are involved in the development of many tumors, but their prognostic significance has not been systematically described in MM. Here, we developed a prognostic signature based on eight RBP-related genes to distinguish MM cohorts with different prognoses.MethodAfter screening the differentially expressed RBPs, univariate Cox regression was performed to evaluate the prognostic relevance of each gene using The Cancer Genome Atlas (TCGA)-Multiple Myeloma Research Foundation (MMRF) dataset. Lasso and stepwise Cox regressions were used to establish a risk prediction model through the training set, and they were validated in three Gene Expression Omnibus (GEO) datasets. We developed a signature based on eight RBP-related genes, which could classify MM patients into high- and low-score groups. The predictive ability was evaluated using bioinformatics methods. Gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment, and gene set enrichment analyses were performed to identify potentially significant biological processes (BPs) in MM.ResultThe prognostic signature performed well in the TCGA-MMRF dataset. The signature includes eight hub genes: HNRNPC, RPLP2, SNRPB, EXOSC8, RARS2, MRPS31, ZC3H6, and DROSHA. Kaplan-Meier survival curves showed that the prognosis of the risk status showed significant differences. A nomogram was constructed with age; B2M, LDH, and ALB levels; and risk status as prognostic parameters. Receiver operating characteristic (ROC) curve, C-index, calibration analysis, and decision curve analysis (DCA) showed that the risk module and nomogram performed well in 1, 3, 5, and 7-year overall survival (OS). Functional analysis suggested that the spliceosome pathway may be a major pathway by which RBPs are involved in myeloma development. Moreover, our signature can improve on the R-International Staging System (ISS)/ISS scoring system (especially for stage II), which may have guiding significance for the future.ConclusionWe constructed and verified the 8-RBP signature, which can effectively predict the prognosis of myeloma patients, and suggested that RBPs are promising biomarkers for MM.

Project description:AimsColon cancer (CRC), with high morbidity and mortality, is a common and highly malignant cancer, which always has a bad prognosis. So it is urgent to employ a reasonable manner to assess the prognosis of patients. We developed and validated a gene model for predicting CRC risk.MethodsThe Gene Expression Omnibus (GEO) database was used to extract the gene expression profiles of CRC patients (N = 181) from GEO to identify genes that were differentially expressed between CRC patients and controls and then stable signature genes by firstly using both robust likelihood-based modeling with 1000 iterations and random survival forest variable hunting algorithms. Cluster analysis using the longest distance method was drawn out, and Kaplan-Meier (KM) survival analysis was used to compare the clusters. Meanwhile, the risk score was evaluated in three independent datasets including the GEO and Illumina HiSeq sequencing platforms. The corresponding risk index was calculated, and samples were clustered into high- and low-risk groups according to the median. And survival ROC analysis was used to evaluate the prognostic model. Finally, the Gene Set Enrichment Analysis (GSEA) was performed for further functional enrichment analyses.ResultsA 10-gene model was obtained, including 7 negative impact factors (SLC39A14, AACS, ERP29, LAMP3, TMEM106C, TMED2, and SLC25A3) and 3 positive ones (CNPY2, GRB10, and PBK), which related with several important oncogenic pathways (KRAS signaling, TNF-α signaling pathway, and WNT signaling pathway) and several cancer-related cellular processes (epithelial mesenchymal transition and cellular apoptosis). By using colon cancer datasets from The Cancer Genome Atlas (TCGA), the model was validated in KM survival analysis (P ≤ 0.001) and significant analysis with recurrence time (P = 0.0018).ConclusionsThis study firstly developed a stable and effective 10-gene model by using novel combined methods, and CRC patients might be able to use it as a prognostic marker for predicting their survival and monitoring their long-term treatment.

Project description:Background: Ferroptosis is a newly discovered form of regulated cell death with distinct properties and recognizing functions involved in physical conditions or various diseases, including cancers. However, the relationship between gliomas and ferroptosis-related lncRNAs (FRLs) remains unclear. Methods: We collected a total of 1850 samples from The Cancer Genome Atlas (TCGA) and Genotype Tissue Expression (GTEX) databases, including 698 tumor and 1,152 normal samples. A list of ferroptosis-related genes was downloaded from the Ferrdb website. Differentially expressed FRLs (DEFRLS) were analyzed using the “limma” package in R software. Subsequently, prognosis-related FRLs were obtained by univariate Cox analysis. Finally, a prognostic model based on the 3 FRLs was constructed using Cox regression analysis with the least absolute shrinkage and selection operator (LASSO) algorithm. The prognostic power of the model was assessed using receiver operating characteristic (ROC) curve analysis and Kaplan-Meier (K-M) survival curve analysis. In addition, we further explored the relationship of the immune landscape and somatic mutations to prognostic model characteristics. Finally, we validated the function of LINC01426 in vitro. Results: We successfully constructed a 3-FRLs signature and classified glioma patients into high-risk and low-risk groups based on the risk score calculated from this signature. Compared with traditional clinicopathological features [age, sex, grade, isocitrate dehydrogenase (IDH) status], the prognostic accuracy of this model is more stable and stronger. Additionally, the model had stable predictive power for overall survival over a 5-year period. In addition, we found significant differences between the two groups in cellular immunity, the numbers of many immune cells, including NK cells, CD4+, CD8+ T-cells, and macrophages, and the expression of many immune-related genes. Finally, the two groups were also significantly different at the level of somatic mutations, especially in glioma prognosis-related genes such as IDH1 and ATRX, with lower mutation rates in the high-risk group leading to poorer prognosis. Finally, we found that the ferroptosis process of glioma cells was inhibited after knocking down the expression of LINC01426. Conclusion: The proposed 3-FRL signature is a promising biomarker for predicting prognostic features in glioma patients.

Project description:BackgroundColorectal cancer (CRC) is one of the most commonly diagnosed cancers in the world. Mitophagy is associated with tumorigenesis and development of malignancy. However, the specific role of mitophagy has yet not been systematically explored in CRC.MethodsThe RNA-sequencing dataset of CRC from The Cancer Genome Atlas (TCGA) and microarray data of gene expression profiles of CRC from Gene Expression Omnibus (GEO) were downloaded. Mitophagy-related gene sets were obtained from the Pathway Unification database. The package "limma" was used for differential gene expression analysis. Kaplan-Meier (KM) survival analyses were utilized to evaluate the prognostic value of the mitophagy regulators. Single-sample gene set enrichment analysis (ssGSEA) was used to estimate the infiltrating immune cells and the activity of immune response. The ConsensusClusterPlus algorithm was used to determine mitophagy-related subtypes. Principal component analysis (PCA) was used to create composite measurement of mitophagy scores. The R packages "survminer" and "ReGlot" were used to plot the nomogram and calibration curves.ResultsIntegrated analysis of the GEO and TCGA databases revealed some common differentially expressed genes (DEGs) in CRC. MFN2, UBB, PINK1, and PRKN were significantly downregulated in CRC samples as compared to normal samples, and other genes were significantly upregulated in CRC samples. KM survival analyses showed that high expression of ATG12 and MAP1LC3B predicted a poor prognosis, whereas high expression of TOMM22 and TOMM40 predicted a better prognosis. Mitophagy showed significant correlation with immune-related pathways in CRC samples. We identified 2 distinct CRC subtypes with different mitophagy accumulation, of which subtype B had better prognosis and immune activity. The mitophagy score may be employed as a new and efficient clinical predictor in conjunction with other clinical indicators to predict the prognosis of CRC patients.ConclusionsWe systematically investigated the CRC heterogeneity with reference to mitophagy based on bioinformatics analyses, and the findings of this study might provide some guidance for future research into potential biomarkers for diagnosis and prognosis prediction of CRC patients.

Project description:BackgroundAccumulating evidence implies that autophagy plays a critical role in breast cancer development and progression. It is crucial to screen out autophagy-related encoding genes (ARGs) with prognostic value in breast cancer and reveal their biological properties in the aggressiveness of breast cancer.MethodsUnivariate and multivariate Cox proportional hazards analyses were used to identify a prognostic risk model of ARGs from The Cancer Genome Atlas (TCGA). Kaplan-Meier analysis, univariate and multivariate Cox regression analyses and receiver operating characteristic (ROC) curve analysis were performed to validate the risk model. Western blot and immunohistochemistry (IHC) were conducted to assess the expression of VPS35 (one of ARGs in risk model). CCK8, Colony formation assay, Transwell migration/invasion assays and autophagy flux assay were used to confirm biological function of VPS35 in breast cancer.ResultsIn this study, the prognostic risk model consisting of six ARGs (VPS35, TRIM21, PRKAB2, RUFY4, MAP1LC3A and LARP1) in breast cancer were identified. The risk model was further verified as a novel independent prognostic factor for breast cancer patients. We also clarified that vacuolar protein sorting-associated protein 35 (VPS35), one of ARGs in the risk model, was upregulated in breast cancer samples and cell lines. VPS35 overexpression was correlated with more aggressive phenotype of breast cancer and indicated worse prognosis in both progression-free survival and overall survival analyses. Meanwhile, VPS35 knockdown inhibited breast cancer cell proliferation, migration and invasion, suggesting that VPS35 promoted the progression of breast cancer. VPS35 silence also influenced autophagy process, indicating that VPS35 was essential for autophagy completion.ConclusionTaken together, the six ARGs risk model has a remarkably prognostic value for breast cancer. Among them, VPS35 might exert as a significant oncogenic and prognostic factor for breast cancer and could be a promising autophagy-related therapeutic target in clinical practice.

Project description:BackgroundThe role of N7-methyladenosine (m7G)-related genes in the progression and prognosis of gastric cancer (GC) remains unclear. This study aimed to explore prognostic biomarkers for GC based on m7G methylation regulators and to construct a prognostic risk model.MethodsRNA sequencing profiles with corresponding clinicopathological information associated with GC of which the histological type was stomach adenocarcinoma (STAD) were obtained from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO), respectively. A total of 29 m7G regulators were extracted from previous studies. According to the expression similarity of m7G regulators, the GC samples obtained from TCGA were further classified into 2 clusters demonstrating different overall survival (OS) rates and genetic heterogeneity, and the differentially expressed genes (DEGs) between these 2 clusters were defined as m7G-related genes. Univariate regression analysis and regression analysis were then used to obtain the prognostic m7G-related genes. The samples in TCGA and Genotype-Tissue Expression (GTEx) were used to verify the differential expression and prognostic value of these m7G-related genes contained in the prognostic model. Subsequently, the risk score was combined with other prognostic factors to develop a nomogram. The predictive ability of the nomogram was evaluated by the standard receiver operating characteristic (ROC) curve. Gene set enrichment analysis (GSEA) was used to identify activation pathways in both groups. Finally, the association between the prognostic model and the immune characteristics of GC were appraised.ResultsA prognostic model consisting of 11 m7G-related genes was constructed. GC patients in the high-risk group were shown to have a poor prognosis and this result was further demonstrated in each group. The risk model can be applied for patients with different clinical features. The results of GSEA showed that cell adhesion, cell junction, and focal adhesion were highly enriched in the high-risk group. In addition, we found that the expression of programmed cell death ligand 1 (PD-L1) was significantly elevated in the low-risk group, whereas programmed cell death ligand 2 (PD-L2) and tumor necrosis factor receptor superfamily member 4 (TNFRSF4) were overexpressed in the high-risk group.ConclusionsWe successfully built and verified a m7G relevant prognostic model for predicting prognosis and providing a new train of thought for improving the treatment of GC.