FlyBase at 25: looking to the future

Tissue-specific modulation of gene expression in response to lowered insulin signalling in Drosophila

Reduced activity of the insulin/IGF signalling network increases health during ageing in multiple species. Diverse and tissue-specific mechanisms drive the health improvement. Here, we performed tissue-specific transcriptional and proteomic profiling of long-lived Drosophila dilp2-3,5 mutants, and identified tissue-specific regulation of >3600 transcripts and >3700 proteins. Most expression changes were regulated post-transcriptionally in the fat body, and only in mutants infected with the endosymbiotic bacteria, Wolbachia pipientis, which increases their lifespan. Bioinformatic analysis identified reduced co-translational ER targeting of secreted and membrane-associated proteins and increased DNA damage/repair response proteins. Accordingly, age-related DNA damage and genome instability were lower in fat body of the mutant, and overexpression of a minichromosome maintenance protein subunit extended lifespan. Proteins involved in carbohydrate metabolism showed altered expression in the mutant intestine, and gut-specific overexpression of a lysosomal mannosidase increased autophagy, gut homeostasis, and lifespan. These processes are candidates for combatting ageing-related decline in other organisms.

Rapid divergence of the male reproductive proteins in the Drosophila dunni group and implications for postmating incompatibilities between species

Proteins involved in post-copulatory interactions between males and females are among the fastest evolving genes in many species, usually attributed to their involvement in reproductive conflict. As a result, these proteins are thought to often be involved in the formation of postmating-prezygotic incompatibilities between species. The Drosophila dunni subgroup consists of a dozen recently diverged species found across the Caribbean islands with varying levels of hybrid incompatibility. We performed experimental crosses between species in the dunni group and see some evidence of hybrid incompatibilities. We also find evidence of reduced survival following hybrid mating, likely due to postmating-prezygotic incompatibilities. We assessed rates of evolution between these species genomes and find evidence of rapid evolution and divergence of some reproductive proteins, specifically the seminal fluid proteins. This work suggests the rapid evolution of seminal fluid proteins may be associated with postmating-prezygotic isolation, which acts as a barrier for gene flow between even the most closely related species.

Genetic Mapping of a new Hippo allele, HpoN.1.2, in Drosophila melanogaster

Genetic screens provide a mechanism to identify genes involved with different cellular and organismal processes. Using a Flp/FRT screen in the Drosophila eye we identified mutations that result in alterations and de-regulation of cell growth and division. From this screen a group of undergraduate researchers part of the Fly-CURE consortium mapped and characterized a new allele of the gene Hippo, Hpo^N.1.2.

Adaptation of codon and amino acid use for translational functions in highly expressed cricket genes

BACKGROUND

For multicellular organisms, much remains unknown about the dynamics of synonymous codon and amino acid use in highly expressed genes, including whether their use varies with expression in different tissue types and sexes. Moreover, specific codons and amino acids may have translational functions in highly transcribed genes, that largely depend on their relationships to tRNA gene copies in the genome. However, these relationships and putative functions are poorly understood, particularly in multicellular systems.

RESULTS

Here, we studied codon and amino acid use in highly expressed genes from reproductive and nervous system tissues (male and female gonad, somatic reproductive system, brain and ventral nerve cord, and male accessory glands) in the cricket Gryllus bimaculatus. We report an optimal codon, defined as the codon preferentially used in highly expressed genes, for each of the 18 amino acids with synonymous codons in this organism. The optimal codons were mostly shared among tissue types and both sexes. However, the frequency of optimal codons was highest in gonadal genes. Concordant with translational selection, a majority of the optimal codons had abundant matching tRNA gene copies in the genome, but sometimes obligately required wobble tRNAs. We suggest the latter may comprise a mechanism for slowing translation of abundant transcripts, particularly for cell-cycle genes. Non-optimal codons, defined as those least commonly used in highly transcribed genes, intriguingly often had abundant tRNAs, and had elevated use in a subset of genes with specialized functions (gametic and apoptosis genes), suggesting their use promotes the translational upregulation of particular mRNAs. In terms of amino acids, we found evidence suggesting that amino acid frequency, tRNA gene copy number, and amino acid biosynthetic costs (size/complexity) had all interdependently evolved in this insect model, potentially for translational optimization.

CONCLUSIONS

Collectively, the results suggest a model whereby codon use in highly expressed genes, including optimal, wobble, and non-optimal codons, and their tRNA abundances, as well as amino acid use, have been influenced by adaptation for various functional roles in translation within this cricket. The effects of expression in different tissue types and the two sexes are discussed.

Glycosyltransferases: the multifaceted enzymatic regulator in insects

Glycosyltransferases (GTs) catalyse the reaction of glyco-conjugation of various biomolecules by transferring the saccharide moieties from an activated nucleotide sugar to nucleophilic glycosyl acceptor. In insects, GTs show diverse temporal and site-specific expression patterns and thus play significant roles in forming the complex biomolecular structures that are necessary for insect survival, growth and development. Several insects exhibit GT-mediated detoxification as a key defence strategy against plant allelochemicals and xenobiotic compounds, as well as a mechanism for pesticide cross-resistance. Also, these enzymes act as crucial effectors and modulators in various developmental processes of insects such as eye development, UV shielding, cuticle formation, epithelial development and other specialized functions. Furthermore, many of the known insect GTs have been shown to play a fundamental role in other physiological processes like body pigmentation, cuticular tanning, chemosensation and stress response. This review provides a detailed overview of the multifaceted functionality of insect GTs and summarizes numerous case studies associated with it.

The nuclear to cytoplasmic ratio directly regulates zygotic transcription in Drosophila through multiple modalities

Early embryos must rapidly generate large numbers of cells to form an organism. Many species accomplish this through a series of rapid, reductive, and transcriptionally silent cleavage divisions. Previous work has demonstrated that the number of divisions before both cell cycle elongation and zygotic genome activation (ZGA) is regulated by the ratio of nuclear content to cytoplasm (N/C). To understand how the N/C ratio affects the timing of ZGA, we directly assayed the behavior of several previously identified N/C ratio-dependent genes using the MS2-MCP reporter system in living Drosophila embryos with altered ploidy and cell cycle durations. For every gene that we examined, we found that nascent RNA output per cycle is delayed in haploid embryos. Moreover, we found that the N/C ratio influences transcription through three overlapping modes of action. For some genes (knirps, fushi tarazu, and snail), the effect of ploidy can be primarily attributed to changes in cell cycle duration. However, additional N/C ratio-mediated mechanisms contribute significantly to transcription delays for other genes. For giant and bottleneck, the kinetics of transcription activation are significantly disrupted in haploids, while for frühstart and Krüppel, the N/C ratio controls the probability of transcription initiation. Our data demonstrate that the regulatory elements of N/C ratio-dependent genes respond directly to the N/C ratio through multiple modes of regulation.

Ultra-conserved sequences in the genomes of highly diverse Anopheles mosquitoes, with implications for malaria vector control

DNA sequences that are exactly conserved over long evolutionary time scales have been observed in a variety of taxa. Such sequences are likely under strong functional constraint and they have been useful in the field of comparative genomics for identifying genome regions with regulatory function. A potential new application for these ultra-conserved elements (UCEs) has emerged in the development of gene drives to control mosquito populations. Many gene drives work by recognizing and inserting at a specific target sequence in the genome, often imposing a reproductive load as a consequence. They can therefore select for target sequence variants that provide resistance to the drive. Focusing on highly conserved, highly constrained sequences lowers the probability that variant, gene drive-resistant alleles can be tolerated. Here, we search for conserved sequences of 18 bp and over in an alignment of 21 Anopheles genomes, spanning an evolutionary timescale of 100 million years, and characterize the resulting sequences according to their location and function. Over 8000 UCEs were found across the alignment, with a maximum length of 164 bp. Length-corrected gene ontology analysis revealed that genes containing Anopheles UCEs were over-represented in categories with structural or nucleotide-binding functions. Known insect transcription factor binding sites were found in 48% of intergenic Anopheles UCEs. When we looked at the genome sequences of 1142 wild-caught mosquitoes, we found that 15% of the Anopheles UCEs contained no polymorphisms. Our list of Anopheles UCEs should provide a valuable starting point for the selection and testing of new targets for gene-drive modification in the mosquitoes that transmit malaria.

