Visualization of individual Scr mRNAs during Drosophila embryogenesis yields evidence for transcriptional bursting

Dynamic regulation of eve stripe 2 expression reveals transcriptional bursts in living Drosophila embryos

We present the use of recently developed live imaging methods to examine the dynamic regulation of even-skipped (eve) stripe 2 expression in the precellular Drosophila embryo. Nascent transcripts were visualized via MS2 RNA stem loops. The eve stripe 2 transgene exhibits a highly dynamic pattern of de novo transcription, beginning with a broad domain of expression during nuclear cycle 12 (nc12), and progressive refinement during nc13 and nc14. The mature stripe 2 pattern is surprisingly transient, constituting just ∼15 min of the ∼90-min period of expression. Nonetheless, this dynamic transcription profile faithfully predicts the limits of the mature stripe visualized by conventional in situ detection methods. Analysis of individual transcription foci reveals intermittent bursts of de novo transcription, with duration cycles of 4-10 min. We discuss a multistate model of transcription regulation and speculate on its role in the dynamic repression of the eve stripe 2 expression pattern during development.

Phylogenetic comparison of oskar mRNA localization signals

As a way to spatially control the expression of genes within cells, RNA localization is being recognized as an important process by which proteins are restricted to specific subcellular domains, which occurs in more diverse types of tissue than previously considered. Although many localized RNAs have been identified, information on cis-acting elements of localization is still limited. As transcripts of oskar (osk) are known to localize to the posterior pole of oocytes, we computationally analyzed a conserved sequence among eight Drosophila species and tested its role as a localization element. Dimerization of osk mRNA did not occur when the motif was deleted, but this did not affect assembly of osk mRNA-containing ribonucleoprotein (RNP) complexes. Without the motif, however, large RNP complex particles accumulated in nurse cells, and only a small fraction of these RNP complexes was transported into oocytes and properly localized to the posterior pole. Therefore, this motif may be required for the early transport of osk mRNA into oocytes. Also, as dimerization of osk mRNA does not seem to be a prerequisite for the assembly of RNP complexes, a dimerization-independent mechanism may also serve to localize osk mRNA to the posterior pole.

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.

Precise developmental gene expression arises from globally stochastic transcriptional activity

Early embryonic patterning events are strikingly precise, a fact that appears incompatible with the stochastic gene expression observed across phyla. Using single-molecule mRNA quantification in Drosophila embryos, we determine the magnitude of fluctuations in the expression of four critical patterning genes. The accumulation of mRNAs is identical across genes and fluctuates by only ∼8% between neighboring nuclei, generating precise protein distributions. In contrast, transcribing loci exhibit an intrinsic noise of ∼45% independent of specific promoter-enhancer architecture or fluctuating inputs. Precise transcript distribution in the syncytium is recovered via straightforward spatiotemporal averaging, i.e., accumulation and diffusion of transcripts during nuclear cycles, without regulatory feedback. Common expression characteristics shared between genes suggest that fluctuations in mRNA production are context independent and are a fundamental property of transcription. The findings shed light on how the apparent paradox between stochastic transcription and developmental precision is resolved.

Quantitative imaging of transcription in living Drosophila embryos links polymerase activity to patterning

Spatiotemporal patterns of gene expression are fundamental to every developmental program. The resulting macroscopic domains have been mainly characterized by their levels of gene products. However, the establishment of such patterns results from differences in the dynamics of microscopic events in individual cells such as transcription. It is unclear how these microscopic decisions lead to macroscopic patterns, as measurements in fixed tissue cannot access the underlying transcriptional dynamics. In vivo transcriptional dynamics have long been approached in single-celled organisms, but never in a multicellular developmental context. Here, we directly address how boundaries of gene expression emerge in the Drosophila embryo by measuring the absolute number of actively transcribing polymerases in real time in individual nuclei. Specifically, we show that the formation of a boundary cannot be quantitatively explained by the rate of mRNA production in each cell, but instead requires amplification of the dynamic range of the expression boundary. This amplification is accomplished by nuclei randomly adopting active or inactive states of transcription, leading to a collective effect where the fraction of active nuclei is modulated in space. Thus, developmental patterns are not just the consequence of reproducible transcriptional dynamics in individual nuclei, but are the result of averaging expression over space and time.

