Unraveling the molecular basis of clinical heterogeneity in Hirschsprung disease through gene expression analysis