biniou (FoxF), a central component in a regulatory network controlling visceral mesoderm development and midgut morphogenesis in Drosophila

The Drosophila hematopoietic niche assembles through collective cell migration controlled by neighbor tissues and Slit-Robo signaling

Niches are often found in specific positions in tissues relative to the stem cells they support. Consistency of niche position suggests that placement is important for niche function. However, the complexity of most niches has precluded a thorough understanding of how their proper placement is established. To address this, we investigated the formation of a genetically tractable niche, the Drosophila Posterior Signaling Center (PSC), the assembly of which had not been previously explored. This niche controls hematopoietic progenitors of the lymph gland (LG). PSC cells were previously shown to be specified laterally in the embryo, but ultimately reside dorsally, at the LG posterior. Here, using live-imaging, we show that PSC cells migrate as a tight collective and associate with multiple tissues during their trajectory to the LG posterior. We find that Slit emanating from two extrinsic sources, visceral mesoderm and cardioblasts, is required for the PSC to remain a collective, and for its attachment to cardioblasts during migration. Without proper Slit-Robo signaling, PSC cells disperse, form aberrant contacts, and ultimately fail to reach their stereotypical position near progenitors. Our work characterizes a novel example of niche formation and identifies an extrinsic signaling relay that controls precise niche positioning.

Role of the Forkhead box family protein FOXF2 in the progression of solid tumor: systematic review

BACKGROUND

FOXF2 was reported to involve in a variety of biological behaviors that include the development of the central nervous system, tissue homeostasis, epithelia-mesenchymal interactions, regulation of embryonic development, and organogenesis.

PURPOSE

Understanding how FOXF2 influences the growth and development of cancer could provide valuable insights for researchers to develop novel therapeutic strategies.

RESULTS

In this review, we investigate the underlying impact of FOXF2 on tumor cells, including the transformation of cellular phenotype, capacity for migration, invasion, and proliferation, colonization of circulating cells, and formation of metastatic nodules. In addition, we discuss the molecular mechanisms of FOXF2 in different cancers, including hepatocellular, esophageal, breast, colon, lung, prostate gland, as well as its role in embryonic development.

CONCLUSION

FOXF2 is a gene encoding a forkhead transcription factor belonging to the Forkhead Box family. The protein functions by recruiting activation transcription factors and basic components to activate the transcription of genes that interact with the complex. This review provides an in-depth analysis of the FOXF2's function and pleiotropic roles in cancer development and progression.

The Drosophila hematopoietic niche assembles through collective cell migration controlled by neighbor tissues and Slit-Robo signaling

Niches are often found in specific positions in tissues relative to the stem cells they support. Consistency of niche position suggests that placement is important for niche function. However, the complexity of most niches has precluded a thorough understanding of how their proper placement is established. To address this, we investigated the formation of a genetically tractable niche, the Drosophila Posterior Signaling Center (PSC), the assembly of which had not been previously explored. This niche controls hematopoietic progenitors of the lymph gland (LG). PSC cells were previously shown to be specified laterally in the embryo, but ultimately reside dorsally, at the LG posterior. Here, using live-imaging, we show that PSC cells migrate as a tight collective and associate with multiple tissues during their trajectory to the LG posterior. We find that Slit emanating from two extrinsic sources, visceral mesoderm and cardioblasts, is required for the PSC to remain a collective, and for its attachment to cardioblasts during migration. Without proper Slit-Robo signaling, PSC cells disperse, form aberrant contacts, and ultimately fail to reach their stereotypical position near progenitors. Our work characterizes a novel example of niche formation and identifies an extrinsic signaling relay that controls precise niche positioning.

Novel Sex-Specific Genes and Diverse Interspecific Expression in the Antennal Transcriptomes of Ithomiine Butterflies

The olfactory sense is crucial for organisms, facilitating environmental recognition and interindividual communication. Ithomiini butterflies exemplify this importance not only because they rely strongly on olfactory cues for both inter- and intra-sexual behaviors, but also because they show convergent evolution of specialized structures within the antennal lobe, called macroglomerular complexes (MGCs). These structures, widely absent in butterflies, are present in moths where they enable heightened sensitivity to, and integration of, information from various types of pheromones. In this study, we investigate chemosensory evolution across six Ithomiini species and identify possible links between expression profiles and neuroanatomical. To enable this, we sequenced four new high-quality genome assemblies and six sex-specific antennal transcriptomes for three of these species with different MGC morphologies. With extensive genomic analyses, we found that the expression of antennal transcriptomes across species exhibit profound divergence, and identified highly expressed ORs, which we hypothesize may be associated to MGCs, as highly expressed ORs are absent in Methona, an Ithomiini lineage which also lacks MGCs. More broadly, we show how antennal sexual dimorphism is prevalent in both chemosensory genes and non-chemosensory genes, with possible relevance for behavior. As an example, we show how lipid-related genes exhibit consistent sexual dimorphism, potentially linked to lipid transport or host selection. In this study, we investigate the antennal chemosensory adaptations, suggesting a link between genetic diversity, ecological specialization, and sensory perception with the convergent evolution of MCGs. Insights into chemosensory gene evolution, expression patterns, and potential functional implications enhance our knowledge of sensory adaptations and sexual dimorphisms in butterflies, laying the foundation for future investigations into the genetic drivers of insect behavior, adaptation, and speciation.

Two sequential gene expression programs bridged by cell division support long-distance collective cell migration

The precise assembly of tissues and organs relies on spatiotemporal regulation of gene expression to coordinate the collective behavior of cells. In Drosophila embryos, the midgut musculature is formed through collective migration of caudal visceral mesoderm (CVM) cells, but how gene expression changes as cells migrate is not well understood. Here, we have focused on ten genes expressed in the CVM and the cis-regulatory sequences controlling their expression. Although some genes are continuously expressed, others are expressed only early or late during migration. Late expression relates to cell cycle progression, as driving string/Cdc25 causes earlier division of CVM cells and accelerates the transition to late gene expression. In particular, we found that the cell cycle effector transcription factor E2F1 is a required input for the late gene CG5080. Furthermore, whereas late genes are broadly expressed in all CVM cells, early gene transcripts are polarized to the anterior or posterior ends of the migrating collective. We show this polarization requires transcription factors Snail, Zfh1 and Dorsocross. Collectively, these results identify two sequential gene expression programs bridged by cell division that support long-distance directional migration of CVM cells.

Construction of single-cell cross-species chromatin accessibility landscapes with combinatorial-hybridization-based ATAC-seq

Despite recent advances in single-cell genomics, the lack of maps for single-cell candidate cis-regulatory elements (cCREs) in non-mammal species has limited our exploration of conserved regulatory programs across vertebrates and invertebrates. Here, we developed a combinatorial-hybridization-based method for single-cell assay for transposase-accessible chromatin using sequencing (scATAC-seq) named CH-ATAC-seq, enabling the construction of single-cell accessible chromatin landscapes for zebrafish, Drosophila, and earthworms (Eisenia andrei). By integrating scATAC censuses of humans, monkeys, and mice, we systematically identified 152 distinct main cell types and around 0.8 million cell-type-specific cCREs. Our analysis provided insights into the conservation of neural, muscle, and immune lineages across species, while epithelial cells exhibited a higher organ-origin heterogeneity. Additionally, a large-scale gene regulatory network (GRN) was constructed in four vertebrates by integrating scRNA-seq censuses. Overall, our study provides a valuable resource for comparative epigenomics, identifying the evolutionary conservation and divergence of gene regulation across different species.

