Interleukin-22 orchestrates a pathological endoplasmic reticulum stress response transcriptional programme in colonic epithelial cells.

OBJECTIVE: The functional role of interleukin-22 (IL22) in chronic inflammation is controversial, and mechanistic insights into how it regulates target tissue are lacking. In this study, we evaluated the functional role of IL22 in chronic colitis and probed mechanisms of IL22-mediated regulation of colonic epithelial cells. DESIGN: To investigate the functional role of IL22 in chronic colitis and how it regulates colonic epithelial cells, we employed a three-dimentional mini-gut epithelial organoid system, in vivo disease models and transcriptomic datasets in human IBD. RESULTS: As well as inducing transcriptional modules implicated in antimicrobial responses, IL22 also coordinated an endoplasmic reticulum (ER) stress response transcriptional programme in colonic epithelial cells. In the colon of patients with active colonic Crohn's disease (CD), there was enrichment of IL22-responsive transcriptional modules and ER stress response modules. Strikingly, in an IL22-dependent model of chronic colitis, targeting IL22 alleviated colonic epithelial ER stress and attenuated colitis. Pharmacological modulation of the ER stress response similarly impacted the severity of colitis. In patients with colonic CD, antibody blockade of IL12p40, which simultaneously blocks IL12 and IL23, the key upstream regulator of IL22 production, alleviated the colonic epithelial ER stress response. CONCLUSIONS: Our data challenge perceptions of IL22 as a predominantly beneficial cytokine in IBD and provide novel insights into the molecular mechanisms of IL22-mediated pathogenicity in chronic colitis. Targeting IL22-regulated pathways and alleviating colonic epithelial ER stress may represent promising therapeutic strategies in patients with colitis. TRIAL REGISTRATION NUMBER: NCT02749630.

Quantitative Systems Toxicology Modeling To Address Key Safety Questions in Drug Development: A Focus of the TransQST Consortium.

Omics data have been increasingly generated with limited demonstrated value in drug safety assessment. The TransQST consortium was launched to use omics and other data in mechanistic-based quantitative systems toxicology (QST) models to evaluate their potential use in species translation.

Host-microbe interaction in the gastrointestinal tract.

The gastrointestinal tract is a highly complex organ in which multiple dynamic physiological processes are tightly coordinated while interacting with a dense and extremely diverse microbial population. From establishment in early life, through to host-microbe symbiosis in adulthood, the gut microbiota plays a vital role in our development and health. The effect of the microbiota on gut development and physiology is highlighted by anatomical and functional changes in germ-free mice, affecting the gut epithelium, immune system and enteric nervous system. Microbial colonisation promotes competent innate and acquired mucosal immune systems, epithelial renewal, barrier integrity, and mucosal vascularisation and innervation. Interacting or shared signalling pathways across different physiological systems of the gut could explain how all these changes are coordinated during postnatal colonisation, or after the introduction of microbiota into germ-free models. The application of cell-based in-vitro experimental systems and mathematical modelling can shed light on the molecular and signalling pathways which regulate the development and maintenance of homeostasis in the gut and beyond.

Cell proliferation within small intestinal crypts is the principal driving force for cell migration on villi.

The functional integrity of the intestinal epithelial barrier relies on tight coordination of cell proliferation and migration, with failure to regulate these processes resulting in disease. It is not known whether cell proliferation is sufficient to drive epithelial cell migration during homoeostatic turnover of the epithelium. Nor is it known precisely how villus cell migration is affected when proliferation is perturbed. Some reports suggest that proliferation and migration may not be related while other studies support a direct relationship. We used established cell-tracking methods based on thymine analog cell labeling and developed tailored mathematical models to quantify cell proliferation and migration under normal conditions and when proliferation is reduced and when it is temporarily halted. We found that epithelial cell migration velocities along the villi are coupled to cell proliferation rates within the crypts in all conditions. Furthermore, halting and resuming proliferation results in the synchronized response of cell migration on the villi. We conclude that cell proliferation within the crypt is the primary force that drives cell migration along the villus. This methodology can be applied to interrogate intestinal epithelial dynamics and characterize situations in which processes involved in cell turnover become uncoupled, including pharmacological treatments and disease models.-Parker, A., Maclaren, O. J., Fletcher, A. G., Muraro, D., Kreuzaler, P. A., Byrne, H. M., Maini, P. K., Watson, A. J. M., Pin, C. Cell proliferation within small intestinal crypts is the principal driving force for cell migration on villi.