GWAS reveal a role for the central nervous system in regulating weight and weight change in response to exercise

Body size and weight show considerable variation both within and between species. This variation is controlled in part by genetics, but also strongly influenced by environmental factors including diet and the level of activity experienced by the individual. Due to the increasing obesity epidemic in much of the world, there is considerable interest in the genetic factors that control body weight and how weight changes in response to exercise treatments. Here, we address this question in the Drosophila model system, utilizing 38 strains of the Drosophila Genetics Reference Panel. We use GWAS to identify the molecular pathways that control weight and weight changes in response to exercise. We find that there is a complex set of molecular pathways controlling weight, with many genes linked to the central nervous system (CNS). The CNS also plays a role in the weight change with exercise, in particular, signaling from the CNS. Additional analyses revealed that weight in Drosophila is driven by two factors, animal size, and body composition, as the amount of fat mass versus lean mass impacts the density. Thus, while the CNS appears to be important for weight and exercise-induced weight change, signaling pathways are particularly important for determining how exercise impacts weight.

Biological impact of mutually exclusive exon switching

Alternative splicing can expand the diversity of proteomes. Homologous mutually exclusive exons (MXEs) originate from the same ancestral exon and result in polypeptides with similar structural properties but altered sequence. Why would some genes switch homologous exons and what are their biological impact? Here, we analyse the extent of sequence, structural and functional variability in MXEs and report the first large scale, structure-based analysis of the biological impact of MXE events from different genomes. MXE-specific residues tend to map to single domains, are highly enriched in surface exposed residues and cluster at or near protein functional sites. Thus, MXE events are likely to maintain the protein fold, but alter specificity and selectivity of protein function. This comprehensive resource of MXE events and their annotations is available at: http://gene3d.biochem.ucl.ac.uk/mxemod/. These findings highlight how small, but significant changes at critical positions on a protein surface are exploited in evolution to alter function.

Control of Cell Growth and Proliferation by the Tribbles Pseudokinase: Lessons from Drosophila

Simple Summary

Tribbles pseudokinases represent a sub-branch of the CAMK (Ca²⁺/calmodulin-dependent protein kinase) subfamily and are associated with disease-associated signaling pathways associated with various cancers, including melanoma, lung, liver, and acute leukemia. The ability of this class of molecules to regulate cell proliferation was first recognized in the model organism Drosophila and the fruit fly genetic model and continues to provide insight into the molecular mechanism by which this family of adapter molecules regulates both normal development and disease associated with corruption of their proper regulation and function.

Abstract

The Tribbles (Trib) family of pseudokinase proteins regulate cell growth, proliferation, and differentiation during normal development and in response to environmental stress. Mutations in human Trib isoforms (Trib1, 2, and 3) have been associated with metabolic disease and linked to leukemia and the formation of solid tumors, including melanomas, hepatomas, and lung cancers. Drosophila Tribbles (Trbl) was the first identified member of this sub-family of pseudokinases and shares a conserved structure and similar functions to bind and direct the degradation of key mediators of cell growth and proliferation. Common Trib targets include Akt kinase (also known as protein kinase B), C/EBP (CAAT/enhancer binding protein) transcription factors, and Cdc25 phosphatases, leading to the notion that Trib family members stand athwart multiple pathways modulating their growth-promoting activities. Recent work using the Drosophila model has provided important insights into novel facets of conserved Tribbles functions in stem cell quiescence, tissue regeneration, metabolism connected to insulin signaling, and tumor formation linked to the Hippo signaling pathway. Here we highlight some of these recent studies and discuss their implications for understanding the complex roles Tribs play in cancers and disease pathologies.

B.2.16 is a non-lethal modifier of the Dark82 mosaic eye phenotype in Drosophila melanogaster

Genetic screens have been used to identify genes involved in the regulation of different biological processes. We identified growth mutants in a Flp/FRT screen using the Drosophila melanogaster eye to identify conditional regulators of cell growth and cell division. One mutant identified from this screen, B.2.16, was mapped and characterized by researchers in undergraduate genetics labs as part of the Fly-CURE. We find that B.2.16 is a non-lethal genetic modifier of the Dark⁸² mosaic eye phenotype.

Adaptation, ancestral variation and gene flow in a 'Sky Island' Drosophila species

Over time, populations of species can expand, contract, fragment and become isolated, creating subpopulations that must adapt to local conditions. Understanding how species maintain variation after divergence as well as adapt to these changes in the face of gene flow is of great interest, especially as the current climate crisis has caused range shifts and frequent migrations for many species. Here, we characterize how a mycophageous fly species, Drosophila innubila, came to inhabit and adapt to its current range which includes mountain forests in south-western USA separated by large expanses of desert. Using population genomic data from more than 300 wild-caught individuals, we examine four populations to determine their population history in these mountain forests, looking for signatures of local adaptation. In this first extensive study, establishing D. innubila as a key genomic "Sky Island" model, we find D. innubila spread northwards during the previous glaciation period (30-100 KYA) and have recently expanded even further (0.2-2 KYA). D. innubila shows little evidence of population structure, consistent with a recent establishment and genetic variation maintained since before geographic stratification. We also find some signatures of recent selective sweeps in chorion proteins and population differentiation in antifungal immune genes suggesting differences in the environments to which flies are adapting. However, we find little support for long-term recurrent selection in these genes. In contrast, we find evidence of long-term recurrent positive selection in immune pathways such as the Toll signalling system and the Toll-regulated antimicrobial peptides.

Dis3L2 regulates cell proliferation and tissue growth through a conserved mechanism

Dis3L2 is a highly conserved 3’-5’ exoribonuclease which is mutated in the human overgrowth disorders Perlman syndrome and Wilms’ tumour of the kidney. Using Drosophila melanogaster as a model system, we have generated a new dis3L2 null mutant together with wild-type and nuclease-dead genetic lines in Drosophila to demonstrate that the catalytic activity of Dis3L2 is required to control cell proliferation. To understand the cellular pathways regulated by Dis3L2 to control proliferation, we used RNA-seq on dis3L2 mutant wing discs to show that the imaginal disc growth factor Idgf2 is responsible for driving the wing overgrowth. IDGFs are conserved proteins homologous to human chitinase-like proteins such as CHI3L1/YKL-40 which are implicated in tissue regeneration as well as cancers including colon cancer and non-small cell lung cancer. We also demonstrate that loss of DIS3L2 in human kidney HEK-293T cells results in cell proliferation, illustrating the conservation of this important cell proliferation pathway. Using these human cells, we show that loss of DIS3L2 results in an increase in the PI3-Kinase/AKT signalling pathway, which we subsequently show to contribute towards the proliferation phenotype in Drosophila. Our work therefore provides the first mechanistic explanation for DIS3L2-induced overgrowth in humans and flies and identifies an ancient proliferation pathway controlled by Dis3L2 to regulate cell proliferation and tissue growth.

Regulation of cell proliferation is not only important during development but also required for repair of damaged tissues and during wound healing. Using human kidney cells as well as the fruit fly Drosophila we have recently discovered that cell proliferation can be regulated by a protein named Dis3L2. Depletion or removal of this protein results in excess proliferation. These results are relevant to human disease as DIS3L2 has been shown to be mutated in an overgrowth syndrome (Perlman syndrome) where affected children have abnormal enlargement of organs (e.g. kidneys) and susceptibility to Wilms’ tumour (a kidney cancer). Dis3L2 is an enzyme known to "chew up" mRNA molecules which instruct the cell to make particular proteins. Using state-of-the-art molecular methods in Drosophila, we have discovered that Dis3L2 targets a small subset of mRNAs, including an mRNA encoding a growth factor named 'imaginal disc growth factor 2' (idgf2). For human kidney cells in culture, we have found that depletion of DIS3L2 results in enhanced proliferation, and that this involves a well-known cellular pathway. Our results mean that we have discovered a new way of controlling cell proliferation, which could, in the future, be used in human therapies.

