Rapid transcription fosters coordinate snail expression in the Drosophila embryo

Transcript accumulation rates in the early Caenorhabditis elegans embryo

Dynamic transcriptional changes are widespread in rapidly dividing developing embryos when cell fate decisions are made quickly. The Caenorhabditis elegans embryo overcomes these constraints partly through the rapid production of high levels of transcription factor mRNAs. Transcript accumulation rates for some developmental genes are known at single-cell resolution, but genome-scale measurements are lacking. We estimate zygotic mRNA accumulation rates from single-cell RNA sequencing data calibrated with single-molecule transcript imaging. Rapid transcription is common in the early C. elegans embryo with rates highest soon after zygotic transcription begins. High-rate genes are enriched for recently duplicated cell-fate regulators and share common genomic features. We identify core promoter elements associated with high rate and measure their contributions for two early endomesodermal genes, ceh-51 and sdz-31. Individual motifs modestly affect accumulation rates, suggesting multifactorial control. These results are a step toward estimating absolute transcription kinetics and understanding how transcript dosage drives developmental decisions.

Combined modelling of mRNA decay dynamics and single-molecule imaging in the Drosophila embryo uncovers a role for P-bodies in 5' to 3' degradation

Regulation of mRNA degradation is critical for a diverse array of cellular processes and developmental cell fate decisions. Many methods for determining mRNA half-lives rely on transcriptional inhibition or metabolic labelling. Here, we use a non-invasive method for estimating half-lives for hundreds of mRNAs in the early Drosophila embryo. This approach uses the intronic and exonic reads from a total RNA-seq time series and Gaussian process regression to model the dynamics of premature and mature mRNAs. We show how regulation of mRNA stability is used to establish a range of mature mRNA dynamics during embryogenesis, despite shared transcription profiles. Using single-molecule imaging, we provide evidence that, for the mRNAs tested, there is a correlation between short half-life and mRNA association with P-bodies. Moreover, we detect an enrichment of mRNA 3' ends in P-bodies in the early embryo, consistent with 5' to 3' degradation occurring in P-bodies for at least a subset of mRNAs. We discuss our findings in relation to recently published data suggesting that the primary function of P-bodies in other biological contexts is mRNA storage.

The scaffolding function of LSD1/KDM1A reinforces a negative feedback loop to repress stem cell gene expression during primitive hematopoiesis

Lsd1/Kdm1a functions both as a histone demethylase enzyme and as a scaffold for assembling chromatin modifier and transcription factor complexes to regulate gene expression. The relative contributions of Lsd1's demethylase and scaffolding functions during embryogenesis are not known. Here, we analyze two independent zebrafish lsd1/kdm1a mutant lines and show Lsd1 is required to repress primitive hematopoietic stem cell gene expression. Lsd1 rescue constructs containing point mutations that selectively abrogate its demethylase or scaffolding capacity demonstrate the scaffolding function of Lsd1, not its demethylase activity, is required for repression of gene expression in vivo. Lsd1's SNAG-binding domain mediates its scaffolding function and reinforces a negative feedback loop to repress the expression of SNAG-domain-containing genes during embryogenesis, including gfi1 and snai1/2. Our findings reveal a model in which the SNAG-binding and scaffolding function of Lsd1, and its associated negative feedback loop, provide transient and reversible regulation of gene expression during hematopoietic development.

Incorporation of CENP-A/CID into centromeres during early Drosophila embryogenesis does not require RNA polymerase II-mediated transcription

In many species, centromere identity is specified epigenetically by special nucleosomes containing a centromere-specific histone H3 variant, designated as CENP-A in humans and CID in Drosophila melanogaster. After partitioning of centromere-specific nucleosomes onto newly replicated sister centromeres, loading of additional CENP-A/CID into centromeric chromatin is required for centromere maintenance in proliferating cells. Analyses with cultured cells have indicated that transcription of centromeric DNA by RNA polymerase II is required for deposition of new CID into centromere chromatin. However, a dependence of centromeric CID loading on transcription is difficult to reconcile with the notion that the initial embryonic stages appear to proceed in the absence of transcription in Drosophila, as also in many other animal species. To address the role of RNA polymerase II–mediated transcription for CID loading in early Drosophila embryos, we have quantified the effects of alpha-amanitin and triptolide on centromeric CID-EGFP levels. Our analyses demonstrate that microinjection of these two potent inhibitors of RNA polymerase II–mediated transcription has at most a marginal effect on centromeric CID deposition during progression through the early embryonic cleavage cycles. Thus, we conclude that at least during early Drosophila embryogenesis, incorporation of CID into centromeres does not depend on RNA polymerase II–mediated transcription.

Oxyresveratrol inhibits human colon cancer cell migration through regulating epithelial-mesenchymal transition and microRNA

The major cause of death in colorectal cancer (CRC) patients is metastasis. Moreover, lots of studies have emphasized that the epithelial-mesenchymal transition (EMT) is a pivotal step in metastasis. Both transforming growth factor beta (TGF-β) and dysregulation of microRNAs (miRNAs) can induce or regulate EMT, promoting the loss of intercellular adhesion and increased motility of cancer cells. Therefore, it is necessary to prevent or inhibit the metastasis of colorectal cancer. Relatively little is known about the anti-metastatic effect of oxyresveratrol (OXY), a natural derivative of resveratrol (RES), compared to RES. Accordingly, RES was used as the positive control to investigate the effects of OXY on colon cancer cell migration. The results showed that OXY could significantly inhibit cell migration (67.17% ± 0.04, 64.89% ± 0.04) compared to RES (84.6% ± 0.07, 76.34% ± 0.08) in HCT116 cells and TGF-β-induced HT-29 cells, respectively, via Snail/E-cadherin expression. In addition, OXY improved EMT-related miRNA expression through, for example, lowering the levels of miR-3687 and miR-301a-3p while upregulating miR-3612 in TGF-β-induced HT-29 cells. In conclusion, OXY inhibits human colon cancer cell migration by regulating EMT and miRNAs. Based on these findings, it can be stated that OXY promotes anti-metastatic properties in CRC.