A hierarchical Bayesian model for understanding the spatiotemporal dynamics of the intestinal epithelium.

Our work addresses two key challenges, one biological and one methodological. First, we aim to understand how proliferation and cell migration rates in the intestinal epithelium are related under healthy, damaged (Ara-C treated) and recovering conditions, and how these relations can be used to identify mechanisms of repair and regeneration. We analyse new data, presented in more detail in a companion paper, in which BrdU/IdU cell-labelling experiments were performed under these respective conditions. Second, in considering how to more rigorously process these data and interpret them using mathematical models, we use a probabilistic, hierarchical approach. This provides a best-practice approach for systematically modelling and understanding the uncertainties that can otherwise undermine the generation of reliable conclusions-uncertainties in experimental measurement and treatment, difficult-to-compare mathematical models of underlying mechanisms, and unknown or unobserved parameters. Both spatially discrete and continuous mechanistic models are considered and related via hierarchical conditional probability assumptions. We perform model checks on both in-sample and out-of-sample datasets and use them to show how to test possible model improvements and assess the robustness of our conclusions. We conclude, for the present set of experiments, that a primarily proliferation-driven model suffices to predict labelled cell dynamics over most time-scales.

Interactions of Salmonella enterica subspecies enterica serovar Typhimurium with gut bacteria.

The aim of this study was to evaluate the impact of the gut microbiota on the growth and survival of S. Typhimurium. This was tested in two-species co-cultures and in mixed cultures with a simplified gut model microbiota. Subsequently, interactions between S. Typhimurium and human faecal bacteria were quantified in both batch and continuous culture systems simulating the human colon. The exponential growth of S. Typhimurium was halted when the population of Escherichia coli reached the maximum population density in a two-compartment co-culture system where the two species were separated by a 0.45 µm pore membrane. Furthermore, the growth of some gut bacteria such as Lactobacillus gasseri and Bifidobacterium bifidum was inhibited by the presence of S. Typhimurium in the other compartment. The survival of S. Typhimurium was severely affected in mixed batch cultures with human faecal samples; a reduction of 10(3)-10(4) cfu/ml in the concentration of S. Typhimurium was observed in these cultures. However, no effect on S. Typhimurium survival was observed in mixed batch cultures with a simplified gut model microbiota under the same conditions. The effect of human faecal samples on S. Typhimurium in a three-stage continuous culture was different to that obtained in batch cultures; its growth rather than survival was affected under these conditions. S. Typhimurium growth was inhibited, and the bacterium was therefore eliminated by the continuous flow of the medium. Depending upon culturing conditions, the gut microbiota caused either growth inhibition, inactivation or did not affect S. Typhimurium.

γδ T-cell-deficient mice show alterations in mucin expression, glycosylation, and goblet cells but maintain an intact mucus layer.

Intestinal homeostasis is maintained by a hierarchy of immune defenses acting in concert to minimize contact between luminal microorganisms and the intestinal epithelial cell surface. The intestinal mucus layer, covering the gastrointestinal tract epithelial cells, contributes to mucosal homeostasis by limiting bacterial invasion. In this study, we used γδ T-cell-deficient (TCRδ(-/-)) mice to examine whether and how γδ T-cells modulate the properties of the intestinal mucus layer. Increased susceptibility of TCRδ(-/-) mice to dextran sodium sulfate (DSS)-induced colitis is associated with a reduced number of goblet cells. Alterations in the number of goblet cells and crypt lengths were observed in the small intestine and colon of TCRδ(-/-) mice compared with C57BL/6 wild-type (WT) mice. Addition of keratinocyte growth factor to small intestinal organoid cultures from TCRδ(-/-) mice showed a marked increase in crypt growth and in both goblet cell number and redistribution along the crypts. There was no apparent difference in the thickness or organization of the mucus layer between TCRδ(-/-) and WT mice, as measured in vivo. However, γδ T-cell deficiency led to reduced sialylated mucins in association with increased gene expression of gel-secreting Muc2 and membrane-bound mucins, including Muc13 and Muc17. Collectively, these data provide evidence that γδ T cells play an important role in the maintenance of mucosal homeostasis by regulating mucin expression and promoting goblet cell function in the small intestine.

Mining the "glycocode"--exploring the spatial distribution of glycans in gastrointestinal mucin using force spectroscopy.