Crosstalk between vrille transcripts, proteins, and regulatory elements controlling circadian rhythms and development in Drosophila

The vrille (vri) gene encodes a transcriptional repressor required for Drosophila development as well as circadian behavior in adults. Alternate first exons produce vri transcripts predicted to produce a short VRI isoform during development and long VRI in adults. A vri mutant (vri^Δ679) lacking long VRI transcripts is viable, confirming that short VRI is sufficient for developmental functions, yet behavioral rhythms in vri^Δ679 flies persist, showing that short VRI is sufficient for clock output. E-box regulatory elements that drive rhythmic long VRI transcript expression are required for developmental expression of short VRI transcripts. Surprisingly, long VRI transcripts primarily produce short VRI in adults, apparently due to a poor Kozak sequence context, demonstrating that short VRI drives circadian behavior. Thus, E-box-driven long VRI transcripts primarily control circadian rhythms via short VRI, whereas the same E-boxes drive short VRI transcripts that control developmental functions using short VRI.

•

vri-E mRNA is sufficient for Drosophila development and circadian behavior

•

E-boxes upstream of the vri-ADF promoter are required for Drosophila development

•

vri-ADF mRNAs primarily produce short VRI protein rather than long VRI protein

•

Short VRI protein primarily controls Drosophila development and circadian behavior

Leveraging Multiple Layers of Data To Predict Drosophila Complex Traits

The ability to accurately predict complex trait phenotypes from genetic and genomic data are critical for the implementation of personalized medicine and precision agriculture; however, prediction accuracy for most complex traits is currently low. Here, we used data on whole genome sequences, deep RNA sequencing, and high quality phenotypes for three quantitative traits in the ∼200 inbred lines of the Drosophila melanogaster Genetic Reference Panel (DGRP) to compare the prediction accuracies of gene expression and genotypes for three complex traits. We found that expression levels (r = 0.28 and 0.38, for females and males, respectively) provided higher prediction accuracy than genotypes (r = 0.07 and 0.15, for females and males, respectively) for starvation resistance, similar prediction accuracy for chill coma recovery (null for both models and sexes), and lower prediction accuracy for startle response (r = 0.15 and 0.14 for female and male genotypes, respectively; and r = 0.12 and 0.11, for females and male transcripts, respectively). Models including both genotype and expression levels did not outperform the best single component model. However, accuracy increased considerably for all the three traits when we included gene ontology (GO) category as an additional layer of information for both genomic variants and transcripts. We found strongly predictive GO terms for each of the three traits, some of which had a clear plausible biological interpretation. For example, for starvation resistance in females, GO:0033500 (r = 0.39 for transcripts) and GO:0032870 (r = 0.40 for transcripts), have been implicated in carbohydrate homeostasis and cellular response to hormone stimulus (including the insulin receptor signaling pathway), respectively. In summary, this study shows that integrating different sources of information improved prediction accuracy and helped elucidate the genetic architecture of three Drosophila complex phenotypes.

Reproductive consequences of an extra long-term sperm storage organ

Background

Sperm storage plays a key role in the reproductive success of many sexually-reproducing organisms, and the capacity of long-term sperm storage varies across species. While there are theoretical explanations for why such variation exists, to date there are no controlled empirical tests of the reproductive consequences of additional long-term sperm storage. While Dipterans ancestrally have three long-term sperm organs, known as the spermathecae, Drosophila contain only two.

Results

We identified a candidate gene, which we call spermathreecae (sp3), in which a disruption cause the development of three functional spermathecae rather than the usual two in Drosophila. We used this disruption to test the reproductive consequences of having an additional long-term sperm storage organ. Compared to females with two spermathecae, females with three spermathecae store a greater total number of sperm and can produce offspring a greater length of time. However, they did not produce a greater total number of offspring.

Conclusions

Thus, additional long-term sperm storage in insects may increase female fitness through extending the range of conditions where she produces offspring, or through increasing the quality of offspring via enhanced local sperm competition at fertilization.

Genes underlying species differences in cuticular hydrocarbon production between Drosophila melanogaster and D. simulans

Surface chemical compounds are key components of survival and reproduction in many species. Cuticular hydrocarbons (CHCs) are chemical compounds produced by all insects that are used for both desiccation resistance and chemical communication, including communication related to mating. In the species pair of Drosophila melanogaster and D. simulans, female CHCs stimulate conspecific males to mate and repel heterospecific males. While CHCs are a critical contributor to both reproductive success within a species and isolation between species, few genes underlying species variation in CHC profiles are known. Here, we use genetic mapping of the 3rd chromosome to test a suite of candidate genes for interspecies variation in CHCs. Candidate gene CG5946 was found to be involved in species differences in the production of 7,11-heptacosadiene and 7-tricosene between D. melanogaster and D. simulans. This is therefore a new candidate locus contributing to species-specific variation in the CHC profile. In the process of mapping genes for CHCs, we also identified 29 candidate genes for the reduced survival or inviability of interspecies hybrids.

A Burst of Genetic Innovation in Drosophila Actin-Related Proteins for Testis-Specific Function

Many cytoskeletal proteins perform fundamental biological processes and are evolutionarily ancient. For example, the superfamily of actin-related proteins (Arps) specialized early in eukaryotic evolution for diverse cellular roles in the cytoplasm and the nucleus. Despite its strict conservation across eukaryotes, we find that the Arp superfamily has undergone dramatic lineage-specific diversification in Drosophila. Our phylogenomic analyses reveal four independent Arp gene duplications that occurred in the common ancestor of the obscura group of Drosophila and have been mostly preserved in this lineage. All four obscura-specific Arp paralogs are predominantly expressed in the male germline and have evolved under positive selection. We focus our analyses on the divergent Arp2D paralog, which arose via a retroduplication event from Arp2, a component of the Arp2/3 complex that polymerizes branched actin networks. Computational modeling analyses suggest that Arp2D can replace Arp2 in the Arp2/3 complex and bind actin monomers. Together with the signature of positive selection, our findings suggest that Arp2D may augment Arp2's functions in the male germline. Indeed, we find that Arp2D is expressed during and following male meiosis, where it localizes to distinct locations such as actin cones-specialized cytoskeletal structures that separate bundled spermatids into individual mature sperm. We hypothesize that this unprecedented burst of genetic innovation in cytoskeletal proteins may have been driven by the evolution of sperm heteromorphism in the obscura group of Drosophila.

A phylogenomic study of Steganinae fruit flies (Diptera: Drosophilidae): strong gene tree heterogeneity and evidence for monophyly

BACKGROUND

The Drosophilidae family is traditionally divided into two subfamilies: Drosophilinae and Steganinae. This division is based on morphological characters, and the two subfamilies have been treated as monophyletic in most of the literature, but some molecular phylogenies have suggested Steganinae to be paraphyletic. To test the paraphyletic-Steganinae hypothesis, here, we used genomic sequences of eight Drosophilidae (three Steganinae and five Drosophilinae) and two Ephydridae (outgroup) species and inferred the phylogeny for the group based on a dataset of 1,028 orthologous genes present in all species (> 1,000,000 bp). This dataset includes three genera that broke the monophyly of the subfamilies in previous works. To investigate possible biases introduced by small sample sizes and automatic gene annotation, we used the same methods to infer species trees from a set of 10 manually annotated genes that are commonly used in phylogenetics.

RESULTS

Most of the 1,028 gene trees depicted Steganinae as paraphyletic with distinct topologies, but the most common topology depicted it as monophyletic (43.7% of the gene trees). Despite the high levels of gene tree heterogeneity observed, species tree inference in ASTRAL, in PhyloNet, and with the concatenation approach strongly supported the monophyly of both subfamilies for the 1,028-gene dataset. However, when using the concatenation approach to infer a species tree from the smaller set of 10 genes, we recovered Steganinae as a paraphyletic group. The pattern of gene tree heterogeneity was asymmetrical and thus could not be explained solely by incomplete lineage sorting (ILS).

CONCLUSIONS

