Early Drosophila Development Research Articles

Global coupling of R-loopdynamics with RNA polymerase II modulates gene expression and early development of Drosophila.

R-loops are involved in many biological processes in cells, yet the regulatory principles for R-loops in vivo and their impact on development remain to be explored. Here, we modified the CUT&Tag strategy to profile R-loops in Drosophila at multiple developmental stages. While high GC content promotes R-loop formation in mammalian cells, it is not required in Drosophila. In contrast, RNAPII abundance appears to be a universal inducing factor for R-loop formation, including active promoters and enhancers, and H3K27me3 decorated repressive regions and intergenic repeat sequences. Importantly, such a regulatory relationship is dynamically maintained throughout development, and development-related transcription factors may regulate RNAPII activation and R-loop dynamics. By ablating Spt6, we further showed the global R-loop induction coupled with RNAPII pausing. Importantly, depending on the gene length, genes underwent up- or down-regulation, both of which were largely reversed by rnh1 overexpression, suggesting that R-loops play a significant role in the divergent regulation of transcription by Spt6 ablation. DNA damage, defects in survival, and cuticle development were similarly alleviated by rnh1 overexpression. Altogether, our findings indicate that dynamic R-loop regulation is dictated by RNAPII pausing and transcription activity, and plays a feedback role in gene regulation, genome stability maintenance, and Drosophila development.

Affinity hierarchies and amphiphilic proteins underlie the co-assembly of nucleolar and heterochromatin condensates.

Nucleoli are surrounded by Pericentromeric Heterochromatin (PCH), reflecting a close spatial association between the two largest biomolecular condensates in eukaryotic nuclei. Nucleoli are the sites of ribosome synthesis, while the repeat-rich PCH is essential for chromosome segregation, genome stability, and transcriptional silencing. How and why these two distinct condensates co-assemble is unclear. Here, using high-resolution live imaging of Drosophila embryogenesis, we find that de novo establishment of PCH around the nucleolus is highly dynamic, transitioning from the nuclear edge to surrounding the nucleolus. Eliminating the nucleolus by removing the ribosomal RNA genes (rDNA) resulted in increased PCH compaction and subsequent reorganization into a toroidal structure. In addition, in embryos lacking rDNA, some nucleolar proteins were redistributed into new bodies or 'neocondensates', including enrichment in the PCH toroidal hole. Combining these observations with physical modeling revealed that nucleolar-PCH associations can be mediated by a hierarchy of interaction strengths between PCH, nucleoli, and 'amphiphilic' protein(s) that have affinities for both nucleolar and PCH components. We validated this model by identifying a candidate amphiphile, a DEAD-Box RNA Helicase called Pitchoune, whose depletion or mutation of its PCH interaction motif disrupted PCH-nucleolar associations. Together, this study unveils a dynamic program for establishing nucleolar-PCH associations during animal development, demonstrates that nucleoli are required for normal PCH organization, and identifies Pitchoune as an amphiphilic molecular link required for PCH-nucleolar associations.

OVO positively regulates essential maternal pathways by binding near the transcriptional start sites in the Drosophila female germline

Differentiation of female germline stem cells into a mature oocyte includes the expression of RNAs and proteins that drive early embryonic development in Drosophila. We have little insight into what activates the expression of these maternal factors. One candidate is the zinc-finger protein OVO. OVO is required for female germline viability and has been shown to positively regulate its own expression, as well as a downstream target, ovarian tumor, by binding to the transcriptional start site (TSS). To find additional OVO targets in the female germline and further elucidate OVO’s role in oocyte development, we performed ChIP-seq to determine genome-wide OVO occupancy, as well as RNA-seq comparing hypomorphic and wild type rescue ovo alleles. OVO preferentially binds in close proximity to target TSSs genome-wide, is associated with open chromatin, transcriptionally active histone marks, and OVO-dependent expression. Motif enrichment analysis on OVO ChIP peaks identified a 5’-TAACNGT-3’ OVO DNA binding motif spatially enriched near TSSs. However, the OVO DNA binding motif does not exhibit precise motif spacing relative to the TSS characteristic of RNA polymerase II complex binding core promoter elements. Integrated genomics analysis showed that 525 genes that are bound and increase in expression downstream of OVO are known to be essential maternally expressed genes. These include genes involved in anterior/posterior/germ plasm specification (bcd, exu, swa, osk, nos, aub, pgc, gcl), egg activation (png, plu, gnu, wisp, C(3)g, mtrm), translational regulation (cup, orb, bru1, me31B), and vitelline membrane formation (fs(1)N, fs(1)M3, clos). This suggests that OVO is a master transcriptional regulator of oocyte development and is responsible for the expression of structural components of the egg as well as maternally provided RNAs that are required for early embryonic development.

