High mobility of bicoid captured by fluorescence correlation spectroscopy: implication for the rapid establishment of its gradient

4D Single-particle tracking with asynchronous read-out single-photon avalanche diode array detector

Single-particle tracking techniques enable investigation of the complex functions and interactions of individual particles in biological environments. Many such techniques exist, each demonstrating trade-offs between spatiotemporal resolution, spatial and temporal range, technical complexity, and information content. To mitigate these trade-offs, we enhanced a confocal laser scanning microscope with an asynchronous read-out single-photon avalanche diode array detector. This detector provides an image of the particle's emission, precisely reflecting its position within the excitation volume. This localization is utilized in a real-time feedback system to drive the microscope scanning mechanism and ensure the particle remains centered inside the excitation volume. As each pixel is an independent single-photon detector, single-particle tracking is combined with fluorescence lifetime measurement. Our system achieves 40 nm lateral and 60 nm axial localization precision with 100 photons and sub-millisecond temporal sampling for real-time tracking. Offline tracking can refine this precision to the microsecond scale. We validated the system's spatiotemporal resolution by tracking fluorescent beads with diffusion coefficients up to 10 μm2/s. Additionally, we investigated the movement of lysosomes in living SK-N-BE cells and measured the fluorescence lifetime of the marker expressed on a membrane protein. We expect that this implementation will open other correlative imaging and tracking studies.

Transfer of polarity information via diffusion of Wnt ligands in C. elegans embryos

Different signaling mechanisms concur to ensure robust tissue patterning and cell fate instruction during animal development. Most of these mechanisms rely on signaling proteins that are produced, transported, and detected. The spatiotemporal dynamics of signaling molecules are largely unknown, yet they determine signal activity's spatial range and time frame. Here, we use the Caenorhabditis elegans embryo to study how Wnt ligands, an evolutionarily conserved family of signaling proteins, dynamically organize to establish cell polarity in a developing tissue. We identify how Wnt ligands, produced in the posterior half of the embryos, spread extracellularly to transmit information to distant target cells in the anterior half. With quantitative live imaging and fluorescence correlation spectroscopy, we show that Wnt ligands diffuse through the embryo over a timescale shorter than the cell cycle, in the intercellular space, and outside the tissue below the eggshell. We extracted diffusion coefficients of Wnt ligands and their receptor Frizzled and characterized their co-localization. Integrating our different measurements and observations in a simple computational framework, we show how fast diffusion in the embryo can polarize individual cells through a time integration of the arrival of the ligands at the target cells. The polarity established at the tissue level by a posterior Wnt source can be transferred to the cellular level. Our results support a diffusion-based long-range Wnt signaling, which is consistent with the dynamics of developing processes.

Deep learning reduces data requirements and allows real-time measurements in imaging FCS

Imaging fluorescence correlation spectroscopy (FCS) is a powerful tool to extract information on molecular mobilities, actions, and interactions in live cells, tissues, and organisms. Nevertheless, several limitations restrict its applicability. First, FCS is data hungry, requiring 50,000 frames at 1-ms time resolution to obtain accurate parameter estimates. Second, the data size makes evaluation slow. Third, as FCS evaluation is model dependent, data evaluation is significantly slowed unless analytic models are available. Here, we introduce two convolutional neural networks-FCSNet and ImFCSNet-for correlation and intensity trace analysis, respectively. FCSNet robustly predicts parameters in 2D and 3D live samples. ImFCSNet reduces the amount of data required for accurate parameter retrieval by at least one order of magnitude and makes correct estimates even in moderately defocused samples. Both convolutional neural networks are trained on simulated data, are model agnostic, and allow autonomous, real-time evaluation of imaging FCS measurements.

Transcription factor clusters as information transfer agents

Deciphering how genes interpret information from the concentration of transcription factors (TFs) within the cell nucleus remains a fundamental question in gene regulation. Recent advancements have unveiled the heterogeneous distribution of TF molecules in the nucleus, posing challenges to the precise decoding of concentration signals. To explore this phenomenon, we employ high-resolution single-cell imaging of a fluorescently tagged TF protein, Bicoid, in living fly embryos. We show that accumulation of Bicoid in submicron clusters preserves the spatial information of the maternal Bicoid gradient, and that cluster intensity, size, and frequency offer remarkably precise spatial cues. We further discover that various known gene targets of Bicoid activation colocalize with clusters and that for the target gene Hunchback, this colocalization is dependent on its enhancer binding affinity. Modeling information transfer through these clusters suggests that clustering offers a more rapid sensing mechanism for global nuclear concentrations than freely diffusing TF molecules detected by simple enhancers.

Long-range formation of the Bicoid gradient requires multiple dynamic modes that spatially vary across the embryo

Morphogen gradients provide essential positional information to gene networks through their spatially heterogeneous distribution, yet how they form is still hotly contested, with multiple models proposed for different systems. Here, we focus on the transcription factor Bicoid (Bcd), a morphogen that forms an exponential gradient across the anterior-posterior (AP) axis of the early Drosophila embryo. Using fluorescence correlation spectroscopy we find there are spatial differences in Bcd diffusivity along the AP axis, with Bcd diffusing more rapidly in the posterior. We establish that such spatially varying differences in Bcd dynamics are sufficient to explain how Bcd can have a steep exponential gradient in the anterior half of the embryo and yet still have an observable fraction of Bcd near the posterior pole. In the nucleus, we demonstrate that Bcd dynamics are impacted by binding to DNA. Addition of the Bcd homeodomain to eGFP::NLS qualitatively replicates the Bcd concentration profile, suggesting this domain regulates Bcd dynamics. Our results reveal how a long-range gradient can form while retaining a steep profile through much of its range.

Revisiting bicoid function: complete inactivation reveals an additional fundamental role in Drosophila egg geometry specification

INTRODUCTION