Quantitative imaging of transcription in living Drosophila embryos reveals the impact of core promoter motifs on promoter state dynamics

Genes are expressed in stochastic transcriptional bursts linked to alternating active and inactive promoter states. A major challenge in transcription is understanding how promoter composition dictates bursting, particularly in multicellular organisms. We investigate two key Drosophila developmental promoter motifs, the TATA box (TATA) and the Initiator (INR). Using live imaging in Drosophila embryos and new computational methods, we demonstrate that bursting occurs on multiple timescales ranging from seconds to minutes. TATA-containing promoters and INR-containing promoters exhibit distinct dynamics, with one or two separate rate-limiting steps respectively. A TATA box is associated with long active states, high rates of polymerase initiation, and short-lived, infrequent inactive states. In contrast, the INR motif leads to two inactive states, one of which relates to promoter-proximal polymerase pausing. Surprisingly, the model suggests pausing is not obligatory, but occurs stochastically for a subset of polymerases. Overall, our results provide a rationale for promoter switching during zygotic genome activation.

Dynamic regulation of anterior-posterior patterning genes in living Drosophila embryos

Expression of the gap and pair-rule genes plays an essential role in body segmentation during Drosophila embryogenesis.^1-5 However, it remains unclear how precise expression patterns of these key developmental genes arise from stochastic transcriptional activation at the single-cell level. Here, I employed genome-editing and live-imaging approaches to comprehensively visualize regulation of the gap and pair-rule genes at the endogenous loci. Quantitative image analysis revealed that the total duration of active transcription (transcription period) is a major determinant of spatial patterning of gene expression in early embryos. The length of the transcription period is determined by the continuity of bursting activities in individual nuclei, with the core expression domain producing more bursts than boundary regions. Each gene exhibits a distinct rate of nascent RNA production during transcriptional bursting, which contributes to gene-to-gene variability in the total output. I also provide evidence for "enhancer interference," wherein a distal weak enhancer interferes with transcriptional activation by a strong proximal enhancer to downregulate the length of the transcription period without changing the transcription rate. Analysis of the endogenous hunchback (hb) locus revealed that the removal of the distal shadow enhancer induces strong ectopic transcriptional activation, which suppresses refinement of the initial broad expression domain into narrower stripe patterns at the anterior part of embryos. This study provides key insights into the link between transcriptional bursting, enhancer-promoter interaction, and spatiotemporal patterning of gene expression during animal development.

Real-time single-cell characterization of the eukaryotic transcription cycle reveals correlations between RNA initiation, elongation, and cleavage

The eukaryotic transcription cycle consists of three main steps: initiation, elongation, and cleavage of the nascent RNA transcript. Although each of these steps can be regulated as well as coupled with each other, their in vivo dissection has remained challenging because available experimental readouts lack sufficient spatiotemporal resolution to separate the contributions from each of these steps. Here, we describe a novel application of Bayesian inference techniques to simultaneously infer the effective parameters of the transcription cycle in real time and at the single-cell level using a two-color MS2/PP7 reporter gene and the developing fruit fly embryo as a case study. Our method enables detailed investigations into cell-to-cell variability in transcription-cycle parameters as well as single-cell correlations between these parameters. These measurements, combined with theoretical modeling, suggest a substantial variability in the elongation rate of individual RNA polymerase molecules. We further illustrate the power of this technique by uncovering a novel mechanistic connection between RNA polymerase density and nascent RNA cleavage efficiency. Thus, our approach makes it possible to shed light on the regulatory mechanisms in play during each step of the transcription cycle in individual, living cells at high spatiotemporal resolution.

Tracing DNA paths and RNA profiles in cultured cells and tissues with ORCA

Chromatin conformation capture (3C) methods and fluorescent in situ hybridization (FISH) microscopy have been used to investigate the spatial organization of the genome. Although powerful, both techniques have limitations. Hi-C is challenging for low cell numbers and requires very deep sequencing to achieve its high resolution. In contrast, FISH can be done on small cell numbers and capture rare cell populations, but typically targets pairs of loci at a lower resolution. Here we detail a protocol for optical reconstruction of chromatin architecture (ORCA), a microscopy approach to trace the 3D DNA path within the nuclei of fixed tissues and cultured cells with a genomic resolution as fine as 2 kb and a throughput of ~10,000 cells per experiment. ORCA can identify structural features with comparable resolution to Hi-C while providing single-cell resolution and multimodal measurements characteristic of microscopy. We describe how to use this DNA labeling in parallel with multiplexed labeling of dozens of RNAs to relate chromatin structure and gene expression in the same cells. Oligopaint probe design, primary probe making, sample collection, cryosectioning and RNA/DNA primary probe hybridization can be completed in 1.5 weeks, while automated RNA/DNA barcode hybridization and RNA/DNA imaging typically takes 2-6 d for data collection and 2-7 d for the automated steps of image analysis.

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.

