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.

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.

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.

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.

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.

Developing in vitro assays to transform gastrointestinal safety assessment: potential for microphysiological systems.

Drug-induced gastrointestinal toxicities (DI-GITs) are among the most common adverse events in clinical trials. High prevalence of DI-GIT has persisted among new drugs due in part to the lack of robust experimental tools to allow early detection or to guide optimization of safer molecules. Developing in vitro assays for the leading GI toxicities (nausea, vomiting, diarrhoea, constipation, and abdominal pain) will likely involve recapitulating complex physiological properties that require contributions from diverse cell/tissue types including epithelial, immune, microbiome, nerve, and muscle. While this stipulation may be beyond traditional 2D monocultures of intestinal cell lines, emerging 3D GI microtissues capture interactions between diverse cell and tissue types. These interactions give rise to microphysiologies fundamental to gut biology. For GI microtissues, organoid technology was the breakthrough that introduced intestinal stem cells with the capability of differentiating into each of the epithelial cell types and that self-organize into a multi-cellular tissue proxy with villus- and crypt-like domains. Recently, GI microtissues generated using miniaturized devices with microfluidic flow and cyclic peristaltic strain were shown to induce Caco2 cells to spontaneously differentiate into each of the principle intestinal epithelial cell types. Second generation models comprised of epithelial organoids or microtissues co-cultured with non-epithelial cell types can successfully reproduce cross-'tissue' functional interactions broadening the potential of these models to accurately study drug-induced toxicities. A new paradigm in which in vitro assays become an early part of GI safety assessment could be realized if microphysiological systems (MPS) are developed in alignment with drug-discovery needs. Herein, approaches for assessing GI toxicity of pharmaceuticals are reviewed and gaps are compared with capabilities of emerging GI microtissues (e.g., organoids, organ-on-a-chip, transwell systems) in order to provide perspective on the assay features needed for MPS models to be adopted for DI-GIT assessment.

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.

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.

Quantifying Drug-Induced Bone Marrow Toxicity Using a Novel Haematopoiesis Systems Pharmacology Model.

Haematological toxicity associated with cancer therapeutics is monitored by changes in blood cell count, and their primary effect is on proliferative progenitors in the bone marrow. Using observations in rat bone marrow and blood, we characterize a mathematical model that comprises cell proliferation and differentiation of the full haematopoietic phylogeny, with interacting feedback loops between lineages in homeostasis as well as following carboplatin exposure. We accurately predicted the temporal dynamics of several mature cell types related to carboplatin-induced bone marrow toxicity and identified novel insights into haematopoiesis. Our model confirms a significant degree of plasticity within bone marrow cells, with the number and type of both early progenitors and circulating cells affecting cell balance, via feedback mechanisms, through fate decisions of the multipotent progenitors. We also demonstrated cross-species translation of our predictions to patients, applying the same core model structure and considering differences in drug-dependent and physiology-dependent parameters.

Elevated apoptosis impairs epithelial cell turnover and shortens villi in TNF-driven intestinal inflammation.

The intestinal epithelial monolayer, at the boundary between microbes and the host immune system, plays an important role in the development of inflammatory bowel disease (IBD), particularly as a target and producer of pro-inflammatory TNF. Chronic overexpression of TNF leads to IBD-like pathology over time, but the mechanisms driving early pathogenesis events are not clear. We studied the epithelial response to inflammation by combining mathematical models with in vivo experimental models resembling acute and chronic TNF-mediated injury. We found significant villus atrophy with increased epithelial cell death along the crypt-villus axis, most dramatically at the villus tips, in both acute and chronic inflammation. In the acute model, we observed overexpression of TNF receptor I in the villus tip rapidly after TNF injection and concurrent with elevated levels of intracellular TNF and rapid shedding at the tip. In the chronic model, sustained villus atrophy was accompanied by a reduction in absolute epithelial cell turnover. Mathematical modelling demonstrated that increased cell apoptosis on the villus body explains the reduction in epithelial cell turnover along the crypt-villus axis observed in chronic inflammation. Cell destruction in the villus was not accompanied by changes in proliferative cell number or division rate within the crypt. Epithelial morphology and immunological changes in the chronic setting suggest a repair response to cell damage although the villus length is not recovered. A better understanding of how this state is further destabilised and results in clinical pathology resembling IBD will help identify suitable pathways for therapeutic intervention.

Human 3D Gastrointestinal Microtissue Barrier Function As a Predictor of Drug-Induced Diarrhea.

Drug-induced gastrointestinal toxicities (GITs) rank among the most common clinical side effects. Preclinical efforts to reduce incidence are limited by inadequate predictivity of in vitro assays. Recent breakthroughs in in vitro culture methods support intestinal stem cell maintenance and continual differentiation into the epithelial cell types resident in the intestine. These diverse cells self-assemble into microtissues with in vivo-like architecture. Here, we evaluate human GI microtissues grown in transwell plates that allow apical and/or basolateral drug treatment and 96-well throughput. Evaluation of assay utility focused on predictivity for diarrhea because this adverse effect correlates with intestinal barrier dysfunction which can be measured in GI microtissues using transepithelial electrical resistance (TEER). A validation set of widely prescribed drugs was assembled and tested for effects on TEER. When the resulting TEER inhibition potencies were adjusted for clinical exposure, a threshold was identified that distinguished drugs that induced clinical diarrhea from those that lack this liability. Microtissue TEER assay predictivity was further challenged with a smaller set of drugs whose clinical development was limited by diarrhea that was unexpected based on 1-month animal studies. Microtissue TEER accurately predicted diarrhea for each of these drugs. The label-free nature of TEER enabled repeated quantitation with sufficient precision to develop a mathematical model describing the temporal dynamics of barrier damage and recovery. This human 3D GI microtissue is the first in vitro assay with validated predictivity for diarrhea-inducing drugs. It should provide a platform for lead optimization and offers potential for dose schedule exploration.