Mucins are the main components of the gastrointestinal mucus layer. Mucin glycosylation is critical to most intermolecular and intercellular interactions. However, due to the highly complex and heterogeneous mucin glycan structures, the encoded biological information remains largely encrypted. Here we have developed a methodology based on force spectroscopy to identify biologically accessible glycoepitopes in purified porcine gastric mucin (pPGM) and purified porcine jejunal mucin (pPJM). The binding specificity of lectins Ricinus communis agglutinin I (RCA), peanut (Arachis hypogaea) agglutinin (PNA), Maackia amurensis lectin II (MALII), and Ulex europaeus agglutinin I (UEA) was utilized in force spectroscopy measurements to quantify the affinity and spatial distribution of their cognate sugars at the molecular scale. Binding energy of 4, 1.6, and 26 aJ was determined on pPGM for RCA, PNA, and UEA. Binding was abolished by competition with free ligands, demonstrating the validity of the affinity data. The distributions of the nearest binding site separations estimated the number of binding sites in a 200-nm mucin segment to be 4 for RCA, PNA, and UEA, and 1.8 for MALII. Binding site separations were affected by partial defucosylation of pPGM. Furthermore, we showed that this new approach can resolve differences between gastric and jejunum mucins.

Loss of culturability of Salmonella enterica subsp. enterica serovar Typhimurium upon cell-cell contact with human fecal bacteria.

Loss of culturability of Salmonella enterica subsp. enterica serovar Typhimurium has been observed in mixed cultures with anaerobic fecal bacteria under conditions that allow local interaction between cells, such as cell contact. A reduction of a population of culturable S. Typhimurium on the order of ∼10(4) to 10(5) CFU/ml was observed in batch anaerobic mixed cultures with fecal samples from different human donors. Culturability was not affected either in supernatants collected at several times from fecal cultures, when separated from fecal bacteria by a membrane of 0.45-µm pore size, or when in contact with inactivated fecal bacterial cells. Loss of culturability kinetics was characterized by a sharp reduction of several logarithmic units followed by a pronounced tail. A mathematical model was developed to describe the rate of loss of culturability as a function of the frequency of encounters between populations and the probability of inactivation after encounter. The model term F(S · F)(1/2) quantifies the effect of the concentration of both populations, fecal bacteria (F) and S. Typhimurium (S), on the loss of culturability of S. Typhimurium by cell contact with fecal bacteria. When the value of F(S · F)(1/2) decreased below ca. 10(15) (CFU/ml)(2), the frequency of encounters sharply decreased, leading to the deceleration of the inactivation rate and to the tailing off of the S. Typhimurium population. The probability of inactivation after encounter, P, was constant, with an estimated value of ∼10(-5) for all data sets. P might be characteristic of the mechanism of growth inhibition after a cell encounter.

Non-essential genes form the hubs of genome scale protein function and environmental gene expression networks in Salmonella enterica serovar Typhimurium.

BACKGROUND: Salmonella Typhimurium is an important pathogen of human and animals. It shows a broad growth range and survives in harsh conditions. The aim of this study was to analyze transcriptional responses to a number of growth and stress conditions as well as the relationship of metabolic pathways and/or cell functions at the genome-scale-level by network analysis, and further to explore whether highly connected genes (hubs) in these networks were essential for growth, stress adaptation and virulence. RESULTS: De novo generated as well as published transcriptional data for 425 selected genes under a number of growth and stress conditions were used to construct a bipartite network connecting culture conditions and significantly regulated genes (transcriptional network). Also, a genome scale network was constructed for strain LT2. The latter connected genes with metabolic pathways and cellular functions. Both networks were shown to belong to the family of scale-free networks characterized by the presence of highly connected nodes or hubs which are genes whose transcription is regulated when responding to many of the assayed culture conditions or genes encoding products involved in a high number of metabolic pathways and cell functions.The five genes with most connections in the transcriptional network (wraB, ygaU, uspA, cbpA and osmC) and in the genome scale network (ychN, siiF (STM4262), yajD, ybeB and dcoC) were selected for mutations, however mutagenesis of ygaU and ybeB proved unsuccessful. No difference between mutants and the wild type strain was observed during growth at unfavorable temperatures, pH values, NaCl concentrations and in the presence of H2O2. Eight mutants were evaluated for virulence in C57/BL6 mice and none differed from the wild type strain. Notably, however, deviations of phenotypes with respect to the wild type were observed when combinations of these genes were deleted. CONCLUSION: Network analysis revealed the presence of hubs in both transcriptional and functional networks of S. Typhimurium. Hubs theoretically confer higher resistance to random mutation but a greater susceptibility to directed attacks, however, we found that genes that formed hubs were dispensable for growth, stress adaptation and virulence, suggesting that evolution favors non-essential genes as main connectors in cellular networks.

