Gap Gene Research Articles

Expression of distal limb patterning genes in Hypsibius exemplaris indicate regionalization and suggest distal identity of tardigrade legs.

Panarthropods, a major group of invertebrate animals comprised of arthropods, onychophorans, and tardigrades, are the only limb-bearing members of Ecdysozoa. The complexity and versatility of panarthropod paired limbs has prompted great interest in their development to better understand the formation of these structures and the genes involved in this process. However, studies of limb patterning and development are overwhelmingly focused on arthropods, followed by select work on onychophorans but almost entirely lacking for tardigrades. This model organism bias is inherently limited and precludes a comparative analysis of how panarthropod legs originated, have evolved, and the likely limb patterning genes present in the earliest panarthropod ancestors. In this study, we investigated tardigrade homologs of seven arthropod distal limb patterning genes (apterous, aristaless, BarH1, clawless, Lim1, rotund, and spineless) to better characterize tardigrade limb development in a comparative context. We detected homologs of all seven genes in the eutardigrade Hypsibius exemplaris and heterotardigrade Echiniscoides cf. sigismundi suggesting their conservation in both tardigrade lineages. Hybridization chain reaction experiments in H. exemplaris reveal a regionalized expression pattern for the genes aristaless, BarH1, clawless, rotund and spineless. The observed regionalized expression of the distal limb patterning genes in H. exemplaris might reflect the external morphological features of tardigrade legs, such as the distal claws, sensory organs in the proximal region, and specific muscle attachment sites. The comparison between the expression of these limb patterning genes in H. exemplaris relative to other panarthropods suggests their conserved role in the last common panarthropod ancestor, such as establishing the distal limb end and the distribution of sensory structures. Our results support the hypothesis that tardigrade legs are homologous to the distal region of other panarthropod limbs, as suggested by previous work on the expression of leg gap genes in H. exemplaris.

Body plan evolvability: The role of variability in gene regulatory networks.

Recent computational modeling of early fruit fly (Drosophila) development has characterized the degree to which gene regulation networks can be robust to natural variability. In the first few hours of development, broad spatial gradients of maternally derived transcription factors activate embryonic gap genes. These gap patterns determine the subsequent segmented insect body plan through pair-rule gene expression. Gap genes are expressed with greater spatial precision than the maternal patterns. Computational modeling of the gap-gap regulatory interactions provides a mechanistic understanding for this robustness to maternal variability in wild-type (WT) patterning. A long-standing question in evolutionary biology has been how a system which is robust, such as the developmental program creating any particular species' body plan, is also evolvable, i.e. how can a system evolve or speciate, if the WT form is strongly buffered and protected? In the present work, we use the WT model to explore the breakdown of such Waddington-type 'canalization'. What levels of variability will push the system out of the WT form; are there particular pathways in the gene regulatory mechanism which are more susceptible to losing the WT form; and when robustness is lost, what types of forms are most likely to occur (i.e. what forms lie near the WT)? Manipulating maternal effects in several different pathways, we find a common gap 'peak-to-step' pattern transition in the loss of WT.We discuss these results in terms of the evolvability of insect segmentation, and in terms of experimental perturbations and mutations which could test the model predictions. We conclude by discussing the prospects for using continuum models of pattern dynamics to investigate a wider range of evo-devo problems.

The gap gene of Rhizobium etli is required for both free life and symbiosis with common beans.

Glyceraldehyde-3-phosphate dehydrogenase (GAPDH or Gap) is a ubiquitous enzyme essential for carbon and energy metabolism in most organisms. Despite its primary role in sugar metabolism, GAPDH is recognized for its involvement in diverse cellular processes, being considered a paradigm among multifunctional/moonlighting proteins. Besides its canonical cytoplasmic location, GAPDH has been detected on cell surfaces or as a secreted protein in prokaryotes, yet little is known about its possible roles in plant symbiotic bacteria. Here we report that Rhizobium etli, a nitrogen-fixing symbiont of common beans, carries a single gap gene responsible for both GAPDH glycolytic and gluconeogenic activities. An active Gap protein is required throughout all stages of the symbiosis between R. etli and its host plant Phaseolus vulgaris. Both glycolytic and gluconeogenic Gap metabolic activities likely contribute to bacterial fitness during early and intermediate stages of the interaction, whereas GAPDH gluconeogenic activity seems critical for nodule invasion and nitrogen fixation. Although the R. etli Gap protein is secreted in a c-di-GMP related manner, no involvement of the R. etli gap gene in c-di-GMP related phenotypes, such as flocculation, biofilm formation or EPS production, was observed. Notably, the R. etli gap gene fully complemented a double gap1/gap2 mutant of Pseudomonas syringae for free life growth, albeit only partially in planta, suggesting potential specific roles for each type of Gap protein. Nevertheless, further research is required to unravel additional functions of the R. etli Gap protein beyond its essential metabolic roles.