MIGGRI: A multi-instance graph neural network model for inferring gene regulatory networks for Drosophila from spatial expression images.

Recent breakthrough in spatial transcriptomics has brought great opportunities for exploring gene regulatory networks (GRNs) from a brand-new perspective. Especially, the local expression patterns and spatio-temporal regulation mechanisms captured by spatial expression images allow more delicate delineation of the interplay between transcript factors and their target genes. However, the complexity and size of spatial image collections pose significant challenges to GRN inference using image-based methods. Extracting regulatory information from expression images is difficult due to the lack of supervision and the multi-instance nature of the problem, where a gene often corresponds to multiple images captured from different views. While graph models, particularly graph neural networks, have emerged as a promising method for leveraging underlying structure information from known GRNs, incorporating expression images into graphs is not straightforward. To address these challenges, we propose a two-stage approach, MIGGRI, for capturing comprehensive regulatory patterns from image collections for each gene and known interactions. Our approach involves a multi-instance graph neural network (GNN) model for GRN inference, which first extracts gene regulatory features from spatial expression images via contrastive learning, and then feeds them to a multi-instance GNN for semi-supervised learning. We apply our approach to a large set of Drosophila embryonic spatial gene expression images. MIGGRI achieves outstanding performance in the inference of GRNs for early eye development and mesoderm development of Drosophila, and shows robustness in the scenarios of missing image information. Additionally, we perform interpretable analysis on image reconstruction and functional subgraphs that may reveal potential pathways or coordinate regulations. By leveraging the power of graph neural networks and the information contained in spatial expression images, our approach has the potential to advance our understanding of gene regulation in complex biological systems.

An improved Erk biosensor detects oscillatory Erk dynamics driven by mitotic erasure during early development

Extracellular signal-regulated kinase (Erk) signaling dynamics elicit distinct cellular responses in a variety of contexts. The early zebrafish embryo is an ideal model to explore the role of Erk signaling dynamics invivo, as a gradient of activated diphosphorylated Erk (P-Erk) is induced by fibroblast growth factor (Fgf) signaling at the blastula margin. Here, we describe an improved Erk-specific biosensor, which we term modified Erk kinase translocation reporter (modErk-KTR). We demonstrate the utility of this biosensor invitro and in developing zebrafish and Drosophila embryos. Moreover, we show that Fgf/Erk signaling is dynamic and coupled to tissue growth during both early zebrafish and Drosophila development. Erk activity is rapidly extinguished just prior to mitosis, which we refer to as mitotic erasure, inducing periods of inactivity, thus providing a source of heterogeneity in an asynchronously dividing tissue. Our modified reporter and transgenic lines represent an important resource for interrogating the role of Erk signaling dynamics invivo.

Localization of the Drosophila pioneer factor GAF to subnuclear foci is driven by DNA binding and required to silence satellite repeat expression

The eukaryotic genome is organized to enable the precise regulation of gene expression. This organization is established as the embryo transitions from a fertilized gamete to a totipotent zygote. To understand the factors and processes that drive genomic organization, we focused on the pioneer factor GAGA factor (GAF) that is required for early development in Drosophila. GAF transcriptionally activates the zygotic genome and is localized to subnuclear foci. This non-uniform distribution is driven by binding to highly abundant GA repeats. At GA repeats, GAF is necessary to form heterochromatin and silence transcription. Thus, GAF is required to establish both active and silent regions. We propose that foci formation enables GAF to have opposing transcriptional roles within a single nucleus. Our data support a model in which the subnuclear concentration of transcription factors acts to organize the nucleus into functionally distinct domains essential for the robust regulation of gene expression.

Studying primary pair-rule gene eve in the early development of Drosophila by using Zelda as approach.