The bicoid (bcd) gene in Drosophila has served as a paradigm for a morphogen in textbooks for decades. Discovered in 1986 as a mutation affecting anterior development in the embryo, its expression pattern as a protein gradient later confirmed the prediction from transplantation experiments. These experiments suggested that the protein fulfills the criteria of a true morphogen, with the existence of a homeodomain crucial for activation of genes along the anterior-posterior axis, based on the concentration of the morphogen. The bcd gene undergoes alternative splicing, resulting in, among other isoforms, a small and often neglected isoform with low abundance, which lacks the homeodomain, termed small bicoid (smbcd). Most importantly, all known classical strong bcd alleles used in the past to determine bcd function apparently do not affect the function of this isoform.

RESULTS

To overcome the uncertainty regarding which isoform regulates what, I removed the bcd locus entirely using CRISPR technology. bcdCRISPR eggs exhibited a short and round appearance. The phenotype could be ascribed to smbcd because all bcd alleles affecting the function of the major transcript, termed large bicoid (lgbcd) showed normally sized eggs. Several patterning genes for the embryo showed expression in the oocyte, and their expression patterns were altered in bcdCRISPR oocytes. In bcdCRISPR embryos, all downstream segmentation genes showed altered expression patterns, consistent with the expression patterns in "classical" alleles; however, due to the altered egg geometry resulting in fewer blastoderm nuclei, additional constraints came into play, further affecting their expression patterns.

CONCLUSIONS

This study unveils a novel and fundamental role of bcd in shaping the egg's geometry. This discovery demands a comprehensive revision of our understanding of this important patterning gene and prompts a reevaluation of past experiments conducted under the assumption that bcd mutants were bcdnull-mutants.

Current capabilities and future perspectives of FCS: super-resolution microscopy, machine learning, and in vivo applications

Fluorescence correlation spectroscopy (FCS) is a single molecule sensitive tool for the quantitative measurement of biomolecular dynamics and interactions. Improvements in biology, computation, and detection technology enable real-time FCS experiments with multiplexed detection even in vivo. These new imaging modalities of FCS generate data at the rate of hundreds of MB/s requiring efficient data processing tools to extract information. Here, we briefly review FCS's capabilities and limitations before discussing recent directions that address these limitations with a focus on imaging modalities of FCS, their combinations with super-resolution microscopy, new evaluation strategies, especially machine learning, and applications in vivo.

Patterning principles of morphogen gradients

Metazoan embryos develop from a single cell into three-dimensional structured organisms while groups of genetically identical cells attain specialized identities. Cells of the developing embryo both create and accurately interpret morphogen gradients to determine their positions and make specific decisions in response. Here, we first cover intellectual roots of morphogen and positional information concepts. Focusing on animal embryos, we then provide a review of current understanding on how morphogen gradients are established and how their spans are controlled. Lastly, we cover how gradients evolve in time and space during development, and how they encode information to control patterning. In sum, we provide a list of patterning principles for morphogen gradients and review recent advances in quantitative methodologies elucidating information provided by morphogens.

How does an organism extract relevant information from transcription factor concentrations?

How does an organism regulate its genes? The involved regulation typically occurs in terms of a signal processing chain: an externally applied stimulus or a maternally supplied transcription factor leads to the expression of some downstream genes, which, in turn, are transcription factors for further genes. Especially during development, these transcription factors are frequently expressed in amounts where noise is still important; yet, the signals that they provide must not be lost in the noise. Thus, the organism needs to extract exactly relevant information in the signal. New experimental approaches involving single-molecule measurements at high temporal precision as well as increased precision in manipulations directly on the genome are allowing us to tackle this question anew. These new experimental advances mean that also from the theoretical side, theoretical advances should be possible. In this review, I will describe, specifically on the example of fly embryo gene regulation, how theoretical approaches, especially from inference and information theory, can help in understanding gene regulation. To do so, I will first review some more traditional theoretical models for gene regulation, followed by a brief discussion of information-theoretical approaches and when they can be applied. I will then introduce early fly development as an exemplary system where such information-theoretical approaches have traditionally been applied and can be applied; I will specifically focus on how one such method, namely the information bottleneck approach, has recently been used to infer structural features of enhancer architecture.

Synthetic reconstruction of the hunchback promoter specifies the role of Bicoid, Zelda and Hunchback in the dynamics of its transcription

For over 40 years, the Bicoid-hunchback (Bcd-hb) system in the fruit fly embryo has been used as a model to study how positional information in morphogen concentration gradients is robustly translated into step-like responses. A body of quantitative comparisons between theory and experiment have since questioned the initial paradigm that the sharp hb transcription pattern emerges solely from diffusive biochemical interactions between the Bicoid transcription factor and the gene promoter region. Several alternative mechanisms have been proposed, such as additional sources of positional information, positive feedback from Hb proteins or out-of-equilibrium transcription activation. By using the MS2-MCP RNA-tagging system and analysing in real time, the transcription dynamics of synthetic reporters for Bicoid and/or its two partners Zelda and Hunchback, we show that all the early hb expression pattern features and temporal dynamics are compatible with an equilibrium model with a short decay length Bicoid activity gradient as a sole source of positional information. Meanwhile, Bicoid’s partners speed-up the process by different means: Zelda lowers the Bicoid concentration threshold required for transcriptional activation while Hunchback reduces burstiness and increases the polymerase firing rate.

Trading bits in the readout from a genetic network

In the regulation of gene expression, information of relevance to the organism is represented by the concentrations of transcription factor molecules. To extract this information the cell must effectively "measure" these concentrations, but there are physical limits to the precision of these measurements. We use the gap gene network in the early fly embryo as an example of the tradeoff between the precision of concentration measurements and the transmission of relevant information. For thresholded measurements we find that lower thresholds are more important, and fine tuning is not required for near-optimal information transmission. We then consider general sensors, constrained only by a limit on their information capacity, and find that thresholded sensors can approach true information theoretic optima. The information theoretic approach allows us to identify the optimal sensor for the entire gap gene network and to argue that the physical limitations of sensing necessitate the observed multiplicity of enhancer elements, with sensitivities to combinations rather than single transcription factors.

Using FCS to accurately measure protein concentration in the presence of noise and photobleaching

