Retrotrapezoid nucleus (RTN) neurons in the brainstem regulate the ventilatory response to hypercarbia. It is unclear how PHOX2B-polyalanine repeat mutations (PHOX2B-PARMs) alter the function of PHOX2B and perturb the formation of RTN neurons. Here, we generated human brainstem organoids (HBSOs) with RTN-like neurons from human pluripotent stem cells. Single-cell transcriptomics revealed that expression of PHOX2B+7Ala PARM alters the differentiation trajectories of the hindbrain neurons and hampers the formation of the RTN-like neurons in HBSOs. With the unguided cerebral organoids (HCOs), PHOX2B+7Ala PARM interrupted the patterning of PHOX2B+ neurons with dysregulation of Hedgehog pathway and HOX genes. With complementary use of HBSOs and HCOs with a patient and two mutant induced pluripotent stem cell lines carrying different polyalanine repetition in PHOX2B, we further defined the association between the length of polyalanine repetition and malformation of RTN-respiratory center and demonstrated the potential toxic gain of function of PHOX2B-PARMs, highlighting the uniqueness of these organoid models for disease modeling.

Hirschsprung disease, or congenital aganglionosis, is believed to be oligogenic-that is, caused by multiple genetic factors. We performed whole-genome sequence analyses of patients with Hirschsprung disease to identify genetic factors that contribute to disease development and analyzed the functional effects of these variants.

Hirschsprung disease is characterized by the absence of enteric neurons caused by the defects of enteric neural crest cells, leading to intestinal obstruction. Here, using induced pluripotent stem cell-based models of Hirschsprung and single-cell transcriptomic analysis, we identify a gene set of 118 genes commonly dysregulated in all patient enteric neural crest cells, and suggest HDAC1 may be a key regulator of these genes. Furthermore, upregulation of RNA splicing mediators and enhanced alternative splicing events are associated with severe form of Hirschsprung. In particular, the higher inclusion rate of exon 9 in PTBP1 and the perturbed expression of a PTBP1-target, PKM, are significantly enriched in these patient cells, and associated with the defective oxidative phosphorylation and impaired neurogenesis. Hedgehog-induced oxidative phosphorylation significantly enhances the survival and differentiation capacity of patient cells. In sum, we define various factors associated with Hirschsprung pathogenesis and demonstrate the implications of oxidative phosphorylation in enteric neural crest development and HSCR pathogenesis.

Gastrointestinal motility disorders occur frequently in patients with ciliopathy, but the underlying genetic link is unclear. The ciliary protein Kif7 can positively or negatively regulate Hedgehog signaling in different cellular contexts. Mice with neural crest cell (NCC)–specific Kif7 deficiency show a marked reduction of enteric NOS+ inhibitory neurons. Malformation of enteric nervous system (ENS) causes growth retardation and gut motility defect in mice. Mechanistically, Kif7 inhibits Gli2 in enteric NCCs (ENCCs), where Gli2 positively regulates the expression of Ezh2 by inhibiting the miR124-mediated suppression. In developing ENCCs, Ezh2 is a master regulator of 102 core genes underlying ENCC differentiation. Deletion of Gli2 or inhibition of Ezh2 favors the neurogenic lineage differentiation of mouse and human ENCCs and rescues the ENS defects of Kif7 mutants. In summary, Hedgehog signal, via Kif7-Gli-Ezh2, controls the timely expressions of the core genes to mediate the differentiation of ENCCs.

It is widely recognized that noncoding genetic variants play important roles in many human diseases, but there are multiple challenges that hinder the identification of functional disease-associated noncoding variants. The number of noncoding variants can be many times that of coding variants; many of them are not functional but in linkage disequilibrium with the functional ones; different variants can have epistatic effects; different variants can affect the same genes or pathways in different individuals; and some variants are related to each other not by affecting the same gene but by affecting the binding of the same upstream regulator. To overcome these difficulties, we propose a novel analysis framework that considers convergent impacts of different genetic variants on protein binding, which provides multiscale information about disease-associated perturbations of regulatory elements, genes, and pathways. Applying it to our whole-genome sequencing data of 918 short-segment Hirschsprung disease patients and matched controls, we identify various novel genes not detected by standard single-variant and region-based tests, functionally centering on neural crest migration and development. Our framework also identifies upstream regulators whose binding is influenced by the noncoding variants. Using human neural crest cells, we confirm cell stage-specific regulatory roles of three top novel regulatory elements on our list, respectively in the RET, RASGEF1A, and PIK3C2B loci. In the PIK3C2B regulatory element, we further show that a noncoding variant found only in the patients affects the binding of the gliogenesis regulator NFIA, with a corresponding up-regulation of multiple genes in the same topologically associating domain.

It has been a challenge to develop fully functioning cells from human pluripotent stem cells (hPSCs). We investigated how activation of hedgehog signaling regulates derivation of enteric neural crest (NC) cells from hPSCs.