Drosophila immune priming to Enterococcus faecalis relies on immune tolerance rather than resistance.

Innate immune priming increases an organism's survival of a second infection after an initial, non-lethal infection. We used Drosophila melanogaster and an insect-derived strain of Enterococcus faecalis to study transcriptional control of priming. In contrast to other pathogens, the enhanced survival in primed animals does not correlate with decreased E. faecalis load. Further analysis shows that primed organisms tolerate, rather than resist infection. Using RNA-seq of immune tissues, we found many genes were upregulated in only primed flies, suggesting a distinct transcriptional program in response to initial and secondary infections. In contrast, few genes continuously express throughout the experiment or more efficiently re-activate upon reinfection. Priming experiments in immune deficient mutants revealed Imd is largely dispensable for responding to a single infection but needed to fully prime. Together, this indicates the fly's innate immune response is plastic-differing in immune strategy, transcriptional program, and pathway use depending on infection history.

Shadow enhancers mediate trade-offs between transcriptional noise and fidelity.

Enhancers are stretches of regulatory DNA that bind transcription factors (TFs) and regulate the expression of a target gene. Shadow enhancers are two or more enhancers that regulate the same target gene in space and time and are associated with most animal developmental genes. These multi-enhancer systems can drive more consistent transcription than single enhancer systems. Nevertheless, it remains unclear why shadow enhancer TF binding sites are distributed across multiple enhancers rather than within a single large enhancer. Here, we use a computational approach to study systems with varying numbers of TF binding sites and enhancers. We employ chemical reaction networks with stochastic dynamics to determine the trends in transcriptional noise and fidelity, two key performance objectives of enhancers. This reveals that while additive shadow enhancers do not differ in noise and fidelity from their single enhancer counterparts, sub- and superadditive shadow enhancers have noise and fidelity trade-offs not available to single enhancers. We also use our computational approach to compare the duplication and splitting of a single enhancer as mechanisms for the generation of shadow enhancers and find that the duplication of enhancers can decrease noise and increase fidelity, although at the metabolic cost of increased RNA production. A saturation mechanism for enhancer interactions similarly improves on both of these metrics. Taken together, this work highlights that shadow enhancer systems may exist for several reasons: genetic drift or the tuning of key functions of enhancers, including transcription fidelity, noise and output.

Longitudinal monitoring of individual infection progression in Drosophila melanogaster.

The innate immune system is critical for infection survival. Drosophila melanogaster is a key model for understanding the evolution and dynamics of innate immunity. Current toolsets for fly infection studies are limited in throughput and, because of their destructive nature, cannot generate longitudinal measurements in individual animals. We report a bioluminescent imaging strategy enabling non-invasive characterization of pathogen load. By using Escherichia coli expressing the ilux operon, we demonstrate that photon flux from autobioluminescent bacteria can be used to monitor pathogen loads in individual, living flies. Because animal sacrifice is not necessary to estimate pathogen load, stochastic responses to infection can be characterized in individuals over time. The high temporal resolution of bioluminescence imaging enables visualization of the dynamics of microbial clearance on the hours time-scale. This non-invasive imaging strategy provides a simple and scalable platform to observe changes in pathogen load in vivo over time.

Two promoters integrate multiple enhancer inputs to drive wild-type knirps expression in the Drosophila melanogaster embryo.

Proper development depends on precise spatiotemporal gene expression patterns. Most developmental genes are regulated by multiple enhancers and often by multiple core promoters that generate similar transcripts. We hypothesize that multiple promoters may be required either because enhancers prefer a specific promoter or because multiple promoters serve as a redundancy mechanism. To test these hypotheses, we studied the expression of the knirps locus in the early Drosophila melanogaster embryo, which is mediated by multiple enhancers and core promoters. We found that one of these promoters resembles a typical "sharp" developmental promoter, while the other resembles a "broad" promoter usually associated with housekeeping genes. Using synthetic reporter constructs, we found that some, but not all, enhancers in the locus show a preference for one promoter, indicating that promoters provide both redundancy and specificity. By analyzing the reporter dynamics, we identified specific burst properties during the transcription process, namely burst size and frequency, that are most strongly tuned by the combination of promoter and enhancer. Using locus-sized reporters, we discovered that enhancers with no promoter preference in a synthetic setting have a preference in the locus context. Our results suggest that the presence of multiple promoters in a locus is due both to enhancer preference and a need for redundancy and that "broad" promoters with dispersed transcription start sites are common among developmental genes. They also imply that it can be difficult to extrapolate expression measurements from synthetic reporters to the locus context, where other variables shape a gene's overall expression pattern.

Mechanistic gene networks inferred from single-cell data with an outlier-insensitive method.