LET-381/FoxF and its target UNC-30/Pitx2 specify and maintain the molecular identity of C. elegans mesodermal glia that regulate motor behavior

While most glial cell types in the central nervous system (CNS) arise from neuroectodermal progenitors, some, like microglia, are mesodermally derived. To understand mesodermal glia development and function, we investigated C. elegans GLR glia, which envelop the brain neuropil and separate it from the circulatory system cavity. Transcriptome analysis shows that GLR glia combine astrocytic and endothelial characteristics, which are relegated to separate cell types in vertebrates. Combined fate acquisition is orchestrated by LET-381/FoxF, a fate-specification/maintenance transcription factor also expressed in glia and endothelia of other animals. Among LET-381/FoxF targets, the UNC-30/Pitx2 transcription factor controls GLR glia morphology and represses alternative mesodermal fates. LET-381 and UNC-30 co-expression in naive cells is sufficient for GLR glia gene expression. GLR glia inactivation by ablation or let-381 mutation disrupts locomotory behavior and promotes salt-induced paralysis, suggesting brain-neuropil activity dysregulation. Our studies uncover mechanisms of mesodermal glia development and show that like neuronal differentiation, glia differentiation requires autoregulatory terminal selector genes that define and maintain the glial fate.

LET-381/FoxF and UNC-30/Pitx2 control the development of C. elegans mesodermal glia that regulate motor behavior

While most CNS glia arise from neuroectodermal progenitors, some, like microglia, are mesodermally derived. To understand mesodermal glia development and function, we investigated C. elegans GLR glia, which ensheath the brain neuropil and separate it from the circulatory-system cavity. Transcriptome analysis suggests GLR glia merge astrocytic and endothelial characteristics relegated to separate cell types in vertebrates. Combined fate acquisition is orchestrated by LET-381/FoxF, a fate-specification/maintenance transcription factor expressed in glia and endothelia of other animals. Among LET-381/FoxF targets, UNC-30/Pitx2 transcription factor controls GLR glia morphology and represses alternative mesodermal fates. LET-381 and UNC-30 co-expression in naïve cells is sufficient for GLR glia gene expression. GLR glia inactivation by ablation or let-381 mutation disrupts locomotory behavior and induces salt hypersensitivity, suggesting brain-neuropil activity dysregulation. Our studies uncover mechanisms of mesodermal glia development and show that like neurons, glia differentiation requires autoregulatory terminal selector genes that define and maintain the glial fate.

The RNF220 domain nuclear factor Teyrha-Meyrha (Tey) regulates the migration and differentiation of specific visceral and somatic muscles in Drosophila

Development of the visceral musculature of the Drosophila midgut encompasses a closely coordinated sequence of migration events of cells from the trunk and caudal visceral mesoderm that underlies the formation of the stereotypic orthogonal pattern of circular and longitudinal midgut muscles. Our study focuses on the last step of migration and morphogenesis of longitudinal visceral muscle precursors and shows that these multinucleated precursors utilize dynamic filopodial extensions to migrate in dorsal and ventral directions over the forming midgut tube. The establishment of maximal dorsoventral distances from one another, and anteroposterior alignments, lead to the equidistant coverage of the midgut with longitudinal muscle fibers. We identify Teyrha-Meyhra (Tey), a tissue-specific nuclear factor related to the RNF220 domain protein family, as a crucial regulator of this process of muscle migration and morphogenesis that is further required for proper differentiation of longitudinal visceral muscles. In addition, Tey is expressed in a single somatic muscle founder cell in each hemisegment, regulates the migration of this founder cell, and is required for proper pathfinding of its developing myotube to specific myotendinous attachment sites.

Aal-circRNA-407 regulates ovarian development of Aedes albopictus, a major arbovirus vector, via the miR-9a-5p/Foxl axis

Aedes albopictus shows a rapid global expansion and dramatic vectorial capacity for various arboviruses, thus posing a severe threat to global health. Although many noncoding RNAs have been confirmed to play functional roles in various biological processes in Ae. albopictus, the roles of circRNA remain a mystery. In the present study, we first performed high-throughput circRNA sequencing in Ae. albopictus. Then, we identified a cysteine desulfurase (CsdA) superfamily gene-originated circRNA, named aal-circRNA-407, which was the third most abundant circRNA in adult females and displayed a fat body highly expressed manifestation and blood feeding-dependent onset. SiRNA-mediated knockdown of circRNA-407 resulted in a decrease in the number of developing follicles and a reduction in follicle size post blood meal. Furthermore, we demonstrated that circRNA-407 can act as a sponge of aal-miR-9a-5p to promote the expression of its target gene Foxl and eventually regulate ovarian development. Our study is the first to report a functional circRNA in mosquitoes, expanding our current understanding of important biological roles in mosquitoes and providing an alternative genetic strategy for mosquito control.

A comprehensive study of arthropod and onychophoran Fox gene expression patterns

Fox genes represent an evolutionary old class of transcription factor encoding genes that evolved in the last common ancestor of fungi and animals. They represent key-components of multiple gene regulatory networks (GRNs) that are essential for embryonic development. Most of our knowledge about the function of Fox genes comes from vertebrate research, and for arthropods the only comprehensive gene expression analysis is that of the fly Drosophila melanogaster. For other arthropods, only selected Fox genes have been investigated. In this study, we provide the first comprehensive gene expression analysis of arthropod Fox genes including representative species of all main groups of arthropods, Pancrustacea, Myriapoda and Chelicerata. We also provide the first comprehensive analysis of Fox gene expression in an onychophoran species. Our data show that many of the Fox genes likely retained their function during panarthropod evolution highlighting their importance in development. Comparison with published data from other groups of animals shows that this high degree of evolutionary conservation often dates back beyond the last common ancestor of Panarthropoda.

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

Transcription 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.

•

Inclusion of indirect DNA binding of transcription factor improves enhancer function prediction

•

Context specific activating or repressive roles of TFs

•

Indirect binding improves fits to experimental TF-DNA binding data

•

Role of Tinman depends on its DNA-binding mode (direct or indirect)

Visceral mesoderm signaling regulates assembly position and function of the Drosophila testis niche

Tissue homeostasis often requires a properly placed niche to support stem cells. Morphogenetic processes that position a niche are just being described. For the Drosophila testis, we recently showed that pro-niche cells, specified at disparate positions during early gonadogenesis, must assemble into one collective at the anterior of the gonad. We now find that Slit and FGF signals emanating from adjacent visceral mesoderm regulate assembly. In response to signaling, niche cells express islet, which we find is also required for niche assembly. Without signaling, niche cells specified furthest from the anterior are unable to migrate, remaining dispersed. The function of such niches is severely disrupted, with niche cells evading cell cycle quiescence, compromised in their ability to signal the incipient stem cell pool, and failing to orient stem cell divisions properly. Our work identifies both extrinsic signaling and intrinsic responses required for proper assembly and placement of the testis niche.

Simultaneous cellular and molecular phenotyping of embryonic mutants using single-cell regulatory trajectories