Steganinae was clearly a monophyletic group in the dataset that we analyzed. In addition to ILS, gene tree discordance was possibly the result of introgression, suggesting complex branching processes during the early evolution of Drosophilidae with short speciation intervals and gene flow. Our study highlights the importance of genomic data in elucidating contentious phylogenetic relationships and suggests that phylogenetic inference for drosophilids based on small molecular datasets should be performed cautiously. Finally, we suggest an approach for the correction and cleaning of BUSCO-derived genomic datasets that will be useful to other researchers planning to use this tool for phylogenomic studies.

Transcriptomic changes across vitellogenesis in the black tiger prawn (Penaeus monodon), neuropeptides and G protein-coupled receptors repertoire curation

The black tiger prawn (Penaeus monodon) is one of the most commercially important prawn species world-wide, yet there are currently key issues that hinder aquaculture of this species, such as low spawning capacity of captive-reared broodstock females and lack of globally available fully domesticated strains. In this study, we analysed the molecular changes that occur from vitellogenesis to spawning of a fully domesticated population of P.monodon (Madagascar) using four tissues [brain and thoracic ganglia (central nervous system - CNS), eyestalks, antennal gland, and ovary] highlighting differentially expressed genes that could be involved in the sexual maturation. In addition, due to their key role in regulating multiple physiological processes including reproduction, transcripts encoding P.monodon neuropeptides and G protein-coupled receptors (GPCRs) were identified and their expression pattern was assessed. A few neuropeptides and their putative GPCRs which were previously implicated in reproduction are discussed. We identified 573 differentially expressed transcripts between previtellogenic and vitellogenic stages, across the four analysed tissues. Multiple transcripts that have been linked to ovarian maturation were highlighted throughout the study, these include vitellogenin, Wnt, heat shock protein 21, heat shock protein 90, teneurin, Fs(1)M3, hemolymph clottable proteins and some other candidates. Seventy neuropeptide transcripts were also characterized from our de novo assembly. In addition, a hybrid approach that involved clustering and phylogenetics analysis was used to annotate all P. monodon GPCRs, revealing 223 Rhodopsin, 100 Secretin and 27 Metabotropic glutamate GPCRs. Given the key commercial significance of P.monodon and the industry requirements for developing better genomic tools to control reproduction in this species, our findings provide a foundation for future gene-based studies, setting the scene for developing innovative tools for reproduction and/or sexual maturation control in P. monodon.

Interactions between the microbiome and mating influence the female's transcriptional profile in Drosophila melanogaster

Drosophila melanogaster females undergo a variety of post-mating changes that influence their activity, feeding behavior, metabolism, egg production and gene expression. These changes are induced either by mating itself or by sperm or seminal fluid proteins. In addition, studies have shown that axenic females-those lacking a microbiome-have altered fecundity compared to females with a microbiome, and that the microbiome of the female's mate can influence reproductive success. However, the extent to which post-mating changes in transcript abundance are affected by microbiome state is not well-characterized. Here we investigated fecundity and the post-mating transcript abundance profile of axenic or control females after mating with either axenic or control males. We observed interactions between the female's microbiome and her mating status: transcripts of genes involved in reproduction and genes with neuronal functions were differentially abundant depending on the females' microbiome status, but only in mated females. In addition, immunity genes showed varied responses to either the microbiome, mating, or a combination of those two factors. We further observed that the male's microbiome status influences the fecundity of both control and axenic females, while only influencing the transcriptional profile of axenic females. Our results indicate that the microbiome plays a vital role in the post-mating switch of the female's transcriptome.

Genome streamlining in a minute herbivore that manipulates its host plant

The tomato russet mite, Aculops lycopersici, is among the smallest animals on earth. It is a worldwide pest on tomato and can potently suppress the host's natural resistance. We sequenced its genome, the first of an eriophyoid, and explored whether there are genomic features associated with the mite's minute size and lifestyle. At only 32.5 Mb, the genome is the smallest yet reported for any arthropod and, reminiscent of microbial eukaryotes, exceptionally streamlined. It has few transposable elements, tiny intergenic regions, and is remarkably intron-poor, as more than 80% of coding genes are intronless. Furthermore, in accordance with ecological specialization theory, this defense-suppressing herbivore has extremely reduced environmental response gene families such as those involved in chemoreception and detoxification. Other losses associate with this species' highly derived body plan. Our findings accelerate the understanding of evolutionary forces underpinning metazoan life at the limits of small physical and genome size.

An early endosome-derived retrograde trafficking pathway promotes secretory granule maturation

Regulated secretion is a fundamental cellular process in which biologically active molecules stored in long-lasting secretory granules (SGs) are secreted in response to external stimuli. Many studies have described mechanisms responsible for biogenesis and secretion of SGs, but how SGs mature remains poorly understood. In a genetic screen, we discovered a large number of endolysosomal trafficking genes required for proper SG maturation, indicating that maturation of SGs might occur in a manner similar to lysosome-related organelles (LROs). CD63, a tetraspanin known to decorate LROs, also decorates SG membranes and facilitates SG maturation. Moreover, CD63-mediated SG maturation requires type II phosphatidylinositol 4 kinase (PI4KII)-dependent early endosomal sorting and accumulation of phosphatidylinositol 4-phosphate (PI4P) on SG membranes. In addition, the PI4P effector Past1 is needed for formation of stable PI4KII-containing endosomal tubules associated with this process. Our results reveal that maturation of post-Golgi-derived SGs requires trafficking via the endosomal system, similar to mechanisms employed by LROs.

Drosophila R8 photoreceptor cell subtype specification requires hibris

Cell differentiation and cell fate determination in sensory systems are essential for stimulus discrimination and coding of environmental stimuli. Color vision is based on the differential color sensitivity of retinal photoreceptors, however the developmental programs that control photoreceptor cell differentiation and specify color sensitivity are poorly understood. In Drosophila melanogaster, there is evidence that the color sensitivity of different photoreceptors in the compound eye is regulated by inductive signals between cells, but the exact nature of these signals and how they are propagated remains unknown. We conducted a genetic screen to identify additional regulators of this process and identified a novel mutation in the hibris gene, which encodes an irre cell recognition module protein (IRM). These immunoglobulin super family cell adhesion molecules include human KIRREL and nephrin (NPHS1). hibris is expressed dynamically in the developing Drosophila melanogaster eye and loss-of-function mutations give rise to a diverse range of mutant phenotypes including disruption of the specification of R8 photoreceptor cell diversity. We demonstrate that hibris is required within the retina, and that hibris over-expression is sufficient to disrupt normal photoreceptor cell patterning. These findings suggest an additional layer of complexity in the signaling process that produces paired expression of opsin genes in adjacent R7 and R8 photoreceptor cells.

Proteome-Wide Analysis of ADAR-Mediated Messenger RNA Editing during Fruit Fly Ontogeny

Adenosine-to-inosine RNA editing is an enzymatic post-transcriptional modification which modulates immunity and neural transmission in multicellular organisms. In particular, it involves editing of mRNA codons with the resulting amino acid substitutions. We identified such sites for developmental proteomes of Drosophila melanogaster at the protein level using available data for 15 stages of fruit fly development from egg to imago and 14 time points of embryogenesis. In total, 40 sites were obtained, each belonging to a unique protein, including four sites related to embryogenesis. The interactome analysis has revealed that the majority of the editing-recoded proteins were associated with synaptic vesicle trafficking and actomyosin organization. Quantitation data analysis suggested the existence of a phase-specific RNA editing regulation with yet unknown mechanisms. These findings supported the transcriptome analysis results, which showed that a burst in the RNA editing occurs during insect metamorphosis from pupa to imago. Finally, targeted proteomic analysis was performed to quantify editing-recoded and genomically encoded versions of five proteins in brains of larvae, pupae, and imago insects, which showed a clear tendency toward an increase in the editing rate for each of them. These results will allow a better understanding of the protein role in physiological effects of RNA editing.

A functionally defined high-density NRF2 interactome reveals new conditional regulators of ARE transactivation