Finding the Last Bits of Positional Information

In a developing embryo, information about the position of cells is encoded in the concentrations of morphogen molecules. In the fruit fly, the local concentrations of just a handful of proteins encoded by the gap genes are sufficient to specify position with a precision comparable to the spacing between cells along the anterior-posterior axis. This matches the precision of downstream events such as the striped patterns of expression in the pair-rule genes, but is not quite sufficient to define unique identities for individual cells. We demonstrate theoretically that this information gap can be bridged if positional errors are spatially correlated, with correlation lengths approximately 20% of the embryo length. We then show experimentally that these correlations are present, with the required strength, in the fluctuating positions of the pair-rule stripes, and this can be traced back to the gap genes. Taking account of these correlations, the available information matches the information needed for unique cellular specification, within error bars of approximately 2%. These observation support a precisionist view of information flow through the underlying genetic networks, in which accurate signals are available from the start and preserved as they are transformed into the final spatial patterns. Published by the American Physical Society 2024

Multifaceted effects on even-skipped transcriptional dynamics upon Krüppel dosage changes.

Although fluctuations in transcription factor (TF) dosage are often well tolerated, TF dosage modulation can change the target gene expression dynamics and result in significant non-lethal developmental phenotypes. Using MS2/MCP-mediated quantitative live imaging in early Drosophila embryos, we analyzed how changing levels of the gap gene Krüppel (Kr) affects transcriptional dynamics of the pair-rule gene even-skipped (eve). Halving the Kr dosage leads to a transient posterior expansion of the eve stripe 2 and an anterior shift of stripe 5. Surprisingly, the most significant changes are observed in eve stripes 3 and 4, the enhancers of which do not contain Kr-binding sites. In Kr heterozygous embryos, both stripes 3 and 4 display narrower widths, anteriorly shifted boundaries and reduced mRNA production levels. We show that Kr dosage indirectly affects stripe 3 and 4 dynamics by modulating other gap gene dynamics. We quantitatively correlate moderate body segment phenotypes of Kr heterozygotes with spatiotemporal changes in eve expression. Our results indicate that nonlinear relationships between TF dosage and phenotypes underlie direct TF-DNA and indirect TF-TF interactions.

Gene regulatory patterning codes in early cell fate specification of the C. elegans embryo.

Pattern formation originates during embryogenesis by a series of symmetry-breaking steps throughout an expanding cell lineage. In Drosophila, classic work has shown that segmentation in the embryo is established by morphogens within a syncytium, and the subsequent action of the gap, pair-rule, and segment polarity genes. This classic model however does not translate directly to species that lack a syncytium - such as Caenorhabditis elegans - where cell fate is specified by cell-autonomous cell lineage programs and their inter-signaling. Previous single-cell RNA-Seq studies in C. elegans have analyzed cells from a mixed suspension of cells from many embryos to study late differentiation stages, or individual early stage embryos to study early gene expression in the embryo. To study the intermediate stages of early and late gastrulation (28- to 102-cells stages) missed by these approaches, here we determine the transcriptomes of the 1- to 102-cell stage to identify 119 embryonic cell states during cell fate specification, including 'equivalence-group' cell identities. We find that gene expression programs are modular according to the sub-cell lineages, each establishing a set of stripes by combinations of transcription factor gene expression across the anterior-posterior axis. In particular, expression of the homeodomain genes establishes a comprehensive lineage-specific positioning system throughout the embryo beginning at the 28-cell stage. Moreover, we find that genes that segment the entire embryo in Drosophila have orthologs in C. elegans that exhibit sub-lineage-specific expression. These results suggest that the C. elegans embryo is patterned by a juxtaposition of distinct lineage-specific gene regulatory programs each with a unique encoding of cell location and fate. This use of homologous gene regulatory patterning codes suggests a deep homology of cell fate specification programs across diverse modes of development.