This paper mainly studies the primary pair-rule gene Even-skipped to reveal its position in the early transcription network established in the Drosophila embryo and gain a deeper view of the network. We used the genome activator Zelda as an approach. We found Zelda’s role in regulating Eve by committing CHIP- seq towards RNA polymerase and Zelda protein in wt and Zld- embryo. With the PWM matrix, we know transcription factors that might bind to eve. By further researching its TFs, we find the relationship between timing and pattern of expression of these genes involved in the early transcriptional network, while Zld is still our approach. Also, by looking into the specific gene, we find the possible timing mechanism on Eve based on the accumulation of earlier transcribed transcription factors. The role of maternally loaded transcription factors in early development is once again shown by affecting the timing and patterning of primary patterning transcription factors.

Determination of Complex Formation between Drosophila Nrf2 and GATA4 Factors at Selective Chromatin Loci Demonstrates Transcription Coactivation.

Nrf2 is the dominant cellular stress response factor that protects cells through transcriptional responses to xenobiotic and oxidative stimuli. Nrf2 malfunction is highly correlated with many human diseases, but the underlying molecular mechanisms remain to be fully uncovered. GATA4 is a conserved GATA family transcription factor that is essential for cardiac and dorsal epidermal development. Here, we describe a novel interaction between Drosophila Nrf2 and GATA4 proteins, i.e., cap'n'collar C (CncC) and Pannier (Pnr), respectively. Using the bimolecular fluorescence complementation (BiFC) assay-a unique imaging tool for probing protein complexes in living cells-we detected CncC-Pnr complexes in the nuclei of Drosophila embryonic and salivary gland cells. Visualization of CncC-Pnr BiFC signals on the polytene chromosome revealed that CncC and Pnr tend to form complexes in euchromatic regions, with a preference for loci that are not highly occupied by CncC or Pnr alone. Most genes within these loci are activated by the CncC-Pnr BiFC, but not by individually expressed CncC or Pnr fusion proteins, indicating a novel mechanism whereby CncC and Pnr interact at specific genomic loci and coactivate genes at these loci. Finally, CncC-induced early lethality can be rescued by Pnr depletion, suggesting that CncC and Pnr function in the same genetic pathway during the early development of Drosophila. Taken together, these results elucidate a novel crosstalk between the Nrf2 xenobiotic/oxidative response factor and GATA factors in the transcriptional regulation of development. This study also demonstrates that the polytene chromosome BiFC assay is a valuable tool for mapping genes that are targeted by specific transcription factor complexes.

Adaptable P body physical states differentially regulate bicoid mRNA storage during early Drosophila development.

SummaryRibonucleoprotein condensates can exhibit diverse physical states in vitro and in vivo. Despite considerable progress, the relevance of condensate physical states for in vivo biological function remains limited. Here, we investigated the physical properties of processing bodies (P bodies) and their impact on mRNA storage in mature Drosophila oocytes. We show that the conserved DEAD-box RNA helicase Me31B forms viscous P body condensates, which adopt an arrested physical state. We demonstrate that structurally distinct proteins and protein-protein interactions, together with RNA, regulate the physical properties of P bodies. Using live imaging and in situ hybridization, we show that the arrested state and integrity of P bodies support the storage of bicoid (bcd) mRNA and that egg activation modulates P body properties, leading to the release of bcd for translation in the early embryo. Together, this work provides an example of how physical states of condensates regulate cellular function in development.

The mRNA decapping protein 2 (DCP2) is a major regulator of developmental events in Drosophila-insights from expression paradigms.

The Drosophila genome codes for two decapping proteins, DCP1 and DCP2, out of which DCP2 is the active decapping enzyme. The present endeavour explores the endogenous promoter firing, transcript and protein expression of DCP2 in Drosophila wherein, besides a ubiquitous expression across development, we identify an active expression paradigm during dorsal closure and a plausible moonlighting expression in the Corazonin neurons of the larval brain. We also demonstrate that the ablation of DCP2 leads to embryonic lethality and defects in vital morphogenetic processes whereas a knockdown of DCP2 in the Corazonin neurons reduces the sensitivity to ethanol in adults, thereby ascribing novel regulatory roles to DCP2. Our findings unravel novel putative roles for DCP2 and identify it as a candidate for studies on the regulated interplay of essential molecules during early development in Drosophila, nay the living world.

