Deviation from baseline mutation burden provides powerful and robust rare-variants association test for complex diseases.

Identifying rare variants that contribute to complex diseases is challenging because of the low statistical power in current tests comparing cases with controls. Here, we propose a novel and powerful rare variants association test based on the deviation of the observed mutation burden of a gene in cases from a baseline predicted by a weighted recursive truncated negative-binomial regression (RUNNER) on genomic features available from public data. Simulation studies show that RUNNER is substantially more powerful than state-of-the-art rare variant association tests and has reasonable type 1 error rates even for stratified populations or in small samples. Applied to real case-control data, RUNNER recapitulates known genes of Hirschsprung disease and Alzheimer's disease missed by current methods and detects promising new candidate genes for both disorders. In a case-only study, RUNNER successfully detected a known causal gene of amyotrophic lateral sclerosis. The present study provides a powerful and robust method to identify susceptibility genes with rare risk variants for complex diseases.

Actionable pharmacogenetic variants in Hong Kong Chinese exome sequencing data and projected prescription impact in the Hong Kong population.

Preemptive pharmacogenetic testing has the potential to improve drug dosing by providing point-of-care patient genotype information. Nonetheless, its implementation in the Chinese population is limited by the lack of population-wide data. In this study, secondary analysis of exome sequencing data was conducted to study pharmacogenomics in 1116 Hong Kong Chinese. We aimed to identify the spectrum of actionable pharmacogenetic variants and rare, predicted deleterious variants that are potentially actionable in Hong Kong Chinese, and to estimate the proportion of dispensed drugs that may potentially benefit from genotype-guided prescription. The projected preemptive pharmacogenetic testing prescription impact was evaluated based on the patient prescription data of the public healthcare system in 2019, serving 7.5 million people. Twenty-nine actionable pharmacogenetic variants/ alleles were identified in our cohort. Nearly all (99.6%) subjects carried at least one actionable pharmacogenetic variant, whereas 93.5% of subjects harbored at least one rare deleterious pharmacogenetic variant. Based on the prescription data in 2019, 13.4% of the Hong Kong population was prescribed with drugs with pharmacogenetic clinical practice guideline recommendations. The total expenditure on actionable drugs was 33,520,000 USD, and it was estimated that 8,219,000 USD (24.5%) worth of drugs were prescribed to patients with an implicated actionable phenotype. Secondary use of exome sequencing data for pharmacogenetic analysis is feasible, and preemptive pharmacogenetic testing has the potential to support prescription decisions in the Hong Kong Chinese population.

Whole-genome analysis of noncoding genetic variations identifies multiscale regulatory element perturbations associated with Hirschsprung disease.

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.

Genetics of Hirschsprung's disease.

Hirschsprung's disease (HSCR) is a classical model of enteric neuropathy, occurring in approximately 2-2.8 in 10,000 newborns. It is the commonest form of congenital bowel obstruction and is characterized by the absence of enteric ganglia in distal colon. Recent advances in genome-wide association analysis (GWAS) and next generation sequencing (NGS) studies have led to the discovery of a number of new HSCR candidate genes, thereby providing new insights into the genetic architecture and molecular mechanisms of the disease. Altogether, these findings indicated that genetic heterogeneity, variable penetrance and expressivity, and genetic interaction are the pervasive characteristics of HSCR genetics. In this review, we will provide an update on the genetic landscape of HSCR and discuss how the common and rare variants may act together to modulate the phenotypic manifestation. Translating the genetic findings to genetic risk prediction and to optimize clinical outcomes are undoubtedly the ultimate goals for genetic studies on HSCR. From this perspective, we will further discuss the major obstacles in the clinical translation of these latest genetic findings. Lastly, new measures to address these clinical challenges are suggested to advance precision medicine and to develop novel alternative therapies.

Developing Biliary Atresia-like Model by Treating Human Liver Organoids with Polyinosinic:Polycytidylic Acid (Poly (I:C)).

