scholarly journals Identification and External Validation of a Transcription Factor-Related Prognostic Signature in Pediatric Neuroblastoma

2021 ◽  
Vol 2021 ◽  
pp. 1-10
Author(s):  
Rujia Wang ◽  
Qian Wang

Background. Neuroblastoma is a common solid tumor originating from the sympathetic nervous system, commonly found in children, and it is one of the leading causes of tumor-related deaths in children. In addition to pathological features, molecular-level features, such as how much gene expression is present and the mutational profile, may provide useful information for the precise treatment of neuroblastoma. Transcription factors (TFs) play an important regulatory role in all aspects of cellular life activities. But there are currently no studies on transcription factor-based biomarkers of neuroblastoma prognosis, and this study is much needed. Methods. We downloaded RNA transcriptome data and clinical data from the TARGET database to construct a prognostic model. The prognostic model was constructed by using univariate Cox analysis, LASSO, and multivariate Cox regression. We divided the patients into low-risk and high-risk groups using the median value of the risk score as the cut-off. Then, we validated the prognostic model with the dataset GSE49710. Results. We constructed a prognostic model consisting of eight genes (SATB1, ZNF564, SOX14, EN1, IKZF2, SLC2A4RG, FOXJ2, and ZNF521). Patients in the high-risk group had a lower survival rate than those in the low-risk group. The area under the 3-year ROC curve of the model reached 0.825, suggesting a good predictive efficacy. We performed target gene prediction for the eight transcription factors in the model using six online databases and found that TUT1 may be a target gene for transcription factor EN1 and is associated with immune infiltration. Conclusion. This prognostic model consisting of eight transcription factor-associated genes demonstrated reliable predictive efficacy. This prediction model may provide new potential targets for the treatment of neuroblastoma and personalized monitoring of neuroblastoma patients with high and low risk.

Author(s):  
Peng Gu ◽  
Lei Zhang ◽  
Ruitao Wang ◽  
Wentao Ding ◽  
Wei Wang ◽  
...  

Background: Female breast cancer is currently the most frequently diagnosed cancer in the world. This study aimed to develop and validate a novel hypoxia-related long noncoding RNA (HRL) prognostic model for predicting the overall survival (OS) of patients with breast cancer.Methods: The gene expression profiles were downloaded from The Cancer Genome Atlas (TCGA) database. A total of 200 hypoxia-related mRNAs were obtained from the Molecular Signatures Database. The co-expression analysis between differentially expressed hypoxia-related mRNAs and lncRNAs based on Spearman’s rank correlation was performed to screen out 166 HRLs. Based on univariate Cox regression and least absolute shrinkage and selection operator Cox regression analysis in the training set, we filtered out 12 optimal prognostic hypoxia-related lncRNAs (PHRLs) to develop a prognostic model. Kaplan–Meier survival analysis, receiver operating characteristic curves, area under the curve, and univariate and multivariate Cox regression analyses were used to test the predictive ability of the risk model in the training, testing, and total sets.Results: A 12-HRL prognostic model was developed to predict the survival outcome of patients with breast cancer. Patients in the high-risk group had significantly shorter median OS, DFS (disease-free survival), and predicted lower chemosensitivity (paclitaxel, docetaxel) compared with those in the low-risk group. Also, the risk score based on the expression of the 12 HRLs acted as an independent prognostic factor. The immune cell infiltration analysis revealed that the immune scores of patients in the high-risk group were lower than those of the patients in the low-risk group. RT-qPCR assays were conducted to verify the expression of the 12 PHRLs in breast cancer tissues and cell lines.Conclusion: Our study uncovered dozens of potential prognostic biomarkers and therapeutic targets related to the hypoxia signaling pathway in breast cancer.


PeerJ ◽  
2019 ◽  
Vol 7 ◽  
pp. e8128 ◽  
Author(s):  
Cheng Yue ◽  
Hongtao Ma ◽  
Yubai Zhou

