Regulation of the BMP Signaling-Responsive Transcriptional Network in the Drosophila Embryo

Genome-wide analysis of Smad and Schnurri transcription factors in C. elegans demonstrates widespread interaction and a function in collagen secretion

Smads and their transcription factor partners mediate the transcriptional responses of target cells to secreted ligands of the transforming growth factor-β (TGF-β) family, including those of the conserved bone morphogenetic protein (BMP) family, yet only a small number of direct target genes have been well characterized. In C. elegans, the BMP2/4 ortholog DBL-1 regulates multiple biological functions, including body size, via a canonical receptor-Smad signaling cascade. Here, we identify functional binding sites for SMA-3/Smad and its transcriptional partner SMA-9/Schnurri based on ChIP-seq peaks (identified by modEncode) and expression differences of nearby genes identified from RNA-seq analysis of corresponding mutants. We found that SMA-3 and SMA-9 have both overlapping and unique target genes. At a genome-wide scale, SMA-3/Smad acts as a transcriptional activator, whereas SMA-9/Schnurri direct targets include both activated and repressed genes. Mutations in sma-9 partially suppress the small body size phenotype of sma-3, suggesting some level of antagonism between these factors and challenging the prevailing model for Schnurri function. Functional analysis of target genes revealed a novel role in body size for genes involved in one-carbon metabolism and in the endoplasmic reticulum (ER) secretory pathway, including the disulfide reductase dpy-11. Our findings indicate that Smads and SMA-9/Schnurri have previously unappreciated complex genetic and genomic regulatory interactions that in turn regulate the secretion of extracellular components like collagen into the cuticle to mediate body size regulation.

Genome-wide analysis of Smad and Schnurri transcription factors in C. elegans demonstrates widespread interaction and a function in collagen secretion

Smads and their transcription factor partners mediate the transcriptional responses of target cells to secreted ligands of the Transforming Growth Factor-β (TGF-β) family, including those of the conserved bone morphogenetic protein (BMP) family, yet only a small number of direct target genes have been well characterized. In C. elegans, the BMP2/4 ortholog DBL-1 regulates multiple biological functions, including body size, via a canonical receptor-Smad signaling cascade. Here, we identify functional binding sites for SMA-3/Smad and its transcriptional partner SMA-9/Schnurri based on ChIP-seq peaks (identified by modEncode) and expression differences of nearby genes identified from RNA-seq analysis of corresponding mutants. We found that SMA-3 and SMA-9 have both overlapping and unique target genes. At a genome-wide scale, SMA-3/Smad acts as a transcriptional activator, whereas SMA-9/Schnurri direct targets include both activated and repressed genes. Mutations in sma-9 partially suppress the small body size phenotype of sma-3, suggesting some level of antagonism between these factors and challenging the prevailing model for Schnurri function. Functional analysis of target genes revealed a novel role in body size for genes involved in one-carbon metabolism and in the endoplasmic reticulum (ER) secretory pathway, including the disulfide reductase dpy-11. Our findings indicate that Smads and SMA-9/Schnurri have previously unappreciated complex genetic and genomic regulatory interactions that in turn regulate the secretion of extracellular components like collagen into the cuticle to mediate body size regulation.

Upregulated bone morphogenetic protein 8A (BMP8A) in triple negative breast cancer (TNBC) and its involvement in the bone metastasis

Objective

The present study aimed to investigate the involvement of aberrant BMP8A expression in TNBC and bone metastasis.

Methods

Aberrant expression of BMP8A in breast cancer was first determined by analyzing The Cancer Genome Atlas breast cancer cohort (TCGA-BRCA) and an immunohistochemical (IHC) staining of BMP8A in a breast cancer tissue microarray (TMA). Clinical relevance of deregulated BMP8A in breast cancer was assessed using Kaplan-Meier online analysis. The influence of BMP8A on cellular functions of two TNBC cell lines was assessed using in vitro assays. Conditional medium (CM) collected from the supernatant of hFOB cells and bone matrix extract (BME) was applied to mimic the bone micro-environment to evaluate the role played by BMP8A in bone metastasis. Correlations with both osteolytic and osteoblastic markers were evaluated in the TCGA-BRCA cohort. Expression of certain responsive genes was quantified in the BMP8A overexpression cell lines. Additionally, signal transduction through both Smad-dependent and independent pathways was evaluated using Western blot assay.

Results

