Microarray comparison of anterior and posterior Drosophila wing imaginal disc cells identifies novel wing genes

Polycomb-mediated transgenerational epigenetic inheritance of Drosophila eye colour is independent of small RNAs

Transgenerational epigenetic inheritance (TEI) describes the process whereby distinct epigenetic states are transmitted between generations, resulting in heritable gene expression and phenotypic differences that are independent of DNA sequence variation. Chromatin modifications have been demonstrated to be important in TEI; however, the extent to which they require other signals to establish and maintain epigenetic states is still unclear. Here we investigate whether small non-coding RNAs contribute to different epigenetic states of the Fab2L transgene in Drosophila triggered by transient long-range chromosomal contacts, which requires Polycomb complex activity to deposit the H3K27me3 modification for long-term TEI. By analysing mutants deficient in small non-coding RNAs, high-throughput sequencing data, long-range chromosomal contacts and gene expression, we demonstrate that small non-coding RNAs do not contribute directly to initiation or maintenance of silencing. However, we uncover an indirect role for microRNA expression in transgene silencing through effects on the Polycomb group gene Pleiohomeotic. Additionally, we show that a commonly used marker gene, Stubble (Sb), affects Pleiohomeotic expression, which may be important in interpreting experiments assaying Polycomb function in Drosophila development. By ruling out a plausible candidate for TEI at the Fab2L transgene, our work highlights the variability in different modes of TEI across species.

Genome resequencing and genome-wide polymorphisms in mosquito vectors Aedes aegypti and Aedes albopictus from south India

Aedes aegypti and Aedes albopictus mosquitoes spread major vector-borne viral diseases in tropical and sub-tropical regions of the globe. In this study, we sequenced the genome of Indian Ae. aegypti and Ae. albopictus and mapped to their reference genomes. Comparative genomics were performed between our strain and the reference strains. A total of 14,416,484 single nucleotide polymorphisms (SNPs) and 156,487 insertions and deletions (InDels) were found in Ae. aegypti, and 28,940,433 SNPs and 188,987 InDels in Ae. albopictus. Particular emphasis was given to gene families involved in mosquito digestion, development, and innate immunity, which could be putative candidates for vector control. Serine protease cascades and their inhibitors called serpins, play a central role in these processes. We extracted high-impact variants in genes associated with serine proteases and serpins. This study reports for the first time a high coverage genome sequence data of an Indian Ae. albopictus mosquito. The results from this study will provide insights into Indian Aedes specific polymorphisms and the evolution of immune related genes in mosquitoes, which can serve as a resource for future comparative genomics and those pursuing the development of targeted biopesticides for effective mosquito control strategies.

The roles of growth regulation and appendage patterning genes in the morphogenesis of treehopper pronota

Treehoppers of the insect family Membracidae have evolved enlarged and elaborate pronotal structures, which is hypothesized to involve co-opted expression of genes that are shared with the wings. Here, we investigate the similarity between the pronotum and wings in relation to growth. Our study reveals that the ontogenetic allometry of the pronotum is similar to that of wings in Membracidae, but not the outgroup. Using transcriptomics, we identify genes related to translation and protein synthesis, which are mutually upregulated. These genes are implicated in the eIF2, eIF4/p70S6K and mTOR pathways, and have known roles in regulating cell growth and proliferation. We find that species-specific differential growth patterning of the pronotum begins as early as the third instar, which suggests that expression of appendage patterning genes occurs long before the metamorphic molt. We propose that a network related to growth and size determination is the more likely mechanism shared with wings. However, regulators upstream of the shared genes in pronotum and wings need to be elucidated to substantiate whether co-option has occurred. Finally, we believe it will be helpful to distinguish the mechanisms leading to pronotal size from those regulating pronotal shape as we make sense of this spectacular evolutionary innovation.

A cell atlas of adult muscle precursors uncovers early events in fibre-type divergence in Drosophila