Background Lung cancer has the highest morbidity and mortality worldwide, and lung adenocarcinoma (LADC) is the most common pathological subtype. Accumulating evidence suggests the tumor microenvironment (TME) is correlated with the tumor progress and the patient’s outcome. As the major components of TME, the tumor-infiltrated immune cells and stromal cells have attracted more and more attention. In this study, differentially expressed immune and stromal signature genes were used to construct a TME-related prognostic model for predicting the outcomes of LADC patients. Methods The expression profiles of LADC samples with clinical information were obtained from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO). The differentially expressed genes (DEGs) related to the TME of LADC were identified using TCGA dataset by Wilcoxon rank sum test. The prognostic effects of TME-related DEGs were analyzed using univariate Cox regression. Then, the least absolute shrinkage and selection operator (LASSO) regression was performed to reduce the overfit and the number of genes for further analysis. Next, the prognostic model was constructed by step multivariate Cox regression and risk score of each sample was calculated. Then, survival and Receiver Operating Characteristic (ROC) analyses were conducted to validate the model using TCGA and GEO datasets, respectively. The Kyoto Encyclopedia of Genes and Genomes analysis of gene signature was performed using Gene Set Enrichment Analysis (GSEA). Finally, the overall immune status, tumor purity and the expression profiles of HLA genes of high- and low-risk samples was further analyzed to reveal the potential mechanisms of prognostic effects of the model. Results A total of 93 TME-related DEGs were identified, of which 23 DEGs were up-regulated and 70 DEGs were down-regulated. The univariate cox analysis indicated that 23 DEGs has the prognostic effects, the hazard ratio ranged from 0.65 to 1.25 (p < 0.05). Then, seven genes were screened out from the 23 DEGs by LASSO regression method and were further analyzed by step multivariate Cox regression. Finally, a three-gene (ADAM12, Bruton Tyrosine Kinase (BTK), ERG) signature was constructed, and ADAM12, BTK can be used as independent prognostic factors. The three-gene signature well stratified the LADC patients in both training (TCGA) and testing (GEO) datasets as high-risk and low-risk groups, the 3-year area under curve (AUC) of ROC curves of three GEO sets were 0.718 (GSE3141), 0.646 (GSE30219) and 0.643 (GSE50081). The GSEA analysis indicated that highly expressed ADAM12, BTK, ERG mainly correlated with the activation of pathways involving in focal adhesion, immune regulation. The immune analysis indicated that the low-risk group has more immune activities and higher expression of HLA genes than that of the high-risk group. In sum, we identified and constructed a three TME-related DEGs signature, which could be used to predict the prognosis of LADC patients.


2020 ◽  
Vol 2020 ◽  
pp. 1-12
Author(s):  
Dakui Luo ◽  
Zezhi Shan ◽  
Qi Liu ◽  
Sanjun Cai ◽  
Qingguo Li ◽  
...  

A metabolic disorder is considered one of the hallmarks of cancer. Multiple differentially expressed metabolic genes have been identified in colon cancer (CC), and their biological functions and prognostic values have been well explored. The purpose of the present study was to establish a metabolic signature to optimize the prognostic prediction in CC. The related data were downloaded from The Cancer Genome Atlas (TCGA), Genotype-Tissue Expression (GTEx) database, and Gene Expression Omnibus (GEO) combined with GSE39582 set, GSE17538 set, GSE33113 set, and GSE37892 set. The differentially expressed metabolic genes were selected for univariate Cox regression and lasso Cox regression analysis using TCGA and GTEx datasets. Finally, a seventeen-gene metabolic signature was developed to divide patients into a high-risk group and a low-risk group. Patients in the high-risk group presented poorer prognosis compared to the low-risk group in both TCGA and GEO datasets. Moreover, gene set enrichment analyses demonstrated multiple significantly enriched metabolism-related pathways. To sum up, our study described a novel seventeen-gene metabolic signature for prognostic prediction of colon cancer.


2020 ◽  
Vol 13 (1) ◽  
Author(s):  
Yinglian Pan ◽  
Li Ping Jia ◽  
Yuzhu Liu ◽  
Yiyu Han ◽  
Qian Li ◽  
...  

