scholarly journals Screening of DNA Damage Repair Genes Involved in the Prognosis of Triple-Negative Breast Cancer Patients Based on Bioinformatics

Author(s):  
Nan Wang ◽  
Yuanting Gu ◽  
Jiangrui Chi ◽  
Xinwei Liu ◽  
Youyi Xiong ◽  
...  

Abstract Background: Triple-negative breast cancer (TNBC) is a special subtype of breast cancer with poor prognosis. DNA damage response (DDR) is one of the hallmarks of this cancer. However, the association of DDR genes with the prognosis of TNBC is still unclear.Methods: We identified differentially expressed genes (DEGs) between normal and TNBC samples from The Cancer Genome Atlas (TCGA). DDR genes were obtained from the Molecular Signatures Database (MSigDB) through six DDR gene sets. We then overlapped the DEGs with DDR genes. Based on univariate and LASSO Cox regression analyses, a prognostic model was constructed to predict overall survival (OS). Kaplan–Meier (K–M) analysis and receiver operating characteristic (ROC) curve were used to assess the performance of the prognostic model. Cox regression analysis was applied to identify independent prognostic factors in TNBC. The prognostic model was validated using an independent dataset. Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis were performed by using gene set enrichment analysis (GSEA). Single-sample gene set enrichment analysis (ssGSEA) was employed to estimate immune cells related to this prognostic model. Finally, we constructed a transcriptional factor (TF) network and a competing endogenous RNA (ceRNA) regulatory network.Results: 23 differentially expressed DDR genes were detected between TNBC and normal samples. The six-gene prognostic model we developed was shown to be related to OS in TNBC using univariate and LASSO Cox regression analyses. By drawing ROC curve and KM curve, we determined the effectiveness of the risk model. The prognostic value of the six-gene prognostic model was further validated using the GSE58812 dataset. The GSEA analysis indicated that the genes in the high-risk group were mainly correlated with leukocyte migration, cytokine interaction with cytokine receptors, oxidative phosphorylation, autoimmune diseases, and coagulation cascade. The mutation data revealed that the mutation frequency of the two groups was the same, while the mutated genes were different. The gene-TF regulatory network showed that Replication Factor C subunit 4 (RFC4) occupied the dominant position.Conclusion: We identified six gene markers related to DDR, which can predict prognosis and serve as an independent biomarker for TNBC patients.

2021 ◽  
Vol 12 ◽  
Author(s):  
Nan Wang ◽  
Yuanting Gu ◽  
Jiangrui Chi ◽  
Xinwei Liu ◽  
Youyi Xiong ◽  
...  

Background: Triple-negative breast cancer (TNBC) is a special subtype of breast cancer with poor prognosis. DNA damage response (DDR) is one of the hallmarks of this cancer. However, the association of DDR genes with the prognosis of TNBC is still unclear.Methods: We identified differentially expressed genes (DEGs) between normal and TNBC samples from The Cancer Genome Atlas (TCGA). DDR genes were obtained from the Molecular Signatures Database through six DDR gene sets. After the expression of six differential genes were verified by quantitative real-time polymerase chain reaction (qRT-PCR), we then overlapped the DEGs with DDR genes. Based on univariate and LASSO Cox regression analyses, a prognostic model was constructed to predict overall survival (OS). Kaplan–Meier analysis and receiver operating characteristic curve were used to assess the performance of the prognostic model. Cox regression analysis was applied to identify independent prognostic factors in TNBC. The Human Protein Atlas was used to study the immunohistochemical data of six DEGs. The prognostic model was validated using an independent dataset. Gene Ontology and the Kyoto Encyclopedia of Genes and Genomes analysis were performed by using gene set enrichment analysis (GSEA). Single-sample gene set enrichment analysis was employed to estimate immune cells related to this prognostic model. Finally, we constructed a transcriptional factor (TF) network and a competing endogenous RNA regulatory network.Results: Twenty-three differentially expressed DDR genes were detected between TNBC and normal samples. The six-gene prognostic model we developed was shown to be related to OS in TNBC using univariate and LASSO Cox regression analyses. All the six DEGs were identified as significantly up-regulated in the tumor samples compared to the normal samples in qRT-PCR. The GSEA analysis indicated that the genes in the high-risk group were mainly correlated with leukocyte migration, cytokine interaction, oxidative phosphorylation, autoimmune diseases, and coagulation cascade. The mutation data revealed the mutated genes were different. The gene-TF regulatory network showed that Replication Factor C subunit 4 occupied the dominant position.Conclusion: We identified six gene markers related to DDR, which can predict prognosis and serve as an independent biomarker for TNBC patients.