HP1 drives de novo 3D genome reorganization in early Drosophila embryos

Fundamental features of 3D genome organization are established de novo in the early embryo, including clustering of pericentromeric regions, the folding of chromosome arms and the segregation of chromosomes into active (A-) and inactive (B-) compartments. However, the molecular mechanisms that drive de novo organization remain unknown1,2. Here, by combining chromosome conformation capture (Hi-C), chromatin immunoprecipitation with high-throughput sequencing (ChIP–seq), 3D DNA fluorescence in situ hybridization (3D DNA FISH) and polymer simulations, we show that heterochromatin protein 1a (HP1a) is essential for de novo 3D genome organization during Drosophila early development. The binding of HP1a at pericentromeric heterochromatin is required to establish clustering of pericentromeric regions. Moreover, HP1a binding within chromosome arms is responsible for overall chromosome folding and has an important role in the formation of B-compartment regions. However, depletion of HP1a does not affect the A-compartment, which suggests that a different molecular mechanism segregates active chromosome regions. Our work identifies HP1a as an epigenetic regulator that is involved in establishing the global structure of the genome in the early embryo.

Cis-regulatory chromatin loops arise before TADs and gene activation, and are independent of cell fate during early Drosophila development.

Acquisition of cell fate is thought to rely on the specific interaction of remote cis-regulatory modules (CRMs), for example, enhancers and target promoters. However, the precise interplay between chromatin structure and gene expression is still unclear, particularly within multicellular developing organisms. In the present study, we employ Hi-M, a single-cell spatial genomics approach, to detect CRM-promoter looping interactions within topologically associating domains (TADs) during early Drosophila development. By comparing cis-regulatory loops in alternate cell types, we show that physical proximity does not necessarily instruct transcriptional states. Moreover, multi-way analyses reveal that multiple CRMs spatially coalesce to form hubs. Loops and CRM hubs are established early during development, before the emergence of TADs. Moreover, CRM hubs are formed, in part, via the action of the pioneer transcription factor Zelda and precede transcriptional activation. Our approach provides insight into the role of CRM-promoter interactions in defining transcriptional states, as well as distinct cell types.

Targeting of the Dosage-Compensated Male X-Chromosome during Early Drosophila Development

SUMMARYDosage compensation, which corrects for the imbalance in X-linked gene expression between XX females and XY males, represents a model for how genes are targeted for coordinated regulation. However, the mechanism by which dosage compensation complexes identify the X chromosome during early development remains unknown because of the difficulty of sexing embryos before zygotic transcription using X- or Y-linked reporter transgenes. We used meiotic drive to sex Drosophila embryos before zygotic transcription and ChIP-seq to measure the dynamics of dosage compensation factor targeting. The Drosophila male-specific lethal dosage compensation complex (MSLc) requires the ubiquitous zinc-finger protein chromatin-linked adaptor for MSL proteins (CLAMP) to identify the X chromosome. We observe a multi-stage process in which MSLc first identifies CLAMP binding sites throughout the genome, followed by concentration at the strongest X-linked MSLc sites. We provide insight into the dynamics of binding site recognition by a large transcription complex during early development.

Meet the First Authors.

