Synchronous and stochastic patterns of gene activation in the Drosophila embryo

The transcriptional and epigenomic foundations of ground state pluripotency

Mouse embryonic stem (ES) cells grown in serum exhibit greater heterogeneity in morphology and expression of pluripotency factors than ES cells cultured in defined medium with inhibitors of two kinases (Mek and GSK3), a condition known as "2i" postulated to establish a naive ground state. We show that the transcriptome and epigenome profiles of serum- and 2i-grown ES cells are distinct. 2i-treated cells exhibit lower expression of lineage-affiliated genes, reduced prevalence at promoters of the repressive histone modification H3K27me3, and fewer bivalent domains, which are thought to mark genes poised for either up- or downregulation. Nonetheless, serum- and 2i-grown ES cells have similar differentiation potential. Precocious transcription of developmental genes in 2i is restrained by RNA polymerase II promoter-proximal pausing. These findings suggest that transcriptional potentiation and a permissive chromatin context characterize the ground state and that exit from it may not require a metastable intermediate or multilineage priming.

Bivalent histone modifications in early embryogenesis

Histone modifications influence the interactions of transcriptional regulators with chromatin. Studies in embryos and embryonic stem (ES) cells have uncovered histone modification patterns that are diagnostic for different cell types and developmental stages. For example, bivalent domains consisting of regions of H3 lysine 27 trimethylation (H3K27me3) and H3 lysine 4 trimethylation (H3K4me3) mark lineage control genes in ES cells and zebrafish blastomeres. Such bivalent domains have garnered attention because the H3K27me3 mark might help repress lineage-regulatory genes during pluripotency while the H3K4me3 mark could poise genes for activation upon differentiation. Despite the prominence of the bivalent domain concept, studies in other model organisms have questioned its universal nature, and the function of bivalent domains has remained unclear. Histone marks are also associated with developmental regulatory genes in sperm. These observations have raised the possibility that specific histone modification patterns might persist from parent to offspring, but it is unclear whether histone marks are inherited or formed de novo. Here, we review the potential roles of H3K4me3 and H3K27me3 marks in embryos and ES cells and discuss how histone marks might be established, maintained and resolved during embryonic development.

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.

Functional association of Gdown1 with RNA polymerase II poised on human genes

Most human genes are loaded with promoter-proximally paused RNA polymerase II (Pol II) molecules that are poised for release into productive elongation by P-TEFb. We present evidence that Gdown1, the product of the POLR2M gene that renders Pol II responsive to Mediator, is involved in Pol II elongation control. During in vitro transcription, Gdown1 specifically blocked elongation stimulation by TFIIF, inhibited the termination activity of TTF2, and influenced pausing factors NELF and DSIF, but did not affect the function of TFIIS or the mRNA capping enzyme. Without P-TEFb, Gdown1 led to the production of stably paused polymerases in the presence of nuclear extract. Supporting these mechanistic insights, ChIP-Seq demonstrated that Gdown1 mapped over essentially all poised polymerases across the human genome. Our results establish that Gdown1 stabilizes poised polymerases while maintaining their responsiveness to P-TEFb and suggest that Mediator overcomes a Gdown1-mediated block of initiation by allowing TFIIF function.

Analytical distribution and tunability of noise in a model of promoter progress

Chromatin template (CT), which accumulates over time until the promoter becomes active, determines upstream dynamics of transcription, but how upstream sequential steps impact downstream dynamics qualitatively and quantitatively is unclear. Here, we analyze a stochastic gene model with a simple yet typical CT that contains one active state and several inactive states of the promoter. We derive the analytical expressions for the noise in mRNA probability distributions governed by master equations. The derived results extend previous work by including the effects of promoter progress on variability and bimodality. Specifically, given a CT for transcription, we analytically demonstrate that inactive phases of the promoter can modulate the noise intensity to the minimum independently of the mean expression of mRNA. If one new inactive state is added to the CT, then the resulting noise will be reduced, implying that the multi-off mechanism plays a role of attenuating the noise. In contrast to the simple on-off mechanism, the multi-off mechanism can also narrow bimodal regions in a certain parameter plane and obscure two peaks, explaining why bimodal distributions are rarely observed in experiments. Our results provide insight into the role of promoter progress in determining the level of cell-to-cell variability in gene expression.

Multiple exposures to drought 'train' transcriptional responses in Arabidopsis

Pre-exposure to stress may alter plants' subsequent responses by producing faster and/or stronger reactions implying that plants exercise a form of 'stress memory'. The mechanisms of plants' stress memory responses are poorly understood leaving this fundamental biological question unanswered. Here we show that during recurring dehydration stresses Arabidopsis plants display transcriptional stress memory demonstrated by an increase in the rate of transcription and elevated transcript levels of a subset of the stress-response genes (trainable genes). During recovery (watered) states, trainable genes produce transcripts at basal (preinduced) levels, but remain associated with atypically high H3K4me3 and Ser5P polymerase II levels, indicating that RNA polymerase II is stalled. This is the first example of a stalled RNA polymerase II and its involvement in transcriptional memory in plants. These newly discovered phenomena might be a general feature of plant stress-response systems and could lead to novel approaches for increasing the flexibility of a plant's ability to respond to the environment.

The Drosophila 7SK snRNP and the essential role of dHEXIM in development