Homeostasis, injury, and recovery dynamics at multiple scales in a self-organizing mouse intestinal crypt.

The maintenance of the functional integrity of the intestinal epithelium requires a tight coordination between cell production, migration, and shedding along the crypt-villus axis. Dysregulation of these processes may result in loss of the intestinal barrier and disease. With the aim of generating a more complete and integrated understanding of how the epithelium maintains homeostasis and recovers after injury, we have built a multi-scale agent-based model (ABM) of the mouse intestinal epithelium. We demonstrate that stable, self-organizing behaviour in the crypt emerges from the dynamic interaction of multiple signalling pathways, such as Wnt, Notch, BMP, ZNRF3/RNF43, and YAP-Hippo pathways, which regulate proliferation and differentiation, respond to environmental mechanical cues, form feedback mechanisms, and modulate the dynamics of the cell cycle protein network. The model recapitulates the crypt phenotype reported after persistent stem cell ablation and after the inhibition of the CDK1 cycle protein. Moreover, we simulated 5-fluorouracil (5-FU)-induced toxicity at multiple scales starting from DNA and RNA damage, which disrupts the cell cycle, cell signalling, proliferation, differentiation, and migration and leads to loss of barrier integrity. During recovery, our in silico crypt regenerates its structure in a self-organizing, dynamic fashion driven by dedifferentiation and enhanced by negative feedback loops. Thus, the model enables the simulation of xenobiotic-, in particular chemotherapy-, induced mechanisms of intestinal toxicity and epithelial recovery. Overall, we present a systems model able to simulate the disruption of molecular events and its impact across multiple levels of epithelial organization and demonstrate its application to epithelial research and drug development.

A dynamic model of the intestinal epithelium integrates multiple sources of preclinical data and enables clinical translation of drug-induced toxicity.

We have built a quantitative systems toxicology modeling framework focused on the early prediction of oncotherapeutic-induced clinical intestinal adverse effects. The model describes stem and progenitor cell dynamics in the small intestinal epithelium and integrates heterogeneous epithelial-related processes, such as transcriptional profiles, citrulline kinetics, and probability of diarrhea. We fitted a mouse-specific version of the model to quantify doxorubicin and 5-fluorouracil (5-FU)-induced toxicity, which included pharmacokinetics and 5-FU metabolism and assumed that both drugs led to cell cycle arrest and apoptosis in stem cells and proliferative progenitors. The model successfully recapitulated observations in mice regarding dose-dependent disruption of proliferation which could lead to villus shortening, decrease of circulating citrulline, increased diarrhea risk, and transcriptional induction of the p53 pathway. Using a human-specific epithelial model, we translated the cytotoxic activity of doxorubicin and 5-FU quantified in mice into human intestinal injury and predicted with accuracy clinical diarrhea incidence. However, for gefitinib, a specific-molecularly targeted therapy, the mice failed to reproduce epithelial toxicity at exposures much higher than those associated with clinical diarrhea. This indicates that, regardless of the translational modeling approach, preclinical experimental settings have to be suitable to quantify drug-induced clinical toxicity with precision at the structural scale of the model. Our work demonstrates the usefulness of translational models at early stages of the drug development pipeline to predict clinical toxicity and highlights the importance of understanding cross-settings differences in toxicity when building these approaches.

Lag phase is a distinct growth phase that prepares bacteria for exponential growth and involves transient metal accumulation.

Lag phase represents the earliest and most poorly understood stage of the bacterial growth cycle. We developed a reproducible experimental system and conducted functional genomic and physiological analyses of a 2-h lag phase in Salmonella enterica serovar Typhimurium. Adaptation began within 4 min of inoculation into fresh LB medium with the transient expression of genes involved in phosphate uptake. The main lag-phase transcriptional program initiated at 20 min with the upregulation of 945 genes encoding processes such as transcription, translation, iron-sulfur protein assembly, nucleotide metabolism, LPS biosynthesis, and aerobic respiration. ChIP-chip revealed that RNA polymerase was not "poised" upstream of the bacterial genes that are rapidly induced at the beginning of lag phase, suggesting a mechanism that involves de novo partitioning of RNA polymerase to transcribe 522 bacterial genes within 4 min of leaving stationary phase. We used inductively coupled plasma mass spectrometry (ICP-MS) to discover that iron, calcium, and manganese are accumulated by S. Typhimurium during lag phase, while levels of cobalt, nickel, and sodium showed distinct growth-phase-specific patterns. The high concentration of iron during lag phase was associated with transient sensitivity to oxidative stress. The study of lag phase promises to identify the physiological and regulatory processes responsible for adaptation to new environments.