Developmental progression and cellular diversity are largely driven by transcription factors (TFs); yet, characterizing their loss-of-function phenotypes remains challenging and often disconnected from their underlying molecular mechanisms. Here, we combine single-cell regulatory genomics with loss-of-function mutants to jointly assess both cellular and molecular phenotypes. Performing sci-ATAC-seq at eight overlapping time points during Drosophila mesoderm development could reconstruct the developmental trajectories of all major muscle types and reveal the TFs and enhancers involved. To systematically assess mutant phenotypes, we developed a single-nucleus genotyping strategy to process embryo pools of mixed genotypes. Applying this to four TF mutants could identify and quantify their characterized phenotypes de novo and discover new ones, while simultaneously revealing their regulatory input and mode of action. Our approach is a general framework to dissect the functional input of TFs in a systematic, unbiased manner, identifying both cellular and molecular phenotypes at a scale and resolution that has not been feasible before.

Simulation and in vivo experimentation predict AdamTS-A location of function during caudal visceral mesoderm (CVM) migration in Drosophila

BACKGROUND

Caudal visceral mesoderm (CVM) cells migrate as a loose collective along the trunk visceral mesoderm (TVM) and are surrounded by extracellular matrix (ECM). In this study, we examined how one extracellular protease, AdamTS-A, facilitates CVM migration.

RESULTS

A comparison of mathematical simulation to experimental results suggests that location of AdamTS-A action in CVM cells is on the sides of the cell not in contact with the TVM, predominantly at the CVM-ECM interface. CVM migration from a top-down view showed CVM cells migrating along the outside of the TVM substrate in the absence of AdamTS-A. Moreover, over-expression of AdamTS-A resulted in similar, but milder, mis-migration of the CVM. These results contrast with the salivary gland where AdamTS-A is proposed to cleave connections at the trailing edge of migrating cells. Subcellular localization of GFP-tagged AdamTS-A suggests this protease is not limited to functioning at the trailing edge of CVM cells.

CONCLUSION

Using both in vivo experimentation and mathematical simulations, we demonstrated that AdamTS-A cleaves connections between CVM cells and the ECM on all sides not attached to the TVM. Clearly, AdamTS-A has a more expansive role around the entire cell in cleaving cell-ECM attachments in cells migrating as a loose collective. This article is protected by copyright. All rights reserved.

DamID transcriptional profiling identifies the Snail/Scratch transcription factor Kahuli as an Alk target in the Drosophila visceral mesoderm

Development of the Drosophila visceral muscle depends on Anaplastic Lymphoma Kinase (Alk) receptor tyrosine kinase (RTK) signaling, which specifies founder cells (FCs) in the circular visceral mesoderm (VM). Although Alk activation by its ligand Jelly Belly (Jeb) is well characterized, few target molecules have been identified. Here, we used targeted DamID (TaDa) to identify Alk targets in embryos overexpressing Jeb versus embryos with abrogated Alk activity, revealing differentially expressed genes, including the Snail/Scratch family transcription factor Kahuli (Kah). We confirmed Kah mRNA and protein expression in the VM, and identified midgut constriction defects in Kah mutants similar to those of pointed (pnt). ChIP and RNA-Seq data analysis defined a Kah target-binding site similar to that of Snail, and identified a set of common target genes putatively regulated by Kah and Pnt during midgut constriction. Taken together, we report a rich dataset of Alk-responsive loci in the embryonic VM and functionally characterize the role of Kah in the regulation of embryonic midgut morphogenesis.

Extremely rapid and reversible optogenetic perturbation of nuclear proteins in living embryos

Many developmental regulators have complex and context-specific roles in different tissues and stages, making the dissection of their function extremely challenging. As regulatory processes often occur within minutes, perturbation methods that match these dynamics are needed. Here, we present the improved light-inducible nuclear export system (iLEXY), an optogenetic loss-of-function approach that triggers translocation of proteins from the nucleus to the cytoplasm. By introducing a series of mutations, we substantially increased LEXY's efficiency and generated variants with different recovery times. iLEXY enables rapid (t1/2 < 30 s), efficient, and reversible nuclear protein depletion in embryos, and is generalizable to proteins of diverse sizes and functions. Applying iLEXY to the Drosophila master regulator Twist, we phenocopy loss-of-function mutants, precisely map the Twist-sensitive embryonic stages, and investigate the effects of timed Twist depletions. Our results demonstrate the power of iLEXY to dissect the function of pleiotropic factors during embryogenesis with unprecedented temporal precision.

FGF signaling acts on different levels of mesoderm development within Spiralia

FGF signaling is involved in mesoderm induction in members of deuterostomes (e.g. tunicates, hemichordates), but not in flies and nematodes, in which it has a role in mesoderm patterning and migration. However, we need comparable studies in other protostome taxa in order to decipher whether this mesoderm-inducing function of FGF extends beyond the lineage of deuterostomes. Here, we investigated the role of FGF signaling in mesoderm development in three species of lophophorates, a clade within the protostome group Spiralia. Our gene expression analyses show that the mesodermal molecular patterning is conserved between brachiopods and phoronids, but the spatial and temporal recruitment of transcription factors differs significantly. Moreover, the use of the inhibitor SU5402 demonstrates that FGF signaling is involved in different steps of mesoderm development, as well as in morphogenetic movements of gastrulation and axial elongation. Our findings suggest that the mesoderm-inducing role of FGF extends beyond the group of deuterostomes.

Multiple Niche Compartments Orchestrate Stepwise Germline Stem Cell Progeny Differentiation

The niche controls stem cell self-renewal and progenitor differentiation for maintaining adult tissue homeostasis in various organisms. However, it remains unclear whether the niche is compartmentalized to control stem cell self-renewal and stepwise progeny differentiation. In the Drosophila ovary, inner germarial sheath (IGS) cells form a niche for controlling germline stem cell (GSC) progeny differentiation. In this study, we have identified four IGS subpopulations, which form linearly arranged niche compartments for controlling GSC maintenance and multi-step progeny differentiation. Single-cell analysis of the adult ovary has identified four IGS subpopulations (IGS1-IGS4), the identities and cellular locations of which have been further confirmed by fluorescent in situ hybridization. IGS1 and IGS2 physically interact with GSCs and mitotic cysts to control GSC maintenance and cyst formation, respectively, whereas IGS3 and IGS4 physically interact with 16-cell cysts to regulate meiosis, oocyte development, and cyst morphological change. Finally, one follicle cell progenitor population has also been transcriptionally defined for facilitating future studies on follicle stem cell regulation. Therefore, this study has structurally revealed that the niche is organized into multiple compartments for orchestrating stepwise adult stem cell development and has also provided useful resources and tools for further functional characterization of the niche in the future.

The lateral plate mesoderm

The lateral plate mesoderm (LPM) forms the progenitor cells that constitute the heart and cardiovascular system, blood, kidneys, smooth muscle lineage and limb skeleton in the developing vertebrate embryo. Despite this central role in development and evolution, the LPM remains challenging to study and to delineate, owing to its lineage complexity and lack of a concise genetic definition. Here, we outline the processes that govern LPM specification, organization, its cell fates and the inferred evolutionary trajectories of LPM-derived tissues. Finally, we discuss the development of seemingly disparate organ systems that share a common LPM origin.

Summary: The lateral plate mesoderm is the origin of several major cell types and organ systems in the vertebrate body plan. How this mesoderm territory emerges and partitions into its downstream fates provides clues about vertebrate development and evolution.

Collective Migrations of Drosophila Embryonic Trunk and Caudal Mesoderm-Derived Muscle Precursor Cells