Quantitative cell biology requires precise and accurate concentration measurements, resolved both in space and time. Fluorescence correlation spectroscopy (FCS) has been held as a promising technique to perform such measurements because the fluorescence fluctuations it relies on are directly dependent on the absolute number of fluorophores in the detection volume. However, the most interesting applications are in cells, where autofluorescence and confinement result in strong background noise and important levels of photobleaching. Both noise and photobleaching introduce systematic bias in FCS concentration measurements and need to be corrected for. Here, we propose to make use of the photobleaching inevitably occurring in confined environments to perform series of FCS measurements at different fluorophore concentration, which we show allows a precise in situ measurement of both background noise and molecular brightness. Such a measurement can then be used as a calibration to transform confocal intensity images into concentration maps. The power of this approach is first illustrated with in vitro measurements using different dye solutions, then its applicability for in vivo measurements is demonstrated in Drosophila embryos for a model nuclear protein and for two morphogens, Bicoid and Capicua.

High rates of calcium-free diffusion in the cytosol of living cells

Calcium (Ca²⁺) is a universal second messenger that participates in the regulation of innumerous physiological processes. The way in which local elevations of the cytosolic Ca²⁺ concentration spread in space and time is key for the versatility of the signals. Ca²⁺ diffusion in the cytosol is hindered by its interaction with proteins that act as buffers. Depending on the concentrations and the kinetics of the interactions, there is a large range of values at which Ca²⁺ diffusion can proceed. Having reliable estimates of this range, particularly of its highest end, which corresponds to the ions free diffusion, is key to understand how the signals propagate. In this work, we present the first experimental results with which the Ca²⁺-free diffusion coefficient is directly quantified in the cytosol of living cells. By means of fluorescence correlation spectroscopy experiments performed in Xenopus laevis oocytes and in cells of Saccharomyces cerevisiae, we show that the ions can freely diffuse in the cytosol at a higher rate than previously thought.

Capicua is a fast-acting transcriptional brake

Even though transcriptional repressors are studied with ever-increasing molecular resolution, the temporal aspects of gene repression remain poorly understood. Here, we address the dynamics of transcriptional repression by Capicua (Cic), which is essential for normal development and is commonly mutated in human cancers and neurodegenerative diseases.^1,2 We report the speed limit for Cic-dependent gene repression based on live imaging and optogenetic perturbations in the early Drosophila embryo, where Cic was originally discovered.³ Our measurements of Cic concentration and intranuclear mobility, along with real-time monitoring of the activity of Cic target genes, reveal remarkably fast transcriptional repression within minutes of removing an optogenetic de-repressive signal. In parallel, quantitative analyses of transcriptional bursting of Cic target genes support a repression mechanism providing a fast-acting brake on burst generation. This work sets quantitative constraints on potential mechanisms for gene regulation by Cic.

Modulating the bicoid gradient in space and time

Background

The formation of the Bicoid (Bcd) gradient in the early Drosophila is one of the most fascinating observations in biology and serves as a paradigm for gradient formation, yet its mechanism is still not fully understood. Two distinct models were proposed in the past, the SDD and the ARTS model.

Results

We define novel cis- and trans-acting factors that are indispensable for gradient formation. The first one is the poly A tail length of the bcd mRNA where we demonstrate that it changes not only in time, but also in space. We show that posterior bcd mRNAs possess a longer poly tail than anterior ones and this elongation is likely mediated by wispy (wisp), a poly A polymerase. Consequently, modulating the activity of Wisp results in changes of the Bcd gradient, in controlling downstream targets such as the gap and pair-rule genes, and also in influencing the cuticular pattern. Attempts to modulate the Bcd gradient by subjecting the egg to an extra nuclear cycle, i.e. a 15^th nuclear cycle by means of the maternal haploid (mh) mutation showed no effect, neither on the appearance of the gradient nor on the control of downstream target. This suggests that the segmental anlagen are determined during the first 14 nuclear cycles. Finally, we identify the Cyclin B (CycB) gene as a trans-acting factor that modulates the movement of Bcd such that Bcd movement is allowed to move through the interior of the egg.

Conclusions

Our analysis demonstrates that Bcd gradient formation is far more complex than previously thought requiring a revision of the models of how the gradient is formed.

The early Drosophila embryo as a model system for quantitative biology

With the rise of new tools, from controlled genetic manipulations and optogenetics to improved microscopy, it is now possible to make clear, quantitative and reproducible measurements of biological processes. The humble fruit fly Drosophila melanogaster, with its ease of genetic manipulation combined with excellent imaging accessibility, has become a major model system for performing quantitative in vivo measurements. Such measurements are driving a new wave of interest from physicists and engineers, who are developing a range of testable dynamic models of active systems to understand fundamental biological processes. The reproducibility of the early Drosophila embryo has been crucial for understanding how biological systems are robust to unavoidable noise during development. Insights from quantitative in vivo experiments in the Drosophila embryo are having an impact on our understanding of critical biological processes, such as how cells make decisions and how complex tissue shape emerges. Here, to highlight the power of using Drosophila embryogenesis for quantitative biology, I focus on three main areas: (1) formation and robustness of morphogen gradients; (2) how gene regulatory networks ensure precise boundary formation; and (3) how mechanical interactions drive packing and tissue folding. I further discuss how such data has driven advances in modelling.

Of numbers and movement - understanding transcription factor pathogenesis by advanced microscopy

Transcription factors (TFs) are life-sustaining and, therefore, the subject of intensive research. By regulating gene expression, TFs control a plethora of developmental and physiological processes, and their abnormal function commonly leads to various developmental defects and diseases in humans. Normal TF function often depends on gene dosage, which can be altered by copy-number variation or loss-of-function mutations. This explains why TF haploinsufficiency (HI) can lead to disease. Since aberrant TF numbers frequently result in pathogenic abnormalities of gene expression, quantitative analyses of TFs are a priority in the field. In vitro single-molecule methodologies have significantly aided the identification of links between TF gene dosage and transcriptional outcomes. Additionally, advances in quantitative microscopy have contributed mechanistic insights into normal and aberrant TF function. However, to understand TF biology, TF-chromatin interactions must be characterised in vivo, in a tissue-specific manner and in the context of both normal and altered TF numbers. Here, we summarise the advanced microscopy methodologies most frequently used to link TF abundance to function and dissect the molecular mechanisms underlying TF HIs. Increased application of advanced single-molecule and super-resolution microscopy modalities will improve our understanding of how TF HIs drive disease.

