Background Neuroblastoma (NBL) is the most common extracranial solid tumor in childhood, and patients with high-risk neuroblastoma had a relatively poor prognosis despite multimodal treatment. To improve immunotherapy efficacy in neuroblastoma, systematic profiling of the immune landscape in neuroblastoma is an urgent need. Methods RNA-seq and according clinical information of neuroblastoma were downloaded from the TARGET database and GEO database (GSE62564). With an immune-related-gene set obtained from the ImmPort database, Immune-related Prognostic Gene Pairs for Neuroblastoma (IPGPN) for overall survival (OS) were established with the TARGET-NBL cohort and then verified with the GEO-NBL cohort. Immune cell infiltration analysis was subsequently performed. The integrated model was established with IPGPN and clinicopathological parameters. Immune cell infiltration was analyzed with the XCELL algorithm. Functional enrichment analysis was performed with clusterProfiler package in R. Results Immune-related Prognostic Gene Pairs for Neuroblastoma was successfully established with seven immune-related gene pairs (IGPs) involving 13 unique genes in the training cohort. In the training cohort, IPGPN successfully stratified neuroblastoma patients into a high and low immune-risk groups with different OS (HR=3.92, P = 2 × 10−8) and event-free survival (HR=3.66, P=2 × 10−8). ROC curve analysis confirmed its predictive power. Consistently, high IPGPN also predicted worse OS (HR=1.84, P = .002) and EFS in validation cohort (HR=1.38, P = .06) Moreover, higher activated dendritic cells, M1 macrophage, Th1 CD4+, and Th2 CD4+ T cell enrichment were evident in low immune-risk group. Further integrating IPGPN with age and stage demonstrated improved predictive performance than IPGPN alone. Conclusion Herein, we presented an immune landscape with IPGPN for prognosis prediction in neuroblastoma, which complements the present understanding of the immune signature in neuroblastoma.