In Drosophila, the wing disc‐associated muscle precursor cells give rise to the fibrillar indirect flight muscles (IFM) and the tubular direct flight muscles (DFM). To understand early transcriptional events underlying this muscle diversification, we performed single‐cell RNA‐sequencing experiments and built a cell atlas of myoblasts associated with third instar larval wing disc. Our analysis identified distinct transcriptional signatures for IFM and DFM myoblasts that underlie the molecular basis of their divergence. The atlas further revealed various states of differentiation of myoblasts, thus illustrating previously unappreciated spatial and temporal heterogeneity among them. We identified and validated novel markers for both IFM and DFM myoblasts at various states of differentiation by immunofluorescence and genetic cell‐tracing experiments. Finally, we performed a systematic genetic screen using a panel of markers from the reference cell atlas as an entry point and found a novel gene, Amalgam which is functionally important in muscle development. Our work provides a framework for leveraging scRNA‐seq for gene discovery and details a strategy that can be applied to other scRNA‐seq datasets.

Single cell transcriptomic landscapes of pattern formation, proliferation and growth in Drosophila wing imaginal discs

Organ formation relies on the orchestration of pattern formation, proliferation and growth during development. How these processes are integrated at the individual cell level remains unclear. In the past decades, studies using Drosophila wing imaginal discs as a model system have provided valuable insights into pattern formation, growth control and regeneration. Here, we provide single cell transcriptomic landscapes of pattern formation, proliferation and growth of wing imaginal discs. We found that patterning information is robustly maintained in the single cell transcriptomic data and can provide reference matrices for computationally mapping single cells into discrete spatial domains. Assignment of wing disc single cells to spatial subregions facilitates examination of patterning refinement processes. We also clustered single cells into different proliferation and growth states and evaluated the correlation between cell proliferation/growth states and spatial patterning. Furthermore, single cell transcriptomic analyses allowed us to quantitatively examine disturbances of differentiation, proliferation and growth in a well-established tumor model. We provide a database to explore these datasets at http://drosophilayanlab-virtual-wingdisc.ust.hk:3838/v2/This article has an associated 'The people behind the papers' interview.

Gene expression atlas of a developing tissue by single cell expression correlation analysis

The Drosophila wing disc has been a fundamental model system for the discovery of key signaling pathways and for our understanding of developmental processes. However, a complete map of gene expression in this tissue is lacking. To obtain a complete gene expression atlas in the wing disc, we employed single-cell sequencing (scRNA-seq) and developed a new method for analyzing scRNA-seq data based on gene expression correlations rather than cell mapping. This enables us to compute expression maps for all detected genes in the wing disc and to discover 824 genes with spatially restricted expression patterns. This approach identifies both known and new clusters of genes with similar expression patterns and functional relevance. As proof of concept, we characterize the previously unstudied gene CG5151 and show that it regulates Wnt signaling. This novel method will enable the leveraging of scRNA-seq data for generating expression atlases of undifferentiated tissues during development.

Ion Channel Contributions to Wing Development in Drosophila melanogaster

During morphogenesis, cells communicate with each other to shape tissues and organs. Several lines of recent evidence indicate that ion channels play a key role in cellular signaling and tissue morphogenesis. However, little is known about the scope of specific ion-channel types that impinge upon developmental pathways. The Drosophila melanogaster wing is an excellent model in which to address this problem as wing vein patterning is acutely sensitive to changes in developmental pathways. We conducted a screen of 180 ion channels expressed in the wing using loss-of-function mutant and RNAi lines. Here we identify 44 candidates that significantly impacted development of the Drosophila melanogaster wing. Calcium, sodium, potassium, chloride, and ligand-gated cation channels were all identified in our screen, suggesting that a wide variety of ion channel types are important for development. Ion channels belonging to the pickpocket family, the ionotropic receptor family, and the bestrophin family were highly represented among the candidates of our screen. Seven new ion channels with human orthologs that have been implicated in human channelopathies were also identified. Many of the human orthologs of the channels identified in our screen are targets of common general anesthetics, anti-seizure and anti-hypertension drugs, as well as alcohol and nicotine. Our results confirm the importance of ion channels in morphogenesis and identify a number of ion channels that will provide the basis for future studies to understand the role of ion channels in development.