Compared to the adjacent normal tissues, BMP8A expression was significantly increased in primary tumors (p < 0.05) which was associated with shorter distant metastasis free survival (DMFS) in TNBC (p < 0.05). BMP8A was observed to enhance cell invasion and migration within TNBC cells. In the simulated bone milieu, both MDA-MB-231BMP8Aexp and BT549BMP8Aexp cells presented enhanced invasiveness. BMP8A level was strongly correlated with most osteolytic and osteoblastic markers, suggesting the potential involvement of BMP8A in bone metastasis in TNBC. Receptor activator of nuclear factor kappa-B ligand (RANKL) expression was significantly increased in BMP8A overexpressed triple-negative cell lines (MDA-MB-231 and BT549). Furthermore, enhanced phosphorylation of Smad3 and increased expression of epidermal growth factor receptor (EGFR) were observed in MDA-MB-231 cells overexpressing BMP8A.

Conclusion

BMP8A was upregulated in TNBC which was associated with poorer DMFS. BMP8A overexpression enhanced the invasion and migration of TNBC cells. With a putative role in osteolytic bone metastasis in TNBC, BMP8A represents a promising candidate for further investigation into its therapeutic potential.

Dietary sucrose determines the regulatory activity of lithium on gene expression and lifespan in Drosophila melanogaster

The amount of dietary sugars and the administration of lithium both impact the lifespan of the fruit fly Drosophila melanogaster. It is noteworthy that lithium is attributed with insulin-like activity as it stimulates protein kinase B/Akt and suppresses the activity of glycogen synthase kinase-3 (GSK-3). However, its interaction with dietary sugar has largely remained unexplored. Therefore, we investigated the effects of lithium supplementation on known lithium-sensitive parameters in fruit flies, such as lifespan, body composition, GSK-3 phosphorylation, and the transcriptome, while varying the dietary sugar concentration. For all these parameters, we observed that the efficacy of lithium was significantly influenced by the sucrose content in the diet. Overall, we found that lithium was most effective in enhancing longevity and altering body composition when added to a low-sucrose diet. Whole-body RNA sequencing revealed a remarkably similar transcriptional response when either increasing dietary sucrose from 1% to 10% or adding 1 mM LiCl to a 1% sucrose diet, characterized by a substantial overlap of nearly 500 differentially expressed genes. Hence, dietary sugar supply is suggested as a key factor in understanding lithium bioactivity, which could hold relevance for its therapeutic applications.

Bone morphogenetic protein signaling: the pathway and its regulation

In the mid-1960s, bone morphogenetic proteins (BMPs) were first identified in the extracts of bone to have the remarkable ability to induce heterotopic bone. When the Drosophila gene decapentaplegic (dpp) was first identified to share sequence similarity with mammalian BMP2/BMP4 in the late-1980s, it became clear that secreted BMP ligands can mediate processes other than bone formation. Following this discovery, collaborative efforts between Drosophila geneticists and mammalian biochemists made use of the strengths of their respective model systems to identify BMP signaling components and delineate the pathway. The ability to conduct genetic modifier screens in Drosophila with relative ease was critical in identifying the intracellular signal transducers for BMP signaling and the related transforming growth factor-beta/activin signaling pathway. Such screens also revealed a host of genes that encode other core signaling components and regulators of the pathway. In this review, we provide a historical account of this exciting time of gene discovery and discuss how the field has advanced over the past 30 years. We have learned that while the core BMP pathway is quite simple, composed of 3 components (ligand, receptor, and signal transducer), behind the versatility of this pathway lies multiple layers of regulation that ensures precise tissue-specific signaling output. We provide a sampling of these discoveries and highlight many questions that remain to be answered to fully understand the complexity of BMP signaling.

Analysis of SMAD1/5 target genes in a sea anemone reveals ZSWIM4-6 as a novel BMP signaling modulator

BMP signaling has a conserved function in patterning the dorsal-ventral body axis in Bilateria and the directive axis in anthozoan cnidarians. So far, cnidarian studies have focused on the role of different BMP signaling network components in regulating pSMAD1/5 gradient formation. Much less is known about the target genes downstream of BMP signaling. To address this, we generated a genome-wide list of direct pSMAD1/5 target genes in the anthozoan Nematostella vectensis, several of which were conserved in Drosophila and Xenopus. Our ChIP-seq analysis revealed that many of the regulatory molecules with documented bilaterally symmetric expression in Nematostella are directly controlled by BMP signaling. We identified several so far uncharacterized BMP-dependent transcription factors and signaling molecules, whose bilaterally symmetric expression may be indicative of their involvement in secondary axis patterning. One of these molecules is zswim4-6, which encodes a novel nuclear protein that can modulate the pSMAD1/5 gradient and potentially promote BMP-dependent gene repression.