Mesoderm migration in the Drosophila embryo is a highly conserved, complex process that is required for the formation of specialized tissues and organs, including the somatic and visceral musculature. In this FlyBook chapter, we will compare and contrast the specification and migration of cells originating from the trunk and caudal mesoderm. Both cell types engage in collective migrations that enable cells to achieve new positions within developing embryos and form distinct tissues. To start, we will discuss specification and early morphogenetic movements of the presumptive mesoderm, then focus on the coordinate movements of the two subtypes trunk mesoderm and caudal visceral mesoderm, ending with a comparison of these processes including general insights gained through study.

The Role of Chromatin Accessibility in cis-Regulatory Evolution

Transcription factor (TF) binding is determined by sequence as well as chromatin accessibility. Although the role of accessibility in shaping TF-binding landscapes is well recorded, its role in evolutionary divergence of TF binding, which in turn can alter cis-regulatory activities, is not well understood. In this work, we studied the evolution of genome-wide binding landscapes of five major TFs in the core network of mesoderm specification, between Drosophila melanogaster and Drosophila virilis, and examined its relationship to accessibility and sequence-level changes. We generated chromatin accessibility data from three important stages of embryogenesis in both Drosophila melanogaster and Drosophila virilis and recorded conservation and divergence patterns. We then used multivariable models to correlate accessibility and sequence changes to TF-binding divergence. We found that accessibility changes can in some cases, for example, for the master regulator Twist and for earlier developmental stages, more accurately predict binding change than is possible using TF-binding motif changes between orthologous enhancers. Accessibility changes also explain a significant portion of the codivergence of TF pairs. We noted that accessibility and motif changes offer complementary views of the evolution of TF binding and developed a combined model that captures the evolutionary data much more accurately than either view alone. Finally, we trained machine learning models to predict enhancer activity from TF binding and used these functional models to argue that motif and accessibility-based predictors of TF-binding change can substitute for experimentally measured binding change, for the purpose of predicting evolutionary changes in enhancer activity.

Combinatorial chromatin dynamics foster accurate cardiopharyngeal fate choices

During embryogenesis, chromatin accessibility profiles control lineage-specific gene expression by modulating transcription, thus impacting multipotent progenitor states and subsequent fate choices. Subsets of cardiac and pharyngeal/head muscles share a common origin in the cardiopharyngeal mesoderm, but the chromatin landscapes that govern multipotent progenitors competence and early fate choices remain largely elusive. Here, we leveraged the simplicity of the chordate model Ciona to profile chromatin accessibility through stereotyped transitions from naive Mesp+ mesoderm to distinct fate-restricted heart and pharyngeal muscle precursors. An FGF-Foxf pathway acts in multipotent progenitors to establish cardiopharyngeal-specific patterns of accessibility, which govern later heart vs. pharyngeal muscle-specific expression profiles, demonstrating extensive spatiotemporal decoupling between early cardiopharyngeal enhancer accessibility and late cell-type-specific activity. We found that multiple cis-regulatory elements, with distinct chromatin accessibility profiles and motif compositions, are required to activate Ebf and Tbx1/10, two key determinants of cardiopharyngeal fate choices. We propose that these 'combined enhancers' foster spatially and temporally accurate fate choices, by increasing the repertoire of regulatory inputs that control gene expression, through either accessibility and/or activity.

The Hox transcription factor Ubx stabilizes lineage commitment by suppressing cellular plasticity in Drosophila

During development cells become restricted in their differentiation potential by repressing alternative cell fates, and the Polycomb complex plays a crucial role in this process. However, how alternative fate genes are lineage-specifically silenced is unclear. We studied Ultrabithorax (Ubx), a multi-lineage transcription factor of the Hox class, in two tissue lineages using sorted nuclei and interfered with Ubx in mesodermal cells. We find that depletion of Ubx leads to the de-repression of genes normally expressed in other lineages. Ubx silences expression of alternative fate genes by retaining the Polycomb Group protein Pleiohomeotic at Ubx targeted genomic regions, thereby stabilizing repressive chromatin marks in a lineage-dependent manner. Our study demonstrates that Ubx stabilizes lineage choice by suppressing the multipotency encoded in the genome via its interaction with Pho. This mechanism may explain why the Hox code is maintained throughout the lifecycle, since it could set a block to transdifferentiation in adult cells.

Distinct roles and requirements for Ras pathway signaling in visceral versus somatic muscle founder specification

Pleiotropic signaling pathways must somehow engender specific cellular responses. In the Drosophila mesoderm, Ras pathway signaling specifies muscle founder cells from among the broader population of myoblasts. For somatic muscles, this is an inductive process mediated by the ETS-domain downstream Ras effectors Pointed and Aop (Yan). We demonstrate here that for the circular visceral muscles, despite superficial similarities, a significantly different specification mechanism is at work. Not only is visceral founder cell specification not dependent on Pointed or Aop, but Ras pathway signaling in its entirety can be bypassed. Our results show that de-repression, not activation, is the predominant role of Ras signaling in the visceral mesoderm and that, accordingly, Ras signaling is not required in the absence of repression. The key repressor acts downstream of the transcription factor Lame duck and is likely a member of the ETS transcription factor family. Our findings fit with a growing body of data that point to a complex interplay between the Ras pathway, ETS transcription factors, and enhancer binding as a crucial mechanism for determining unique responses to Ras signaling.

foxF-1 Controls Specification of Non-body Wall Muscle and Phagocytic Cells in Planarians

Planarians are flatworms capable of regenerating any missing body part in a process requiring stem cells and positional information. Muscle is a major source of planarian positional information and consists of several types of fibers with distinct regulatory roles in regeneration. The transcriptional regulatory programs used to specify different muscle fibers are poorly characterized. Using single-cell RNA sequencing, we define the transcriptomes of planarian dorsal-ventral muscle (DVM), intestinal muscle (IM), and pharynx muscle. This analysis identifies foxF-1, which encodes a broadly conserved Fox-family transcription factor, as a master transcriptional regulator of all non-body wall muscle. The transcription factors encoded by nk4 and gata4/5/6-2 specify two different subsets of DVM, lateral and medial, respectively, whereas gata4/5/6-3 specifies IM. These muscle types all express planarian patterning genes. Both lateral and medial DVM are required for medial-lateral patterning in regeneration, whereas medial DVM and IM have a role in maintaining and regenerating intestine morphology. In addition to the role in muscle, foxF-1 is required for the specification of multiple cell types with transcriptome similarities, including high expression levels of cathepsin genes. These cells include pigment cells, glia, and several other cells with unknown function. cathepsin⁺ cells phagocytose E. coli, suggesting these are phagocytic cells. In conclusion, we describe a regulatory program for planarian muscle cell subsets and phagocytic cells, both driven by foxF-1. FoxF proteins specify different mesoderm-derived tissues in other organisms, suggesting that FoxF regulates formation of an ancient and broadly conserved subset of mesoderm derivatives in the Bilateria.

Gene Regulation of BMP Ligands in Drosophila

Drosophila is a valuable system to study bone morphogenetic proteins (BMPs) due to the high functional conservation of the pathway and the molecular genetic tools available. Drosophila has three BMP ligands, decapentaplegic (BMP2/4), screw, and glass bottom boat (BMP5/6/7/8). Of these genes, the transcriptional regulation of decapentaplegic has been studied, and some of the enhancers directing its spatially specific gene expression have been described. These analyses have used many of the standard tools of molecular biology, but a valuable method of analysis often used in Drosophila is the creation of patches of mutant tissue at any stage and in any location by induced somatic recombination. The ability to create transgenic flies and manipulate the Drosophila genome with recombinases is well established. This method can be used to evaluate the requirements for specific transcription factors to act on enhancer elements in vivo, in stage- and tissue-specific manners. The yeast FLP/FRT recombination system facilitates experiments to interrogate loss- or gain-of-function for transcription factor activity on known enhancers. This chapter will outline the necessary steps to create the tools needed and conduct somatic cell recombination experiments to interrogate the function of transcription factors on BMP enhancers.

