Unbiased characterization of the larval zebrafish enteric nervous system at a single cell transcriptomic level.

The enteric nervous system (ENS) regulates many gastrointestinal functions including peristalsis, immune regulation and uptake of nutrients. Defects in the ENS can lead to severe enteric neuropathies such as Hirschsprung disease (HSCR). Zebrafish have proven to be fruitful in the identification of genes involved in ENS development and HSCR pathogenesis. However, composition and specification of enteric neurons and glial subtypes at larval stages, remains mainly unexplored. Here, we performed single cell RNA sequencing of zebrafish ENS at 5 days post-fertilization. We identified vagal neural crest progenitors, Schwann cell precursors, and four clusters of differentiated neurons. In addition, a previously unrecognized elavl3+/phox2bb-population of neurons and cx43+/phox2bb-enteric glia was found. Pseudotime analysis supported binary neurogenic branching of ENS differentiation, driven by a notch-responsive state. Taken together, we provide new insights on ENS development and specification, proving that the zebrafish is a valuable model for the study of congenital enteric neuropathies.

A combinatorial panel for flow cytometry-based isolation of enteric nervous system cells from human intestine.

Efficient isolation of neurons and glia from the human enteric nervous system (ENS) is challenging because of their rare and fragile nature. Here, we describe a staining panel to enrich ENS cells from the human intestine by fluorescence-activated cell sorting (FACS). We find that CD56/CD90/CD24 co-expression labels ENS cells with higher specificity and resolution than previous methods. Surprisingly, neuronal (CD24, TUBB3) and glial (SOX10) selective markers appear co-expressed by all ENS cells. We demonstrate that this contradictory staining pattern is mainly driven by neuronal fragments, either free or attached to glial cells, which are the most abundant cell types. Live neurons can be enriched by the highest CD24 and CD90 levels. By applying our protocol to isolate ENS cells for single-cell RNA sequencing, we show that these cells can be obtained with high quality, enabling interrogation of the human ENS transcriptome. Taken together, we present a selective FACS protocol that allows enrichment and discrimination of human ENS cells, opening up new avenues to study this complex system in health and disease.

The long Filamin-A isoform is required for intestinal development and motility: implications for chronic intestinal pseudo-obstruction.

Filamin A (FLNA) is a cytoplasmic actin binding protein, recently shown to be expressed as a long and short isoform. Mutations in FLNA are associated with a wide spectrum of disorders, including an X-linked form of chronic intestinal pseudo-obstruction (CIPO). However, the role of FLNA in intestinal development and function is largely unknown. In this study, we show that FLNA is expressed in the muscle layer of the small intestine from early human fetal stages. Expression of FLNA variants associated with CIPO, blocked expression of the long flna isoform and led to an overall reduction of RNA and protein levels. As a consequence, contractility of human intestinal smooth muscle cells was affected. Lastly, our transgenic zebrafish line showed that the flna long isoform is required for intestinal elongation and peristalsis. Histological analysis revealed structural and architectural changes in the intestinal smooth muscle of homozygous fish, likely triggered by the abnormal expression of intestinal smooth muscle markers. No defect in the localization or numbers of enteric neurons was observed. Taken together, our study demonstrates that the long FLNA isoform contributes to intestinal development and function. Since loss of the long FLNA isoform does not seem to affect the enteric nervous system, it likely results in a myopathic form of CIPO, bringing new insights to disease pathogenesis.

TFAP2B Haploinsufficiency Impacts Gastrointestinal Function and Leads to Pediatric Intestinal Pseudo-obstruction.