Neofunctionalization of Toll Signaling in Insects: From Immunity to Dorsoventral Patterning

Toll signaling plays a crucial role in pathogen defense throughout the animal kingdom. It was discovered, however, for its function in dorsoventral (DV) axis formation in Drosophila. In all other insects studied so far, but not outside the insects, Toll is also required for DV patterning. However, in insects more distantly related to Drosophila, Toll's patterning role is frequently reduced and substituted by an expanded influence of BMP signaling, the pathway implicated in DV axis formation in all major metazoan lineages. This suggests that Toll was integrated into an ancestral BMP-based patterning system at the base of the insects or during insect evolution. The observation that Toll signaling has an immune function in the extraembryonic serosa, an early differentiating tissue of most insect embryos, suggests a scenario of how Toll was co-opted from an ancestral immune function for its new role in axis formation.

Fragile X mental retardation protein coordinates neuron-to-glia communication for clearance of developmentally transient brain neurons

In the developmental remodeling of brain circuits, neurons are removed by glial phagocytosis to optimize adult behavior. Fragile X mental retardation protein (FMRP) regulates neuron-to-glia signaling to drive glial phagocytosis for targeted neuron pruning. We find that FMRP acts in a mothers against decapentaplegic (Mad)-insulin receptor (InR)-protein kinase B (Akt) pathway to regulate pretaporter (Prtp) and amyloid precursor protein-like (APPL) signals directing this glial clearance. Neuronal RNAi of Drosophila fragile X mental retardation 1 (dfmr1) elevates mad transcript levels and increases pMad signaling. Neuronal dfmr1 and mad RNAi both elevate phospho-protein kinase B (pAkt) and delay neuron removal but cause opposite effects on InR expression. Genetically correcting pAkt levels in the mad RNAi background restores normal remodeling. Consistently, neuronal dfmr1 and mad RNAi both decrease Prtp levels, whereas neuronal InR and akt RNAi increase Prtp levels, indicating FMRP works with pMad and insulin signaling to tightly regulate Prtp signaling and thus control glial phagocytosis for correct circuit remodeling. Neuronal dfmr1 and mad and akt RNAi all decrease APPL levels, with the pathway signaling higher glial endolysosome activity for phagocytosis. These findings reveal a FMRP-dependent control pathway for neuron-to-glia communication in neuronal pruning, identifying potential molecular mechanisms for devising fragile X syndrome treatments.

Combined modelling of mRNA decay dynamics and single-molecule imaging in the Drosophila embryo uncovers a role for P-bodies in 5' to 3' degradation

Regulation of mRNA degradation is critical for a diverse array of cellular processes and developmental cell fate decisions. Many methods for determining mRNA half-lives rely on transcriptional inhibition or metabolic labelling. Here, we use a non-invasive method for estimating half-lives for hundreds of mRNAs in the early Drosophila embryo. This approach uses the intronic and exonic reads from a total RNA-seq time series and Gaussian process regression to model the dynamics of premature and mature mRNAs. We show how regulation of mRNA stability is used to establish a range of mature mRNA dynamics during embryogenesis, despite shared transcription profiles. Using single-molecule imaging, we provide evidence that, for the mRNAs tested, there is a correlation between short half-life and mRNA association with P-bodies. Moreover, we detect an enrichment of mRNA 3' ends in P-bodies in the early embryo, consistent with 5' to 3' degradation occurring in P-bodies for at least a subset of mRNAs. We discuss our findings in relation to recently published data suggesting that the primary function of P-bodies in other biological contexts is mRNA storage.

Germline/soma distinction in Drosophila embryos requires regulators of zygotic genome activation

In Drosophila melanogaster embryos, somatic versus germline identity is the first cell fate decision. Zygotic genome activation (ZGA) orchestrates regionalized gene expression, imparting specific identity on somatic cells. ZGA begins with a minor wave that commences at nuclear cycle (NC)8 under the guidance of chromatin accessibility factors (Zelda, CLAMP, GAF), followed by the major wave during NC14. By contrast, primordial germ cell (PGC) specification requires maternally deposited and posteriorly anchored germline determinants. This is accomplished by a centrosome coordinated release and sequestration of germ plasm during the precocious cellularization of PGCs in NC10. Here, we report a novel requirement for Zelda and CLAMP during the establishment of the germline/soma distinction. When their activity is compromised, PGC determinants are not properly sequestered, and specification is disrupted. Conversely, the spreading of PGC determinants from the posterior pole adversely influences transcription in the neighboring somatic nuclei. These reciprocal aberrations can be correlated with defects in centrosome duplication/separation that are known to induce inappropriate transmission of the germ plasm. Interestingly, consistent with the ability of bone morphogenetic protein (BMP) signaling to influence specification of embryonic PGCs, reduction in the transcript levels of a BMP family ligand, decapentaplegic (dpp), is exacerbated at the posterior pole.