2021 ◽  
Vol 11 ◽  
Author(s):  
Libo Yang ◽  
Chunyan Li ◽  
Yang Qin ◽  
Guoying Zhang ◽  
Bin Zhao ◽  
...  

BackgroundBladder cancer (BC) is a molecular heterogeneous malignant tumor; the treatment strategies for advanced-stage patients were limited. Therefore, it is vital for improving the clinical outcome of BC patients to identify key biomarkers affecting prognosis. Ferroptosis is a newly discovered programmed cell death and plays a crucial role in the occurrence and progression of tumors. Ferroptosis-related genes (FRGs) can be promising candidate biomarkers in BC. The objective of our study was to construct a prognostic model to improve the prognosis prediction of BC.MethodsThe mRNA expression profiles and corresponding clinical data of bladder urothelial carcinoma (BLCA) patients were downloaded from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. FRGs were identified by downloading data from FerrDb. Differential analysis was performed to identify differentially expressed genes (DEGs) related to ferroptosis. Univariate and multivariate Cox regression analyses were conducted to establish a prognostic model in the TCGA cohort. BLCA patients from the GEO cohort were used for validation. Gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and single-sample gene set enrichment analysis (ssGSEA) were used to explore underlying mechanisms.ResultsNine genes (ALB, BID, FADS2, FANCD2, IFNG, MIOX, PLIN4, SCD, and SLC2A3) were identified to construct a prognostic model. Patients were classified into high-risk and low-risk groups according to the signature-based risk score. Receiver operating characteristic (ROC) and Kaplan–Meier (K–M) survival analysis confirmed the superior predictive performance of the novel survival model based on the nine-FRG signature. Multivariate Cox regression analyses showed that risk score was an independent risk factor associated with overall survival (OS). GO and KEGG enrichment analysis indicated that apart from ferroptosis-related pathways, immune-related pathways were significantly enriched. ssGSEA analysis indicated that the immune status was different between the two risk groups.ConclusionThe results of our study indicated that a novel prognostic model based on the nine-FRG signature can be used for prognostic prediction in BC patients. FRGs are potential prognostic biomarkers and therapeutic targets.


2021 ◽  
Author(s):  
Feng Wang ◽  
Cheng Chen ◽  
Wei-Peng Chen ◽  
Zu-Ling Li ◽  
Hui Cheng

Abstract Background Ferroptosis is a mode of regulated cell death that depends on iron, plays pivotal roles in regulating various biological process in human cancers. However, the role of ferroptosis in Gastric cancer (GC) remains unclear. Methods A total of 2721 differentially expressed genes (DEGs) were filtered base on The Cancer Genome Atlas (TCGA) (n = 375) dataset. Gene modules were identified based on co-expression network analysis (WGCNA). Functional analysis was performed to explore the biological function. Lasso-penalized and univariate Cox regression (UCR) analysis, survival genes were screened out to construct a prognostic model, which validated by the GSE43292 dataset. Gene set enrichment analysis (GSEA) for prognostic index was performed. Finally, the correlations of ferroptosis and immune cells were assessed through the TIMER database. Results Compared to normal specimens, 1063 highly upregulated and 1658 downregulated genes respectively and their normal counterparts in GC specimens were screened. WGCNA analysis was used and identified 7 modules, of which, blue module with the most significant enrichment result was selected. By taking intersections of blue module and differentially expressed ferroptosis-related genes (DEFRGs), we got 23 common genes. Functional analysis was performed to explore the biological function of the interested genes, and with the consequences Lasso-penalized and univariate Cox regression (UCR) analysis, survival genes were screened out to construct a prognostic model based on 3 genes (SLC1A5, ANGPTL4, and CGAS), which could play a role in predicting the survival of GC patients. UCR and multivariate Cox regression (MCR) analysis revealed that the prognostic index could be used as independent prognostic indicators and validated using another GSE84437 dataset. Notably, patients in high-risk groups had higher levels of higher mutation frequencies such as TTN and TP53.Mechanistically. Gene set enrichment analysis (GSEA) unveiled several significant and pathways involved in GC. TIMER analysis demonstrated that risk score strongly correlated with Macrophage and CD4 + T cells infiltration. In addition, high- and low-risk group illustrated different distributions in different immune status. Conclusions In this study, a novel FRGs signature was built. It could accurately predict GC prognosis and pave the new way for diagnosis and therapy strategy. May reflect the status of tumor immune microenvironment.


2021 ◽  
Author(s):  
Zhian Ling ◽  
Yuting Liang ◽  
Suping Wei ◽  
Yuanming Chen ◽  
Jinmin Zhao

