Drosophila Embryogenesis Research Articles

Confinement promotes nematic alignment of spindle-shaped cells during Drosophila embryogenesis.

Tissue morphogenesis is often controlled by actomyosin networks pulling on adherens junctions (AJs), but junctional myosin levels vary. At an extreme, the Drosophila embryo amnioserosa forms a horseshoe-shaped strip of aligned, spindle-shaped cells lacking junctional myosin. What are the bases of amnioserosal cell interactions and alignment? Compared with surrounding tissue, we find that amnioserosal AJ continuity has lesser dependence on α-catenin, the mediator of AJ-actomyosin association, and greater dependence on Bazooka/Par-3, a junction-associated scaffold protein. Microtubule bundles also run along amnioserosal AJs and support their long-range curvilinearity. Amnioserosal confinement is apparent from partial overlap of its spindle-shaped cells, its outward bulging from surrounding tissue and from compressive stress detected within the amnioserosa. Genetic manipulations that alter amnioserosal confinement by surrounding tissue also result in amnioserosal cells losing alignment and gaining topological defects characteristic of nematically ordered systems. With Bazooka depletion, confinement by surrounding tissue appears to be relatively normal and amnioserosal cells align despite their AJ fragmentation. Overall, the fully elongated amnioserosa appears to form through tissue-autonomous generation of spindle-shaped cells that nematically align in response to confinement by surrounding tissue.

Organogenetic transcriptomes of the Drosophila embryo at single cell resolution.

To gain insight into the transcription programs activated during the formation of Drosophila larval structures, we carried out single cell RNA sequencing during two periods of Drosophila embryogenesis: stages 10-12, when most organs are first specified and initiate morphological and physiological specialization; and stages 13-16, when organs achieve their final mature architectures and begin to function. Our data confirm previous findings with regards to functional specialization of some organs - the salivary gland and trachea - and clarify the embryonic functions of another - the plasmatocytes. We also identify two early developmental trajectories in germ cells and uncover a potential role for proteolysis during germline stem cell specialization. We identify the likely cell type of origin for key components of the Drosophila matrisome and several commonly used Drosophila embryonic cell culture lines. Finally, we compare our findings with other recent related studies and with other modalities for identifying tissue-specific gene expression patterns. These data provide a useful community resource for identifying many new players in tissue-specific morphogenesis and functional specialization of developing organs.

Chromatin gene-gene loops support the cross-regulation of genes with related function

Chromatin loops between gene pairs have been observed in diverse contexts in both flies and vertebrates. Combining high-resolution Capture-C, DNA fluorescence in situ hybridization, and genetic perturbations, we dissect the functional role of three loops between genes with related function during Drosophila embryogenesis. By mutating the loop anchor (but not the gene) or the gene (but not loop anchor), we disentangle loop formation and gene expression and show that the 3D proximity of paralogous gene loci supports their co-regulation. Breaking the loop leads to either an attenuation or enhancement of expression and perturbs their relative levels of expression and cross-regulation. Although many loops appear constitutive across embryogenesis, their function can change in different developmental contexts. Taken together, our results indicate that chromatin gene-gene loops act as architectural scaffolds that can be used in different ways in different contexts to fine-tune the coordinated expression of genes with related functions and sustain their cross-regulation.

Ectopic Expression of Human, Mouse, and Fruit Fly Aryl Hydrocarbon Receptor in <i>Drosophila</i> <i>melanogaster</i> Tissues

The aryl hydrocarbon receptor (AHR) is a ligand-dependent transcription factor and its target genes play a fundamental role in detoxification, regulation of developmental processes, maintenance of homeostasis, as well as in the occurrence of oncological and autoimmune diseases and drug metabolism. The high conservatism of vertebrate AHR allowed us to study its functions in vivo using transformed Drosophila melanogaster fruit flies with human or mouse AHR gene and compare the ectopic effect of their expression with the expression of spineless gene, Drosophila AHR homologue. This work demonstrates for the first time that vertebrate AHR exhibits its functional activity in Drosophila embryogenesis, in leg imaginal discs and in somatic cells of female reproductive system in the absence of exogenous ligands.