NRF2 (NFE2L2) is a cytoprotective transcription factor associated with >60 human diseases, adverse drug reactions and therapeutic resistance. To provide insight into the complex regulation of NRF2 responses, 1962 predicted NRF2-partner interactions were systematically tested to generate an experimentally defined high-density human NRF2 interactome. Verification and conditional stratification of 46 new NRF2 partners was achieved by co-immunoprecipitation and the novel integration of quantitative data from dual luminescence-based co-immunoprecipitation (DULIP) assays and live-cell fluorescence cross-correlation spectroscopy (FCCS). The functional impact of new partners was then assessed in genetically edited loss-of-function (NRF2-/-) and disease-related gain-of-function (NRF2T80K and KEAP1-/-) cell-lines. Of the new partners investigated >77% (17/22) modified NRF2 responses, including partners that only exhibited effects under disease-related conditions. This experimentally defined binary NRF2 interactome provides a new vision of the complex molecular networks that govern the modulation and consequence of NRF2 activity in health and disease.

Functional Transcription Factor Target Networks Illuminate Control of Epithelial Remodelling

Cell identity is governed by gene expression, regulated by transcription factor (TF) binding at cis-regulatory modules. Decoding the relationship between TF binding patterns and gene regulation is nontrivial, remaining a fundamental limitation in understanding cell decision-making. We developed the NetNC software to predict functionally active regulation of TF targets; demonstrated on nine datasets for the TFs Snail, Twist, and modENCODE Highly Occupied Target (HOT) regions. Snail and Twist are canonical drivers of epithelial to mesenchymal transition (EMT), a cell programme important in development, tumour progression and fibrosis. Predicted "neutral" (non-functional) TF binding always accounted for the majority (50% to 95%) of candidate target genes from statistically significant peaks and HOT regions had higher functional binding than most of the Snail and Twist datasets examined. Our results illuminated conserved gene networks that control epithelial plasticity in development and disease. We identified new gene functions and network modules including crosstalk with notch signalling and regulation of chromatin organisation, evidencing networks that reshape Waddington's epigenetic landscape during epithelial remodelling. Expression of orthologous functional TF targets discriminated breast cancer molecular subtypes and predicted novel tumour biology, with implications for precision medicine. Predicted invasion roles were validated using a tractable cell model, supporting our approach.

Topology-driven protein-protein interaction network analysis detects genetic sub-networks regulating reproductive capacity

Understanding the genetic regulation of organ structure is a fundamental problem in developmental biology. Here, we use egg-producing structures of insect ovaries, called ovarioles, to deduce systems-level gene regulatory relationships from quantitative functional genetic analysis. We previously showed that Hippo signalling, a conserved regulator of animal organ size, regulates ovariole number in Drosophila melanogaster. To comprehensively determine how Hippo signalling interacts with other pathways in this regulation, we screened all known signalling pathway genes, and identified Hpo-dependent and Hpo-independent signalling requirements. Network analysis of known protein-protein interactions among screen results identified independent gene regulatory sub-networks regulating one or both of ovariole number and egg laying. These sub-networks predict involvement of previously uncharacterised genes with higher accuracy than the original candidate screen. This shows that network analysis combining functional genetic and large-scale interaction data can predict function of novel genes regulating development.

Validating the Potential of Double-Stranded RNA Targeting Colorado Potato Beetle Mesh Gene in Laboratory and Field Trials

Colorado potato beetle (CPB) is an agricultural pest of solanaceous crops, notorious for its rapid resistance development to chemical pesticides. Foliar spraying of dsRNA formulations is a promising innovative technology providing highly specific and environmentally acceptable option for CPB management. We designed dsRNA to silence CPB mesh gene (dsMESH) and performed laboratory feeding trials to assess impacts on beetle survival and development. We compared the effectiveness of in vivo and in vitro produced dsRNA in a series of laboratory experiments. We additionally performed a field trial in which the efficacy of dsRNA sprayed onto potato foliage was compared to a spinosad-based insecticide. We showed that dsMESH ingestion consistently and significantly impaired larval growth and decreased larval survival in laboratory feeding experiments. In vivo produced dsRNA performed similarly as in vitro synthesized dsRNA in laboratory settings. In the field trial, dsMESH was as effective in controlling CPB larvae as a commercial spinosad insecticide, its activity was however slower. We discuss limitations and benefits of a potential dsMESH-based CPB management strategy and list some important RNAi based CPB research topics, which will have to be addressed in future.

Out of the testis, into the ovary: biased outcomes of gene duplication and deletion in Drosophila

Gene turnover is a key source of adaptive variation. Yet most evolutionary studies have focused on gene duplication, dismissing gene deletion as a mechanism that simply eradicates redundancy. Here, I use genome-scale sequence and multi-tissue expression data from Drosophila melanogaster and Drosophila pseudoobscura to simultaneously assess the evolutionary outcomes of gene duplication and deletion in Drosophila. I find that gene duplication is more frequent than gene deletion in both species, indicating that it may play a more important role in Drosophila evolution. However, examination of several genic properties reveals that genes likely possess distinct functions after duplication that diverge further before deletion, suggesting that loss of redundancy cannot explain a majority of gene deletion events in Drosophila. Moreover, in addition to providing support for the well-known "out of the testis" origin of young duplicate genes, analyses of gene expression profiles uncover a preferential bias against deletion of old ovary-expressed genes. Therefore, I propose a novel "into the ovary" hypothesis for gene deletion in Drosophila, in which gene deletion may promote adaptation by salvaging genes that contribute to the evolution of female reproductive phenotypes. Under this combined "out of the testis, into the ovary" evolutionary model, gene duplication and deletion work in concert to generate and maintain a balanced repertoire of genes that promote sex-specific adaptation in Drosophila.

JNK pathway restricts DENV2, ZIKV and CHIKV infection by activating complement and apoptosis in mosquito salivary glands

Arbovirus infection of Aedes aegypti salivary glands (SGs) determines transmission. However, there is a dearth of knowledge on SG immunity. Here, we characterized SG immune response to dengue, Zika and chikungunya viruses using high-throughput transcriptomics. We also describe a transcriptomic response associated to apoptosis, blood-feeding and lipid metabolism. The three viruses differentially regulate components of Toll, Immune deficiency (IMD) and c-Jun N- terminal Kinase (JNK) pathways. However, silencing of the Toll and IMD pathway components showed variable effects on SG infection by each virus. In contrast, regulation of the JNK pathway produced consistent responses in both SGs and midgut. Infection by the three viruses increased with depletion of the activator Kayak and decreased with depletion of the negative regulator Puckered. Virus-induced JNK pathway regulates the complement factor, Thioester containing protein-20 (TEP20), and the apoptosis activator, Dronc, in SGs. Individual and co-silencing of these genes demonstrate their antiviral effects and that both may function together. Co-silencing either TEP20 or Dronc with Puckered annihilates JNK pathway antiviral effect. Upon infection in SGs, TEP20 induces antimicrobial peptides (AMPs), while Dronc is required for apoptosis independently of TEP20. In conclusion, we revealed the broad antiviral function of JNK pathway in SGs and showed that it is mediated by a TEP20 complement and Dronc-induced apoptosis response. These results expand our understanding of the immune arsenal that blocks arbovirus transmission.

Arboviral diseases caused by dengue (DENV), Zika (ZIKV) and chikungunya (CHIKV) viruses are responsible for large number of death and debilitation around the world. These viruses are transmitted to humans by the mosquito vector, Aedes aegypti. During the bites, infected salivary glands (SGs) release saliva containing viruses, which initiate human infection. As the tissue where transmitted viruses are produced, SG infection is a key determinant of transmission. To bridge the knowledge gap in vector-virus molecular interactions in SGs, we describe the transcriptome after DENV, ZIKV and CHIKV infection using RNA-sequencing and characterized the immune response in this tissue. Our study reveals the broad antiviral function of c-Jun N-terminal kinase (JNK) pathway against DENV, ZIKV and CHIKV in SGs. We further show that it is mediated by the complement system and apoptosis, identifying the mechanism. Our study adds the JNK pathway to the immune arsenal that can be harnessed to engineer refractory vectors.