Background: Pediatric Intestinal Pseudo-obstruction (PIPO) is a congenital enteric disorder characterized by severe gastrointestinal (GI) dysmotility, without mechanical obstruction. Although several genes have been described to cause this disease, most patients do not receive a genetic diagnosis. Here, we aim to identify the genetic cause of PIPO in a patient diagnosed with severe intestinal dysmotility shortly after birth. Methods: Whole exome sequencing (WES) was performed in the patient and unaffected parents, in a diagnostic setting. After identification of the potential disease-causing variant, its functional consequences were determined in vitro and in vivo. For this, expression constructs with and without the causing variant, were overexpressed in HEK293 cells. To investigate the role of the candidate gene in GI development and function, a zebrafish model was generated where its expression was disrupted using CRISPR/Cas9 editing. Results: WES analysis identified a de novo heterozygous deletion in TFAP2B (NM_003221.4:c.602-5_606delTCTAGTTCCA), classified as a variant of unknown significance. In vitro studies showed that this deletion affects RNA splicing and results in loss of exon 4, leading to the appearance of a premature stop codon and absence of TFAP2B protein. Disruption of tfap2b in zebrafish led to decreased enteric neuronal numbers and delayed transit time. However, no defects in neuronal differentiation were detected. tfap2b crispants also showed decreased levels of ednrbb mRNA, a downstream target of tfap2b. Conclusion: We showed that TFAP2B haploinsufficiency leads to reduced neuronal numbers and GI dysmotility, suggesting for the first time, that this gene is involved in PIPO pathogenesis.

Corrigendum: TALPID3/KIAA0586 Regulates Multiple Aspects of Neuromuscular Patterning During Gastrointestinal Development in Animal Models and Human.

[This corrects the article DOI: 10.3389/fnmol.2021.757646.].

Universal Immunohistochemistry for Lynch Syndrome: A Systematic Review and Meta-analysis of 58,580 Colorectal Carcinomas.

BACKGROUND & AIMS: Lynch syndrome is a form of hereditary colorectal cancer (CRC) caused by pathogenic germline variants (PV) in DNA mismatch repair (MMR) genes. Currently, many Western countries perform universal immunohistochemistry testing on CRC to increase the identification of Lynch syndrome patients and their relatives. For a clear understanding of health benefits and costs, data on its outcomes are required: proportions of Lynch syndrome, sporadic MMR-deficient (MMRd) cases, and unexplained MMRd cases. METHODS: Ovid Medline, Embase, and Cochrane CENTRAL were searched for studies reporting on universal MMR immunohistochemistry, followed by MMR germline analysis, until March 20, 2020. Proportions were calculated, subgroup analyses were performed based on age and diagnostics used, and random effects meta-analyses were conducted. Quality was assessed using the Joanna Briggs Critical Appraisal Tool for Prevalence Studies. RESULTS: Of 2723 identified articles, 56 studies covering 58,580 CRCs were included. In 6.22% (95% CI, 5.08%-7.61%; I2 = 96%) MMRd was identified. MMR germline PV was present in 2.00% (95% CI, 1.59%-2.50%; I2 = 92%), ranging from 1.80% to 7.27% based on completeness of diagnostics and age restriction. Immunohistochemistry outcomes were missing in 11.81%, and germline testing was performed in 76.30% of eligible patients. In 7 studies, including 6848 CRCs completing all diagnostic stages, germline PV and biallelic somatic MMR inactivation were found in 3.01% and 1.75%, respectively; 0.61% remained unexplained MMRd. CONCLUSIONS: Age, completeness, and type of diagnostics affect the percentage of MMR PV and unexplained MMRd percentages. Complete diagnostics explain almost all MMRd CRCs, reducing the amount of subsequent multigene panel testing. This contributes to optimizing testing and surveillance in MMRd CRC patients and relatives.

The Somatic Mutation Paradigm in Congenital Malformations: Hirschsprung Disease as a Model.

Patients with Hirschsprung disease (HSCR) do not always receive a genetic diagnosis after routine screening in clinical practice. One of the reasons for this could be that the causal mutation is not present in the cell types that are usually tested-whole blood, dermal fibroblasts or saliva-but is only in the affected tissue. Such mutations are called somatic, and can occur in a given cell at any stage of development after conception. They will then be present in all subsequent daughter cells. Here, we investigated the presence of somatic mutations in HSCR patients. For this, whole-exome sequencing and copy number analysis were performed in DNA isolated from purified enteric neural crest cells (ENCCs) and blood or fibroblasts of the same patient. Variants identified were subsequently validated by Sanger sequencing. Several somatic variants were identified in all patients, but causative mutations for HSCR were not specifically identified in the ENCCs of these patients. Larger copy number variants were also not found to be specific to ENCCs. Therefore, we believe that somatic mutations are unlikely to be identified, if causative for HSCR. Here, we postulate various modes of development following the occurrence of a somatic mutation, to describe the challenges in detecting such mutations, and hypothesize how somatic mutations may contribute to 'missing heritability' in developmental defects.