Regulation of the positive transcription elongation factor, P-TEFb, plays a major role in controlling mammalian transcription and this is accomplished in part by controlled release of P-TEFb from the 7SK snRNP that sequesters the kinase in an inactive state. We demonstrate here that a similar P-TEFb control system exists in Drosophila. We show that an RNA previously suggested to be a 7SK homolog is, in fact, associated with P-TEFb, through the action of a homolog of the human HEXIM1/2 proteins (dHEXIM). In addition, a Drosophila La related protein (now called dLARP7) is shown to be the functional homolog of human LARP7. The Drosophila 7SK snRNP (d7SK snRNP) responded to treatment of cells with P-TEFb inhibitors and to nuclease treatment of cell lysates by releasing P-TEFb. Supporting a critical role for the d7SK snRNP in Drosophila development, dLARP7 and dHEXIM were found to be ubiquitously expressed throughout embryos and tissues at all stages. Importantly, knockdown of dHEXIM was embryonic lethal, and reduction of dHEXIM in specific tissues led to serious developmental defects. Our results suggest that regulation of P-TEFb by the d7SK snRNP is essential for the growth and differentiation of tissues required during Drosophila development.

Dorsal-ventral gene expression in the Drosophila embryo reflects the dynamics and precision of the dorsal nuclear gradient

Patterning of the dorsal-ventral axis in the early Drosophila embryo depends on the nuclear distribution of the Dorsal transcription factor. Using live two-photon light-sheet microscopy, we quantified the nuclear Dorsal gradient in space and time and found that its amplitude and basal levels display oscillations throughout early embryonic development. These dynamics raise questions regarding how cells can reproducibly establish patterns of gene expression from a rapidly varying signal. We therefore quantified domains of Dorsal target genes, discovering their expression patterns are also dynamic. Computational modeling of this system reveals a correlation between Dorsal gradient dynamics and changes in target gene expression and suggests that these dynamics, together with time averaging of noise, results in the formation of graded gene expression borders in regions where the gradient is nearly flat. We propose that mRNA levels remain plastic during transient signaling events, allowing tissues to refine patterns in the face of genetic or environmental variation.

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.

Consequences of eukaryotic enhancer architecture for gene expression dynamics, development, and fitness

The regulatory logic of time- and tissue-specific gene expression has mostly been dissected in the context of the smallest DNA fragments that, when isolated, recapitulate native expression in reporter assays. It is not known if the genomic sequences surrounding such fragments, often evolutionarily conserved, have any biological function or not. Using an enhancer of the even-skipped gene of Drosophila as a model, we investigate the functional significance of the genomic sequences surrounding empirically identified enhancers. A 480 bp long "minimal stripe element" is able to drive even-skipped expression in the second of seven stripes but is embedded in a larger region of 800 bp containing evolutionarily conserved binding sites for required transcription factors. To assess the overall fitness contribution made by these binding sites in the native genomic context, we employed a gene-replacement strategy in which whole-locus transgenes, capable of rescuing even-skipped(-) lethality to adulthood, were substituted for the native gene. The molecular phenotypes were characterized by tagging Even-skipped with a fluorescent protein and monitoring gene expression dynamics in living embryos. We used recombineering to excise the sequences surrounding the minimal enhancer and site-specific transgenesis to create co-isogenic strains differing only in their stripe 2 sequences. Remarkably, the flanking sequences were dispensable for viability, proving the sufficiency of the minimal element for biological function under normal conditions. These sequences are required for robustness to genetic and environmental perturbation instead. The mutant enhancers had measurable sex- and dose-dependent effects on viability. At the molecular level, the mutants showed a destabilization of stripe placement and improper activation of downstream genes. Finally, we demonstrate through live measurements that the peripheral sequences are required for temperature compensation. These results imply that seemingly redundant regulatory sequences beyond the minimal enhancer are necessary for robust gene expression and that "robustness" itself must be an evolved characteristic of the wild-type enhancer.

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.

Dynamic transcriptional events in embryonic stem cells mediated by the super elongation complex (SEC)

Transcriptional regulation of developmentally controlled genes is at the heart of differentiation and organogenesis. In this study, we performed global genomic analyses in murine embryonic stem (ES) cells and in human cells in response to activation signals. We identified an essential role for the ELL (eleven-nineteen lysine-rich leukemia gene)/P-TEFb (positive transcription elongation factor)-containing super elongation complex (SEC) in the regulation of gene expression, including several genes bearing paused RNA polymerase II (Pol II). Paused Pol II has been proposed to be associated with loci that respond rapidly to environmental stimuli. However, our studies in ES cells also identified a requirement for SEC at genes without paused Pol II, which also respond dynamically to differentiation signals. Our findings suggest that SEC is a major class of active P-TEFb-containing complexes required for transcriptional activation in response to environmental cues such as differentiation signals.

Gene length may contribute to graded transcriptional responses in the Drosophila embryo