Does Gene Expression Noise Play a Functional Role in Plants?

Gene expression in individual cells can be surprisingly noisy. In unicellular organisms this noise can be functional; for example, by allowing a subfraction of the population to prepare for environmental stress. The role of gene expression noise in multicellular organisms has, however, remained unclear. In this review, we discuss how new techniques are revealing an unexpected level of variability in gene expression between and within genetically identical plants. We describe recent progress as well as speculate on the function of transcriptional noise as a mechanism for generating functional phenotypic diversity in plants.

Modulation of the Promoter Activation Rate Dictates the Transcriptional Response to Graded BMP Signaling Levels in the Drosophila Embryo

Morphogen gradients specify cell fates during development, with a classic example being the bone morphogenetic protein (BMP) gradient's conserved role in embryonic dorsal-ventral axis patterning. Here, we elucidate how the BMP gradient is interpreted in the Drosophila embryo by combining live imaging with computational modeling to infer transcriptional burst parameters at single-cell resolution. By comparing burst kinetics in cells receiving different levels of BMP signaling, we show that BMP signaling controls burst frequency by regulating the promoter activation rate. We provide evidence that the promoter activation rate is influenced by both enhancer and promoter sequences, whereas Pol II loading rate is primarily modulated by the enhancer. Consistent with BMP-dependent regulation of burst frequency, the numbers of BMP target gene transcripts per cell are graded across their expression domains. We suggest that graded mRNA output is a general feature of morphogen gradient interpretation and discuss how this can impact on cell-fate decisions.

Twist-dependent ratchet functioning downstream from Dorsal revealed using a light-inducible degron

Graded transcription factors are pivotal regulators of embryonic patterning, but whether their role changes over time is unclear. A light-regulated protein degradation system was used to assay temporal dependence of the transcription factor Dorsal in dorsal-ventral axis patterning of Drosophila embryos. Surprisingly, the high-threshold target gene snail only requires Dorsal input early but not late when Dorsal levels peak. Instead, late snail expression can be supported by action of the Twist transcription factor, specifically, through one enhancer, sna.distal This study demonstrates that continuous input is not required for some Dorsal targets and downstream responses, such as twist, function as molecular ratchets.

A technical review and guide to RNA fluorescence in situ hybridization

RNA-fluorescence in situ hybridization (FISH) is a powerful tool to visualize target messenger RNA transcripts in cultured cells, tissue sections or whole-mount preparations. As the technique has been developed over time, an ever-increasing number of divergent protocols have been published. There is now a broad selection of options available to facilitate proper tissue preparation, hybridization, and post-hybridization background removal to achieve optimal results. Here we review the technical aspects of RNA-FISH, examining the most common methods associated with different sample types including cytological preparations and whole-mounts. We discuss the application of commonly used reagents for tissue preparation, hybridization, and post-hybridization washing and provide explanations of the functional roles for each reagent. We also discuss the available probe types and necessary controls to accurately visualize gene expression. Finally, we review the most recent advances in FISH technology that facilitate both highly multiplexed experiments and signal amplification for individual targets. Taken together, this information will guide the methods development process for investigators that seek to perform FISH in organisms that lack documented or optimized protocols.

Ordered patterning of the sensory system is susceptible to stochastic features of gene expression

Sensory neuron numbers and positions are precisely organized to accurately map environmental signals in the brain. This precision emerges from biochemical processes within and between cells that are inherently stochastic. We investigated impact of stochastic gene expression on pattern formation, focusing on senseless (sens), a key determinant of sensory fate in Drosophila. Perturbing microRNA regulation or genomic location of sens produced distinct noise signatures. Noise was greatly enhanced when both sens alleles were present in homologous loci such that each allele was regulated in trans by the other allele. This led to disordered patterning. In contrast, loss of microRNA repression of sens increased protein abundance but not sensory pattern disorder. This suggests that gene expression stochasticity is a critical feature that must be constrained during development to allow rapid yet accurate cell fate resolution.

Direct and simultaneous observation of transcription and chromosome architecture in single cells with Hi-M

Simultaneous observation of 3D chromatin organization and transcription at the single-cell level and with high spatial resolution may hold the key to unveiling the mechanisms regulating embryonic development, cell differentiation and even disease. We recently developed Hi-M, a technology that enables the sequential labeling, 3D imaging and localization of multiple genomic DNA loci, together with RNA expression, in single cells within whole, intact Drosophila embryos. Importantly, Hi-M enables simultaneous detection of RNA expression and chromosome organization without requiring sample unmounting and primary probe rehybridization. Here, we provide a step-by-step protocol describing the design of probes, the preparation of samples, the stable immobilization of embryos in microfluidic chambers, and the complete procedure for image acquisition. The combined RNA/DNA fluorescence in situ hybridization procedure takes 4-5 d, including embryo collection. In addition, we describe image analysis software to segment nuclei, detect genomic spots, correct for drift and produce Hi-M matrices. A typical Hi-M experiment takes 1-2 d to complete all rounds of labeling and imaging and 4 additional days for image analysis. This technology can be easily expanded to investigate cell differentiation in cultured cells or organization of chromatin within complex tissues.

Multimodal transcriptional control of pattern formation in embryonic development

