Genome-wide measurement of spatial expression in patterning mutants of Drosophila melanogaster

Multi-level regulation of even-skipped stripes by the ubiquitous factor Zelda

The zinc-finger protein Zelda (Zld) is a key activator of zygotic transcription in early Drosophila embryos. Here, we study Zld-dependent regulation of the seven-striped pattern of the pair-rule gene even-skipped (eve). Individual stripes are regulated by discrete enhancers that respond to broadly distributed activators; stripe boundaries are formed by localized repressors encoded by the gap genes. The strongest effects of Zld are on stripes 2, 3 and 7, which are regulated by two enhancers in a 3.8 kb genomic fragment that includes the eve basal promoter. We show that Zld facilitates binding of the activator Bicoid and the gap repressors to this fragment, consistent with its proposed role as a pioneer protein. To test whether the effects of Zld are direct, we mutated all canonical Zld sites in the 3.8 kb fragment, which reduced expression but failed to phenocopy the abolishment of stripes caused by removing Zld in trans. We show that Zld also indirectly regulates the eve stripes by establishing specific gap gene expression boundaries, which provides the embryonic spacing required for proper stripe activation.

Chromatin accessibility in the Drosophila embryo is determined by transcription factor pioneering and enhancer activation

Chromatin accessibility is integral to the process by which transcription factors (TFs) read out cis-regulatory DNA sequences, but it is difficult to differentiate between TFs that drive accessibility and those that do not. Deep learning models that learn complex sequence rules provide an unprecedented opportunity to dissect this problem. Using zygotic genome activation in Drosophila as a model, we analyzed high-resolution TF binding and chromatin accessibility data with interpretable deep learning and performed genetic validation experiments. We identify a hierarchical relationship between the pioneer TF Zelda and the TFs involved in axis patterning. Zelda consistently pioneers chromatin accessibility proportional to motif affinity, whereas patterning TFs augment chromatin accessibility in sequence contexts where they mediate enhancer activation. We conclude that chromatin accessibility occurs in two tiers: one through pioneering, which makes enhancers accessible but not necessarily active, and the second when the correct combination of TFs leads to enhancer activation.

Genome-wide analysis of anterior-posterior mRNA localization in Stentor coeruleus reveals a role for the microtubule cytoskeleton

Cells have complex and beautiful structures that are important for their function, but understanding the molecular mechanisms that produce these structures is a challenging problem due to the gap in size scales between molecular interactions and cellular structures. The giant ciliate Stentor coeruleus is a unicellular model organism whose large size, reproducible structure, and ability to heal wounds and regenerate has historically allowed the formation of structure in a single cell to be addressed using methods of experimental embryology. Such studies have shown that specific cellular structures, such as the oral apparatus, always form in specific regions of the cell, which raises the question: what is the source of positional information within this organism? By analogy with embryonic development, in which localized mRNA is often used to mark position, we asked whether position along the anterior-posterior axis of Stentor might be marked by specific regionalized mRNAs. By physically bisecting cells and conducting half-cell RNA sequencing, we were able to identify sets of messages enriched in either the anterior or posterior half. We repeated this analysis in cells in which a set of longitudinal microtubule bundles running down the whole length of the cell, known as KM-fibers, were disrupted by RNAi of b-tubulin. We found that many messages either lost their regionalized distribution or switched to an opposite distribution, such that anterior-enriched messages in control became posterior-enriched in the RNAi cells, or vice versa. This study indicates that mRNA can be regionalized within a single giant cell and that microtubules may play a role, possibly by serving as tracks for the movement of the messages.

An integrated microfluidic device for multiplexed imaging of spatial gene expression patterns of Drosophila embryos

To reveal the underlying mechanism of the biological function of multicellular systems, it is important to obtain comprehensive spatial gene expression profiles. Among the emerging single-cell spatial-omics techniques, immunofluorescence (IF)-based iterative multiplexed imaging is a promising approach. However, the conventional method is usually costly, time-consuming, labor-intensive, and has low throughput. Moreover, it has yet to be demonstrated in intact multicellular organisms. Here, we developed an integrated microfluidic system to overcome these challenges for quantitatively measuring multiple protein profiles sequentially in situ in the same Drosophila embryo. We designed an array of hydrodynamic trapping sites to automatically capture over ten Drosophila embryos with orientation selectivity at more than 90% trapping rates. We also optimized the geometry of confinement and the on-chip IF protocol to achieve the same high signal-to-noise ratio as the off-chip traditional IF experiments. Moreover, we developed an efficient de-staining protocol by combining on-chip antibody stripping and fluorophore bleaching. Using the same secondary antibody to sequentially stain different genes, we confirmed that the de-stained genes have no detectable interference with the subsequently stained genes, and the gene expression profiles are preserved after multiple cycles of staining and de-staining processes. This preliminary test shows that our newly developed integrated microfluidic system can be a powerful tool for multiplexed imaging of Drosophila embryos. Our work opens a new avenue to design microfluidic chips for multicellular organisms and single-cell spatial-omics techniques.

Single-nucleus RNA-sequencing in pre-cellularization Drosophila melanogaster embryos