SPOCD1 is an essential executor of piRNA-directed de novo DNA methylation

In mammals, the acquisition of the germline from the soma provides the germline with an essential challenge: the need to erase and reset genomic methylation1. In the male germline, RNA-directed DNA methylation silences young, active transposable elements2-4. The PIWI protein MIWI2 (PIWIL4) and its associated PIWI-interacting RNAs (piRNAs) instruct DNA methylation of transposable elements3,5. piRNAs are proposed to tether MIWI2 to nascent transposable element transcripts; however, the mechanism by which MIWI2 directs the de novo methylation of transposable elements is poorly understood, although central to the immortality of the germline. Here we define the interactome of MIWI2 in mouse fetal gonocytes undergoing de novo genome methylation and identify a previously unknown MIWI2-associated factor, SPOCD1, that is essential for the methylation and silencing of young transposable elements. The loss of Spocd1 in mice results in male-specific infertility but does not affect either piRNA biogenesis or the localization of MIWI2 to the nucleus. SPOCD1 is a nuclear protein whose expression is restricted to the period of de novo genome methylation. It co-purifies in vivo with DNMT3L and DNMT3A, components of the de novo methylation machinery, as well as with constituents of the NURD and BAF chromatin remodelling complexes. We propose a model whereby tethering of MIWI2 to a nascent transposable element transcript recruits repressive chromatin remodelling activities and the de novo methylation apparatus through SPOCD1. In summary, we have identified a previously unrecognized and essential executor of mammalian piRNA-directed DNA methylation.

Transcytosis via the late endocytic pathway as a cell morphogenetic mechanism

Plasma membranes fulfil many physiological functions. In polarized cells, different membrane compartments take on specialized roles, each being allocated correct amounts of membrane. The Drosophila tracheal system, an established tubulogenesis model, contains branched terminal cells with subcellular tubes formed by apical plasma membrane invagination. We show that apical endocytosis and late endosome‐mediated trafficking are required for membrane allocation to the apical and basal membrane domains. Basal plasma membrane growth stops if endocytosis is blocked, whereas the apical membrane grows excessively. Plasma membrane is initially delivered apically and then continuously endocytosed, together with apical and basal cargo. We describe an organelle carrying markers of late endosomes and multivesicular bodies (MVBs) that is abolished by inhibiting endocytosis and which we suggest acts as transit station for membrane destined to be redistributed both apically and basally. This is based on the observation that disrupting MVB formation prevents growth of both compartments.

Comprehensive Genetic Characterization of Mitochondrial Ca2+ Uniporter Components Reveals Their Different Physiological Requirements In Vivo

Mitochondrial Ca²⁺ uptake is an important mediator of metabolism and cell death. Identification of components of the highly conserved mitochondrial Ca²⁺ uniporter has opened it up to genetic analysis in model organisms. Here, we report a comprehensive genetic characterization of all known uniporter components conserved in Drosophila. While loss of pore-forming MCU or EMRE abolishes fast mitochondrial Ca²⁺ uptake, this results in only mild phenotypes when young, despite shortened lifespans. In contrast, loss of the MICU1 gatekeeper is developmentally lethal, consistent with unregulated Ca²⁺ uptake. Mutants for the neuronally restricted regulator MICU3 are viable with mild neurological impairment. Genetic interaction analyses reveal that MICU1 and MICU3 are not functionally interchangeable. More surprisingly, loss of MCU or EMRE does not suppress MICU1 mutant lethality, suggesting that this results from uniporter-independent functions. Our data reveal the interplay among components of the mitochondrial Ca²⁺ uniporter and shed light on their physiological requirements in vivo.

•

MCU or EMRE loss blocks fast mitochondrial calcium uptake but are relatively benign

•

MCU knockout flies are short lived compared to EMRE or MICU3 mutants

•

MICU1 mutants are developmentally lethal, and MCU or EMRE knockout fails to rescue

•

MICU1 and MICU3 are not functionally interchangeable

Overlapping but Distinct Sequences Play Roles in the Insulator and Promoter Activities of the Drosophila BEAF-Dependent scs' Insulator

Chromatin domain insulators are thought to help partition the genome into genetic units called topologically associating domains (TADs). In Drosophila, TADs are often separated by inter-TAD regions containing active housekeeping genes and associated insulator binding proteins. This raises the question of whether insulator binding proteins are involved primarily in chromosomal TAD architecture or gene activation, or if these two activities are linked. The Boundary Element-Associated Factor of 32 kDa (BEAF-32, or BEAF for short) is usually found in inter-TADs. BEAF was discovered based on binding to the scs' insulator, and is important for the insulator activity of scs' and other BEAF binding sites. There are divergent promoters in scs' with a BEAF binding site by each. Here, we dissect the scs' insulator to identify DNA sequences important for insulator and promoter activity, focusing on the half of scs' with a high affinity BEAF binding site. We find that the BEAF binding site is important for both insulator and promoter activity, as is another sequence we refer to as LS4. Aside from that, different sequences play roles in insulator and promoter activity. So while there is overlap and BEAF is important for both, insulator and promoter activity can be separated.

Survey for positively selected coding regions in the genome of the hematophagous tsetse fly Glossina morsitans identifies candidate genes associated with feeding habits and embryonic development

Tsetse flies are responsible for the transmission of Trypanossoma sp. to vertebrate animals in Africa causing huge health issues and economic loss. The availability of the genome sequence of Glossina morsitans enabled the discovery of several genes related to medically important phenotypes and novel physiological features. However, a genome-wide scan for coding regions that underwent positive selection is still missing, which is surprising given the evolution of traits associated with the hematophagy in this lineage. In this study, we employed an experimental design that controlled for the rate of false positives and we performed a scan of 3,318 G. morsitans genes. We found 145 genes with significant historical signal of positive selection. These genes were categorized into 18 functional classes after careful manual annotation. Based on their attributed functions, we identified candidate genes related with feeding habits and embryonic development. When our results were contrasted with gene expression data, we confirmed that most genes that underwent adaptive molecular evolution were frequently expressed in organs associated with key physiological evolutionary innovations in the G. morsitans lineage, namely, the salivary gland, the midgut, fat body tissue, and in the spermatophore.

A Nucleolar Isoform of the Drosophila Ubiquitin Specific Protease dUSP36 Regulates MYC-Dependent Cell Growth

The c-Myc oncogene is a transcription factor that regulates the expression of a very large set of genes mainly involved in cell growth and proliferation. It is overexpressed in more than 70% of human cancers, illustrating the importance of keeping its levels and activity under control. The ubiquitin proteasome system is a major regulator of MYC levels in humans as well as in model organisms such as Drosophila melanogaster. Although the E3 ligases that promote MYC ubiquitination have been largely investigated, the identity and the role of the deubiquitinating enzymes, which counteract their action is only beginning to be unraveled. Using isoform-specific CRISPR-Cas9 mutagenesis, we show that the Drosophila homolog of the Ubiquitin Specific Protease USP36 has different isoforms with specific sub-cellular localizations and that the nucleolar dUSP36-D isoform is specifically required for cell and organismal growth. We also demonstrate that this isoform interacts with dMYC and the E3 ligase AGO and regulates their stability and ubiquitination levels. Furthermore, we show that dUSP36 is ubiquitinated by AGO and is able to self-deubiquitinate. Finally, we provide in vivo evidence supporting the functional relevance of these regulatory relationships. Together these results reveal that dMYC, AGO and dUSP36 form a tripartite, evolutionary conserved complex that acts as a regulatory node to control dMYC protein levels.

Amyloid-like Assembly Activates a Phosphatase in the Developing Drosophila Embryo

Prion-like proteins can assume distinct conformational and physical states in the same cell. Sequence analysis suggests that prion-like proteins are prevalent in various species; however, it remains unclear what functional space they occupy in multicellular organisms. Here, we report the identification of a prion-like protein, Herzog (CG5830), through a multimodal screen in Drosophila melanogaster. Herzog functions as a membrane-associated phosphatase and controls embryonic patterning, likely being involved in TGF-β/BMP and FGF/EGF signaling pathways. Remarkably, monomeric Herzog is enzymatically inactive and becomes active upon amyloid-like assembly. The prion-like domain of Herzog is necessary for both its assembly and membrane targeting. Removal of the prion-like domain impairs activity, while restoring assembly on the membrane using a heterologous prion-like domain and membrane-targeting motif can restore phosphatase activity. This study provides an example of a prion-like domain that allows an enzyme to gain essential functionality via amyloid-like assembly to control animal development.