Quantification of fluctuations from fluorescence correlation spectroscopy experiments in reaction-diffusion systems

Fluorescence correlation spectroscopy (FCS) is commonly used to estimate diffusion and reaction rates. In FCS the fluorescence coming from a small volume is recorded and the autocorrelation function (ACF) of the fluorescence fluctuations is computed. Scaling out the fluctuations due to the emission process, this ACF can be related to the ACF of the fluctuations in the number of observed fluorescent molecules. In this paper the ACF of the molecule number fluctuations is studied theoretically, with no approximations, for a reaction-diffusion system in which the fluorescence changes with binding and unbinding. Theoretical ACFs are usually derived assuming that fluctuations in the number of molecules of one species are instantaneously uncorrelated to those of the others and obey Poisson statistics. Under these assumptions, the ACF derived in this paper is characterized only by the diffusive timescale of the fluorescent species and its total weight is the inverse of the mean number of observed fluorescent molecules. The theory is then scrutinized in view of previous experimental results which, for a similar system, gave a different total weight and correct estimates of other diffusive timescales. The total weight mismatch is corrected by assuming that the variance of the number of fluorescent molecules depends on the variance of the particle numbers of the other species, as in the variance decomposition formula. Including the finite acquisition time in its computation, it is shown that the ACF depends on various timescales of the system and that its total weight coincides with the one obtained with the variance decomposition formula. This calculation implies that diffusion coefficients of nonobservable species can be estimated with FCS experiments performed in reaction-diffusion systems. Ways to proceed in future experiments are also discussed.

Phase separation at the nanoscale quantified by dcFCCS

Liquid–liquid phase separation causes formation of membraneless condensates containing concentrated biomolecules, which play essential roles in diverse cellular processes. Currently, micrometer-scale phase-separated condensates are examined routinely to elucidate their functions and mechanisms in details. However, limited by current commonly used methods, the transition process from miscible individual molecules to micrometer-scale condensates is mostly unknown. Herein, we captured formation of nanoscale condensates, whose size, growth rate, molecular stoichiometry, and binding affinity to recruit client molecules were all quantified. To achieve this, we developed a dual-color fluorescence cross-correlation spectroscopy (dcFCCS) method, and we expect that dcFCCS can be widely applied to investigate phase separation of other biomolecules at the nanoscale.

Liquid–liquid phase separation, driven by multivalent macromolecular interactions, causes formation of membraneless compartments, which are biomolecular condensates containing concentrated macromolecules. These condensates are essential in diverse cellular processes. Formation and dynamics of micrometer-scale phase-separated condensates are examined routinely. However, limited by commonly used methods which cannot capture small-sized free-diffusing condensates, the transition process from miscible individual molecules to micrometer-scale condensates is mostly unknown. Herein, with a dual-color fluorescence cross-correlation spectroscopy (dcFCCS) method, we captured formation of nanoscale condensates beyond the detection limit of conventional fluorescence microscopy. In addition, dcFCCS is able to quantify size and growth rate of condensates as well as molecular stoichiometry and binding affinity of client molecules within condensates. The critical concentration to form nanoscale condensates, identified by our experimental measurements and Monte Carlo simulations, is at least several fold lower than the detection limit of conventional fluorescence microscopy. Our results emphasize that, in addition to micrometer-scale condensates, nanoscale condensates are likely to play important roles in various cellular processes and dcFCCS is a simple and powerful quantitative tool to examine them in detail.

Quantifying Temperature Compensation of Bicoid Gradients with a Fast T-Tunable Microfluidic Device

As a reaction-diffusion system strongly affected by temperature, early fly embryos surprisingly show highly reproducible and accurate developmental patterns during embryogenesis under temperature perturbations. To reveal the underlying temperature compensation mechanism, it is important to overcome the challenge in quantitative imaging on fly embryos under temperature perturbations. Inspired by microfluidics generating temperature steps on fly embryos, here we design a microfluidic device capable of ensuring the normal development of multiple fly embryos as well as achieving real-time temperature control and fast temperature switches for quantitative live imaging with a home-built two-photon microscope. We apply this system to quantify the temperature compensation of the morphogen Bicoid (Bcd) gradient in fly embryos. The length constant of the exponential Bcd gradient reaches the maximum at 25°C within the measured temperatures of 18-29°C and gradually adapts to the corresponding value at new temperatures upon a fast temperature switch. The relaxation time of such an adaptation becomes longer if the temperature is switched in a later developmental stage. This age-dependent temperature compensation could be explained if the traditional synthesis-diffusion-degradation model is extended to incorporate the dynamic change of the parameters controlling the formation of Bcd gradients.

A mechanism for hunchback promoters to readout morphogenetic positional information in less than a minute

Cell fate decisions in the fly embryo are rapid: hunchback genes decide in minutes whether nuclei follow the anterior/posterior developmental blueprint by reading out positional information in the Bicoid morphogen. This developmental system is a prototype of regulatory decision processes that combine speed and accuracy. Traditional arguments based on fixed-time sampling of Bicoid concentration indicate that an accurate readout is impossible within the experimental times. This raises the general issue of how speed-accuracy tradeoffs are achieved. Here, we compare fixed-time to on-the-fly decisions, based on comparing the likelihoods of anterior/posterior locations. We found that these more efficient schemes complete reliable cell fate decisions within the short embryological timescales. We discuss the influence of promoter architectures on decision times and error rates, present concrete examples that rapidly readout the morphogen, and predictions for new experiments. Lastly, we suggest a simple mechanism for RNA production and degradation that approximates the log-likelihood function.

Modeling of Wnt-mediated tissue patterning in vertebrate embryogenesis

During embryogenesis, morphogens form a concentration gradient in responsive tissue, which is then translated into a spatial cellular pattern. The mechanisms by which morphogens spread through a tissue to establish such a morphogenetic field remain elusive. Here, we investigate by mutually complementary simulations and in vivo experiments how Wnt morphogen transport by cytonemes differs from typically assumed diffusion-based transport for patterning of highly dynamic tissue such as the neural plate in zebrafish. Stochasticity strongly influences fate acquisition at the single cell level and results in fluctuating boundaries between pattern regions. Stable patterning can be achieved by sorting through concentration dependent cell migration and apoptosis, independent of the morphogen transport mechanism. We show that Wnt transport by cytonemes achieves distinct Wnt thresholds for the brain primordia earlier compared with diffusion-based transport. We conclude that a cytoneme-mediated morphogen transport together with directed cell sorting is a potentially favored mechanism to establish morphogen gradients in rapidly expanding developmental systems.