Quantitative assessment of ratiometric bimolecular beacons as a tool for imaging single engineered RNA transcripts and measuring gene expression in living cells

Recently, we developed an oligonucleotide-based probe, ratiometric bimolecular beacon (RBMB), which generates a detectable fluorescent signal in living cells that express the target RNA. Here, we show that RBMBs can also be used to image single RNA transcripts in living cells, when the target RNA is engineered to contain as few as four hybridization sites. Moreover, comparison with single-molecule fluorescence in situ hybridization confirmed that RBMBs could be used to accurately quantify the number of RNA transcripts within individual cells. Measurements of gene expression could be acquired within 30 min and using a wide range of RBMB concentrations. The ability to acquire accurate measurements of RNA copy number in both HT-1080 cells and CHO cells also suggests that RBMBs can be used to image and quantify single RNA transcripts in a wide range of cell lines. Overall, these findings highlight the robustness and versatility of RBMBs as a tool for imaging RNA in live cells. We envision that the unique capabilities of RBMBs will open up new avenues for RNA research.

Are gene loops the cause of transcriptional noise?

Expression levels of the same mRNA or protein vary significantly among the cells of an otherwise identical population. Such biological noise has great functional implications and is largely due to transcriptional bursting, the episodic production of mRNAs in short, intense bursts, interspersed by periods of transcriptional inactivity. Bursting has been demonstrated in a wide range of pro- and eukaryotic species, attesting to its universal importance. However, the mechanistic origins of bursting remain elusive. A different type of phenomenon, which has also been suggested to be widespread, is the physical interaction between the promoter and 3' end of a gene. Several functional roles have been proposed for such gene loops, including the facilitation of transcriptional reinitiation. Here, I discuss the most recent findings related to these subjects and argue that gene loops are a likely cause of transcriptional bursting and, thus, biological noise.

Transcription stochasticity of complex gene regulation models

Transcription is regulated by a multitude of factors that concertedly induce genes to switch between activity states. Eukaryotic transcription involves a multitude of complexes that sequentially assemble on chromatin under the influence of transcription factors and the dynamic state of chromatin. Prokaryotic transcription depends on transcription factors, sigma-factors, and, in some cases, on DNA looping. We present a stochastic model of transcription that considers these complex regulatory mechanisms. We coarse-grain the molecular details in such a way that the model can describe a broad class of gene-regulation mechanisms. We solve this model analytically for various measures of stochastic transcription and compare alternative gene-regulation designs. We find that genes with complex multiprotein regulation can have peaked burst-size distributions in contrast to the geometric distributions found for simple models of transcription regulation. Burst-size distributions are, in addition, shaped by mRNA degradation during transcription bursts. We derive the stochastic properties of genes in the limit of deterministic switch times. These genes typically have reduced transcription noise. Severe timescale separation between gene regulation and transcription initiation enhances noise and leads to bimodal mRNA copy number distributions. In general, complex mechanisms for gene regulation lead to nonexponential waiting-time distributions for gene switching and transcription initiation, which typically reduce noise in mRNA copy numbers and burst size. Finally, we discuss that qualitatively different gene regulation models can often fit the same experimental data on single-cell mRNA abundance even though they have qualitatively different burst-size statistics and regulatory parameters.

Rapid transcription fosters coordinate snail expression in the Drosophila embryo

