Integrative Single-Cell and Multi-Cohort Analysis of the Netrin-1 Signaling Pathway Reveals Divergent Prognostic Trends and Sex-Dimorphic Associations in Glioblastoma

Abstract

BackgroundThe Netrin-1 dependence receptor pathway plays critical roles in neural development, but its expression landscape and prognostic significance in glioblastoma (GBM) remain poorly characterized. MethodsSingle-cell RNA-seq data from 148,019 cells across 34 tumors (Neftel et al., 2019) were analyzed to map Netrin-1 pathway gene expression across GBM cellular states. Differential gene expression and pathway enrichment analyses were performed on NEO1-defined subpopulations. Bulk RNA-seq survival analysis was conducted across three independent GBM cohorts TCGA (n=106), CGGA mRNAseq_325 (n=137), and CGGA mRNAseq_693 (n=237), totaling 480 patients. Primary analysis used continuous Cox regression (per-SD hazard ratios); meta-analysis employed fixed-effects inverse-variance weighting. ResultsIn GBM single-cell data, Netrin-1 pathway genes showed state-specific enrichment --NEO1, DCC, NTN1, and RGMB were predominantly expressed in oligodendrocyte-precursor (OPC) and neural-progenitor (NPC) states. Cells positive for NEO1 were enriched for neural differentiation programs (nervous system development, p=9.6x10-; Axon Guidance, p=2.8x10-), whereas NEO1-negative cells were dominated by ribosomal/translational and immune activation programs. In the 3-cohort survival meta-analysis, NTN1 (Netrin-1 ligand) emerged as the sole gene reaching meta-analytic significance as a risk factor (Meta HR=1.163 per SD, 95% CI 1.056-1.281, p=0.0021, I{superscript 2}=0%, 3/3 cohorts concordant), while DCC and RGMB showed directionally consistent protective trends (DCC: Meta HR=0.938, 95% CI 0.858-1.025, p=0.156; RGMB: Meta HR=0.979, 95% CI 0.881-1.087, p=0.686; both 3/3 cohorts concordant). NEO1 itself did not independently predict survival (Meta HR=1.008, 95% CI 0.885-1.147, p=0.910). After Bonferroni correction for 10 genes tested (threshold p<0.005), only NTN1 met strict significance. In exploratory sex-stratified analysis of a single cohort (CGGA 693, n=237), NEO1 and NTN1 exhibited female-specific risk enhancement (NEO1: HR=1.417, p=0.014; NTN1: HR=1.249, p=0.019), with minimal effects in males. UNC5B showed context-dependent risk in MGMT-unmethylated tumors (HR=1.331, p=0.037). These sex-dimorphic findings require independent validation. ConclusionsThe Netrin-1 pathway exhibits divergent prognostic trends in GBM, with NTN1 as a risk factor and DCC trending toward protection--consistent with the dependence receptor model. These findings, which should be interpreted as hypothesis-generating, nominate NTN1 as a candidate therapeutic target and highlight the potential importance of sex-stratified evaluation in future Netrin-1-directed trials. Independent replication in larger cohorts is warranted.