HomeCirculation ResearchVol. 125, No. 10Meet the First Authors Free AccessIn BriefPDF/EPUBAboutView PDFView EPUBSections ToolsAdd to favoritesDownload citationsTrack citationsPermissions ShareShare onFacebookTwitterLinked InMendeleyReddit Jump toFree AccessIn BriefPDF/EPUBMeet the First Authors Originally published24 Oct 2019https://doi.org/10.1161/RES.0000000000000306Circulation Research. 2019;125:870is related toBlood Flow Suppresses Vascular Anomalies in a Zebrafish Model of Cerebral Cavernous MalformationsDiversification and CXCR4-Dependent Establishment of the Bone Marrow B-1a Cell Pool Governs Atheroprotective IgM Production Linked to Human Coronary AtherosclerosisPeripheral Vasoconstriction During Mental Stress and Adverse Cardiovascular Outcomes in Patients With Coronary Artery DiseaseBlood Flow Suppresses CCM Phenotypes in Zebrafish (p e43)Download figureDownload PowerPointDr Claudia J. Rödel is a molecular biologist interested in cardiovascular development and disease. She earned her undergraduate degree in Biology from the Georg-August University of Göttingen, as well as her PhD. As a graduate student, she investigated the evolution of post-transcriptional control during early Drosophila development in the lab of Michalis Averof (mentor). In the lab of Dimitris Beis at the Biomedical Research Foundation Academy of Athens, she discovered the beauty of the zebrafish vertebrate model, and switched to the field of cardiovascular research. She then joined the lab of Salim Seyfried at the University of Potsdam, where she began her work on cerebral cavernous malformation (CCM), an enigmatic disease with vascular anomalies that affect almost exclusively the cerebral vasculature. Combining zebrafish genetics and cell biology, Claudia and her colleagues were able to address the impact of biomechanical stimuli that affect CCM angiogenesis, revealing that blood flow takes on a protective role in CCM. For her, the zebrafish system will continue to provide essential and versatile tools in the quest for a better understanding of the complex mechanisms that contribute to this rare but devastating disease. Outside the lab, the science never stops but she also enjoys spending time with her family and friends.Vasoconstriction During Mental Stress and Outcomes (p 874)Download figureDownload PowerPointDr Jeong Hwan Kim is a cardiology fellow at Emory University School of Medicine. He earned his BA in Biology and Economics from Swarthmore College, and he obtained his MD from the University of Chicago Pritzker School of Medicine in 2014. He subsequently completed his internal medicine residency at Emory University, where he stayed for his cardiology fellowship under the mentorship of Dr Arshed Quyyumi. Dr Kim is interested in studying the mechanisms by which psychological stress affects patients with cardiovascular disease. He also has interests in biomarkers, vascular biology, risk prediction, and social determinants of cardiovascular disease, and plans to pursue a career as a clinical investigator in cardiology. He can be found on twitter, @JeongHwanKimMD.CXCR4 Sustains Atheroprotective IgM in Bone Marrow (p e55)Download figureDownload PowerPointAditi Upadhye is a PhD candidate at the University of Virginia in the laboratory of Dr Coleen McNamara. She earned her BS in Biochemistry from Lyman Briggs College at Michigan State University, where her coursework and exposure to translational research led her to pursue graduate studies in biomedical sciences. Aditi’s thesis project in Dr McNamara’s lab focuses on B-1 cell-mediated protection against atherosclerosis. In the future, she plans to continue pursuing research on immune mechanisms present in chronic inflammation and opportunities to translate scientific discoveries to human health and disease. Outside of the laboratory, Aditi enjoys staying active through yoga, running, and hiking, and trying new baking recipes to make up for it. Previous Back to top Next FiguresReferencesRelatedDetailsRelated articlesBlood Flow Suppresses Vascular Anomalies in a Zebrafish Model of Cerebral Cavernous MalformationsClaudia Jasmin Rödel, et al. Circulation Research. 2019;125:e43-e54Diversification and CXCR4-Dependent Establishment of the Bone Marrow B-1a Cell Pool Governs Atheroprotective IgM Production Linked to Human Coronary AtherosclerosisAditi Upadhye, et al. Circulation Research. 2019;125:e55-e70Peripheral Vasoconstriction During Mental Stress and Adverse Cardiovascular Outcomes in Patients With Coronary Artery DiseaseJeong Hwan Kim, et al. Circulation Research. 2019;125:874-883 October 25, 2019Vol 125, Issue 10 Advertisement Article InformationMetrics © 2019 American Heart Association, Inc.https://doi.org/10.1161/RES.0000000000000306PMID: 31647769 Originally publishedOctober 24, 2019 PDF download Advertisement

Nature and Functions of Telomeric Transcripts.

Telomeres are complex and dynamic structures whose functions and composition change during the cell cycle and development. Telomeric transcripts are essential components of telomeres. Transcription regulation and cellular levels of telomeric RNAs are closely associated with the control of telomere length, formation of telomeric chromatin, telomere replication, and regulation of non-telomeric gene transcription, which indicates a critical regulatory role of telomeric RNAs in telomere protection and transmission of signals about the state of telomeres to cellular genes. The studies of telomeric transcriptome in early Drosophila development have revealed a new level of genomic stability regulation involving telomeric RNAs. Due to their ability to interact with multiple proteins and to translocate in the cell, telomeric transcripts are important participants of telomeric signaling pathways, whose mechanisms are still to be understood at the organism level.

