Dynamic Maternal Gradients Control Timing and Shift-Rates for Drosophila Gap Gene Expression

PLoS Comput Biol. 2017 Feb 3;13(2):e1005285. doi: 10.1371/journal.pcbi.1005285. eCollection 2017 Feb.

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.

MeSH terms

  • Animals
  • Body Patterning / physiology*
  • Computer Simulation
  • Drosophila / embryology
  • Drosophila / physiology*
  • Drosophila Proteins / metabolism*
  • Embryonic Development / physiology
  • Female
  • Gene Expression Regulation, Developmental / physiology*
  • Homeodomain Proteins / metabolism*
  • Models, Biological*
  • Trans-Activators / metabolism*
  • Transcription Factors / metabolism*

Substances

  • Drosophila Proteins
  • Homeodomain Proteins
  • Trans-Activators
  • Transcription Factors
  • bcd protein, Drosophila
  • cad protein, Drosophila

Grants and funding

The authors thankfully acknowledge the computer resources, technical expertise and assistance provided by the Barcelona Supercomputing Center—Centro Nacional de Supercomputación. BV was supported by a “la Caixa” fellowship, and a Writing-Up Fellowship at the KLI Klosterneuburg. The research group of JJ was supported by the MEC-EMBL agreement for the EMBL/CRG Research Unit in Systems Biology, European Commission grant FP7-KBBE-2011-5/289434 (BioPreDyn), and grants BFU2009-10184 and BFU2012-33775 from MINECO. JJ thanks the Wissenschaftskolleg zu Berlin (Wiko) for a 10-month fellowship in 2014/15 and AC, for a half-year fellowship in 2014/15. The Centre for Genomic Regulation (CRG) acknowledges support from MINECO, “Centro de Excelencia Severo Ochoa 2013–2017” SEV-2012-0208. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.