Transcription is commonly held to be a highly stochastic process, resulting in considerable heterogeneity of gene expression among the different cells in a population. Here, we employ quantitative in situ hybridization methods coupled with high-resolution imaging assays to measure the expression of snail, a developmental patterning gene necessary for coordinating the invagination of the mesoderm during gastrulation of the Drosophila embryo. Our measurements of steady-state mRNAs suggest that there is very little variation in snail expression across the different cells that make up the mesoderm and that synthesis approaches the kinetic limits of Pol II processivity. We propose that rapid transcription kinetics and negative autoregulation are responsible for the remarkable homogeneity of snail expression and the coordination of mesoderm invagination.

Quantitative biology of single neurons

The building blocks of complex biological systems are single cells. Fundamental insights gained from single-cell analysis promise to provide the framework for understanding normal biological systems development as well as the limits on systems/cellular ability to respond to disease. The interplay of cells to create functional systems is not well understood. Until recently, the study of single cells has concentrated primarily on morphological and physiological characterization. With the application of new highly sensitive molecular and genomic technologies, the quantitative biochemistry of single cells is now accessible.

Live imaging of nascent RNA dynamics reveals distinct types of transcriptional pulse regulation

Transcription of genes can be discontinuous, occurring in pulses or bursts. It is not clear how properties of transcriptional pulses vary between different genes. We compared the pulsing of five housekeeping and five developmentally induced genes by direct imaging of single gene transcriptional events in individual living Dictyostelium cells. Each gene displayed its own transcriptional signature, differing in probability of firing and pulse duration, frequency, and intensity. In contrast to the prevailing view from both prokaryotes and eukaryotes that transcription displays binary behavior, strongly expressed housekeeping genes altered the magnitude of their transcriptional pulses during development. These nonbinary "tunable" responses may be better suited than stochastic switch behavior for housekeeping functions. Analysis of RNA synthesis kinetics using fluorescence recovery after photobleaching implied modulation of housekeeping-gene pulse strength occurs at the level of transcription initiation rather than elongation. In addition, disparities between single cell and population measures of transcript production suggested differences in RNA stability between gene classes. Analysis of stability using RNAseq revealed no major global differences in stability between developmental and housekeeping transcripts, although strongly induced RNAs showed unusually rapid decay, indicating tight regulation of expression.

Transcription goes digital

Transcription is a complex process that integrates the state of the cell and its environment to generate adequate responses for cell fitness and survival. Recent microscopy experiments have been able to monitor transcription from single genes in individual cells. These observations have revealed two striking features: transcriptional activity can vary markedly from one cell to another, and is subject to large changes over time, sometimes within minutes. How the chromatin structure, transcription machinery assembly and signalling networks generate such patterns is still unclear. In this review, we present the techniques used to investigate transcription from single genes, introduce quantitative modelling tools, and discuss transcription mechanisms and their implications for gene expression regulation.

Three Drosophila Hox complex microRNAs do not have major effects on expression of evolutionarily conserved Hox gene targets during embryogenesis

The discovery of microRNAs has resulted in a major expansion of the number of molecules known to be involved in gene regulation. Elucidating the functions of animal microRNAs has posed a significant challenge as their target interactions with messenger RNAs do not adhere to simple rules. Of the thousands of known animal microRNAs, relatively few microRNA:messenger RNA regulatory interactions have been biologically validated in an normal organismal context. Here we present evidence that three microRNAs from the Hox complex in Drosophila (miR-10-5p, miR-10-3p, miR-iab-4-5p) do not have significant effects during embryogenesis on the expression of Hox genes that contain high confidence microRNAs target sites in the 3′ untranslated regions of their messenger RNAs. This is significant, in that it suggests that many predicted microRNA-target interactions may not be biologically relevant, or that the outcomes of these interactions may be so subtle that mutants may only show phenotypes in specific contexts, such as in environmental stress conditions, or in combinations with other microRNA mutations.

Single-molecule assays reveal that RNA localization signals regulate dynein-dynactin copy number on individual transcript cargoes