How two extraembryonic epithelia became one: serosa and amnion features and functions of Drosophila's amnioserosa

The conservation of gene networks that specify and differentiate distinct tissues has long been a subject of great interest to evolutionary developmental biologists, but the question of how pre-existing tissue-specific developmental trajectories merge is rarely asked. During the radiation of flies, two extraembryonic epithelia, known as serosa and amnion, evolved into one, called amnioserosa. This unique extraembryonic epithelium is found in fly species of the group Schizophora, including the genetic model organism Drosophila melanogaster, and has been studied in depth. Close relatives of this group develop a serosa and a rudimentary amnion. The scuttle fly Megaselia abdita has emerged as an excellent model organism to study this extraembryonic tissue organization. In this review, development and functions of the extraembryonic tissue complements of Drosophila and Megaselia are compared. It is concluded that the amnioserosa combines cells, genetic pathway components and functions that were previously associated either with serosa development or amnion development. The composite developmental trajectory of the amnioserosa raises the question of whether merging tissue-specific gene networks is a common evolutionary process. This article is part of the theme issue 'Extraembryonic tissues: exploring concepts, definitions and functions across the animal kingdom'.

Chromatin accessibility-based characterisation of brain gene regulatory networks in three distinct honey bee polyphenisms

The honey bee genome has the capacity to produce three phenotypically distinct organisms (two diploid female castes: queen and worker, and a haploid male drone). Previous studies have implicated metabolic flux acting via epigenetic regulation in directing nutrition-driven phenotypic plasticity in the honey bee. However, the cis-acting DNA regulatory elements that establish tissue and polyphenism -specific epigenomes and gene expression programmes, remain unclear. Using a high resolution multiomic approach including assay for transposase-accessible chromatin by sequencing (ATAC-seq), RNA-seq and ChIP-seq, we produce the first genome-wide maps of the regulatory landscape across all three adult honey bee phenotypes identifying > 5000 regulatory regions in queen, 7500 in worker and 6500 in drone, with the vast majority of these sites located within intronic regions. These regions are defined by positive enrichment of H3K27ac and depletion of H3K4me3 and show a positive correlation with gene expression. Using ATAC-seq footprinting we determine queen, worker and drone -specific transcription factor occupancy and uncover novel phenotype-specific regulatory networks identifying two key nuclear receptors that have previously been implicated in caste-determination and adult behavioural maturation in honey bees; ecdysone receptor and ultraspiracle. Collectively, this study provides novel insights into key gene regulatory networks that are associated with these distinct polyphenisms in the honey bee.

Dichotomous cis-regulatory motifs mediate the maturation of the neuromuscular junction by retrograde BMP signaling

Retrograde bone morphogenetic protein (BMP) signaling at the Drosophila neuromuscular junction (NMJ) has served as a paradigm to study TGF-β-dependent synaptic function and maturation. Yet, how retrograde BMP signaling transcriptionally regulates these functions remains unresolved. Here, we uncover a gene network, enriched for neurotransmission-related genes, that is controlled by retrograde BMP signaling in motor neurons through two Smad-binding cis-regulatory motifs, the BMP-activating (BMP-AE) and silencer (BMP-SE) elements. Unpredictably, both motifs mediate direct gene activation, with no involvement of the BMP derepression pathway regulators Schnurri and Brinker. Genome editing of candidate BMP-SE and BMP-AE within the locus of the active zone gene bruchpilot, and a novel Ly6 gene witty, demonstrated the role of these motifs in upregulating genes required for the maturation of pre- and post-synaptic NMJ compartments. Our findings uncover how Smad-dependent transcriptional mechanisms specific to motor neurons directly orchestrate a gene network required for synaptic maturation by retrograde BMP signaling.

Evolution of the dorsoventral axis in insects: the changing role of Bone Morphogenetic Proteins