Predicting how interactions between transcription factors and regulatory DNA sequence dictate rates of transcription and, ultimately, drive developmental outcomes remains an open challenge in physical biology. Using stripe 2 of the even-skipped gene in Drosophila embryos as a case study, we dissect the regulatory forces underpinning a key step along the developmental decision-making cascade: the generation of cytoplasmic mRNA patterns via the control of transcription in individual cells. Using live imaging and computational approaches, we found that the transcriptional burst frequency is modulated across the stripe to control the mRNA production rate. However, we discovered that bursting alone cannot quantitatively recapitulate the formation of the stripe and that control of the window of time over which each nucleus transcribes even-skipped plays a critical role in stripe formation. Theoretical modeling revealed that these regulatory strategies (bursting and the time window) respond in different ways to input transcription factor concentrations, suggesting that the stripe is shaped by the interplay of 2 distinct underlying molecular processes.

Spatiotemporal control of gene expression boundaries using a feedforward loop

BACKGROUND

A feedforward loop (FFL) is commonly observed in several biological networks. The FFL network motif has been mostly studied with respect to variation of the input signal in time, with only a few studies of FFL activity in a spatially distributed system such as morphogen-mediated tissue patterning. However, most morphogen gradients also evolve in time.

RESULTS

We studied the spatiotemporal behavior of a coherent FFL in two contexts: (a) a generic, oscillating morphogen gradient and (b) the dorsal-ventral patterning of the early Drosophila embryo by a gradient of the NF-κB homolog dorsal with its early target Twist. In both models, we found features in the dynamics of the intermediate node-phase difference and noise filtering-that were largely independent of the parameterization of the models, and thus were functions of the structure of the FFL itself. In the dorsal gradient model, we also found that proper target gene expression was not possible without including the effect of maternal pioneer factor Zelda.

CONCLUSIONS

An FFL buffers fluctuation to changes in the morphogen signal ensuring stable gene expression boundaries.

Noise in the Vertebrate Segmentation Clock Is Boosted by Time Delays but Tamed by Notch Signaling

Taming cell-to-cell variability in gene expression is critical for precise pattern formation during embryonic development. To investigate the source and buffering mechanism of expression variability, we studied a biological clock, the vertebrate segmentation clock, controlling the precise spatiotemporal patterning of the vertebral column. By counting single transcripts of segmentation clock genes in zebrafish, we show that clock genes have low RNA amplitudes and expression variability is primarily driven by gene extrinsic sources, which is suppressed by Notch signaling. We further show that expression noise surprisingly increases from the posterior progenitor zone to the anterior segmentation and differentiation zone. Our computational model reproduces the spatial noise profile by incorporating spatially increasing time delays in gene expression. Our results, suggesting that expression variability is controlled by the balance of time delays and cell signaling in a vertebrate tissue, will shed light on the accuracy of natural clocks in multi-cellular systems and inspire engineering of robust synthetic oscillators.

Lighting up the central dogma for predictive developmental biology

Although the last 30years have witnessed the mapping of the wiring diagrams of the gene regulatory networks that dictate cell fate and animal body plans, specific understanding building on such network diagrams that shows how DNA regulatory regions control gene expression lags far behind. These networks have yet to yield the predictive power necessary to, for example, calculate how the concentration dynamics of input transcription factors and DNA regulatory sequence prescribes output patterns of gene expression that, in turn, determine body plans themselves. Here, we argue that reaching a predictive understanding of developmental decision-making calls for an interplay between theory and experiment aimed at revealing how the regulation of the processes of the central dogma dictate network connections and how network topology guides cells toward their ultimate developmental fate. To make this possible, it is crucial to break free from the snapshot-based understanding of embryonic development facilitated by fixed-tissue approaches and embrace new technologies that capture the dynamics of developmental decision-making at the single cell level, in living embryos.

Control of Hox transcription factor concentration and cell-to-cell variability by an auto-regulatory switch

The variability in transcription factor concentration among cells is an important developmental determinant, yet how variability is controlled remains poorly understood. Studies of variability have focused predominantly on monitoring mRNA production noise. Little information exists about transcription factor protein variability, as this requires the use of quantitative methods with single-molecule sensitivity. Using Fluorescence Correlation Spectroscopy (FCS), we have characterized the concentration and variability of 14 endogenously tagged TFs in live Drosophila imaginal discs. For the Hox TF Antennapedia, we investigated whether protein variability results from random stochastic events or is developmentally regulated. We found that Antennapedia transitioned from low concentration/high variability early, to high concentration/low variability later, in development. FCS and temporally resolved genetic studies uncovered that Antennapedia itself is necessary and sufficient to drive a developmental regulatory switch from auto-activation to auto-repression, thereby reducing variability. This switch is controlled by progressive changes in relative concentrations of preferentially activating and repressing Antennapedia isoforms, which bind chromatin with different affinities. Mathematical modeling demonstrated that the experimentally supported auto-regulatory circuit can explain the increase of Antennapedia concentration and suppression of variability over time.

Summary: Preferentially repressing and activating isoforms of the Hox transcription factor Antennapedia elicit a developmental regulatory switch from auto-activation to auto-repression that increases concentration and suppresses cell-to-cell variability over time.

Detection and Automated Analysis of Single Transcripts at Subcellular Resolution in Zebrafish Embryos

Single molecule fluorescence in situ hybridization (smFISH) is a method to visualize single mRNA molecules. When combined with cellular and nuclear segmentation, transcripts can be assigned to different cellular compartments resulting in quantitative information on transcript levels at subcellular resolution. The use of smFISH in zebrafish has been limited by the lack of protocols and an automated image analysis pipeline for samples of multicellular organisms. Here we present a protocol for smFISH on zebrafish cryosections. The protocol includes a method to obtain high-quality sections of zebrafish embryos, an smFISH protocol optimized for zebrafish cryosections, and a user-friendly, automated analysis pipeline for cell segmentation and transcript detection. The software is freely available and can be used to analyze sections of any multicellular organism.