Synaptojanin regulates Hedgehog signalling by modulating phosphatidylinositol 4-phosphate levels

In Hedgehog (Hh) signalling, Hh ligand concentration gradient is effectively translated into a spatially distinct transcriptional program to give precisely controlled context dependent developmental outcomes. In the absence of Hh, the receptor Patched (Ptc) inhibits the signal transducer Smoothened (Smo) by maintaining low phosphatidylinositol 4-phosphate (PI(4)P) levels. Binding of Hh to its receptor Ptc promotes PI(4)P production, which in turn activates Smo. Using wingdiscs of Drosophila melanogaster, this study shows that Synaptojanin (Synj), a dual phosphatase, modulates PI(4)P levels and affects Smo activation, and thereby functions as an additional regulatory step in the Hh pathway. Reducing the levels of Synj in the wing-discs caused enhancement of a Hh dominant gain-of-function Moonrat phenotype in the adult wings. Synj downregulation augmented Hh signalling, which was associated with elevated PI(4)P levels and Smo activation. Synj did not control the absolute pathway activity but rather fine-tuned the response since its downregulation increased expression of decapentaplegic (dpp), a low-threshold target of the pathway while the high-threshold targets remained unaffected. This is the first report that identifies Synj as a negative regulator of Hh signalling, implying its importance and an additional regulatory step in Hh signal transduction.

Wound healing, calcium signaling, and other novel pathways are associated with the formation of butterfly eyespots

Background

One hypothesis surrounding the origin of novel traits is that they originate from the co-option of pre-existing genes or larger gene regulatory networks into novel developmental contexts. Insights into a trait’s evolutionary origins can, thus, be gained via identification of the genes underlying trait development, and exploring whether those genes also function in other developmental contexts. Here we investigate the set of genes associated with the development of eyespot color patterns, a trait that originated once within the Nymphalid family of butterflies. Although several genes associated with eyespot development have been identified, the eyespot gene regulatory network remains largely unknown.

Results

In this study, next-generation sequencing and transcriptome analyses were used to identify a large set of genes associated with eyespot development of Bicyclus anynana butterflies, at 3-6 h after pupation, prior to the differentiation of the color rings. Eyespot-associated genes were identified by comparing the transcriptomes of homologous micro-dissected wing tissues that either develop or do not develop eyespots in wild-type and a mutant line of butterflies, Spotty, with extra eyespots. Overall, 186 genes were significantly up and down-regulated in wing tissues that develop eyespots compared to wing tissues that do not. Many of the differentially expressed genes have yet to be annotated. New signaling pathways, including the Toll, Fibroblast Growth Factor (FGF), extracellular signal–regulated kinase (ERK) and/or Jun N-terminal kinase (JNK) signaling pathways are associated for the first time with eyespot development. In addition, several genes involved in wound healing and calcium signaling were also found to be associated with eyespots.

Conclusions

Overall, this study provides the identity of many new genes and signaling pathways associated with eyespots, and suggests that the ancient wound healing gene regulatory network may have been co-opted to cells at the center of the pattern to aid in eyespot origins. New transcription factors that may be providing different identities to distinct wing sectors, and genes with sexually dimorphic expression in the eyespots were also identified.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4175-7) contains supplementary material, which is available to authorized users.

A clip domain serine protease involved in moulting in the silkworm, Bombyx mori: cloning, characterization, expression patterns and functional analysis