BELMM: Bayesian model selection and random walk smoothing in time-series clustering.

Due to advances in measuring technology, many new phenotype, gene expression, and other omics time-course datasets are now commonly available. Cluster analysis may provide useful information about the structure of such data. In this work, we propose BELMM (Bayesian Estimation of Latent Mixture Models): a flexible framework for analysing, clustering, and modelling time-series data in a Bayesian setting. The framework is built on mixture modelling: first, the mean curves of the mixture components are assumed to follow random walk smoothing priors. Second, we choose the most plausible model and the number of mixture components using the Reversible-jump Markov chain Monte Carlo. Last, we assign the individual time series into clusters based on the similarity to the cluster-specific trend curves determined by the latent random walk processes. We demonstrate the use of fast and slow implementations of our approach on both simulated and real time-series data using widely available software R, Stan, and CU-MSDSp. The French mortality dataset is available at http://www.mortality.org, the Drosophila melanogaster embryogenesis gene expression data at https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE121160. Details on our simulated datasets are available in the Supplementary Material, and R scripts and a detailed tutorial on GitHub at https://github.com/ollisa/BELMM. The software CU-MSDSp is available on GitHub at https://github.com/jtchavisIII/CU-MSDSp.

Adhesion and Polarity protein distribution-regulates hexagon dominated plasma membrane organization in Drosophila blastoderm embryos.

Epithelial cells contain polarity complexes on the lateral membrane and are organized in a hexagon-dominated polygonal array. The mechanisms regulating the organization of polygonal architecture in metazoan embryogenesis are not completely understood. Drosophila embryogenesis enables mechanistic analysis of epithelial polarity formation and its impact on polygonal organization. The plasma membrane (PM) of syncytial Drosophila blastoderm embryos is organized as a polygonal array with pseudocleavage furrow formation during the almost synchronous cortical division cycles. We find that polygonal (PM) organization arises in the metaphase (MP) of division cycle 11, and hexagon dominance occurs with an increase in furrow length in the metaphase of cycle 12. There is a decrease in cell shape index in metaphase from cycles 11 to 13. This coincides with Drosophila E-cad (DE-cadherin) and Bazooka enrichment at the edges and the septin, Peanut at the vertices of the furrow. We further assess the role of polarity and adhesion proteins in pseudocleavage furrow formation and its organization as a polygonal array. We find that DE-cadherin depletion leads to decreased furrow length, loss of hexagon dominance, and increased cell shape index. Bazooka and Peanut depletion lead to decreased furrow length, delay in onset of hexagon dominance from cycle 12 to 13, and increased cell shape index. Hexagon dominance occurs with an increase in furrow length in cycle 13 and increased DE-cadherin, possibly due to the inhibition of endocytosis. We conclude that polarity protein recruitment and regulation of endocytic pathways enable pseudocleavage furrow stability and the formation of a hexagon-dominated polygon array.

Collection, Fixation, and Antibody Staining of Drosophila Embryos.

The Drosophila embryonic central nervous system has been used for decades as a model for understanding the genetic regulation of axon guidance and other aspects of neural development. Foundational studies using antibody staining to examine the embryonic ventral nerve cord in wild-type and mutant animals led to the discovery of evolutionarily conserved genes that regulate fundamental aspects of axon guidance, including midline crossing of axons. The development of the regular, segmentally repeating structure of axon pathways in the ventral nerve cord can illustrate basic principles of axon guidance to beginning students and can also be used by expert researchers to characterize new mutants, detect genetic interactions between known genes, and precisely quantify variations in gene function in engineered mutant lines. Here, we describe a protocol for collecting and fixing Drosophila embryos and visualizing axon pathways in the embryonic ventral nerve cord using immunofluorescence or immunohistochemical staining methods. As embryogenesis in Drosophila takes ∼24 h to complete, a 1-d collection yields embryos representing all stages of development from newly fertilized through ready-to-hatch larvae, allowing investigation of multiple developmental events within a single batch of collected embryos. The methods described in this protocol should be accessible to introductory laboratory courses as well as seasoned investigators in established research laboratories.