Embryonic dorsal-ventral (DV) patterning by Bone Morphogenetic Proteins (BMPs) is a conserved feature of Bilateria, based on graded BMP activity set up by diffusible BMP ligands and Chordin/Sog antagonists. In the fly Drosophila melanogaster BMP function is secondary to patterning by the Toll pathway, suggesting a more restricted role for BMPs in insects. With widespread genome sequencing technologies allied to functional analysis in a growing number of species, recent work has shown that BMP's role in DV patterning relative to Toll varies among insect orders. Further, the role of BMP antagonists to set up BMP gradients is also greatly diversified. Here we review the recent findings concerning the role of BMP in the DV patterning of insects and address the potential aspects that may have co-evolved with BMPs to attain this functional divergence.

Scalable inference of transcriptional kinetic parameters from MS2 time series data

MOTIVATION

The MS2-MCP (MS2 coat protein) live imaging system allows for visualization of transcription dynamics through the introduction of hairpin stem-loop sequences into a gene. A fluorescent signal at the site of nascent transcription in the nucleus quantifies mRNA production. Computational modelling can be used to infer the promoter states along with the kinetic parameters governing transcription, such as promoter switching frequency and polymerase loading rate. However, modelling of the fluorescent trace presents a challenge due its persistence; the observed fluorescence at a given time point depends on both current and previous promoter states. A compound state Hidden Markov Model (cpHMM) was recently introduced to allow inference of promoter activity from MS2-MCP data. However, the computational time for inference scales exponentially with gene length and the cpHMM is therefore not currently practical for application to many eukaryotic genes.

RESULTS

We present a scalable implementation of the cpHMM for fast inference of promoter activity and transcriptional kinetic parameters. This new method can model genes of arbitrary length through the use of a time-adaptive truncated compound state space. The truncated state space provides a good approximation to the full state space by retaining the most likely set of states at each time during the forward pass of the algorithm. Testing on MS2-MCP fluorescent data collected from early Drosophila melanogaster embryos indicates that the method provides accurate inference of kinetic parameters within a computationally feasible timeframe. The inferred promoter traces generated by the model can also be used to infer single-cell transcriptional parameters.

AVAILABILITY AND IMPLEMENTATION

Python implementation is available at https://github.com/ManchesterBioinference/burstInfer, along with code to reproduce the examples presented here.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

BMP Signaling: Lighting up the Way for Embryonic Dorsoventral Patterning

One of the most significant events during early embryonic development is the establishment of a basic embryonic body plan, which is defined by anteroposterior, dorsoventral (DV), and left-right axes. It is well-known that the morphogen gradient created by BMP signaling activity is crucial for DV axis patterning across a diverse set of vertebrates. The regulation of BMP signaling during DV patterning has been strongly conserved across evolution. This is a remarkable regulatory and evolutionary feat, as the BMP gradient has been maintained despite the tremendous variation in embryonic size and shape across species. Interestingly, the embryonic DV axis exhibits robust stability, even in face of variations in BMP signaling. Multiple lines of genetic, molecular, and embryological evidence have suggested that numerous BMP signaling components and their attendant regulators act in concert to shape the developing DV axis. In this review, we summarize the current knowledge of the function and regulation of BMP signaling in DV patterning. Throughout, we focus specifically on popular model animals, such as Xenopus and zebrafish, highlighting the similarities and differences of the regulatory networks between species. We also review recent advances regarding the molecular nature of DV patterning, including the initiation of the DV axis, the formation of the BMP gradient, and the regulatory molecular mechanisms behind BMP signaling during the establishment of the DV axis. Collectively, this review will help clarify our current understanding of the molecular nature of DV axis formation.

Preformation and epigenesis converge to specify primordial germ cell fate in the early Drosophila embryo

A critical step in animal development is the specification of primordial germ cells (PGCs), the precursors of the germline. Two seemingly mutually exclusive mechanisms are implemented across the animal kingdom: epigenesis and preformation. In epigenesis, PGC specification is non-autonomous and depends on extrinsic signaling pathways. The BMP pathway provides the key PGC specification signals in mammals. Preformation is autonomous and mediated by determinants localized within PGCs. In Drosophila, a classic example of preformation, constituents of the germ plasm localized at the embryonic posterior are thought to be both necessary and sufficient for proper determination of PGCs. Contrary to this longstanding model, here we show that these localized determinants are insufficient by themselves to direct PGC specification in blastoderm stage embryos. Instead, we find that the BMP signaling pathway is required at multiple steps during the specification process and functions in conjunction with components of the germ plasm to orchestrate PGC fate.