Quantifying robustness of the gap gene network

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.

Identification and Characterization of Arcanobacterium canis from Companion Animals in Germany and The United Kingdom

Arcanobacterium canis is a novel species of the Arcanobacterium most closely related to A. haemolyticum. This study aims to characterize two A. canis isolates recovered from companion animals, specifically the claw of a cat and a vaginal swab from a dog. This study used real-time PCR to characterize A. canis isolated from companion animals. Two isolates of A. canis were recovered from purulent material from the claw of an 11-year-old cat in Germany and a vaginal swab of a dog in the United Kingdom. The samples were characterized phenotypically and genotypically. Both isolates were analyzed using culture methods, biochemical analysis, MALDI-TOF MS, real-time PCR amplification and sequencing of the 16S rRNA gene, and rpoB, gap, and tuf genes. The findings showed that the isolates P5197-15 and M214-96-1 obtained from companion animals were successfully characterized and confirmed to species level by real-time PCR amplification and sequencing of the 16S rRNA gene, as well as the genes of rpoB, gap, and tuf. This study seeks to comprehensively understand the characteristics of A. canis isolates obtained from companion animals. Such knowledge is essential for accurate diagnosis, treatment, and control of infections caused by this pathogen in veterinary medicine. Additionally, it contributes to the broader understanding of the genetic diversity and characteristics of A. canis, which can have implications for public health and animal well-being.

Multi-level regulation of even-skipped stripes by the ubiquitous factor Zelda.

The zinc-finger protein Zelda (Zld) is a key activator of zygotic transcription in early Drosophila embryos. Here, we study Zld-dependent regulation of the seven-striped pattern of the pair-rule gene even-skipped (eve). Individual stripes are regulated by discrete enhancers that respond to broadly distributed activators; stripe boundaries are formed by localized repressors encoded by the gap genes. The strongest effects of Zld are on stripes 2, 3 and 7, which are regulated by two enhancers in a 3.8 kb genomic fragment that includes the eve basal promoter. We show that Zld facilitates binding of the activator Bicoid and the gap repressors to this fragment, consistent with its proposed role as a pioneer protein. To test whether the effects of Zld are direct, we mutated all canonical Zld sites in the 3.8 kb fragment, which reduced expression but failed to phenocopy the abolishment of stripes caused by removing Zld in trans. We show that Zld also indirectly regulates the eve stripes by establishing specific gap gene expression boundaries, which provides the embryonic spacing required for proper stripe activation.

Functional analysis of the Drosophila eve locus in response to non-canonical combinations of gap gene expression levels

Transcription factor combinations play a key role in shaping cellular identity. However, the precise relationship between specific combinations and downstream effects remains elusive. Here, we investigate this relationship within the context of the Drosophila eve locus, which is controlled by gap genes. We measure spatiotemporal levels of four gap genes in heterozygous and homozygous gap mutant embryos and correlate them with the striped eve activity pattern. Although changes in gap gene expression extend beyond the manipulated gene, the spatial patterns of Eve expression closely mirror canonical activation levels in wild type. Interestingly, some combinations deviate from the wild-type repertoire but still drive eve activation. Although in homozygous mutants some Eve stripes exhibit partial penetrance, stripes consistently emerge at reproducible positions, even with varying gap gene levels. Our findings suggest a robust molecular canalization of cell fates in gap mutants and provide insights into the regulatory constraints governing multi-enhancer gene loci.

Global repression by tailless during segmentation