Odd-Paired: The Drosophila Zic Gene

Zinc finger in the cerebellum (Zic) proteins are a family of transcription factors with multiple roles during development, particularly in neural tissues. The founding member of the Zic family is the Drosophila odd-paired (opa) gene. The Opa protein has a DNA binding domain containing five Cys2His2-type zinc fingers and has been shown to act as a sequence-specific DNA binding protein. Opa has significant homology to mammalian Zic1, Zic2, and Zic3 within the zinc finger domain and in two other conserved regions outside that domain. opa was initially identified as a pair-rule gene, part of the hierarchy of genes that establish the segmental body plan of the early Drosophila embryo. However, its wide expression pattern during embryogenesis indicates it plays additional roles. Embryos deficient in opa die before hatching with aberrant segmentation but also with defects in larval midgut formation. Post-embryonically, opa plays important roles in adult head development and circadian rhythm. Based on extensive neural expression, opa is predicted to be involved in many aspects of neural development and behavior, like other proteins of the Zic family. Consensus DNA binding sites have been identified for Opa and have been shown to activate transcription in vivo. However, there is evidence Opa may serve as a transcriptional regulator in the absence of direct DNA binding, as has been seen for other Zic proteins.

Uncoupling evolutionary changes in DNA sequence, transcription factor occupancy and enhancer activity

Sequence variation within enhancers plays a major role in both evolution and disease, yet its functional impact on transcription factor (TF) occupancy and enhancer activity remains poorly understood. Here, we assayed the binding of five essential TFs over multiple stages of embryogenesis in two distant Drosophila species (with 1.4 substitutions per neutral site), identifying thousands of orthologous enhancers with conserved or diverged combinatorial occupancy. We used these binding signatures to dissect two properties of developmental enhancers: (1) potential TF cooperativity, using signatures of co-associations and co-divergence in TF occupancy. This revealed conserved combinatorial binding despite sequence divergence, suggesting protein-protein interactions sustain conserved collective occupancy. (2) Enhancer in-vivo activity, revealing orthologous enhancers with conserved activity despite divergence in TF occupancy. Taken together, we identify enhancers with diverged motifs yet conserved occupancy and others with diverged occupancy yet conserved activity, emphasising the need to functionally measure the effect of divergence on enhancer activity.

Transcription Factor Forkhead Regulates Expression of Antimicrobial Peptides in the Tobacco Hornworm, Manduca sexta

Antimicrobial peptides (AMPs) play an important role in defense against microbial infections in insects. Expression of AMPs is regulated mainly by NF-κB factors Dorsal, Dif and Relish. Our previous study showed that both NF-κB and GATA-1 factors are required for activation of moricin promoter in the tobacco hornworm, Manduca sexta, and a 140-bp region in the moricin promoter contains binding sites for additional transcription factors. In this study, we identified three forkhead (Fkh)-binding sites in the 140-bp region of the moricin promoter and several Fkh-binding sites in the lysozyme promoter, and demonstrated that Fkh-binding sites are required for activation of both moricin and lysozyme promoters by Fkh factors. In addition, we found that Fkh mRNA was undetectable in Drosophila S2 cells, and M. sexta Fkh (MsFkh) interacted with Relish-Rel-homology domain (RHD) but not with Dorsal-RHD. Dual luciferase assays with moricin mutant promoters showed that co-expression of MsFkh with Relish-RHD did not have an additive effect on the activity of moricin promoter, suggesting that MsFkh and Relish regulate moricin activation independently. Our results suggest that insect AMPs can be activated by Fkh factors under non-infectious conditions, which may be important for protection of insects from microbial infection during molting and metamorphosis.

Cleavage modification did not alter blastomere fates during bryozoan evolution

BACKGROUND

Stereotypic cleavage patterns play a crucial role in cell fate determination by precisely positioning early embryonic blastomeres. Although misplaced cell divisions can alter blastomere fates and cause embryonic defects, cleavage patterns have been modified several times during animal evolution. However, it remains unclear how evolutionary changes in cleavage impact the specification of blastomere fates. Here, we analyze the transition from spiral cleavage - a stereotypic pattern remarkably conserved in many protostomes - to a biradial cleavage pattern, which occurred during the evolution of bryozoans.

RESULTS

Using 3D-live imaging time-lapse microscopy (4D-microscopy), we characterize the cell lineage, MAPK signaling, and the expression of 16 developmental genes in the bryozoan Membranipora membranacea. We found that the molecular identity and the fates of early bryozoan blastomeres are similar to the putative homologous blastomeres in spiral-cleaving embryos.

CONCLUSIONS

Our work suggests that bryozoans have retained traits of spiral development, such as the early embryonic fate map, despite the evolution of a novel cleavage geometry. These findings provide additional support that stereotypic cleavage patterns can be modified during evolution without major changes to the molecular identity and fate of embryonic blastomeres.

Comparative analysis of gene expression profiles for several migrating cell types identifies cell migration regulators

Cell migration is an instrumental process that ensures cells are properly positioned to support the specification of distinct tissue types during development. To provide insight, we used fluorescence activated cell sorting (FACS) to isolate two migrating cell types from the Drosophila embryo: caudal visceral mesoderm (CVM) cells, precursors of longitudinal muscles of the gut, and hemocytes (HCs), the Drosophila equivalent of blood cells. ~350 genes were identified from each of the sorted samples using RNA-seq, and in situ hybridization was used to confirm expression within each cell type or, alternatively, within other interacting, co-sorted cell types. To start, the two gene expression profiling datasets were compared to identify cell migration regulators that are potentially generally-acting. 73 genes were present in both CVM cell and HC gene expression profiles, including the transcription factor zinc finger homeodomain-1 (zfh1). Comparisons with gene expression profiles of Drosophila border cells that migrate during oogenesis had a more limited overlap, with only the genes neyo (neo) and singed (sn) found to be expressed in border cells as well as CVM cells and HCs, respectively. Neo encodes a protein with Zona pellucida domain linked to cell polarity, while sn encodes an actin binding protein. Tissue specific RNAi expression coupled with live in vivo imaging was used to confirm cell-autonomous roles for zfh1 and neo in supporting CVM cell migration, whereas previous studies had demonstrated a role for Sn in supporting HC migration. In addition, comparisons were made to migrating cells from vertebrates. Seven genes were found expressed by chick neural crest cells, CVM cells, and HCs including extracellular matrix (ECM) proteins and proteases. In summary, we show that genes shared in common between CVM cells, HCs, and other migrating cell types can help identify regulators of cell migration. Our analyses show that neo in addition to zfh1 and sn studied previously impact cell migration. This study also suggests that modification of the extracellular milieu may be a fundamental requirement for cells that undergo cell streaming migratory behaviors.

The Zic family homologue Odd-paired regulates Alk expression in Drosophila