Systems Modeling to Quantify Safety Risks in Early Drug Development: Using Bifurcation Analysis and Agent-Based Modeling as Examples.

Quantitative Systems Toxicology (QST) models, recapitulating pharmacokinetics and mechanism of action together with the organic response at multiple levels of biological organization, can provide predictions on the magnitude of injury and recovery dynamics to support study design and decision-making during drug development. Here, we highlight the application of QST models to predict toxicities of cancer treatments, such as cytopenia(s) and gastrointestinal adverse effects, where narrow therapeutic indexes need to be actively managed. The importance of bifurcation analysis is demonstrated in QST models of hematologic toxicity to understand how different regions of the parameter space generate different behaviors following cancer treatment, which results in asymptotically stable predictions, yet highly irregular for specific schedules, or oscillating predictions of blood cell levels. In addition, an agent-based model of the intestinal crypt was used to simulate how the spatial location of the injury within the crypt affects the villus disruption severity. We discuss the value of QST modeling approaches to support drug development and how they align with technological advances impacting trial design including patient selection, dose/regimen selection, and ultimately patient safety.

Light-Induced Smooth Endoplasmic Reticulum Rearrangement in a Unique Interlaced Compartmental Pattern in Macaca mulatta RPE.

Purpose: To investigate light-induced modifications of the smooth endoplasmic reticulum of the RPE in primates. Methods: Eyes of three terminally anesthetized Rhesus monkeys were exposed to 5000 lux for 10 minutes or kept in the dark. Transmission electron microscopy and electron tomography were conducted on small fragments of retina sampled from different regions of the retina. Results: RPE cells smooth endoplasmic reticulum shows a previously unknown arrangement characterized by an interlaced compartmental pattern (ICP). Electron tomograms and 3D-modelling demonstrated that the smooth endoplasmic reticulum with an ICP (ICPSER) consisted of four parallel, independent and interwoven networks of tubules arranged as interconnected coiled coils. Its architecture realized a compact labyrinthine structure of tightly packed tubules stabilized by intertubular filamentous tethers. On average, the ICPSER is present in about 14.6% of RPE cells. Although ICPSER was preferentially found in cells located in the peripheral and in the para/perifoveal retina, ICPSER cells significantly increased in number upon light exposure in the para/perifovea and in the fovea. Conclusions: An ICPSER is apparently a unique feature to primate RPE. Its rapid appearance in the area centralis of the retina upon light exposure suggests a function related to the foveate structure of primate retina or to the diurnal habits of animals that may require additional protection from photo-oxidation or enhanced requests of visual pigments regeneration.

Modeling the lag period and exponential growth of Listeria monocytogenes under conditions of fluctuating temperature and water activity values.

The dynamic model for the growth of a bacterial population described by Baranyi and Roberts (J. Baranyi and T. A. Roberts, Int. J. Food Microbiol. 23:277-294, 1994) was applied to model the lag period and exponential growth of Listeria monocytogenes under conditions of fluctuating temperature and water activity (a(w)) values. To model the duration of the lag phase, the dependence of the parameter h(0), which quantifies the amount of work done during the lag period, on the previous and current environmental conditions was determined experimentally. This parameter depended not only on the magnitude of the change between the previous and current environmental conditions but also on the current growth conditions. In an exponentially growing population, any change in the environment requiring a certain amount of work to adapt to the new conditions initiated a lag period that lasted until that work was finished. Observations for several scenarios in which exponential growth was halted by a sudden change in the temperature and/or a(w) were in good agreement with predictions. When a population already in a lag period was subjected to environmental fluctuations, the system was reset with a new lag phase. The work to be done during the new lag phase was estimated to be the workload due to the environmental change plus the unfinished workload from the uncompleted previous lag phase.

Effects of carbon dioxide on growth of proteolytic Clostridium botulinum, its ability to produce neurotoxin, and its transcriptome.