Proper specification of primordial germ cells (PGCs) is crucial as PGCs serve as the precursors of germline stem cells. To specify PGC fate, invertebrates rely upon cell autonomous preformation involving maternally deposited germ plasm. In Drosophila melanogaster, to insulate newly formed PGCs from the adverse effects of the cell-cell signaling pathways, germ plasm determinants silence transcription and attenuate the cell cycle. However, our data on the BMP signaling pathway challenge this long-held view of PGC specification and suggest that appropriate specification of embryonic PGCs is sensitive to the BMP ligand, decapentaplegic (dpp), and its cognate receptor, thickveins. We find that PGCs are not only capable of responding to BMP signals from the soma, but also that these signals impact the proper determination of the germ cells. Based on these unanticipated similarities between mammals and flies, we propose a model integrating contribution of both the cell-autonomous (preformation) and non-autonomous (epigenesis) pathways during PGC determination. Consistent with the model, we have observed dominant genetic interactions between, oskar, the maternal determinant of PGC fate, and the BMP pathway ligand dpp.

Cis-regulatory chromatin loops arise before TADs and gene activation, and are independent of cell fate during early Drosophila development

Acquisition of cell fate is thought to rely on the specific interaction of remote cis-regulatory modules (CRMs), for example, enhancers and target promoters. However, the precise interplay between chromatin structure and gene expression is still unclear, particularly within multicellular developing organisms. In the present study, we employ Hi-M, a single-cell spatial genomics approach, to detect CRM-promoter looping interactions within topologically associating domains (TADs) during early Drosophila development. By comparing cis-regulatory loops in alternate cell types, we show that physical proximity does not necessarily instruct transcriptional states. Moreover, multi-way analyses reveal that multiple CRMs spatially coalesce to form hubs. Loops and CRM hubs are established early during development, before the emergence of TADs. Moreover, CRM hubs are formed, in part, via the action of the pioneer transcription factor Zelda and precede transcriptional activation. Our approach provides insight into the role of CRM-promoter interactions in defining transcriptional states, as well as distinct cell types.

The BMP signaling gradient is interpreted through concentration thresholds in dorsal-ventral axial patterning

Bone Morphogenetic Protein (BMP) patterns the dorsal–ventral (DV) embryonic axis in all vertebrates, but it is unknown how cells along the DV axis interpret and translate the gradient of BMP signaling into differential gene activation that will give rise to distinct cell fates. To determine the mechanism of BMP morphogen interpretation in the zebrafish gastrula, we identified 57 genes that are directly activated by BMP signaling. By using Seurat analysis of single-cell RNA sequencing (scRNA-seq) data, we found that these genes are expressed in at least 3 distinct DV domains of the embryo. We distinguished between 3 models of BMP signal interpretation in which cells activate distinct gene expression through interpretation of thresholds of (1) the BMP signaling gradient slope; (2) the BMP signal duration; or (3) the level of BMP signal activation. We tested these 3 models using quantitative measurements of phosphorylated Smad5 (pSmad5) and by examining the spatial relationship between BMP signaling and activation of different target genes at single-cell resolution across the embryo. We found that BMP signaling gradient slope or BMP exposure duration did not account for the differential target gene expression domains. Instead, we show that cells respond to 3 distinct levels of BMP signaling activity to activate and position target gene expression. Together, we demonstrate that distinct pSmad5 threshold levels activate spatially distinct target genes to pattern the DV axis.

This study tested three models of how a BMP morphogen gradient is translated into differential gene activation that specifies distinct cell fates, finding that BMP signal concentration thresholds, not gradient shape or signal duration, position three distinct gene activation domains.

Modulation of the Promoter Activation Rate Dictates the Transcriptional Response to Graded BMP Signaling Levels in the Drosophila Embryo

Morphogen gradients specify cell fates during development, with a classic example being the bone morphogenetic protein (BMP) gradient's conserved role in embryonic dorsal-ventral axis patterning. Here, we elucidate how the BMP gradient is interpreted in the Drosophila embryo by combining live imaging with computational modeling to infer transcriptional burst parameters at single-cell resolution. By comparing burst kinetics in cells receiving different levels of BMP signaling, we show that BMP signaling controls burst frequency by regulating the promoter activation rate. We provide evidence that the promoter activation rate is influenced by both enhancer and promoter sequences, whereas Pol II loading rate is primarily modulated by the enhancer. Consistent with BMP-dependent regulation of burst frequency, the numbers of BMP target gene transcripts per cell are graded across their expression domains. We suggest that graded mRNA output is a general feature of morphogen gradient interpretation and discuss how this can impact on cell-fate decisions.

Coacting enhancers can have complementary functions within gene regulatory networks and promote canalization