The Anaplastic Lymphoma Kinase (Alk) receptor tyrosine kinase (RTK) plays a critical role in the specification of founder cells (FCs) in the Drosophila visceral mesoderm (VM) during embryogenesis. Reporter gene and CRISPR/Cas9 deletion analysis reveals enhancer regions in and upstream of the Alk locus that influence tissue-specific expression in the amnioserosa (AS), the VM and the epidermis. By performing high throughput yeast one-hybrid screens (Y1H) with a library of Drosophila transcription factors (TFs) we identify Odd-paired (Opa), the Drosophila homologue of the vertebrate Zic family of TFs, as a novel regulator of embryonic Alk expression. Further characterization identifies evolutionarily conserved Opa-binding cis-regulatory motifs in one of the Alk associated enhancer elements. Employing Alk reporter lines as well as CRISPR/Cas9-mediated removal of regulatory elements in the Alk locus, we show modulation of Alk expression by Opa in the embryonic AS, epidermis and VM. In addition, we identify enhancer elements that integrate input from additional TFs, such as Binou (Bin) and Bagpipe (Bap), to regulate VM expression of Alk in a combinatorial manner. Taken together, our data show that the Opa zinc finger TF is a novel regulator of embryonic Alk expression.

The Alk receptor tyrosine kinase is employed repeatedly during Drosophila development to drive signaling events in a variety of tissues. The spatial and temporal expression pattern of the Alk gene is tightly regulated. Identifying factors that influence the expression of Alk is important to better understand how Alk signaling is controlled. In this paper we characterize cis-regulatory sequences in the Alk locus and the transcription factors that bind them to govern Alk expression in the Drosophila embryo. Using a robotic protein-DNA interaction assay, we identified the Zic family transcription factor Odd-paired as a factor that binds to regulatory elements in the Alk locus. Binding of Odd-paired to Alkcis-regulatory elements varies spatially, revealing a requirement for additional transcription factors such as the NK3 and FoxF orthologues Bagpipe and Biniou in a subset of Alk-expressing tissues. Our findings provide new insight into the dynamics underlying temporal and spatial regulation of the Alk receptor during embryogenesis.

The evolutionary origin of bilaterian smooth and striated myocytes

The dichotomy between smooth and striated myocytes is fundamental for bilaterian musculature, but its evolutionary origin is unsolved. In particular, interrelationships of visceral smooth muscles remain unclear. Absent in fly and nematode, they have not yet been characterized molecularly outside vertebrates. Here, we characterize expression profile, ultrastructure, contractility and innervation of the musculature in the marine annelid Platynereis dumerilii and identify smooth muscles around the midgut, hindgut and heart that resemble their vertebrate counterparts in molecular fingerprint, contraction speed and nervous control. Our data suggest that both visceral smooth and somatic striated myocytes were present in the protostome-deuterostome ancestor and that smooth myocytes later co-opted the striated contractile module repeatedly – for example, in vertebrate heart evolution. During these smooth-to-striated myocyte conversions, the core regulatory complex of transcription factors conveying myocyte identity remained unchanged, reflecting a general principle in cell type evolution.

DOI:http://dx.doi.org/10.7554/eLife.19607.001

Qualitative Dynamical Modelling Can Formally Explain Mesoderm Specification and Predict Novel Developmental Phenotypes

Given the complexity of developmental networks, it is often difficult to predict the effect of genetic perturbations, even within coding genes. Regulatory factors generally have pleiotropic effects, exhibit partially redundant roles, and regulate highly interconnected pathways with ample cross-talk. Here, we delineate a logical model encompassing 48 components and 82 regulatory interactions involved in mesoderm specification during Drosophila development, thereby providing a formal integration of all available genetic information from the literature. The four main tissues derived from mesoderm correspond to alternative stable states. We demonstrate that the model can predict known mutant phenotypes and use it to systematically predict the effects of over 300 new, often non-intuitive, loss- and gain-of-function mutations, and combinations thereof. We further validated several novel predictions experimentally, thereby demonstrating the robustness of model. Logical modelling can thus contribute to formally explain and predict regulatory outcomes underlying cell fate decisions.

Genome-Wide Approaches to Drosophila Heart Development

The development of the dorsal vessel in Drosophila is one of the first systems in which key mechanisms regulating cardiogenesis have been defined in great detail at the genetic and molecular level. Due to evolutionary conservation, these findings have also provided major inputs into studies of cardiogenesis in vertebrates. Many of the major components that control Drosophila cardiogenesis were discovered based on candidate gene approaches and their functions were defined by employing the outstanding genetic tools and molecular techniques available in this system. More recently, approaches have been taken that aim to interrogate the entire genome in order to identify novel components and describe genomic features that are pertinent to the regulation of heart development. Apart from classical forward genetic screens, the availability of the thoroughly annotated Drosophila genome sequence made new genome-wide approaches possible, which include the generation of massive numbers of RNA interference (RNAi) reagents that were used in forward genetic screens, as well as studies of the transcriptomes and proteomes of the developing heart under normal and experimentally manipulated conditions. Moreover, genome-wide chromatin immunoprecipitation experiments have been performed with the aim to define the full set of genomic binding sites of the major cardiogenic transcription factors, their relevant target genes, and a more complete picture of the regulatory network that drives cardiogenesis. This review will give an overview on these genome-wide approaches to Drosophila heart development and on computational analyses of the obtained information that ultimately aim to provide a description of this process at the systems level.

E2F function in muscle growth is necessary and sufficient for viability in Drosophila

The E2F transcription factor is a key cell cycle regulator. However, the inactivation of the entire E2F family in Drosophila is permissive throughout most of animal development until pupation when lethality occurs. Here we show that E2F function in the adult skeletal muscle is essential for animal viability since providing E2F function in muscles rescues the lethality of the whole-body E2F-deficient animals. Muscle-specific loss of E2F results in a significant reduction in muscle mass and thinner myofibrils. We demonstrate that E2F is dispensable for proliferation of muscle progenitor cells, but is required during late myogenesis to directly control the expression of a set of muscle-specific genes. Interestingly, E2f1 provides a major contribution to the regulation of myogenic function, while E2f2 appears to be less important. These findings identify a key function of E2F in skeletal muscle required for animal viability, and illustrate how the cell cycle regulator is repurposed in post-mitotic cells.

Hox control of Drosophila larval anatomy; The Alary and Thoracic Alary-Related Muscles

The body plan of arthropods and vertebrates involves the formation of repetitive segments, which subsequently diversify to give rise to different body parts along the antero-posterior/rostro-caudal body axis. Anatomical variations between body segments are crucial for organ function and organismal fitness. Pioneering work in Drosophila has established that Hox transcription factors play key roles both in endowing initially identical segments with distinct identities and organogenesis. The focus of this review is on Alary Muscles (AMs) and the newly discovered Thoracic Alary-Related Muscles (TARMs). AMs and TARMs are thin muscles which together connect the circulatory system and different midgut regions to the exoskeleton, while intertwining with the respiratory tubular network. They were hypothesized to represent a new type of muscles with spring-like properties, maintaining internal organs in proper anatomical positions during larval locomotion. Both the morphology of TARMs relative to AMs, and morphogenesis of connected tissues is under Hox control, emphasizing the key role of Hox proteins in coordinating the anatomical development of the larva.

Mesodermal gene expression during the embryonic and larval development of the articulate brachiopod Terebratalia transversa

Background

Brachiopods undergo radial cleavage, which is distinct from the stereotyped development of closely related spiralian taxa. The mesoderm has been inferred to derive from the archenteron walls following gastrulation, and the primary mesoderm derivative in the larva is a complex musculature. To investigate the specification and differentiation of the mesoderm in the articulate brachiopod Terebratalia transversa, we have identified orthologs of genes involved in mesoderm development in other taxa and investigated their spatial and temporal expression during the embryonic and larval development of T. transversa.

