FlyBase: establishing a Gene Group resource for Drosophila melanogaster

Drosophila models used to simulate human ATP1A1 gene mutations that cause Charcot-Marie-Tooth type 2 disease and refractory seizures

JOURNAL/nrgr/04.03/01300535-202501000-00034/figure1/v/2024-05-14T021156Z/r/image-tiff Certain amino acids changes in the human Na+/K+-ATPase pump, ATPase Na+/K+ transporting subunit alpha 1 (ATP1A1), cause Charcot-Marie-Tooth disease type 2 (CMT2) disease and refractory seizures. To develop in vivo models to study the role of Na+/K+-ATPase in these diseases, we modified the Drosophila gene homolog, Atpα, to mimic the human ATP1A1 gene mutations that cause CMT2. Mutations located within the helical linker region of human ATP1A1 (I592T, A597T, P600T, and D601F) were simultaneously introduced into endogenous DrosophilaAtpα by CRISPR/Cas9-mediated genome editing, generating the AtpαTTTF model. In addition, the same strategy was used to generate the corresponding single point mutations in flies (AtpαI571T, AtpαA576T, AtpαP579T, and AtpαD580F). Moreover, a deletion mutation (Atpαmut) that causes premature termination of translation was generated as a positive control. Of these alleles, we found two that could be maintained as homozygotes (AtpαI571T and AtpαP579T). Three alleles (AtpαA576T, AtpαP579 and AtpαD580F) can form heterozygotes with the Atpαmut allele. We found that the Atpα allele carrying these CMT2-associated mutations showed differential phenotypes in Drosophila. Flies heterozygous for AtpαTTTF mutations have motor performance defects, a reduced lifespan, seizures, and an abnormal neuronal morphology. These Drosophila models will provide a new platform for studying the function and regulation of the sodium-potassium pump.

EMC1 Is Required for the Sarcoplasmic Reticulum and Mitochondrial Functions in the Drosophila Muscle

EMC1 is part of the endoplasmic reticulum (ER) membrane protein complex, whose functions include the insertion of transmembrane proteins into the ER membrane, ER-mitochondria contact, and lipid exchange. Here, we show that the Drosophila melanogaster EMC1 gene is expressed in the somatic musculature and the protein localizes to the sarcoplasmic reticulum (SR) network. Muscle-specific EMC1 RNAi led to severe motility defects and partial late pupae/early adulthood lethality, phenotypes that are rescued by co-expression with an EMC1 transgene. Motility impairment in EMC1-depleted flies was associated with aberrations in muscle morphology in embryos, larvae, and adults, including tortuous and misaligned fibers with reduced size and weakness. They were also associated with an altered SR network, cytosolic calcium overload, and mitochondrial dysfunction and dysmorphology that impaired membrane potential and oxidative phosphorylation capacity. Genes coding for ER stress sensors, mitochondrial biogenesis/dynamics, and other EMC components showed altered expression and were mostly rescued by the EMC1 transgene expression. In conclusion, EMC1 is required for the SR network's mitochondrial integrity and influences underlying programs involved in the regulation of muscle mass and shape. We believe our data can contribute to the biology of human diseases caused by EMC1 mutations.

A genome-wide CRISPR screen in Anopheles mosquito cells identifies essential genes and required components of clodronate liposome function

Anopheles mosquitoes are the sole vector of human malaria, the most burdensome vector-borne disease worldwide. Strategies aimed at reducing mosquito populations and limiting their ability to transmit disease show the most promise for disease control. Therefore, gaining an improved understanding of mosquito biology, and specifically that of the immune response, can aid efforts to develop new approaches that limit malaria transmission. Here, we use a genome-wide CRISPR screening approach for the first time in mosquito cells to identify essential genes in Anopheles and identify genes for which knockout confers resistance to clodronate liposomes, which have been widely used in mammals and arthropods to ablate immune cells. In the essential gene screen, we identified a set of 1280 Anopheles genes that are highly enriched for genes involved in fundamental cell processes. For the clodronate liposome screen, we identified several candidate resistance factors and confirm their roles in the uptake and processing of clodronate liposomes through in vivo validation in Anopheles gambiae, providing new mechanistic detail of phagolysosome formation and clodronate liposome function. In summary, we demonstrate the application of a genome-wide CRISPR knockout platform in a major malaria vector and the identification of genes that are important for fitness and immune-related processes.

Long non-coding RNAs involved in Drosophila development and regeneration

The discovery of functional long non-coding RNAs (lncRNAs) changed their initial concept as transcriptional noise. LncRNAs have been identified as regulators of multiple biological processes, including chromatin structure, gene expression, splicing, mRNA degradation, and translation. However, functional studies of lncRNAs are hindered by the usual lack of phenotypes upon deletion or inhibition. Here, we used Drosophila imaginal discs as a model system to identify lncRNAs involved in development and regeneration. We examined a subset of lncRNAs expressed in the wing, leg, and eye disc development. Additionally, we analyzed transcriptomic data from regenerating wing discs to profile the expression pattern of lncRNAs during tissue repair. We focused on the lncRNA CR40469, which is upregulated during regeneration. We generated CR40469 mutant flies that developed normally but showed impaired wing regeneration upon cell death induction. The ability of these mutants to regenerate was restored by the ectopic expression of CR40469. Furthermore, we found that the lncRNA CR34335 has a high degree of sequence similarity with CR40469 and can partially compensate for its function during regeneration in the absence of CR40469. Our findings point to a potential role of the lncRNA CR40469 in trans during the response to damage in the wing imaginal disc.

Misregulation of bromotyrosine compromises fertility in male Drosophila

Biological regulation often depends on reversible reactions such as phosphorylation, acylation, methylation, and glycosylation, but rarely halogenation. A notable exception is the iodination and deiodination of thyroid hormones. Here, we report detection of bromotyrosine and its subsequent debromination during Drosophila spermatogenesis. Bromotyrosine is not evident when Drosophila express a native flavin-dependent dehalogenase that is homologous to the enzyme responsible for iodide salvage from iodotyrosine in mammals. Deletion or suppression of the dehalogenase-encoding condet (cdt) gene in Drosophila allows bromotyrosine to accumulate with no detectable chloro- or iodotyrosine. The presence of bromotyrosine in the cdt mutant males disrupts sperm individualization and results in decreased fertility. Transgenic expression of the cdt gene in late-staged germ cells rescues this defect and enhances tolerance of male flies to bromotyrosine. These results are consistent with reversible halogenation affecting Drosophila spermatogenesis in a process that had previously eluded metabolomic, proteomic, and genomic analyses.

The alpha-ketoacid dehydrogenase complexes of Drosophila melanogaster

The conserved family of alpha-ketoacid dehydrogenase complexes (AKDHCs) catalyze essential reactions in central metabolism and their dysregulation is implicated in several human diseases. Drosophila melanogaster provides an excellent model system to study the genetics and functions of these complexes. However, a systematic account of Drosophila AKDHCs and their composition has been lacking. Here, I identify and classify the genes encoding all Drosophila AKDHC subunits, update their functional annotations and integrate this work into the FlyBase database.

A new experimental evidence-weighted signaling pathway resource in FlyBase