Developmental genes are often regulated by multiple enhancers exhibiting similar spatiotemporal outputs, which are generally considered redundantly acting though few have been studied functionally. Using CRISPR-Cas9, we created deletions of two enhancers, brk5' and brk3', that drive similar but not identical expression of the gene brinker (brk) in early Drosophila embryos. Utilizing both in situ hybridization and quantitative mRNA analysis, we investigated the changes in the gene network state caused by the removal of one or both of the early acting enhancers. brk5' deletion generally phenocopied the gene mutant, including expansion of the BMP ligand decapentaplegic (dpp) as well as inducing variability in amnioserosa tissue cell number suggesting a loss of canalization. In contrast, brk3' deletion presented unique phenotypes including dorsal expansion of several ventrally expressed genes and a decrease in amnioserosa cell number. Similarly, deletions were made for two enhancers associated with the gene short-gastrulation (sog), sog.int and sog.dist, demonstrating that they also exhibit distinct patterning phenotypes and affect canalization. In summary, this study shows that similar gene expression driven by coacting enhancers can support distinct, and sometimes complementary, functions within gene regulatory networks and, moreover, that phenotypes associated with individual enhancer deletion mutants can provide insight into new gene functions.

Retrograde BMP signaling activates neuronal gene expression through widespread deployment of a conserved BMP-responsive cis-regulatory activation element

Retrograde Bone Morphogenetic Protein (BMP) signaling in neurons is essential for the differentiation and synaptic function of many neuronal subtypes. BMP signaling regulates these processes via Smad transcription factor activity, yet the scope and nature of Smad-dependent gene regulation in neurons are mostly unknown. Here, we applied a computational approach to predict Smad-binding cis-regulatory BMP-Activating Elements (BMP-AEs) in Drosophila, followed by transgenic in vivo reporter analysis to test their neuronal subtype enhancer activity in the larval central nervous system (CNS). We identified 34 BMP-AE-containing genomic fragments that are responsive to BMP signaling in neurons, and showed that the embedded BMP-AEs are required for this activity. RNA-seq analysis identified BMP-responsive genes in the CNS and revealed that BMP-AEs selectively enrich near BMP-activated genes. These data suggest that functional BMP-AEs control nearby BMP-activated genes, which we validated experimentally. Finally, we demonstrated that the BMP-AE motif mediates a conserved Smad-responsive function in the Drosophila and vertebrate CNS. Our results provide evidence that BMP signaling controls neuronal function by directly coordinating the expression of a battery of genes through widespread deployment of a conserved Smad-responsive cis-regulatory motif.

A Developmental Program Truncates Long Transcripts to Temporally Regulate Cell Signaling

Rapid mitotic divisions and a fixed transcription rate limit the maximal length of transcripts in early Drosophila embryos. Previous studies suggested that transcription of long genes is initiated but aborted, as early nuclear divisions have short interphases. Here, we identify long genes that are expressed during short nuclear cycles as truncated transcripts. The RNA binding protein Sex-lethal physically associates with transcripts for these genes and is required to support early termination to specify shorter transcript isoforms in early embryos of both sexes. In addition, one truncated transcript for the gene short-gastrulation encodes a product in embryos that functionally relates to a previously characterized dominant-negative form, which maintains TGF-β signaling in the off-state. In summary, our results reveal a developmental program of short transcripts functioning to help temporally regulate Drosophila embryonic development, keeping cell signaling at early stages to a minimum in order to support its proper initiation at cellularization.

Phenotypically distinct female castes in honey bees are defined by alternative chromatin states during larval development

The capacity of the honey bee to produce three phenotypically distinct organisms (two female castes; queens and sterile workers, and haploid male drones) from one genotype represents one of the most remarkable examples of developmental plasticity in any phylum. The queen-worker morphological and reproductive divide is environmentally controlled during post-embryonic development by differential feeding. Previous studies implicated metabolic flux acting via epigenetic regulation, in particular DNA methylation and microRNAs, in establishing distinct patterns of gene expression underlying caste-specific developmental trajectories. We produce the first genome-wide maps of chromatin structure in the honey bee at a key larval stage in which developmental canalization into queen or worker is virtually irreversible. We find extensive genome-wide differences in H3K4me3, H3K27ac, and H3K36me3, many of which correlate with caste-specific transcription. Furthermore, we identify H3K27ac as a key chromatin modification, with caste-specific regions of intronic H3K27ac directing the worker caste. These regions may harbor the first examples of caste-specific enhancer elements in the honey bee. Our results demonstrate a key role for chromatin modifications in the establishment and maintenance of caste-specific transcriptional programs in the honey bee. We show that at 96 h of larval growth, the queen-specific chromatin pattern is already established, whereas the worker determination is not, thus providing experimental support for the perceived timing of this critical point in developmental heterochrony in two types of honey bee females. In a broader context, our study provides novel data on environmentally regulated organismal plasticity and the molecular foundation of the evolutionary origins of eusociality.

