The glycosphingolipid GD2 as an effective but enigmatic target of passive immunotherapy in children with aggressive neuroblastoma (HR-NBL)

Molecular profiling of core immune-escape genes highlights TNFAIP3 as an immune-related prognostic biomarker in neuroblastoma

BACKGROUND

Neuroblastoma (NB) is the most prevalent and deadliest pediatric solid tumor. With of over 50% of high-risk neuroblastoma cases relapse, the imperative for novel drug targets and therapeutic strategies is accentuated. In neuroblastoma, the existence of tumor-associated macrophages (TAMs) correlates with an unfavorable patient prognosis. However, the clinical relevance and prognostic implications of regulatory genes linked to TAMs infiltration in neuroblastoma remain unclear, and further study is required.

METHODS

We conducted a comprehensive analysis utilizing transcriptome expression profiles from three primary datasets associated with neuroblastoma (GSE45547, GSE49710, TARGET) to identify hub genes implicated in immune evasion within neuroblastoma. Subsequently, we utilized single-cell RNA sequencing analysis on 17 clinical neuroblastoma samples to investigate the expression and distribution of these hub genes, leading to the identification of TNFAIP3. The above three public databases were merged to allowed for the validation of TNFAIP3's molecular functions through GO and KEGG analysis. Furthermore, we assessed TNFAIP3's correlation with immune infiltration and its potential immunotherapeutic impact by multiple algorithms. Our single-cell transcriptome data revealed the role of TNFAIP3 in macrophage polarization. Finally, preliminary experimental verifications to confirm the biological functions of TNFAIP3-mediated TAMs in NB.

RESULTS

A total of 6 genes related to immune evasion were screened and we found that TNFAIP3 exhibited notably higher expression in macrophages than other immune cell types, based on the scRNA-sequencing data. GO and KEGG analysis showed that low expression of TNFAIP3 significantly correlated with the activation of multiple oncogenic pathways as well as immune-related pathways. Then validation affirmed that individuals within the TNFAIP3 high-expression cohort could potentially derive greater advantages from immunotherapeutic interventions, alongside exhibiting heightened immune responsiveness. Deciphering the pseudotime trajectory of macrophages, we revealed the potential of TNFAIP3 in inducing the polarization of macrophages towards the M1 phenotype. Finally, we confirmed that patients in the TNFAIP3 high expression group might benefit more from immunotherapy or chemotherapy as substantiated by RT-qPCR and immunofluorescence examinations. Moreover, the role of TNFAIP3 in macrophage polarization was validated. Preliminary experiment showed that TNFAIP3-mediated TAMs inhibit the proliferation, migration and invasion capabilities of NB cells.

CONCLUSIONS

Our results suggest that TNFAIP3 was first identified as a promising biomarker for immunotherapy and potential molecular target in NB. Besides, the presence of TNFAIP3 within TAMs may offer a novel therapeutic strategy for NB.

Identification of a Glycosyltransferase Signature for Predicting Prognosis and Immune Microenvironment in Neuroblastoma

Neuroblastoma (NB) is one of the most common solid tumors in children. Glycosyltransferases (GTs) play a crucial role in tumor development and immune escape and have been used as prognostic biomarkers in various tumors. However, the biological functions and prognostic significance of GTs in NB remain poorly understood. The expression data from Gene Expression Omnibus (GEO) and Therapeutically Applicable Research to Generate Effective Treatments (TARGET) were collected as training and testing data. Based on a progression status, differentially expressed GTs were identified. We constructed a GTscore through support vector machine, least absolute shrinkage and selection operator, and Cox regression in NB, which included four prognostic GTs and was an independent prognostic risk factor for NB. Patients in the high GTscore group had an older age, MYCN amplification, advanced International Neuroblastoma Staging System stage, and high risk. Samples with high GTscores revealed high disialoganglioside (GD2) and neuron-specific enolase expression levels. In addition, a lack of immune cell infiltration was observed in the high GTscore group. This GTscore was also associated with the expression of chemokines (CCL2, CXCL9, and CXCL10) and immune checkpoint genes (cytotoxic T-lymphocyte–associated protein 4, granzyme H, and granzyme K). A low GTscore was also linked to an enhanced response to anti–PD-1 immunotherapy in melanoma patients, and one type of tumor was also derived from neuroectodermal cells such as NB. In conclusion, the constructed GTscore revealed the relationship between GT expression and the NB outcome, GD2 phenotype, and immune infiltration and provided novel clues for the prediction of prognosis and immunotherapy response in NB.