The orphan nuclear receptor Tailless (Tll) exhibits conserved roles in brain formation and maintenance that are shared, for example, with vertebrate orthologous forms (Tlx). However, the early expression of tll in two gap domains in the segmentation cascade of Drosophila is unusual even for most other insects. Here we investigate tll regulation on pair-rule stripes. With ectopic misexpression of tll we detected unexpected repression of almost all pair-rule stripes of hairy (h), even-skipped (eve), runt (run), and fushi-tarazu (ftz). Examining Tll embryonic ChIP-chip data with regions mapped as Cis-Regulatory Modules (CRMs) of pair-rule stripes we verified Tll interactions to these regions. With the ChIP-chip data we also verified Tll interactions to the CRMs of gap domains and in the misexpression assay, Tll-mediated repression on Kruppel (Kr), kni (kni) and giant (gt) according to their differential sensitivity to Tll. These results with gap genes confirmed previous data from the literature and argue against indirect repression roles of Tll in the striped pattern. Moreover, the prediction of Tll binding sites in the CRMs of eve stripes and the mathematical modeling of their removal using an experimentally validated theoretical framework shows effects on eve stripes compatible with the absence of a repressor binding to the CRMs. In addition, modeling increased tll levels in the embryo results in the differential repression of eve stripes, agreeing well with the results of the misexpression assay. In genetic assays we investigated eve 5, that is strongly repressed by the ectopic domain and representative of more central stripes not previously implied to be under direct regulation of tll. While this stripe is little affected in tll-, its posterior border is expanded in gt- but detected with even greater expansion in gt-;tll-. We end up by discussing tll with key roles in combinatorial repression mechanisms to contain the expression of medial patterns of the segmentation cascade in the extremities of the embryo.

Genotypic identification and antimicrobial susceptibility profiles of coagulase-negative staphylococci isolated from sheep’s milk

Mastitis affect milk production and dairy products in sheep farming. Its main etiological agents are coagulase-negative staphylococci (CoNS). Identification based on DNA sequencing is considered the gold-standard method for these organisms. The present study identify CoNS species from 32 sheep’s milk isolates through polymerase chain reaction (PCR) with restriction fragment length polymorphism (RFLP) of the gap gene and to ascertain their susceptibility to nine selected antimicrobial agents by the agar disk diffusion technique. Four PCR-RFLP profiles were obtained, corresponding to three CoNS species: S. epidermidis (59.4%), S. chromogenes (34.4%), and S. devriesei (6.2%). All isolates were susceptible to cefalotin, ceftiofur, and penicillin/novobiocin. Twenty-one isolates (65.6%) were resistant to at least one antimicrobial agent, and 11 isolates (34.4%) were susceptible to all tested agents. Four isolates (12.5%) were multidrug-resistant, including an S. epidermidis isolate which exhibited resistance to oxacillin, confirmed by amplification of the mecA gene. In addition to constituting an animal health problem and causing economic losses to sheep farmers, CoNS have zoonotic potential

In-silico identification and comparison of transcription factor binding sites cluster in anterior-posterior patterning genes in Drosophila melanogaster and Tribolium castaneum.

The cis-regulatory data that help in transcriptional regulation is arranged into modular pieces of a few hundred base pairs called CRMs (cis-regulatory modules) and numerous binding sites for multiple transcription factors are prominent characteristics of these cis-regulatory modules. The present study was designed to localize transcription factor binding site (TFBS) clusters on twelve Anterior-posterior (A-P) genes in Tribolium castaneum and compare them to their orthologous gene enhancers in Drosophila melanogaster. Out of the twelve A-P patterning genes, six were gap genes (Kruppel, Knirps, Tailless, Hunchback, Giant, and Caudal) and six were pair rule genes (Hairy, Runt, Even-skipped, Fushi-tarazu, Paired, and Odd-skipped). The genes along with 20 kb upstream and downstream regions were scanned for TFBS clusters using the Motif Cluster Alignment Search Tool (MCAST), a bioinformatics tool that looks for set of nucleotide sequences for statistically significant clusters of non-overlapping occurrence of a given set of motifs. The motifs used in the current study were Hunchback, Caudal, Giant, Kruppel, Knirps, and Even-skipped. The results of the MCAST analysis revealed the maximum number of TFBS for Hunchback, Knirps, Caudal, and Kruppel in both D. melanogaster and T. castaneum, while Bicoid TFBS clusters were found only in D. melanogaster. The size of all the predicted TFBS clusters was less than 1kb in both insect species. These sequences revealed more transversional sites (Tv) than transitional sites (Ti) and the average Ti/Tv ratio was 0.75.

Shaping the scaling characteristics of gap gene expression patterns in Drosophila