An important question in developmental biology is how relatively shallow gradients of morphogens can reliably establish a series of distinct transcriptional readouts. Current models emphasize interactions between transcription factors binding in distinct modes to cis-acting sequences of target genes. Another recent idea is that the cis-acting interactions may amplify preexisting biases or prepatterns to establish robust transcriptional responses. In this study, we examine the possible contribution of one such source of prepattern, namely gene length. We developed quantitative imaging tools to measure gene expression levels for several loci at a time on a single-cell basis and applied these quantitative imaging tools to dissect the establishment of a gene expression border separating the mesoderm and neuroectoderm in the early Drosophila embryo. We first characterized the formation of a transient ventral-to-dorsal gradient of the Snail (Sna) repressor and then examined the relationship between this gradient and repression of neural target genes in the mesoderm. We found that neural genes are repressed in a nested pattern within a zone of the mesoderm abutting the neuroectoderm, where Sna levels are graded. While several factors may contribute to the transient graded response to the Sna gradient, our analysis suggests that gene length may play an important, albeit transient, role in establishing these distinct transcriptional responses. One prediction of the gene-length-dependent transcriptional patterning model is that the co-regulated genes knirps (a short gene) and knirps-related (a long gene) should be transiently expressed in domains of differing widths, which we confirmed experimentally. These findings suggest that gene length may contribute to establishing graded responses to morphogen gradients by providing transient prepatterns that are subsequently amplified and stabilized by traditional cis-regulatory interactions.

Negative elongation factor accelerates the rate at which heat shock genes are shut off by facilitating dissociation of heat shock factor

Promoter-proximal pausing of RNA polymerase II (Pol II) occurs on thousands of genes in animal cells. This pausing often correlates with the rapid induction of genes, but direct tests of the relationship between pausing and induction rates are lacking. hsp70 and hsp26 in Drosophila are rapidly induced by heat shock. Contrary to current expectations, depletion of negative elongation factor (NELF), a key factor in setting up paused Pol II, reduced pausing but did not interfere with rapid induction. Instead, depletion of NELF delayed the time taken for these genes to shut off during recovery from heat shock. NELF depletion also delayed the dissociation of HSF from hsp70 and hsp26, and a similar delay was observed when cells were depleted of the histone acetyltransferase CBP. CBP has been reported to associate with Pol II, and acetylation of HSF by CBP has been implicated in inhibiting the DNA-binding activity of HSF. We propose that NELF-mediated pausing allows Pol II to direct CBP-mediated acetylation of HSF, thus causing HSF to dissociate from the gene. Activators are typically viewed as controlling Pol II. Our results reveal a possible reciprocal relationship in which paused Pol II influences the activator.

Multiple enhancers ensure precision of gap gene-expression patterns in the Drosophila embryo

Segmentation of the Drosophila embryo begins with the establishment of spatially restricted gap gene-expression patterns in response to broad gradients of maternal transcription factors, such as Bicoid. Numerous studies have documented the fidelity of these expression patterns, even when embryos are subjected to genetic or environmental stress, but the underlying mechanisms for this transcriptional precision are uncertain. Here we present evidence that every gap gene contains multiple enhancers with overlapping activities to produce authentic patterns of gene expression. For example, a recently identified hunchback (hb) enhancer (located 5-kb upstream of the classic enhancer) ensures repression at the anterior pole. The combination of intronic and 5' knirps (kni) enhancers produces a faithful expression pattern, even though the intronic enhancer alone directs an abnormally broad expression pattern. We present different models for "enhancer synergy," whereby two enhancers with overlapping activities produce authentic patterns of gene expression.

Morphogen gradients: from generation to interpretation

Morphogens are long-range signaling molecules that pattern developing tissues in a concentration-dependent manner. The graded activity of morphogens within tissues exposes cells to different signal levels and leads to region-specific transcriptional responses and cell fates. In its simplest incarnation, a morphogen signal forms a gradient by diffusion from a local source and clearance in surrounding tissues. Responding cells often transduce morphogen levels in a linear fashion, which results in the graded activation of transcriptional effectors. The concentration-dependent expression of morphogen target genes is achieved by their different binding affinities for transcriptional effectors as well as inputs from other transcriptional regulators. Morphogen distribution and interpretation are the result of complex interactions between the morphogen and responding tissues. The response to a morphogen is dependent not simply on morphogen concentration but also on the duration of morphogen exposure and the state of the target cells. In this review, we describe the morphogen concept and discuss the mechanisms that underlie the generation, modulation, and interpretation of morphogen gradients.

Paused RNA polymerase II as a developmental checkpoint

The textbook view of gene activation is that the rate-limiting step is the interaction of RNA polymerase II (Pol II) with the gene's promoter. However, studies in a variety of systems, including human embryonic stem cells and the early Drosophila embryo, have begun to challenge this view. There is increasing evidence that differential gene expression often depends on the regulation of transcription elongation via the release of Pol II from the proximal promoter. I review the implications of this mechanism of gene activation with respect to the orderly unfolding of complex gene networks governing animal development.

Exact time-dependent solutions for a self-regulating gene

The exact time-dependent solution for the stochastic equations governing the behavior of a binary self-regulating gene is presented. Using the generating function technique to rephrase the master equations in terms of partial differential equations, we show that the model is totally integrable and the analytical solutions are the celebrated confluent Heun functions. Self-regulation plays a major role in the control of gene expression, and it is remarkable that such a microscopic model is completely integrable in terms of well-known complex functions.

Rapid activity-induced transcription of Arc and other IEGs relies on poised RNA polymerase II

