Daniele Muraro

Integration of hormonal signaling networks and mobile microRNAs is required for vascular patterning in Arabidopsis roots

Significance The vascular tissues form a continuous network providing the long-distance transport of water and nutrients in all higher plants (tracheophytes). To incorporate separate organs into this network, it is essential that the position of different vascular cell types is tightly regulated. Several factors required for root vascular patterning (including hormones and gene products) have previously been identified in the model plant Arabidopsis. We have now established a mathematical model formulizing the interaction between these factors, allowing us to identify a minimal regulatory network capable of maintaining a stable vascular pattern in Arabidopsis roots. We envisage that this model will help future researchers understand how similar regulatory units can be applied to create alternative patterns in other species. As multicellular organisms grow, positional information is continually needed to regulate the pattern in which cells are arranged. In the Arabidopsis root, most cell types are organized in a radially symmetric pattern; however, a symmetry-breaking event generates bisymmetric auxin and cytokinin signaling domains in the stele. Bidirectional cross-talk between the stele and the surrounding tissues involving a mobile transcription factor, SHORT ROOT (SHR), and mobile microRNA species also determines vascular pattern, but it is currently unclear how these signals integrate. We use a multicellular model to determine a minimal set of components necessary for maintaining a stable vascular pattern. Simulations perturbing the signaling network show that, in addition to the mutually inhibitory interaction between auxin and cytokinin, signaling through SHR, microRNA165/6, and PHABULOSA is required to maintain a stable bisymmetric pattern. We have verified this prediction by observing loss of bisymmetry in shr mutants. The model reveals the importance of several features of the network, namely the mutual degradation of microRNA165/6 and PHABULOSA and the existence of an additional negative regulator of cytokinin signaling. These components form a plausible mechanism capable of patterning vascular tissues in the absence of positional inputs provided by the transport of hormones from the shoot.

The influence of cytokinin-auxin cross-regulation on cell-fate determination in Arabidopsis thaliana root development

Root growth and development in Arabidopsis thaliana are sustained by a specialised zone termed the meristem, which contains a population of dividing and differentiating cells that are functionally analogous to a stem cell niche in animals. The hormones auxin and cytokinin control meristem size antagonistically. Local accumulation of auxin promotes cell division and the initiation of a lateral root primordium. By contrast, high cytokinin concentrations disrupt the regular pattern of divisions that characterises lateral root development, and promote differentiation. The way in which the hormones interact is controlled by a genetic regulatory network. In this paper, we propose a deterministic mathematical model to describe this network and present model simulations that reproduce the experimentally observed effects of cytokinin on the expression of auxin regulated genes. We show how auxin response genes and auxin efflux transporters may be affected by the presence of cytokinin. We also analyse and compare the responses of the hormones auxin and cytokinin to changes in their supply with the responses obtained by genetic mutations of SHY2, which encodes a protein that plays a key role in balancing cytokinin and auxin regulation of meristem size. We show that although shy2 mutations can qualitatively reproduce the effect of varying auxin and cytokinin supply on their response genes, some elements of the network respond differently to changes in hormonal supply and to genetic mutations, implying a different, general response of the network. We conclude that an analysis based on the ratio between these two hormones may be misleading and that a mathematical model can serve as a useful tool for stimulate further experimental work by predicting the response of the network to changes in hormone levels and to other genetic mutations.

The role of auxin and cytokinin signalling in specifying the root architecture of Arabidopsis thaliana

Auxin and cytokinin are key hormonal signals that control the cellular architecture of the primary root and the initiation of new lateral root organs in the plant Arabidopsis thaliana. Both developmental processes are regulated by cross-talk between these hormones and their signalling pathways. In this paper, sub-cellular and multi-cellular mathematical models are developed to investigate how interactions between auxin and cytokinin influence the size and location of regions of division and differentiation within the primary root, and describe how their cross-regulation may cause periodic branching of lateral roots. We show how their joint activity may influence tissue-specific oscillations in gene expression, as shown in Moreno-Risueno et al. (2010) and commented upon in Traas and Vernoux (2010), and we propose mechanisms that may generate synchronisation of such periodic behaviours inside a cell and with its neighbours. Using a multi-cellular model, we also analyse the roles of cytokinin and auxin in specifying the three main regions of the primary root (elongation, transition and division zones), our simulation results being in good agreement with independent experimental observations. We then use our model to generate testable predictions concerning the effect of varying the concentrations of the auxin efflux transporters on the sizes of the different root regions. In particular, we predict that over-expression of the transporters will generate a longer root with a longer elongation zone and a smaller division zone than that of a wild type root. This root will contain fewer cells than its wild type counterpart. We conclude that our model can provide a useful tool for investigating the response of cell division and elongation to perturbations in hormonal signalling.

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. FASEB J. 31, 636–649 (2017). www.fasebj.org

Validation of a Morphogenesis Model of Drosophila Early Development by a Multi-objective Evolutionary Optimization Algorithm

We apply evolutionary computation to calibrate the parameters of a morphogenesis model of Drosophila early development. The model aims to describe the establishment of the steady gradients of Bicoid and Caudal proteins along the antero-posterior axis of the embryo of Drosophila . The model equations consist of a system of non-linear parabolic partial differential equations with initial and zero flux boundary conditions. We compare the results of single- and multi-objective variants of the CMA-ES algorithm for the model the calibration with the experimental data. Whereas the multi-objective algorithm computes a full approximation of the Pareto front, repeated runs of the single-objective algorithm give solutions that dominate (in the Pareto sense) the results of the multi-objective approach. We retain as best solutions those found by the latter technique. From the biological point of view, all such solutions are all equally acceptable, and for our test cases, the relative error between the experimental data and validated model solutions on the Pareto front are in the range 3% *** 6%. This technique is general and can be used as a generic tool for parameter calibration problems.