Subcellular localization of mRNAs by cytoskeletal motors plays critical roles in the spatial control of protein function¹. However, optical limitations of studying mRNA transport in vivo mean that there is little mechanistic insight into how transcripts are packaged and linked to motors, and how the movement of mRNA:motor complexes on the cytoskeleton is orchestrated. Here, we have reconstituted transport of mRNPs containing specific RNAs in vitro. We show directly that mRNAs that are either apically localized or non-localized in Drosophila embryos associate with the dynein motor and move bidirectionally on individual microtubules, with localizing mRNPs exhibiting a strong minus-end-directed bias. Single-molecule fluorescence measurements reveal that RNA localization signals increase the average number of dynein and dynactin components recruited to individual mRNPs. We find that, surprisingly, individual RNA molecules are present in motile mRNPs in vitro and present evidence that this is also the case in vivo. Thus, RNA oligomerization is not obligatory for transport. Our findings lead to a model in which RNA localization signals produce highly polarized distributions of transcript populations through modest changes in motor copy number on single mRNA molecules.

Evaluating the Drosophila Bicoid morphogen gradient system through dissecting the noise in transcriptional bursts

MOTIVATION

We describe a statistical model to dissect the noise in transcriptional bursts in a developmental system.

RESULTS

We assume that, at any given moment of time, each copy of a native gene inside a cell can exist in either a bursting (active) or non-bursting (inactive) state. The experimentally measured total noise in the transcriptional states of a gene in a population of cells can be mathematically dissected into two contributing components: internal and external. While internal noise quantifies the stochastic nature of transcriptional bursts, external noise is caused by cell-to-cell differences including fluctuations in activator concentration. We use our developed methods to analyze the Drosophila Bicoid (Bcd) morphogen gradient system. For its target gene hunchback (hb), the noise properties can be recapitulated by a simplified gene regulatory model in which Bcd acts as the only input, suggesting that the external noise in hb transcription is primarily derived from fluctuations in the Bcd activator input. However, such a simplified gene regulatory model is insufficient to predict the noise properties of another Bcd target gene, orthodenticle (otd), suggesting that otd transcription is sensitive to additional external fluctuations beyond those in Bcd. Our results show that analysis of the relationship between input and output noise can reveal important insights into how a morphogen gradient system works. Our study also advances the knowledge about transcription at a fundamental level.

CONTACT

jun.ma@cchmc.org

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Transcriptional activators and activation mechanisms

Transcriptional activators are required to turn on the expression of genes in a eukaryotic cell. Activators bound to the enhancer can facilitate either the recruitment of RNA polymerase II to the promoter or its elongation. This article examines a few selected issues in understanding activator functions and activation mechanisms.

Pulsatile patterns of pituitary hormone gene expression change during development

Important questions in biology have emerged recently concerning the timing of transcription in living cells. Studies on clonal cell lines have shown that transcription is often pulsatile and stochastic, with implications for cellular differentiation. Currently, information regarding transcriptional activity at cellular resolution within a physiological context remains limited. To investigate single-cell transcriptional activity in real-time in living tissue we used bioluminescence imaging of pituitary tissue from transgenic rats in which luciferase gene expression is driven by a pituitary hormone gene promoter. We studied fetal and neonatal pituitary tissue to assess whether dynamic patterns of transcription change during tissue development. We show that gene expression in single cells is highly pulsatile at the time endocrine cells first appear but becomes stabilised as the tissue develops in early neonatal life. This stabilised transcription pattern might depend upon tissue architecture or paracrine signalling, as isolated cells, generated from enzymatic dispersion of the tissue, display pulsatile luminescence. Nascent cells in embryonic tissue also showed coordinated transcription activity over short distances further indicating that cellular context is important for transcription activity. Overall, our data show that cells alter their patterns of gene expression according to their context and developmental stage, with important implications for cellular differentiation.

Nuclear organization and transcriptional dynamics in Dictyostelium