Chronic TNFα-driven injury delays cell migration to villi in the intestinal epithelium.

The intestinal epithelium is a single layer of cells which provides the first line of defence of the intestinal mucosa to bacterial infection. Cohesion of this physical barrier is supported by renewal of epithelial stem cells, residing in invaginations called crypts, and by crypt cell migration onto protrusions called villi; dysregulation of such mechanisms may render the gut susceptible to chronic inflammation. The impact that excessive or misplaced epithelial cell death may have on villus cell migration is currently unknown. We integrated cell-tracking methods with computational models to determine how epithelial homeostasis is affected by acute and chronic TNFα-driven epithelial cell death. Parameter inference reveals that acute inflammatory cell death has a transient effect on epithelial cell dynamics, whereas cell death caused by chronic elevated TNFα causes a delay in the accumulation of labelled cells onto the villus compared to the control. Such a delay may be reproduced by using a cell-based model to simulate the dynamics of each cell in a crypt-villus geometry, showing that a prolonged increase in cell death slows the migration of cells from the crypt to the villus. This investigation highlights which injuries (acute or chronic) may be regenerated and which cause disruption of healthy epithelial homeostasis.

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.

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.

Lineage-dependent spatial and functional organization of the mammalian enteric nervous system.

The enteric nervous system (ENS) is essential for digestive function and gut homeostasis. Here we show that the amorphous neuroglia networks of the mouse ENS are composed of overlapping clonal units founded by postmigratory neural crest-derived progenitors. The spatial configuration of ENS clones depends on proliferation-driven local interactions of ENS progenitors with lineally unrelated neuroectodermal cells, the ordered colonization of the serosa-mucosa axis by clonal descendants, and gut expansion. Single-cell transcriptomics and mutagenesis analysis delineated dynamic molecular states of ENS progenitors and identified RET as a regulator of neurogenic commitment. Clonally related enteric neurons exhibit synchronous activity in response to network stimulation. Thus, lineage relationships underpin the organization of the peripheral nervous system.

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.

Evidence of metabolically active but non-culturable Listeria monocytogenes in long-term growth at 10 °C.

Cultures of Listeria monocytogenes at low temperatures (10 °C) in a broth model revealed long-term survival at about 0.04% cell density in relation to the log phase. In contrast, direct viable counts and PMA real-time PCR data suggested that 50% and 1% of the population retain membrane integrity, respectively. To elucidate the observed difference, the metabolic activity of the bacterial population was investigated by multiparametric flow cytometry, including the assessment of membrane integrity, reductase activity, as well as forward and side scatter properties. These analyses were complemented by 16S rRNA real-time PCR. The majority of the cells retained their membrane integrity and reductase activity until day 29. On day 42, 48.00 ± 4.00% (L. monocytogenes strain 3251) and 68.67 ± 3.74% (L. monocytogenes strain 535) of the cells had intact membranes, whereas 57.23 ± 1.85% (strain 3251) and 74.97 ± 3.01% (strain 535) exhibited high reductase activity. On day 42, mean 16S rRNA copy numbers of 3.98 ± 1.37 (membrane integrity) and 3.86 ± 1.32 (reductase activity) remained per intact or active cell. Our data suggest the transition of L. monocytogenes into a state of metabolic dormancy during long-term culture at low temperature.

Quantifying Lgr5-positive stem cell behaviour in the pyloric epithelium.

Using in-vivo lineage tracing data we quantified clonal expansion as well as proliferation and differentiation of the Lgr5-positive stem cell population in pyloric gastric glands. Fitting clone expansion models, we estimated that there are five effective Lgr5-positive cells able to give rise to monoclonal glands by replacing each other following a pattern of neutral drift dynamics. This analysis is instrumental to assess stem cell performance; however, stem cell proliferation is not quantified by clone expansion analysis. We identified a suitable mathematical model to quantify proliferation and differentiation of the Lgr5-positive population. As expected for populations in steady-state, the proliferation rate of the Lgr5-positive population was equal to its rate of differentiation. This rate was significantly faster than the rate at which effective cells are replaced, estimated by modelling clone expansion/contraction. This suggests that the majority of Lgr5-positive cell divisions serve to renew epithelial cells and only few result in the effective replacement of a neighbour to effect expansion to the entire gland. The application of the model under altered situations with uncoupled differentiation and proliferation was demonstrated. This methodology represents a valuable tool for quantifying stem cell performance in homeostasis and importantly for deciphering altered stem cell behaviour in disease.

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.