Research in model organisms is central to the characterization of signaling pathways in multicellular organisms. Here, we present the comprehensive and systematic curation of 17 Drosophila signaling pathways using the Gene Ontology framework to establish a dynamic resource that has been incorporated into FlyBase, providing visualization and data integration tools to aid research projects. By restricting to experimental evidence reported in the research literature and quantifying the amount of such evidence for each gene in a pathway, we captured the landscape of empirical knowledge of signaling pathways in Drosophila.

Revisiting bicoid function: complete inactivation reveals an additional fundamental role in Drosophila egg geometry specification

INTRODUCTION

The bicoid (bcd) gene in Drosophila has served as a paradigm for a morphogen in textbooks for decades. Discovered in 1986 as a mutation affecting anterior development in the embryo, its expression pattern as a protein gradient later confirmed the prediction from transplantation experiments. These experiments suggested that the protein fulfills the criteria of a true morphogen, with the existence of a homeodomain crucial for activation of genes along the anterior-posterior axis, based on the concentration of the morphogen. The bcd gene undergoes alternative splicing, resulting in, among other isoforms, a small and often neglected isoform with low abundance, which lacks the homeodomain, termed small bicoid (smbcd). Most importantly, all known classical strong bcd alleles used in the past to determine bcd function apparently do not affect the function of this isoform.

RESULTS

To overcome the uncertainty regarding which isoform regulates what, I removed the bcd locus entirely using CRISPR technology. bcdCRISPR eggs exhibited a short and round appearance. The phenotype could be ascribed to smbcd because all bcd alleles affecting the function of the major transcript, termed large bicoid (lgbcd) showed normally sized eggs. Several patterning genes for the embryo showed expression in the oocyte, and their expression patterns were altered in bcdCRISPR oocytes. In bcdCRISPR embryos, all downstream segmentation genes showed altered expression patterns, consistent with the expression patterns in "classical" alleles; however, due to the altered egg geometry resulting in fewer blastoderm nuclei, additional constraints came into play, further affecting their expression patterns.

CONCLUSIONS

This study unveils a novel and fundamental role of bcd in shaping the egg's geometry. This discovery demands a comprehensive revision of our understanding of this important patterning gene and prompts a reevaluation of past experiments conducted under the assumption that bcd mutants were bcdnull-mutants.

Finding information about uncharacterized Drosophila melanogaster genes

Genes that have been identified in the genome but remain uncharacterized with regards to function offer an opportunity to uncover novel biological information. Novelty is exciting but can also be a barrier. If nothing is known, how does one start planning and executing experiments? Here, we provide a recommended information-mining workflow and a corresponding guide to accessing information about uncharacterized Drosophila melanogaster genes, such as those assigned only a systematic coding gene identifier. The available information can provide insights into where and when the gene is expressed, what the function of the gene might be, whether there are similar genes in other species, whether there are known relationships to other genes, and whether any other features have already been determined. In addition, available information about relevant reagents can inspire and facilitate experimental studies. Altogether, mining available information can help prioritize genes for further study, as well as provide starting points for experimental assays and other analyses.

A new experimental evidence-weighted signaling pathway resource in FlyBase

Research in model organisms is central to the characterization of signaling pathways in multicellular organisms. Here, we present the systematic curation of 17 Drosophila signaling pathways using the Gene Ontology framework to establish a comprehensive and dynamic resource that has been incorporated into FlyBase, providing visualization and data integration tools to aid research projects. By restricting to experimental evidence reported in the research literature and quantifying the amount of such evidence for each gene in a pathway, we captured the landscape of empirical knowledge of signaling pathways in Drosophila .

Summary statement

Comprehensive curation of Drosophila signaling pathways and new visual displays of the pathways provides a new FlyBase resource for researchers, and new insights into signaling pathway architecture.

A new metabolic model of Drosophila melanogaster and the integrative analysis of Parkinson's disease