With advances in single-cell techniques, measuring gene dynamics at cellular resolution has become practicable. In contrast, the increased complexity of data has made it more challenging computationally to unravel underlying biological mechanisms. Thus, it is critical to develop novel computational methods capable of dealing with such complexity and of providing predictive deductions from such data. Many methods have been developed to address such challenges, each with its own advantages and limitations. We present an iterative regression algorithm for inferring a mechanistic gene network from single-cell data, especially suited to overcoming problems posed by measurement outliers. Using this regression, we infer a developmental model for the gene dynamics in Drosophila melanogaster blastoderm embryo. Our results show that the predictive power of the inferred model is higher than that of other models inferred with least squares and ridge regressions. As a baseline for how well a mechanistic model should be expected to perform, we find that model predictions of the gene dynamics are more accurate than predictions made with neural networks of varying architectures and complexity. This holds true even in the limit of small sample sizes. We compare predictions for various gene knockouts with published experimental results, finding substantial qualitative agreement. We also make predictions for gene dynamics under various gene network perturbations, impossible in non-mechanistic models.

Molecular competition can shape enhancer activity in the Drosophila embryo.

Transgenic reporters allow the measurement of regulatory DNA activity in vivo and consequently have long been useful tools for studying enhancers. Despite their utility, few studies have investigated the effects these reporters may have on the expression of other genes. Understanding these effects is required to accurately interpret reporter data and characterize gene regulatory mechanisms. By measuring the expression of Kruppel (Kr) enhancer reporters in live Drosophila embryos, we find reporters inhibit one another's expression and that of a nearby endogenous gene. Using synthetic transcription factor (TF) binding site arrays, we present evidence that competition for TFs is partially responsible for the observed transcriptional inhibition. We develop a simple thermodynamic model that predicts competition of the measured magnitude specifically when TF binding is restricted to distinct nuclear subregions. Our findings underline an unexpected role of the non-homogenous nature of the nucleus in regulating gene expression.

The mode of expression divergence in Drosophila fat body is infection-specific.

Transcription is controlled by interactions of cis-acting DNA elements with diffusible trans-acting factors. Changes in cis or trans factors can drive expression divergence within and between species, and their relative prevalence can reveal the evolutionary history and pressures that drive expression variation. Previous work delineating the mode of expression divergence in animals has largely used whole-body expression measurements in one condition. Because cis-acting elements often drive expression in a subset of cell types or conditions, these measurements may not capture the complete contribution of cis-acting changes. Here, we quantify the mode of expression divergence in the Drosophila fat body, the primary immune organ, in several conditions, using two geographically distinct lines of D. melanogaster and their F1 hybrids. We measured expression in the absence of infection and in infections with Gram-negative S. marcescens or Gram-positive E. faecalis bacteria, which trigger the two primary signaling pathways in the Drosophila innate immune response. The mode of expression divergence strongly depends on the condition, with trans-acting effects dominating in response to Gram-negative infection and cis-acting effects dominating in Gram-positive and preinfection conditions. Expression divergence in several receptor proteins may underlie the infection-specific trans effects. Before infection, when the fat body has a metabolic role, there are many compensatory effects, changes in cis and trans that counteract each other to maintain expression levels. This work shows that within a single tissue, the mode of expression divergence varies between conditions and suggests that these differences reflect the diverse evolutionary histories of host-pathogen interactions.

Enhancer redundancy in development and disease.

Shadow enhancers are seemingly redundant transcriptional cis-regulatory elements that regulate the same gene and drive overlapping expression patterns. Recent studies have shown that shadow enhancers are remarkably abundant and control most developmental gene expression in both invertebrates and vertebrates, including mammals. Shadow enhancers might provide an important mechanism for buffering gene expression against mutations in non-coding regulatory regions of genes implicated in human disease. Technological advances in genome editing and live imaging have shed light on how shadow enhancers establish precise gene expression patterns and confer phenotypic robustness. Shadow enhancers can interact in complex ways and may also help to drive the formation of transcriptional hubs within the nucleus. Despite their apparent redundancy, the prevalence and evolutionary conservation of shadow enhancers underscore their key role in emerging metazoan gene regulatory networks.

A Mutation in the Drosophila melanogaster eve Stripe 2 Minimal Enhancer Is Buffered by Flanking Sequences.

Enhancers are DNA sequences composed of transcription factor binding sites that drive complex patterns of gene expression in space and time. Until recently, studying enhancers in their genomic context was technically challenging. Therefore, minimal enhancers, the shortest pieces of DNA that can drive an expression pattern that resembles a gene's endogenous pattern, are often used to study features of enhancer function. However, evidence suggests that some enhancers require sequences outside the minimal enhancer to maintain function under environmental perturbations. We hypothesized that these additional sequences also prevent misexpression caused by a transcription factor binding site mutation within a minimal enhancer. Using the Drosophila melanogastereven-skipped stripe 2 enhancer as a case study, we tested the effect of a Giant binding site mutation (gt-2) on the expression patterns driven by minimal and extended enhancer reporter constructs. We found that, in contrast to the misexpression caused by the gt-2 binding site deletion in the minimal enhancer, the same gt-2 binding site deletion in the extended enhancer did not have an effect on expression. The buffering of expression levels, but not expression pattern, is partially explained by an additional Giant binding site outside the minimal enhancer. Deleting the gt-2 binding site in the endogenous locus had no significant effect on stripe 2 expression. Our results indicate that rules derived from mutating enhancer reporter constructs may not represent what occurs in the endogenous context.