Uniform gene expression in embryos is achieved by temporal averaging of transcription noise

Here, Stapel et al. investigated how stochastic transcription activation can result in uniform expression during development of cellularized embryos. Using smFISH and mathematical modeling, the authors show that during zebrafish embryogenesis, transcription activation is stochastic due to (1) genes acquiring transcriptional competence at different times in different cells, (2) differences in cell cycle stage between cells, and (3) the stochastic nature of transcription.

Transcription is often stochastic. This is seemingly incompatible with the importance of gene expression during development. Here we show that during zebrafish embryogenesis, transcription activation is stochastic due to (1) genes acquiring transcriptional competence at different times in different cells, (2) differences in cell cycle stage between cells, and (3) the stochastic nature of transcription. Initially, stochastic transcription causes large cell-to-cell differences in transcript levels. However, variability is reduced by lengthening cell cycles and the accumulation of transcription events in each cell. Temporal averaging might provide a general context in which to understand how embryos deal with stochastic transcription.

Shaping development by stochasticity and dynamics in gene regulation

Animal development is orchestrated by spatio-temporal gene expression programmes that drive precise lineage commitment, proliferation and migration events at the single-cell level, collectively leading to large-scale morphological change and functional specification in the whole organism. Efforts over decades have uncovered two 'seemingly contradictory' mechanisms in gene regulation governing these intricate processes: (i) stochasticity at individual gene regulatory steps in single cells and (ii) highly coordinated gene expression dynamics in the embryo. Here we discuss how these two layers of regulation arise from the molecular and the systems level, and how they might interplay to determine cell fate and to control the complex body plan. We also review recent technological advancements that enable quantitative analysis of gene regulation dynamics at single-cell, single-molecule resolution. These approaches outline next-generation experiments to decipher general principles bridging gaps between molecular dynamics in single cells and robust gene regulations in the embryo.

Enhancer Control of Transcriptional Bursting

Transcription is episodic, consisting of a series of discontinuous bursts. Using live-imaging methods and quantitative analysis, we examine transcriptional bursting in living Drosophila embryos. Different developmental enhancers positioned downstream of synthetic reporter genes produce transcriptional bursts with similar amplitudes and duration but generate very different bursting frequencies, with strong enhancers producing more bursts than weak enhancers. Insertion of an insulator reduces the number of bursts and the corresponding level of gene expression, suggesting that enhancer regulation of bursting frequency is a key parameter of gene control in development. We also show that linked reporter genes exhibit coordinated bursting profiles when regulated by a shared enhancer, challenging conventional models of enhancer-promoter looping.

Noise propagation with interlinked feed-forward pathways

Functionally similar pathways are often seen in biological systems, forming feed-forward controls. The robustness in network motifs such as feed-forward loops (FFLs) has been reported previously. In this work, we studied noise propagation in a development network that has multiple interlinked FFLs. A FFL has the potential of asymmetric noise-filtering (i.e., it works at either the “ON” or the “OFF” state in the target gene). With multiple, interlinked FFLs, we show that the propagated noises are largely filtered regardless of the states in the input genes. The noise-filtering property of an interlinked FFL can be largely derived from that of the individual FFLs, and with interlinked FFLs, it is possible to filter noises in both “ON” and “OFF” states in the output. We demonstrated the noise filtering effect in the developmental regulatory network of Caenorhabditis elegans that controls the timing of distal tip cell (DTC) migration. The roles of positive feedback loops involving blmp-1 and the degradation regulation of DRE-1 also studied. Our analyses allow for better inference from network structures to noise-filtering properties, and provide insights into the mechanisms behind the precise DTC migration controls in space and time.

New methods to image transcription in living fly embryos: the insights so far, and the prospects

The regulation of transcription is a fundamental process underlying the determination of cell identity and its maintenance during development. In the last decades, most of the transcription factors, which have to be expressed at the right place and at the right time for the proper development of the fly embryo, have been identified. However, mostly because of the lack of methods to visualize transcription as the embryo develops, their coordinated spatiotemporal dynamics remains largely unexplored. Efforts have been made to decipher the transcription process with single molecule resolution at the single cell level. Recently, the fluorescent labeling of nascent RNA in developing fly embryos allowed the direct visualization of ongoing transcription at single loci within each nucleus. Together with powerful imaging and quantitative data analysis, these new methods provide unprecedented insights into the temporal dynamics of the transcription process and its intrinsic noise. Focusing on the Drosophila embryo, we discuss how the detection of single RNA molecules enhanced our comprehension of the transcription process and we outline the potential next steps made possible by these new imaging tools. In combination with genetics and theoretical analysis, these new imaging methods will aid the search for the mechanisms responsible for the robustness of development. For further resources related to this article, please visit the WIREs website.

Quantitative Single-Embryo Profile of Drosophila Genome Activation and the Dorsal-Ventral Patterning Network