Our current understanding of the regulation of gene expression in the early Drosophila melanogaster embryo comes from observations of a few genes at a time, as with in situ hybridizations, or observation of gene expression levels without regards to patterning, as with RNA-sequencing. Single-nucleus RNA-sequencing however, has the potential to provide new insights into the regulation of gene expression for many genes at once while simultaneously retaining information regarding the position of each nucleus prior to dissociation based on patterned gene expression. In order to establish the use of single-nucleus RNA sequencing in Drosophila embryos prior to cellularization, here we look at gene expression in control and insulator protein, dCTCF, maternal null embryos during zygotic genome activation at nuclear cycle 14. We find that early embryonic nuclei can be grouped into distinct clusters according to gene expression. From both virtual and published in situ hybridizations, we also find that these clusters correspond to spatial regions of the embryo. Lastly, we provide a resource of candidate differentially expressed genes that might show local changes in gene expression between control and maternal dCTCF null nuclei with no detectable differential expression in bulk. These results highlight the potential for single-nucleus RNA-sequencing to reveal new insights into the regulation of gene expression in the early Drosophila melanogaster embryo.

Transcriptome-scale methods for uncovering subcellular RNA localization mechanisms

Across a variety of systems, thousands of RNAs are localized to specific subcellular locations. However, for the vast majority of these RNAs, the mechanisms that underlie their transport are unknown. Historically, these mechanisms were uncovered for a single transcript at a time by laboriously testing the ability of RNA fragments to direct transcript localization. Recently developed methods profile the content of subcellular transcriptomes using high-throughput sequencing, allowing the analysis of the localization of thousands of transcripts at once. By identifying commonalities shared among multiple localized transcripts, these methods have the potential to rapidly expand our understanding of RNA localization mechanisms.

When Bigger Is Better: 3D RNA Profiling of the Developing Head in the Catshark Scyliorhinus canicula

We report the adaptation of RNA tomography, a technique allowing spatially resolved, genome-wide expression profiling, to a species occupying a key phylogenetic position in gnathostomes, the catshark Scyliorhinus canicula. We focused analysis on head explants at an embryonic stage, shortly following neural tube closure and of interest for a number of developmental processes, including early brain patterning, placode specification or the establishment of epithalamic asymmetry. As described in the zebrafish, we have sequenced RNAs extracted from serial sections along transverse, horizontal and sagittal planes, mapped the data onto a gene reference taking advantage of the high continuity genome recently released in the catshark, and projected read counts onto a digital model of the head obtained by confocal microscopy. This results in the generation of a genome-wide 3D atlas, containing expression data for most protein-coding genes in a digital model of the embryonic head. The digital profiles obtained for candidate forebrain regional markers along antero-posterior, dorso-ventral and left-right axes reproduce those obtained by in situ hybridization (ISH), with expected relative organizations. We also use spatial autocorrelation and correlation as measures to analyze these data and show that they provide adequate statistical tools to extract novel expression information from the model. These data and tools allow exhaustive searches of genes exhibiting any predefined expression characteristic, such a restriction to a territory of interest, thus providing a reference for comparative analyses across gnathostomes. This methodology appears best suited to species endowed with large embryo or organ sizes and opens novel perspectives to a wide range of evo-devo model organisms, traditionally counter-selected on size criterion.

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

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

From whole-mount to single-cell spatial assessment of gene expression in 3D

Unravelling spatio-temporal patterns of gene expression is crucial to understanding core biological principles from embryogenesis to disease. Here we review emerging technologies, providing automated, high-throughput, spatially resolved quantitative gene expression data. Novel techniques expand on current benchmark protocols, expediting their incorporation into ongoing research. These approaches digitally reconstruct patterns of embryonic expression in three dimensions, and have successfully identified novel domains of expression, cell types, and tissue features. Such technologies pave the way for unbiased and exhaustive recapitulation of gene expression levels in spatial and quantitative terms, promoting understanding of the molecular origin of developmental defects, and improving medical diagnostics.

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.

Quantitative Comparison of the Anterior-Posterior Patterning System in the Embryos of Five Drosophila Species

Complex spatiotemporal gene expression patterns direct the development of the fertilized egg into an adult animal. Comparisons across species show that, in spite of changes in the underlying regulatory DNA sequence, developmental programs can be maintained across millions of years of evolution. Reciprocally, changes in gene expression can be used to generate morphological novelty. Distinguishing between changes in regulatory DNA that lead to changes in gene expression and those that do not is therefore a central goal of evolutionary developmental biology. Quantitative, spatially-resolved measurements of developmental gene expression patterns play a crucial role in this goal, enabling the detection of subtle phenotypic differences between species and the development of computations models that link the sequence of regulatory DNA to expression patterns. Here we report the generation of two atlases of cellular resolution gene expression measurements for the primary anterior-posterior patterning genes in Drosophila simulans and Drosophila virilis By combining these data sets with existing atlases for three other Drosophila species, we detect subtle differences in the gene expression patterns and dynamics driving the highly conserved axis patterning system and delineate inter-species differences in the embryonic morphology. These data sets will be a resource for future modeling studies of the evolution of developmental gene regulatory networks.