Abstract Background In this study we aimed to identify a prognostic signature in BRCA1/2 mutations to predict disease progression and the efficiency of chemotherapy ovarian cancer (OV), the second most common cause of death from gynecologic cancer in women worldwide. Methods Univariate Cox proportional-hazards and multivariate Cox regression analyses were used to identifying prognostic factors from data obtained from The Cancer Genome Atlas (TCGA) database. The area under the curve of the receiver operating characteristic curve was assessed, and the sensitivity and specificity of the prediction model were determined. Results A signature consisting of two long noncoding RNAs(lncRNAs), Z98885.2 and AC011601.1, was selected as the basis for classifying patients into high and low-risk groups (median survival: 7.2 years vs. 2.3 years). The three-year overall survival (OS) rates for the high- and low-risk group were approximately 38 and 100%, respectively. Chemotherapy treatment survival rates indicated that the high-risk group had significantly lower OS rates with adjuvant chemotherapy than the low-risk group. The one-, three-, and five-year OS were 100, 40, and 15% respectively in the high-risk group. The survival rate of the high-risk group declined rapidly after 2 years of OV chemotherapy treatment. Multivariate Cox regression associated with other traditional clinical factors showed that the 2-lncRNA model could be used as an independent OV prognostic factor. Analyses of data from the Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) indicated that these signatures are pivotal to cancer development. Conclusion In conclusion, Z98885.2 and AC011601.1 comprise a novel prognostic signature for OV patients with BRCA1/2 mutations, and can be used to predict prognosis and the efficiency of chemotherapy.


2020 ◽  
Author(s):  
Li Liu ◽  
She Tian ◽  
Zhu Li ◽  
Yongjun Gong ◽  
Hao Zhang

Abstract Background : Hepatocellular carcinoma (HCC) is one of the most common clinical malignant tumors, resulting in high mortality and poor prognosis. Studies have found that LncRNA plays an important role in the onset, metastasis and recurrence of hepatocellular carcinoma. The immune system plays a vital role in the development, progression, metastasis and recurrence of cancer. Therefore, immune-related lncRNA can be used as a novel biomarker to predict the prognosis of hepatocellular carcinoma. Methods : The transcriptome data and clinical data of HCC patients were obtained by using The Cancer Genome Atlas-Liver Hepatocellular Carcinoma (TCGA‑LIHC), and immune-related genes were extracted from the Molecular Signatures Database (IMMUNE RESPONSE M19817 and IMMUNE SYSTEM PROCESS M13664). By constructing the co-expression network and Cox regression analysis, 13 immune-lncRNAs was identified to predict the prognosis of HCC patients. Patients were divided into high risk group and low risk group by using the risk score formula, and the difference in overall survival (OS) between the two groups was reflected by Kaplan-Meier survival curve. The time - dependent receiver operating characteristics (ROC) analysis and principal component analysis (PCA) were used to evaluate 13 immune -lncRNAs signature. Results : Through TCGA - LIHC extracted from 343 cases of patients with hepatocellular carcinoma RNA - Seq data and clinical data, 331 immune-related genes were extracted from the Molecular Signatures Database , co-expression networks and Cox regression analysis were constructed, 13 immune-lncRNAs signature was identified as biomarkers to predict the prognosis of patients. At the same time using the risk score median divided the patients into high risk and low risk groups, and through the Kaplan-Meier survival curve analysis found that high-risk group of patients' overall survival (OS) less low risk group of patients. The AUC value of the ROC curve is 0.828, and principal component analysis (PCA) results showed that patients could be clearly divided into two parts by immune-lncRNAs, which provided evidence for the use of 13 immune-lncRNAs signature as prognostic markers. Conclusion : Our study identified 13 immune-lncRNAs signature that can effectively predict the prognosis of HCC patients, which may be a new prognostic indicator for predicting clinical outcomes.


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.


