Abstract:Objective: This study aims to identify differentially expressed genes (DEGs) in neuroblastoma (NB) through comprehensive bioinformatics analysis and machine learning techniques. We seek to elucidate these DEGs' biological functions and associated signaling pathways. Furthermore, our objective extends to predicting upstream microRNAs (miRNAs) and relevant transcription factors of pivotal genes, with the ultimate goal of guiding clinical diagnostics and informing future treatment strategies for Neuroblastoma. Methods: In this study, we sourced datasets GSE49710 and TARGET from the GEO and UCSC-XENA databases, respectively. Differentially expressed genes (DEGs) were identified using the R language "limma" package. Subsequent Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses of these DEGs were conducted using the "clusterProfiler" package. We employed Weighted Gene Co-expression Network Analysis (WGCNA) to isolate the most significant modules associated with death and MYCN amplification, specifically MEpink and MEbrown modules. These modules were then cross-referenced with the DEGs for further GO and KEGG pathway analyses. LASSO regression analysis, facilitated by the "glmnet" package, was utilized to pinpoint three hub genes. We performed differential analysis on these genes and constructed Receiver Operating Characteristic (ROC) curves for disease diagnosis purposes. Immune infiltration analysis was conducted using the "GSVA" package's ssGSEA function. Additionally, single-gene Gene Set Enrichment Analysis (GSEA) on the hub gene was carried out based on Reactome and KEGG databases. Upstream miRNA and transcription factors associated with the hub gene were predicted using RegNetwork, with visual representations created in Cytoscape. Furthermore, to validate the three identified markers in neuroblastoma tissues, quantitative Real-Time Polymerase Chain Reaction (qRT-PCR) analysis was conducted. Results: We identified 483 differentially expressed genes (DEGs) in neuroblastoma. These genes predominantly function in protein translation, membrane composition, and RNA transcription regulation, and are implicated in multiple signaling pathways relevant to neurodegenerative diseases. Utilizing LASSO regression analysis, we pinpointed three hub genes: VGF, DGKD, and C19orf52. The Receiver Operating Characteristic (ROC) curve analysis yielded Area Under Curve (AUC) values of 0.751 and 0.722 for VGF, 0.79 and 0.656 for DGKD, and 0.8 and 0.753 for C19orf52, respectively. Our immune infiltration analysis revealed significant correlations among monocytes, follicular helper T cells, and CD4+ T cells. Notably, in the death group, we observed heightened infiltration levels of activated CD4+ T cells, macrophages, and Th2 cells. C19orf52 exhibited a close association with the infiltration of monocytes, CD4+ T cells, and Th2 cells, with P-values less than 0.05. Furthermore, qRT-PCR analysis corroborated the upregulation of VGF in neuroblastoma tissues, further validating our findings. Conclusion: The hub genes (VGF, DGKD, and C19orf52) of neuroblastoma are screened. VGF, one of the hub genes, may have a high diagnostic value and is involved in the immune cell infiltration in neuroblastoma tissue, which may be used as a biomarker for the diagnosis of neuroblastoma and provides a new direction for clinical prognosis prediction and management improvement.

Clinical Significance of a CD3/CD8-Based Immunoscore in Neuroblastoma Patients Using Digital Pathology

Prognostic Value of Immune Score in Nasopharyngeal Carcinoma Using Digital Pathology

Predictive Value Of Immunoscore In Patients With Colorectal Cancer

Cross-cohort analysis identified an immune checkpoint-based signature to predict the clinical outcomes of neuroblastoma

CD20highCD138low Tumor-Infiltrating Lymphocytes Predominantly Related to Cytokine‒cytokine Receptor Interactions Are Associated with Favorable Outcomes in Neuroblastoma Patients

Cellular components in tumor microenvironment of neuroblastoma and the prognostic value.

Prognostic Impact of the Immunoscore Based on Whole-Slide Image Analysis of CD3+ Tumor-Infiltrating Lymphocytes in Diffuse Large B-Cell Lymphoma

Establishment and Validation of a Prognostic Immune Signature in Neuroblastoma

Prognostic Signature of Immune Genes and Immune-Related LncRNAs in Neuroblastoma: A Study Based on GEO and TARGET Datasets

Tumor Microenvironment Profiling Identifies Prognostic Signatures and Suggests Immunotherapeutic Benefits in Neuroblastoma

Prognostic and Predictive Value of an Immunoscore Signature in Glioblastoma Multiform

Exploration Of The Molecular Characteristics Of The Tumor-Immune Interaction And The Development Of An Individualized Immune Prognostic Signature For Neuroblastoma

The relationship between circulating tumor cells in peripheral blood and clinical characteristics of pediatric neuroblastoma and prognostic evaluation

Association of image-defined risk factors, tumor resectability, and prognosis in children with localized neuroblastoma

Clinical Impact of the Tumor Immune Microenvironment in Completely Resected Stage IIIA(N2) Non-Small Cell Lung Cancer Based on an Immunoscore Approach

Immunoscore Combining CD8, FoxP3, and CD68-Positive Cells Density and Distribution Predicts the Prognosis of Head and Neck Cancer Patients

Identification of an Autophagy-Related Risk Signature Correlates With Immunophenotype and Predicts Immune Checkpoint Blockade Efficacy of Neuroblastoma

Bioinformatic Identification of Neuroblastoma Microenvironment-Associated Biomarkers with Prognostic Value

Bioinformatics-Based Identification of Key Prognostic Genes in Neuroblastoma with a Focus on Immune Cell Infiltration and Diagnostic Potential of VGF

A signature based on five immune-related genes to predict the survival and immune characteristics of neuroblastoma

Clinical Immunoprofiling Reveals that High Numbers of CD8+ and PD-1+ Cells Predict Superior Patient Survival Across Major Cancer Types Independent of Major Risk Factors