Size matters: Large copy number losses in Hirschsprung disease patients reveal genes involved in enteric nervous system development.

Hirschsprung disease (HSCR) is a complex genetic disease characterized by absence of ganglia in the intestine. HSCR etiology can be explained by a unique combination of genetic alterations: rare coding variants, predisposing haplotypes and Copy Number Variation (CNV). Approximately 18% of patients have additional anatomical malformations or neurological symptoms (HSCR-AAM). Pinpointing the responsible culprits within a CNV is challenging as often many genes are affected. Therefore, we selected candidate genes based on gene enrichment strategies using mouse enteric nervous system transcriptomes and constraint metrics. Next, we used a zebrafish model to investigate whether loss of these genes affects enteric neuron development in vivo. This study included three groups of patients, two groups without coding variants in disease associated genes: HSCR-AAM and HSCR patients without associated anomalies (HSCR-isolated). The third group consisted of all HSCR patients in which a confirmed pathogenic rare coding variant was identified. We compared these patient groups to unaffected controls. Predisposing haplotypes were determined, confirming that every HSCR subgroup had increased contributions of predisposing haplotypes, but their contribution was highest in isolated HSCR patients without RET coding variants. CNV profiling proved that specifically HSCR-AAM patients had larger Copy Number (CN) losses. Gene enrichment strategies using mouse enteric nervous system transcriptomes and constraint metrics were used to determine plausible candidate genes located within CN losses. Validation in zebrafish using CRISPR/Cas9 targeting confirmed the contribution of UFD1L, TBX2, SLC8A1, and MAPK8 to ENS development. In addition, we revealed epistasis between reduced Ret and Gnl1 expression and between reduced Ret and Tubb5 expression in vivo. Rare large CN losses-often de novo-contribute to HSCR in HSCR-AAM patients. We proved the involvement of six genes in enteric nervous system development and Hirschsprung disease.

Intestinal multicellular organoids to study colorectal cancer.

Modeling colorectal cancer (CRC) using organoids has burgeoned in the last decade, providing enhanced in vitro models to study the development and possible treatment options for this type of cancer. In this review, we describe both normal and CRC intestinal organoid models and their utility in the cancer research field. Besides highlighting studies that develop epithelial CRC organoid models, i.e. organoids without tumor microenvironment (TME) cellular components, we emphasize on the need for TME in CRC modeling, to help reduce translational disparities in this area. Also, we discuss the utilization of CRC organoids derived from pluripotent stem cells, as well as their potential to be used in cancer research. Finally, limitations and challenges in the current CRC organoids field, are discussed.

Loss of enteric neuronal Ndrg4 promotes colorectal cancer via increased release of Nid1 and Fbln2.

The N-Myc Downstream-Regulated Gene 4 (NDRG4), a prominent biomarker for colorectal cancer (CRC), is specifically expressed by enteric neurons. Considering that nerves are important members of the tumor microenvironment, we here establish different Ndrg4 knockout (Ndrg4-/- ) CRC models and an indirect co-culture of primary enteric nervous system (ENS) cells and intestinal organoids to identify whether the ENS, via NDRG4, affects intestinal tumorigenesis. Linking immunostainings and gastrointestinal motility (GI) assays, we show that the absence of Ndrg4 does not trigger any functional or morphological GI abnormalities. However, combining in vivo, in vitro, and quantitative proteomics data, we uncover that Ndrg4 knockdown is associated with enlarged intestinal adenoma development and that organoid growth is boosted by the Ndrg4-/- ENS cell secretome, which is enriched for Nidogen-1 (Nid1) and Fibulin-2 (Fbln2). Moreover, NID1 and FBLN2 are expressed in enteric neurons, enhance migration capacities of CRC cells, and are enriched in human CRC secretomes. Hence, we provide evidence that the ENS, via loss of Ndrg4, is involved in colorectal pathogenesis and that ENS-derived Nidogen-1 and Fibulin-2 enhance colorectal carcinogenesis.