Transcription of immediate early genes (IEGs) in neurons is highly sensitive to neuronal activity, but the mechanism underlying these early transcription events is largely unknown. We found that several IEGs, such as Arc (also known as Arg3.1), are poised for near-instantaneous transcription by the stalling of RNA polymerase II (Pol II) just downstream of the transcription start site in rat neurons. Depletion through RNA interference of negative elongation factor, a mediator of Pol II stalling, reduced the Pol II occupancy of the Arc promoter and compromised the rapid induction of Arc and other IEGs. In contrast, reduction of Pol II stalling did not prevent transcription of IEGs that were expressed later and largely lacked promoter-proximal Pol II stalling. Together, our data strongly indicate that the rapid induction of neuronal IEGs requires poised Pol II and suggest a role for this mechanism in a wide variety of transcription-dependent processes, including learning and memory.

The super elongation complex (SEC) and MLL in development and disease

Transcriptional regulation at the level of elongation is vital for the control of gene expression and metazoan development. The mixed lineage leukemia (MLL) protein and its Drosophila homolog, Trithorax, which exist within COMPASS (complex of proteins associated with Set1)-like complexes, are master regulators of development. They are required for proper homeotic gene expression, in part through methylation of histone H3 on Lys 4. In humans, the MLL gene is involved in a large number of chromosomal translocations that create chimeric proteins, fusing the N terminus of MLL to several proteins that share little sequence similarity. Several frequent translocation partners of MLL were found recently to coexist in a super elongation complex (SEC) that includes known transcription elongation factors such as eleven-nineteen lysine-rich leukemia (ELL) and P-TEFb. Importantly, the SEC is required for HOX gene expression in leukemic cells, suggesting that chromosomal translocations involving MLL could lead to the overexpression of HOX and other genes through the involvement of the SEC. Here, we review the normal developmental roles of MLL and the SEC, and how MLL fusion proteins can mediate leukemogenesis.

A rapid, extensive, and transient transcriptional response to estrogen signaling in breast cancer cells

We report the immediate effects of estrogen signaling on the transcriptome of breast cancer cells using global run-on and sequencing (GRO-seq). The data were analyzed using a new bioinformatic approach that allowed us to identify transcripts directly from the GRO-seq data. We found that estrogen signaling directly regulates a strikingly large fraction of the transcriptome in a rapid, robust, and unexpectedly transient manner. In addition to protein-coding genes, estrogen regulates the distribution and activity of all three RNA polymerases and virtually every class of noncoding RNA that has been described to date. We also identified a large number of previously undetected estrogen-regulated intergenic transcripts, many of which are found proximal to estrogen receptor binding sites. Collectively, our results provide the most comprehensive measurement of the primary and immediate estrogen effects to date and a resource for understanding rapid signal-dependent transcription in other systems.

Transcriptional regulation: effects of promoter proximal pausing on speed, synchrony and reliability

Recent whole genome polymerase binding assays in the Drosophila embryo have shown that a substantial proportion of uninduced genes have pre-assembled RNA polymerase-II transcription initiation complex (PIC) bound to their promoters. These constitute a subset of promoter proximally paused genes for which mRNA elongation instead of promoter access is regulated. This difference can be described as a rearrangement of the regulatory topology to control the downstream transcriptional process of elongation rather than the upstream transcriptional initiation event. It has been shown experimentally that genes with the former mode of regulation tend to induce faster and more synchronously, and that promoter-proximal pausing is observed mainly in metazoans, in accord with a posited impact on synchrony. However, it has not been shown whether or not it is the change in the regulated step per se that is causal. We investigate this question by proposing and analyzing a continuous-time Markov chain model of PIC assembly regulated at one of two steps: initial polymerase association with DNA, or release from a paused, transcribing state. Our analysis demonstrates that, over a wide range of physical parameters, increased speed and synchrony are functional consequences of elongation control. Further, we make new predictions about the effect of elongation regulation on the consistent control of total transcript number between cells. We also identify which elements in the transcription induction pathway are most sensitive to molecular noise and thus possibly the most evolutionarily constrained. Our methods produce symbolic expressions for quantities of interest with reasonable computational effort and they can be used to explore the interplay between interaction topology and molecular noise in a broader class of biochemical networks. We provide general-purpose code implementing these methods.

Gene activation is an inherently random process because numerous diffusing proteins and DNA must first interact by random association before transcription can begin. For many genes the necessary protein–DNA associations only begin after activation, but it has recently been noted that a large class of genes in multicellular organisms can assemble the initiation complex of proteins on the core promoter prior to activation. For these genes, activation merely releases polymerase from the preassembled complex to transcribe the gene. It has been proposed on the basis of experiments that such a mechanism, while possibly costly, increases both the speed and the synchrony of the process of gene transcription. We study a realistic model of gene transcription, and show that this conclusion holds for all but a tiny fraction of the space of physical rate parameters that govern the process. The improved control of cell-to-cell variations afforded by regulation through a paused polymerase may help multicellular organisms achieve the high degree of coordination required for development. Our approach has also generated tools with which one can study the effects of analogous changes in other molecular networks and determine the relative importance of various molecular binding rates to particular system properties.

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.

Drosophila heat shock system as a general model to investigate transcriptional regulation

