Natural variation of the expression pattern of the segmentation gene even-skipped in melanogaster

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.

Robust morphogenesis by chaotic dynamics

This research illustrates that complex dynamics of gene products enable the creation of any prescribed cellular differentiation patterns. These complex dynamics can take the form of chaotic, stochastic, or noisy chaotic dynamics. Based on this outcome and previous research, it is established that a generic open chemical reactor can generate an exceptionally large number of different cellular patterns. The mechanism of pattern generation is robust under perturbations and it is based on a combination of Turing's machines, Turing instability and L. Wolpert's gradients. These results can help us to explain the formidable adaptive capacities of biochemical systems.

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.

REDfly: An Integrated Knowledgebase for Insect Regulatory Genomics

We provide here an updated description of the REDfly (Regulatory Element Database for Fly) database of transcriptional regulatory elements, a unique resource that provides regulatory annotation for the genome of Drosophila and other insects. The genomic sequences regulating insect gene expression-transcriptional cis-regulatory modules (CRMs, e.g., "enhancers") and transcription factor binding sites (TFBSs)-are not currently curated by any other major database resources. However, knowledge of such sequences is important, as CRMs play critical roles with respect to disease as well as normal development, phenotypic variation, and evolution. Characterized CRMs also provide useful tools for both basic and applied research, including developing methods for insect control. REDfly, which is the most detailed existing platform for metazoan regulatory-element annotation, includes over 40,000 experimentally verified CRMs and TFBSs along with their DNA sequences, their associated genes, and the expression patterns they direct. Here, we briefly describe REDfly's contents and data model, with an emphasis on the new features implemented since 2020. We then provide an illustrated walk-through of several common REDfly search use cases.

The rate of thermodynamic cost against adiabatic and nonadiabatic fluctuations of a single gene circuit in Drosophila embryos

We study the stochastic dynamics of the externally regulating　gene circuit as an example of eve-skipped gene stripe in the development of Drosophila. Three gene regulation regimes are　considered: adiabatic phase when the switching rate of the gene from the OFF to the ON state is faster than the rate of mRNA degradation; nonadiabatic phase when the switching rate from the OFF to the ON state is slower than that of the mRNA degradation; the bursting phase when the gene switching is fast and transcription is very fast, while the ON state probability is very low. We found that the rate of　thermodynamic cost quantified by the entropy production rate can suppress the fluctuations of the gene circuit. Higher(lower) rate of thermodynamic cost leads to reduced (increased) fluctuations on the number of gene products in the adiabatic (nonadiabatic) regime. We also found that higher thermodynamic cost is often required to sustain the emergence of more gene states and therefore more heterogeneity coming from genetic mutations or epigenetics. We also study the stability of the gene state using the mean first passage time from one state to another. We found the monotonic decrease in time, i.e. on the stability of the state, in the transition from the nonadiabatic to the adiabatic regimes. Therefore, as the higher rate of thermodynamic cost suppresses the fluctuations, higher stability requires higher thermodynamics cost to maintain.

Annotating the Insect Regulatory Genome

An ever-growing number of insect genomes is being sequenced across the evolutionary spectrum. Comprehensive annotation of not only genes but also regulatory regions is critical for reaping the full benefits of this sequencing. Driven by developments in sequencing technologies and in both empirical and computational discovery strategies, the past few decades have witnessed dramatic progress in our ability to identify cis-regulatory modules (CRMs), sequences such as enhancers that play a major role in regulating transcription. Nevertheless, providing a timely and comprehensive regulatory annotation of newly sequenced insect genomes is an ongoing challenge. We review here the methods being used to identify CRMs in both model and non-model insect species, and focus on two tools that we have developed, REDfly and SCRMshaw. These resources can be paired together in a powerful combination to facilitate insect regulatory annotation over a broad range of species, with an accuracy equal to or better than that of other state-of-the-art methods.

Kinetic sculpting of the seven stripes of the Drosophila even-skipped gene

We used live imaging to visualize the transcriptional dynamics of the Drosophila melanogaster even-skipped gene at single-cell and high-temporal resolution as its seven stripe expression pattern forms, and developed tools to characterize and visualize how transcriptional bursting varies over time and space. We find that despite being created by the independent activity of five enhancers, even-skipped stripes are sculpted by the same kinetic phenomena: a coupled increase of burst frequency and amplitude. By tracking the position and activity of individual nuclei, we show that stripe movement is driven by the exchange of bursting nuclei from the posterior to anterior stripe flanks. Our work provides a conceptual, theoretical and computational framework for dissecting pattern formation in space and time, and reveals how the coordinated transcriptional activity of individual nuclei shapes complex developmental patterns.

REDfly: the transcriptional regulatory element database for Drosophila