How patterns are formed to scale with tissue size remains an unresolved problem. Here we investigate embryonic patterns of gap gene expression along the anterior-posterior (AP) axis in Drosophila. We use embryos that greatly differ in length and, importantly, possess distinct length-scaling characteristics of the Bicoid (Bcd) gradient. We systematically analyze the dynamic movements of gap gene expression boundaries in relation to both embryo length and Bcd input as a function of time. We document the process through which such dynamic movements drive both an emergence of a global scaling landscape and evolution of boundary-specific scaling characteristics. We show that, despite initial differences in pattern scaling characteristics that mimic those of Bcd in the anterior, such characteristics of final patterns converge. Our study thus partitions the contributions of Bcd input and regulatory dynamics inherent to the AP patterning network in shaping embryonic pattern's scaling characteristics.

A timer gene network is spatially regulated by the terminal system in the Drosophila embryo.

In insect embryos, anteroposterior patterning is coordinated by the sequential expression of the 'timer' genes caudal, Dichaete, and odd-paired, whose expression dynamics correlate with the mode of segmentation. In Drosophila, the timer genes are expressed broadly across much of the blastoderm, which segments simultaneously, but their expression is delayed in a small 'tail' region, just anterior to the hindgut, which segments during germband extension. Specification of the tail and the hindgut depends on the terminal gap gene tailless, but beyond this the regulation of the timer genes is poorly understood. We used a combination of multiplexed imaging, mutant analysis, and gene network modelling to resolve the regulation of the timer genes, identifying 11 new regulatory interactions and clarifying the mechanism of posterior terminal patterning. We propose that a dynamic Tailless expression gradient modulates the intrinsic dynamics of a timer gene cross-regulatory module, delineating the tail region and delaying its developmental maturation.

Kruppel-like Factors in Skeletal Physiology and Pathologies

Kruppel-like factors (KLFs) belong to a large group of zinc finger-containing transcription factors with amino acid sequences resembling the Drosophila gap gene Krüppel. Since the first report of molecular cloning of the KLF family gene, the number of KLFs has increased rapidly. Currently, 17 murine and human KLFs are known to play crucial roles in the regulation of transcription, cell proliferation, cellular differentiation, stem cell maintenance, and tissue and organ pathogenesis. Recent evidence has shown that many KLF family molecules affect skeletal cells and regulate their differentiation and function. This review summarizes the current understanding of the unique roles of each KLF in skeletal cells during normal development and skeletal pathologies.

An additive repression mechanism sets the anterior limits of anterior pair-rule stripes 1.

Segments are repeated anatomical units forming the body of insects. In Drosophila, the specification of the body takes place during the blastoderm through the segmentation cascade. Pair-rule genes such as hairy (h), even-skipped (eve), runt (run), and fushi-tarazu (ftz) are of the intermediate level of the cascade and each pair-rule gene is expressed in seven transversal stripes along the antero-posterior axis of the embryo. Stripes are formed by independent cis-regulatory modules (CRMs) under the regulation of transcription factors of maternal source and of gap proteins of the first level of the cascade. The initial blastoderm of Drosophila is a syncytium and it also coincides with the mid-blastula transition when thousands of zygotic genes are transcribed and their products are able to diffuse in the cytoplasm. Thus, we anticipated a complex regulation of the CRMs of the pair-rule stripes. The CRMs of h 1, eve 1, run 1, ftz 1 are able to be activated by bicoid (bcd) throughout the anterior blastoderm and several lines of evidence indicate that they are repressed by the anterior gap genes slp1 (sloppy-paired 1), tll (tailless) and hkb (huckebein). The modest activity of these repressors led to the premise of a combinatorial mechanism regulating the expression of the CRMs of h 1, eve 1, run 1, ftz 1 in more anterior regions of the embryo. We tested this possibility by progressively removing the repression activities of slp1, tll and hkb. In doing so, we were able to expose a mechanism of additive repression limiting the anterior borders of stripes 1. Stripes 1 respond depending on their distance from the anterior end and repressors operating at different levels.

Identification, Superantigen Toxin Gene Profile and Antimicrobial Resistance of Staphylococci Isolated from Polish Primitive Sheep Breeds.