Deep Low-Shot Learning for Biological Image Classification and Visualization From Limited Training Samples.

Predictive modeling is useful but very challenging in biological image analysis due to the high cost of obtaining and labeling training data. For example, in the study of gene interaction and regulation in Drosophila embryogenesis, the analysis is most biologically meaningful when in situ hybridization (ISH) gene expression pattern images from the same developmental stage are compared. However, labeling training data with precise stages is very time-consuming even for developmental biologists. Thus, a critical challenge is how to build accurate computational models for precise developmental stage classification from limited training samples. In addition, identification and visualization of developmental landmarks are required to enable biologists to interpret prediction results and calibrate models. To address these challenges, we propose a deep two-step low-shot learning framework to accurately classify ISH images using limited training images. Specifically, to enable accurate model training on limited training samples, we formulate the task as a deep low-shot learning problem and develop a novel two-step learning approach, including data-level learning and feature-level learning. We use a deep residual network as our base model and achieve improved performance in the precise stage prediction task of ISH images. Furthermore, the deep model can be interpreted by computing saliency maps, which consists of pixel-wise contributions of an image to its prediction result. In our task, saliency maps are used to assist the identification and visualization of developmental landmarks. Our experimental results show that the proposed model can not only make accurate predictions but also yield biologically meaningful interpretations. We anticipate our methods to be easily generalizable to other biological image classification tasks with small training datasets. Our open-source code is available at https://github.com/divelab/lsl-fly.

Stress‐sensitive dynamics of miRNAs and Elba1 in Drosophila embryogenesis

Early‐life stress can result in life‐long effects that impact adult health and disease risk, but little is known about how such programming is established and maintained. Here, we show that such epigenetic memories can be initiated in the Drosophila embryo before the major wave of zygotic transcription, and higher‐order chromatin structures are established. An early short heat shock results in elevated levels of maternal miRNA and reduced levels of a subgroup of zygotic genes in stage 5 embryos. Using a Dicer‐1 mutant, we show that the stress‐induced decrease in one of these genes, the insulator‐binding factor Elba1, is dependent on functional miRNA biogenesis. Reduction in Elba1 correlates with the upregulation of early developmental genes and promotes a sustained weakening of heterochromatin in the adult fly as indicated by an increased expression of the PEV wm4h reporter. We propose that maternal miRNAs, retained in response to an early embryonic heat shock, shape the subsequent de novo heterochromatin establishment that occurs during early development via direct or indirect regulation of some of the earliest expressed genes, including Elba1.

Advanced maternal age has negative multigenerational impacts during Drosophila melanogaster embryogenesis

Increasing maternal age is commonly accompanied by decreased fitness in offspring. In Drosophila melanogaster, maternal senescence negatively affects multiple facets of offspring phenotype and fitness. These maternal effects are particularly large on embryonic viability. Identifying which embryonic stages are disrupted can indicate mechanisms of maternal effect senescence. Some maternal effects can also carry-over to subsequent generations. We examined potential multi- and transgenerational effects maternal senescence on embryonic development in two laboratory strains of D. melanogaster. We categorized the developmental stages of embryos from every combination of old and young mother, grandmother and great grandmother. We then modelled embryonic survival across the stages and compared these models among the multigenerational maternal age groups in order to identify which developmental processes were most sensitive to the effects of maternal effect senescence. Maternal effect senescence has negative multigenerational effects on multiple embryonic stages, indicating that maternal provisioning and, possibly epigenetics, but not mutation accumulation, contribute to decreased offspring survival. This study shows the large, early and multi-faceted nature of maternal effects senescence in an insect population.

Patterned mechanical feedback establishes a global myosin gradient