Background: We explored the feasibility of creating BA-like organoids by treating human liver organoids with Polyinosinic:Polycytidylic acid (Poly I:C). Methods: Organoids were developed from the liver parenchyma collected during Kasai portoenterostomy (BA) and surgery for other liver disorders (non-BA). The non-BA organoids were co-cultured with poly I:C (40 µg/mL). The organoid morphology from both samples was compared on day 17. RNA-sequencing was performed to examine the transcriptomic differences. Results: Non-BA liver organoids developed into well-expanded spherical organoids with a single-cell layer of epithelial cells and a single vacuole inside. After poly I:C treatment, the majority of these organoids developed into an aberrant morphology with a high index of similarity to BA organoids which are multi-vacuoled and/or unexpanded. RNA-sequencing analysis revealed that 19 inflammatory genes were commonly expressed in both groups. Conditional cluster analysis revealed several genes (SOCS6, SOCS6.1, ARAF, CAMK2G, GNA1C, ITGA2, PRKACA, PTEN) that are involved in immune-mediated signaling pathway had a distinct pattern of expression in the poly I:C treated organoids. This resembled the expression pattern in BA organoids (p < 0.05). Conclusions: Poly I:C treated human liver organoids exhibit morphology and genetic signature highly compatible to organoids developed from BA liver samples. They are potential research materials to study immune-mediated inflammation in BA.

Immunomodulatory Mechanisms of Mesenchymal Stem Cells and Their Potential Clinical Applications.

Mesenchymal stem cells (MSCs) are multipotent stem cells with the capacity of self-renewal, homing, and low immunogenicity. These distinct biological characteristics have already shown immense potential in regenerative medicine. MSCs also possess immunomodulatory properties that can maintain immune homeostasis when the immune response is over-activated or under-activated. The secretome of MSCs consists of cytokines, chemokines, signaling molecules, and growth factors, which effectively contribute to the regulation of immune and inflammatory responses. The immunomodulatory effects of MSCs can also be achieved through direct cell contact with microenvironmental factors and immune cells. Furthermore, preconditioned and engineered MSCs can specifically improve the immunomodulation effects in diverse clinical applications. These multifunctional properties of MSCs enable them to be used as a prospective therapeutic strategy to treat immune disorders, including autoimmune diseases and incurable inflammatory diseases. Here we review the recent exploration of immunomodulatory mechanisms of MSCs and briefly discuss the promotion of the genetically engineered MSCs. Additionally, we review the potential clinical applications of MSC-mediated immunomodulation in four types of immune diseases, including systemic lupus erythematosus, Crohn's disease, graft-versus-host disease, and COVID-19.

Current Understanding in the Clinical Characteristics and Molecular Mechanisms in Different Subtypes of Biliary Atresia.

Biliary atresia is a severe obliterative cholangiopathy in early infancy that is by far the most common cause of surgical jaundice and the most common indicator for liver transplantation in children. With the advanced knowledge gained from different clinical trials and the development of research models, a more precise clinical classification of BA (i.e., isolated BA (IBA), cystic BA (CBA), syndromic BA (SBA), and cytomegalovirus-associated BA (CMVBA)) is proposed. Different BA subtypes have similar yet distinguishable clinical manifestations. The clinical and etiological heterogeneity leads to dramatically different prognoses; hence, treatment needs to be specific. In this study, we reviewed the clinical characteristics of different BA subtypes and revealed the molecular mechanisms of their developmental contributors. We aimed to highlight the differences among these various subtypes of BA which ultimately contribute to the development of a specific management protocol for each subtype.

Actionable secondary findings in 1116 Hong Kong Chinese based on exome sequencing data.

The use of exome and genome sequencing has increased rapidly nowadays. After primary analysis, further analysis can be performed to identify secondary findings that offer medical benefit for patient care. Multiple studies have been performed to evaluate secondary findings in different ethnicities. However, relevant data are limited in Chinese. Here, with the use of a cohort consisted of 1116 Hong Kong Chinese exome sequencing data, we evaluated the secondary findings in the 59 genes recommended by the American College of Medical Genetics and Genomics. Fifteen unique pathogenic or likely pathogenic variants in 17 individuals were identified, representing a frequency of 1.52% in our cohort. Although 20 individuals harboured pathogenic or likely pathogenic variants in recessive conditions, none carried bi-allelic mutations in the same gene. Our finding was in accordance with the estimation from the American College of Medical Genetics and Genomics that about 1% individuals harbour secondary findings.