Juvenile hormone signaling in short germ-band hemimetabolan embryos

The role of juvenile hormone (JH) in insect embryos is far from understood, especially in short germ-band hemimetabolan species. To shed light on this issue, we depleted the mRNA levels of Krüppel homolog 1, Methoprene-tolerant and JH acid O-methyltransferase, key elements of JH signaling, in embryos of the short germ-band hemimetabolan species Blattella germanica This precluded the formation of the germ-band anlage in a group of embryos. Hatchability was also reduced, which might have been caused by premature upregulation of laccase 2, a promoter of cuticle tanning. In other cases, development was interrupted in mid embryogenesis, involving defects related to dorsal closure and appendage formation. These phenotypes possibly result from the low levels of Broad-complex (BR-C) produced under JH-depleted conditions. This contrasts with holometabolan species, in which JH does not promote BR-C expression, which remains low during embryo development. Possibly, the stimulatory role of JH on BR-C expression and the morphogenetic functions of BR-C in hemimetabolan embryos were lost in holometabolan species. If so, this might have been a key driver for the evolution of holometabolan metamorphosis.

FlyExpress 7: An Integrated Discovery Platform To Study Coexpressed Genes Using in Situ Hybridization Images in Drosophila

Gene expression patterns assayed across development can offer key clues about a gene's function and regulatory role. Drosophila melanogaster is ideal for such investigations as multiple individual and high-throughput efforts have captured the spatiotemporal patterns of thousands of embryonic expressed genes in the form of in situ images. FlyExpress (www.flyexpress.net), a knowledgebase based on a massive and unique digital library of standardized images and a simple search engine to find coexpressed genes, was created to facilitate the analytical and visual mining of these patterns. Here, we introduce the next generation of FlyExpress resources to facilitate the integrative analysis of sequence data and spatiotemporal patterns of expression from images. FlyExpress 7 now includes over 100,000 standardized in situ images and implements a more efficient, user-defined search algorithm to identify coexpressed genes via Genomewide Expression Maps (GEMs). Shared motifs found in the upstream 5' regions of any pair of coexpressed genes can be visualized in an interactive dotplot. Additional webtools and link-outs to assist in the downstream validation of candidate motifs are also provided. Together, FlyExpress 7 represents our largest effort yet to accelerate discovery via the development and dispersal of new webtools that allow researchers to perform data-driven analyses of coexpression (image) and genomic (sequence) data.

Modulation of BMP signalling by integrins

The bone morphogenetic protein (BMP) pathway is a major conserved signalling pathway with diverse roles in development and homeostasis. Given that cells exist in three-dimensional environments, one important area is to understand how the BMP pathway operates within such complex cellular environments. The extracellular matrix contains information regarding tissue architecture and its mechanical properties that is transmitted to the cell via integrin receptors. In this review, I describe various examples of modulation of the BMP pathway by integrins. In the case of the Drosophila embryo and some cell line-based studies, integrins have been found to enhance BMP responses through different mechanisms, such as enhancement of BMP ligand–receptor binding and effects on Smad phosphorylation or stability. In these contexts, BMP-dependent activation of integrins is a common theme. However, I also discuss examples where integrins inhibit the BMP pathway, highlighting the context-dependent nature of integrin–BMP cross-talk.

The Functionality and Evolution of Eukaryotic Transcriptional Enhancers

Enhancers regulate precise spatial and temporal patterns of gene expression in eukaryotes and, moreover, evolutionary changes in these modular cis-regulatory elements may represent the predominant genetic basis for phenotypic evolution. Here, we review approaches to identify and functionally analyze enhancers and their transcription factor binding sites, including assay for transposable-accessible chromatin-sequencing (ATAC-Seq) and clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9, respectively. We also explore enhancer functionality, including how transcription factor binding sites combine to regulate transcription, as well as research on shadow and super enhancers, and how enhancers can act over great distances and even in trans. Finally, we discuss recent theoretical and empirical data on how transcription factor binding sites and enhancers evolve. This includes how the function of enhancers is maintained despite the turnover of transcription factor binding sites as well as reviewing studies where mutations in enhancers have been shown to underlie morphological change.