Clip domain serine proteases (CLIPs), characterized by one or more conserved clip domains, are essential components of extracellular signalling cascades in various biological processes, especially in innate immunity and the embryonic development of insects. Additionally, CLIPs may have additional non-immune functions in insect development. In the present study, the clip domain serine protease gene Bombyx mori serine protease 95 (BmSP95), which encodes a 527-residue protein, was cloned from the integument of B. mori. Bioinformatics analysis indicated that BmSP95 is a typical CLIP of the subfamily D and possesses a clip domain at the N terminus, a trypsin-like serine protease (tryp_spc) domain at the C terminus and a conserved proline-rich motif between these two domains. At the transcriptional level, BmSP95 is expressed in the integument during moulting and metamorphosis, and the expression pattern is consistent with the fluctuating 20-hydroxyecdysone (20E) titre in B. mori. At the translational level, BmSP95 protein is synthesized in the epidermal cells, secreted as a zymogen and activated in the moulting fluid. Immunofluorescence revealed that BmSP95 is distributed into the old endocuticle in the moulting stage. The expression of BmSP95 was upregulated by 20E. Moreover, expression of BmSP95 was downregulated by pathogen infection. RNA interference-mediated silencing of BmSP95 led to delayed moulting from pupa to moth. These results suggest that BmSP95 is involved in integument remodelling during moulting and metamorphosis.

Molecular mapping and characterization of the silkworm apodal mutant

The morphological diversity of insects is important for their survival; in essence, it results from the differential expression of genes during development of the insect body. The silkworm apodal (ap) mutant has degraded thoracic legs making crawling and eating difficult and the female is sterile, which is an ideal subject for studying the molecular mechanisms of morphogenesis. Here, we confirmed that the infertility of ap female moths is a result of the degradation of the bursa copulatrix. Positional cloning of ap locus and expression analyses reveal that the Bombyx mori sister of odd and bowl (Bmsob) gene is a strong candidate for the ap mutant. The expression of Bmsob is down-regulated, while the corresponding Hox genes are up-regulated in the ap mutant compared to the wild type. Analyses with the dual luciferase assay present a declined activity of the Bmsob promoter in the ap mutant. Furthermore, we demonstrate that Bmsob can inhibit Hox gene expression directly and by suppressing the expression of other genes, including the BmDsp gene. The results of this study are an important contribution to our understanding of the diversification of insect body plan.

Drosophila melanogaster clip-domain serine proteases: Structure, function and regulation

Mammalian chymotrypsin-like serine proteases (SPs) are one of the best-studied family of enzymes with roles in a wide range of physiological processes, including digestion, blood coagulation, fibrinolysis and humoral immunity. Extracellular SPs can form cascades, in which one protease activates the zymogen of the next protease in the chain, to amplify physiological or pathological signals. These extracellular SPs are generally multi-domain proteins, with pro-domains that are involved in protein-protein interactions critical for the sequential organization of the cascades, the control of their intensity and their proper localization. Far less is known about invertebrate SPs than their mammalian counterparts. In insect genomes, SPs and their proteolytically inactive homologs (SPHs) constitute large protein families. In addition to the chymotrypsin fold, many of these proteins contain additional structural domains, often with conserved mammalian orthologues. However, the largest group of arthropod SP regulatory modules is the clip domains family, which has only been identified in arthropods. The clip-domain SPs are extracellular and have roles in the immune response and embryonic development. The powerful reverse-genetics tools in Drosophila melanogaster have been essential to identify the functions of clip-SPs and their organization in sequential cascades. This review focuses on the current knowledge of Drosophila clip-SPs and presents, when necessary, data obtained in other insect models. We will first cover the biochemical and structural features of clip domain SPs and SPHs. Clip-SPs are implicated in three main biological processes: the control of the dorso-ventral patterning during embryonic development; the activation of the Toll-mediated response to microbial infections and the prophenoloxydase cascade, which triggers melanization. Finally, we review the regulation of SPs and SPHs, from specificity of activation to inhibition by endogenous or pathogen-encoded inhibitors.

Tumor Suppressor WWOX Contributes to the Elimination of Tumorigenic Cells in Drosophila melanogaster

