Combinatorial multiomic analysis from a pedigree of Sox10Dom Hirschsprung mice identifies multiple high confidence...

Hirschsprung disease (HSCR) is characterized by absence of enteric ganglia (aganglionosis) along variable lengths of the distal intestine. This disorder results from deficient colonization of fetal intestine by enteric neural crest-derived cells (ENCDCs). HSCR exhibits complex, multifactorial inheritance with penetrance and severity varying widely even within families. SOX10 is among causal genes that predispose to aganglionosis. Yet, how gene interactions influence severity of HSCR aganglionosis is not understood. Prior mapping of aganglionosis modifiers was achieved in a standard F1-intercross utilizing the Sox10Dom HSCR mouse model. Here we deploy a novel strategy of genotyping an extended pedigree pedigree of Sox10Dom mice on a mixed genetic background. GWAS in this pedigree points to novel aganglionosis modifier intervals with replication and refinement of prior modifier regions. Complementary omics analysis of the developing Enteric Nervous System (ENS) enabled identification of multiple high-priority candidate genes within these modifier intervals based on gene expression, chromatin accessibility, and presence of conserved SOX10 binding motifs. We implemented a prioritization pipeline for ranking potential modifiers that generated candidate lists including several well-known for effects on ENS development as well as multiple novel genes. Among the novel genes, Dach1 ranked as a top priority candidate gene for modifying migration of ENCDCs and thus influencing aganglionosis severity. The results identify genome intervals with intrinsic genes that are logical candidates for modifying Sox10Dom aganglionosis severity. We also note that several human orthologs to aganglionosis modifier candidate genes are within linkage disequilibrium blocks containing genetic variants associated with human gut motility disorders, which offers opportunity for gaining biological insight into human HSCR severity.