The Dictyostelium model has a set of features uniquely well-suited to developing our understanding of transcriptional control. The complete Dictyostelium discoideum genome sequence has revealed that many of the molecular components regulating transcription in larger eukaryotes are conserved in Dictyostelium, from transcription factors and chromatin components to the enzymes and signals that regulate them. In addition, the system permits visualization of single gene firing events in living cells, which provides a more detailed view of transcription and its relationships to cell and developmental processes. This review will bring together the available knowledge of the structure and dynamics of the Dictyostelium nucleus and discuss recent transcription imaging studies and their implications for stability and accuracy of cell decisions.

The snail repressor inhibits release, not elongation, of paused Pol II in the Drosophila embryo

The development of the precellular Drosophila embryo is characterized by exceptionally rapid transitions in gene activity, with broadly distributed maternal regulatory gradients giving way to precise on/off patterns of gene expression within a one-hour window, between two and three hours after fertilization [1]. Transcriptional repression plays a pivotal role in this process, delineating sharp expression patterns (e.g., pair-rule stripes) within broad domains of gene activation. As many as 20 different sequence-specific repressors have been implicated in this process, yet the mechanisms by which they silence gene expression have remained elusive [2]. Here we report the development of a method for the quantitative visualization of transcriptional repression. We focus on the Snail repressor, which establishes the boundary between the presumptive mesoderm and neurogenic ectoderm [3]. We find that elongating Pol II complexes complete transcription after the onset of Snail repression. As a result, moderately sized genes (e.g., the 22 kb sog locus) are fully silenced only after tens of minutes of repression. We propose that this "repression lag" imposes a severe constraint on the regulatory dynamics of embryonic patterning and further suggest that posttranscriptional regulators, like microRNAs, are required to inhibit unwanted transcripts produced during protracted periods of gene silencing.

Pabp binds to the osk 3'UTR and specifically contributes to osk mRNA stability and oocyte accumulation

RNA localization is tightly coordinated with RNA stability and translation control. Bicaudal-D (Bic-D), Egalitarian (Egl), microtubules and their motors are part of a Drosophila transport machinery that localizes mRNAs to specific cellular regions during oogenesis and embryogenesis. We identified the Poly(A)-binding protein (Pabp) as a protein that forms an RNA-dependent complex with Bic-D in embryos and ovaries. pabp also interacts genetically with Bic-D and, similar to Bic-D, pabp is essential in the germline for oocyte growth and accumulation of osk mRNA in the oocyte. In the absence of pabp, reduced stability of osk mRNA and possibly also defects in osk mRNA transport prevent normal oocyte localization of osk mRNA. pabp also interacts genetically with osk and lack of one copy of pabp(+) causes osk to become haploinsufficient. Moreover, pointing to a poly(A)-independent role, Pabp binds to A-rich sequences (ARS) in the osk 3'UTR and these turned out to be required in vivo for osk function during early oogenesis. This effect of pabp on osk mRNA is specific for this RNA and other tested mRNAs localizing to the oocyte are less and more indirectly affected by the lack of pabp.

Uniform distribution of elongating RNA polymerase II complexes in transcribed gene locus

The intensity of gene transcription is generally reflected by the level of RNA polymerase II (RNAPII) recruitment to the gene. However, genome-wide studies of polymerase occupancy indicate that RNAPII distribution varies among genes. In some loci more polymerases are found in the 5′ region, whereas in other loci, in the 3′ region of the gene. We studied the distribution of elongating RNAPII complexes at highly transcribed GAL-VPS13 locus in Saccharomyces cerevisiae and found that in the cell population the amount of polymerases gradually decreased toward the 3′ end of the gene. However, the conventional chromatin immunoprecipitation assay averages the signal from the cell population, and no data on single cell level can be gathered. To study the spacing of elongating polymerases on single chromosomes, we used a sequential chromatin immunoprecipitation assay for the detection of multiple RNAPII complexes on the same DNA fragment. Our results demonstrate uniform distribution of elongating polymerases throughout all regions of the GAL-VPS13 gene.

General properties of transcriptional time series in Escherichia coli

