Identifying prognostic hub genes and key pathways in pediatric adrenocortical tumors through RNA sequencing and Co-expression analysis

Pediatric adrenocortical tumors (ACTs), rare conditions with uncertain prognoses, have high incidence in southern and southeastern Brazil. Pediatric ACTs are highly heterogeneous, so establishing prognostic markers for these tumors is challenging. We have conducted transcriptomic analysis on 14 pediatric ACT samples and compared cases with favorable and unfavorable clinical outcomes to identify prognostically significant genes. This comparison showed 1257 differentially expressed genes in favorable and unfavorable cases. Among these genes, 15 out of 60 hub genes were significantly associated with five-year event-free survival (EFS), and 10 had significant diagnostic value for predicting ACT outcomes in an independent microarray dataset of pediatric adrenocortical carcinomas (GSE76019). Overexpression of N4BP2, HSPB6, JUN, APBB1IP, STK17B, CSNK1D, and KDM3A was associated with poorer EFS, whereas lower expression of ISCU, PTPR, PRKAB2, CD48, PRF1, ITGAL, KLK15, and HIST1H3J was associated with worse outcomes. Collectively, these findings underscore the prognostic significance of these hub genes and suggest that they play a potential role in pediatric ACT progression and are useful predictors of clinical outcomes.