Roles of Enteric Neural Stem Cell Niche and Enteric Nervous System Development in Hirschsprung Disease.

The development of the enteric nervous system (ENS) is highly modulated by the synchronized interaction between the enteric neural crest cells (ENCCs) and the neural stem cell niche comprising the gut microenvironment. Genetic defects dysregulating the cellular behaviour(s) of the ENCCs result in incomplete innervation and hence ENS dysfunction. Hirschsprung disease (HSCR) is a rare complex neurocristopathy in which the enteric neural crest-derived cells fail to colonize the distal colon. In addition to ENS defects, increasing evidence suggests that HSCR patients may have intrinsic defects in the niche impairing the extracellular matrix (ECM)-cell interaction and/or dysregulating the cellular niche factors necessary for controlling stem cell behaviour. The niche defects in patients may compromise the regenerative capacity of the stem cell-based therapy and advocate for drug- and niche-based therapies as complementary therapeutic strategies to alleviate/enhance niche-cell interaction. Here, we provide a summary of the current understandings of the role of the enteric neural stem cell niche in modulating the development of the ENS and in the pathogenesis of HSCR. Deciphering the contribution of the niche to HSCR may provide important implications to the development of regenerative medicine for HSCR.

Beta-amyloid deposition around hepatic bile ducts is a novel pathobiological and diagnostic feature of biliary atresia.

BACKGROUND AND AIMS: Biliary atresia (BA) is a poorly understood and devastating obstructive bile duct disease of newborns. It is often diagnosed late, is incurable and frequently requires liver transplantation. In this study, we aimed to investigate the underlying pathogenesis and molecular signatures associated with BA. METHODS: We combined organoid and transcriptomic analysis to gain new insights into BA pathobiology using patient samples and a mouse model of BA. RESULTS: Liver organoids derived from patients with BA and a rhesus rotavirus A-infected mouse model of BA, exhibited aberrant morphology and disturbed apical-basal organization. Transcriptomic analysis of BA organoids revealed a shift from cholangiocyte to hepatocyte transcriptional signatures and altered beta-amyloid-related gene expression. Beta-amyloid accumulation was observed around the bile ducts in BA livers and exposure to beta-amyloid induced the aberrant morphology in control organoids. CONCLUSION: The novel observation that beta-amyloid accumulates around bile ducts in the livers of patients with BA has important pathobiological implications, as well as diagnostic potential. LAY SUMMARY: Biliary atresia is a poorly understood and devastating obstructive bile duct disease of newborns. It is often diagnosed late, is incurable and frequently requires liver transplantation. Using human and mouse 'liver mini-organs in the dish', we unexpectedly identified beta-amyloid deposition - the main pathological feature of Alzheimer's disease and cerebral amyloid angiopathy - around bile ducts in livers from patients with biliary atresia. This finding reveals a novel pathogenic mechanism that could have important diagnostic and therapeutic implications.

Risk factors for the development of post-operative enterocolitis in short segment Hirschsprung's disease.

AIM OF THE STUDY: The objective of this study is to identify risk factors associated with the development of post-operative enterocolitis (HAEC), in short segment Hirschsprung's disease (HSCR-S). METHODS: A retrospective study was carried out for post-operative patients with HSCR-S from 1997 to 2017. HSCR-S was defined as the most proximal extension of aganglionosis limited to the sigmoid colon. An episode of HAEC was defined as the presence of (1) vomiting or explosive diarrhea; (2) abdominal distension; (3) fever and (4) leukocytosis. Risk factors for the development of HACE were determined using multivariate logistic regression. MAIN RESULTS: The medical records of 96 patients were reviewed. The overall incidence of HAEC was 20.8% (n = 20) and 65.0% (n = 13) of HAEC occurred within the first year of operation. After a univariate logistic regression analysis, three risk factors for HAEC were identified: (1) presence of other major anomalies [OR: 1.43 (1.12-2.32), p = 0.041]; (2) creation of pre-operative defunctioning stoma [OR: 2.28 (1.47-3.23), p = 0.035]; (3) extension of aganglionosis to the sigmoid colon [OR: 1.89 (1.05-3.19), p = 0.049]. After multivariate logistic regression analysis, a significant association was demonstrated for creation of pre-operative defunctioning stoma [OR: 1.81 (1.08-3.22), p = 0.045] and extension of aganglionosis to the sigmoid colon [OR: 1.91 (1.37-2.98), p = 0.038]. CONCLUSIONS: The requirement of pre-operative defunctioning stoma and a more proximal extension of aganglionosis are risk factors for the development of post-operative HAEC in HSCR-S. Patients with these risk factors should be closely followed up especially during the first year after the operation.