Dynamic multifactor hubs interact transiently with sites of active transcription in Drosophila embryos

The regulation of transcription requires the coordination of numerous activities on DNA, yet how transcription factors mediate these activities remains poorly understood. Here, we use lattice light-sheet microscopy to integrate single-molecule and high-speed 4D imaging in developing Drosophila embryos to study the nuclear organization and interactions of the key transcription factors Zelda and Bicoid. In contrast to previous studies suggesting stable, cooperative binding, we show that both factors interact with DNA with surprisingly high off-rates. We find that both factors form dynamic subnuclear hubs, and that Bicoid binding is enriched within Zelda hubs. Remarkably, these hubs are both short lived and interact only transiently with sites of active Bicoid-dependent transcription. Based on our observations, we hypothesize that, beyond simply forming bridges between DNA and the transcription machinery, transcription factors can organize other proteins into hubs that transiently drive multiple activities at their gene targets.

Editorial note

This article has been through an editorial process in which the authors decide how to respond to the issues raised during peer review. The Reviewing Editor's assessment is that all the issues have been addressed (see decision letter).

Spatially varying cis-regulatory divergence in Drosophila embryos elucidates cis-regulatory logic

Spatial patterning of gene expression is a key process in development, yet how it evolves is still poorly understood. Both cis- and trans-acting changes could participate in complex interactions, so to isolate the cis-regulatory component of patterning evolution, we measured allele-specific spatial gene expression patterns in D. melanogaster × simulans hybrid embryos. RNA-seq of cryo-sectioned slices revealed 66 genes with strong spatially varying allele-specific expression. We found that hunchback, a major regulator of developmental patterning, had reduced expression of the D. simulans allele specifically in the anterior tip of hybrid embryos. Mathematical modeling of hunchback cis-regulation suggested a candidate transcription factor binding site variant, which we verified as causal using CRISPR-Cas9 genome editing. In sum, even comparing morphologically near-identical species we identified surprisingly extensive spatial variation in gene expression, suggesting not only that development is robust to many such changes, but also that natural selection may have ample raw material for evolving new body plans via changes in spatial patterning.

Dosage effects in morphogenetic gradients of transcription factors: insights from a simple mathematical model

The classical Hill equation is the simplest way to model sharp transcriptional responses (TR) to changing concentrations of a regulator. Such steep sigmoidal transitions of the TR are often involved in the creation of boundaries or domains during development. Here, I use the Hill function to explore the response of a promoter to a gradient of a transcription factor (TF). Specifically, I examine how various sources of nonlinearity (such as that imposed by the number of TF-binding sites) and the nonlinearity of the gradient itself influence TR and modulate its robustness to variations in the dosage of morphogen. This qualitative analysis underscores how nonlinearities can modulate the robustness of TR and is a convenient way to convey concepts in a teaching setting.

Patterns of chromatin accessibility along the anterior-posterior axis in the early Drosophila embryo

As the Drosophila embryo transitions from the use of maternal RNAs to zygotic transcription, domains of open chromatin, with relatively low nucleosome density and specific histone marks, are established at promoters and enhancers involved in patterned embryonic transcription. However it remains unclear how regions of activity are established during early embryogenesis, and if they are the product of spatially restricted or ubiquitous processes. To shed light on this question, we probed chromatin accessibility across the anterior-posterior axis (A-P) of early Drosophila melanogaster embryos by applying a transposon based assay for chromatin accessibility (ATAC-seq) to anterior and posterior halves of hand-dissected, cellular blastoderm embryos. We find that genome-wide chromatin accessibility is highly similar between the two halves, with regions that manifest significant accessibility in one half of the embryo almost always accessible in the other half, even for promoters that are active in exclusively one half of the embryo. These data support previous studies that show that chromatin accessibility is not a direct result of activity, and point to a role for ubiquitous factors or processes in establishing chromatin accessibility at promoters in the early embryo. However, in concordance with similar works, we find that at enhancers active exclusively in one half of the embryo, we observe a significant skew towards greater accessibility in the region of their activity, highlighting the role of patterning factors such as Bicoid in this process.

Dense Bicoid hubs accentuate binding along the morphogen gradient

Morphogen gradients direct the spatial patterning of developing embryos; however, the mechanisms by which these gradients are interpreted remain elusive. Here we used lattice light-sheet microscopy to perform in vivo single-molecule imaging in early Drosophila melanogaster embryos of the transcription factor Bicoid that forms a gradient and initiates patterning along the anteroposterior axis. In contrast to canonical models, we observed that Bicoid binds to DNA with a rapid off rate throughout the embryo such that its average occupancy at target loci is on-rate-dependent. We further observed Bicoid forming transient "hubs" of locally high density that facilitate binding as factor levels drop, including in the posterior, where we observed Bicoid binding despite vanishingly low protein levels. We propose that localized modulation of transcription factor on rates via clustering provides a general mechanism to facilitate binding to low-affinity targets and that this may be a prevalent feature of other developmental transcription factors.