Shadow enhancers can suppress input transcription factor noise through distinct regulatory logic.

Shadow enhancers, groups of seemingly redundant enhancers, are found in a wide range of organisms and are critical for robust developmental patterning. However, their mechanism of action is unknown. We hypothesized that shadow enhancers drive consistent expression levels by buffering upstream noise through a separation of transcription factor (TF) inputs at the individual enhancers. By measuring the transcriptional dynamics of several Kruppel shadow enhancer configurations in live Drosophila embryos, we showed that individual member enhancers act largely independently. We found that TF fluctuations are an appreciable source of noise that the shadow enhancer pair can better buffer than duplicated enhancers. The shadow enhancer pair is also uniquely able to maintain low levels of expression noise across a wide range of temperatures. A stochastic model demonstrated the separation of TF inputs is sufficient to explain these findings. Our results suggest the widespread use of shadow enhancers is partially due to their noise suppressing ability.

In all higher organisms, life begins with a single cell. During the early stages of development, this single cell grows and divides multiple times to develop into the many different kinds of cells that make up an organism. This is a highly regulated process during which cells receive instructions telling them what kind of cell to become. These instructions are relayed via genes, and a particular combination of activated genes determines the cell's fate. Specific pieces of DNA, known as enhancers, act as switches that control when and where genes are active, while so-called shadow enhancers are found in groups and work together to turn on the same gene in a similar way. Shadow enhancers are often active during the early stages of life to direct the formation of specialized cells in different parts of the body. But so far, it has been unclear why it is beneficial to the divide the role of activating genes across several shadow enhancers rather than a single one. Here, Waymack et al. examined shadow enhancers around a gene called Kruppel in embryos of the fruit fly Drosophila melanogaster. Manipulating the shadow enhancers showed that they help to make gene activity more resistant to changes. Factors such as fluctuations in temperature have different effects on each shadow enhancer. Having several shadow enhancers working together ensures that, whatever happens, the right genes still get activated. For genes like Kruppel, which are key for healthy development, the ability to withstand unexpected changes is a valuable evolutionary benefit. The study of Waymack et al. reveals why shadow enhancers are involved in the regulation of many genes, which may help to better understand developmental defects. Many conditions caused by such defects are influenced by both genetics and the environment. Genetic illnesses can vary in severity, which may be related to the roles of shadow enhancers. As such, studying shadow enhancers could lead to new approaches for treating genetic diseases.

Dissecting the sharp response of a canonical developmental enhancer reveals multiple sources of cooperativity.

Developmental enhancers integrate graded concentrations of transcription factors (TFs) to create sharp gene expression boundaries. Here we examine the hunchback P2 (HbP2) enhancer which drives a sharp expression pattern in the Drosophila blastoderm embryo in response to the transcriptional activator Bicoid (Bcd). We systematically interrogate cis and trans factors that influence the shape and position of expression driven by HbP2, and find that the prevailing model, based on pairwise cooperative binding of Bcd to HbP2 is not adequate. We demonstrate that other proteins, such as pioneer factors, Mediator and histone modifiers influence the shape and position of the HbP2 expression pattern. Comparing our results to theory reveals how higher-order cooperativity and energy expenditure impact boundary location and sharpness. Our results emphasize that the bacterial view of transcription regulation, where pairwise interactions between regulatory proteins dominate, must be reexamined in animals, where multiple molecular mechanisms collaborate to shape the gene regulatory function.

Quantitative Comparison of the Anterior-Posterior Patterning System in the Embryos of Five Drosophila Species.

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.

Hunchback is counter-repressed to regulate even-skipped stripe 2 expression in Drosophila embryos.

Hunchback is a bifunctional transcription factor that can activate and repress gene expression in Drosophila development. We investigated the regulatory DNA sequence features that control Hunchback function by perturbing enhancers for one of its target genes, even-skipped (eve). While Hunchback directly represses the eve stripe 3+7 enhancer, we found that in the eve stripe 2+7 enhancer, Hunchback repression is prevented by nearby sequences-this phenomenon is called counter-repression. We also found evidence that Caudal binding sites are responsible for counter-repression, and that this interaction may be a conserved feature of eve stripe 2 enhancers. Our results alter the textbook view of eve stripe 2 regulation wherein Hb is described as a direct activator. Instead, to generate stripe 2, Hunchback repression must be counteracted. We discuss how counter-repression may influence eve stripe 2 regulation and evolution.