Correction of Hirschsprung-Associated Mutations in Human Induced Pluripotent Stem Cells Via Clustered Regularly Interspaced Short Palindromic Repeats/Cas9, Restores Neural Crest Cell Function.

BACKGROUND & AIMS: Hirschsprung disease is caused by failure of enteric neural crest cells (ENCCs) to fully colonize the bowel, leading to bowel obstruction and megacolon. Heterozygous mutations in the coding region of the RET gene cause a severe form of Hirschsprung disease (total colonic aganglionosis). However, 80% of HSCR patients have short-segment Hirschsprung disease (S-HSCR), which has not been associated with genetic factors. We sought to identify mutations associated with S-HSCR, and used the clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9 gene editing system to determine how mutations affect ENCC function. METHODS: We created induced pluripotent stem cell (iPSC) lines from 1 patient with total colonic aganglionosis (with the G731del mutation in RET) and from 2 patients with S-HSCR (without a RET mutation), as well as RET+/- and RET-/- iPSCs. IMR90-iPSC cells were used as the control cell line. Migration and differentiation capacities of iPSC-derived ENCCs were analyzed in differentiation and migration assays. We searched for mutation(s) associated with S-HSCR by combining genetic and transcriptome data from patient blood- and iPSC-derived ENCCs, respectively. Mutations in the iPSCs were corrected using the CRISPR/Cas9 system. RESULTS: ENCCs derived from all iPSC lines, but not control iPSCs, had defects in migration and neuronal lineage differentiation. RET mutations were associated with differentiation and migration defects of ENCCs in vitro. Genetic and transcriptome analyses associated a mutation in the vinculin gene (VCL M209L) with S-HSCR. CRISPR/Cas9 correction of the RET G731del and VCL M209L mutations in iPSCs restored the differentiation and migration capacities of ENCCs. CONCLUSIONS: We identified mutations in VCL associated with S-HSCR. Correction of this mutation in iPSC using CRISPR/Cas9 editing, as well as the RET G731del mutation that causes Hirschsprung disease with total colonic aganglionosis, restored ENCC function. Our study demonstrates how human iPSCs can be used to identify disease-associated mutations and determine how they affect cell functions and contribute to pathogenesis.

Long-term Results and Quality of Life Assessment in Biliary Atresia Patients: A 35-Year Experience in a Tertiary Hospital.

OBJECTIVES: To review long-term transplant-free survival and quality of life (QOL) of patients with biliary atresia (BA). METHODS: A retrospective study reviewing all patients with Kasai operation between January 1, 1980 and December 31, 2015 was performed to evaluate the transplant-free survival. Subgroup analysis of patients older than 20 years was carried out to assess the QOL using the Short Form-36 Health Survey and incidences of disease-related complications. Comparison between patients with native and transplanted liver was performed using two-tailed independent samples t-test (P value < 0.05, significant). RESULTS: The 20-year Kaplan-Meier transplant-free survival of the 141 patients in our study was 51%. The subgroup analysis of long-term survivors revealed a trend of increased prevalence of complications like esophageal varices, portal hypertension, and recurrent admissions in the patient groups with raised serum bilirubin (SB).Thirty-one patients were successfully contacted for QOL assessment, 26 (16 with native liver and 10 with transplanted liver) responded (76.5%). BA patients who were documented to have active complications have a significantly lower vitality score (50.7 vs 57.5, P = 0.015). There was no statistically significant difference in the scores between the transplanted group and the disease-free control group. However, the native liver group achieved a lower score in both the general health section (42.9 vs 49.6, P = 0.029) and the overall physical component (49.6 vs 54.4, P = 0.037). CONCLUSIONS: A significant proportion of our patients survive with their native liver for more than 20 years. These long-term survivors may suffer from complications that impair their QOL. They require continuous life-long care.