Results

Orthologs of 17 developmental regulatory genes with roles in the development of the mesoderm in other bilaterian animals were found to be expressed in the developing mesoderm of T. transversa. Five genes, Tt.twist, Tt.GATA456, Tt.dachshund, Tt.mPrx, and Tt.NK1, were found to have expression throughout the archenteron wall at the radial gastrula stage, shortly after the initiation of gastrulation. Three additional genes, Tt.Pax1/9, Tt.MyoD, and Tt.Six1/2, showed expression at this stage in only a portion of the archenteron wall. Tt.eya, Tt.FoxC, Tt.FoxF, Tt.Mox, Tt.paraxis, Tt.Limpet, and Tt.Mef2 all showed initial mesodermal expression during later gastrula or early larval stages. At the late larval stage, Tt.dachshund, Tt.Limpet, and Tt.Mef2 showed expression in nearly all mesoderm cells, while all other genes were localized to specific regions of the mesoderm. Tt.FoxD and Tt.noggin both showed expression in the ventral mesoderm at the larval stages, with gastrula expression patterns in the archenteron roof and blastopore lip, respectively.

Conclusions

Expression analyses support conserved roles for developmental regulators in the specification and differentiation of the mesoderm during the development of T. transversa. Expression of multiple mesodermal factors in the archenteron wall during gastrulation supports previous morphological observations that this region gives rise to larval mesoderm. Localized expression domains during gastrulation and larval development evidence early regionalization of the mesoderm and provide a basis for hypotheses regarding the molecular regulation underlying the complex system of musculature observed in the larva.

Electronic supplementary material

The online version of this article (doi:10.1186/s13227-015-0004-8) contains supplementary material, which is available to authorized users.

Distinct genetic programs guide Drosophila circular and longitudinal visceral myoblast fusion

BACKGROUND

The visceral musculature of Drosophila larvae comprises circular visceral muscles tightly interwoven with longitudinal visceral muscles. During myogenesis, the circular muscles arise by one-to-one fusion of a circular visceral founder cell (FC) with a visceral fusion-competent myoblast (FCM) from the trunk visceral mesoderm, and longitudinal muscles arise from FCs of the caudal visceral mesoderm. Longitudinal FCs migrate anteriorly under guidance of fibroblast growth factors during embryogenesis; it is proposed that they fuse with FCMs from the trunk visceral mesoderm to give rise to syncytia containing up to six nuclei.

RESULTS

Using fluorescence in situ hybridization and immunochemical analyses, we investigated whether these fusion events during migration use the same molecular repertoire and cellular components as fusion-restricted myogenic adhesive structure (FuRMAS), the adhesive signaling center that mediates myoblast fusion in the somatic mesoderm. Longitudinal muscles were formed by the fusion of one FC with Sns-positive FCMs, and defects in FCM specification led to defects in longitudinal muscle formation. At the fusion sites, Duf/Kirre and the adaptor protein Rols7 accumulated in longitudinal FCs, and Blow and F-actin accumulated in FCMs. The accumulation of these four proteins at the fusion sites argues for FuRMAS-like adhesion and signaling centers. Longitudinal fusion was disturbed in rols and blow single, and scar wip double mutants. Mutants of wasp or its interaction partner wip had no defects in longitudinal fusion.

CONCLUSIONS

Our results indicated that all embryonic fusion events depend on the same cell-adhesion molecules, but that the need for Rols7 and regulators of F-actin distinctly differs. Rols7 was required for longitudinal visceral and somatic myoblast fusion but not for circular visceral fusion. Importantly, longitudinal fusion depended on Kette and SCAR/Wave but was independent of WASp-dependent Arp2/3 activation. Thus, the complexity of the players involved in muscle formation increases from binucleated circular muscles to longitudinal visceral muscles to somatic muscles.

Subtle changes in motif positioning cause tissue-specific effects on robustness of an enhancer's activity

Deciphering the specific contribution of individual motifs within cis-regulatory modules (CRMs) is crucial to understanding how gene expression is regulated and how this process is affected by sequence variation. But despite vast improvements in the ability to identify where transcription factors (TFs) bind throughout the genome, we are limited in our ability to relate information on motif occupancy to function from sequence alone. Here, we engineered 63 synthetic CRMs to systematically assess the relationship between variation in the content and spacing of motifs within CRMs to CRM activity during development using Drosophila transgenic embryos. In over half the cases, very simple elements containing only one or two types of TF binding motifs were capable of driving specific spatio-temporal patterns during development. Different motif organizations provide different degrees of robustness to enhancer activity, ranging from binary on-off responses to more subtle effects including embryo-to-embryo and within-embryo variation. By quantifying the effects of subtle changes in motif organization, we were able to model biophysical rules that explain CRM behavior and may contribute to the spatial positioning of CRM activity in vivo. For the same enhancer, the effects of small differences in motif positions varied in developmentally related tissues, suggesting that gene expression may be more susceptible to sequence variation in one tissue compared to another. This result has important implications for human eQTL studies in which many associated mutations are found in cis-regulatory regions, though the mechanism for how they affect tissue-specific gene expression is often not understood.

Jeb/Alk signalling regulates the Lame duck GLI family transcription factor in the Drosophila visceral mesoderm

The Jelly belly (Jeb)/Anaplastic Lymphoma Kinase (Alk) signalling pathway regulates myoblast fusion in the circular visceral mesoderm (VM) of Drosophila embryos via specification of founder cells. However, only a limited number of target molecules for this pathway are described. We have investigated the role of the Lame Duck (Lmd) transcription factor in VM development in relationship to Jeb/Alk signal transduction. We show that Alk signalling negatively regulates Lmd activity post-transcriptionally through the MEK/MAPK (ERK) cascade resulting in a relocalisation of Lmd protein from the nucleus to cytoplasm. It has previously been shown that downregulation of Lmd protein is necessary for the correct specification of founder cells. In the visceral mesoderm of lmd mutant embryos, fusion-competent myoblasts seem to be converted to 'founder-like' cells that are still able to build a gut musculature even in the absence of fusion. The ability of Alk signalling to downregulate Lmd protein requires the N-terminal 140 amino acids, as a Lmd(141-866) mutant remains nuclear in the presence of active ALK and is able to drive robust expression of the Lmd downstream target Vrp1 in the developing VM. Our results suggest that Lmd is a target of Jeb/Alk signalling in the VM of Drosophila embryos.

Org-1 is required for the diversification of circular visceral muscle founder cells and normal midgut morphogenesis

The T-Box family of transcription factors plays fundamental roles in the generation of appropriate spatial and temporal gene expression profiles during cellular differentiation and organogenesis in animals. In this study we report that the Drosophila Tbx1 orthologue optomotor-blind-related-gene-1 (org-1) exerts a pivotal function in the diversification of circular visceral muscle founder cell identities in Drosophila. In embryos mutant for org-1, the specification of the midgut musculature per se is not affected, but the differentiating midgut fails to form the anterior and central midgut constrictions and lacks the gastric caeca. We demonstrate that this phenotype results from the nearly complete loss of the founder cell specific expression domains of several genes known to regulate midgut morphogenesis, including odd-paired (opa), teashirt (tsh), Ultrabithorax (Ubx), decapentaplegic (dpp) and wingless (wg). To address the mechanisms that mediate the regulatory inputs from org-1 towards Ubx, dpp, and wg in these founder cells we genetically dissected known visceral mesoderm specific cis-regulatory-modules (CRMs) of these genes. The analyses revealed that the activities of the dpp and wg CRMs depend on org-1, the CRMs are bound by Org-1 in vivo and their T-Box binding sites are essential for their activation in the visceral muscle founder cells. We conclude that Org-1 acts within a well-defined signaling and transcriptional network of the trunk visceral mesoderm as a crucial founder cell-specific competence factor, in concert with the general visceral mesodermal factor Biniou. As such, it directly regulates several key genes involved in the establishment of morphogenetic centers along the anteroposterior axis of the visceral mesoderm, which subsequently organize the formation of midgut constrictions and gastric caeca and thereby determine the morphology of the midgut.