Whereas the regulation of a gene is uniquely tailored to respond to specific biological needs, general transcriptional mechanisms are used by diversely regulated genes within and across species. The primary mode of regulation is achieved by modulating specific steps in the transcription cycle of RNA polymerase II (Pol II). Pol II "pausing" has recently been identified as a prevalent rate-limiting and regulated step in the transcription cycle. Many sequence-specific transcription factors (TFs) modulate the duration of the pause by directly or indirectly recruiting positive transcription elongation factor b (P-TEFb) kinase, which promotes escape of Pol II from the pause into productive elongation. These specialized TFs find their target-binding sites by discriminating between DNA sequence elements based on the chromatin context in which these elements reside and can result in productive changes in gene expression or nonfunctional "promiscuous" binding. The binding of a TF can precipitate drastic changes in chromatin architecture that can be both dependent and independent of active Pol II transcription. Here, we highlight heat-shock-mediated gene transcription as a model system in which to study common mechanistic features of gene regulation.

Regulating RNA polymerase pausing and transcription elongation in embryonic stem cells

Transitions between pluripotent stem cells and differentiated cells are executed by key transcription regulators. Comparative measurements of RNA polymerase distribution over the genome's primary transcription units in different cell states can identify the genes and steps in the transcription cycle that are regulated during such transitions. To identify the complete transcriptional profiles of RNA polymerases with high sensitivity and resolution, as well as the critical regulated steps upon which regulatory factors act, we used genome-wide nuclear run-on (GRO-seq) to map the density and orientation of transcriptionally engaged RNA polymerases in mouse embryonic stem cells (ESCs) and mouse embryonic fibroblasts (MEFs). In both cell types, progression of a promoter-proximal, paused RNA polymerase II (Pol II) into productive elongation is a rate-limiting step in transcription of ∼40% of mRNA-encoding genes. Importantly, quantitative comparisons between cell types reveal that transcription is controlled frequently at paused Pol II's entry into elongation. Furthermore, "bivalent" ESC genes (exhibiting both active and repressive histone modifications) bound by Polycomb group complexes PRC1 (Polycomb-repressive complex 1) and PRC2 show dramatically reduced levels of paused Pol II at promoters relative to an average gene. In contrast, bivalent promoters bound by only PRC2 allow Pol II pausing, but it is confined to extremely 5' proximal regions. Altogether, these findings identify rate-limiting targets for transcription regulation during cell differentiation.

Cellular decision making and biological noise: from microbes to mammals

Cellular decision making is the process whereby cells assume different, functionally important and heritable fates without an associated genetic or environmental difference. Such stochastic cell fate decisions generate nongenetic cellular diversity, which may be critical for metazoan development as well as optimized microbial resource utilization and survival in a fluctuating, frequently stressful environment. Here, we review several examples of cellular decision making from viruses, bacteria, yeast, lower metazoans, and mammals, highlighting the role of regulatory network structure and molecular noise. We propose that cellular decision making is one of at least three key processes underlying development at various scales of biological organization.

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

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.

Promoter proximal pausing and the control of gene expression

The advent of methods for mapping the location of specific proteins across genomes is substantially enlightening our understanding of gene regulation. One recent discovery is that Pol II is concentrated at the 5' end of thousands of genes in mammalian and Drosophila cells. Before this, much research had focused on understanding how sequence-specific, DNA-binding proteins orchestrate the actions of regulators of chromatin structure and the general transcriptional machinery to control transcription initiation. The concentration of Pol II at the 5' ends of genes indicates that key steps regulating transcription occur after Pol II has associated with a gene's promoter.

Gene regulatory networks and the role of robustness and stochasticity in the control of gene expression

In any given cell, thousands of genes are expressed and work in concert to ensure the cell's function, fitness, and survival. Each gene, in turn, must be expressed at the proper time and in the proper amounts to ensure the appropriate functional outcome. The regulation and expression of some genes are highly robust; their expression is controlled by invariable expression programs. For instance, developmental gene expression is extremely similar in a given cell type from one individual to another. The expression of other genes is more variable: Their levels are noisy and are different from cell to cell and from individual to individual. This can be highly beneficial in physiological responses to outside cues and stresses. Recent advances have enabled the analysis of differential gene expression at a systems level. Gene regulatory networks (GRNs) involving interactions between large numbers of genes and their regulators have been mapped onto graphic diagrams that are used to visualize the regulatory relationships. The further characterization of GRNs has already uncovered global principles of gene regulation. Together with synthetic network biology, such studies are starting to provide insights into the transcriptional mechanisms that cause robust versus stochastic gene expression and their relationships to phenotypic robustness and variability. Here, we discuss GRNs and their topological properties in relation to transcriptional and phenotypic outputs in development and organismal physiology.

Gene expression noise in spatial patterning: hunchback promoter structure affects noise amplitude and distribution in Drosophila segmentation

Positional information in developing embryos is specified by spatial gradients of transcriptional regulators. One of the classic systems for studying this is the activation of the hunchback (hb) gene in early fruit fly (Drosophila) segmentation by the maternally-derived gradient of the Bicoid (Bcd) protein. Gene regulation is subject to intrinsic noise which can produce variable expression. This variability must be constrained in the highly reproducible and coordinated events of development. We identify means by which noise is controlled during gene expression by characterizing the dependence of hb mRNA and protein output noise on hb promoter structure and transcriptional dynamics. We use a stochastic model of the hb promoter in which the number and strength of Bcd and Hb (self-regulatory) binding sites can be varied. Model parameters are fit to data from WT embryos, the self-regulation mutant hb(14F), and lacZ reporter constructs using different portions of the hb promoter. We have corroborated model noise predictions experimentally. The results indicate that WT (self-regulatory) Hb output noise is predominantly dependent on the transcription and translation dynamics of its own expression, rather than on Bcd fluctuations. The constructs and mutant, which lack self-regulation, indicate that the multiple Bcd binding sites in the hb promoter (and their strengths) also play a role in buffering noise. The model is robust to the variation in Bcd binding site number across a number of fly species. This study identifies particular ways in which promoter structure and regulatory dynamics reduce hb output noise. Insofar as many of these are common features of genes (e.g. multiple regulatory sites, cooperativity, self-feedback), the current results contribute to the general understanding of the reproducibility and determinacy of spatial patterning in early development.