Dysregulation of the NRG1/ERBB pathway causes a developmental disorder with gastrointestinal dysmotility in humans.

Hirschsprung disease (HSCR) is the most frequent developmental anomaly of the enteric nervous system, with an incidence of 1 in 5000 live births. Chronic intestinal pseudo-obstruction (CIPO) is less frequent and classified as neurogenic or myogenic. Isolated HSCR has an oligogenic inheritance with RET as the major disease-causing gene, while CIPO is genetically heterogeneous, caused by mutations in smooth muscle-specific genes. Here, we describe a series of patients with developmental disorders including gastrointestinal dysmotility, and investigate the underlying molecular bases. Trio-exome sequencing led to the identification of biallelic variants in ERBB3 and ERBB2 in 8 individuals variably associating HSCR, CIPO, peripheral neuropathy, and arthrogryposis. Thorough gut histology revealed aganglionosis, hypoganglionosis, and intestinal smooth muscle abnormalities. The cell type-specific ErbB3 and ErbB2 function was further analyzed in mouse single-cell RNA sequencing data and in a conditional ErbB3-deficient mouse model, revealing a primary role for ERBB3 in enteric progenitors. The consequences of the identified variants were evaluated using quantitative real-time PCR (RT-qPCR) on patient-derived fibroblasts or immunoblot assays on Neuro-2a cells overexpressing WT or mutant proteins, revealing either decreased expression or altered phosphorylation of the mutant receptors. Our results demonstrate that dysregulation of ERBB3 or ERBB2 leads to a broad spectrum of developmental anomalies, including intestinal dysmotility.

TALPID3/KIAA0586 Regulates Multiple Aspects of Neuromuscular Patterning During Gastrointestinal Development in Animal Models and Human.

TALPID3/KIAA0586 is an evolutionary conserved protein, which plays an essential role in protein trafficking. Its role during gastrointestinal (GI) and enteric nervous system (ENS) development has not been studied previously. Here, we analyzed chicken, mouse and human embryonic GI tissues with TALPID3 mutations. The GI tract of TALPID3 chicken embryos was shortened and malformed. Histologically, the gut smooth muscle was mispatterned and enteric neural crest cells were scattered throughout the gut wall. Analysis of the Hedgehog pathway and gut extracellular matrix provided causative reasons for these defects. Interestingly, chicken intra-species grafting experiments and a conditional knockout mouse model showed that ENS formation did not require TALPID3, but was dependent on correct environmental cues. Surprisingly, the lack of TALPID3 in enteric neural crest cells (ENCC) affected smooth muscle and epithelial development in a non-cell-autonomous manner. Analysis of human gut fetal tissues with a KIAA0586 mutation showed strikingly similar findings compared to the animal models demonstrating conservation of TALPID3 and its necessary role in human GI tract development and patterning.

Using Out-of-Batch Reference Populations to Improve Untargeted Metabolomics for Screening Inborn Errors of Metabolism.

Untargeted metabolomics is an emerging technology in the laboratory diagnosis of inborn errors of metabolism (IEM). Analysis of a large number of reference samples is crucial for correcting variations in metabolite concentrations that result from factors, such as diet, age, and gender in order to judge whether metabolite levels are abnormal. However, a large number of reference samples requires the use of out-of-batch samples, which is hampered by the semi-quantitative nature of untargeted metabolomics data, i.e., technical variations between batches. Methods to merge and accurately normalize data from multiple batches are urgently needed. Based on six metrics, we compared the existing normalization methods on their ability to reduce the batch effects from nine independently processed batches. Many of those showed marginal performances, which motivated us to develop Metchalizer, a normalization method that uses 10 stable isotope-labeled internal standards and a mixed effect model. In addition, we propose a regression model with age and sex as covariates fitted on reference samples that were obtained from all nine batches. Metchalizer applied on log-transformed data showed the most promising performance on batch effect removal, as well as in the detection of 195 known biomarkers across 49 IEM patient samples and performed at least similar to an approach utilizing 15 within-batch reference samples. Furthermore, our regression model indicates that 6.5-37% of the considered features showed significant age-dependent variations. Our comprehensive comparison of normalization methods showed that our Log-Metchalizer approach enables the use out-of-batch reference samples to establish clinically-relevant reference values for metabolite concentrations. These findings open the possibilities to use large scale out-of-batch reference samples in a clinical setting, increasing the throughput and detection accuracy.