High-fidelity promoter profiling reveals widespread alternative promoter usage and transposon-driven developmental gene expression

Many eukaryotic genes possess multiple alternative promoters with distinct expression specificities. Therefore, comprehensively annotating promoters and deciphering their individual regulatory dynamics is critical for gene expression profiling applications and for our understanding of regulatory complexity. We introduce RAMPAGE, a novel promoter activity profiling approach that combines extremely specific 5'-complete cDNA sequencing with an integrated data analysis workflow, to address the limitations of current techniques. RAMPAGE features a streamlined protocol for fast and easy generation of highly multiplexed sequencing libraries, offers very high transcription start site specificity, generates accurate and reproducible promoter expression measurements, and yields extensive transcript connectivity information through paired-end cDNA sequencing. We used RAMPAGE in a genome-wide study of promoter activity throughout 36 stages of the life cycle of Drosophila melanogaster, and describe here a comprehensive data set that represents the first available developmental time-course of promoter usage. We found that >40% of developmentally expressed genes have at least two promoters and that alternative promoters generally implement distinct regulatory programs. Transposable elements, long proposed to play a central role in the evolution of their host genomes through their ability to regulate gene expression, contribute at least 1300 promoters shaping the developmental transcriptome of D. melanogaster. Hundreds of these promoters drive the expression of annotated genes, and transposons often impart their own expression specificity upon the genes they regulate. These observations provide support for the theory that transposons may drive regulatory innovation through the distribution of stereotyped cis-regulatory modules throughout their host genomes.

Predicting spatial and temporal gene expression using an integrative model of transcription factor occupancy and chromatin state

Precise patterns of spatial and temporal gene expression are central to metazoan complexity and act as a driving force for embryonic development. While there has been substantial progress in dissecting and predicting cis-regulatory activity, our understanding of how information from multiple enhancer elements converge to regulate a gene's expression remains elusive. This is in large part due to the number of different biological processes involved in mediating regulation as well as limited availability of experimental measurements for many of them. Here, we used a Bayesian approach to model diverse experimental regulatory data, leading to accurate predictions of both spatial and temporal aspects of gene expression. We integrated whole-embryo information on transcription factor recruitment to multiple cis-regulatory modules, insulator binding and histone modification status in the vicinity of individual gene loci, at a genome-wide scale during Drosophila development. The model uses Bayesian networks to represent the relation between transcription factor occupancy and enhancer activity in specific tissues and stages. All parameters are optimized in an Expectation Maximization procedure providing a model capable of predicting tissue- and stage-specific activity of new, previously unassayed genes. Performing the optimization with subsets of input data demonstrated that neither enhancer occupancy nor chromatin state alone can explain all gene expression patterns, but taken together allow for accurate predictions of spatio-temporal activity. Model predictions were validated using the expression patterns of more than 600 genes recently made available by the BDGP consortium, demonstrating an average 15-fold enrichment of genes expressed in the predicted tissue over a naïve model. We further validated the model by experimentally testing the expression of 20 predicted target genes of unknown expression, resulting in an accuracy of 95% for temporal predictions and 50% for spatial. While this is, to our knowledge, the first genome-wide approach to predict tissue-specific gene expression in metazoan development, our results suggest that integrative models of this type will become more prevalent in the future.

Development is a complex process in which a single cell gives rise to a multi-cellular organism comprised of diverse cell types and well-organized tissues. This transformation requires tightly coordinated expression, both spatially and temporally, of hundreds to thousands of genes specific to any given tissue. To orchestrate these patterns, gene expression is regulated at multiple steps, from TF binding to cis-regulatory modules, general transcription factor and RNA polymerase II recruitment to promoters, chromatin remodeling, and three-dimensional looping interactions. Despite this level of complexity, the regulation of gene expression is typically modeled in the context of transcription factor binding and a single enhancer's activity as this is where the majority of experimental data is available. Recent advances in the measurement of chromatin modifications and insulator binding during embryogenesis provide new datasets that can be used for modeling gene expression. Here we use a Bayesian approach to integrate all three levels of information to combine the activity of multiple regulatory elements into a single model of a gene's expression, implementing an expectation maximization strategy to overcome the problem of missing data. Importantly, while the data for histone modifications and insulator binding represents merged signals from all cells in the embryo, the model can extract cell type specific and stage-specific predictions on gene expression for hundreds of genes of unknown expression.

A graphical modelling approach to the dissection of highly correlated transcription factor binding site profiles

Inferring the combinatorial regulatory code of transcription factors (TFs) from genome-wide TF binding profiles is challenging. A major reason is that TF binding profiles significantly overlap and are therefore highly correlated. Clustered occurrence of multiple TFs at genomic sites may arise from chromatin accessibility and local cooperation between TFs, or binding sites may simply appear clustered if the profiles are generated from diverse cell populations. Overlaps in TF binding profiles may also result from measurements taken at closely related time intervals. It is thus of great interest to distinguish TFs that directly regulate gene expression from those that are indirectly associated with gene expression. Graphical models, in particular Bayesian networks, provide a powerful mathematical framework to infer different types of dependencies. However, existing methods do not perform well when the features (here: TF binding profiles) are highly correlated, when their association with the biological outcome is weak, and when the sample size is small. Here, we develop a novel computational method, the Neighbourhood Consistent PC (NCPC) algorithms, which deal with these scenarios much more effectively than existing methods do. We further present a novel graphical representation, the Direct Dependence Graph (DDGraph), to better display the complex interactions among variables. NCPC and DDGraph can also be applied to other problems involving highly correlated biological features. Both methods are implemented in the R package ddgraph, available as part of Bioconductor (http://bioconductor.org/packages/2.11/bioc/html/ddgraph.html). Applied to real data, our method identified TFs that specify different classes of cis-regulatory modules (CRMs) in Drosophila mesoderm differentiation. Our analysis also found depletion of the early transcription factor Twist binding at the CRMs regulating expression in visceral and somatic muscle cells at later stages, which suggests a CRM-specific repression mechanism that so far has not been characterised for this class of mesodermal CRMs.

Transcription factors (TFs) are proteins that bind to DNA and regulate gene expression. Recent technological advances make it possible to map TF binding patterns across the whole genome. Multiple single-gene studies showed that combinatorial binding of multiple transcription factors determines the gene transcriptional output. A common naive assumption is that correlated binding profiles may indicate combinatorial binding. However, it has been found that many TFs bind to distinct hotspots whose role is currently unclear. It is thus of great interest to find transcription factor combinations whose correlated binding is causally most immediate to gene expression. Building upon theories of statistical dependence and causality, we develop novel graphical modelbased algorithms that handle highly correlated transcription factor binding profiles more efficiently and reliably than existing algorithms do. These algorithms can also be applied to other biological areas involving highly correlated variables, such as the analysis of high-throughput gene knock-down experiments.