Gene activity is described by the time-series of discrete, stochastic mRNA production events. This transcriptional time-series exhibits intermittent, bursty behavior. One consequence of this temporal intricacy is that gene expression can be tuned by varying different features of the time-series. What schemes for varying the transcriptional time-series are observed in the cell? Are the observed properties of these time-series optimized for cellular function? To address these questions, we characterize mRNA copy-number statistics at single-molecule resolution from multiple Escherichia coli promoters. We find that the degree of burstiness depends only on the gene expression level, while being independent of the details of gene regulation. The observed behavior is explained by the underlying variation in the duration of bursting events. Using information theory, we find that the properties of the transcriptional time series allow the cell to efficiently map the extracellular concentration of inducer molecules to intracellular levels of mRNA and proteins.

A multiscale investigation of bicoid-dependent transcriptional events in Drosophila embryos

Background

Morphogen molecules form concentration gradients to provide spatial information to cells in a developing embryo. Precisely how cells decode such information to form patterns with sharp boundaries remains an open question. For example, it remains controversial whether the Drosophila morphogenetic protein Bicoid (Bcd) plays a transient or sustained role in activating its target genes to establish sharp expression boundaries during development.

Methodology/Principal Findings

In this study, we describe a method to simultaneously detect Bcd and the nascent transcripts of its target genes in developing embryos. This method allows us to investigate the relationship between Bcd and the transcriptional status of individual copies of its target genes on distinct scales. We show that, on three scales analyzed concurrently—embryonic, nuclear and local, the actively-transcribing gene copies are associated with high Bcd concentrations. These results underscore the importance of Bcd as a sustained input for transcriptional decisions of individual copies of its target genes during development. We also show that the Bcd-dependent transcriptional decisions have a significantly higher noise than Bcd-dependent gene products, suggesting that, consistent with theoretical studies, time and/or space averaging reduces the noise of Bcd-activated transcriptional output. Finally, our analysis of an X-linked Bcd target gene reveals that Bcd-dependent transcription bursts at twice the frequency in males as in females, providing a mechanism for dosage compensation in early Drosophila embryos.

Conclusion/Significance

Our study represents a first experimental uncovering of the actions of Bcd in controlling the actual transcriptional events while its positional information is decoded during development. It establishes a sustained role of Bcd in transcriptional decisions of individual copies of its target genes to generate sharp expression boundaries. It also provides an experimental evaluation of the effect of time and/or space averaging on Bcd-dependent transcriptional output, and establishes a dosage compensation mechanism in early Drosophila embryos.

Validating transcripts with probes and imaging technology

High-throughput gene expression screens provide a quantitative picture of the average expression signature of biological samples. However, the analysis of spatial gene expression patterns with single-cell resolution requires quantitative in situ measurement techniques. Here we describe recent technological advances in RNA fluorescence in situ hybridization (FISH) techniques that facilitate detection of individual fluorescently labeled mRNA molecules of practically any endogenous gene. These methods, which are based on advances in probe design, imaging technology and image processing, enable the absolute measurement of transcript abundance in individual cells with single-molecule resolution.

Mammalian genes are transcribed with widely different bursting kinetics

In prokaryotes and eukaryotes, most genes appear to be transcribed during short periods called transcriptional bursts, interspersed by silent intervals. We describe how such bursts generate gene-specific temporal patterns of messenger RNA (mRNA) synthesis in mammalian cells. To monitor transcription at high temporal resolution, we established various gene trap cell lines and transgenic cell lines expressing a short-lived luciferase protein from an unstable mRNA, and recorded bioluminescence in real time in single cells. Mathematical modeling identified gene-specific on- and off-switching rates in transcriptional activity and mean numbers of mRNAs produced during the bursts. Transcriptional kinetics were markedly altered by cis-regulatory DNA elements. Our analysis demonstrated that bursting kinetics are highly gene-specific, reflecting refractory periods during which genes stay inactive for a certain time before switching on again.

The formation of the Bicoid morphogen gradient requires protein movement from anteriorly localized mRNA