Morphogenesis, the coordinated execution of developmental programs that shape embryos, raises many fundamental questions at the interface between physics and biology. In particular, how the dynamics of active cytoskeletal processes are coordinated across the surface of entire embryos to generate global cell flows is poorly understood. Two distinct regulatory principles have been identified: genetic programs and dynamic response to mechanical stimuli. Despite progress, disentangling these two contributions remains challenging. Here, we combine in toto light sheet microscopy with genetic and optogenetic perturbations of tissue mechanics to examine theoretically predicted dynamic recruitment of non-muscle myosin II to cell junctions during Drosophila embryogenesis. We find dynamic recruitment has a long-range impact on global myosin configuration, and the rate of junction deformation sets the rate of myosin recruitment. Mathematical modeling and high frequency analysis reveal myosin fluctuations on junctions around a mean value set by mechanical feedback. Our model accounts for the early establishment of the global myosin pattern at 80% fidelity. Taken together our results indicate spatially modulated mechanical feedback as a key regulatory input in the establishment of long-range gradients of cytoskeletal configurations and global tissue flow patterns.

Manipulating the nature of embryonic mitotic waves

Early embryogenesis is characterized by rapid and synchronous cleavage divisions, which are often controlled by wave-like patterns of Cdk1 activity. Two mechanisms have been proposed for mitotic waves: sweep and trigger waves.1,2 The two mechanisms give rise to different wave speeds, dependencies on physical and molecular parameters, and spatial profiles of Cdk1 activity: upward sweeping gradients versus traveling wavefronts. Both mechanisms hinge on the transient bistability governing the cell cycle and are differentiated by the speed of the cell-cycle progression: sweep and trigger waves arise for rapid and slow drives, respectively. Here, using quantitative imaging of Cdk1 activity and theory, we illustrate that sweep waves are the dominant mechanism in Drosophila embryos and test two fundamental predictions on the transition from sweep to trigger waves. We demonstrate that sweep waves can be turned into trigger waves if the cell cycle is slowed down genetically or if significant delays in the cell-cycle progression are introduced across the embryo by altering nuclear density. Our genetic experiments demonstrate that Polo kinase is a major rate-limiting regulator of the blastoderm divisions, and genetic perturbations reducing its activity can induce the transition from sweep to trigger waves. Furthermore, we show that changes in temperature cause an essentially uniform slowdown of interphase and mitosis. That results in sweep waves being observed across a wide temperature range despite the cell-cycle durations being significantly different. Collectively, our combination of theory and experiments elucidates the nature of mitotic waves in Drosophila embryogenesis, their control mechanisms, and their mutual transitions.

Translation and natural selection of micropeptides from long non-canonical RNAs

Long noncoding RNAs (lncRNAs) are transcripts longer than 200 nucleotides but lacking canonical coding sequences. Apparently unable to produce peptides, lncRNA function seems to rely only on RNA expression, sequence and structure. Here, we exhaustively detect in-vivo translation of small open reading frames (small ORFs) within lncRNAs using Ribosomal profiling during Drosophila melanogaster embryogenesis. We show that around 30% of lncRNAs contain small ORFs engaged by ribosomes, leading to regulated translation of 100 to 300 micropeptides. We identify lncRNA features that favour translation, such as cistronicity, Kozak sequences, and conservation. For the latter, we develop a bioinformatics pipeline to detect small ORF homologues, and reveal evidence of natural selection favouring the conservation of micropeptide sequence and function across evolution. Our results expand the repertoire of lncRNA biochemical functions, and suggest that lncRNAs give rise to novel coding genes throughout evolution. Since most lncRNAs contain small ORFs with as yet unknown translation potential, we propose to rename them “long non-canonical RNAs”.

RNA Polymerase II "Pause" Prepares Promoters for Upcoming Transcription during Drosophila Development.