2020 ◽  
Vol 40 (6) ◽  
Author(s):  
Xin Zhu ◽  
Qian Zhao ◽  
Xiaoyu Su ◽  
Jinming Ke ◽  
Yunyun Yi ◽  
...  

Abstract The identification of effective signatures is crucial to predict the prognosis of acute myeloid leukemia (AML). The investigation aimed to identify a new signature for AML prognostic prediction by using the three-gene expression (octamer-binding transcription factor 4 (OCT4), POU domain type 5 transcription factor 1B (POU5F1B) and B-cell-specific Moloney murine leukemia virus integration site-1 pseudogene 1 (BMI1P1). The expressions of genes were obtained from our previous study. Only the specimens in which three genes were all expressed were included in this research. A three-gene signature was constructed by the multivariate Cox regression analyses to divide patients into high-risk and low-risk groups. Receiver operating characteristic (ROC) analysis of the three-gene signature (area under ROC curve (AUC) = 0.901, 95% CI: 0.821–0.981, P&lt;0.001) indicated that it was a more valuable signature for distinguishing between patients and controls than any of the three genes. Moreover, white blood cells (WBCs, P=0.004), platelets (PLTs, P=0.017), percentage of blasts in bone marrow (BM) (P=0.011) and complete remission (CR, P=0.027) had significant differences between two groups. Furthermore, high-risk group had shorter leukemia-free survival (LFS) and overall survival (OS) than low-risk group (P=0.026; P=0.006), and the three-gene signature was a prognostic factor. Our three-gene signature for prognosis prediction in AML may serve as a prognostic biomarker.


2021 ◽  
Vol 12 ◽  
Author(s):  
Qian Zhang ◽  
Liping Lv ◽  
Ping Ma ◽  
Yangyang Zhang ◽  
Jiang Deng ◽  
...  

BackgroundPancreatic adenocarcinoma (PAAD) spreads quickly and has a poor prognosis. Autophagy research on PAAD could reveal new biomarkers and targets for diagnosis and treatment.MethodsAutophagy-related genes were translated into autophagy-related gene pairs, and univariate Cox regression was performed to obtain overall survival (OS)-related IRGPs (P&lt;0.001). LASSO Cox regression analyses were performed to construct an autophagy-related gene pair (ARGP) model for predicting OS. The Cancer Genome Atlas (TCGA)-PAAD cohort was set as the training group for model construction. The model predictive value was validated in multiple external datasets. Receiver operating characteristic (ROC) curves were used to evaluate model performance. Tumor microenvironments and immune infiltration were compared between low- and high-risk groups with ESTIMATE and CIBERSORT. Differentially expressed genes (DEGs) between the groups were further analyzed by Gene Ontology biological process (GO-BP) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses and used to identify potential small-molecule compounds in L1000FWD.ResultsRisk scores were calculated as follows: ATG4B|CHMP4C×(-0.31) + CHMP2B|MAP1LC3B×(0.30) + CHMP6|RIPK2 ×(-0.33) + LRSAM1|TRIM5×(-0.26) + MAP1LC3A|PAFAH1B2×(-0.15) + MAP1LC3A|TRIM21×(-0.08) + MET|MFN2×(0.38) + MET|MTDH×(0.47) + RASIP1|TRIM5×(-0.23) + RB1CC1|TPCN1×(0.22). OS was significantly shorter in the high-risk group than the low-risk group in each PAAD cohort. The ESTIMATE analysis showed no difference in stromal scores but a significant difference in immune scores (p=0.0045) and ESTIMATE scores (p=0.014) between the groups. CIBERSORT analysis showed higher naive B cell, Treg cell, CD8 T cell, and plasma cell levels in the low-risk group and higher M1 and M2 macrophage levels in the high-risk group. In addition, the results showed that naive B cells (r=-0.32, p&lt;0.001), Treg cells (r=-0.31, p&lt;0.001), CD8 T cells (r=-0.24, p=0.0092), and plasma cells (r=-0.2, p&lt;0.026) were statistically correlated with the ARGP risk score. The top 3 enriched GO-BPs were signal release, regulation of transsynaptic signaling, and modulation of chemical synaptic transmission, and the top 3 enriched KEGG pathways were the insulin secretion, dopaminergic synapse, and NF-kappa B signaling pathways. Several potential small-molecule compounds targeting ARGs were also identified.ConclusionOur results demonstrate that the ARGP-based model may be a promising prognostic indicator for identifying drug targets in patients with PAAD.


2021 ◽  
Author(s):  
Chen-jie Qiu ◽  
Xue-bing Wang ◽  
Zi-ruo Zheng ◽  
Chao-zhi Yang ◽  
Kai Lin ◽  
...  

Abstract Background: The purpose of this study was to identify ferroptosis-related genes (FRGs) associated with the prognosis of pancreatic cancer and to construct a prognostic model based on FRGs. Methods: Based on pancreatic cancer data obtained from The Cancer Genome Atlas database, we established the prognostic model from 232 FRGs. A nomogram was constructed by combining the prognostic model and clinicopathological features. Gene Expression Omnibus datasets and tissue samples obtained from our center were utilized to validate the model. Relationship between risk score and immune cell infiltration was explored by CIBERSORT and TIMER.Results: The prognostic model was established based on four FRGs (ENPP2, ATG4D, SLC2A1 and MAP3K5) and can be an independent risk factor in pancreatic cancer (HR 1.648, 95% CI 1.335-2.035, p < 0.001). Based on the median risk score, patients were divided into a high-risk group and a low-risk group. The prognosis of the low-risk group was significantly better than that of the high-risk group. In the high-risk group, patients treated with chemotherapy had a better prognosis. The nomogram showed that the model was the most important element. Gene set enrichment analysis identified three key pathways, namely, TGFβ signaling, HIF signaling pathway and adherens junction. The prognostic model can also affect the immune cell infiltration, such as macrophages M0, M1, CD4+T cell and CD8+T cell. Conclusion: A ferroptosis-related prognostic model can be employed to predict the prognosis of pancreatic cancer. Ferroptosis can be an important marker and immunotherapy can be a potential therapeutic target for pancreatic cancer.


2021 ◽  
Vol 11 ◽  
Author(s):  
Xiaoxia Tong ◽  
Xiaofei Qu ◽  
Mengyun Wang

BackgroundCutaneous melanoma (CM) is one of the most aggressive cancers with highly metastatic ability. To make things worse, there are limited effective therapies to treat advanced CM. Our study aimed to investigate new biomarkers for CM prognosis and establish a novel risk score system in CM.MethodsGene expression data of CM from Gene Expression Omnibus (GEO) datasets were downloaded and analyzed to identify differentially expressed genes (DEGs). The overlapped DEGs were then verified for prognosis analysis by univariate and multivariate COX regression in The Cancer Genome Atlas (TCGA) datasets. Based on the gene signature of multiple survival associated DEGs, a risk score model was established, and its prognostic and predictive role was estimated through Kaplan-Meier (K-M) analysis and log-rank test. Furthermore, the correlations between prognosis related genes expression and immune infiltrates were analyzed via Tumor Immune Estimation Resource (TIMER) site.ResultsA total of 103 DEGs were obtained based on GEO cohorts, and four genes were verified in TCGA datasets. Subsequently, four genes (ADAMDEC1, GNLY, HSPA13, and TRIM29) model was developed by univariate and multivariate Cox regression analyses. The K-M plots showed that the high-risk group was associated with shortened survival than that in the low-risk group (P &lt; 0.0001). Multivariate analysis suggested that the model was an independent prognostic factor (high-risk vs. low-risk, HR= 2.06, P &lt; 0.001). Meanwhile, the high-risk group was prone to have larger breslow depth (P&lt; 0.001) and ulceration (P&lt; 0.001).ConclusionsThe four-gene risk score model functions well in predicting the prognosis and treatment response in CM and will be useful for guiding therapeutic strategies for CM patients. Additional clinical trials are needed to verify our findings.


Sign in / Sign up

Export Citation Format

Share Document