Fluorescence fluctuation analysis reveals PpV dependent Cdc25 protein dynamics in living embryos

The protein phosphatase Cdc25 is a key regulator of the cell cycle by activating Cdk-cyclin complexes. Cdc25 is regulated by its expression levels and post-translational mechanisms. In early Drosophila embryogenesis, Cdc25/Twine drives the fast and synchronous nuclear cycles. A pause in the cell cycle and the remodeling to a more generic cell cycle mode with a gap phase are determined by Twine inactivation and destruction in early interphase 14, in response to zygotic genome activation. Although the pseudokinase Tribbles contributes to the timely degradation of Twine, Twine levels are controlled by additional yet unknown post-translational mechanisms. Here, we apply a non-invasive method based on fluorescence fluctuation analysis (FFA) to record the absolute concentration profiles of Twine with minute-scale resolution in single living embryos. Employing this assay, we found that Protein phosphatase V (PpV), the homologue of the catalytic subunit of human PP6, ensures appropriately low Twine protein levels at the onset of interphase 14. PpV controls directly or indirectly the phosphorylation of Twine at multiple serine and threonine residues as revealed by phosphosite mapping. Mutational analysis confirmed that these sites are involved in control of Twine protein dynamics, and cell cycle remodeling is delayed in a fraction of the phosphosite mutant embryos. Our data reveal a novel mechanism for control of Twine protein levels and their significance for embryonic cell cycle remodeling.

Constraints and limitations on the transcriptional response downstream of the Bicoid morphogen gradient

The regulation of the hunchback promoter expression by the maternal Bicoid gradient has been studied as a model system in development for many years. Yet, at the level of quantitative agreement between data and theoretical models, even the first step of this regulation, transcription, continues to be challenging. This situation is slowly progressing, thanks to quantitative live-imaging techniques coupled to advanced statistical data analysis and modeling. Here, we outline the current state of our knowledge of this apparently "simple" step, highlighting the newly appreciated role of bursty transcription dynamics and its regulation.

A matter of time: Formation and interpretation of the Bicoid morphogen gradient

Spatially distributed signaling molecules, known as morphogens, provide spatial information during development. A host of different morphogens have now been identified, from subcellular gradients through to morphogens that act across a whole embryo. These gradients form over a wide-range of timescales, from seconds to hours, and their time windows for interpretation are also highly variable; the processes of morphogen gradient formation and interpretation are highly dynamic. The morphogen Bicoid (Bcd), present in the early Drosophila embryo, is essential for setting up the future Drosophila body segments. Due to its accessibility for both genetic perturbations and imaging, this system has provided key insights into how precise patterning can occur within a highly dynamic system. Here, we review the temporal scales of Bcd gradient formation and interpretation. In particular, we discuss the quantitative evidence for different models of Bcd gradient formation, outline the time windows for Bcd interpretation, and describe how Bcd temporally adapts its own ability to be interpreted. The utilization of temporal information in morphogen readout may provide crucial inputs to ensure precise spatial patterning, particularly in rapidly developing systems.

bicoid RNA localization requires the trans-Golgi network

Background

The formation of the bicoid (bcd) mRNA gradient is a crucial step for Bcd protein gradient formation in Drosophila. In the past, a microtubule (MT)-based cortical network had been shown to be indispensable for bcd mRNA transport to the posterior.

Results

We report the identification of a MT-binding protein CLASP/Chb as the first component associated with this cortical MT network. Since CLASPs in vertebrates were shown to serve as an acentriolar microtubule organization center (aMTOC) in concert with trans-Golgi proteins, we examined the effect of the Drosophila trans-Golgins on bcd localization and gradient formation. Using a genetic approach, we demonstrate that the Drosophila trans-Golgins dGCC88, dGolgin97 and dGCC185 indeed affect bcd mRNA localization during oocyte development. Consequently, the bcd mRNA is already mislocalized before the egg is fertilized. The expression domains of genes downstream of the hierarchy of bcd, e.g. of the gap gene empty spiracles or of the pair-rule gene even-skipped are changed, indicating an altered segmental anlagen, due to a faulty bcd gradient. Thus, at the end of embryogenesis, trans-Golgin mutants show bcd-like cuticle phenotypes.

Conclusions

Our data provides evidence that the Golgi as a cellular member of the secretory pathway exerts control on bcd localization which indicates that bcd gradient formation is probably more intricate than previously presumed.

Imaging, Visualization, and Computation in Developmental Biology

Embryonic development is highly complex and dynamic, requiring the coordination of numerous molecular and cellular events at precise times and places. Advances in imaging technology have made it possible to follow developmental processes at cellular, tissue, and organ levels over time as they take place in the intact embryo. Parallel innovations of in vivo probes permit imaging to report on molecular, physiological, and anatomical events of embryogenesis, but the resulting multidimensional data sets pose significant challenges for extracting knowledge. In this review, we discuss recent and emerging advances in imaging technologies, in vivo labeling, and data processing that offer the greatest potential for jointly deciphering the intricate cellular dynamics and the underlying molecular mechanisms. Our discussion of the emerging area of “image-omics” highlights both the challenges of data analysis and the promise of more fully embracing computation and data science for rapidly advancing our understanding of biology.

Fluorescence Correlation and Cross-Correlation Spectroscopy in Zebrafish