Drosophila melanogaster research: history, breakthrough and perspectives

Abstract The common fruit fly, or Drosophila Melanogaster, has been used as an object of biomedicals studies for over a century. It has been mostly employed in genetic research, as it exhibits several advantages which make its use relatively easy and cheap, with the results widely translatable into further vertebrate studies. This model been the basis of the work of Christiane Nusslein-Volhard, who together with Eric Wieschaus unravelled much of the mystery surrounding early drosophila development in the 1970s-1980s, laying foundations for broader understanding of multicellular organism embryogenesis, which brought them a Nobel prize in Physiology and Medicine in 1995. The knowledge gained from drosophila studies improves the basic understanding of developmental processes, while the model itself is relatively easy to maintain, analyse and translate the results onto other species. While models such as Zebrafish present better with other vertebrates, drosophila remains a very important element of genetic research, finding even more applications with the development of current science and medicine. Hence, in this short review, the outline of the history, breakthroughs and perspectives of the drosophila research has been presented.

Robustness of expression pattern formation due to dynamic equilibrium in gap gene system of an early Drosophila embryo

The first manifestation of a segmentation pattern in the early Drosophila development is the formation of expression domains of genes belonging to the gap class. In our previous research the phenomenon of the gap system’s robustness, exhibited as the ability to reduce highly variable gene expression in the course of development, was explained as a result of gene cross-regulation. In this paper we formulate the rigorous robustness conditions using the inherent properties of gap gene family.We propose an approach based on the observation that for the formation of a pattern with well-established domain borders it is necessary that there exist a stationary nucleus in which the gene expression level is almost constant in time. The dynamics of gap gene expression is described by a gene circuit model that correctly reproduces the observed principles of the border positioning.We take an advantage of the fact that the system behavior in a stationary nucleus is described by an algebraic equation and hence can be easily handled analytically. This enables us to explicitly characterize the gene cross-regulation properties guaranteeing the system robustness through the spatial behavior of the patterning system along the main embryo axis. In particular, it is proved that if the total regulatory action of all the genes acting in the border formation area changes synchronously with the initial gene gradient, the system will filtrate the initial positioning error and action of highly variable external input factors.We will now show how this approach can be applied as an instrument for the analysis of the robustness mechanism and revealing new biological findings.

LiveFly: A Toolbox for the Analysis of Transcription Dynamics in Live Drosophila Embryos.

We present the LiveFly toolbox for quantitative analysis of transcription dynamics in live Drosophila embryos. The toolbox allows users to process two-color 3D confocal movies acquired using nuclei-labeling and the fluorescent RNA-tagging system described in the previous chapter and export the nuclei's position as a function of time, their lineages and the intensity traces of the active loci. The toolbox, which is tailored for the context of Drosophila early development, is semiautomatic, and requires minimal user intervention. It also includes a tool to combine data from multiple movies and visualize several features of the intensity traces and the expression pattern.

Conserved noncoding transcription and core promoter regulatory code in early Drosophila development.

Multicellular development is driven by regulatory programs that orchestrate the transcription of protein-coding and noncoding genes. To decipher this genomic regulatory code, and to investigate the developmental relevance of noncoding transcription, we compared genome-wide promoter activity throughout embryogenesis in 5 Drosophila species. Core promoters, generally not thought to play a significant regulatory role, in fact impart restrictions on the developmental timing of gene expression on a global scale. We propose a hierarchical regulatory model in which core promoters define broad windows of opportunity for expression, by defining a range of transcription factors from which they can receive regulatory inputs. This two-tiered mechanism globally orchestrates developmental gene expression, including extremely widespread noncoding transcription. The sequence and expression specificity of noncoding RNA promoters are evolutionarily conserved, implying biological relevance. Overall, this work introduces a hierarchical model for developmental gene regulation, and reveals a major role for noncoding transcription in animal development.

Translating natural genetic variation to gene expression in a computational model of the Drosophila gap gene regulatory network.

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.