The antimicrobial gas carbon dioxide is frequently used in modified atmosphere packaging. In the present study, the effects of CO2 (10 to 70%, vol/vol) on gene expression (measured using quantitative reverse transcription-PCR and a whole-genome DNA microarray) and neurotoxin formation (measured using an enzyme-linked immunosorbent assay [ELISA]) by proteolytic Clostridium botulinum type A1 strain ATCC 3502 were studied during the growth cycle. Interestingly, in marked contrast to the situation with nonproteolytic C. botulinum types B and E, CO2 had little effect on any of these parameters. At all CO2 concentrations, relative expression of neurotoxin cluster genes peaked in the transition between exponential and stationary phases, with evidence of a second rise in expression in late stationary phase. Microarray analysis enabled identification of coding sequences whose expression profiles matched those of the neurotoxin cluster. Further research is needed to determine whether these are connected to neurotoxin formation or are merely growth phase associated.

Macrophage migration inhibitory factor plays a role in the regulation of microfold (M) cell-mediated transport in the gut.

It has been shown previously that certain bacteria rapidly (3 h) up-regulated in vivo microfold cell (M cell)-mediated transport of Ag across the follicle-associated epithelium of intestinal Peyer's patch. Our aim was to determine whether soluble mediators secreted following host-bacteria interaction were involved in this event. A combination of proteomics and immunohistochemical analyses was used to identify molecules produced in the gut in response to bacterial challenge in vivo; their effects were then tested on human intestinal epithelial cells in vitro. Macrophage migration inhibitory factor (MIF) was the only cytokine produced rapidly after in vivo bacterial challenge by CD11c(+) cells located beneath the M cell-rich area of the follicle-associated epithelium of the Peyer's patch. Subsequently, in vitro experiments conducted using human Caco-2 cells showed that, within hours, MIF induced the appearance of cells that showed temperature-dependent transport of microparticles and M cell-specific bacterium Vibrio cholerae, and acquired biochemical features of M cells. Furthermore, using an established in vitro human M cell model, we showed that anti-MIF Ab blocked Raji B cell-mediated conversion of Caco-2 cells into Ag-sampling cells. Finally, we report that MIF(-/-) mice, in contrast to wild-type mice, failed to show increased M cell-mediated transport following in vivo bacterial challenge. These data show that MIF plays a role in M cell-mediated transport, and cross-talk between bacteria, gut epithelium, and immune system is instrumental in regulating key functions of the gut, including M cell-mediated Ag sampling.

Importance of Stability Analysis When Using Nonlinear Semimechanistic Models to Describe Drug-Induced Hematotoxicity.

Stability analysis, often overlooked in pharmacometrics, is essential to explore dynamical systems. The model developed by Friberg et al.1 to describe drug-induced hematotoxicity is widely used to support decisions across drug development, and parameter values are often identified from observed blood counts. We use stability analysis to study the parametric dependence of stable and unstable solutions of several Friberg-type models and highlight the risks associated with system instability in the context of nonlinear mixed effects modeling. We emphasize the consequences of unstable solutions on prediction performance by demonstrating nonbiological system behaviors in a real case study of drug-induced thrombocytopenia. Ultimately, we provide simple criteria for identifying parameters associated with stable solutions of Friberg-type models. For instance, in the original Friberg model, we find that stability depends only on the parameter that governs the feedback from peripheral cells to progenitors and provide the exact range of values that results in stable solutions.

Enteric glia as a source of neural progenitors in adult zebrafish.

The presence and identity of neural progenitors in the enteric nervous system (ENS) of vertebrates is a matter of intense debate. Here, we demonstrate that the non-neuronal ENS cell compartment of teleosts shares molecular and morphological characteristics with mammalian enteric glia but cannot be identified by the expression of canonical glial markers. However, unlike their mammalian counterparts, which are generally quiescent and do not undergo neuronal differentiation during homeostasis, we show that a relatively high proportion of zebrafish enteric glia proliferate under physiological conditions giving rise to progeny that differentiate into enteric neurons. We also provide evidence that, similar to brain neural stem cells, the activation and neuronal differentiation of enteric glia are regulated by Notch signalling. Our experiments reveal remarkable similarities between enteric glia and brain neural stem cells in teleosts and open new possibilities for use of mammalian enteric glia as a potential source of neurons to restore the activity of intestinal neural circuits compromised by injury or disease.