Identification of a Five-Gene Signature Derived From MYCN Amplification and Establishment of a Nomogram for Predicting the Prognosis of Neuroblastoma

Identification and validation of a novel five-gene signature in high-risk MYCN-not-amplified neuroblastoma

OBJECTIVE

High-risk neuroblastoma patients often have poor outcomes despite multi-treatment options. The risk stratification of high-risk MYCN-not-amplified (HR-MYCN-NA) patients remains difficult. This study aims to identify a gene set signature that can help further stratify HR-MYCN-NA patients for a potential personalized therapeutic strategy.

METHODS

Three microarrays and one single-cell RNA sequence dataset were acquired and analyzed. Firstly, the prognostic-related genes (PRGs) in HR-MYCN-NA tumor cells were identified using TARGET-NB and GSE137804 datasets. Then, the prognostic model was established by LASSO-Cox regression, and verified in external cohort (GSE49710, GSE45547). Moreover, a time-dependent receiver operating characteristic curve (ROC) and area under the ROC (AUC) was used to assess survival prediction. A nomogram was established to predict the 1-, 3- and 5-year overall survival (OS) of HR-MYCN-NA patients.

RESULTS

In the training set, a five-PRGs signature, which include GAL, GFRA3, MARCKS, PSMD13, and ZNHIT3 genes, was identified and successfully stratified HR-MYCN-NA patients into ultra-high risk (UHR) and high-risk (HR) subtypes (HR = 4.29, P < 0.001). ROC curve analysis confirmed its predictive power (AUC = 0.74-0.82), suggesting a good predictive efficacy. Consistently, high-risk scores also predicted worse OS (HR = 2, P = 0.033) in the external validation dataset (AUC = 0.67-0.71). Moreover, the overall C-index of the nomogram was 0.75 (P < 0.001), which indicated good agreement between the observed and predicted survival rates. Further integrating the five PRGs signature with clinical factors, these 5 gene signature (HR = 4.45, P < 0.001) and tumor grade (HR = 4.15, P = 0.02) were found to be independent prognostic factors for HR-MYCN-NA patients.

CONCLUSION

The novel five PRGs signature could well predict the survival of HR-MYCN-NA patients, which may provide constructive information for these subsets.

High expression of NXPH4 correlates with poor prognosis, metabolic reprogramming, and immune infiltration in colon adenocarcinoma

Background

Colon adenocarcinoma (COAD) is a prevalent gastrointestinal malignant disease with high mortality rate, and identification of novel prognostic biomarkers and therapeutic targets is urgently needed. Although neurexophilin 4 (NXPH4) has been investigated in several tumors, its role in COAD remains unclear. The aim of this study was to explore the prognostic value and potential functions of NXPH4 in COAD.

Methods

The expression of NXPH4 in COAD were analyzed using The Cancer Genome Atlas (TCGA) and datasets from the Gene Expression Omnibus (GEO) database. The prognostic value of NXPH4 was determined using Kaplan-Meier analysis and Cox regression analysis. To investigate the possible mechanism underlying the role of NXPH4 in COAD, Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and gene set enrichment analysis (GSEA) were employed. The correlation between NXPH4 expression and immune cell infiltration levels was examined thorough single-sample gene set enrichment analysis (ssGSEA). Furthermore, the competing endogenous RNA (ceRNA) regulatory network that may be involved in NXPH4 in COAD was predicted and constructed through a variety of databases.

Results

NXPH4 expression was significantly higher in COAD tissue compared with normal colon tissues. Meanwhile, high expression of NXPH4 was associated with poor prognosis in COAD patients. GO-KEGG and GSEA analyses indicated that NXPH4 was associated with glycolysis and hypoxia pathway, and may promote COAD progression and metastasis by modulating metabolic reprogramming. ssGSEA analysis demonstrated that NXPH4 expression also associated with immune infiltration. Furthermore, we identified various microRNAs (miRNAs) and long noncoding RNAs (lncRNAs) as upstream regulators of NXPH4 in COAD.

Conclusions

The present study revealed that high expression of NXPH4 is associated with tumor progression, metabolic reprogramming, and immune infiltration. These findings suggest that NXPH4 could serve as a reliable prognostic biomarker and a promising therapeutic target in COAD.

Development and validation of a novel stemness-related prognostic model for neuroblastoma using integrated machine learning and bioinformatics analyses

Background

Neuroblastoma (NB) is a common solid tumor in children, with a dismal prognosis in high-risk cases. Despite advancements in NB treatment, the clinical need for precise prognostic models remains critical, particularly to address the heterogeneity of cancer stemness which plays a pivotal role in tumor aggressiveness and patient outcomes. By utilizing machine learning (ML) techniques, we aimed to explore the cancer stemness features in NB and identify stemness-related hub genes for future investigation and potential targeted therapy.

Methods