There has been increasing interest in biophysical studies on live organisms to gain better insights into physiologically relevant biological events at the molecular level. Zebrafish (Danio rerio) is a viable vertebrate model to study such events due to its genetic and evolutionary similarities to humans, amenability to less invasive fluorescence techniques owing to its transparency and well-characterized genetic manipulation techniques. Fluorescence techniques used to probe biomolecular dynamics and interactions of molecules in live zebrafish embryos are therefore highly sought-after to bridge molecular and developmental events. Fluorescence correlation and cross-correlation spectroscopy (FCS and FCCS) are two robust techniques that provide molecular level information on dynamics and interactions respectively. Here, we detail the steps for applying confocal FCS and FCCS, in particular single-wavelength FCCS (SW-FCCS), in live zebrafish embryos, beginning with sample preparation, instrumentation, calibration, and measurements on the FCS/FCCS instrument and ending with data analysis.

Temporal control of gene expression by the pioneer factor Zelda through transient interactions in hubs

Pioneer transcription factors can engage nucleosomal DNA, which leads to local chromatin remodeling and to the establishment of transcriptional competence. However, the impact of enhancer priming by pioneer factors on the temporal control of gene expression and on mitotic memory remains unclear. Here we employ quantitative live imaging methods and mathematical modeling to test the effect of the pioneer factor Zelda on transcriptional dynamics and memory in Drosophila embryos. We demonstrate that increasing the number of Zelda binding sites accelerates the kinetics of nuclei transcriptional activation regardless of their transcriptional past. Despite its known pioneering activities, we show that Zelda does not remain detectably associated with mitotic chromosomes and is neither necessary nor sufficient to foster memory. We further reveal that Zelda forms sub-nuclear dynamic hubs where Zelda binding events are transient. We propose that Zelda facilitates transcriptional activation by accumulating in microenvironments where it could accelerate the duration of multiple pre-initiation steps.

Precision in a rush: Trade-offs between reproducibility and steepness of the hunchback expression pattern

Fly development amazes us by the precision and reproducibility of gene expression, especially since the initial expression patterns are established during very short nuclear cycles. Recent live imaging of hunchback promoter dynamics shows a stable steep binary expression pattern established within the three minute interphase of nuclear cycle 11. Considering expression models of different complexity, we explore the trade-off between the ability of a regulatory system to produce a steep boundary and minimize expression variability between different nuclei. We show how a limited readout time imposed by short developmental cycles affects the gene’s ability to read positional information along the embryo’s anterior posterior axis and express reliably. Comparing our theoretical results to real-time monitoring of the hunchback transcription dynamics in live flies, we discuss possible regulatory strategies, suggesting an important role for additional binding sites, gradients or non-equilibrium binding and modified transcription factor search strategies.

Despite very limited time, organisms develop in reproducible ways. In the early stages of fly development the information about maternal signals is read out in a few minutes to produce steep and precise gene expression patterns. Motivated by recent live imaging experiments in fly embryos, we explore the consequences of the trade-off between a rushed but reproducible readout and a steep expression pattern on the regulatory modules of gene expression. We show that the current view of one anterior gradient morphogen binding to six binding sites is quantitatively inconsistent with the experimental data given the short readout time, suggesting other regulatory features.

3 minutes to precisely measure morphogen concentration

Morphogen gradients provide concentration-dependent positional information along polarity axes. Although the dynamics of the establishment of these gradients is well described, precision and noise in the downstream activation processes remain elusive. A simple paradigm to address these questions is the Bicoid morphogen gradient that elicits a rapid step-like transcriptional response in young fruit fly embryos. Focusing on the expression of the major Bicoid target, hunchback (hb), at the onset of zygotic transcription, we used the MS2-MCP approach which combines fluorescent labeling of nascent mRNA with live imaging at high spatial and temporal resolution. Removing 36 putative Zelda binding sites unexpectedly present in the original MS2 reporter, we show that the 750 bp of the hb promoter are sufficient to recapitulate endogenous expression at the onset of zygotic transcription. After each mitosis, in the anterior, expression is turned on to rapidly reach a plateau with all nuclei expressing the reporter. Consistent with a Bicoid dose-dependent activation process, the time period required to reach the plateau increases with the distance to the anterior pole. Despite the challenge imposed by frequent mitoses and high nuclei-to-nuclei variability in transcription kinetics, it only takes 3 minutes at each interphase for the MS2 reporter loci to distinguish subtle differences in Bicoid concentration and establish a steadily positioned and steep (Hill coefficient ~ 7) expression boundary. Modeling based on the cooperativity between the 6 known Bicoid binding sites in the hb promoter region, assuming rate limiting concentrations of the Bicoid transcription factor at the boundary, is able to capture the observed dynamics of pattern establishment but not the steepness of the boundary. This suggests that a simple model based only on the cooperative binding of Bicoid is not sufficient to describe the spatiotemporal dynamics of early hb expression.

Seeing is believing: the Bicoid protein reveals its path

In this commentary, I will review the latest findings on the Bicoid (Bcd) morphogen in Drosophila, a paradigm for gradient formation taught to biology students for more than two decades. “Seeing is believing” also summarizes the erroneous steps that were needed to elucidate the mechanisms of gradient formation and the path of movement of Bcd. Initially proclaimed as a dogma in 1988 and later incorporated into the SDD model where the broad diffusion of Bcd throughout the embryo was the predominant step leading to gradient formation, the SDD model was irrefutable for more than two decades until first doubts were raised in 2007 regarding the diffusion properties of Bcd associated with the SDD model. This led to re-thinking of the issue and the definition of a new model, termed the ARTS model which could explain most of the physical constraints that were inherently associated with the SDD model. In the ARTS model, gradient formation is mediated by the mRNA which is redistributed along cortical microtubules to form a mRNA gradient which is translated to form the protein gradient. Contrary to the SDD model, there is no Bcd diffusion from the tip. The ARTS model is also compatible with the observed cortical movement of Bcd. I will critically compare the SDD and the ARTS models as well as other models, analyze the major differences, and highlight the path where Bcd is localized during early nuclear cycles.

Bicoid gradient formation mechanism and dynamics revealed by protein lifetime analysis