A transgenic mouse for in vivo detection of endogenous labeled mRNA

Live-cell single mRNA imaging is a powerful tool but has been restricted in higher eukaryotes to artificial cell lines and reporter genes. We describe an approach that enables live-cell imaging of single endogenous labeled mRNA molecules transcribed in primary mammalian cells and tissue. We generated a knock-in mouse line with an MS2 binding site (MBS) cassette targeted to the 3' untranslated region of the essential β-actin gene. As β-actin-MBS was ubiquitously expressed, we could uniquely address endogenous mRNA regulation in any tissue or cell type. We simultaneously followed transcription from the β-actin alleles in real time and observed transcriptional bursting in response to serum stimulation with precise temporal resolution. We tracked single endogenous labeled mRNA particles being transported in primary hippocampal neurons. The MBS cassette also enabled high-sensitivity fluorescence in situ hybridization (FISH), allowing detection and localization of single β-actin mRNA molecules in various mouse tissues.

Pol II waiting in the starting gates: Regulating the transition from transcription initiation into productive elongation

Proper regulation of gene expression is essential for the differentiation, development and survival of all cells and organisms. Recent work demonstrates that transcription of many genes, including key developmental and stimulus-responsive genes, is regulated after the initiation step, by pausing of RNA polymerase II during elongation through the promoter-proximal region. Thus, there is great interest in better understanding the events that follow transcription initiation and the ways in which the efficiency of early elongation can be modulated to impact expression of these highly regulated genes. Here we describe our current understanding of the steps involved in the transition from an unstable initially transcribing complex into a highly stable and processive elongation complex. We also discuss the interplay between factors that affect early transcript elongation and the potential physiological consequences for genes that are regulated through transcriptional pausing.

The meaning of pausing

Pausing of RNA polymerase II (RNAPII) at the 5' end of genes is a widespread phenomenon in metazoans, but the role of this event in gene regulation is poorly understood. Gilchrist et al. (2010) now demonstrate that RNAPII pausing counteracts DNA-influenced nucleosome organization to allow precise gene activation.

Ct3d: tracking microglia motility in 3D using a novel cosegmentation approach

MOTIVATION

Cell tracking is an important method to quantitatively analyze time-lapse microscopy data. While numerous methods and tools exist for tracking cells in 2D time-lapse images, only few and very application-specific tracking tools are available for 3D time-lapse images, which is of high relevance in immunoimaging, in particular for studying the motility of microglia in vivo.

RESULTS

We introduce a novel algorithm for tracking cells in 3D time-lapse microscopy data, based on computing cosegmentations between component trees representing individual time frames using the so-called tree-assignments. For the first time, our method allows to track microglia in three dimensional confocal time-lapse microscopy images. We also evaluate our method on synthetically generated data, demonstrating that our algorithm is robust even in the presence of different types of inhomogeneous background noise.

AVAILABILITY

Our algorithm is implemented in the ct3d package, which is available under http://www.picb.ac.cn/patterns/Software/ct3d; supplementary videos are available from http://www.picb.ac.cn/patterns/Supplements/ct3d.

Recruitment timing and dynamics of transcription factors at the Hsp70 loci in living cells

Chromatin immunoprecipitation (ChIP) studies provide snapshots of factors on chromatin in cell populations. Here, we use live-cell imaging to examine at high temporal resolution the recruitment and dynamics of transcription factors to the inducible Hsp70 loci in individual Drosophila salivary gland nuclei. Recruitment of the master regulator, HSF, is first detected within 20 s of gene activation; the timing of its recruitment resolves from RNA polymerase II and P-TEFb, and these factors resolve from Spt6 and Topo I. Remarkably, the recruitment of each factor is highly synchronous between different cells. In addition, fluorescence recovery after photobleaching (FRAP) analyses show that the entry and exit of multiple factors are progressively constrained upon gene activation, suggesting the gradual formation of a transcription compartment. Furthermore, we demonstrate that poly(ADP-ribose) (PAR) polymerase activity is required to maintain the transcription compartment. We propose that PAR polymers locally retain factors in a transcription compartment.

Pausing of RNA polymerase II disrupts DNA-specified nucleosome organization to enable precise gene regulation

Metazoan transcription is controlled through either coordinated recruitment of transcription machinery to the gene promoter or regulated pausing of RNA polymerase II (Pol II) in early elongation. We report that a striking difference between genes that use these distinct regulatory strategies lies in the "default" chromatin architecture specified by their DNA sequences. Pol II pausing is prominent at highly regulated genes whose sequences inherently disfavor nucleosome formation within the gene but favor occlusion of the promoter by nucleosomes. In contrast, housekeeping genes that lack pronounced Pol II pausing show higher nucleosome occupancy downstream, but their promoters are deprived of nucleosomes regardless of polymerase binding. Our results indicate that a key role of paused Pol II is to compete with nucleosomes for occupancy of highly regulated promoters, thereby preventing the formation of repressive chromatin architecture to facilitate further or future gene activation.