The REDfly database provides a comprehensive curation of experimentally-validated Drosophila transcriptional cis-regulatory elements and includes information on DNA sequence, experimental evidence, patterns of regulated gene expression, and more. Now in its thirteenth year, REDfly has grown to over 23 000 records of tested reporter gene constructs and 2200 tested transcription factor binding sites. Recent developments include the start of curation of predicted cis-regulatory modules in addition to experimentally-verified ones, improved search and filtering, and increased interaction with the authors of curated papers. An expanded data model that will capture information on temporal aspects of gene regulation, regulation in response to environmental and other non-developmental cues, sexually dimorphic gene regulation, and non-endogenous (ectopic) aspects of reporter gene expression is under development and expected to be in place within the coming year. REDfly is freely accessible at http://redfly.ccr.buffalo.edu, and news about database updates and new features can be followed on Twitter at @REDfly_database.

Evolution of larval segment position across 12 Drosophila species

Many developmental traits that are critical to the survival of the organism are also robust. These robust traits are resistant to phenotypic change in the face of variation. This presents a challenge to evolution. In this article, we asked whether and how a well‐established robust trait, Drosophila segment patterning, changed over the evolutionary history of the genus. We compared segment position scaled to body length at the first‐instar larval stage among 12 Drosophila species. We found that relative segment position has changed many times across the phylogeny. Changes were frequent, but primarily small in magnitude. Phylogenetic analysis demonstrated that rates of change in segment position are variable along the Drosophila phylogenetic tree, and that these changes can occur in short evolutionary timescales. Correlation between position shifts of segments decreased as the distance between two segments increased, suggesting local control of segment position. The posterior‐most abdominal segment showed the highest magnitude of change on average, had the highest rate of evolution between species, and appeared to be evolving more independently as compared to the rest of the segments. This segment was exceptionally elongated in the cactophilic species in our dataset, raising questions as to whether this change may be adaptive.

Lessons and perspectives for applications of stochastic models in biological and cancer research

The effects of randomness, an unavoidable feature of intracellular environments, are observed at higher hierarchical levels of living matter organization, such as cells, tissues, and organisms. Additionally, the many compounds interacting as a well-orchestrated network of reactions increase the difficulties of assessing these systems using only experiments. This limitation indicates that elucidation of the dynamics of biological systems is a complex task that will benefit from the establishment of principles to help describe, categorize, and predict the behavior of these systems. The theoretical machinery already available, or ones to be discovered to help solve biological problems, might play an important role in these processes. Here, we demonstrate the application of theoretical tools by discussing some biological problems that we have approached mathematically: fluctuations in gene expression and cell proliferation in the context of loss of contact inhibition. We discuss the methods that have been employed to provide the reader with a biologically motivated phenomenological perspective of the use of theoretical methods. Finally, we end this review with a discussion of new research perspectives motivated by our results.

Cooperative recruitment of Yan via a high-affinity ETS supersite organizes repression to confer specificity and robustness to cardiac cell fate specification

Here, Boisclair Lachance et al. investigate how the cis-regulatory logic of a tissue-specific cis-regulatory module (CRM) responsible for even-skipped (eve) induction during cardiogenesis organizes the competing inputs of two ETS members: the activator Pointed (Pnt) and the repressor Yan. Their findings provide insight into a novel mechanism by which differential interpretation of CRM syntax by a competing repressor–activator pair can confer both specificity and robustness to developmental transitions.

Cis-regulatory modules (CRMs) are defined by unique combinations of transcription factor-binding sites. Emerging evidence suggests that the number, affinity, and organization of sites play important roles in regulating enhancer output and, ultimately, gene expression. Here, we investigate how the cis-regulatory logic of a tissue-specific CRM responsible for even-skipped (eve) induction during cardiogenesis organizes the competing inputs of two E-twenty-six (ETS) members: the activator Pointed (Pnt) and the repressor Yan. Using a combination of reporter gene assays and CRISPR–Cas9 gene editing, we suggest that Yan and Pnt have distinct syntax preferences. Not only does Yan prefer high-affinity sites, but an overlapping pair of such sites is necessary and sufficient for Yan to tune Eve expression levels in newly specified cardioblasts and block ectopic Eve induction and cell fate specification in surrounding progenitors. Mechanistically, the efficient Yan recruitment promoted by this high-affinity ETS supersite not only biases Yan–Pnt competition at the specific CRM but also organizes Yan-repressive complexes in three dimensions across the eve locus. Taken together, our results uncover a novel mechanism by which differential interpretation of CRM syntax by a competing repressor–activator pair can confer both specificity and robustness to developmental transitions.

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.

A sequence level model of an intact locus predicts the location and function of nonadditive enhancers

Metazoan gene expression is controlled through the action of long stretches of noncoding DNA that contain enhancers-shorter sequences responsible for controlling a single aspect of a gene's expression pattern. Models built on thermodynamics have shown how enhancers interpret protein concentration in order to determine specific levels of gene expression, but the emergent regulatory logic of a complete regulatory locus shows qualitative and quantitative differences from isolated enhancers. Such differences may arise from steric competition limiting the quantity of DNA that can simultaneously influence the transcription machinery. We incorporated this competition into a mechanistic model of gene regulation, generated efficient algorithms for this computation, and applied it to the regulation of Drosophila even-skipped (eve). This model finds the location of enhancers and identifies which factors control the boundaries of eve expression. This model predicts a new enhancer that, when assayed in vivo, drives expression in a non-eve pattern. Incorporation of chromatin accessibility eliminates this inconsistency.