Embryogenesis relies on instructions provided by spatially organized signaling molecules known as morphogens. Understanding the principles behind morphogen distribution and how cells interpret locally this information remains a major challenge in developmental biology. Here, we introduce morphogen-age measurements as a novel approach to test models of morphogen gradient formation. Using a tandem fluorescent timer as a protein age sensor, we find a gradient of increasing age of Bicoid along the anterior-posterior axis in the early Drosophila embryo. Quantitative analysis of the protein age distribution across the embryo reveals that the synthesis-diffusion-degradation model is the most likely model underlying Bicoid gradient formation, and rules out other hypotheses for gradient formation. Moreover, we show that the timer can detect transitions in the dynamics associated with syncytial cellularization. Our results provide new insight into Bicoid gradient formation and demonstrate how morphogen-age information can complement knowledge about movement, abundance, and distribution, which should be widely applicable to other systems.

LlamaTags: A Versatile Tool to Image Transcription Factor Dynamics in Live Embryos

Embryonic cell fates are defined by transcription factors that are rapidly deployed, yet attempts to visualize these factors in vivo often fail because of slow fluorescent protein maturation. Here, we pioneer a protein tag, LlamaTag, which circumvents this maturation limit by binding mature fluorescent proteins, making it possible to visualize transcription factor concentration dynamics in live embryos. Implementing this approach in the fruit fly Drosophila melanogaster, we discovered stochastic bursts in the concentration of transcription factors that are correlated with bursts in transcription. We further used LlamaTags to show that the concentration of protein in a given nucleus heavily depends on transcription of that gene in neighboring nuclei; we speculate that this inter-nuclear signaling is an important mechanism for coordinating gene expression to delineate straight and sharp boundaries of gene expression. Thus, LlamaTags now make it possible to visualize the flow of information along the central dogma in live embryos.

Spatially Different Tissue-Scale Diffusivity Shapes ANGUSTIFOLIA3 Gradient in Growing Leaves

The spatial gradient of signaling molecules is pivotal for establishing developmental patterns of multicellular organisms. It has long been proposed that these gradients could arise from the pure diffusion process of signaling molecules between cells, but whether this simplest mechanism establishes the formation of the tissue-scale gradient remains unclear. Plasmodesmata are unique channel structures in plants that connect neighboring cells for molecular transport. In this study, we measured cellular- and tissue-scale kinetics of molecular transport through plasmodesmata in Arabidopsis thaliana developing leaf primordia by fluorescence recovery assays. These trans-scale measurements revealed biophysical properties of diffusive molecular transport through plasmodesmata and revealed that the tissue-scale diffusivity, but not the cellular-scale diffusivity, is spatially different along the leaf proximal-to-distal axis. We found that the gradient in cell size along the developmental axis underlies this spatially different tissue-scale diffusivity. We then asked how this diffusion-based framework functions in establishing a signaling gradient of endogenous molecules. ANGUSTIFOLIA3 (AN3) is a transcriptional co-activator, and as we have shown here, it forms a long-range signaling gradient along the leaf proximal-to-distal axis to determine a cell-proliferation domain. By genetically engineering AN3 mobility, we assessed each contribution of cell-to-cell movement and tissue growth to the distribution of the AN3 gradient. We constructed a diffusion-based theoretical model using these quantitative data to analyze the AN3 gradient formation and demonstrated that it could be achieved solely by the diffusive molecular transport in a growing tissue. Our results indicate that the spatially different tissue-scale diffusivity is a core mechanism for AN3 gradient formation. This provides evidence that the pure diffusion process establishes the formation of the long-range signaling gradient in leaf development.

On the importance of protein diffusion in biological systems: The example of the Bicoid morphogen gradient

Morphogens are proteins that form concentration gradients in embryos and developing tissues, where they act as postal codes, providing cells with positional information and allowing them to behave accordingly. Bicoid was the first discovered morphogen, and remains one of the most studied. It regulates segmentation in flies, forming a striking exponential gradient along the anterior-posterior axis of early Drosophila embryos, and activating the transcription of multiple target genes in a concentration-dependent manner. In this review, the work done by us and by others to characterize the mobility of Bicoid in D. melanogaster embryos is presented. The central role played by the diffusion of Bicoid in both the establishment of the gradient and the activation of target genes is discussed, and placed in the context of the need for these processes to be all at once rapid, precise and robust. The Bicoid system, and morphogen gradients in general, remain amongst the most amazing examples of the coexistence, often observed in living systems, of small-scale disorder and large-scale spatial order. This article is part of a Special Issue entitled: Biophysics in Canada, edited by Lewis Kay, John Baenziger, Albert Berghuis and Peter Tieleman.

Fluorescence correlation spectroscopy experiments to quantify free diffusion coefficients in reaction-diffusion systems: The case of Ca^{2+} and its dyes

Many cell signaling pathways involve the diffusion of messengers that bind and unbind to and from intracellular components. Quantifying their net transport rate under different conditions then requires having separate estimates of their free diffusion coefficient and binding or unbinding rates. In this paper, we show how performing sets of fluorescence correlation spectroscopy (FCS) experiments under different conditions, it is possible to quantify free diffusion coefficients and on and off rates of reaction-diffusion systems. We develop the theory and present a practical implementation for the case of the universal second messenger, calcium (Ca^{2+}) and single-wavelength dyes that increase their fluorescence upon Ca^{2+} binding. We validate the approach with experiments performed in aqueous solutions containing Ca^{2+} and Fluo4 dextran (both in its high and low affinity versions). Performing FCS experiments with tetramethylrhodamine-dextran in Xenopus laevis oocytes, we infer the corresponding free diffusion coefficients in the cytosol of these cells. Our approach can be extended to other physiologically relevant reaction-diffusion systems to quantify biophysical parameters that determine the dynamics of various variables of interest.

Maternal AP determinants in the Drosophila oocyte and embryo

An animal embryo cannot initiate its journey of forming a new life on its own. It must rely on maternally provided resources and inputs to kick-start its developmental process. In Drosophila, the initial polarities of the embryo along both the anterior-posterior (AP) and dorsal-ventral (DV) axes are also specified by maternal determinants. Over the past several decades, genetic and molecular studies have identified and characterized such determinants, as well as the zygotic genetic regulatory networks that control patterning in the early embryo. Extensive studies of oogenesis have also led to a detailed knowledge of the cellular and molecular interactions that control the formation of a mature egg. Despite these efforts, oogenesis and embryogenesis have been studied largely as separate problems, except for qualitative aspects with regard to maternal regulation of the asymmetric localization of maternal determinants. Can oogenesis and embryogenesis be viewed from a unified perspective at a quantitative level, and can that improve our understanding of how robust embryonic patterning is achieved? Here, we discuss the basic knowledge of the regulatory mechanisms controlling oogenesis and embryonic patterning along the AP axis. We explore properties of the maternal Bicoid gradient in relation to embryo size in search for a unified framework for robust AP patterning. WIREs Dev Biol 2016, 5:562-581. doi: 10.1002/wdev.235 For further resources related to this article, please visit the WIREs website.