Are all patients with short segment Hirschsprung's disease equal? A retrospective multicenter study.

BACKGROUND/PURPOSE: Short segment Hirschsprung's disease (HSCR) carries a better prognosis than long segment disease, but the definition of short is controversial. The objective of this study is to determine anatomically the extent of disease involvement that would be associated with a better functional outcome. METHODS: This is a retrospective multicenter (n = 3) study with patients (≥ 3 years) who had transanal pullthrough operation done for aganglionosis limited to the recto-sigmoid colon were reviewed. The extent of disease involvement and bowel resection was retrieved by reviewing the operative records as well as histopathological reports of the resected specimens. Clinical assessment was performed according to the criteria of a seven-itemed bowel function score (BFS) (maximum score = 20). Manometric assessment was performed with anorectal manometry. RESULTS: The study period started from 2003 to 45 patients were studied with median age at assessment = 52.0 months and operation = 3.0 months. The disease involvement was categorized into upper sigmoid-descending colon (DC) (n = 8), sigmoid colon (SC) (n = 12), upper rectum (UR) (n = 14) and lower rectum (LR) (n = 11) according to the level of normal biopsy result. There was no significant difference in the age of assessment between the four groups. The median BFSs in the DC, SC, UR and LR were 13, 15, 17 and 17, respectively (p = 0.01). Nine patients from the DC and SC groups reported soiling for more than twice per week. Sub-group analysis comparing patients with and without the entire sigmoid colon resected revealed worse functional outcomes in terms of the incidence of soiling (40.7 vs 22.2%, p = 0.05) and the BFS (14 vs 18, p = 0.04) in the former group. Anorectal manometry did not reveal any significant difference between the four groups, but a higher proportion of patients in the UR and LR groups appeared to have a normal sphincter resting pressure (DC vs SC vs UR vs LR = 62.5 vs 75.0 vs 85.7 vs 80.0%, p = 0.10). CONCLUSION: Patients with short segment HSCR are not equal at all. HSCR patients with aganglionosis limited to the rectum without the need of removing the entire sigmoid colon have a better bowel control and overall functional score. Less bowel loss and colonic dissection maybe the underlying reasons. Although future studies with a larger sample size and a longer follow-up period are required to validate the results of this study, it has provided a new insight to the current understanding of short segment disease in HSCR.

Meta-analysis of GWAS on two Chinese populations followed by replication identifies novel genetic variants on the X chromosome associated with systemic lupus erythematosus.

Systemic lupus erythematosus (SLE) is a prototypic autoimmune disease that affects mainly females. What role the X chromosome plays in the disease has always been an intriguing question. In this study, we examined the genetic variants on the X chromosome through meta-analysis of two genome-wide association studies (GWAS) on SLE on Chinese Han populations. Prominent association signals from the meta-analysis were replicated in 4 additional Asian cohorts, with a total of 5373 cases and 9166 matched controls. We identified a novel variant in PRPS2 on Xp22.3 as associated with SLE with genome-wide significance (rs7062536, OR = 0.84, P = 1.00E-08). Association of the L1CAM-MECP2 region with SLE was reported previously. In this study, we identified independent contributors in this region in NAA10 (rs2071128, OR = 0.81, P = 2.19E-13) and TMEM187 (rs17422, OR = 0.75, P = 1.47E-15), in addition to replicating the association from IRAK1-MECP2 region (rs1059702, OR = 0.71, P = 2.40E-18) in Asian cohorts. The X-linked susceptibility variants showed higher effect size in males than that in females, similar to results from a genome-wide survey of associated SNPs on the autosomes. These results suggest that susceptibility genes identified on the X chromosome, while contributing to disease predisposition, might not contribute significantly to the female predominance of this prototype autoimmune disease.

Somatic PIK3CA mutations in seven patients with PIK3CA-related overgrowth spectrum.