A Comprehensive Drosophila melanogaster Transcription Factor Interactome

Combinatorial interactions among transcription factors (TFs) play essential roles in generating gene expression specificity and diversity in metazoans. Using yeast 2-hybrid (Y2H) assays on nearly all sequence-specific Drosophila TFs, we identified 1,983 protein-protein interactions (PPIs), more than doubling the number of currently known PPIs among Drosophila TFs. For quality assessment, we validated a subset of our interactions using MITOMI and bimolecular fluorescence complementation assays. We combined our interactome with prior PPI data to generate an integrated Drosophila TF-TF binary interaction network. Our analysis of ChIP-seq data, integrating PPI and gene expression information, uncovered different modes by which interacting TFs are recruited to DNA. We further demonstrate the utility of our Drosophila interactome in shedding light on human TF-TF interactions. This study reveals how TFs interact to bind regulatory elements in vivo and serves as a resource of Drosophila TF-TF binary PPIs for understanding tissue-specific gene regulation.

Combinatorial regulation by transcription factors (TFs) is one mechanism for achieving condition and tissue-specific gene regulation. Shokri et al. mapped TF-TF interactions between most Drosophila TFs, reporting a comprehensive TF-TF network integrated with previously known interactions. They used this network to discern distinct TF-DNA binding modes.

Adaptive evolution among cytoplasmic piRNA proteins leads to decreased genomic auto-immunity

In metazoan germlines, the piRNA pathway acts as a genomic immune system, employing small RNA-mediated silencing to defend host DNA from the harmful effects of transposable elements (TEs). Expression of genomic TEs is proposed to initiate self regulation by increasing the production of repressive piRNAs, thereby “adapting” piRNA-mediated control to the most active TE families. Surprisingly, however, piRNA pathway proteins, which execute piRNA biogenesis and enforce silencing of targeted sequences, evolve rapidly and adaptively in animals. If TE silencing is ensured through piRNA biogenesis, what necessitates changes in piRNA pathway proteins? Here we used interspecific complementation to test for functional differences between Drosophila melanogaster and D. simulans alleles of three adaptively evolving piRNA pathway proteins: Armitage, Aubergine and Spindle-E. In contrast to piRNA-mediated transcriptional regulators examined in previous studies, these three proteins have cytoplasmic functions in piRNA maturation and post-transcriptional silencing. Across all three proteins we observed interspecific divergence in the regulation of only a handful of TE families, which were more robustly silenced by the heterospecific piRNA pathway protein. This unexpected result suggests that unlike transcriptional regulators, positive selection has not acted on cytoplasmic piRNA effector proteins to enhance their function in TE repression. Rather, TEs may evolve to “escape” silencing by host proteins. We further discovered that D. simulans alleles of aub and armi exhibit enhanced off-target effects on host transcripts in a D. melanogaster background, as well as modest reductions in the efficiency of piRNA biogenesis, suggesting that promiscuous binding of D. simulans Aub and Armi proteins to host transcripts reduces their participation in piRNA production. Avoidance of genomic auto-immunity may therefore be a critical target of selection. Our observations suggest that piRNA effector proteins are subject to an evolutionary trade-off between defending the host genome from the harmful effect of TEs while also minimizing collateral damage to host genes.

Transposable elements are mobile fragments of selfish DNA that burden host genomes with deleterious mutations and incite genome instability. Host cells employ a specialized small-RNA mediated silencing pathway, the piRNA pathway, to act as a genomic immune system suppressing the mobilization of TEs. Changes in genomic TE content are met with rapid changes in the piRNA pool, thereby maintaining host control over transposition. However, piRNA pathway proteins—which enact piRNA biogenesis and silence target TEs—also evolve adaptively. To isolate forces that underlie this adaptive evolution, we examined functional divergence between two Drosophila species for three adaptively evolving piRNA pathway proteins. To our surprise, we found very few differences in TE regulation, suggesting that evolution has not generally acted to enhance control of TE parasites. Rather, we discovered interspecific differences in the regulation of host mRNAs for two proteins, which suggested that proteins evolve to avoid off-target silencing of host transcripts. We propose that the avoidance of such “genomic autoimmunity” is an important and underappreciated force driving the adaptive evolution of piRNA proteins.

Gephebase, a database of genotype-phenotype relationships for natural and domesticated variation in Eukaryotes

Gephebase is a manually-curated database compiling our accumulated knowledge of the genes and mutations that underlie natural, domesticated and experimental phenotypic variation in all Eukaryotes—mostly animals, plants and yeasts. Gephebase aims to compile studies where the genotype–phenotype association (based on linkage mapping, association mapping or a candidate gene approach) is relatively well supported. Human clinical traits and aberrant mutant phenotypes in laboratory organisms are not included and can be found in other databases (e.g. OMIM, OMIA, Monarch Initiative). Gephebase contains more than 1700 entries. Each entry corresponds to an allelic difference at a given gene and its associated phenotypic change(s) between two species or two individuals of the same species, and is enriched with molecular details, taxonomic information, and bibliographic information. Users can easily browse entries and perform searches at various levels using boolean operators (e.g. transposable elements, snakes, carotenoid content, Doebley). Data is exportable in spreadsheet format. This database allows to perform meta-analyses to extract global trends about the living world and the research fields. Gephebase should also help breeders, conservationists and others to identify promising target genes for crop improvement, parasite/pest control, bioconservation and genetic diagnostic. It is freely available at www.gephebase.org.

Neuronal Function and Dopamine Signaling Evolve at High Temperature in Drosophila

Neuronal activity is temperature sensitive and affects behavioral traits important for individual fitness, such as locomotion and courtship. Yet, we do not know enough about the evolutionary response of neuronal phenotypes in new temperature environments. Here, we use long-term experimental evolution of Drosophila simulans populations exposed to novel temperature regimes. Here, we demonstrate a direct relationship between thermal selective pressure and the evolution of neuronally expressed molecular and behavioral phenotypes. Several essential neuronal genes evolve lower expression at high temperatures and higher expression at low temperatures, with dopaminergic neurons standing out by displaying the most consistent expression change across independent replicates. We functionally validate the link between evolved gene expression and behavioral changes by pharmacological intervention in the experimentally evolved D. simulans populations as well as by genetically triggered expression changes of key genes in D. melanogaster. As natural temperature clines confirm our results for Drosophila and Anopheles populations, we conclude that neuronal dopamine evolution is a key factor for temperature adaptation.

The RNAcentral Consortium http://orcid.org/0000-0002-6497-2883SweeneyBlake Ahttp://orcid.org/0000-0001-7279-2682PetrovAnton IBurkovBorishttp://orcid.org/0000-0001-8626-2148FinnRobert Dhttp://orcid.org/0000-0002-6982-4660BatemanAlexSzymanskiMaciejKarlowskiWojciech MGorodkinJanSeemannStefan ECannoneJamie JGutellRobin RFeyPetraBasuSiddharthaKaySimonhttp://orcid.org/0000-0001-7954-7057CochraneGuyBillisKostantinosEmmertDavidMarygoldSteven Jhttp://orcid.org/0000-0001-6718-3559HuntleyRachael Phttp://orcid.org/0000-0002-9791-0064LoveringRuth CFrankishAdamChanPatricia Phttp://orcid.org/0000-0003-3253-6021LoweTodd Mhttp://orcid.org/0000-0002-8380-5247BrufordElspethSealRuthhttp://orcid.org/0000-0001-6274-0184VandesompeleJohttp://orcid.org/0000-0002-2685-2637VoldersPieter-JanParaskevopoulouMariaMaLinaZhangZhangGriffiths-JonesSamBujnickiJanusz MBoccalettoPietrohttp://orcid.org/0000-0001-8522-334XBlakeJudith ABultCarol JChenRunshengZhaoYiWoodValerieRutherfordKimhttp://orcid.org/0000-0002-2084-269XRivasElenaColeJameshttp://orcid.org/0000-0001-5356-4174LaulederkindStanley J FShimoyamaMaryGillespieMarc EOrlic-MilacicMarijahttp://orcid.org/0000-0001-9424-9197KalvariIoannahttp://orcid.org/0000-0002-2497-3427NawrockiEricEngelStacia Rhttp://orcid.org/0000-0001-9163-5180CherryJ MichaelTeamSILVABerardiniTanya ZHatzigeorgiouArtemisKaragkouniDimitrahttp://orcid.org/0000-0002-1751-9226HoweKevinDavisPaulDingerMarcelhttp://orcid.org/0000-0002-7294-0865HeShunminYoshihamaMakiKenmochiNaoyaStadlerPeter FWilliamsKelly P. RNAcentral: a hub of information for non-coding RNA sequences