A complementary study approach unravels novel players in the pathoetiology of Hirschsprung disease.

Hirschsprung disease (HSCR, OMIM 142623) involves congenital intestinal obstruction caused by dysfunction of neural crest cells and their progeny during enteric nervous system (ENS) development. HSCR is a multifactorial disorder; pathogenetic variants accounting for disease phenotype are identified only in a minority of cases, and the identification of novel disease-relevant genes remains challenging. In order to identify and to validate a potential disease-causing relevance of novel HSCR candidate genes, we established a complementary study approach, combining whole exome sequencing (WES) with transcriptome analysis of murine embryonic ENS-related tissues, literature and database searches, in silico network analyses, and functional readouts using candidate gene-specific genome-edited cell clones. WES datasets of two patients with HSCR and their non-affected parents were analysed, and four novel HSCR candidate genes could be identified: ATP7A, SREBF1, ABCD1 and PIAS2. Further rare variants in these genes were identified in additional HSCR patients, suggesting disease relevance. Transcriptomics revealed that these genes are expressed in embryonic and fetal gastrointestinal tissues. Knockout of these genes in neuronal cells demonstrated impaired cell differentiation, proliferation and/or survival. Our approach identified and validated candidate HSCR genes and provided further insight into the underlying pathomechanisms of HSCR.

Goldberg-Shprintzen syndrome is determined by the absence, or reduced expression levels, of KIFBP.

Goldberg-Shprintzen syndrome (GOSHS) is caused by loss of function variants in the kinesin binding protein gene (KIFBP). However, the phenotypic range of this syndrome is wide, indicating that other factors may play a role. To date, 37 patients with GOSHS have been reported. Here, we document nine new patients with variants in KIFBP: seven with nonsense variants and two with missense variants. To our knowledge, this is the first time that missense variants have been reported in GOSHS. We functionally investigated the effect of the variants identified, in an attempt to find a genotype-phenotype correlation. We also determined whether common Hirschsprung disease (HSCR)-associated single nucleotide polymorphisms (SNPs), could explain the presence of HSCR in GOSHS. Our results showed that the missense variants led to reduced expression of KIFBP, while the truncating variants resulted in lack of protein. However, no correlation was found between the severity of GOSHS and the location of the variants. We were also unable to find a correlation between common HSCR-associated SNPs, and HSCR development in GOSHS. In conclusion, we show that reduced, as well as lack of KIFBP expression can lead to GOSHS, and our results suggest that a threshold expression of KIFBP may modulate phenotypic variability of the disease.

Inhibition of ROCK signaling pathway accelerates enteric neural crest cell-based therapy after transplantation in a rat hypoganglionic model.

BACKGROUND: Hirschsprung's disease (HSCR) is a congenital gastrointestinal disorder, characterized by enteric ganglia absence in part or entire of the colon, due to abnormal colonization and migration of enteric neural crest cells (ENCCs) during development. Currently, besides surgery which is the main therapy for HSCR, the potential of stem cell-based transplantation was investigated as an alternative option. Although promising, it has limitations, including poor survival, differentiation, and migration of the grafted cells. We hypothesized that modulation of extracellular factors during transplantation could promote ENCCs proliferation and migration, leading to increased transplantation efficiency. Considering that the RhoA/ROCK pathway is highly involved in cytoskeletal dynamics and neurite growth, our study explored the effect of inhibition of this pathway to improve the success of ENCCs transplantation. METHODS: Enteric neural crest cells were isolated from rat embryos and labeled with a GFP-tag. Cell viability, apoptosis, differentiation, and migration assays were performed with and without RhoA/ROCK inhibition. Labeled ENCCs were transplanted into the muscle layer of an induced hypoganglionic rat model followed by intraperitoneal injections of ROCK inhibitor. The transplanted segments were collected 3 weeks after for histological analysis. KEY RESULTS: Our results showed that inhibition of ROCK increased viable cell number, differentiation, and migration of ENCCs in vitro. Moreover, transplantation of labeled ENCCs into the hypoganglionic model showed enhanced distribution of grafted ENCCs, upon treatment with ROCK inhibitor. CONCLUSIONS AND INFERENCES: ROCK inhibitors influence ENCCs growth and migration in vitro and in vivo, and should be considered to improve the efficiency of ENCCs transplantation.