During embryonic development of Drosophila melanogaster, the maternal-to-zygotic transition (MZT) marks a significant and rapid turning point when zygotic transcription begins and control of development is transferred from maternally deposited transcripts. Characterizing the sequential activation of the genome during the MZT requires precise timing and a sensitive assay to measure changes in expression. We utilized the NanoString nCounter instrument, which directly counts messenger RNA transcripts without reverse transcription or amplification, to study >70 genes expressed along the dorsal-ventral (DV) axis of early Drosophila embryos, dividing the MZT into 10 time points. Transcripts were quantified for every gene studied at all time points, providing the first dataset of absolute numbers of transcripts during Drosophila development. We found that gene expression changes quickly during the MZT, with early nuclear cycle 14 (NC14) the most dynamic time for the embryo. twist is one of the most abundant genes in the entire embryo and we use mutants to quantitatively demonstrate how it cooperates with Dorsal to activate transcription and is responsible for some of the rapid changes in transcription observed during early NC14. We also uncovered elements within the gene regulatory network that maintain precise transcript levels for sets of genes that are spatiotemporally cotranscribed within the presumptive mesoderm or dorsal ectoderm. Using these new data, we show that a fine-scale, quantitative analysis of temporal gene expression can provide new insights into developmental biology by uncovering trends in gene networks, including coregulation of target genes and specific temporal input by transcription factors.

Bicoid gradient formation and function in the Drosophila pre-syncytial blastoderm

Bicoid (Bcd) protein distributes in a concentration gradient that organizes the anterior/posterior axis of the Drosophila embryo. It has been understood that bcd RNA is sequestered at the anterior pole during oogenesis, is not translated until fertilization, and produces a protein gradient that functions in the syncytial blastoderm after 9–10 nuclear divisions. However, technical issues limited the sensitivity of analysis of pre-syncytial blastoderm embryos and precluded studies of oocytes after stage 13. We developed methods to analyze stage 14 oocytes and pre-syncytial blastoderm embryos, and found that stage 14 oocytes make Bcd protein, that bcd RNA and Bcd protein distribute in matching concentration gradients in the interior of nuclear cycle 2–6 embryos, and that Bcd regulation of target gene expression is apparent at nuclear cycle 7, two cycles prior to syncytial blastoderm. We discuss the implications for the generation and function of the Bcd gradient.

DOI:http://dx.doi.org/10.7554/eLife.13222.001

As an embryo develops, a single cell transforms into a collection of different types of cells. One protein that is crucial for this process in fruit fly embryos is Bicoid. Thirty years ago, scientists discovered that Bicoid protein is concentrated at the head end of the embryo and gradually decreases in amount towards the rear end. This concentration gradient of Bicoid protein organizes the embryo body and regulates the expression of many genes, thus directing the cells to develop different identities.

Several assumptions had been made about how this gradient is established. It was thought that in the unfertilized egg, the mRNA molecules that will be translated to produce Bicoid proteins are stored in an inactive state in the region of the egg that later develops into the embryo’s head. In the embryo, the mRNA molecules were believed to remain in the head region while being translated, with the newly formed proteins then gradually spreading from this site to create the Bicoid gradient. It was also thought that no Bicoid proteins are stored in the unfertilized egg. However, no known methods were sensitive enough to investigate these assumptions.

Now, using newer and more sensitive methods, Ali-Murthy and Kornberg show that Bicoid protein is present in the unfertilized fruit fly egg in the same region as the mRNA molecules that make Bicoid. Furthermore, the Bicoid gradient forms when the embryo has fewer than 32 nuclei, much earlier in development than previously thought. The Bicoid protein also does not appear to spread passively towards the rear of the embryo, but is transported in a more orchestrated manner.

Overall, Ali-Murthy and Kornberg’s results suggest that the early fruit fly embryo is more organized and actively regulated than had been previously understood. This paves the way for further studies that use sensitive techniques to investigate this early stage of development.

DOI:http://dx.doi.org/10.7554/eLife.13222.002

Transcriptional Memory in the Drosophila Embryo

Transmission of active transcriptional states from mother to daughter cells has the potential to foster precision in the gene expression programs underlying development. Such transcriptional memory has been specifically proposed to promote rapid reactivation of complex gene expression profiles after successive mitoses in Drosophila development [1]. By monitoring transcription in living Drosophila embryos, we provide the first evidence for transcriptional memory in animal development. We specifically monitored the activities of stochastically expressed transgenes in order to distinguish active and inactive mother cells and the behaviors of their daughter nuclei after mitosis. Quantitative analyses reveal that there is a 4-fold higher probability for rapid reactivation after mitosis when the mother experienced transcription. Moreover, memory nuclei activate transcription twice as fast as neighboring inactive mothers, thus leading to augmented levels of gene expression. We propose that transcriptional memory is a mechanism of precision, which helps coordinate gene activity during embryogenesis.

Automated detection and quantification of single RNAs at cellular resolution in zebrafish embryos

Analysis of differential gene expression is crucial for the study of cell fate and behavior during embryonic development. However, automated methods for the sensitive detection and quantification of RNAs at cellular resolution in embryos are lacking. With the advent of single-molecule fluorescence in situ hybridization (smFISH), gene expression can be analyzed at single-molecule resolution. However, the limited availability of protocols for smFISH in embryos and the lack of efficient image analysis pipelines have hampered quantification at the (sub)cellular level in complex samples such as tissues and embryos. Here, we present a protocol for smFISH on zebrafish embryo sections in combination with an image analysis pipeline for automated transcript detection and cell segmentation. We use this strategy to quantify gene expression differences between different cell types and identify differences in subcellular transcript localization between genes. The combination of our smFISH protocol and custom-made, freely available, analysis pipeline will enable researchers to fully exploit the benefits of quantitative transcript analysis at cellular and subcellular resolution in tissues and embryos.