Characterizing anomalous diffusion in crowded polymer solutions and gels over five decades in time with variable-lengthscale fluorescence correlation spectroscopy

The diffusion of macromolecules in cells and in complex fluids is often found to deviate from simple Fickian diffusion. One explanation offered for this behavior is that molecular crowding renders diffusion anomalous, where the mean-squared displacement of the particles scales as 〈r(2)〉∝t(α) with α < 1. Unfortunately, methods such as fluorescence correlation spectroscopy (FCS) or fluorescence recovery after photobleaching (FRAP) probe diffusion only over a narrow range of lengthscales and cannot directly test the dependence of the mean-squared displacement (MSD) on time. Here we show that variable-lengthscale FCS (VLS-FCS), where the volume of observation is varied over several orders of magnitude, combined with a numerical inversion procedure of the correlation data, allows retrieving the MSD for up to five decades in time, bridging the gap between diffusion experiments performed at different lengthscales. In addition, we show that VLS-FCS provides a way to assess whether the propagator associated with the diffusion is Gaussian or non-Gaussian. We used VLS-FCS to investigate two systems where anomalous diffusion had been previously reported. In the case of dense cross-linked agarose gels, the measured MSD confirmed that the diffusion of small beads was anomalous at short lengthscales, with a cross-over to simple diffusion around ≈1 μm, consistent with a caged diffusion process. On the other hand, for solutions crowded with marginally entangled dextran molecules, we uncovered an apparent discrepancy between the MSD, found to be linear, and the propagators at short lengthscales, found to be non-Gaussian. These contradicting features call to mind the "anomalous, yet Brownian" diffusion observed in several biological systems, and the recently proposed "diffusing diffusivity" model.

Fluctuations, Correlations and the Estimation of Concentrations inside Cells

Information transmission in cells occurs quite accurately even when concentration changes are “read” by individual binding sites. In this paper we study ligand number and site occupancy fluctuations when ligands diffuse and react going beyond the analyses that focus on their asymptotic decay. In this way we show that, for immobile binding sites, fluctuations in the number of bound molecules decay on a relatively fast scale before the asymptotic behavior kicks in. This result can explain the observed co-existence of highly fluctuating instantaneous transcriptional activities with accumulated mRNA concentrations that have relatively small noise levels. We also show that the initial stages of the decay in the bound molecule number fluctuations have one or two characteristic timescales depending on the concentration of free molecules. This transition can explain the changes in enzyme activity observed at the single molecule level.

Kinetics of Reactive Modules Adds Discriminative Dimensions for Selective Cell Imaging

Living cells are chemical mixtures of exceptional interest and significance, whose investigation requires the development of powerful analytical tools fulfilling the demanding constraints resulting from their singular features. In particular, multiplexed observation of a large number of molecular targets with high spatiotemporal resolution appears highly desirable. One attractive road to address this analytical challenge relies on engaging the targets in reactions and exploiting the rich kinetic signature of the resulting reactive module, which originates from its topology and its rate constants. This review explores the various facets of this promising strategy. We first emphasize the singularity of the content of a living cell as a chemical mixture and suggest that its multiplexed observation is significant and timely. Then, we show that exploiting the kinetics of analytical processes is relevant to selectively detect a given analyte: upon perturbing the system, the kinetic window associated to response read-out has to be matched with that of the targeted reactive module. Eventually, we introduce the state-of-the-art of cell imaging exploiting protocols based on reaction kinetics and draw some promising perspectives.

Positional Information and Cell Fate Determination in the Early Drosophila Embryo

During early development in Drosophila, the spatial information of maternal gradients is translated into discrete transcriptional states determining cell fate. Information transfer depends on reproducibility of the gradients themselves, as well as the ability of cells to accurately measure and utilize morphogen concentrations in biologically meaning ways. In the following essay, these issues are discussed in context of the Bicoid gradient.

Combining protein and mRNA quantification to decipher transcriptional regulation

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

Bayesian model selection applied to the analysis of fluorescence correlation spectroscopy data of fluorescent proteins in vitro and in vivo

Fluorescence correlation spectroscopy (FCS) is a powerful technique to investigate molecular dynamics with single molecule sensitivity. In particular, in the life sciences it has found widespread application using fluorescent proteins as molecularly specific labels. However, FCS data analysis and interpretation using fluorescent proteins remains challenging due to typically low signal-to-noise ratio of FCS data and correlated noise in autocorrelated data sets. As a result, naive fitting procedures that ignore these important issues typically provide similarly good fits for multiple competing models without clear distinction of which model is preferred given the signal-to-noise ratio present in the data. Recently, we introduced a Bayesian model selection procedure to overcome this issue with FCS data analysis. The method accounts for the highly correlated noise that is present in FCS data sets and additionally penalizes model complexity to prevent over interpretation of FCS data. Here, we apply this procedure to evaluate FCS data from fluorescent proteins assayed in vitro and in vivo. Consistent with previous work, we demonstrate that model selection is strongly dependent on the signal-to-noise ratio of the measurement, namely, excitation intensity and measurement time, and is sensitive to saturation artifacts. Under fixed, low intensity excitation conditions, physical transport models can unambiguously be identified. However, at excitation intensities that are considered moderate in many studies, unwanted artifacts are introduced that result in nonphysical models to be preferred. We also determined the appropriate fitting models of a GFP tagged secreted signaling protein, Wnt3, in live zebrafish embryos, which is necessary for the investigation of Wnt3 expression and secretion in development. Bayes model selection therefore provides a robust procedure to determine appropriate transport and photophysical models for fluorescent proteins when appropriate models are provided, to help detect and eliminate experimental artifacts in solution, cells, and in living organisms.