According to previous studies, during Drosophila embryogenesis, the recruitment of RNA polymerase II precedes active gene transcription. This work is aimed at exploring whether this mechanism is used during Drosophila metamorphosis. In addition, the composition of the RNA polymerase II “paused” complexes associated with promoters at different developmental stages are described in detail. For this purpose, we performed ChIP-Seq analysis using antibodies for various modifications of RNA polymerase II (total, Pol II CTD Ser5P, and Pol II CTD Ser2P) as well as for subunits of the NELF, DSIF, and PAF complexes and Brd4/Fs(1)h that control transcription elongation. We found that during metamorphosis, similar to mid-embryogenesis, the promoters were bound by RNA polymerase II in the “paused” state, preparing for activation at later stages of development. During mid-embryogenesis, RNA polymerase II in a “pause” state was phosphorylated at Ser5 and Ser2 of Pol II CTD and bound the NELF, DSIF, and PAF complexes, but not Brd4/Fs(1)h. During metamorphosis, the “paused” RNA polymerase II complex included Brd4/Fs(1)h in addition to NELF, DSIF, and PAF. The RNA polymerase II in this complex was phosphorylated at Ser5 of Pol II CTD, but not at Ser2. These results indicate that, during mid-embryogenesis, RNA polymerase II stalls in the “post-pause” state, being phosphorylated at Ser2 of Pol II CTD (after the stage of p-TEFb action). During metamorphosis, the “pause” mechanism is closer to classical promoter-proximal pausing and is characterized by a low level of Pol II CTD Ser2P.

R-loop Mapping and Characterization During Drosophila Embryogenesis Reveals Developmental Plasticity in R-loop Signatures

R-loops are involved in transcriptional regulation, DNA and histone post-translational modifications, genome replication and genome stability. To what extent R-loop abundance and genome-wide localization is actively regulated during metazoan embryogenesis is unknown. Drosophila embryogenesis provides a powerful system to address these questions due to its well-characterized developmental program, the sudden onset of zygotic transcription and available genome-wide data sets. Here, we measure the overall abundance and genome localization of R-loops in early and late-stage embryos relative to Drosophila cultured cells. We demonstrate that absolute R-loop levels change during embryogenesis and that RNaseH1 catalytic activity is critical for embryonic development. R-loop mapping by strand-specific DRIP-seq reveals that R-loop localization is plastic across development, both in the genes which form R-loops and where they localize relative to gene bodies. Importantly, these changes are not driven by changes in the transcriptional program. Negative GC skew and absolute changes in AT skew are associated with R-loop formation in Drosophila. Furthermore, we demonstrate that while some chromatin binding proteins and histone modifications such as H3K27me3 are associated with R-loops throughout development, other chromatin factors associated with R-loops in a developmental specific manner. Our findings highlight the importance and developmental plasticity of R-loops during Drosophila embryogenesis.

The Lateral Epidermis Actively Counteracts Pulling by the Amnioserosa During Dorsal Closure.

Dorsal closure is a prominent morphogenetic process during Drosophila embryogenesis, which involves two epithelial tissues, that is, the squamous amnioserosa and the columnar lateral epidermis. Non-muscle myosin II-driven constriction in the amnioserosa leads to a decrease in the apical surface area and pulls on the adjacent lateral epidermis, which subsequently moves dorsally. The pull by the amnioserosa becomes obvious in an elongation of the epidermal cells, especially of those in the first row. The contribution of the epidermal cell elongation has remained unclear to dorsal closure. Cell elongation may be a mere passive consequence or an active response to the pulling by the amnioserosa. Here, we found that the lateral epidermis actively responds. We analyzed tensions within tissues and cell junctions by laser ablation before and during dorsal closure, the elliptical and dorsal closure stages, respectively. Furthermore, we genetically and optochemically induced chronic and acute cell contraction, respectively. In this way, we found that tension in the epidermis increased during dorsal closure. A correspondingly increased tension was not observed at individual junctions, however. Junctional tension even decreased during dorsal closure in the epidermis. We strikingly observed a strong increase of the microtubule amount in the epidermis, while non-muscle myosin II increased in both tissues. Our data suggest that the epidermis actively antagonizes the pull from the amnioserosa during dorsal closure and the increased microtubules might help the epidermis bear part of the mechanical force.

A non-canonical Raf function is required for dorsal–ventral patterning during Drosophila embryogenesis