Abstract Background N6-methylandenosine (m6A) methylation is one of the most common methylation modifications in RNA. At present, a large number of studies have found that m6A methylation can regulate the occurrence and development of tumors by modifying mRNA. However, it is still unclear how m6A modifies Long non-coding RNA (lncRNA) that regulates mRNA expression by interacting with miRNA to affect the occurrence and development of osteosarcoma(OS). Therefore, exploring the lncRNAs related to m6A methylation and identifying lncRNAs that have both prognostic effects and immune functions are things that need to be solved urgently. Methods The published gene expression data of OS and complete clinical annotation files were obtained from the TARGET database. LncRNAs with P <0.001 from the results of Pearson correlation coefficient analysis as m6A-related lncRNAs were screened. Single-factor Cox regression analysis was used to screening prognostic- related lncRNA combined with the clinical information of patients and constructed a prognostic model based on lasso regression analysis. Then we explored the differences in survival and immune function of different subtypes that be obtained using the Consensus Cluster. The enrichment of differential genes between high and low risk groups in the KEGG pathway is achieved through Gene set enrichment analysis(GSEA). Results We obtained 706 lncRNAs in the TARGET database. Consensus clustering method were used to divide patients with OS into subgroups based on the expression of 26 prognostic-related lncRNAs. Through Kaplan-Meier survival analysis, there are significant differences between the two subgroups. The average immune score (P = 0.02), stromal score(P =0.027), and estimate score༈P = 0.015༉were higher in cluster 1 than in cluster 2. We found that compared with cluster 2, SIGLEC15, HAVCR2, LAG3, and PDCD1 were highly expressed in cluster 1.We obtain a prognostic model by lasso regression analysis. In the training group and the text group, the OS curve showed that patients in the high-risk group had a poorer prognosis than those in the low-risk group. In the training set, univariate Cox regression analysis and multivariate Cox regression analysis showed that the risk score was correlated with the prognosis of OS patients. In the high-risk group, the Linoleic acid metabolism and the Glycine, serine and threonine metabolism pathway were mainly involved by Gene Set Enrichment analysis. The abundance of Mast cells activated (P ≦0.024) and T cells CD4 (P ≦0.0044) naive were positively association the risk score. Conclusions This study clarified the important role of m6A-related lncRNAs in the prognosis and immune microenvironment of patients with OS, and indicate that m6A-related prognostic lncRNA signals may provide new targets for the diagnosis and treatment of OS.


2021 ◽  
Vol 44 (3) ◽  
pp. E32-44
Author(s):  
Jia Shen ◽  
Ming Shu ◽  
Shujie Xie ◽  
Jia Yan ◽  
Kaile Pan ◽  
...  

Purpose: This study aimed to screen hepatitis B virus (HBV)-associated hepatocellular carcinoma (HCC)-related feature ribonucleic acids (RNAs) and to establish a prognostic model. Methods: The transcriptome expression data of HBV-associated HCC were downloaded from The Cancer Genome Atlas (TCGA) database and Gene Expression Omnibus database. Differential RNAs between HBV-associated HCC and normal controls were identified by a meta-analysis of TCGA, GSE55092 and GSE121248. Weighted gene co-expression network analysis was performed to identify key RNAs and modules. A prognostic score model was established using TCGA as a training set by Cox regression analysis and was validated in E-TABM-36 dataset. Additionally, independent prognostic clinical factors were screened, and the function of lncRNAs waspredicted through Gene Set Enrichment Analysis. Results: A total of 710 consistent differential RNAs between HBV-associated HCC and normal controls were obtained, including five lncRNAs and 705 mRNAs. An optimized combination of six differential RNAs (DSCR4, DBH, ECM1, GDAP1, MATR3 and RFC4) was selected and a prognostic score model was constructed. Kaplan-Meier analysis demonstrated that the prognosis of the high-risk and low-risk groups separated by this model was significantly different in the training set and the validation set. Gene Set Enrichment Analysis showed that the co-expression genes of DSCR4 were significantly correlated with neuroactive ligand receptor interactionpathway. Conclusion: A prognostic model based on DSCR4, DBH, ECM1, GDAP1, MATR3 and RFC4 was developed that can accurately predict the prognosis of patients with HBV-associated HCC. These genes, as well as histologic grade, may serve as independent prognostic factors in HBV-associated HCC.


2021 ◽  
Author(s):  
Yanjia Hu ◽  
Jing Zhang ◽  
Jing Chen