Licensed to elongate: a molecular mechanism for MLL-based leukaemogenesis

The RNA polymerase II (Pol II) elongation factor (ELL) was the first translocation partner of mixed lineage leukaemia (MLL) for which a biochemical function was determined. It was therefore proposed that the regulation of the elongation stage of transcription could be fundamental to MLL-based leukaemogenesis. Recent studies have identified ELL complexed with several of the translocation partners of MLL in a transcriptional super elongation complex (SEC). These studies provide evidence for the importance of the regulation of Pol II elongation in disease pathogenesis and suggest that MLL chimaeras function by licensing Pol II transcription elongation without the appropriate checkpoints.

Shadow enhancers foster robustness of Drosophila gastrulation

Critical developmental control genes sometimes contain "shadow" enhancers that can be located in remote positions, including the introns of neighboring genes [1]. They nonetheless produce patterns of gene expression that are the same as or similar to those produced by more proximal primary enhancers. It was suggested that shadow enhancers help foster robustness in gene expression in response to environmental or genetic perturbations [2, 3]. We critically tested this hypothesis by employing a combination of bacterial artificial chromosome (BAC) recombineering and quantitative confocal imaging methods [2, 4]. Evidence is presented that the snail gene is regulated by a distal shadow enhancer located within a neighboring locus. Removal of the proximal primary enhancer does not significantly perturb snail function, including the repression of neurogenic genes and formation of the ventral furrow during gastrulation at normal temperatures. However, at elevated temperatures, there is sporadic loss of snail expression and coincident disruptions in gastrulation. Similar defects are observed at normal temperatures upon reductions in the levels of Dorsal, a key activator of snail expression (reviewed in [5]). These results suggest that shadow enhancers represent a novel mechanism of canalization whereby complex developmental processes "bring about one definite end-result regardless of minor variations in conditions" [6].

The time to measure positional information: maternal hunchback is required for the synchrony of the Bicoid transcriptional response at the onset of zygotic transcription

It is widely accepted that morphogenetic gradients determine cell identity by concentration-dependent activation of target genes. How precise is each step in the gene expression process that acts downstream of morphogens, however, remains unclear. The Bicoid morphogen is a transcription factor directly activating its target genes and provides thus a simple system to address this issue in a quantitative manner. Recent studies indicate that the Bicoid gradient is precisely established in Drosophila embryos after eight nuclear divisions (cycle 9) and that target protein expression is specified five divisions later (cycle 14), with a precision that corresponds to a relative difference of Bicoid concentration of 10%. To understand how such precision was achieved, we directly analyzed nascent transcripts of the hunchback target gene at their site of synthesis. Most anterior nuclei in cycle 11 interphasic embryos exhibit efficient biallelic transcription of hunchback and this synchronous expression is specified within a 10% difference of Bicoid concentration. The fast diffusion of Bcd-EGFP (7.7 mum(2)/s) that we captured by fluorescent correlation spectroscopy in the nucleus is consistent with this robust expression at cycle 11. However, given the interruption of transcription during mitosis, it remains too slow to be consistent with precise de novo reading of Bicoid concentration at each interphase, suggesting the existence of a memorization process that recalls this information from earlier cycles. The two anterior maternal morphogens, Bicoid and Hunchback, contribute differently to this early response: whereas Bicoid provides dose-dependent positional information along the axis, maternal Hunchback is required for the synchrony of the response and is therefore likely to be involved in this memorization process.

beta-Catenin primes organizer gene expression by recruiting a histone H3 arginine 8 methyltransferase, Prmt2

An emerging concept in development is that transcriptional poising presets patterns of gene expression in a manner that reflects a cell's developmental potential. However, it is not known how certain loci are specified in the embryo to establish poised chromatin architecture as the developmental program unfolds. We find that, in the context of transcriptional quiescence prior to the midblastula transition in Xenopus, dorsal specification by the Wnt/beta-catenin pathway is temporally uncoupled from the onset of dorsal target gene expression, and that beta-catenin establishes poised chromatin architecture at target promoters. beta-catenin recruits the arginine methyltransferase Prmt2 to target promoters, thereby establishing asymmetrically dimethylated H3 arginine 8 (R8). Recruitment of Prmt2 to beta-catenin target genes is necessary and sufficient to establish the dorsal developmental program, indicating that Prmt2-mediated histone H3(R8) methylation plays a critical role downstream of beta-catenin in establishing poised chromatin architecture and marking key organizer genes for later expression.

Phenotypic robustness conferred by apparently redundant transcriptional enhancers