Proper embryonic development requires directional axes to pattern cells into embryonic structures. In Drosophila, spatially discrete expression of transcription factors determines the anterior to posterior organization of the early embryo, while the Toll and TGFβ signalling pathways determine the early dorsal to ventral pattern. Embryonic MAPK/ERK signaling contributes to both anterior to posterior patterning in the terminal regions and to dorsal to ventral patterning during oogenesis and embryonic stages. Here we describe a novel loss of function mutation in the Raf kinase gene, which leads to loss of ventral cell fates as seen through the loss of the ventral furrow, the absence of Dorsal/NFκB nuclear localization, the absence of mesoderm determinants Twist and Snail, and the expansion of TGFβ. Gene expression analysis showed cells adopting ectodermal fates much like loss of Toll signaling. Our results combine novel mutants, live imaging, optogenetics and transcriptomics to establish a novel role for Raf, that appears to be independent of the MAPK cascade, in embryonic patterning.

Condensation of the Drosophila nerve cord is oscillatory and depends on coordinated mechanical interactions

During development, organs reach precise shapes and sizes. Organ morphology is not always obtained through growth; a classic counterexample is the condensation of the nervous system during Drosophila embryogenesis. The mechanics underlying such condensation remain poorly understood. Here, we characterize the condensation of the embryonic ventral nerve cord (VNC) at both subcellular and tissue scales. This analysis reveals that condensation is not a unidirectional continuous process but instead occurs through oscillatory contractions. The VNC mechanical properties spatially and temporally vary, and forces along its longitudinal axis are spatially heterogeneous. We demonstrate that the process of VNC condensation is dependent on the coordinated mechanical activities of neurons and glia. These outcomes are consistent with a viscoelastic model of condensation, which incorporates time delays and effective frictional interactions. In summary, we have defined the progressive mechanics driving VNC condensation, providing insights into how a highly viscous tissue can autonomously change shape and size.

Cullin-5 mutants reveal collective sensing of the nucleocytoplasmic ratio in Drosophila embryogenesis

In most metazoans, early embryonic development is characterized by rapid division cycles that pause before gastrulation at the midblastula transition (MBT).1 These cleavage divisions are accompanied by cytoskeletal rearrangements that ensure proper nuclear positioning. However, the molecular mechanisms controlling nuclear positioning are not fully elucidated. In Drosophila, early embryogenesis unfolds in a multinucleated syncytium. Nuclei rapidly move across the anterior-posterior (AP) axis at cell cycles 4-6 in a process driven by actomyosin contractility and cytoplasmic flows.2,3 In shackleton (shkl) mutants, this axial spreading is impaired.4 Here, we show that shkl mutants carry mutations in the cullin-5 (cul-5) gene. Live imaging experiments show that Cul-5 is downstream of the cell cycle but is required for cortical actomyosin contractility. The nuclear spreading phenotype of cul-5 mutants can be rescued by reducing Src activity, suggesting that a major target of cul-5 is Src kinase. cul-5 mutants display gradients of nuclear density across the AP axis that we exploit to study cell-cycle control as a function of the N/C ratio. We found that the N/C ratio is sensed collectively in neighborhoods of about 100μm, and such collective sensing is required for a precise MBT, in which all the nuclei in the embryo pause their division cycle. Moreover, we found that the response to the N/C ratio is slightly graded along the AP axis. These two features can be linked to Cdk1 dynamics. Collectively, we reveal a new pathway controlling nuclear positioning and provide a dissection of how nuclear cyclesrespond to the N/C ratio.

Thermodynamics-based modeling reveals regulatory effects of indirect transcription factor-DNA binding

SummaryTranscription factors (TFs) influence gene expression by binding to DNA, yet experimental data suggests that they also frequently bind regulatory DNA indirectly by interacting with other DNA-bound proteins. Here, we used a data modeling approach to test if such indirect binding by TFs plays a significant role in gene regulation. We first incorporated regulatory function of indirectly bound TFs into a thermodynamics-based model for predicting enhancer-driven expression from its sequence. We then fit the new model to a rich data set comprising hundreds of enhancers and their regulatory activities during mesoderm specification in Drosophila embryogenesis and showed that the newly incorporated mechanism results in significantly better agreement with data. In the process, we derived the first sequence-level model of this extensively characterized regulatory program. We further showed that allowing indirect binding of a TF explains its localization at enhancers more accurately than with direct binding only. Our model also provided a simple explanation of how a TF may switch between activating and repressive roles depending on context.