WWOX is a >1Mb gene spanning FRA16D Common Chromosomal Fragile Site, a region of DNA instability in cancer. Consequently, altered WWOX levels have been observed in a wide variety of cancers. In vitro studies have identified a large number and variety of potential roles for WWOX. Although its normal role in vivo and functional contribution to cancer have not been fully defined, WWOX does have an integral role in metabolism and can suppress tumor growth. Using Drosophila melanogaster as an in vivo model system, we find that WWOX is a modulator of TNFα/Egr-mediated cell death. We found that altered levels of WWOX can modify phenotypes generated by low level ectopic expression of TNFα/Egr and this corresponds to altered levels of Caspase 3 activity. These results demonstrate an in vivo role for WWOX in promoting cell death. This form of cell death is accompanied by an increase in levels of reactive oxygen species, the regulation of which we have previously shown can also be modified by altered WWOX activity. We now hypothesise that, through regulation of reactive oxygen species, WWOX constitutes a link between alterations in cellular metabolism observed in cancer cells and their ability to evade normal cell death pathways. We have further shown that WWOX activity is required for the efficient removal of tumorigenic cells from a developing epithelial tissue. Together these results provide a molecular basis for the tumor suppressor functions of WWOX and the better prognosis observed in cancer patients with higher levels of WWOX activity. Understanding the conserved cellular pathways to which WWOX contributes provides novel possibilities for the development of therapeutic approaches to restore WWOX function in cancer.

A large-scale, in vivo transcription factor screen defines bivalent chromatin as a key property of regulatory factors mediating Drosophila wing development

Transcription factors (TFs) are key regulators of cell fate. The estimated 755 genes that encode DNA binding domain-containing proteins comprise ∼ 5% of all Drosophila genes. However, the majority has remained uncharacterized so far due to the lack of proper genetic tools. We generated 594 site-directed transgenic Drosophila lines that contain integrations of individual UAS-TF constructs to facilitate spatiotemporally controlled misexpression in vivo. All transgenes were expressed in the developing wing, and two-thirds induced specific phenotypic defects. In vivo knockdown of the same genes yielded a phenotype for 50%, with both methods indicating a great potential for misexpression to characterize novel functions in wing growth, patterning, and development. Thus, our UAS-TF library provides an important addition to the genetic toolbox of Drosophila research, enabling the identification of several novel wing development-related TFs. In parallel, we established the chromatin landscape of wing imaginal discs by ChIP-seq analyses of five chromatin marks and RNA Pol II. Subsequent clustering revealed six distinct chromatin states, with two clusters showing enrichment for both active and repressive marks. TFs that carry such "bivalent" chromatin are highly enriched for causing misexpression phenotypes in the wing, and analysis of existing expression data shows that these TFs tend to be differentially expressed across the wing disc. Thus, bivalently marked chromatin can be used as a marker for spatially regulated TFs that are functionally relevant in a developing tissue.

Large-scale imaginal disc sorting: A protocol for "omics"-approaches

Imaginal discs, especially the wing imaginal disc, are powerful model systems to study organ development. The traditional methods to analyze wing imaginal discs depend on the laborious and time-consuming dissection of larvae. "Omics"-based approaches, such as RNA-seq, ChIP-seq, proteomics and lipidomics, offer new opportunities for the systems-level investigation of organ development. However, it is impractical to manually isolate the required starting material. This is even more problematic when experiments strive for enhanced temporal and spatial resolution. The mass isolation workflow discussed in this review, solves this problem. The semi-automated sorting of 1000 wing imaginal discs in less than 3h forms the basis of a workflow that can be connected to biochemical analyses of organ patterning and growth. In addition to the mass isolation workflow we briefly describe key "omics" technologies and their applications. The combination of mass isolation and "omics"-approaches ensures that the wing imaginal disc will continue to be a key model organ for studying developmental processes, both on the genetic, but increasingly also on the biochemical level.