The expression of the T-box selector gene midline in the leg imaginal disc is controlled by both transcriptional regulation and cell lineage

The Drosophila Tbx20 homologs midline and H15 act as selector genes for ventral fate in Drosophila legs. midline and H15 expression defines the ventral domain of the leg and the two genes are necessary and sufficient for the development of ventral fate. Ventral-specific expression of midline and H15 is activated by Wingless (Wg) and repressed by Decapentaplegic (Dpp). Here we identify VLE, a 5 kb enhancer that drives ventral specific expression in the leg disc that is very similar to midline expression. Subdivision of VLE identifies two regions that mediate both activation and repression and third region that only mediates repression. Loss- and gain-of-function genetic mosaic analysis shows that the activating and repressing regions respond to Wg and Dpp signaling respectively. All three repression regions depend on the activity of Mothers-against-decapentaplegic, a Drosophila r-Smad that mediates Dpp signaling, and respond to ectopic expression of the Dpp target genes optomoter-blind and Dorsocross 3. However, only one repression region is responsive to loss of schnurri, a co-repressor required for direct repression by Dpp-signaling. Thus, Dpp signaling restricts midline expression through both direct repression and through the activation of downstream repressors. We also find that midline and H15 expression are both subject to cross-repression and feedback inhibition. Finally, a lineage analysis indicates that ventral midline-expressing cells and dorsal omb-expressing cells do not mix during development. Together this data indicates that the ventral-specific expression of midline results from both transcriptional regulation and from a lack of cell-mixing between dorsal and ventral cells.

Deciphering Transcriptional Dynamics In Vivo by Counting Nascent RNA Molecules

Deciphering how the regulatory DNA sequence of a gene dictates its expression in response to intra and extracellular cues is one of the leading challenges in modern genomics. The development of novel single-cell sequencing and imaging techniques, as well as a better exploitation of currently available single-molecule imaging techniques, provides an avenue to interrogate the process of transcription and its dynamics in cells by quantifying the number of RNA polymerases engaged in the transcription of a gene (or equivalently the number of nascent RNAs) at a given moment in time. In this paper, we propose that measurements of the cell-to-cell variability in the number of nascent RNAs provide a mostly unexplored method for deciphering mechanisms of transcription initiation in cells. We propose a simple kinetic model of transcription initiation and elongation from which we calculate nascent RNA copy-number fluctuations. To demonstrate the usefulness of this approach, we test our theory against published nascent RNA data for twelve constitutively expressed yeast genes. Rather than transcription being initiated through a single rate limiting step, as it had been previously proposed, our single-cell analysis reveals the presence of at least two rate limiting steps. Surprisingly, half of the genes analyzed have nearly identical rates of transcription initiation, suggesting a common mechanism. Our analytical framework can be used to extract quantitative information about dynamics of transcription from single-cell sequencing data, as well as from single-molecule imaging and electron micrographs of fixed cells, and provides the mathematical means to exploit the quantitative power of these technologies.

Combining protein and mRNA quantification to decipher transcriptional regulation

We combine immunofluorescence and single-molecule fluorescence in situ hybridization (smFISH), followed by automated image analysis, to quantify the concentration of nuclear transcription factors, number of transcription factors bound, and number of nascent mRNAs synthesized at individual gene loci. A theoretical model is used to decipher how transcription-factor binding modulates the stochastic kinetics of mRNA production. We demonstrate this approach by examining the regulation of hunchback in the early Drosophila embryo.

Bursty gene expression in the intact mammalian liver

Bursts of nascent mRNA have been shown to lead to substantial cell-cell variation in unicellular organisms, facilitating diverse responses to environmental challenges. It is unknown whether similar bursts and gene-expression noise occur in mammalian tissues. To address this, we combine single molecule transcript counting with dual-color labeling and quantification of nascent mRNA to characterize promoter states, transcription rates, and transcript lifetimes in the intact mouse liver. We find that liver gene expression is highly bursty, with promoters stochastically switching between transcriptionally active and inactive states. Promoters of genes with short mRNA lifetimes are active longer, facilitating rapid response while reducing burst-associated noise. Moreover, polyploid hepatocytes exhibit less noise than diploid hepatocytes, suggesting a possible benefit to liver polyploidy. Thus, temporal averaging and liver polyploidy dampen the intrinsic variability associated with transcriptional bursts. Our approach can be used to study transcriptional bursting in diverse mammalian tissues.

Hormone-induced and DNA demethylation-induced relief of a tissue-specific and developmentally regulated block in transcriptional elongation

Genome-wide studies have revealed that genes commonly have a high density of RNA polymerase II just downstream of the transcription start site. This has raised the possibility that genes are commonly regulated by transcriptional elongation, but this remains largely untested in vivo, particularly in vertebrates. Here, we show that the proximal promoter from the Rhox5 homeobox gene recruits polymerase II and begins elongating in all tissues and cell lines that we tested, but it only completes elongation in a tissue-specific and developmentally regulated manner. Relief of the elongation block is associated with recruitment of the elongation factor P-TEFb, the co-activator GRIP1, the chromatin remodeling factor BRG1, and specific histone modifications. We provide evidence that two mechanisms relieve the elongation block at the proximal promoter: demethylation and recruitment of androgen receptor. Together, our findings support a model in which promoter proximal pausing helps confer tissue-specific and developmental gene expression through a mechanism regulated by DNA demethylation-dependent nuclear hormone receptor recruitment.

Analytic approaches to stochastic gene expression in multicellular systems