The public dataset GSE49710 was employed as the training set for acquire gene expression data and NB sample information, including age, stage, and MYCN amplification status and survival. The messenger RNA (mRNA) expression-based stemness index (mRNAsi) was calculated and patients were grouped according to their mRNAsi value. Stemness-related hub genes were identified from the differentially expressed genes (DEGs) to construct a gene signature. This was followed by evaluating the relationship between cancer stemness and the NB immune microenvironment, and the development of a predictive nomogram. We assessed the prognostic outcomes including overall survival (OS) and event-free survival, employing machine learning methods to measure predictive accuracy through concordance indices and validation in an independent cohort E-MTAB-8248.

Results

Based on mRNAsi, we categorized NB patients into two groups to explore the association between varying levels of stemness and their clinical outcomes. High mRNAsi was linked to the advanced International Neuroblastoma Staging System (INSS) stage, amplified MYCN, and elder age. High mRNAsi patients had a significantly poorer prognosis than low mRNAsi cases. According to the multivariate Cox analysis, the mRNAsi was an independent risk factor of prognosis in NB patients. After least absolute shrinkage and selection operator (LASSO) regression analysis, four key genes (ERCC6L, DUXAP10, NCAN, DIRAS3) most related to mRNAsi scores were discovered and a risk model was built. Our model demonstrated a significant prognostic capacity with hazard ratios (HR) ranging from 18.96 to 41.20, P values below 0.0001, and area under the receiver operating characteristic curve (AUC) values of 0.918 in the training set, suggesting high predictive accuracy which was further confirmed by external verification. Individuals with a low four-gene signature score had a favorable outcome and better immune responses. Finally, a nomogram for clinical practice was constructed by integrating the four-gene signature and INSS stage.

Conclusions

Our findings confirm the influence of CSC features in NB prognosis. The newly developed NB stemness-related four-gene signature prognostic signature could facilitate the prognostic prediction, and the identified hub genes may serve as promising targets for individualized treatments.

Plasma cell-free DNA methylome profiling in pre- and post-surgery oral cavity squamous cell carcinoma

Head and neck squamous cell carcinoma (HNSCC), a highly heterogeneous disease that involves multiple anatomic sites, is a leading cause of cancer-related mortality worldwide. Although the utility of noninvasive biomarkers based on circulating cell-free DNA (cfDNA) methylation profiling has been widely recognized, limited studies have been reported so far regarding the dynamics of cfDNA methylome in oral cavity squamous cell carcinoma (OCSCC). It is hypothesized in this study that comparison of methylation profiles in pre- and postsurgery plasma samples will reveal OCSCC-specific prognostic and diagnostic biomarkers. As a strategy to further prioritize tumor-specific targets, top differential methylated regions (DMRs) were called by reanalyzing methylation data from paired tumor and normal tissue collected in the the cancer genome atlas head-neck squamous cell carcinoma (TCGA) head and neck cancer cohort. Matched plasma samples from eight patients with OCSCC were collected at Moffitt Cancer Center before and after surgical resection. Plasma-derived cfDNA was analyzed by cfMBD-seq, which is a high-sensitive methylation profiling assay. Differential methylation analysis was then performed based on the matched samples profiled. In the top 200 HNSCC-specific DMRs detected based on the TCGA data set, a total of 23 regions reached significance in the plasma-based DMR test. The top five validated DMR regions (ranked by the significance in the plasma study) are located in the promoter regions of genes PENK, NXPH1, ZIK1, TBXT, and CDO1, respectively. The genome-wide cfDNA DMR analysis further highlighted candidate biomarkers located in genes SFRP4, SOX1, IRF4, and PCDH17. The prognostic relevance of candidate genes was confirmed by survival analysis using the TCGA data. This study supports the utility of cfDNA-based methylome profiling as a promising noninvasive biomarker source for OCSCC and HNSCC.

GDPD5 Related to Lipid Metabolism Is a Potential Prognostic Biomarker in Neuroblastoma

Neuroblastoma (NB) is an extracranial solid tumor in children with poor prognosis in high-risk patients and its pathogenesis and prognostic markers urgently need to be explored. This study aimed to explore potential biomarkers related to NB from the aspect of lipid metabolism. Fifty-eight lipid metabolism-related differentially expressed genes between high-risk NB and non-high-risk NB in the GSE49710 dataset were analyzed using bioinformatics, including 45 down-regulated genes and 13 up-regulated genes. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis identified steroid hormone biosynthesis as an abnormal metabolic pathway in high-risk NB. Survival analysis established a three-gene prognostic model, including ACHE, GDPD5 and PIK3R1. In the test data, the AUCs of the established prognostic models used to predict patient survival at 1, 3 and 5 years were 0.84, 0.90 and 0.91, respectively. Finally, in the SH-SY5Y cell line, it was verified that overexpression of GDPD5 can inhibit cell proliferation and migration, as well as affect the lipid metabolism of SH-SY5Y, but not the sugar metabolism. hsa-miR-592 was predicted to be a potential target miRNA of GDPD5 by bioinformatics. In conclusion, this study develops a lipid-metabolism-related gene-based prognostic model for NB and demonstrates that GDPD5 inhibits SH-SY5Y proliferation and migration and may be targeted by hsa-miR-592 and inhibit SH-SY5Y fat synthesis.