1
|
Andreas E, Cummins B, Gedeon T. Quantifying robustness of the gap gene network. J Theor Biol 2024; 580:111720. [PMID: 38211890 DOI: 10.1016/j.jtbi.2024.111720] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/02/2023] [Revised: 12/28/2023] [Accepted: 12/30/2023] [Indexed: 01/13/2024]
Abstract
Early development of Drosophila melanogaster (fruit fly) facilitated by the gap gene network has been shown to be incredibly robust, and the same patterns emerge even when the process is seriously disrupted. We investigate this robustness using a previously developed computational framework called DSGRN (Dynamic Signatures Generated by Regulatory Networks). Our mathematical innovations include the conceptual extension of this established modeling technique to enable modeling of spatially monotone environmental effects, as well as the development of a collection of graph theoretic robustness scores for network models. This allows us to rank order the robustness of network models of cellular systems where each cell contains the same genetic network topology but operates under a parameter regime that changes continuously from cell to cell. We demonstrate the power of this method by comparing the robustness of two previously introduced network models of gap gene expression along the anterior-posterior axis of the fruit fly embryo, both to each other and to a random sample of networks with same number of nodes and edges. We observe that there is a substantial difference in robustness scores between the two models. Our biological insight is that random network topologies are in general capable of reproducing complex patterns of expression, but that using measures of robustness to rank order networks permits a large reduction in hypothesis space for highly conserved systems such as developmental networks.
Collapse
Affiliation(s)
- Elizabeth Andreas
- Department of Mathematical Sciences, Montana State University, Bozeman, 59718, MT, USA.
| | - Breschine Cummins
- Department of Mathematical Sciences, Montana State University, Bozeman, 59718, MT, USA
| | - Tomáš Gedeon
- Department of Mathematical Sciences, Montana State University, Bozeman, 59718, MT, USA
| |
Collapse
|
2
|
Santos-Moreno J, Tasiudi E, Kusumawardhani H, Stelling J, Schaerli Y. Robustness and innovation in synthetic genotype networks. Nat Commun 2023; 14:2454. [PMID: 37117168 PMCID: PMC10147661 DOI: 10.1038/s41467-023-38033-3] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/11/2022] [Accepted: 04/13/2023] [Indexed: 04/30/2023] Open
Abstract
Genotype networks are sets of genotypes connected by small mutational changes that share the same phenotype. They facilitate evolutionary innovation by enabling the exploration of different neighborhoods in genotype space. Genotype networks, first suggested by theoretical models, have been empirically confirmed for proteins and RNAs. Comparative studies also support their existence for gene regulatory networks (GRNs), but direct experimental evidence is lacking. Here, we report the construction of three interconnected genotype networks of synthetic GRNs producing three distinct phenotypes in Escherichia coli. Our synthetic GRNs contain three nodes regulating each other by CRISPR interference and governing the expression of fluorescent reporters. The genotype networks, composed of over twenty different synthetic GRNs, provide robustness in face of mutations while enabling transitions to innovative phenotypes. Through realistic mathematical modeling, we quantify robustness and evolvability for the complete genotype-phenotype map and link these features mechanistically to GRN motifs. Our work thereby exemplifies how GRN evolution along genotype networks might be driving evolutionary innovation.
Collapse
Affiliation(s)
- Javier Santos-Moreno
- Department of Fundamental Microbiology, University of Lausanne, Biophore Building, 1015, Lausanne, Switzerland
- Department of Medicine and Life Sciences, Pompeu Fabra University, 00803, Barcelona, Spain
| | - Eve Tasiudi
- Department of Biosystems Science and Engineering, ETH Zurich and SIB Swiss Institute of Bioinformatics, Basel, Switzerland
| | - Hadiastri Kusumawardhani
- Department of Fundamental Microbiology, University of Lausanne, Biophore Building, 1015, Lausanne, Switzerland
| | - Joerg Stelling
- Department of Biosystems Science and Engineering, ETH Zurich and SIB Swiss Institute of Bioinformatics, Basel, Switzerland.
| | - Yolanda Schaerli
- Department of Fundamental Microbiology, University of Lausanne, Biophore Building, 1015, Lausanne, Switzerland.
| |
Collapse
|
3
|
Hawkes WL, Sivell O, Wotton KR. The genome sequence of the Marmalade Hoverfly, Episyrphus balteatus (De Geer, 1776). Wellcome Open Res 2023. [DOI: 10.12688/wellcomeopenres.19073.1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/26/2023] Open
Abstract
We present a genome assembly from an individual female Episyrphus balteatus (the Marmalade Hoverfly; Arthropoda; Insecta; Diptera; Syrphidae). The genome sequence is 535 megabases in span. Most of the assembly is scaffolded into five chromosomal pseudomolecules, including the assembled X sex chromosome. The mitochondrial genome has also been assembled and is 16.8 kilobases in length.
Collapse
|
4
|
Perkins ML, Gandara L, Crocker J. A synthetic synthesis to explore animal evolution and development. Philos Trans R Soc Lond B Biol Sci 2022; 377:20200517. [PMID: 35634925 PMCID: PMC9149795 DOI: 10.1098/rstb.2020.0517] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
Abstract
Identifying the general principles by which genotypes are converted into phenotypes remains a challenge in the post-genomic era. We still lack a predictive understanding of how genes shape interactions among cells and tissues in response to signalling and environmental cues, and hence how regulatory networks generate the phenotypic variation required for adaptive evolution. Here, we discuss how techniques borrowed from synthetic biology may facilitate a systematic exploration of evolvability across biological scales. Synthetic approaches permit controlled manipulation of both endogenous and fully engineered systems, providing a flexible platform for investigating causal mechanisms in vivo. Combining synthetic approaches with multi-level phenotyping (phenomics) will supply a detailed, quantitative characterization of how internal and external stimuli shape the morphology and behaviour of living organisms. We advocate integrating high-throughput experimental data with mathematical and computational techniques from a variety of disciplines in order to pursue a comprehensive theory of evolution. This article is part of the theme issue ‘Genetic basis of adaptation and speciation: from loci to causative mutations’.
Collapse
Affiliation(s)
- Mindy Liu Perkins
- Developmental Biology Unit, European Molecular Biology Laboratory, 69117 Heidelberg, Germany
| | - Lautaro Gandara
- Developmental Biology Unit, European Molecular Biology Laboratory, 69117 Heidelberg, Germany
| | - Justin Crocker
- Developmental Biology Unit, European Molecular Biology Laboratory, 69117 Heidelberg, Germany
| |
Collapse
|
5
|
Nogueira Alves A, Oliveira MM, Koyama T, Shingleton A, Mirth CK. Ecdysone coordinates plastic growth with robust pattern in the developing wing. eLife 2022; 11:72666. [PMID: 35261337 PMCID: PMC8947767 DOI: 10.7554/elife.72666] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2021] [Accepted: 03/07/2022] [Indexed: 11/25/2022] Open
Abstract
Animals develop in unpredictable, variable environments. In response to environmental change, some aspects of development adjust to generate plastic phenotypes. Other aspects of development, however, are buffered against environmental change to produce robust phenotypes. How organ development is coordinated to accommodate both plastic and robust developmental responses is poorly understood. Here, we demonstrate that the steroid hormone ecdysone coordinates both plasticity of organ size and robustness of organ pattern in the developing wings of the fruit fly Drosophila melanogaster. Using fed and starved larvae that lack prothoracic glands, which synthesize ecdysone, we show that nutrition regulates growth both via ecdysone and via an ecdysone-independent mechanism, while nutrition regulates patterning only via ecdysone. We then demonstrate that growth shows a graded response to ecdysone concentration, while patterning shows a threshold response. Collectively, these data support a model where nutritionally regulated ecdysone fluctuations confer plasticity by regulating disc growth in response to basal ecdysone levels and confer robustness by initiating patterning only once ecdysone peaks exceed a threshold concentration. This could represent a generalizable mechanism through which hormones coordinate plastic growth with robust patterning in the face of environmental change.
Collapse
Affiliation(s)
| | | | | | - Alexander Shingleton
- Department of Biological Sciences, University of Illinois at Chicago, Chicago, United States
| | - Christen K Mirth
- School of Biological Sciences, Monash University, Melbourne, Australia
| |
Collapse
|
6
|
Abstract
Even if a species' phenotype does not change over evolutionary time, the underlying mechanism may change, as distinct molecular pathways can realize identical phenotypes. Here we use linear system theory to explore the consequences of this idea, describing how a gene network underlying a conserved phenotype evolves, as the genetic drift of small changes to these molecular pathways causes a population to explore the set of mechanisms with identical phenotypes. To do this, we model an organism's internal state as a linear system of differential equations for which the environment provides input and the phenotype is the output, in which context there exists an exact characterization of the set of all mechanisms that give the same input-output relationship. This characterization implies that selectively neutral directions in genotype space should be common and that the evolutionary exploration of these distinct but equivalent mechanisms can lead to the reproductive incompatibility of independently evolving populations. This evolutionary exploration, or system drift, is expected to proceed at a rate proportional to the amount of intrapopulation genetic variation divided by the effective population size ( Ne$N_e$ ). At biologically reasonable parameter values this could lead to substantial interpopulation incompatibility, and thus speciation, on a time scale of Ne$N_e$ generations. This model also naturally predicts Haldane's rule, thus providing a concrete explanation of why heterogametic hybrids tend to be disrupted more often than homogametes during the early stages of speciation.
Collapse
Affiliation(s)
- Joshua S. Schiffman
- New York Genome CenterNew YorkNew York 10013,Weill Cornell MedicineNew YorkNew York 10065,Department of Molecular and Computational BiologyUniversity of Southern CaliforniaLos AngelesCalifornia 90089
| | - Peter L. Ralph
- Department of Molecular and Computational BiologyUniversity of Southern CaliforniaLos AngelesCalifornia 90089,Department of Mathematics, Institute of Ecology and EvolutionUniversity of OregonEugeneOregon 97403,Department of Biology, Institute of Ecology and EvolutionUniversity of OregonEugeneOregon 97403
| |
Collapse
|
7
|
Chipman AD. The evolution of the gene regulatory networks patterning the Drosophila Blastoderm. Curr Top Dev Biol 2021; 139:297-324. [PMID: 32450964 DOI: 10.1016/bs.ctdb.2020.02.004] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
The Drosophila blastoderm gene regulatory network is one of the best studied networks in biology. It is composed of a series of tiered sub-networks that act sequentially to generate a primary segmental pattern. Many of these sub-networks have been studied in other arthropods, allowing us to reconstruct how each of them evolved over the transition from the arthropod ancestor to the situation seen in Drosophila today. I trace the evolution of each of these networks, showing how some of them have been modified significantly in Drosophila relative to the ancestral state while others are largely conserved across evolutionary timescales. I compare the putative ancestral arthropod segmentation network with that found in Drosophila and discuss how and why it has been modified throughout evolution, and to what extent this modification is unusual.
Collapse
Affiliation(s)
- Ariel D Chipman
- The Department of Ecology, Evolution & Behavior, The Hebrew University of Jerusalem, Edmond J. Safra Campus, Givat Ram, Jerusalem, Israel.
| |
Collapse
|
8
|
DiFrisco J, Jaeger J. Homology of process: developmental dynamics in comparative biology. Interface Focus 2021; 11:20210007. [PMID: 34055306 PMCID: PMC8086918 DOI: 10.1098/rsfs.2021.0007] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 02/22/2021] [Indexed: 12/14/2022] Open
Abstract
Comparative biology builds up systematic knowledge of the diversity of life, across evolutionary lineages and levels of organization, starting with evidence from a sparse sample of model organisms. In developmental biology, a key obstacle to the growth of comparative approaches is that the concept of homology is not very well defined for levels of organization that are intermediate between individual genes and morphological characters. In this paper, we investigate what it means for ontogenetic processes to be homologous, focusing specifically on the examples of insect segmentation and vertebrate somitogenesis. These processes can be homologous without homology of the underlying genes or gene networks, since the latter can diverge over evolutionary time, while the dynamics of the process remain the same. Ontogenetic processes like these therefore constitute a dissociable level and distinctive unit of comparison requiring their own specific criteria of homology. In addition, such processes are typically complex and nonlinear, such that their rigorous description and comparison requires not only observation and experimentation, but also dynamical modelling. We propose six criteria of process homology, combining recognized indicators (sameness of parts, morphological outcome and topological position) with novel ones derived from dynamical systems modelling (sameness of dynamical properties, dynamical complexity and evidence for transitional forms). We show how these criteria apply to animal segmentation and other ontogenetic processes. We conclude by situating our proposed dynamical framework for homology of process in relation to similar research programmes, such as process structuralism and developmental approaches to morphological homology.
Collapse
Affiliation(s)
- James DiFrisco
- Institute of Philosophy, KU Leuven, 3000 Leuven, Belgium
| | - Johannes Jaeger
- Complexity Science Hub (CSH) Vienna, Josefstädter Strasse 39, 1080 Vienna, Austria
| |
Collapse
|
9
|
Abstract
Arthropod segmentation and vertebrate somitogenesis are leading fields in the experimental and theoretical interrogation of developmental patterning. However, despite the sophistication of current research, basic conceptual issues remain unresolved. These include: (i) the mechanistic origins of spatial organization within the segment addition zone (SAZ); (ii) the mechanistic origins of segment polarization; (iii) the mechanistic origins of axial variation; and (iv) the evolutionary origins of simultaneous patterning. Here, I explore these problems using coarse-grained models of cross-regulating dynamical processes. In the morphogenetic framework of a row of cells undergoing axial elongation, I simulate interactions between an 'oscillator', a 'switch' and up to three 'timers', successfully reproducing essential patterning behaviours of segmenting systems. By comparing the output of these largely cell-autonomous models to variants that incorporate positional information, I find that scaling relationships, wave patterns and patterning dynamics all depend on whether the SAZ is regulated by temporal or spatial information. I also identify three mechanisms for polarizing oscillator output, all of which functionally implicate the oscillator frequency profile. Finally, I demonstrate significant dynamical and regulatory continuity between sequential and simultaneous modes of segmentation. I discuss these results in the context of the experimental literature.
Collapse
Affiliation(s)
- Erik Clark
- Department of Systems Biology, Harvard Medical School, 210 Longwood Ave, Boston, MA 02115, USA
- Trinity College Cambridge, University of Cambridge, Trinity Street, Cambridge CB2 1TQ, UK
| |
Collapse
|
10
|
Perkins ML. Implications of diffusion and time-varying morphogen gradients for the dynamic positioning and precision of bistable gene expression boundaries. PLoS Comput Biol 2021; 17:e1008589. [PMID: 34061823 PMCID: PMC8195430 DOI: 10.1371/journal.pcbi.1008589] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/09/2020] [Revised: 06/11/2021] [Accepted: 05/11/2021] [Indexed: 11/19/2022] Open
Abstract
The earliest models for how morphogen gradients guide embryonic patterning failed to account for experimental observations of temporal refinement in gene expression domains. Following theoretical and experimental work in this area, dynamic positional information has emerged as a conceptual framework to discuss how cells process spatiotemporal inputs into downstream patterns. Here, we show that diffusion determines the mathematical means by which bistable gene expression boundaries shift over time, and therefore how cells interpret positional information conferred from morphogen concentration. First, we introduce a metric for assessing reproducibility in boundary placement or precision in systems where gene products do not diffuse, but where morphogen concentrations are permitted to change in time. We show that the dynamics of the gradient affect the sensitivity of the final pattern to variation in initial conditions, with slower gradients reducing the sensitivity. Second, we allow gene products to diffuse and consider gene expression boundaries as propagating wavefronts with velocity modulated by local morphogen concentration. We harness this perspective to approximate a PDE model as an ODE that captures the position of the boundary in time, and demonstrate the approach with a preexisting model for Hunchback patterning in fruit fly embryos. We then propose a design that employs antiparallel morphogen gradients to achieve accurate boundary placement that is robust to scaling. Throughout our work we draw attention to tradeoffs among initial conditions, boundary positioning, and the relative timescales of network and gradient evolution. We conclude by suggesting that mathematical theory should serve to clarify not just our quantitative, but also our intuitive understanding of patterning processes.
Collapse
Affiliation(s)
- Melinda Liu Perkins
- Developmental Biology Unit, European Molecular Biology Laboratory, Heidelberg, Germany
- * E-mail:
| |
Collapse
|
11
|
Manrubia S, Cuesta JA, Aguirre J, Ahnert SE, Altenberg L, Cano AV, Catalán P, Diaz-Uriarte R, Elena SF, García-Martín JA, Hogeweg P, Khatri BS, Krug J, Louis AA, Martin NS, Payne JL, Tarnowski MJ, Weiß M. From genotypes to organisms: State-of-the-art and perspectives of a cornerstone in evolutionary dynamics. Phys Life Rev 2021; 38:55-106. [PMID: 34088608 DOI: 10.1016/j.plrev.2021.03.004] [Citation(s) in RCA: 36] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/24/2020] [Accepted: 03/01/2021] [Indexed: 12/21/2022]
Abstract
Understanding how genotypes map onto phenotypes, fitness, and eventually organisms is arguably the next major missing piece in a fully predictive theory of evolution. We refer to this generally as the problem of the genotype-phenotype map. Though we are still far from achieving a complete picture of these relationships, our current understanding of simpler questions, such as the structure induced in the space of genotypes by sequences mapped to molecular structures, has revealed important facts that deeply affect the dynamical description of evolutionary processes. Empirical evidence supporting the fundamental relevance of features such as phenotypic bias is mounting as well, while the synthesis of conceptual and experimental progress leads to questioning current assumptions on the nature of evolutionary dynamics-cancer progression models or synthetic biology approaches being notable examples. This work delves with a critical and constructive attitude into our current knowledge of how genotypes map onto molecular phenotypes and organismal functions, and discusses theoretical and empirical avenues to broaden and improve this comprehension. As a final goal, this community should aim at deriving an updated picture of evolutionary processes soundly relying on the structural properties of genotype spaces, as revealed by modern techniques of molecular and functional analysis.
Collapse
Affiliation(s)
- Susanna Manrubia
- Department of Systems Biology, Centro Nacional de Biotecnología (CSIC), Madrid, Spain; Grupo Interdisciplinar de Sistemas Complejos (GISC), Madrid, Spain.
| | - José A Cuesta
- Grupo Interdisciplinar de Sistemas Complejos (GISC), Madrid, Spain; Departamento de Matemáticas, Universidad Carlos III de Madrid, Leganés, Spain; Instituto de Biocomputación y Física de Sistemas Complejos (BiFi), Universidad de Zaragoza, Spain; UC3M-Santander Big Data Institute (IBiDat), Getafe, Madrid, Spain
| | - Jacobo Aguirre
- Grupo Interdisciplinar de Sistemas Complejos (GISC), Madrid, Spain; Centro de Astrobiología, CSIC-INTA, ctra. de Ajalvir km 4, 28850 Torrejón de Ardoz, Madrid, Spain
| | - Sebastian E Ahnert
- Department of Chemical Engineering and Biotechnology, University of Cambridge, Philippa Fawcett Drive, Cambridge CB3 0AS, UK; The Alan Turing Institute, British Library, 96 Euston Road, London NW1 2DB, UK
| | | | - Alejandro V Cano
- Institute of Integrative Biology, ETH Zurich, Zurich, Switzerland; Swiss Institute of Bioinformatics, Lausanne, Switzerland
| | - Pablo Catalán
- Grupo Interdisciplinar de Sistemas Complejos (GISC), Madrid, Spain; Departamento de Matemáticas, Universidad Carlos III de Madrid, Leganés, Spain
| | - Ramon Diaz-Uriarte
- Department of Biochemistry, Universidad Autónoma de Madrid, Madrid, Spain; Instituto de Investigaciones Biomédicas "Alberto Sols" (UAM-CSIC), Madrid, Spain
| | - Santiago F Elena
- Instituto de Biología Integrativa de Sistemas, I(2)SysBio (CSIC-UV), València, Spain; The Santa Fe Institute, Santa Fe, NM, USA
| | | | - Paulien Hogeweg
- Theoretical Biology and Bioinformatics Group, Utrecht University, the Netherlands
| | - Bhavin S Khatri
- The Francis Crick Institute, London, UK; Department of Life Sciences, Imperial College London, London, UK
| | - Joachim Krug
- Institute for Biological Physics, University of Cologne, Köln, Germany
| | - Ard A Louis
- Rudolf Peierls Centre for Theoretical Physics, University of Oxford, Oxford, UK
| | - Nora S Martin
- Theory of Condensed Matter Group, Cavendish Laboratory, University of Cambridge, Cambridge, UK; Sainsbury Laboratory, University of Cambridge, Cambridge, UK
| | - Joshua L Payne
- Institute of Integrative Biology, ETH Zurich, Zurich, Switzerland; Swiss Institute of Bioinformatics, Lausanne, Switzerland
| | | | - Marcel Weiß
- Theory of Condensed Matter Group, Cavendish Laboratory, University of Cambridge, Cambridge, UK; Sainsbury Laboratory, University of Cambridge, Cambridge, UK
| |
Collapse
|
12
|
Abstract
Diverse cellular phenotypes are determined by groups of transcription factors (TFs) and other regulators that influence each others' gene expression, forming transcriptional gene regulatory networks (GRNs). In many biological contexts, especially in development and associated diseases, the expression of the genes in GRNs is not static but evolves in time. Modeling the dynamics of GRN state is an important approach for understanding diverse cellular phenomena such as cell-fate specification, pluripotency and cell-fate reprogramming, oncogenesis, and tissue regeneration. In this protocol, we describe how to model GRNs using a data-driven dynamic modeling methodology, gene circuits. Gene circuits do not require knowledge of the GRN topology and connectivity but instead learn them from training data, making them very general and applicable to diverse biological contexts. We utilize the MATLAB-based gene circuit modeling software Fast Inference of Gene Regulation (FIGR) for training the model on quantitative gene expression data and simulating the GRN. We describe all the steps in the modeling life cycle, from formulating the model, training the model using FIGR, simulating the GRN, to analyzing and interpreting the model output. This protocol highlights these steps with the example of a dynamical model of the gap gene GRN involved in Drosophila segmentation and includes example MATLAB statements for each step.
Collapse
Affiliation(s)
- Joanna E Handzlik
- Department of Biology, University of North Dakota, Grand Forks, ND, 58202, USA
| | - Yen Lee Loh
- Department of Physics and Astrophysics, University of North Dakota, Grand Forks, ND, 58202, USA
| | - Manu
- Department of Biology, University of North Dakota, Grand Forks, ND, 58202, USA.
| |
Collapse
|
13
|
Oates AC. Waiting on the Fringe: cell autonomy and signaling delays in segmentation clocks. Curr Opin Genet Dev 2020; 63:61-70. [PMID: 32505051 DOI: 10.1016/j.gde.2020.04.008] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/16/2020] [Revised: 04/19/2020] [Accepted: 04/23/2020] [Indexed: 12/16/2022]
Abstract
The rhythmic and sequential segmentation of the vertebrate body axis into somites during embryogenesis is governed by a multicellular, oscillatory patterning system called the segmentation clock. Despite many overt similarities between vertebrates, differences in genetic and dynamic regulation have been reported, raising intriguing questions about the evolution and conservation of this fundamental patterning process. Recent studies have brought insights into two important and related issues: (1) whether individual cells of segmentation clocks are autonomous oscillators or require cell-cell communication for their rhythm; and (2) the role of delays in the cell-cell communication that synchronizes the population of genetic oscillators. Although molecular details differ between species, conservation may exist at the level of the dynamics, hinting at rules for evolutionary trajectories in the system.
Collapse
Affiliation(s)
- Andrew C Oates
- Institute of Bioengineering, School of Life Sciences and School of Engineering, Ecole Polytechnique Fédéral de Lausanne (EPFL), CH-1015, Switzerland.
| |
Collapse
|
14
|
Huang A, Rupprecht JF, Saunders TE. Embryonic geometry underlies phenotypic variation in decanalized conditions. eLife 2020; 9:e47380. [PMID: 32048988 PMCID: PMC7032927 DOI: 10.7554/elife.47380] [Citation(s) in RCA: 21] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/03/2019] [Accepted: 02/11/2020] [Indexed: 11/13/2022] Open
Abstract
During development, many mutations cause increased variation in phenotypic outcomes, a phenomenon termed decanalization. Phenotypic discordance is often observed in the absence of genetic and environmental variations, but the mechanisms underlying such inter-individual phenotypic discordance remain elusive. Here, using the anterior-posterior (AP) patterning of the Drosophila embryo, we identified embryonic geometry as a key factor predetermining patterning outcomes under decanalizing mutations. With the wild-type AP patterning network, we found that AP patterning is robust to variations in embryonic geometry; segmentation gene expression remains reproducible even when the embryo aspect ratio is artificially reduced by more than twofold. In contrast, embryonic geometry is highly predictive of individual patterning defects under decanalized conditions of either increased bicoid (bcd) dosage or bcd knockout. We showed that the phenotypic discordance can be traced back to variations in the gap gene expression, which is rendered sensitive to the geometry of the embryo under mutations.
Collapse
Affiliation(s)
- Anqi Huang
- Mechanobiology Institute, National University of SingaporeSingaporeSingapore
| | - Jean-François Rupprecht
- Mechanobiology Institute, National University of SingaporeSingaporeSingapore
- CNRS and Turing Center for Living Systems, Centre de Physique Théorique, Aix-Marseille UniversitéMarseilleFrance
| | - Timothy E Saunders
- Mechanobiology Institute, National University of SingaporeSingaporeSingapore
- Department of Biological Sciences, National University of SingaporeSingaporeSingapore
- Institute of Molecular and Cell Biology, Proteos, A*StarSingaporeSingapore
| |
Collapse
|
15
|
An Atlas of Transcription Factors Expressed in Male Pupal Terminalia of Drosophila melanogaster. G3-GENES GENOMES GENETICS 2019; 9:3961-3972. [PMID: 31619460 PMCID: PMC6893207 DOI: 10.1534/g3.119.400788] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
During development, transcription factors and signaling molecules govern gene regulatory networks to direct the formation of unique morphologies. As changes in gene regulatory networks are often implicated in morphological evolution, mapping transcription factor landscapes is important, especially in tissues that undergo rapid evolutionary change. The terminalia (genital and anal structures) of Drosophila melanogaster and its close relatives exhibit dramatic changes in morphology between species. While previous studies have identified network components important for patterning the larval genital disc, the networks governing adult structures during pupal development have remained uncharted. Here, we performed RNA-seq in whole Drosophila melanogaster male terminalia followed by in situ hybridization for 100 highly expressed transcription factors during pupal development. We find that the male terminalia are highly patterned during pupal stages and that specific transcription factors mark separate structures and substructures. Our results are housed online in a searchable database (https://flyterminalia.pitt.edu/) as a resource for the community. This work lays a foundation for future investigations into the gene regulatory networks governing the development and evolution of Drosophila terminalia.
Collapse
|
16
|
Classification-Based Inference of Dynamical Models of Gene Regulatory Networks. G3-GENES GENOMES GENETICS 2019; 9:4183-4195. [PMID: 31624138 PMCID: PMC6893186 DOI: 10.1534/g3.119.400603] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 12/05/2022]
Abstract
Cell-fate decisions during development are controlled by densely interconnected gene regulatory networks (GRNs) consisting of many genes. Inferring and predictively modeling these GRNs is crucial for understanding development and other physiological processes. Gene circuits, coupled differential equations that represent gene product synthesis with a switch-like function, provide a biologically realistic framework for modeling the time evolution of gene expression. However, their use has been limited to smaller networks due to the computational expense of inferring model parameters from gene expression data using global non-linear optimization. Here we show that the switch-like nature of gene regulation can be exploited to break the gene circuit inference problem into two simpler optimization problems that are amenable to computationally efficient supervised learning techniques. We present FIGR (Fast Inference of Gene Regulation), a novel classification-based inference approach to determining gene circuit parameters. We demonstrate FIGR’s effectiveness on synthetic data generated from random gene circuits of up to 50 genes as well as experimental data from the gap gene system of Drosophila melanogaster, a benchmark for inferring dynamical GRN models. FIGR is faster than global non-linear optimization by a factor of 600 and its computational complexity scales much better with GRN size. On a practical level, FIGR can accurately infer the biologically realistic gap gene network in under a minute on desktop-class hardware instead of requiring hours of parallel computing. We anticipate that FIGR would enable the inference of much larger biologically realistic GRNs than was possible before.
Collapse
|
17
|
Coronado-Zamora M, Salvador-Martínez I, Castellano D, Barbadilla A, Salazar-Ciudad I. Adaptation and Conservation throughout the Drosophila melanogaster Life-Cycle. Genome Biol Evol 2019; 11:1463-1482. [PMID: 31028390 PMCID: PMC6535812 DOI: 10.1093/gbe/evz086] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 04/16/2019] [Indexed: 01/09/2023] Open
Abstract
Previous studies of the evolution of genes expressed at different life-cycle stages of Drosophila melanogaster have not been able to disentangle adaptive from nonadaptive substitutions when using nonsynonymous sites. Here, we overcome this limitation by combining whole-genome polymorphism data from D. melanogaster and divergence data between D. melanogaster and Drosophila yakuba. For the set of genes expressed at different life-cycle stages of D. melanogaster, as reported in modENCODE, we estimate the ratio of substitutions relative to polymorphism between nonsynonymous and synonymous sites (α) and then α is discomposed into the ratio of adaptive (ωa) and nonadaptive (ωna) substitutions to synonymous substitutions. We find that the genes expressed in mid- and late-embryonic development are the most conserved, whereas those expressed in early development and postembryonic stages are the least conserved. Importantly, we found that low conservation in early development is due to high rates of nonadaptive substitutions (high ωna), whereas in postembryonic stages it is due, instead, to high rates of adaptive substitutions (high ωa). By using estimates of different genomic features (codon bias, average intron length, exon number, recombination rate, among others), we also find that genes expressed in mid- and late-embryonic development show the most complex architecture: they are larger, have more exons, more transcripts, and longer introns. In addition, these genes are broadly expressed among all stages. We suggest that all these genomic features are related to the conservation of mid- and late-embryonic development. Globally, our study supports the hourglass pattern of conservation and adaptation over the life-cycle.
Collapse
Affiliation(s)
- Marta Coronado-Zamora
- Genomics, Bioinformatics and Evolution, Departament de Genètica i de Microbiologia, Universitat Autònoma de Barcelona, Cerdanyola del Vallès, Spain
| | - Irepan Salvador-Martínez
- Evo-Devo Helsinki Community, Centre of Excellence in Experimental and Computational Developmental Biology, Institute of Biotechnology, University of Helsinki, Finland.,Department of Genetics, Evolution and Environment, University College London, United Kingdom
| | | | - Antonio Barbadilla
- Genomics, Bioinformatics and Evolution, Departament de Genètica i de Microbiologia, Universitat Autònoma de Barcelona, Cerdanyola del Vallès, Spain
| | - Isaac Salazar-Ciudad
- Genomics, Bioinformatics and Evolution, Departament de Genètica i de Microbiologia, Universitat Autònoma de Barcelona, Cerdanyola del Vallès, Spain.,Evo-Devo Helsinki Community, Centre of Excellence in Experimental and Computational Developmental Biology, Institute of Biotechnology, University of Helsinki, Finland.,Centre de Recerca Matemàtica, Cerdanyola del Vallès, Spain
| |
Collapse
|
18
|
Khatri BS, Goldstein RA. Biophysics and population size constrains speciation in an evolutionary model of developmental system drift. PLoS Comput Biol 2019; 15:e1007177. [PMID: 31335870 PMCID: PMC6677325 DOI: 10.1371/journal.pcbi.1007177] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/30/2018] [Revised: 08/02/2019] [Accepted: 06/13/2019] [Indexed: 02/06/2023] Open
Abstract
Developmental system drift is a likely mechanism for the origin of hybrid incompatibilities between closely related species. We examine here the detailed mechanistic basis of hybrid incompatibilities between two allopatric lineages, for a genotype-phenotype map of developmental system drift under stabilising selection, where an organismal phenotype is conserved, but the underlying molecular phenotypes and genotype can drift. This leads to number of emergent phenomenon not obtainable by modelling genotype or phenotype alone. Our results show that: 1) speciation is more rapid at smaller population sizes with a characteristic, Orr-like, power law, but at large population sizes slow, characterised by a sub-diffusive growth law; 2) the molecular phenotypes under weakest selection contribute to the earliest incompatibilities; and 3) pair-wise incompatibilities dominate over higher order, contrary to previous predictions that the latter should dominate. The population size effect we find is consistent with previous results on allopatric divergence of transcription factor-DNA binding, where smaller populations have common ancestors with a larger drift load because genetic drift favours phenotypes which have a larger number of genotypes (higher sequence entropy) over more fit phenotypes which have far fewer genotypes; this means less substitutions are required in either lineage before incompatibilities arise. Overall, our results indicate that biophysics and population size provide a much stronger constraint to speciation than suggested by previous models, and point to a general mechanistic principle of how incompatibilities arise the under stabilising selection for an organismal phenotype. The process of speciation is of fundamental importance to the field of evolution as it is intimately connected to understanding the immense bio-diversity of life. There is still relatively little understanding of the underlying genetic mechanisms that give rise to hybrid incompatibilities with results suggesting that divergence in transcription factor DNA binding and gene expression play an important role. A key finding from the field of evo-devo is that organismal phenotypes show developmental system drift, where species maintain the same phenotype, but diverge in developmental pathways; this is an important potential source of hybrid incompatibilities. Here, we explore a theoretical framework to understand how incompatibilities arise due to developmental system drift, using a tractable biophysically inspired genotype-phenotype for spatial gene expression. Modelling the evolution of phenotypes in this way has the key advantage that it mirrors how selection works in nature, i.e. that selection acts on phenotypes, but variation (mutation) arise at the level of genotypes. This results, as we demonstrate, in a number of non-trivial and testable predictions concerning speciation due to developmental system drift, which would not be obtainable by modelling evolution of genotypes or phenotypes alone.
Collapse
Affiliation(s)
| | - Richard A. Goldstein
- Division of Infection & Immunity, University College London, London, United Kingdom
| |
Collapse
|
19
|
Wunderlich Z, Fowlkes CC, Eckenrode KB, Bragdon MDJ, Abiri A, DePace AH. Quantitative Comparison of the Anterior-Posterior Patterning System in the Embryos of Five Drosophila Species. G3 (BETHESDA, MD.) 2019; 9:2171-2182. [PMID: 31048401 PMCID: PMC6643877 DOI: 10.1534/g3.118.200953] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/07/2018] [Accepted: 05/01/2019] [Indexed: 11/18/2022]
Abstract
Complex spatiotemporal gene expression patterns direct the development of the fertilized egg into an adult animal. Comparisons across species show that, in spite of changes in the underlying regulatory DNA sequence, developmental programs can be maintained across millions of years of evolution. Reciprocally, changes in gene expression can be used to generate morphological novelty. Distinguishing between changes in regulatory DNA that lead to changes in gene expression and those that do not is therefore a central goal of evolutionary developmental biology. Quantitative, spatially-resolved measurements of developmental gene expression patterns play a crucial role in this goal, enabling the detection of subtle phenotypic differences between species and the development of computations models that link the sequence of regulatory DNA to expression patterns. Here we report the generation of two atlases of cellular resolution gene expression measurements for the primary anterior-posterior patterning genes in Drosophila simulans and Drosophila virilis By combining these data sets with existing atlases for three other Drosophila species, we detect subtle differences in the gene expression patterns and dynamics driving the highly conserved axis patterning system and delineate inter-species differences in the embryonic morphology. These data sets will be a resource for future modeling studies of the evolution of developmental gene regulatory networks.
Collapse
Affiliation(s)
- Zeba Wunderlich
- Department of Developmental and Cell Biology, University of California, Irvine, CA, 92697
| | - Charless C Fowlkes
- Department of Computer Science, University of California, Irvine, CA, 92697
| | - Kelly B Eckenrode
- Department of Systems Biology, Harvard Medical School, Boston, MA, 20115
| | - Meghan D J Bragdon
- Department of Systems Biology, Harvard Medical School, Boston, MA, 20115
| | - Arash Abiri
- Department of Developmental and Cell Biology, University of California, Irvine, CA, 92697
| | - Angela H DePace
- Department of Systems Biology, Harvard Medical School, Boston, MA, 20115
| |
Collapse
|
20
|
Abstract
Proteins and RNA molecules are deposited into the developing egg by the mother. These gene products will drive the first stages of development and are coded by maternal genes. Maternal genes are essential, yet, despite their importance, their evolutionary dynamics is largely unknown. Here I review the current knowledge of maternal gene evolution. The evolutionary origin of maternal genes tends to be more recent than that of zygotic genes. Some studies support the theoretical prediction that maternal genes evolve faster than zygotic genes. However, most studies were done on a limited set of species and genes. I also discuss the way forward to understand the evolution of maternal genes by combining high-throughput genomics and theoretical evolutionary approaches.
Collapse
|
21
|
Verd B, Monk NAM, Jaeger J. Modularity, criticality, and evolvability of a developmental gene regulatory network. eLife 2019; 8:e42832. [PMID: 31169494 PMCID: PMC6645726 DOI: 10.7554/elife.42832] [Citation(s) in RCA: 47] [Impact Index Per Article: 9.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2019] [Accepted: 06/05/2019] [Indexed: 01/16/2023] Open
Abstract
The existence of discrete phenotypic traits suggests that the complex regulatory processes which produce them are functionally modular. These processes are usually represented by networks. Only modular networks can be partitioned into intelligible subcircuits able to evolve relatively independently. Traditionally, functional modularity is approximated by detection of modularity in network structure. However, the correlation between structure and function is loose. Many regulatory networks exhibit modular behaviour without structural modularity. Here we partition an experimentally tractable regulatory network-the gap gene system of dipteran insects-using an alternative approach. We show that this system, although not structurally modular, is composed of dynamical modules driving different aspects of whole-network behaviour. All these subcircuits share the same regulatory structure, but differ in components and sensitivity to regulatory interactions. Some subcircuits are in a state of criticality, while others are not, which explains the observed differential evolvability of the various expression features in the system.
Collapse
Affiliation(s)
- Berta Verd
- EMBL/CRG Systems Biology Research Unit, Centre for Genomic Regulation (CRG)The Barcelona Institute of Science and TechnologyBarcelonaSpain
- Universitat Pompeu Fabra (UPF)BarcelonaSpain
- Konrad Lorenz Institute for Evolution and Cognition Research (KLI)KlosterneuburgAustria
- Department of GeneticsUniversity of CambridgeCambridgeUnited Kingdom
| | - Nicholas AM Monk
- School of Mathematics and StatisticsUniversity of SheffieldSheffieldUnited States
| | - Johannes Jaeger
- EMBL/CRG Systems Biology Research Unit, Centre for Genomic Regulation (CRG)The Barcelona Institute of Science and TechnologyBarcelonaSpain
- Universitat Pompeu Fabra (UPF)BarcelonaSpain
- Konrad Lorenz Institute for Evolution and Cognition Research (KLI)KlosterneuburgAustria
- School of Mathematics and StatisticsUniversity of SheffieldSheffieldUnited States
- Wissenschaftskolleg zu BerlinBerlinGermany
- Center for Systems Biology Dresden (CSBD)DresdenGermany
- Complexity Science Hub (CSH)ViennaAustria
- Centre de Recherches Interdisciplinaires (CRI)ParisFrance
| |
Collapse
|
22
|
Ali S, Signor SA, Kozlov K, Nuzhdin SV. Novel approach to quantitative spatial gene expression uncovers genetic stochasticity in the developing Drosophila eye. Evol Dev 2019; 21:157-171. [PMID: 30756455 DOI: 10.1111/ede.12283] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Abstract
Robustness in development allows for the accumulation of genetically based variation in expression. However, this variation is usually examined in response to large perturbations, and examination of this variation has been limited to being spatial, or quantitative, but because of technical restrictions not both. Here we bridge these gaps by investigating replicated quantitative spatial gene expression using rigorous statistical models, in different genotypes, sexes, and species (Drosophila melanogaster and D. simulans). Using this type of quantitative approach with molecular developmental data allows for comparison among conditions, such as different genetic backgrounds. We apply this approach to the morphogenetic furrow, a wave of differentiation that patterns the developing eye disc. Within the morphogenetic furrow, we focus on four genes, hairy, atonal, hedgehog, and Delta. Hybridization chain reaction quantitatively measures spatial gene expression, co-staining for all four genes simultaneously. We find considerable variation in the spatial expression pattern of these genes in the eye between species, genotypes, and sexes. We also find that there has been evolution of the regulatory relationship between these genes, and that their spatial interrelationships have evolved between species. This variation has no phenotypic effect, and could be buffered by network thresholds or compensation from other genes. Both of these mechanisms could potentially be contributing to long term developmental systems drift.
Collapse
Affiliation(s)
- Sammi Ali
- Molecular and Computational Biology, University of Southern California, Los Angeles, California
| | - Sarah A Signor
- Molecular and Computational Biology, University of Southern California, Los Angeles, California
| | - Konstantin Kozlov
- Department of Applied Mathematics, St. Petersburg State Polytechnic University, St. Petersburg, Russia
| | - Sergey V Nuzhdin
- Molecular and Computational Biology, University of Southern California, Los Angeles, California.,Department of Applied Mathematics, St. Petersburg State Polytechnic University, St. Petersburg, Russia
| |
Collapse
|
23
|
Developmental Bias and Evolution: A Regulatory Network Perspective. Genetics 2018; 209:949-966. [PMID: 30049818 DOI: 10.1534/genetics.118.300995] [Citation(s) in RCA: 97] [Impact Index Per Article: 16.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2017] [Accepted: 04/19/2018] [Indexed: 01/12/2023] Open
Abstract
Phenotypic variation is generated by the processes of development, with some variants arising more readily than others-a phenomenon known as "developmental bias." Developmental bias and natural selection have often been portrayed as alternative explanations, but this is a false dichotomy: developmental bias can evolve through natural selection, and bias and selection jointly influence phenotypic evolution. Here, we briefly review the evidence for developmental bias and illustrate how it is studied empirically. We describe recent theory on regulatory networks that explains why the influence of genetic and environmental perturbation on phenotypes is typically not uniform, and may even be biased toward adaptive phenotypic variation. We show how bias produced by developmental processes constitutes an evolving property able to impose direction on adaptive evolution and influence patterns of taxonomic and phenotypic diversity. Taking these considerations together, we argue that it is not sufficient to accommodate developmental bias into evolutionary theory merely as a constraint on evolutionary adaptation. The influence of natural selection in shaping developmental bias, and conversely, the influence of developmental bias in shaping subsequent opportunities for adaptation, requires mechanistic models of development to be expanded and incorporated into evolutionary theory. A regulatory network perspective on phenotypic evolution thus helps to integrate the generation of phenotypic variation with natural selection, leaving evolutionary biology better placed to explain how organisms adapt and diversify.
Collapse
|
24
|
|
25
|
Pers D, Lynch JA. Ankyrin domain encoding genes from an ancient horizontal transfer are functionally integrated into Nasonia developmental gene regulatory networks. Genome Biol 2018; 19:148. [PMID: 30266092 PMCID: PMC6161386 DOI: 10.1186/s13059-018-1526-x] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/13/2018] [Accepted: 09/05/2018] [Indexed: 02/06/2023] Open
Abstract
BACKGROUND How regulatory networks incorporate additional components and how novel genes are functionally integrated into well-established developmental processes are two important and intertwined questions whose answers have major implications for understanding the evolution of development. We recently discovered a set of lineage-restricted genes with strong and specific expression patterns along the dorsal-ventral (DV) axis of the embryo of the wasp Nasonia that may serve as a powerful system for addressing these questions. We sought to both understand the evolutionary history of these genes and to determine their functions in the Nasonia DV patterning system. RESULTS We have found that the novel DV genes are part of a large family of rapidly duplicating and diverging ankyrin domain-encoding genes that originated most likely by horizontal transfer from a prokaryote in a common ancestor of the wasp superfamily Chalcidoidea. We tested the function of those ankyrin-encoding genes expressed along the DV axis and found that they participate in early embryonic DV patterning. We also developed a new wasp model system (Melittobia) and found that some functional integration of ankyrin genes have been preserved for over 90 million years. CONCLUSIONS Our results indicate that regulatory networks can incorporate novel genes that then become necessary for stable and repeatable outputs. Even a modest role in developmental networks may be enough to allow novel or duplicate genes to be maintained in the genome and become fully integrated network components.
Collapse
Affiliation(s)
- Daniel Pers
- Department of Biological Sciences, University of Illinois at Chicago, MBRB 4020, 900 S. Ashland Ave, Chicago, IL, 60607, USA
| | - Jeremy A Lynch
- Department of Biological Sciences, University of Illinois at Chicago, MBRB 4020, 900 S. Ashland Ave, Chicago, IL, 60607, USA.
| |
Collapse
|
26
|
Jiménez-Guri E, Wotton KR, Jaeger J. tarsal-less is expressed as a gap gene but has no gap gene phenotype in the moth midge Clogmia albipunctata. ROYAL SOCIETY OPEN SCIENCE 2018; 5:180458. [PMID: 30225035 PMCID: PMC6124123 DOI: 10.1098/rsos.180458] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 04/03/2018] [Accepted: 07/18/2018] [Indexed: 06/08/2023]
Abstract
Gap genes are involved in segment determination during early development of the vinegar fly Drosophila melanogaster and other dipteran insects (flies, midges and mosquitoes). They are expressed in overlapping domains along the antero-posterior (A-P) axis of the blastoderm embryo. While gap domains cover the entire length of the A-P axis in Drosophila, there is a region in the blastoderm of the moth midge Clogmia albipunctata, which lacks canonical gap gene expression. Is a non-canonical gap gene functioning in this area? Here, we characterize tarsal-less (tal) in C. albipunctata. The homologue of tal in the flour beetle Tribolium castaneum (called milles-pattes, mlpt) is a bona fide gap gene. We find that Ca-tal is expressed in the region previously reported as lacking gap gene expression. Using RNA interference, we study the interaction of Ca-tal with gap genes. We show that Ca-tal is regulated by gap genes, but only has a very subtle effect on tailless (Ca-tll), while not affecting other gap genes at all. Moreover, cuticle phenotypes of Ca-tal depleted embryos do not show any gap phenotype. We conclude that Ca-tal is expressed and regulated like a gap gene, but does not function as a gap gene in C. albipunctata.
Collapse
Affiliation(s)
- Eva Jiménez-Guri
- EMBL/CRG Research Unit in Systems Biology, Centre de Regulació Genòmica (CRG), The Barcelona Institute of Science and Technology (BIST), Dr. Aiguader 88, 08003 Barcelona, Spain
- Universitat Pompeu Fabra (UPF), Barcelona, Spain
- Centre for Ecology and Conservation, College of Life and Environmental Sciences, University of Exeter, Penryn, Cornwall TR10 9EZ, UK
| | - Karl R. Wotton
- EMBL/CRG Research Unit in Systems Biology, Centre de Regulació Genòmica (CRG), The Barcelona Institute of Science and Technology (BIST), Dr. Aiguader 88, 08003 Barcelona, Spain
- Universitat Pompeu Fabra (UPF), Barcelona, Spain
| | - Johannes Jaeger
- EMBL/CRG Research Unit in Systems Biology, Centre de Regulació Genòmica (CRG), The Barcelona Institute of Science and Technology (BIST), Dr. Aiguader 88, 08003 Barcelona, Spain
- Universitat Pompeu Fabra (UPF), Barcelona, Spain
| |
Collapse
|
27
|
Fraire-Zamora JJ, Jaeger J, Solon J. Two consecutive microtubule-based epithelial seaming events mediate dorsal closure in the scuttle fly Megaselia abdita. eLife 2018. [PMID: 29537962 PMCID: PMC5851697 DOI: 10.7554/elife.33807] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022] Open
Abstract
Evolution of morphogenesis is generally associated with changes in genetic regulation. Here, we report evidence indicating that dorsal closure, a conserved morphogenetic process in dipterans, evolved as the consequence of rearrangements in epithelial organization rather than signaling regulation. In Drosophila melanogaster, dorsal closure consists of a two-tissue system where the contraction of extraembryonic amnioserosa and a JNK/Dpp-dependent epidermal actomyosin cable result in microtubule-dependent seaming of the epidermis. We find that dorsal closure in Megaselia abdita, a three-tissue system comprising serosa, amnion and epidermis, differs in morphogenetic rearrangements despite conservation of JNK/Dpp signaling. In addition to an actomyosin cable, M. abdita dorsal closure is driven by the rupture and contraction of the serosa and the consecutive microtubule-dependent seaming of amnion and epidermis. Our study indicates that the evolutionary transition to a reduced system of dorsal closure involves simplification of the seaming process without changing the signaling pathways of closure progression.
Collapse
Affiliation(s)
- Juan Jose Fraire-Zamora
- Cell and Developmental Biology Programme, Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology, Barcelona, Spain.,Universitat Pompeu Fabra, Barcelona, Spain
| | - Johannes Jaeger
- Universitat Pompeu Fabra, Barcelona, Spain.,System Biology Programme, Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology, Barcelona, Spain.,Konrad Lorenz Institute for Evolution and Cognition Research (KLI), Klosterneuburg, Austria
| | - Jérôme Solon
- Cell and Developmental Biology Programme, Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology, Barcelona, Spain.,Universitat Pompeu Fabra, Barcelona, Spain
| |
Collapse
|
28
|
A damped oscillator imposes temporal order on posterior gap gene expression in Drosophila. PLoS Biol 2018; 16:e2003174. [PMID: 29451884 PMCID: PMC5832388 DOI: 10.1371/journal.pbio.2003174] [Citation(s) in RCA: 38] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2017] [Revised: 03/01/2018] [Accepted: 01/31/2018] [Indexed: 12/21/2022] Open
Abstract
Insects determine their body segments in two different ways. Short-germband insects, such as the flour beetle Tribolium castaneum, use a molecular clock to establish segments sequentially. In contrast, long-germband insects, such as the vinegar fly Drosophila melanogaster, determine all segments simultaneously through a hierarchical cascade of gene regulation. Gap genes constitute the first layer of the Drosophila segmentation gene hierarchy, downstream of maternal gradients such as that of Caudal (Cad). We use data-driven mathematical modelling and phase space analysis to show that shifting gap domains in the posterior half of the Drosophila embryo are an emergent property of a robust damped oscillator mechanism, suggesting that the regulatory dynamics underlying long- and short-germband segmentation are much more similar than previously thought. In Tribolium, Cad has been proposed to modulate the frequency of the segmentation oscillator. Surprisingly, our simulations and experiments show that the shift rate of posterior gap domains is independent of maternal Cad levels in Drosophila. Our results suggest a novel evolutionary scenario for the short- to long-germband transition and help explain why this transition occurred convergently multiple times during the radiation of the holometabolan insects. Different insect species exhibit one of two distinct modes of determining their body segments (known as segmentation) during development: they either use a molecular oscillator to position segments sequentially, or they generate segments simultaneously through a hierarchical gene-regulatory cascade. The sequential mode is ancestral, while the simultaneous mode has been derived from it independently several times during evolution. In this paper, we present evidence suggesting that simultaneous segmentation also involves an oscillator in the posterior end of the embryo of the vinegar fly, Drosophila melanogaster. This surprising result indicates that both modes of segment determination are much more similar than previously thought. Such similarity provides an important step towards our understanding of the frequent evolutionary transitions observed between sequential and simultaneous segmentation.
Collapse
|
29
|
Valfort AC, Launay C, Sémon M, Delattre M. Evolution of mitotic spindle behavior during the first asymmetric embryonic division of nematodes. PLoS Biol 2018; 16:e2005099. [PMID: 29357348 PMCID: PMC5794175 DOI: 10.1371/journal.pbio.2005099] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/13/2017] [Revised: 02/01/2018] [Accepted: 01/03/2018] [Indexed: 01/08/2023] Open
Abstract
Asymmetric cell division is essential to generate cellular diversity. In many animal cells, the cleavage plane lies perpendicular to the mitotic spindle, and it is the spindle positioning that dictates the size of the daughter cells. Although some properties of spindle positioning are conserved between distantly related model species and different cell types, little is known of the evolutionary robustness of the mechanisms underlying this event. We recorded the first embryonic division of 42 species of nematodes closely related to Caenorhabditis elegans, which is an excellent model system to study the biophysical properties of asymmetric spindle positioning. Our recordings, corresponding to 128 strains from 27 Caenorhabditis and 15 non-Caenorhabditis species (accessible at http://www.ens-lyon.fr/LBMC/NematodeCell/videos/), constitute a powerful collection of subcellular phenotypes to study the evolution of various cellular processes across species. In the present work, we analyzed our collection to the study of asymmetric spindle positioning. Although all the strains underwent an asymmetric first cell division, they exhibited large intra- and inter-species variations in the degree of cell asymmetry and in several parameters controlling spindle movement, including spindle oscillation, elongation, and displacement. Notably, these parameters changed frequently during evolution with no apparent directionality in the species phylogeny, with the exception of spindle transverse oscillations, which were an evolutionary innovation at the base of the Caenorhabditis genus. These changes were also unrelated to evolutionary variations in embryo size. Importantly, spindle elongation, displacement, and oscillation each evolved independently. This finding contrasts starkly with expectations based on C. elegans studies and reveals previously unrecognized evolutionary changes in spindle mechanics. Collectively, these data demonstrate that, while the essential process of asymmetric cell division has been conserved over the course of nematode evolution, the underlying spindle movement parameters can combine in various ways. Like other developmental processes, asymmetric cell division is subject to system drift.
Collapse
Affiliation(s)
- Aurore-Cécile Valfort
- Department of Pharmacology & Physiology (Colin Flaveny lab), Saint Louis University School of Medicine, Saint Louis, Missouri, United States of America
| | - Caroline Launay
- UnivLyon, ENS de Lyon, Univ Claude Bernard, Laboratory of Biology and Modelling of the Cell, Lyon University, Lyon, France
| | - Marie Sémon
- UnivLyon, ENS de Lyon, Univ Claude Bernard, Laboratory of Biology and Modelling of the Cell, Lyon University, Lyon, France
| | - Marie Delattre
- UnivLyon, ENS de Lyon, Univ Claude Bernard, Laboratory of Biology and Modelling of the Cell, Lyon University, Lyon, France
| |
Collapse
|
30
|
Jutras-Dubé L, Henry A, François P. Modelling Time-Dependent Acquisition of Positional Information. Methods Mol Biol 2018; 1863:281-301. [PMID: 30324604 DOI: 10.1007/978-1-4939-8772-6_16] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
Theoretical and computational modelling are crucial to understand dynamics of embryonic development. In this tutorial chapter, we describe two models of gene networks performing time-dependent acquisition of positional information under control of a dynamic morphogen: a toy-model of a bistable gene under control of a morphogen, allowing for the numerical computation of a simple Waddington's epigenetic landscape, and a recently published model of gap genes in Tribolium under control of multiple enhancers. We present detailed commented implementations of the models using python and jupyter notebooks.
Collapse
Affiliation(s)
- Laurent Jutras-Dubé
- McGill University, Ernest Rutherford Physics Building, 3600 rue University, H3A2T8, Montreal, QC, Canada
| | - Adrien Henry
- McGill University, Ernest Rutherford Physics Building, 3600 rue University, H3A2T8, Montreal, QC, Canada
| | - Paul François
- McGill University, Ernest Rutherford Physics Building, 3600 rue University, H3A2T8, Montreal, QC, Canada
| |
Collapse
|
31
|
Wotton KR, Alcaine-Colet A, Jaeger J, Jiménez-Guri E. Non-canonical dorsoventral patterning in the moth midge Clogmia albipunctata. EvoDevo 2017; 8:20. [PMID: 29158889 PMCID: PMC5683363 DOI: 10.1186/s13227-017-0083-9] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/24/2017] [Accepted: 11/03/2017] [Indexed: 11/20/2022] Open
Abstract
Background Bone morphogenetic proteins (BMPs) are of central importance for dorsal–ventral (DV) axis specification. They are core components of a signalling cascade that includes the BMP ligand decapentaplegic (DPP) and its antagonist short gastrulation (SOG) in Drosophila melanogaster. These components are very ancient, with orthologs involved in DV patterning in both protostomes and deuterostomes. Despite such strong conservation, recent comparative work in insects has revealed interesting differences in the way the patterning function of the DV system is achieved in different species. Results In this paper, we characterise the expression patterns of the principal components of the BMP DV patterning system, as well as its signalling outputs and downstream targets, in the non-cyclorrhaphan moth midge Clogmia albipunctata (Diptera: Psychodidae). We previously reported ventral expression patterns of dpp in the pole regions of C. albipunctata blastoderm embryos. Strikingly, we also find ventral sog and posteriorly restricted tkv expression, as well as expanded polar activity of pMad. We use our results from gene knock-down by embryonic RNA interference to propose a mechanism of polar morphogen shuttling in C. albipunctata. We compare these results to available data from other species and discuss scenarios for the evolution of DV signalling in the holometabolan insects. Conclusions A comparison of gene expression patterns across hemipteran and holometabolan insects reveals that expression of upstream signalling factors in the DV system is very variable, while signalling output is highly conserved. This has two major implications: first, as long as ligand shuttling and other upstream regulatory mechanisms lead to an appropriately localised activation of BMP signalling at the dorsal midline, it is of less importance exactly where the upstream components of the DV system are expressed. This, in turn, explains why the early-acting components of the DV patterning system in insects exhibit extensive amounts of developmental systems drift constrained by highly conserved downstream signalling output.
Collapse
Affiliation(s)
- Karl R Wotton
- EMBL/CRG Research Unit in Systems Biology, Centre de Regulació Genòmica (CRG), The Barcelona Institute of Science and Technology (BIST), Dr. Aiguader 88, 08003 Barcelona, Spain.,Universitat Pompeu Fabra (UPF), Barcelona, Spain.,Present Address: Centre for Ecology and Conservation, College of Life and Environmental Sciences, University of Exeter, Penryn, Cornwall TR10 9EZ UK
| | - Anna Alcaine-Colet
- EMBL/CRG Research Unit in Systems Biology, Centre de Regulació Genòmica (CRG), The Barcelona Institute of Science and Technology (BIST), Dr. Aiguader 88, 08003 Barcelona, Spain.,Universitat Pompeu Fabra (UPF), Barcelona, Spain
| | - Johannes Jaeger
- EMBL/CRG Research Unit in Systems Biology, Centre de Regulació Genòmica (CRG), The Barcelona Institute of Science and Technology (BIST), Dr. Aiguader 88, 08003 Barcelona, Spain.,Universitat Pompeu Fabra (UPF), Barcelona, Spain.,Present Address: Complexity Science Hub Vienna, Josefstädter Straße 39, 1080 Vienna, Austria
| | - Eva Jiménez-Guri
- EMBL/CRG Research Unit in Systems Biology, Centre de Regulació Genòmica (CRG), The Barcelona Institute of Science and Technology (BIST), Dr. Aiguader 88, 08003 Barcelona, Spain.,Universitat Pompeu Fabra (UPF), Barcelona, Spain.,Present Address: Centre for Ecology and Conservation, College of Life and Environmental Sciences, University of Exeter, Penryn, Cornwall TR10 9EZ UK
| |
Collapse
|
32
|
Myasnikova EM, Spirov AV. A Method for Estimating the Predictive Power in a Model of a Biological System with Low Sensitivity to Parameters. Biophysics (Nagoya-shi) 2017. [DOI: 10.1134/s0006350917060185] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022] Open
|
33
|
Gursky VV, Kozlov KN, Kulakovskiy IV, Zubair A, Marjoram P, Lawrie DS, Nuzhdin SV, Samsonova MG. Translating natural genetic variation to gene expression in a computational model of the Drosophila gap gene regulatory network. PLoS One 2017; 12:e0184657. [PMID: 28898266 PMCID: PMC5595321 DOI: 10.1371/journal.pone.0184657] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/05/2017] [Accepted: 08/28/2017] [Indexed: 11/18/2022] Open
Abstract
Annotating the genotype-phenotype relationship, and developing a proper quantitative description of the relationship, requires understanding the impact of natural genomic variation on gene expression. We apply a sequence-level model of gap gene expression in the early development of Drosophila to analyze single nucleotide polymorphisms (SNPs) in a panel of natural sequenced D. melanogaster lines. Using a thermodynamic modeling framework, we provide both analytical and computational descriptions of how single-nucleotide variants affect gene expression. The analysis reveals that the sequence variants increase (decrease) gene expression if located within binding sites of repressors (activators). We show that the sign of SNP influence (activation or repression) may change in time and space and elucidate the origin of this change in specific examples. The thermodynamic modeling approach predicts non-local and non-linear effects arising from SNPs, and combinations of SNPs, in individual fly genotypes. Simulation of individual fly genotypes using our model reveals that this non-linearity reduces to almost additive inputs from multiple SNPs. Further, we see signatures of the action of purifying selection in the gap gene regulatory regions. To infer the specific targets of purifying selection, we analyze the patterns of polymorphism in the data at two phenotypic levels: the strengths of binding and expression. We find that combinations of SNPs show evidence of being under selective pressure, while individual SNPs do not. The model predicts that SNPs appear to accumulate in the genotypes of the natural population in a way biased towards small increases in activating action on the expression pattern. Taken together, these results provide a systems-level view of how genetic variation translates to the level of gene regulatory networks via combinatorial SNP effects.
Collapse
Affiliation(s)
- Vitaly V. Gursky
- Theoretical Department, Ioffe Institute, Saint Petersburg, Russia
- Systems Biology and Bioinformatics Laboratory, Peter the Great Saint Petersburg Polytechnic University, Saint Petersburg, Russia
- * E-mail:
| | - Konstantin N. Kozlov
- Systems Biology and Bioinformatics Laboratory, Peter the Great Saint Petersburg Polytechnic University, Saint Petersburg, Russia
| | - Ivan V. Kulakovskiy
- Engelhardt Institute of Molecular Biology, Moscow, Russia
- Vavilov Institute of General Genetics, Moscow, Russia
- Center for Data-Intensive Biomedicine and Biotechnology, Skolkovo Institute of Science and Technology, Moscow, Russia
| | - Asif Zubair
- Molecular and Computational Biology, University of Southern California, Los Angeles, California, United States of America
| | - Paul Marjoram
- Molecular and Computational Biology, University of Southern California, Los Angeles, California, United States of America
| | - David S. Lawrie
- Molecular and Computational Biology, University of Southern California, Los Angeles, California, United States of America
| | - Sergey V. Nuzhdin
- Molecular and Computational Biology, University of Southern California, Los Angeles, California, United States of America
| | - Maria G. Samsonova
- Systems Biology and Bioinformatics Laboratory, Peter the Great Saint Petersburg Polytechnic University, Saint Petersburg, Russia
| |
Collapse
|
34
|
Clark E. Dynamic patterning by the Drosophila pair-rule network reconciles long-germ and short-germ segmentation. PLoS Biol 2017; 15:e2002439. [PMID: 28953896 PMCID: PMC5633203 DOI: 10.1371/journal.pbio.2002439] [Citation(s) in RCA: 71] [Impact Index Per Article: 10.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/13/2017] [Revised: 10/09/2017] [Accepted: 09/07/2017] [Indexed: 02/07/2023] Open
Abstract
Drosophila segmentation is a well-established paradigm for developmental pattern formation. However, the later stages of segment patterning, regulated by the "pair-rule" genes, are still not well understood at the system level. Building on established genetic interactions, I construct a logical model of the Drosophila pair-rule system that takes into account the demonstrated stage-specific architecture of the pair-rule gene network. Simulation of this model can accurately recapitulate the observed spatiotemporal expression of the pair-rule genes, but only when the system is provided with dynamic "gap" inputs. This result suggests that dynamic shifts of pair-rule stripes are essential for segment patterning in the trunk and provides a functional role for observed posterior-to-anterior gap domain shifts that occur during cellularisation. The model also suggests revised patterning mechanisms for the parasegment boundaries and explains the aetiology of the even-skipped null mutant phenotype. Strikingly, a slightly modified version of the model is able to pattern segments in either simultaneous or sequential modes, depending only on initial conditions. This suggests that fundamentally similar mechanisms may underlie segmentation in short-germ and long-germ arthropods.
Collapse
Affiliation(s)
- Erik Clark
- Laboratory for Development and Evolution, Department of Zoology, University of Cambridge, Cambridge, United Kingdom
| |
Collapse
|
35
|
Verd B, Crombach A, Jaeger J. Dynamic Maternal Gradients Control Timing and Shift-Rates for Drosophila Gap Gene Expression. PLoS Comput Biol 2017; 13:e1005285. [PMID: 28158178 PMCID: PMC5291410 DOI: 10.1371/journal.pcbi.1005285] [Citation(s) in RCA: 45] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/22/2016] [Accepted: 12/06/2016] [Indexed: 11/18/2022] Open
Abstract
Pattern formation during development is a highly dynamic process. In spite of this, few experimental and modelling approaches take into account the explicit time-dependence of the rules governing regulatory systems. We address this problem by studying dynamic morphogen interpretation by the gap gene network in Drosophila melanogaster. Gap genes are involved in segment determination during early embryogenesis. They are activated by maternal morphogen gradients encoded by bicoid (bcd) and caudal (cad). These gradients decay at the same time-scale as the establishment of the antero-posterior gap gene pattern. We use a reverse-engineering approach, based on data-driven regulatory models called gene circuits, to isolate and characterise the explicitly time-dependent effects of changing morphogen concentrations on gap gene regulation. To achieve this, we simulate the system in the presence and absence of dynamic gradient decay. Comparison between these simulations reveals that maternal morphogen decay controls the timing and limits the rate of gap gene expression. In the anterior of the embyro, it affects peak expression and leads to the establishment of smooth spatial boundaries between gap domains. In the posterior of the embryo, it causes a progressive slow-down in the rate of gap domain shifts, which is necessary to correctly position domain boundaries and to stabilise the spatial gap gene expression pattern. We use a newly developed method for the analysis of transient dynamics in non-autonomous (time-variable) systems to understand the regulatory causes of these effects. By providing a rigorous mechanistic explanation for the role of maternal gradient decay in gap gene regulation, our study demonstrates that such analyses are feasible and reveal important aspects of dynamic gene regulation which would have been missed by a traditional steady-state approach. More generally, it highlights the importance of transient dynamics for understanding complex regulatory processes in development. Animal development is a highly dynamic process. Biochemical or environmental signals can cause the rules that shape it to change over time. We know little about the effects of such changes. For the sake of simplicity, we usually leave them out of our models and experimental assays. Here, we do exactly the opposite. We characterise precisely those aspects of pattern formation caused by changing signalling inputs to a gene regulatory network, the gap gene system of Drosophila melanogaster. Gap genes are involved in determining the body segments of flies and other insects during early development. Gradients of maternal morphogens activate the expression of the gap genes. These gradients are highly dynamic themselves, as they decay while being read out. We show that this decay controls the peak concentration of gap gene products, produces smooth boundaries of gene expression, and slows down the observed positional shifts of gap domains in the posterior of the embryo, thereby stabilising the spatial pattern. Our analysis demonstrates that the dynamics of gene regulation not only affect the timing, but also the positioning of gene expression. This suggests that we must pay closer attention to transient dynamic aspects of development than is currently the case.
Collapse
Affiliation(s)
- Berta Verd
- EMBL/CRG Systems Biology Research Unit, Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology, Barcelona, Spain
- Universitat Pompeu Fabra (UPF), Barcelona, Spain
- KLI Klosterneuburg, Klosterneuburg, Austria
- * E-mail: (BV); (JJ)
| | - Anton Crombach
- EMBL/CRG Systems Biology Research Unit, Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology, Barcelona, Spain
- Universitat Pompeu Fabra (UPF), Barcelona, Spain
- Center for Interdisciplinary Research in Biology (CIRB), Collège de France, CNRS, INSERM, PSL Research University, Paris, France
| | - Johannes Jaeger
- EMBL/CRG Systems Biology Research Unit, Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology, Barcelona, Spain
- Universitat Pompeu Fabra (UPF), Barcelona, Spain
- KLI Klosterneuburg, Klosterneuburg, Austria
- Wissenschaftskolleg zu Berlin, Berlin, Germany
- * E-mail: (BV); (JJ)
| |
Collapse
|
36
|
Surkova SY, Golubkova EV, Mamon LA, Samsonova MG. Morphogenetic networks which determine the spatial expression of zygotic genes in early Drosophila embryo. Russ J Dev Biol 2016. [DOI: 10.1134/s1062360416040093] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
37
|
Rothschild JB, Tsimiklis P, Siggia ED, François P. Predicting Ancestral Segmentation Phenotypes from Drosophila to Anopheles Using In Silico Evolution. PLoS Genet 2016; 12:e1006052. [PMID: 27227405 PMCID: PMC4882032 DOI: 10.1371/journal.pgen.1006052] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/24/2016] [Accepted: 04/23/2016] [Indexed: 12/23/2022] Open
Abstract
Molecular evolution is an established technique for inferring gene homology but regulatory DNA turns over so rapidly that inference of ancestral networks is often impossible. In silico evolution is used to compute the most parsimonious path in regulatory space for anterior-posterior patterning linking two Dipterian species. The expression pattern of gap genes has evolved between Drosophila (fly) and Anopheles (mosquito), yet one of their targets, eve, has remained invariant. Our model predicts that stripe 5 in fly disappears and a new posterior stripe is created in mosquito, thus eve stripe modules 3+7 and 4+6 in fly are homologous to 3+6 and 4+5 in mosquito. We can place Clogmia on this evolutionary pathway and it shares the mosquito homologies. To account for the evolution of the other pair-rule genes in the posterior we have to assume that the ancestral Dipterian utilized a dynamic method to phase those genes in relation to eve. The last common ancestor of the fruit fly (Drosophila) and mosquito (Anopheles) lived more than 200 Million years ago. Can we use available data on insects alive today to infer what their ancestor looked like? In this manuscript, we focus on early embryonic development, when stripes of genetic expression appear and define the location of insect segments (“segmentation”). We use an evolutionary algorithm to reconstruct and predict dynamics of genes controlling stripes in the last common ancestor of fly and mosquito. We predict a new and different combinatorial logic of stripe formation in mosquito compared to fly, which is fully consistent with development of intermediate species such as moth-fly (Clogmia). Our simulations further suggest that the dynamics of gene expression in this last common ancestor were similar to other insects, such as wasps (Nasonia). Our method illustrates how computational methods inspired by machine learning and non-linear physics can be used to infer gene dynamics in species that disappeared millions of years ago.
Collapse
Affiliation(s)
- Jeremy B. Rothschild
- Physics Department, McGill University, Ernest Rutherford Physics Building, Montreal, Quebec, Canada
| | - Panagiotis Tsimiklis
- Physics Department, McGill University, Ernest Rutherford Physics Building, Montreal, Quebec, Canada
| | - Eric D. Siggia
- Center for Studies in Physics and Biology, The Rockefeller University, New York, New York, United States of America
| | - Paul François
- Physics Department, McGill University, Ernest Rutherford Physics Building, Montreal, Quebec, Canada
- * E-mail:
| |
Collapse
|
38
|
Crombach A, Wotton KR, Jiménez-Guri E, Jaeger J. Gap Gene Regulatory Dynamics Evolve along a Genotype Network. Mol Biol Evol 2016; 33:1293-307. [PMID: 26796549 PMCID: PMC4839219 DOI: 10.1093/molbev/msw013] [Citation(s) in RCA: 51] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Abstract
Developmental gene networks implement the dynamic regulatory mechanisms that pattern and shape the organism. Over evolutionary time, the wiring of these networks changes, yet the patterning outcome is often preserved, a phenomenon known as “system drift.” System drift is illustrated by the gap gene network—involved in segmental patterning—in dipteran insects. In the classic model organism Drosophila melanogaster and the nonmodel scuttle fly Megaselia abdita, early activation and placement of gap gene expression domains show significant quantitative differences, yet the final patterning output of the system is essentially identical in both species. In this detailed modeling analysis of system drift, we use gene circuits which are fit to quantitative gap gene expression data in M. abdita and compare them with an equivalent set of models from D. melanogaster. The results of this comparative analysis show precisely how compensatory regulatory mechanisms achieve equivalent final patterns in both species. We discuss the larger implications of the work in terms of “genotype networks” and the ways in which the structure of regulatory networks can influence patterns of evolutionary change (evolvability).
Collapse
Affiliation(s)
- Anton Crombach
- EMBL/CRG Systems Biology Research Unit, Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology, Dr. Aiguader 88, Barcelona 08003, Spain Universitat Pompeu Fabra (UPF), Barcelona, Spain
| | - Karl R Wotton
- EMBL/CRG Systems Biology Research Unit, Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology, Dr. Aiguader 88, Barcelona 08003, Spain Universitat Pompeu Fabra (UPF), Barcelona, Spain
| | - Eva Jiménez-Guri
- EMBL/CRG Systems Biology Research Unit, Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology, Dr. Aiguader 88, Barcelona 08003, Spain Universitat Pompeu Fabra (UPF), Barcelona, Spain
| | - Johannes Jaeger
- EMBL/CRG Systems Biology Research Unit, Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology, Dr. Aiguader 88, Barcelona 08003, Spain Universitat Pompeu Fabra (UPF), Barcelona, Spain
| |
Collapse
|
39
|
Wotton KR, Jiménez-Guri E, Jaeger J. Maternal co-ordinate gene regulation and axis polarity in the scuttle fly Megaselia abdita. PLoS Genet 2015; 11:e1005042. [PMID: 25757102 PMCID: PMC4355411 DOI: 10.1371/journal.pgen.1005042] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2014] [Accepted: 01/30/2015] [Indexed: 02/01/2023] Open
Abstract
Axis specification and segment determination in dipteran insects are an excellent model system for comparative analyses of gene network evolution. Antero-posterior polarity of the embryo is established through systems of maternal morphogen gradients. In Drosophila melanogaster, the anterior system acts through opposing gradients of Bicoid (Bcd) and Caudal (Cad), while the posterior system involves Nanos (Nos) and Hunchback (Hb) protein. These systems act redundantly. Both Bcd and Hb need to be eliminated to cause a complete loss of polarity resulting in mirror-duplicated abdomens, so-called bicaudal phenotypes. In contrast, knock-down of bcd alone is sufficient to induce double abdomens in non-drosophilid cyclorrhaphan dipterans such as the hoverfly Episyrphus balteatus or the scuttle fly Megaselia abdita. We investigate conserved and divergent aspects of axis specification in the cyclorrhaphan lineage through a detailed study of the establishment and regulatory effect of maternal gradients in M. abdita. Our results show that the function of the anterior maternal system is highly conserved in this species, despite the loss of maternal cad expression. In contrast, hb does not activate gap genes in this species. The absence of this activatory role provides a precise genetic explanation for the loss of polarity upon bcd knock-down in M. abdita, and suggests a general scenario in which the posterior maternal system is increasingly replaced by the anterior one during the evolution of the cyclorrhaphan dipteran lineage. The basic head-to-tail polarity of an animal is established very early in development. In dipteran insects (flies, midges, and mosquitoes), polarity is established with the help of so-called morphogen gradients. Morphogens are regulatory proteins that are distributed as a concentration gradient, often involving diffusion from a localised source. This graded distribution then leads to the concentration-dependent activation of different target genes along the embryo’s axis. We examine this process, which differs to a surprising extent between dipteran species, in the scuttle fly Megaselia abdita, and compare our results to the model organism Drosophila melanogaster. In this way, we not only gain insights into how the mechanisms that establish polarity function differently in different species, but also how the system has evolved since these two flies shared a common ancestor. Specifically, we pin down the main difference between Drosophila and Megaselia in the altered function of the maternal Hunchback morphogen gradient, which activates target genes in the former, but not the latter species, where it has been completely replaced by the Bicoid morphogen during evolution.
Collapse
Affiliation(s)
- Karl R. Wotton
- EMBL/CRG Research Unit in Systems Biology, Centre for Genomic Regulation (CRG), Barcelona, Spain
- Universitat Pompeu Fabra (UPF), Barcelona, Spain
- * E-mail: (KW); (JJ)
| | - Eva Jiménez-Guri
- EMBL/CRG Research Unit in Systems Biology, Centre for Genomic Regulation (CRG), Barcelona, Spain
- Universitat Pompeu Fabra (UPF), Barcelona, Spain
| | - Johannes Jaeger
- EMBL/CRG Research Unit in Systems Biology, Centre for Genomic Regulation (CRG), Barcelona, Spain
- Universitat Pompeu Fabra (UPF), Barcelona, Spain
- * E-mail: (KW); (JJ)
| |
Collapse
|
40
|
Wotton KR, Jiménez-Guri E, Crombach A, Cicin-Sain D, Jaeger J. High-resolution gene expression data from blastoderm embryos of the scuttle fly Megaselia abdita. Sci Data 2015; 2:150005. [PMID: 25977812 PMCID: PMC4423355 DOI: 10.1038/sdata.2015.5] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/19/2014] [Accepted: 02/03/2015] [Indexed: 12/21/2022] Open
Abstract
Gap genes are involved in segment determination during early development in dipteran insects (flies, midges, and mosquitoes). We carried out a systematic quantitative comparative analysis of the gap gene network across different dipteran species. Our work provides mechanistic insights into the evolution of this pattern-forming network. As a central component of our project, we created a high-resolution quantitative spatio-temporal data set of gap and maternal co-ordinate gene expression in the blastoderm embryo of the non-drosophilid scuttle fly, Megaselia abdita. Our data include expression patterns in both wild-type and RNAi-treated embryos. The data-covering 10 genes, 10 time points, and over 1,000 individual embryos-consist of original embryo images, quantified expression profiles, extracted positions of expression boundaries, and integrated expression patterns, plus metadata and intermediate processing steps. These data provide a valuable resource for researchers interested in the comparative study of gene regulatory networks and pattern formation, an essential step towards a more quantitative and mechanistic understanding of developmental evolution.
Collapse
Affiliation(s)
- Karl R Wotton
- EMBL/CRG Research Unit in Systems Biology, Centre for Genomic Regulation (CRG), 08003 Barcelona, Spain
- Universitat Pompeu Fabra (UPF), 08002 Barcelona, Spain
| | - Eva Jiménez-Guri
- EMBL/CRG Research Unit in Systems Biology, Centre for Genomic Regulation (CRG), 08003 Barcelona, Spain
- Universitat Pompeu Fabra (UPF), 08002 Barcelona, Spain
| | - Anton Crombach
- EMBL/CRG Research Unit in Systems Biology, Centre for Genomic Regulation (CRG), 08003 Barcelona, Spain
- Universitat Pompeu Fabra (UPF), 08002 Barcelona, Spain
| | - Damjan Cicin-Sain
- EMBL/CRG Research Unit in Systems Biology, Centre for Genomic Regulation (CRG), 08003 Barcelona, Spain
- Universitat Pompeu Fabra (UPF), 08002 Barcelona, Spain
| | - Johannes Jaeger
- EMBL/CRG Research Unit in Systems Biology, Centre for Genomic Regulation (CRG), 08003 Barcelona, Spain
- Universitat Pompeu Fabra (UPF), 08002 Barcelona, Spain
| |
Collapse
|