Deterministic thermodynamic models of the complex systems, which control gene expression in metazoa, are helping researchers identify fundamental themes in the regulation of transcription. However, quantitative single cell studies are increasingly identifying regulatory mechanisms that control variability in expression. Such behaviors cannot be captured by deterministic models and are poorly suited to contemporary stochastic approaches that rely on continuum approximations, such as Langevin methods. Fortunately, theoretical advances in the modeling of transcription have assembled some general results that can be readily applied to systems being explored only through a deterministic approach. Here, I review some of the recent experimental evidence for the importance of genetically regulating stochastic effects during embryonic development and discuss key results from Markov theory that can be used to model this regulation. I then discuss several pairs of regulatory mechanisms recently investigated through a Markov approach. In each case, a deterministic treatment predicts no difference between the mechanisms, but the statistical treatment reveals the potential for substantially different distributions of transcriptional activity. In this light, features of gene regulation that seemed needlessly complex evolutionary baggage may be appreciated for their key contributions to reliability and precision of gene expression.

The embryo as a laboratory: quantifying transcription in Drosophila

Transcriptional regulation of gene expression is fundamental to most cellular processes, including determination of cellular fates. Quantitative studies of transcription in cultured cells have led to significant advances in identifying mechanisms underlying transcriptional control. Recent progress allowed implementation of these same quantitative methods in multicellular organisms to ask how transcriptional regulation unfolds both in vivo and at the single molecule level in the context of embryonic development. Here we review some of these advances in early Drosophila development, which bring the embryo on par with its single celled counterparts. In particular, we discuss progress in methods to measure mRNA and protein distributions in fixed and living embryos, and we highlight some initial applications that lead to fundamental new insights about molecular transcription processes. We end with an outlook on how to further exploit the unique advantages that come with investigating transcriptional control in the multicellular context of development.

Feedback control of gene expression variability in the Caenorhabditis elegans Wnt pathway

Variability in gene expression contributes to phenotypic heterogeneity even in isogenic populations. Here, we used the stereotyped, Wnt signaling-dependent development of the Caenorhabditis elegans Q neuroblast to probe endogenous mechanisms that control gene expression variability. We found that the key Hox gene that orients Q neuroblast migration exhibits increased gene expression variability in mutants in which Wnt pathway activity has been perturbed. Distinct features of the gene expression distributions prompted us on a systematic search for regulatory interactions, revealing a network of interlocked positive and negative feedback loops. Interestingly, positive feedback appeared to cooperate with negative feedback to reduce variability while keeping the Hox gene expression at elevated levels. A minimal model correctly predicts the increased gene expression variability across mutants. Our results highlight the influence of gene network architecture on expression variability and implicate feedback regulation as an effective mechanism to ensure developmental robustness.

Comparing mRNA levels using in situ hybridization of a target gene and co-stain

In situ hybridization is an important technique for measuring the spatial expression patterns of mRNA in cells, tissues, and whole animals. However, mRNA levels cannot be compared across experiments using typical protocols. Here we present a semi-quantitative method to compare mRNA levels of a gene across multiple samples. This method yields an estimate of the error in the measurement to allow statistical comparison. Our method uses a typical in situ hybridization protocol to stain for a target gene and an internal standard, which we refer to as a co-stain. As a proof of concept, we apply this method to multiple lines of transgenic Drosophila embryos, harboring constructs that express reporter genes to different levels. We generated this test set by mutating enhancer sequences to contain different numbers of binding sites for Zelda, a transcriptional activator. We demonstrate that using a co-stain with in situ hybridization is an effective method to compare mRNA levels across samples. This method requires only minor modifications to existing in situ hybridization protocols and uses straightforward analysis techniques. This strategy can be broadly applied to detect quantitative, spatially resolved changes in mRNA levels.

Genetic determinants and cellular constraints in noisy gene expression

In individual cells, transcription is a random process obeying single-molecule kinetics. Often, it occurs in a bursty, intermittent manner. The frequency and size of these bursts affect the magnitude of temporal fluctuations in messenger RNA and protein content within a cell, creating variation or "noise" in gene expression. It is still unclear to what degree transcriptional kinetics are specific to each gene and determined by its promoter sequence. Alternative scenarios have been proposed, in which the kinetics of transcription are governed by cellular constraints and follow universal rules across the genome. Evidence from genome-wide noise studies and from systematic perturbations of promoter sequences suggest that both scenarios-namely gene-specific versus genome-wide regulation of transcription kinetics-may be present to different degrees in bacteria, yeast, and animal cells.

Paused Pol II coordinates tissue morphogenesis in the Drosophila embryo

Paused RNA polymerase (Pol II) is a pervasive feature of Drosophila embryos and mammalian stem cells, but its role in development is uncertain. Here, we demonstrate that a spectrum of paused Pol II determines the "time to synchrony"-the time required to achieve coordinated gene expression across the cells of a tissue. To determine whether synchronous patterns of gene activation are significant in development, we manipulated the timing of snail expression, which controls the coordinated invagination of ∼1,000 mesoderm cells during gastrulation. Replacement of the strongly paused snail promoter with moderately paused or nonpaused promoters causes stochastic activation of snail expression and increased variability of mesoderm invagination. Computational modeling of the dorsal-ventral patterning network recapitulates these variable and bistable gastrulation profiles and emphasizes the importance of timing of gene activation in development. We conclude that paused Pol II and transcriptional synchrony are essential for coordinating cell behavior during morphogenesis.