New quantitative data show that the Bicoid morphogen gradient is generated from a dynamic localized source and that protein gradient formation requires protein movement along the anterior-posterior axis.

The Bicoid morphogen gradient directs the patterning of cell fates along the anterior-posterior axis of the syncytial Drosophila embryo and serves as a paradigm of morphogen-mediated patterning. The simplest models of gradient formation rely on constant protein synthesis and diffusion from anteriorly localized source mRNA, coupled with uniform protein degradation. However, currently such models cannot account for all known gradient characteristics. Recent work has proposed that bicoid mRNA spatial distribution is sufficient to produce the observed protein gradient, minimizing the role of protein transport. Here, we adapt a novel method of fluorescent in situ hybridization to quantify the global spatio-temporal dynamics of bicoid mRNA particles. We determine that >90% of all bicoid mRNA is continuously present within the anterior 20% of the embryo. bicoid mRNA distribution along the body axis remains nearly unchanged despite dynamic mRNA translocation from the embryo core to the cortex. To evaluate the impact of mRNA distribution on protein gradient dynamics, we provide detailed quantitative measurements of nuclear Bicoid levels during the formation of the protein gradient. We find that gradient establishment begins 45 minutes after fertilization and that the gradient requires about 50 minutes to reach peak levels. In numerical simulations of gradient formation, we find that incorporating the actual bicoid mRNA distribution yields a closer prediction of the observed protein dynamics compared to modeling protein production from a point source at the anterior pole. We conclude that the spatial distribution of bicoid mRNA contributes to, but cannot account for, protein gradient formation, and therefore that protein movement, either active or passive, is required for gradient formation.

The Bicoid protein gradient plays a crucial role in determining the anterior body pattern of Drosophila embryos. This gradient is the classic example of morphogen-mediated patterning of a developing metazoan and serves as a major topic for mathematical modeling. Accurate modeling of the gradient requires a detailed account of the underlying bicoid mRNA distribution. The classic model holds that mRNA protein gradient arises via protein diffusion from mRNA localized at the anterior of the developing egg. In contrast, recent proposals suggest that an mRNA gradient generates the protein gradient without protein movement. In this study, we introduce a novel mRNA quantification method for Drosophila embryos, which allows us to visualize each individual mRNA particle accurately in whole embryos. We demonstrate that all but a few mRNA particles are confined to the anterior 20% of the egg, and consequently that the protein must move in order to establish a gradient. We further report that the mRNA distribution is highly dynamic during the time of protein synthesis. In numerical simulations, we show that incorporating realistic spatial locations of the individual source mRNA molecules throughout the developmental period is necessary to accurately model the experimentally observed protein gradient dynamics.

Experimental approaches to evaluate the contributions of candidate cis-regulatory mutations to phenotypic evolution

Elucidating the molecular bases by which phenotypic traits have evolved provides a glimpse into the past, allowing the characterization of genetic changes that cumulatively contribute to evolutionary innovations. Historically, much of the experimental attention has been focused on changes in protein-coding regions that can readily be identified by the genetic code for translating gene coding sequences into proteins. Resultantly, the role of noncoding sequences in trait evolution has remained more mysterious. In recent years, several studies have reached an unprecedented level of detail in describing how noncoding mutations in gene cis-regulatory elements contribute to morphological evolution. Based on these and other studies, we describe an experimental framework and some of the genetic and molecular methods to connect a particular cis-regulatory mutation to the evolution of any phenotypic trait.

Identifying fluorescently labeled single molecules in image stacks using machine learning