ERNICA guidelines for the management of rectosigmoid Hirschsprung's disease.

BACKGROUND: Hirschsprung's disease (HSCR) is a serious congenital bowel disorder with a prevalence of 1/5000. Currently, there is a lack of systematically developed guidelines to assist clinical decision-making regarding diagnostics and management. AIMS: This guideline aims to cover the diagnostics and management of rectosigmoid HSCR up to adulthood. It aims to describe the preferred approach of ERNICA, the European Reference Network for rare inherited and congenital digestive disorders. METHODS: Recommendations within key topics covering the care pathway for rectosigmoid HSCR were developed by an international workgroup of experts from 8 European countries within ERNICA European Reference Network from the disciplines of surgery, medicine, histopathology, microbiology, genetics, and patient organization representatives. Recommendation statements were based on a comprehensive review of the available literature and expert consensus. AGREE II and GRADE approaches were used during development. Evidence levels and levels of agreement are noted. RESULTS: Thirty-three statements within 9 key areas were generated. Most recommendations were based on expert opinion. CONCLUSION: In rare or low-prevalence diseases such as HSCR, there remains limited availability of high-quality clinical evidence. Consensus-based guidelines for care are presented.

Infantile hypertrophic pyloric stenosis in patients with esophageal atresia.

BACKGROUND: Patients born with esophageal atresia (EA) have a higher incidence of infantile hypertrophic pyloric stenosis (IHPS), suggestive of a relationship. A shared etiology makes sense from a developmental perspective as both affected structures are foregut derived. A genetic component has been described for both conditions as single entities and EA and IHPS are variable components in several monogenetic syndromes. We hypothesized that defects disturbing foregut morphogenesis are responsible for this combination of malformations. METHODS: We investigated the genetic variation of 15 patients with both EA and IHPS with unaffected parents using exome sequencing and SNP array-based genotyping, and compared the results to mouse transcriptome data of the developing foregut. RESULTS: We did not identify putatively deleterious de novo mutations or recessive variants. However, we detected rare inherited variants in EA or IHPS disease genes or in genes important in foregut morphogenesis, expressed at the proper developmental time-points. Two pathways were significantly enriched (p < 1 × 10-5 ): proliferation and differentiation of smooth muscle cells and self-renewal of satellite cells. CONCLUSIONS: None of our findings could fully explain the combination of abnormalities on its own, which makes complex inheritance the most plausible genetic explanation, most likely in combination with mechanical and/or environmental factors. As we did not find one defining monogenetic cause for the EA/IHPS phenotype, the impact of the corrective surgery could should be further investigated.

Zebrafish: A Model Organism for Studying Enteric Nervous System Development and Disease.

The Enteric Nervous System (ENS) is a large network of enteric neurons and glia that regulates various processes in the gastrointestinal tract including motility, local blood flow, mucosal transport and secretion. The ENS is derived from stem cells coming from the neural crest that migrate into and along the primitive gut. Defects in ENS establishment cause enteric neuropathies, including Hirschsprung disease (HSCR), which is characterized by an absence of enteric neural crest cells in the distal part of the colon. In this review, we discuss the use of zebrafish as a model organism to study the development of the ENS. The accessibility of the rapidly developing gut in zebrafish embryos and larvae, enables in vivo visualization of ENS development, peristalsis and gut transit. These properties make the zebrafish a highly suitable model to bring new insights into ENS development, as well as in HSCR pathogenesis. Zebrafish have already proven fruitful in studying ENS functionality and in the validation of novel HSCR risk genes. With the rapid advancements in gene editing techniques and their unique properties, research using zebrafish as a disease model, will further increase our understanding on the genetics underlying HSCR, as well as possible treatment options for this disease.