Analysis of Genetic Variation Indicates DNA Shape Involvement in Purifying Selection.

Noncoding DNA sequences, which play various roles in gene expression and regulation, are under evolutionary pressure. Gene regulation requires specific protein-DNA binding events, and our previous studies showed that both DNA sequence and shape readout are employed by transcription factors (TFs) to achieve DNA binding specificity. By investigating the shape-disrupting properties of single nucleotide polymorphisms (SNPs) in human regulatory regions, we established a link between disruptive local DNA shape changes and loss of specific TF binding. Furthermore, we described cases where disease-associated SNPs may alter TF binding through DNA shape changes. This link led us to hypothesize that local DNA shape within and around TF binding sites is under selection pressure. To verify this hypothesis, we analyzed SNP data derived from 216 natural strains of Drosophila melanogaster. Comparing SNPs located in functional and nonfunctional regions within experimentally validated cis-regulatory modules (CRMs) from D. melanogaster that are active in the blastoderm stage of development, we found that SNPs within functional regions tended to cause smaller DNA shape variations. Furthermore, SNPs with higher minor allele frequency were more likely to result in smaller DNA shape variations. The same analysis based on a large number of SNPs in putative CRMs of the D. melanogaster genome derived from DNase I accessibility data confirmed these observations. Taken together, our results indicate that common SNPs in functional regions tend to maintain DNA shape, whereas shape-disrupting SNPs are more likely to be eliminated through purifying selection.

Inference of Transcription Factor Regulation Patterns Using Gene Expression Covariation in Natural Populations of Drosophila melanogaster.

Gene regulatory networks control the complex programs that drive development. Deciphering the connections between transcription factors (TFs) and target genes is challenging, in part because TFs bind to thousands of places in the genome but control expression through a subset of these binding events. We hypothesize that we can combine natural variation of expression levels and predictions of TF binding sites to identify TF targets. We gather RNA-seq data from 71 genetically distinct F1 Drosophila melanogaster embryos and calculate the correlations between TF and potential target genes' expression levels, which we call "regulatory strength." To separate direct and indirect TF targets, we hypothesize that direct TF targets will have a preponderance of binding sites in their upstream regions. Using 14 TFs active during embryogenesis, we find that 12 TFs showed a significant correlation between their binding strength and regulatory strength on downstream targets, and 10 TFs showed a significant correlation between the number of binding sites and the regulatory effect on target genes. The general roles, e.g. bicoid's role as an activator, and the particular interactions we observed between our TFs, e.g. twist's role as a repressor of sloppy paired and odd paired, generally coincide with the literature.

Quantitative Measurement and Thermodynamic Modeling of Fused Enhancers Support a Two-Tiered Mechanism for Interpreting Regulatory DNA.

Computational models of enhancer function generally assume that transcription factors (TFs) exert their regulatory effects independently, modeling an enhancer as a "bag of sites." These models fail on endogenous loci that harbor multiple enhancers, and a "two-tier" model appears better suited: in each enhancer TFs work independently, and the total expression is a weighted sum of their expression readouts. Here, we test these two opposing views on how cis-regulatory information is integrated. We fused two Drosophila blastoderm enhancers, measured their readouts, and applied the above two models to these data. The two-tier mechanism better fits these readouts, suggesting that these fused enhancers comprise multiple independent modules, despite having sequence characteristics typical of single enhancers. We show that short-range TF-TF interactions are not sufficient to designate such modules, suggesting unknown underlying mechanisms. Our results underscore that mechanisms of how modules are defined and how their outputs are combined remain to be elucidated.

An Enhancer's Length and Composition Are Shaped by Its Regulatory Task.

Enhancers drive the gene expression patterns required for virtually every process in metazoans. We propose that enhancer length and transcription factor (TF) binding site composition-the number and identity of TF binding sites-reflect the complexity of the enhancer's regulatory task. In development, we define regulatory task complexity as the number of fates specified in a set of cells at once. We hypothesize that enhancers with more complex regulatory tasks will be longer, with more, but less specific, TF binding sites. Larger numbers of binding sites can be arranged in more ways, allowing enhancers to drive many distinct expression patterns, and therefore cell fates, using a finite number of TF inputs. We compare ~100 enhancers patterning the more complex anterior-posterior (AP) axis and the simpler dorsal-ventral (DV) axis in Drosophila and find that the AP enhancers are longer with more, but less specific binding sites than the (DV) enhancers. Using a set of ~3,500 enhancers, we find enhancer length and TF binding site number again increase with increasing regulatory task complexity. Therefore, to be broadly applicable, computational tools to study enhancers must account for differences in regulatory task.