In the past several years, a host of new technologies have made it possible to visualize single molecules within cells and organisms (Raj et al., Nat Methods 5:877-879, 2008; Paré et al., Curr Biol 19:2037-2042, 2009; Lu and Tsourkas, Nucleic Acids Res 37:e100, 2009; Femino et al., Science 280:585-590, 1998; Rodriguez et al., Semin Cell Dev Biol 18:202-208, 2007; Betzig et al., Science 313:1642-1645, 2006; Rust et al., Nat Methods 3:793-796, 2006; Fusco et al., Curr Biol 13:161-167, 2003). Many of these are based on fluorescence, either fluorescent proteins or fluorescent dyes coupled to a molecule of interest. In many applications, the fluorescent signal is limited to a few pixels, which poses a classic signal processing problem: how can actual signal be distinguished from background noise? In this chapter, I present a MATLAB (MathWorks (2010) MATLAB. Retrieved from http://www.mathworks.com) software suite designed to work with these single-molecule visualization technologies (Rifkin (2010) spotFinding Suite. http://www.biology.ucsd.edu/labs/rifkin/software.html). It takes images or image stacks from a fluorescence microscope as input and outputs locations of the molecules. Although the software was developed for the specific application of identifying single mRNA transcripts in fixed specimens, it is more general than this and can be used and/or customized for other applications that produce localized signals embedded in a potentially noisy background. The analysis pipeline consists of the following steps: (a) create a gold-standard dataset, (b) train a machine-learning algorithm to classify image features as signal or noise depending upon user defined statistics, (c) run the machine-learning algorithm on a new dataset to identify mRNA locations, and (d) visually inspect and correct the results.

Functional roles for noise in genetic circuits

The genetic circuits that regulate cellular functions are subject to stochastic fluctuations, or 'noise', in the levels of their components. Noise, far from just a nuisance, has begun to be appreciated for its essential role in key cellular activities. Noise functions in both microbial and eukaryotic cells, in multicellular development, and in evolution. It enables coordination of gene expression across large regulons, as well as probabilistic differentiation strategies that function across cell populations. At the longest timescales, noise may facilitate evolutionary transitions. Here we review examples and emerging principles that connect noise, the architecture of the gene circuits in which it is present, and the biological functions it enables. We further indicate some of the important challenges and opportunities going forward.

Inferring the ancestral function of the posterior Hox gene within the bilateria: controlling the maintenance of reproductive structures, the musculature and the nervous system in the acoel flatworm Isodiametra pulchra

Molecular phylogenies place the acoel flatworms as the sister-group to the remaining Bilateria, a position that should prove very valuable when trying to understand the evolutionary origins of the bilaterian body plan. A major feature characterizing Bilateria is the presence of two, orthogonal, body axis. In this article we aim at tackling the problem of how the bilaterian anterior-posterior (AP) axis is organized, and how this axis have been established over evolutionary time. To this purpose we have studied the role of some key regulatory genes involved in the control of the AP axis, the Hox family of transcription factors. All acoels studied to date contain a minimal complement of three Hox genes that are all expressed in nested domains along this major axis, providing the oldest evidence for a "Hox vectorial system" working in Bilateria. However, this proposition is not based in the analysis of Hox functions. Here we document the specific roles of one posterior Hox gene, IpHoxPost, in the postembryonic development of the acoel Isodiametra pulchra. The analysis has been done using RNA interference technologies, for the first time in acoels, and we demonstrate that the functions of this gene are restricted to the posterior region of the animal, within the muscular and neural tissues. We conclude, therefore, that the posterior Hox genes were used to specify and maintain defined anatomical regions within the AP axis of animals since the beginning of bilaterian evolution.

Bursts and pulses: insights from single cell studies into transcriptional mechanisms

With a developing appreciation of how noisy gene expression can be, and difficulties in deciphering conventional gene expression data into cell control mechanisms, it has become clear that single cell techniques for measuring transcription are necessary to illuminate basic cell regulation strategies. The resultant use of in situ hybridisation and live cell RNA visualisation approaches in single cells revealed transcription is not adequately reflected by the smooth, seamless process we tend to infer from standard measures of RNA level. When RNA production is measured in single cells, the process of transcription has been shown to occur in bursts, or pulses. This review will highlight the evidence for these phenomena, the proposed mechanisms underlying discontinuity, and the biological implications of such behaviour.