Abstract Background Hypoxia-related long non-coding RNAs (lncRNAs) have been proven to play a role in multiple cancers and can serve as prognostic markers. Lower-grade gliomas (LGGs) are characterized by large heterogeneity. Methods This study aimed to construct a hypoxia-related lncRNA signature for predicting the prognosis of LGG patients. Transcriptome and clinical data of LGG patients were obtained from The Cancer Genome Atlas (TCGA) and the Chinese Glioma Genome Atlas (CGGA). LGG cohort in TCGA was chosen as training set and LGG cohorts in CGGA served as validation sets. A prognostic signature consisting of fourteen hypoxia-related lncRNAs was constructed using univariate and LASSO Cox regression. A risk score formula involving the fourteen lncRNAs was developed to calculate the risk score and patients were classified into high- and low-risk groups based on cutoff. Kaplan-Meier survival analysis was used to compare the survival between two groups. Cox regression analysis was used to determine whether risk score was an independent prognostic factor. A nomogram was then constructed based on independent prognostic factors and assessed by C-index and calibration plot. Gene set enrichment analysis and immune cell infiltration analysis were performed to uncover further mechanisms of this lncRNA signature. Results LGG patients with high risk had poorer prognosis than those with low risk in both training and validation sets. Recipient operating characteristic curves showed good performance of the prognostic signature. Univariate and multivariate Cox regression confirmed that the established lncRNA signature was an independent prognostic factor. C-index and calibration plots showed good predictive performance of nomogram. Gene set enrichment analysis showed that genes in the high-risk group were enriched in apoptosis, cell adhesion, pathways in cancer, hypoxia etc. Immune cells were higher in high-risk group. Conclusion The present study showed the value of the 14-lncRNA signature in predicting survival of LGGs and these 14 lncRNAs could be further investigated to reveal more mechanisms involved in gliomas.


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) &lt; 1), and HLA-F was risky (HR &gt; 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.


2021 ◽  
Vol 11 ◽  
Author(s):  
Wenxing Qin ◽  
Feng Qi ◽  
Jia Li ◽  
Ping Li ◽  
Yuan-Sheng Zang

The objective of this study was to construct a competitive endogenous RNA (ceRNA) regulatory network using differentially expressed long noncoding RNAs (lncRNAs), microRNAs (miRNAs), and mRNAs in patients with triple-negative breast cancer (TNBC) and to construct a prognostic model for predicting overall survival (OS) in patients with TNBC. Differentially expressed lncRNAs, miRNAs, and mRNAs in TNBC patients from the TCGA and Metabric databases were examined. A prognostic model based on prognostic scores (PSs) was established for predicting OS in TNBC patients, and the performance of the model was assessed by a recipient that operated on a distinctive curve. A total of 874 differentially expressed RNAs (DERs) were screened, among which 6 lncRNAs, 295 miRNAs and 573 mRNAs were utilized to construct targeted and coexpression ceRNA regulatory networks. Eight differentially expressed genes (DEGs) associated with survival prognosis, DBX2, MYH7, TARDBP, POU4F1, ABCB11, LHFPL5, TRHDE and TIMP4, were identified by multivariate Cox regression and then used to establish a prognostic model. Our study shows that the ceRNA network has a critical role in maintaining the aggressiveness of TNBC and provides comprehensive molecular-level insight for predicting individual mortality hazards for TNBC patients. Our data suggest that these prognostic mRNAs from the ceRNA network are promising therapeutic targets for clinical intervention.


Author(s):  
Bo Xiao ◽  
Liyan Liu ◽  
Zhuoyuan Chen ◽  
Aoyu Li ◽  
Pingxiao Wang ◽  
...  

Melanoma is the most common cancer of the skin, associated with a worse prognosis and distant metastasis. Epithelial–mesenchymal transition (EMT) is a reversible cellular biological process that plays significant roles in diverse tumor functions, and it is modulated by specific genes and transcription factors. The relevance of EMT-related lncRNAs in melanoma has not been determined. Therefore, RNA expression data and clinical features were collected from the TCGA database (N = 447). Melanoma samples were randomly assigned into the training (315) and testing sets (132). An EMT-related lncRNA signature was constructed via comprehensive analyses of lncRNA expression level and corresponding clinical data. The Kaplan-Meier analysis showed significant differences in overall survival in patients with melanoma in the low and high-risk groups in two sets. Receiver operating characteristic (ROC) curves were used to measure the performance of the model. Cox regression analysis indicated that the risk score was an independent prognostic factor in two sets. Besides, a nomogram was constructed based on the independent variables. Gene Set Enrichment Analysis (GSEA) was applied to evaluate the potential biological functions in the two risk groups. Furthermore, the melanoma microenvironment was evaluated using ESTIMATE and CIBERSORT algorithms in the risk groups. This study indicates that EMT-related lncRNAs can function as potential independent prognostic biomarkers for melanoma survival.


Sign in / Sign up

Export Citation Format

Share Document