Somatic mutations in PIK3CA cause many overgrowth syndromes that have been recently coined the "PIK3CA-Related Overgrowth Spectrum." Here, we present seven molecularly confirmed patients with PIK3CA-Related Overgrowth Spectrum, including patients with Congenital Lipomatous Overgrowth, Vascular Malformations, Epidermal Nevi, Scoliosis/Skeletal and Spinal syndrome, Klippel-Trenaunay syndrome, lymphatic malformation and two with atypical phenotypes that cannot be classified into existing disease categories. The literature on PIK3CA-Related Overgrowth Spectrum, suggests that PIK3CA c.1258T>C; p.(Cys420Arg), c.1624G>A; p.(Glu542Lys), c.1633G>A; p.(Glu545Lys), c.3140A>G; p.(His1047Arg), and c.3140A>T; p.(His1047Leu) can be identified in approximately 90% of patients without brain overgrowth. Therefore, droplet digital polymerase chain reaction targeting these mutation hotspots could be used as the first-tier genetic test on patients with PIK3CA-Related Overgrowth Spectrum who do not have signs of overgrowth in their central nervous system. © 2017 Wiley Periodicals, Inc.

Three SNPs in chromosome 11q23.3 are independently associated with systemic lupus erythematosus in Asians.

Systemic lupus erythematosus (SLE) has a complex etiology and is affected by both genetic and environmental factors. Although more than 40 loci have shown robust association with SLE, the details of these loci, such as the independent contributors and the genes involved, are still unclear. In this study, we performed meta-analysis of two existing genome-wide association studies (GWASs) on Chinese Han populations from Hong Kong and Anhui, China, and followed the findings by further replication on three additional Chinese and Thailand cohorts with a total of 4254 cases and 6262 controls matched geographically and ethnically. We discovered multiple susceptibility variants for SLE in the 11q23.3 region, including variants in/near PHLDB1 (rs11603023, P(_combined) = 1.25E-08, OR = 1.20), DDX6 (rs638893, P(_combined) = 5.19E-07, OR = 1.22) and CXCR5 (rs10892301, P(_combined) = 2.51E-08, OR = 0.85). Genetic contributions from the newly identified variants were all independent of SNP rs4639966, whose association was reported from the previous GWAS. In addition, the three newly identified variants all showed independent association with the disease through modeling by both stepwise and conditional logistic regression. The presence of multiple independent variants in this region emphasizes its role in SLE susceptibility, and also hints the possibility that distinct biological mechanisms might be involved in the disease involving this genomic region.

Identification of GLI Mutations in Patients With Hirschsprung Disease That Disrupt Enteric Nervous System Development in Mice.

BACKGROUND & AIMS: Hirschsprung disease is characterized by a deficit in enteric neurons, which are derived from neural crest cells (NCCs). Aberrant hedgehog signaling disrupts NCC differentiation and might cause Hirschsprung disease. We performed genetic analyses to determine whether hedgehog signaling is involved in pathogenesis. METHODS: We performed deep-target sequencing of DNA from 20 patients with Hirschsprung disease (16 men, 4 women), and 20 individuals without (controls), and searched for mutation(s) in GLI1, GLI2, GLI3, SUFU, and SOX10. Biological effects of GLI mutations were tested in luciferase reporter assays using HeLa or neuroblastoma cell lines. Development of the enteric nervous system was studied in Sufu(f/f), Gli3(Δ699), Wnt1-Cre, and Sox10(NGFP) mice using immunohistochemical and whole-mount staining procedures to quantify enteric neurons and glia and analyze axon fasciculation, respectively. NCC migration was studied using time-lapse imaging. RESULTS: We identified 3 mutations in GLI in 5 patients with Hirschsprung disease but no controls; all lead to increased transcription of SOX10 in cell lines. SUFU, GLI, and SOX10 form a regulatory loop that controls the neuronal vs glial lineages and migration of NCCs. Sufu mutants mice had high Gli activity, due to loss of Sufu, disrupting the regulatory loop and migration of enteric NCCs, leading to defective axonal fasciculation, delayed gut colonization, or intestinal hypoganglionosis. The ratio of enteric neurons to glia correlated inversely with Gli activity. CONCLUSIONS: We identified mutations that increase GLI activity in patients with Hirschsprung disease. Disruption of the SUFU-GLI-SOX10 regulatory loop disrupts migration of NCCs and development of the enteric nervous system in mice.