RNAcentral is a comprehensive database of non-coding RNA (ncRNA) sequences, collating information on ncRNA sequences of all types from a broad range of organisms. We have recently added a new genome mapping pipeline that identifies genomic locations for ncRNA sequences in 296 species. We have also added several new types of functional annotations, such as tRNA secondary structures, Gene Ontology annotations, and miRNA-target interactions. A new quality control mechanism based on Rfam family assignments identifies potential contamination, incomplete sequences, and more. The RNAcentral database has become a vital component of many workflows in the RNA community, serving as both the primary source of sequence data for academic and commercial groups, as well as a source of stable accessions for the annotation of genomic and functional features. These examples are facilitated by an improved RNAcentral web interface, which features an updated genome browser, a new sequence feature viewer, and improved text search functionality. RNAcentral is freely available at https://rnacentral.org.

Genenames.org: the HGNC and VGNC resources in 2019

The HUGO Gene Nomenclature Committee (HGNC) based at EMBL's European Bioinformatics Institute (EMBL-EBI) assigns unique symbols and names to human genes. There are over 40 000 approved gene symbols in our current database of which over 19 000 are for protein-coding genes. The Vertebrate Gene Nomenclature Committee (VGNC) was established in 2016 to assign standardized nomenclature in line with human for vertebrate species that lack their own nomenclature committees. The VGNC initially assigned nomenclature for over 15000 protein-coding genes in chimpanzee. We have extended this process to other vertebrate species, naming over 14000 protein-coding genes in cow and dog and over 13 000 in horse to date. Our HGNC website https://www.genenames.org has undergone a major design update, simplifying the homepage to provide easy access to our search tools and making the site more mobile friendly. Our gene families pages are now known as 'gene groups' and have increased in number to over 1200, with nearly half of all named genes currently assigned to at least one gene group. This article provides an overview of our online data and resources, focusing on our work over the last two years.

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.

REDfly: the transcriptional regulatory element database for Drosophila

The REDfly database provides a comprehensive curation of experimentally-validated Drosophila transcriptional cis-regulatory elements and includes information on DNA sequence, experimental evidence, patterns of regulated gene expression, and more. Now in its thirteenth year, REDfly has grown to over 23 000 records of tested reporter gene constructs and 2200 tested transcription factor binding sites. Recent developments include the start of curation of predicted cis-regulatory modules in addition to experimentally-verified ones, improved search and filtering, and increased interaction with the authors of curated papers. An expanded data model that will capture information on temporal aspects of gene regulation, regulation in response to environmental and other non-developmental cues, sexually dimorphic gene regulation, and non-endogenous (ectopic) aspects of reporter gene expression is under development and expected to be in place within the coming year. REDfly is freely accessible at http://redfly.ccr.buffalo.edu, and news about database updates and new features can be followed on Twitter at @REDfly_database.

The metabolome as a link in the genotype-phenotype map for peroxide resistance in the fruit fly, Drosophila melanogaster

BACKGROUND

Genetic association studies that seek to explain the inheritance of complex traits typically fail to explain a majority of the heritability of the trait under study. Thus, we are left with a gap in the map from genotype to phenotype. Several approaches have been used to fill this gap, including those that attempt to map endophenotype such as the transcriptome, proteome or metabolome, that underlie complex traits. Here we used metabolomics to explore the nature of genetic variation for hydrogen peroxide (H2O2) resistance in the sequenced inbred Drosophila Genetic Reference Panel (DGRP).

RESULTS

We first studied genetic variation for H2O2 resistance in 179 DGRP lines and along with identifying the insulin signaling modulator u-shaped and several regulators of feeding behavior, we estimate that a substantial amount of phenotypic variation can be explained by a polygenic model of genetic variation. We then profiled a portion of the aqueous metabolome in subsets of eight 'high resistance' lines and eight 'low resistance' lines. We used these lines to represent collections of genotypes that were either resistant or sensitive to the stressor, effectively modeling a discrete trait. Across the range of genotypes in both populations, flies exhibited surprising consistency in their metabolomic signature of resistance. Importantly, the resistance phenotype of these flies was more easily distinguished by their metabolome profiles than by their genotypes. Furthermore, we found a metabolic response to H2O2 in sensitive, but not in resistant genotypes. Metabolomic data further implicated at least two pathways, glycogen and folate metabolism, as determinants of sensitivity to H2O2. We also discovered a confounding effect of feeding behavior on assays involving supplemented food.

CONCLUSIONS

This work suggests that the metabolome can be a point of convergence for genetic variation influencing complex traits, and can efficiently elucidate mechanisms underlying trait variation.

Identification of genomic enhancers through spatial integration of single-cell transcriptomics and epigenomics

Single-cell technologies allow measuring chromatin accessibility and gene expression in each cell, but jointly utilizing both layers to map bona fide gene regulatory networks and enhancers remains challenging. Here, we generate independent single-cell RNA-seq and single-cell ATAC-seq atlases of the Drosophila eye-antennal disc and spatially integrate the data into a virtual latent space that mimics the organization of the 2D tissue using ScoMAP (Single-Cell Omics Mapping into spatial Axes using Pseudotime ordering). To validate spatially predicted enhancers, we use a large collection of enhancer-reporter lines and identify ~ 85% of enhancers in which chromatin accessibility and enhancer activity are coupled. Next, we infer enhancer-to-gene relationships in the virtual space, finding that genes are mostly regulated by multiple, often redundant, enhancers. Exploiting cell type-specific enhancers, we deconvolute cell type-specific effects of bulk-derived chromatin accessibility QTLs. Finally, we discover that Prospero drives neuronal differentiation through the binding of a GGG motif. In summary, we provide a comprehensive spatial characterization of gene regulation in a 2D tissue.

Promoter-Proximal Chromatin Domain Insulator Protein BEAF Mediates Local and Long-Range Communication with a Transcription Factor and Directly Activates a Housekeeping Promoter in Drosophila

BEAF (Boundary Element-Associated Factor) was originally identified as a Drosophila melanogaster chromatin domain insulator-binding protein, suggesting a role in gene regulation through chromatin organization and dynamics. Genome-wide mapping found that BEAF usually binds near transcription start sites, often of housekeeping genes, suggesting a role in promoter function. This would be a nontraditional role for an insulator-binding protein. To gain insight into molecular mechanisms of BEAF function, we identified interacting proteins using yeast two-hybrid assays. Here, we focus on the transcription factor Serendipity δ (Sry-δ). Interactions were confirmed in pull-down experiments using bacterially expressed proteins, by bimolecular fluorescence complementation, and in a genetic assay in transgenic flies. Sry-δ interacted with promoter-proximal BEAF both when bound to DNA adjacent to BEAF or > 2-kb upstream to activate a reporter gene in transient transfection experiments. The interaction between BEAF and Sry-δ was detected using both a minimal developmental promoter (y) and a housekeeping promoter (RpS12), while BEAF alone strongly activated the housekeeping promoter. These two functions for BEAF implicate it in playing a direct role in gene regulation at hundreds of BEAF-associated promoters.