Root development: cytokinin transport matters, too!

Unlike the plant hormone auxin, the mechanism and function of cytokinin transport is poorly characterised. Two new studies now demonstrate that cytokinins transported from shoot to roots via the phloem are critical for creating mutually exclusive auxin and cytokinin signalling domains that control root vascular patterning.

A Software Tool to Model Genetic Regulatory Networks. Applications to the Modeling of Threshold Phenomena and of Spatial Patterning in Drosophila

We present a general methodology in order to build mathematical models of genetic regulatory networks. This approach is based on the mass action law and on the Jacob and Monod operon model. The mathematical models are built symbolically by the Mathematica software package GeneticNetworks. This package accepts as input the interaction graphs of the transcriptional activators and repressors of a biological process and, as output, gives the mathematical model in the form of a system of ordinary differential equations. All the relevant biological parameters are chosen automatically by the software. Within this framework, we show that concentration dependent threshold effects in biology emerge from the catalytic properties of genes and its associated conservation laws. We apply this methodology to the segment patterning in Drosophila early development and we calibrate the genetic transcriptional network responsible for the patterning of the gap gene proteins Hunchback and Knirps, along the antero-posterior axis of the Drosophila embryo. In this approach, the zygotically produced proteins Hunchback and Knirps do not diffuse along the antero-posterior axis of the embryo of Drosophila, developing a spatial pattern due to concentration dependent thresholds. This shows that patterning at the gap genes stage can be explained by the concentration gradients along the embryo of the transcriptional regulators.

An integrative analysis of gene expression and molecular interaction data to identify dys-regulated sub-networks in inflammatory bowel disease

BackgroundInflammatory bowel disease (IBD) consists of two main disease-subtypes, Crohn’s disease (CD) and ulcerative colitis (UC); these subtypes share overlapping genetic and clinical features. Genome-wide microarray data enable unbiased documentation of alterations in gene expression that may be disease-specific. As genetic diseases are believed to be caused by genetic alterations affecting the function of signalling pathways, module-centric optimisation algorithms, whose aim is to identify sub-networks that are dys-regulated in disease, are emerging as promising approaches.ResultsIn order to account for the topological structure of molecular interaction networks, we developed an optimisation algorithm that integrates databases of known molecular interactions with gene expression data; such integration enables identification of differentially regulated network modules. We verified the performance of our algorithm by testing it on simulated networks; we then applied the same method to study experimental data derived from microarray analysis of CD and UC biopsies and human interactome databases. This analysis allowed the extraction of dys-regulated subnetworks under different experimental conditions (inflamed and uninflamed tissues in CD and UC). Optimisation was performed to highlight differentially expressed network modules that may be common or specific to the disease subtype.ConclusionsWe show that the selected subnetworks include genes and pathways of known relevance for IBD; in particular, the solutions found highlight cross-talk among enriched pathways, mainly the JAK/STAT signalling pathway and the EGF receptor signalling pathway. In addition, integration of gene expression with molecular interaction data highlights nodes that, although not being differentially expressed, interact with differentially expressed nodes and are part of pathways that are relevant to IBD. The method proposed here may help identifying dys-regulated sub-networks that are common in different diseases and sub-networks whose dys-regulation is specific to a particular disease.

A parallel multi-objective optimization algorithm for the calibration of mathematical models

Abstract Based on evolutionary computation techniques, we present a parallel, globally convergent, multi-objective optimization algorithm which extends the Covariance Matrix Adaptation Evolutionary Strategy (CMA-ES). This approach enables identifying multiple global optima and multiple discontinuous Pareto set solutions of the optimization problem in a compact search space. After evaluating the algorithm with test functions, we apply our method to the identification of the parameters of a reaction–diffusion model of a genetic regulatory mechanism during Drosophila early development, our simulations being in agreement with the experimental data. Comparisons with a multi-objective version of the CMA-ES (MO-CMA-ES) on this dataset show that our algorithm outperforms largely the speed of convergence of MO-CMA-ES. We have identified an infinite number of accurate solutions of the model equations, associated with the Pareto set of the optimization problem. This non-unicity property of a biological developmental process explains phenotypic plasticity and resilience in biological systems.

Snapin promotes HIV-1 transmission from dendritic cells by dampening TLR8 signaling.

HIV‐1 traffics through dendritic cells (DCs) en route to establishing a productive infection in T lymphocytes but fails to induce an innate immune response. Within DC endosomes, HIV‐1 somehow evades detection by the pattern‐recognition receptor (PRR) Toll‐like receptor 8 (TLR8). Using a phosphoproteomic approach, we identified a robust and diverse signaling cascade triggered by HIV‐1 upon entry into human DCs. A secondary siRNA screen of the identified signaling factors revealed several new mediators of HIV‐1 trans‐infection of CD4+ T cells in DCs, including the dynein motor protein Snapin. Inhibition of Snapin enhanced localization of HIV‐1 with TLR8+ early endosomes, triggered a pro‐inflammatory response, and inhibited trans‐infection of CD4+ T cells. Snapin inhibited TLR8 signaling in the absence of HIV‐1 and is a general regulator of endosomal maturation. Thus, we identify a new mechanism of innate immune sensing by TLR8 in DCs, which is exploited by HIV‐1 to promote transmission.