Genes include cis-regulatory regions that contain transcriptional enhancers. Recent reports have shown that developmental genes often possess multiple discrete enhancer modules that drive transcription in similar spatio-temporal patterns: primary enhancers located near the basal promoter and secondary, or 'shadow', enhancers located at more remote positions. It has been proposed that the seemingly redundant activity of primary and secondary enhancers contributes to phenotypic robustness. We tested this hypothesis by generating a deficiency that removes two newly discovered enhancers of shavenbaby (svb, a transcript of the ovo locus), a gene encoding a transcription factor that directs development of Drosophila larval trichomes. At optimal temperatures for embryonic development, this deficiency causes minor defects in trichome patterning. In embryos that develop at both low and high extreme temperatures, however, absence of these secondary enhancers leads to extensive loss of trichomes. These temperature-dependent defects can be rescued by a transgene carrying a secondary enhancer driving transcription of the svb cDNA. Finally, removal of one copy of wingless, a gene required for normal trichome patterning, causes a similar loss of trichomes only in flies lacking the secondary enhancers. These results support the hypothesis that secondary enhancers contribute to phenotypic robustness in the face of environmental and genetic variability.

Moonshine illuminates a developmental role for regulated transcription elongation

Mutations in the zebrafish gene moonshine, encoding the ortholog of TIF1 gamma, cause profound anemia and embryonic lethality. In a recent issue of Cell, Bai et al. provide evidence that these defects arise from inefficient transcription elongation, implicating elongation as an important point of regulation during cell differentiation and development.

Interactions between DSIF (DRB sensitivity inducing factor), NELF (negative elongation factor), and the Drosophila RNA polymerase II transcription elongation complex

Negative elongation factor (NELF) and 5,6-dichloro-1-beta-D-ribofuranosylbenzimidazole sensitivity-inducing factor (DSIF) are involved in pausing RNA Polymerase II (Pol II) in the promoter-proximal region of the hsp70 gene in Drosophila, before heat shock induction. Such blocks in elongation are widespread in the Drosophila genome. However, the mechanism by which DSIF and NELF participate in setting up the paused Pol II remains unclear. We analyzed the interactions among DSIF, NELF, and a reconstituted Drosophila Pol II elongation complex to gain insight into the mechanism of pausing. Our results show that DSIF and NELF require a nascent transcript longer than 18 nt to stably associate with the Pol II elongation complex. Protein-RNA cross-linking reveals that Spt5, the largest subunit of DSIF, contacts the nascent RNA as the RNA emerges from the elongation complex. Taken together, these results provide a possible model by which DSIF binds the elongation complex via association with the nascent transcript and subsequently recruits NELF. Although DSIF and NELF were both required for inhibition of transcription, we did not detect a NELF-RNA contact when the nascent transcript was between 22 and 31 nt long, which encompasses the region where promoter-proximal pausing occurs on many genes in Drosophila. This raises the possibility that RNA binding by NELF is not necessary in promoter-proximal pausing.

Genome architecture and the role of transcription

During development or in response to environmental stimuli, eukaryotic genes change both their expression and position in 3D nuclear space. Then, is a gene transcribed because of its position, or is position determined by transcription? Are genes stochastically or deterministically engaged in transcription cycles? Recent results confirm that RNA polymerases and their transcription factors play central roles in genome organization, and that stochastic events can give rise to apparently deterministic expression. As is so often the case in biology, structure both determines function and is influenced by it.

Regulation of RNA polymerase II elongation by c-Myc

The genome-wide mapping presented in this issue by Rahl et al. (2010) demonstrates how the transcription factor c-Myc interfaces with the machinery that controls RNA polymerase II elongation. These findings account for the ability of c-Myc to activate the expression of a large number of genes that promote rapid proliferation.

Inducible gene expression: diverse regulatory mechanisms

The rapid activation of gene expression in response to stimuli occurs largely through the regulation of RNA polymerase II-dependent transcription. In this Review, we discuss events that occur during the transcription cycle in eukaryotes that are important for the rapid and specific activation of gene expression in response to external stimuli. In addition to regulated recruitment of the transcription machinery to the promoter, it has now been shown that control steps can include chromatin remodelling and the release of paused polymerase. Recent work suggests that some components of signal transduction cascades also play an integral part in activating transcription at target genes.

Chromatin signature of embryonic pluripotency is established during genome activation

After fertilization the embryonic genome is inactive until transcription is initiated during the maternal-zygotic transition. This transition coincides with the formation of pluripotent cells, which in mammals can be used to generate embryonic stem cells. To study the changes in chromatin structure that accompany pluripotency and genome activation, we mapped the genomic locations of histone H3 molecules bearing lysine trimethylation modifications before and after the maternal-zygotic transition in zebrafish. Histone H3 lysine 27 trimethylation (H3K27me3), which is repressive, and H3K4me3, which is activating, were not detected before the transition. After genome activation, more than 80% of genes were marked by H3K4me3, including many inactive developmental regulatory genes that were also marked by H3K27me3. Sequential chromatin immunoprecipitation demonstrated that the same promoter regions had both trimethylation marks. Such bivalent chromatin domains also exist in embryonic stem cells and are thought to poise genes for activation while keeping them repressed. Furthermore, we found many inactive genes that were uniquely marked by H3K4me3. Despite this activating modification, these monovalent genes were neither expressed nor stably bound by RNA polymerase II. Inspection of published data sets revealed similar monovalent domains in embryonic stem cells. Moreover, H3K4me3 marks could form in the absence of both sequence-specific transcriptional activators and stable association of RNA polymerase II, as indicated by the analysis of an inducible transgene. These results indicate that bivalent and monovalent domains might poise embryonic genes for activation and that the chromatin profile associated with pluripotency is established during the maternal-zygotic transition.