Simple SummaryCurrent knowledge regarding the occurrence, differentiation, and pathogenicity of staphylococci isolated from animals is mainly based on clinical isolates, but little is known about staphylococci colonizing healthy animals. Even less is known about staphylococci that colonize sheep. While considering the microbiota inhabiting the primitive native breeds of sheep from the ecological/agritourism farm, there is no such information at all. The current study aimed to identify staphylococcal species isolated from healthy sheep of two native breeds (Świniarka and Wrzosówka) kept on the same ecological/agritourism farm. Furthermore, the prevalence of selected toxin genes and antimicrobial resistance was also determined. A total of 127 coagulase-negative staphylococci were identified from 61 sheep (33 of Świniarka) and (28 of Wrzosówka). From 1 to 4 staphylococcal isolates were obtained from one sheep. Staphylococcus aureus was not identified in any of the samples. In total, seven different species were identified. Over 60% of staphylococci were resistant to at least one antimicrobial of therapeutic importance and over 77% possessed from 1 to 7 genes encoding different toxins. The findings of this study prove that toxigenic and antimicrobial-resistant coagulase-negative staphylococci can colonize the nasal cavity of healthy sheep of primitive, native breeds.The study aimed to analyze staphylococcal microbiota of the nasal cavity of the primitive sheep breeds Polish Świniarka and Wrzosówka kept on the same ecological farm. The research included the identification of staphylococcal species, evaluation of the prevalence of genes encoding enterotoxins, staphylococcal enterotoxin-like proteins, exfoliative toxins, toxic shock syndrome toxin 1, and detection of antimicrobial resistance. From 61 swab samples gathered from Świniarka (33) and Wrzosówka (28) healthy sheep, 127 coagulase-negative staphylococci (CoNS) were isolated. Based on PCR-RFLP analysis of the gap gene using AluI and HpyCH4V enzymes, the isolates were identified as: Staphylococcus xylosus (33.9%), S. equorum (29.1%), S. arlettae (15%), S. warneri (9.4%), S. lentus (7.9%), S. succinus (3.9%) and S. sciuri (0.8%). Three of these species, S. lentus, S. succinus, and S. sciuri, were detected only from the Świniarka breed. It was found that 77.2% of isolates harbored from 1 to 7 out of 21 analyzed genes for superantigenic toxins. The greatest diversity of toxin genes was recorded for S. equorum (16 different genes). The most prevalent gene was ser (40.2%). The incidence and number of resistances to antimicrobials were found to be bacterial species but not sheep breed dependent. The highest percentage of resistance was found for S. sciuri. The most frequent resistance was observed to clindamycin (45.7%). The findings of this study prove that toxigenic and antimicrobial resistant CoNS can colonize the nasal cavity of healthy sheep.

Latent space of a small genetic network: Geometry of dynamics and information

The high-dimensional character of most biological systems presents genuine challenges for modeling and prediction. Here we propose a neural network-based approach for dimensionality reduction and analysis of biological gene expression data, using, as a case study, a well-known genetic network in the early Drosophila embryo, the gap gene patterning system. We build an autoencoder compressing the dynamics of spatial gap gene expression into a two-dimensional (2D) latent map. The resulting 2D dynamics suggests an almost linear model, with a small bare set of essential interactions. Maternally defined spatial modes control gap genes positioning, without the classically assumed intricate set of repressive gap gene interactions. This, surprisingly, predicts minimal changes of neighboring gap domains when knocking out gap genes, consistent with previous observations. Latent space geometries in maternal mutants are also consistent with the existence of such spatial modes. Finally, we show how positional information is well defined and interpretable as a polar angle in latent space. Our work illustrates how optimization of small neural networks on medium-sized biological datasets is sufficiently informative to capture essential underlying mechanisms of network function.

Scaling dictates the decoder structure

Despite fluctuations in embryo size within a species, the spatial gene expression pattern and hence the embryonic structure can nonetheless maintain the correct proportion to the embryo size. This is known as the scaling phenomenon. For morphogen-induced patterning of gene expression, the positional information encoded in the local morphogen concentrations is decoded by the downstream genetic network (the decoder). In this paper, we show that the requirement of scaling sets severe constraints on the geometric structure of such a local decoder, which in turn enables deduction of mutants’ behavior and extraction of regulation information without going into any molecular details. We demonstrate that the Drosophila gap gene system achieves scaling in the way consistent with our theory—the decoder geometry required by scaling correctly accounts for the observed gap gene expression pattern in nearly all maternal morphogen mutants. Furthermore, the regulation logic and the coding/decoding strategy of the gap gene system can also be revealed from the decoder geometry. Our work provides a general theoretical framework for a large class of problems where scaling output is achieved by non-scaling inputs and a local decoder, as well as a unified understanding of scaling, mutants’ behavior, and gene regulation for the Drosophila gap gene system.