High conservation of the disease-associated genes between flies and humans facilitates the common use of Drosophila melanogaster to study metabolic disorders under controlled laboratory conditions. However, metabolic modeling studies are highly limited for this organism. We here report a comprehensively curated genome-scale metabolic network model of Drosophila using an orthology-based approach. The gene coverage and metabolic information of the draft model derived from a reference human model were expanded via Drosophila-specific KEGG and MetaCyc databases, with several curation steps to avoid metabolic redundancy and stoichiometric inconsistency. Furthermore, we performed literature-based curations to improve gene-reaction associations, subcellular metabolite locations, and various metabolic pathways. The performance of the resulting Drosophila model (8,230 reactions, 6,990 metabolites, and 2,388 genes), iDrosophila1 (https://github.com/SysBioGTU/iDrosophila), was assessed using flux balance analysis in comparison with the other currently available fly models leading to superior or comparable results. We also evaluated the transcriptome-based prediction capacity of iDrosophila1, where differential metabolic pathways during Parkinson's disease could be successfully elucidated. Overall, iDrosophila1 is promising to investigate system-level metabolic alterations in response to genetic and environmental perturbations.

Pervasive SUMOylation of heterochromatin and piRNA pathway proteins

Genome regulation involves complex protein interactions that are often mediated through post-translational modifications (PTMs). SUMOylation-modification by the small ubiquitin-like modifier (SUMO)-has been implicated in numerous essential processes in eukaryotes. In Drosophila, SUMO is required for viability and fertility, with its depletion from ovaries leading to heterochromatin loss and ectopic transposon and gene activation. Here, we developed a proteomics-based strategy to uncover the Drosophila ovarian "SUMOylome," which revealed that SUMOylation is widespread among proteins involved in heterochromatin regulation and different aspects of the Piwi-interacting small RNA (piRNA) pathway that represses transposons. Furthermore, we show that SUMOylation of several piRNA pathway proteins occurs in a Piwi-dependent manner. Together, these data highlight broad implications of protein SUMOylation in epigenetic regulation and indicate novel roles of this modification in the cellular defense against genomic parasites. Finally, this work provides a resource for the study of SUMOylation in other biological contexts in the Drosophila model.

PANGEA: a new gene set enrichment tool for Drosophila and common research organisms

Gene set enrichment analysis (GSEA) plays an important role in large-scale data analysis, helping scientists discover the underlying biological patterns over-represented in a gene list resulting from, for example, an 'omics' study. Gene Ontology (GO) annotation is the most frequently used classification mechanism for gene set definition. Here we present a new GSEA tool, PANGEA (PAthway, Network and Gene-set Enrichment Analysis; https://www.flyrnai.org/tools/pangea/), developed to allow a more flexible and configurable approach to data analysis using a variety of classification sets. PANGEA allows GO analysis to be performed on different sets of GO annotations, for example excluding high-throughput studies. Beyond GO, gene sets for pathway annotation and protein complex data from various resources as well as expression and disease annotation from the Alliance of Genome Resources (Alliance). In addition, visualizations of results are enhanced by providing an option to view network of gene set to gene relationships. The tool also allows comparison of multiple input gene lists and accompanying visualisation tools for quick and easy comparison. This new tool will facilitate GSEA for Drosophila and other major model organisms based on high-quality annotated information available for these species.

Tropical super flies: Integrating Cas9 into Drosophila ananassae and its phenotypic effects

Ectotherms such as insects are animals whose body temperature largely depends on ambient temperature and temperature variations provide a selection pressure affecting the geographical distribution of these species. However, over the course of evolution, some insect species managed to colonize environments characterized by various temperature ranges. Therefore, insects provide an excellent study system to investigate the basis of adaptation to temperature changes and extremes. We are generally using the vinegar fly Drosophila ananassae as a model system to investigate the genetic basis of cold tolerance. This species has expanded from its tropical ancestral range to more temperate regions resulting in a cosmopolitan, domestic distribution. Previously, we identified candidate genes significantly associated with cold tolerance in this species. We now established molecular genetic tools to assess the function of these genes. Using CRISPR/Cas9 methodology for genome editing and the PiggyBac system, the Cas9 enzyme was successfully integrated into the genome of three fly strains with different levels of cold tolerance. We further report on preliminary findings that the Cas9 integration itself did not have a consistent effect on tolerance to cold. In conclusion, we offer with our study the molecular tools that allow studying stress-related candidate genes in D. ananassae in the future. In addition, we point out and provide guidance on the challenges that come with genome editing in a non-model species.

On the identification of potential novel therapeutic targets for spinocerebellar ataxia type 1 (SCA1) neurodegenerative disease using EvoPPI3

EvoPPI (http://evoppi.i3s.up.pt), a meta-database for protein-protein interactions (PPI), has been upgraded (EvoPPI3) to accept new types of data, namely, PPI from patients, cell lines, and animal models, as well as data from gene modifier experiments, for nine neurodegenerative polyglutamine (polyQ) diseases caused by an abnormal expansion of the polyQ tract. The integration of the different types of data allows users to easily compare them, as here shown for Ataxin-1, the polyQ protein involved in spinocerebellar ataxia type 1 (SCA1) disease. Using all available datasets and the data here obtained for Drosophila melanogaster wt and exp Ataxin-1 mutants (also available at EvoPPI3), we show that, in humans, the Ataxin-1 network is much larger than previously thought (380 interactors), with at least 909 interactors. The functional profiling of the newly identified interactors is similar to the ones already reported in the main PPI databases. 16 out of 909 interactors are putative novel SCA1 therapeutic targets, and all but one are already being studied in the context of this disease. The 16 proteins are mainly involved in binding and catalytic activity (mainly kinase activity), functional features already thought to be important in the SCA1 disease.

PANGEA: A New Gene Set Enrichment Tool for Drosophila and Common Research Organisms

Gene set enrichment analysis (GSEA) plays an important role in large-scale data analysis, helping scientists discover the underlying biological patterns over-represented in a gene list resulting from, for example, an 'omics' study. Gene Ontology (GO) annotation is the most frequently used classification mechanism for gene set definition. Here we present a new GSEA tool, PANGEA (PAthway, Network and Gene-set Enrichment Analysis; https://www.flyrnai.org/tools/pangea/ ), developed to allow a more flexible and configurable approach to data analysis using a variety of classification sets. PANGEA allows GO analysis to be performed on different sets of GO annotations, for example excluding high-throughput studies. Beyond GO, gene sets for pathway annotation and protein complex data from various resources as well as expression and disease annotation from the Alliance of Genome Resources (Alliance). In addition, visualisations of results are enhanced by providing an option to view network of gene set to gene relationships. The tool also allows comparison of multiple input gene lists and accompanying visualisation tools for quick and easy comparison. This new tool will facilitate GSEA for Drosophila and other major model organisms based on high-quality annotated information available for these species.

Exploring FlyBase Data Using QuickSearch

FlyBase (www.flybase.org) is the primary online database of genetic, genomic, and functional information about Drosophila melanogaster. The long and rich history of Drosophila research, combined with recent surges in genomic-scale and high-throughput technologies, means that FlyBase now houses a huge quantity of data. Researchers need to be able to query these data rapidly and intuitively, and the QuickSearch tool has been designed to meet these needs. This tool is conveniently located on the FlyBase homepage and is organized into a series of simple tabbed interfaces that cover the major data and annotation classes within the database. This article describes the functionality of all aspects of the QuickSearch tool. With this knowledge, FlyBase users will be equipped to take full advantage of all QuickSearch features and thereby gain improved access to data relevant to their research. © 2023 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: Using the "Search FlyBase" tab of QuickSearch Basic Protocol 2: Using the "Data Class" tab of QuickSearch Basic Protocol 3: Using the "References" tab of QuickSearch Basic Protocol 4: Using the "Gene Groups" tab of QuickSearch Basic Protocol 5: Using the "Pathways" tab of QuickSearch Basic Protocol 6: Using the "GO" tab of QuickSearch Basic Protocol 7: Using the "Protein Domains" tab of QuickSearch Basic Protocol 8: Using the "Expression" tab of QuickSearch Basic Protocol 9: Using the "GAL4 etc" tab of QuickSearch Basic Protocol 10: Using the "Phenotype" tab of QuickSearch Basic Protocol 11: Using the "Human Disease" tab of QuickSearch Basic Protocol 12: Using the "Homologs" tab of QuickSearch Support Protocol 1: Managing FlyBase hit lists.

Canalization of Phenotypes-When the Transcriptome is Constantly but Weakly Perturbed

Recent studies have increasingly pointed to microRNAs (miRNAs) as the agent of gene regulatory network (GRN) stabilization as well as developmental canalization against constant but small environmental perturbations. To analyze mild perturbations, we construct a Dicer-1 knockdown line (dcr-1 KD) in Drosophila that modestly reduces all miRNAs by, on average, ∼20%. The defining characteristic of stabilizers is that, when their capacity is compromised, GRNs do not change their short-term behaviors. Indeed, even with such broad reductions across all miRNAs, the changes in the transcriptome are very modest during development in stable environment. By comparison, broad knockdowns of other regulatory genes (esp. transcription factors) by the same method should lead to drastic changes in the GRNs. The consequence of destabilization may thus be in long-term development as postulated by the theory of canalization. Flies with modest miRNA reductions may gradually deviate from the developmental norm, resulting in late-stage failures such as shortened longevity. In the optimal culture condition, the survival to adulthood is indeed normal in the dcr-1 KD line but, importantly, adult longevity is reduced by ∼90%. When flies are stressed by high temperature, dcr-1 KD induces lethality earlier in late pupation and, as the perturbations are shifted earlier, the affected stages are shifted correspondingly. Hence, in late stages of development with deviations piling up, GRN would be increasingly in need of stabilization. In conclusion, miRNAs appear to be a solution to weak but constant environmental perturbations.

Identification and Characterization of Small RNA Markers of Age in the Blow Fly Cochliomyia macellaria (Fabricius) (Diptera: Calliphoridae)

Blow fly development is important in decomposition ecology, agriculture, and forensics. Much of the impact of these species is from immature samples, thus knowledge of their development is important to enhance or ameliorate their effects. One application of this information is the estimation of immature insect age to provide temporal information for death investigations. While traditional markers of age such as stage and size are generally accurate, they lack precision in later developmental stages. We used miRNA sequencing to measure miRNA expression, throughout development, of the secondary screwworm, Cochliomyia macellaria (Fabricius) (Diptera: Calliphoridae) and identified 217 miRNAs present across the samples. Ten were identified to be significantly differentially expressed in larval samples and seventeen were found to be significantly differentially expressed in intrapuparial samples. Twenty-eight miRNAs were identified to be differentially expressed between sexes. Expression patterns of two miRNAs, miR-92b and bantam, were qPCR-validated in intrapuparial samples; these and likely food-derived miRNAs appear to be stable markers of age in C. macellaria. Our results support the use of miRNAs for developmental markers of age and suggest further investigations across species and under a range of abiotic and biotic conditions.

Systematic identification of tRNA genes in Drosophila melanogaster

Transfer RNAs (tRNAs) are ubiquitous adapter molecules that link specific codons in messenger RNA (mRNA) with their corresponding amino acids during protein synthesis. The tRNA genes of Drosophila have been investigated for over half a century but have lacked systematic identification and nomenclature. Here, we review and integrate data within FlyBase and the Genomic tRNA Database (GtRNAdb) to identify the full complement of tRNA genes in the D. melanogaster nuclear and mitochondrial genomes. We apply a logical and informative nomenclature to all tRNA genes, and provide an overview of their characteristics and genomic features.

Phylogenetic annotation of Drosophila melanogaster heat shock protein 70 genes

The traditional classification of stress-inducible 70 kDa heat shock protein (Hsp70) and heat shock cognate (Hsc70) requires a revision because of a recent finding that neither of them constitutes a monophyletic gene family. Here we inferred a phylogenetic relationship among Drosophila melanogaster Hsp70 family members. D. melanogaster Hsp70 family members were separated into four known metazoan Hsp70 lineages: cytosolic A, cytosolic B, endoplasmic reticulum, and mitochondria. Hsc70s sporadically distributed in the phylogenetic tree, indicating their paraphyletic origin. Detailed sequence inspection found several motifs that support the phylogenetic analysis. Taken together, we propose new aliases of D. melanogaster Hsp70 family members based on their evolutionary history.

FlyBase: a guided tour of highlighted features

FlyBase provides a centralized resource for the genetic and genomic data of Drosophila melanogaster. As FlyBase enters our fourth decade of service to the research community, we reflect on our unique aspects and look forward to our continued collaboration with the larger research and model organism communities. In this study, we emphasize the dedicated reports and tools we have constructed to meet the specialized needs of fly researchers but also to facilitate use by other research communities. We also highlight ways that we support the fly community, including an external resources page, help resources, and multiple avenues by which researchers can interact with FlyBase.

Invertebrate Model Organisms as a Platform to Investigate Rare Human Neurological Diseases

Patients suffering from rare human diseases often go through a painful journey for finding a definite molecular diagnosis prerequisite of appropriate cures. With a novel variant isolated from a single patient, determination of its pathogenicity to end such "diagnostic odyssey" requires multi-step processes involving experts in diverse areas of interest, including clinicians, bioinformaticians and research scientists. Recent efforts in building large-scale genomic databases and in silico prediction platforms have facilitated identification of potentially pathogenic variants causative of rare human diseases of a Mendelian basis. However, the functional significance of individual variants remains elusive in many cases, thus requiring incorporation of versatile and rapid model organism (MO)-based platforms for functional analyses. In this review, the current scope of rare disease research is briefly discussed. In addition, an overview of invertebrate MOs for their key features relevant to rare neurological diseases is provided, with the characteristics of two representative invertebrate MOs, Drosophila melanogaster and Caenorhabditis elegans, as well as the challenges against them. Finally, recently developed research networks integrating these MOs in collaborative research are portraited with an array of bioinformatical analyses embedded. A comprehensive survey of MO-based research activities provided in this review will help us to design a wellstructured analysis of candidate genes or potentially pathogenic variants for their roles in rare neurological diseases in future.

A spatiotemporal atlas of the lepidopteran pest Helicoverpa armigera midgut provides insights into nutrient processing and pH regulation

Background

Caterpillars from the insect order Lepidoptera are some of the most widespread and destructive agricultural pests. Most of their impact is at the larval stage, where the midgut epithelium mediates the digestion and absorption of an astonishing amount of food. Although this tissue has been the subject of frequent investigation in Lepidoptera, a comprehensive expression atlas has yet to be generated.

Results

Here, we perform RNA-sequencing and proteomics on the gut of the polyphagous pest Helicoverpa armigera across, life stages, diet types, and compartments of the anterior-posterior axis. A striking relationship between the structural homology and expression pattern of a group of sugar transporters was observed in the early larval stages. Further comparisons were made among the spatial compartments of the midgut, which suggested a putative role for vATPases and SLC9 transporters in the generation of alkaline conditions in the H. armigera midgut.

Conclusions

This comprehensive resource will aid the scientific community in understanding lepidopteran gut physiology in unprecedented resolution. It is hoped that this study advances the understanding of the lepidopteran midgut and also facilitates functional work in this field.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08274-x.

FlyAtlas 2 in 2022: enhancements to the Drosophila melanogaster expression atlas

FlyAtlas 2 (flyatlas2.org) is a database and web application for studying the expression of the genes of Drosophila melanogaster in different tissues of adults and larvae. It is based on RNA-Seq data, and incorporates both genes encoding proteins and microRNAs. We have now completed the population of the database with 13 tissues from both male and female adults, five sex-specific tissues, and eight larval tissues. Larval garland cell nephrocytes have also been included. Major enhancements have been made to the application. First, a facility has been added for a 'Profile' search for genes with a similar pattern of tissue expression as a query gene. This may help establish the function of genes for which this is currently unknown. Second, a facility has been added dedicated to the larval midgut, where the difference in gene expression in the five regions of different pH can be explored. A variety of further improvements to the interface are described.

Using FlyBase: A Database of Drosophila Genes and Genetics

Since 1992, FlyBase has provided a freely available online database of information about the model organism Drosophila melanogaster. Data in FlyBase is curated manually from research papers as well as computationally from a variety of relevant sources, to serve as an information hub that enables and accelerates research discovery. This chapter aims to give users new to the database an overview of the layout and types of data available, as well as introducing some tools with which to access the data. More experienced users will find useful information about recent improvements and descriptions to enable more efficient navigation of the database.

Nuclear transport genes recurrently duplicate by means of RNA intermediates in Drosophila but not in other insects

Background

The nuclear transport machinery is involved in a well-known male meiotic drive system in Drosophila. Fast gene evolution and gene duplications have been major underlying mechanisms in the evolution of meiotic drive systems, and this might include some nuclear transport genes in Drosophila. So, using a comprehensive, detailed phylogenomic study, we examined 51 insect genomes for the duplication of the same nuclear transport genes.

Results

We find that most of the nuclear transport duplications in Drosophila are of a few classes of nuclear transport genes, RNA mediated and fast evolving. We also retrieve many pseudogenes for the Ran gene. Some of the duplicates are relatively young and likely contributing to the turnover expected for genes under strong but changing selective pressures. These duplications are potentially revealing what features of nuclear transport are under selection. Unlike in flies, we find only a few duplications when we study the Drosophila duplicated nuclear transport genes in dipteran species outside of Drosophila, and none in other insects.

Conclusions

These findings strengthen the hypothesis that nuclear transport gene duplicates in Drosophila evolve either as drivers or suppressors of meiotic drive systems or as other male-specific adaptations circumscribed to flies and involving a handful of nuclear transport functions.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08170-4.

Up-regulation of apoptotic- and cell survival-related gene pathways following exposures of western corn rootworm to B. thuringiensis crystalline pesticidal proteins in transgenic maize roots

BACKGROUND

Resistance of pest insect species to insecticides, including B. thuringiensis (Bt) pesticidal proteins expressed by transgenic plants, is a threat to global food security. Despite the western corn rootworm, Diabrotica virgifera virgifera, being a major pest of maize and having populations showing increasing levels of resistance to hybrids expressing Bt pesticidal proteins, the cell mechanisms leading to mortality are not fully understood.

RESULTS

Twenty unique RNA-seq libraries from the Bt susceptible D. v. virgifera inbred line Ped12, representing all growth stages and a range of different adult and larval exposures, were assembled into a reference transcriptome. Ten-day exposures of Ped12 larvae to transgenic Bt Cry3Bb1 and Gpp34/Tpp35Ab1 maize roots showed significant differential expression of 1055 and 1374 transcripts, respectively, compared to cohorts on non-Bt maize. Among these, 696 were differentially expressed in both Cry3Bb1 and Gpp34/Tpp35Ab1 maize exposures. Differentially-expressed transcripts encoded protein domains putatively involved in detoxification, metabolism, binding, and transport, were, in part, shared among transcripts that changed significantly following exposures to the entomopathogens Heterorhabditis bacteriophora and Metarhizium anisopliae. Differentially expressed transcripts in common between Bt and entomopathogen treatments encode proteins in general stress response pathways, including putative Bt binding receptors from the ATP binding cassette transporter superfamily. Putative caspases, pro- and anti-apoptotic factors, as well as endoplasmic reticulum (ER) stress-response factors were identified among transcripts uniquely up-regulated following exposure to either Bt protein.

CONCLUSIONS

Our study suggests that the up-regulation of genes involved in ER stress management and apoptotic progression may be important in determining cell fate following exposure of susceptible D. v. virgifera larvae to Bt maize roots. This study provides novel insights into insect response to Bt intoxication, and a possible framework for future investigations of resistance mechanisms.

E2F/Dp inactivation in fat body cells triggers systemic metabolic changes

The E2F transcription factors play a critical role in controlling cell fate. In Drosophila, the inactivation of E2F in either muscle or fat body results in lethality, suggesting an essential function for E2F in these tissues. However, the cellular and organismal consequences of inactivating E2F in these tissues are not fully understood. Here, we show that the E2F loss exerts both tissue-intrinsic and systemic effects. The proteomic profiling of E2F-deficient muscle and fat body revealed that E2F regulates carbohydrate metabolism, a conclusion further supported by metabolomic profiling. Intriguingly, animals with E2F-deficient fat body had a lower level of circulating trehalose and reduced storage of fat. Strikingly, a sugar supplement was sufficient to restore both trehalose and fat levels, and subsequently rescued animal lethality. Collectively, our data highlight the unexpected complexity of E2F mutant phenotype, which is a result of combining both tissue-specific and systemic changes that contribute to animal development.

Anti-Fibrotic Activity of an Antimicrobial Peptide in a Drosophila Model

Fibrotic lesions accompany several pathological conditions, including tumors. We show that expression of a dominant-active form of the Ras oncogene in Drosophila salivary glands (SGs) leads to redistribution of components of the basement membrane (BM) and fibrotic lesions. Similar to several types of mammalian fibrosis, the disturbed BM attracts clot components, including insect transglutaminase and phenoloxidase. SG epithelial cells show reduced apicobasal polarity accompanied by a loss of secretory activity. Both the fibrotic lesions and the reduced cell polarity are alleviated by ectopic expression of the antimicrobial peptide drosomycin (Drs), which also restores the secretory activity of the SGs. In addition to extracellular matrix components, both Drs and F-actin localize to fibrotic lesions.

Run or Die in the Evolution of New MicroRNAs-Testing the Red Queen Hypothesis on De Novo New Genes

The Red Queen hypothesis depicts evolution as the continual struggle to adapt. According to this hypothesis, new genes, especially those originating from nongenic sequences (i.e., de novo genes), are eliminated unless they evolve continually in adaptation to a changing environment. Here, we analyze two Drosophila de novo miRNAs that are expressed in a testis-specific manner with very high rates of evolution in their DNA sequence. We knocked out these miRNAs in two sibling species and investigated their contributions to different fitness components. We observed that the fitness contributions of miR-975 in Drosophila simulans seem positive, in contrast to its neutral contributions in D. melanogaster, whereas miR-983 appears to have negative contributions in both species, as the fitness of the knockout mutant increases. As predicted by the Red Queen hypothesis, the fitness difference of these de novo miRNAs indicates their different fates.

Subunits of the PBAP Chromatin Remodeler Are Capable of Mediating Enhancer-Driven Transcription in Drosophila

The chromatin remodeler SWI/SNF is an important participant in gene activation, functioning predominantly by opening the chromatin structure on promoters and enhancers. Here, we describe its novel mode of action in which SWI/SNF factors mediate the targeted action of an enhancer. We studied the functions of two signature subunits of PBAP subfamily, BAP170 and SAYP, in Drosophila. These subunits were stably tethered to a transgene reporter carrying the hsp70 core promoter. The tethered subunits mediate transcription of the reporter in a pattern that is generated by enhancers close to the insertion site in multiple loci throughout the genome. Both tethered SAYP and BAP170 recruit the whole PBAP complex to the reporter promoter. However, we found that BAP170-dependent transcription is more resistant to the depletion of other PBAP subunits, suggesting that BAP170 may play a more critical role in establishing enhancer-dependent transcription.

Testing Implications of the Omnigenic Model for the Genetic Analysis of Loci Identified through Genome-wide Association

Organismal phenotypes usually have a quantitative distribution, and their genetic architecture can be studied by genome-wide association (GWA) mapping approaches. In most of such studies, it has become clear that many genes of moderate or small effects contribute to the phenotype.^1-4 Hence, the attention has turned toward the loci falling below the GWA cut-off, which may contribute to the phenotype through modifier interactions with a set of core genes, as proposed in the omnigenic model.⁵ One can thus predict that both moderate effect GWA-derived candidate genes and randomly chosen genes should have a similar likelihood to affect a given phenotype when they are analyzed via gene disruption assays. We have tested this hypothesis by using an automated phenotyping system for Drosophila pupal phenotypes.^6,7 We first identified candidate genes for pupal length in a GWA based on the Drosophila Genetic Reference Panel (DGRP)^8,9 and showed that most of these candidate genes are indeed involved in the phenotype. We then randomly chose genes below a GWA significance threshold and found that three-quarters of them had also an effect on the trait with comparable effect sizes as the GWA candidate genes. We further tested the effects of these knockout lines on an independent behavioral pupal trait (pupation site choice) and found that a similar fraction had a significant effect as well. Our data thus confirm the implication that a large number of genes can influence independent quantitative traits.

Wolbachia affects reproduction in the spider mite Tetranychus truncatus (Acari: Tetranychidae) by regulating chorion protein S38-like and Rop

Wolbachia-induced reproductive regulation in hosts has been used to control pest populations, but little is known about the molecular mechanism underlying Wolbachia regulation of host genes. Here, reproductive regulation by Wolbachia in the spider mite Tetranychus truncatus was studied at the molecular level. Infection with Wolbachia resulted in decreasing oviposition and cytoplasmic incompatibility in T. truncatus. Further RNA-seq revealed genes regulated by Wolbachia in T. truncatus. Real-time quantitative polymerase chain reaction (qPCR) showed that genes, including chorion protein S38-like and Rop were down-regulated by Wolbachia. RNA interference (RNAi) of chorion protein S38-like and Rop in Wolbachia-uninfected T. truncatus decreased oviposition, which was consistent with Wolbachia-induced oviposition decrease. Interestingly, suppressing Rop in Wolbachia-infected T. truncatus led to increased Wolbachia titres in eggs; however, this did not occur after RNAi of chorion protein S38-like. This is the first study to show that chorion protein S38-like and Rop facilitate Wolbachia-mediated changes in T. truncatus fertility. In addition, RNAi of Rop turned the body colour of Wolbachia-uninfected T. truncatus black, which indicates that the role of Rop is not limited to the reproductive regulation of T. truncatus.

FlyBase: updates to the Drosophila melanogaster knowledge base

FlyBase (flybase.org) is an essential online database for researchers using Drosophila melanogaster as a model organism, facilitating access to a diverse array of information that includes genetic, molecular, genomic and reagent resources. Here, we describe the introduction of several new features at FlyBase, including Pathway Reports, paralog information, disease models based on orthology, customizable tables within reports and overview displays ('ribbons') of expression and disease data. We also describe a variety of recent important updates, including incorporation of a developmental proteome, upgrades to the GAL4 search tab, additional Experimental Tool Reports, migration to JBrowse for genome browsing and improvements to batch queries/downloads and the Fast-Track Your Paper tool.

A functional genomics screen identifying blood cell development genes in Drosophila by undergraduates participating in a course-based research experience

Undergraduate students participating in the UCLA Undergraduate Research Consortium for Functional Genomics (URCFG) have conducted a two-phased screen using RNA interference (RNAi) in combination with fluorescent reporter proteins to identify genes important for hematopoiesis in Drosophila. This screen disrupted the function of approximately 3500 genes and identified 137 candidate genes for which loss of function leads to observable changes in the hematopoietic development. Targeting RNAi to maturing, progenitor, and regulatory cell types identified key subsets that either limit or promote blood cell maturation. Bioinformatic analysis reveals gene enrichment in several previously uncharacterized areas, including RNA processing and export and vesicular trafficking. Lastly, the participation of students in this course-based undergraduate research experience (CURE) correlated with increased learning gains across several areas, as well as increased STEM retention, indicating that authentic, student-driven research in the form of a CURE represents an impactful and enriching pedagogical approach.

Widely Used Mutants of eiger, Encoding the Drosophila Tumor Necrosis Factor, Carry Additional Mutations in the NimrodC1 Phagocytosis Receptor

The process of apoptosis in epithelia involves activation of caspases, delamination of cells, and degradation of cellular components. Corpses and cellular debris are then rapidly cleared from the tissue by phagocytic blood cells. In studies of the Drosophila TNF, Eiger (Egr) and cell death in wing imaginal discs, the epithelial primordia of fly wings, we noticed that dying cells appeared to transiently accumulate in egr3 mutant wing discs, raising the possibility that their phagocytic engulfment by hemocytes was impaired. Further investigation revealed that lymph glands and circulating hemocytes from egr3 mutant larvae were completely devoid of NimC1 staining, a marker of phagocytic hemocytes. Genome sequencing uncovered mutations in the NimC1 coding region that are predicted to truncate the NimC1 protein before its transmembrane domain, and provide an explanation for the lack of NimC staining. The work that we report here demonstrates the presence of these NimC1 mutations in the widely used egr3 mutant, its sister allele, egr1 , and its parental strain, Regg1GS9830 As the egr3 and egr1 alleles have been used in numerous studies of immunity and cell death, it may be advisable to re-evaluate their associated phenotypes.

bsAS, an antisense long non-coding RNA, essential for correct wing development through regulation of blistered/DSRF isoform usage

Natural Antisense Transcripts (NATs) are long non-coding RNAs (lncRNAs) that overlap coding genes in the opposite strand. NATs roles have been related to gene regulation through different mechanisms, including post-transcriptional RNA processing. With the aim to identify NATs with potential regulatory function during fly development, we generated RNA-Seq data in Drosophila developing tissues and found bsAS, one of the most highly expressed lncRNAs in the fly wing. bsAS is antisense to bs/DSRF, a gene involved in wing development and neural processes. bsAS plays a crucial role in the tissue specific regulation of the expression of the bs/DSRF isoforms. This regulation is essential for the correct determination of cell fate during Drosophila development, as bsAS knockouts show highly aberrant phenotypes. Regulation of bs isoform usage by bsAS is mediated by specific physical interactions between the promoters of these two genes, which suggests a regulatory mechanism involving the collision of RNA polymerases transcribing in opposite directions. Evolutionary analysis suggests that bsAS NAT emerged simultaneously to the long-short isoform structure of bs, preceding the emergence of wings in insects.

Magnesium efflux from Drosophila Kenyon cells is critical for normal and diet-enhanced long-term memory

Dietary magnesium (Mg²⁺) supplementation can enhance memory in young and aged rats. Memory-enhancing capacity was largely ascribed to increases in hippocampal synaptic density and elevated expression of the NR2B subunit of the NMDA-type glutamate receptor. Here we show that Mg²⁺ feeding also enhances long-term memory in Drosophila. Normal and Mg²⁺-enhanced fly memory appears independent of NMDA receptors in the mushroom body and instead requires expression of a conserved CNNM-type Mg²⁺-efflux transporter encoded by the unextended (uex) gene. UEX contains a putative cyclic nucleotide-binding homology domain and its mutation separates a vital role for uex from a function in memory. Moreover, UEX localization in mushroom body Kenyon cells (KCs) is altered in memory-defective flies harboring mutations in cAMP-related genes. Functional imaging suggests that UEX-dependent efflux is required for slow rhythmic maintenance of KC Mg²⁺. We propose that regulated neuronal Mg²⁺ efflux is critical for normal and Mg²⁺-enhanced memory.

The proverbial saying ‘you are what you eat’ perfectly summarizes the concept that our diet can influence both our mental and physical health. We know that foods that are good for the heart, such as nuts, oily fish and berries, are also good for the brain. We know too that vitamins and minerals are essential for overall good health. But is there any evidence that increasing your intake of specific vitamins or minerals could help boost your brain power?

While it might sound almost too good to be true, there is some evidence that this is the case for at least one mineral, magnesium. Studies in rodents have shown that adding magnesium supplements to food improves how well the animals perform on memory tasks. Both young and old animals benefit from additional magnesium. Even elderly rodents with a condition similar to Alzheimer’s disease show less memory loss when given magnesium supplements. But what about other species?

Wu et al. now show that magnesium supplements also boost memory performance in fruit flies. One group of flies was fed with standard cornmeal for several days, while the other group received cornmeal supplemented with magnesium. Both groups were then trained to associate an odor with a food reward. Flies that had received the extra magnesium showed better memory for the odor when tested 24 hours after training.

Wu et al. show that magnesium improves memory in the flies via a different mechanism to that reported previously for rodents. In rodents, magnesium increased levels of a receptor protein for a brain chemical called glutamate. In fruit flies, by contrast, the memory boost depended on a protein that transports magnesium out of neurons. Mutant flies that lacked this transporter showed memory impairments. Unlike normal flies, those without the transporter showed no memory improvement after eating magnesium-enriched food. The results suggest that the transporter may help adjust magnesium levels inside brain cells in response to neural activity.

Humans produce four variants of this magnesium transporter, each encoded by a different gene. One of these transporters has already been implicated in brain development. The findings of Wu et al. suggest that the transporters may also act in the adult brain to influence cognition. Further studies are needed to test whether targeting the magnesium transporter could ultimately hold promise for treating memory impairments.

The genomes of a monogenic fly: views of primitive sex chromosomes

The production of male and female offspring is often determined by the presence of specific sex chromosomes which control sex-specific expression, and sex chromosomes evolve through reduced recombination and specialized gene content. Here we present the genomes of Chrysomya rufifacies, a monogenic blow fly (females produce female or male offspring, exclusively) by separately sequencing and assembling each type of female and the male. The genomes (> 25X coverage) do not appear to have any sex-linked Muller F elements (typical for many Diptera) and exhibit little differentiation between groups supporting the morphological assessments of C. rufifacies homomorphic chromosomes. Males in this species are associated with a unimodal coverage distribution while females exhibit bimodal coverage distributions, suggesting a potential difference in genomic architecture. The presence of the individual-sex draft genomes herein provides new clues regarding the origination and evolution of the diverse sex-determining mechanisms observed within Diptera. Additional genomic analysis of sex chromosomes and sex-determining genes of other blow flies will allow a refined evolutionary understanding of how flies with a typical X/Y heterogametic amphogeny (male and female offspring in similar ratios) sex determination systems evolved into one with a dominant factor that results in single sex progeny in a chromosomally monomorphic system.

Geographical distribution and molecular insights into abamectin and milbemectin cross-resistance in European field populations of Tetranychus urticae

BACKGROUND

Milbemectin and abamectin are frequently used to control the spider mite Tetranychus urticae. The development of abamectin resistance in this major pest has become an increasing problem worldwide, potentially compromising the use of milbemectin. In this study, a large collection of European field populations was screened for milbemectin and abamectin resistance, and both target-site and metabolic (cross-)resistance mechanisms were investigated.

RESULTS

High to very high levels of abamectin resistance were found in one third of all populations, while milbemectin resistance levels were low for most populations. The occurrence of well-known target-site resistance mutations in glutamate-gated chloride channels (G314D in GluCl1 and G326E in GluCl3) was documented in the most resistant populations. However, a new mutation, I321T in GluCl3, was also uncovered in three resistant populations, while a V327G and L329F mutation was found in GluCl3 of one resistant population. A differential gene-expression analysis revealed the overexpression of detoxification genes, more specifically cytochrome P450 monooxygenase (P450) and UDP-glycosyltransferase (UGT) genes. Multiple UGTs were functionally expressed, and their capability to glycosylate abamectin and milbemectin, was tested and confirmed.

CONCLUSIONS

We found a clear correlation between abamectin and milbemectin resistance in European T. urticae populations, but as milbemectin resistance levels were low, the observed cross-resistance is probably not of operational importance. The presence of target-site resistance mutations in GluCl genes was confirmed in most but not all resistant populations. Gene-expression analysis and functional characterization of P450s and UGTs suggests that also metabolic abamectin resistance mechanisms are common in European T. urticae populations. © 2020 Society of Chemical Industry.

The Identification and Evolutionary Trends of the Solute Carrier Superfamily in Arthropods

The solute carrier (SLC) transporter superfamily comprises an ancient and ubiquitous group of proteins capable of translocating a range of nutrients, endogenous molecules, and xenobiotics. Although the group has been the subject of intense investigation in both bacteria and mammals, its systematic identification in arthropods has not yet been undertaken. Here, we present a genome-wide identification of all 66 human SLC families in 174 arthropod species. A pipeline (SLC_id) was constructed to identify and group SLCs using a combination of hidden Markov model and BLAST searches followed by filtering based on polypeptide length and the number of transmembrane domains. Comparative analysis of the number of transporters in each family across diverse arthropod lineages was accomplished using one-way analysis of variance (ANOVA) and the Computational Analysis of gene Family Evolution (CAFE). These results suggested that many SLC families have undergone expansions or contractions in particular evolutionary lineages. Notably, the sugar transporting SLC2 family was significantly larger in insects compared with arachnids. This difference may have been complemented by a rapid expansion of the SLC60 family in arachnids which also acts on dietary sugars. Furthermore, the SLC33 family underwent a recent and drastic expansion in aphids, although the biological relevance of this expansion was not possible to infer. Information on specific SLC transporter families across arthropod species can be accessed through an R shiny web application at http://chrysalida.imbb.forth.gr : 3838/Arthropod_SLC_Database/. The present study greatly facilitates further investigation of the diverse group of SLC transporters in arthropods.

FlyBase 2.0: the next generation

FlyBase (flybase.org) is a knowledge base that supports the community of researchers that use the fruit fly, Drosophila melanogaster, as a model organism. The FlyBase team curates and organizes a diverse array of genetic, molecular, genomic, and developmental information about Drosophila. At the beginning of 2018, ‘FlyBase 2.0’ was released with a significantly improved user interface and new tools. Among these important changes are a new organization of search results into interactive lists or tables (hitlists), enhanced reference lists, and new protein domain graphics. An important new data class called ‘experimental tools’ consolidates information on useful fly strains and other resources related to a specific gene, which significantly enhances the ability of the Drosophila researcher to design and carry out experiments. With the release of FlyBase 2.0, there has also been a restructuring of backend architecture and a continued development of application programming interfaces (APIs) for programmatic access to FlyBase data. In this review, we describe these major new features and functionalities of the FlyBase 2.0 site and how they support the use of Drosophila as a model organism for biological discovery and translational research.

Behavioral and Transcriptional Response to Selection for Olfactory Behavior in Drosophila

The detection, discrimination, and behavioral responses to chemical cues in the environment can have marked effects on organismal survival and reproduction, eliciting attractive or aversive behavior. To gain insight into mechanisms mediating this hedonic valence, we applied thirty generations of divergent artificial selection for Drosophila melanogaster olfactory behavior. We independently selected for positive and negative behavioral responses to two ecologically relevant chemical compounds: 2,3-butanedione and cyclohexanone. We also tested the correlated responses to selection by testing behavioral responses to other odorants and life history traits. Measurements of behavioral responses of the selected lines and unselected controls to additional odorants showed that the mechanisms underlying responses to these odorants are, in some cases, differentially affected by selection regime and generalization of the response to other odorants was only detected in the 2,3-butanedione selection lines. Food consumption and lifespan varied with selection regime and, at times, sex. An analysis of gene expression of both selection regimes identified multiple differentially expressed genes. New genes and genes previously identified in mediating olfactory behavior were identified. In particular, we found functional enrichment of several gene ontology terms, including cell-cell adhesion and sulfur compound metabolic process, the latter including genes belonging to the glutathione S-transferase family. These findings highlight a potential role for glutathione S-transferases in the evolution of hedonic valence to ecologically relevant volatile compounds and set the stage for a detailed investigation into mechanisms by which these genes mediate attraction and aversion.

Macdonald S. Dissecting the Genetic Basis of Variation in Drosophila Sleep Using a Multiparental QTL Mapping Resource

There is considerable variation in sleep duration, timing and quality in human populations, and sleep dysregulation has been implicated as a risk factor for a range of health problems. Human sleep traits are known to be regulated by genetic factors, but also by an array of environmental and social factors. These uncontrolled, non-genetic effects complicate powerful identification of the loci contributing to sleep directly in humans. The model system, Drosophila melanogaster, exhibits a behavior that shows the hallmarks of mammalian sleep, and here we use a multitiered approach, encompassing high-resolution QTL mapping, expression QTL data, and functional validation with RNAi to investigate the genetic basis of sleep under highly controlled environmental conditions. We measured a battery of sleep phenotypes in >750 genotypes derived from a multiparental mapping panel and identified several, modest-effect QTL contributing to natural variation for sleep. Merging sleep QTL data with a large head transcriptome eQTL mapping dataset from the same population allowed us to refine the list of plausible candidate causative sleep loci. This set includes genes with previously characterized effects on sleep and circadian rhythms, in addition to novel candidates. Finally, we employed adult, nervous system-specific RNAi on the Dopa decarboxylase, dyschronic, and timeless genes, finding significant effects on sleep phenotypes for all three. The genes we resolve are strong candidates to harbor causative, regulatory variation contributing to sleep.

Detoxification Genes Differ Between Cactus-, Fruit-, and Flower-Feeding Drosophila

We use annotated genomes of 14 Drosophila species covering diverse host use phenotypes to test whether 4 gene families that often have detoxification functions are associated with host shifts among species. Bark, slime flux, flower, and generalist necrotic fruit-feeding species all have similar numbers of carboxyl/cholinesterase, glutathione S-transferase, cytochrome P450, and UDP-glucuronosyltransferase genes. However, species feeding on toxic Morinda citrifolia fruit and the fresh fruit-feeding Drosophila suzukii have about 30 and 60 more, respectively. ABC transporters show a different pattern, with the flower-feeding D. elegans and the generalist necrotic fruit and cactus feeder D. hydei having about 20 and >100 more than the other species, respectively. Surprisingly, despite the complex secondary chemistry we find that 3 cactophilic specialists in the mojavensis species cluster have variably fewer genes than any of the other species across all 4 families. We also find 82 positive selection events across the 4 families, with the terminal D. suzukii and M. citrifolia-feeding D. sechellia branches again having the highest number of such events in proportion to their respective branch lengths. Many of the genes involved in these host-use-specific gene number differences or positive selection events lie in specific clades of the gene families that have been recurrently associated with detoxification. Several genes are also found to be involved in multiple duplication and/or positive selection events across the species studied regardless of their host use phenotypes; the most frequently involved are the ABC transporter CG1718, which is not in a specific clade associated with detoxification, and the α-esterase gene cluster, which is.

A transcriptomic and proteomic atlas of expression in the Nezara viridula (Heteroptera: Pentatomidae) midgut suggests the compartmentalization of xenobiotic metabolism and nutrient digestion

BACKGROUND

Stink bugs are an emerging threat to crop security in many parts of the globe, but there are few genetic resources available to study their physiology at a molecular level. This is especially true for tissues such as the midgut, which forms the barrier between ingested material and the inside of the body.

RESULTS

Here, we focus on the midgut of the southern green stink bug Nezara viridula and use both transcriptomic and proteomic approaches to create an atlas of expression along the four compartments of the anterior-posterior axis. Estimates of the transcriptome completeness were high, which led us to compare our predicted gene set to other related stink bugs and Hemiptera, finding a high number of species-specific genes in N. viridula. To understand midgut function, gene ontology and gene family enrichment analyses were performed for the most highly expressed and specific genes in each midgut compartment. These data suggested a role for the anterior midgut (regions M1-M3) in digestion and xenobiotic metabolism, while the most posterior compartment (M4) was enriched in transmembrane proteins. A more detailed characterization of these findings was undertaken by identifying individual members of the cytochrome P450 superfamily and nutrient transporters thought to absorb amino acids or sugars.

CONCLUSIONS

These findings represent an initial step to understand the compartmentalization and physiology of the N. viridula midgut at a genetic level. Future studies will be able to build on this work and explore the molecular physiology of the stink bug midgut.

Building a pipeline to solicit expert knowledge from the community to aid gene summary curation

Brief summaries describing the function of each gene's product(s) are of great value to the research community, especially when interpreting genome-wide studies that reveal changes to hundreds of genes. However, manually writing such summaries, even for a single species, is a daunting task; for example, the Drosophila melanogaster genome contains almost 14 000 protein-coding genes. One solution is to use computational methods to generate summaries, but this often fails to capture the key functions or express them eloquently. Here, we describe how we solicited help from the research community to generate manually written summaries of D. melanogaster gene function. Based on the data within the FlyBase database, we developed a computational pipeline to identify researchers who have worked extensively on each gene. We e-mailed these researchers to ask them to draft a brief summary of the main function(s) of the gene's product, which we edited for consistency to produce a 'gene snapshot'. This approach yielded 1800 gene snapshot submissions within a 3-month period. We discuss the general utility of this strategy for other databases that capture data from the research literature. Database URL: https://flybase.org/.

Microarray-Based Quality Assessment as a Supporting Criterion for de novo Transcriptome Assembly Selection

RNA-Sequencing and de novo assembly have enabled the analysis of species with non-available reference transcriptomes, although intrinsic features (biological and technical) induce errors in the reconstruction. A strategy to resolve these errors consists of varying assembling process parameters to generate multiple reconstructions. However, the best assembly selection remains a challenge. Quantitative metrics for quality assessment have been inconsistent when compared with pertinent references. In this paper, a criterion for supporting assembly selection based on mapping DNA microarray hybridized probes to assembly sets is proposed. Mouse and fruit fly RNA-Seq datasets were assembled with standard de novo procedures. Quality assessment was estimated using quantitative metrics and the proposed criterion. The assembly that best mapped to the available reference transcriptomes of these model species provided the highest quality assembly. The hybridized probes identified the best assemblies, whereas quantitative metrics remained inconsistent. For example, subtle probe mapping difference of 0.25 percent, but statistically significant (ANOVA, p < 0.05), enabled the assembly selection that led to identify 3,719 more contigs and led to 1,049 further mapped contigs to the mouse reference transcriptome. The microarray data availability for non-model species makes the proposed criterion suitable for quality assessment of multiple de novo assembly strategies.