FlyBase at 25: looking to the future

CapsNet-TIS: Predicting translation initiation site based on multi-feature fusion and improved capsule network

Genes are the basic units of protein synthesis in organisms, and accurately identifying the translation initiation site (TIS) of genes is crucial for understanding the regulation, transcription, and translation processes of genes. However, the existing models cannot adequately extract the feature information in TIS sequences, and they also inadequately capture the complex hierarchical relationships among features. Therefore, a novel predictor named CapsNet-TIS is proposed in this paper. CapsNet-TIS first fully extracts the TIS sequence information using four encoding methods, including One-hot encoding, physical structure property (PSP) encoding, nucleotide chemical property (NCP) encoding, and nucleotide density (ND) encoding. Next, multi-scale convolutional neural networks are used to perform feature fusion of the encoded features to enhance the comprehensiveness of the feature representation. Finally, the fused features are classified using capsule network as the main network of the classification model to capture the complex hierarchical relationships among the features. Moreover, we improve the capsule network by introducing residual block, channel attention, and BiLSTM to enhance the model's feature extraction and sequence data modeling capabilities. In this paper, the performance of CapsNet-TIS is evaluated using TIS datasets from four species: human, mouse, bovine, and fruit fly, and the effectiveness of each part is demonstrated by performing ablation experiments. By comparing the experimental results with models proposed by other researchers, the results demonstrate the superior performance of CapsNet-TIS.

Temporal genomics in Hawaiian crickets reveals compensatory intragenomic coadaptation during adaptive evolution

Theory predicts that compensatory genetic changes reduce negative indirect effects of selected variants during adaptive evolution, but evidence is scarce. Here, we test this in a wild population of Hawaiian crickets using temporal genomics and a high-quality chromosome-level cricket genome. In this population, a mutation, flatwing, silences males and rapidly spread due to an acoustically-orienting parasitoid. Our sampling spanned a social transition during which flatwing fixed and the population went silent. We find long-range linkage disequilibrium around the putative flatwing locus was maintained over time, and hitchhiking genes had functions related to negative flatwing-associated effects. We develop a combinatorial enrichment approach using transcriptome data to test for compensatory, intragenomic coevolution. Temporal changes in genomic selection were distributed genome-wide and functionally associated with the population's transition to silence, particularly behavioural responses to silent environments. Our results demonstrate how 'adaptation begets adaptation'; changes to the sociogenetic environment accompanying rapid trait evolution can generate selection provoking further, compensatory adaptation.

Using EMBL-EBI Services via Web Interface and Programmatically via Web Services

The European Bioinformatics Institute (EMBL-EBI)'s Job Dispatcher framework provides access to a wide range of core databases and analysis tools that are of key importance in bioinformatics. As well as providing web interfaces to these resources, web services are available using REST and SOAP protocols that enable programmatic access and allow their integration into other applications and analytical workflows and pipelines. This article describes the various options available to researchers and bioinformaticians who would like to use our resources via the web interface employing RESTful web services clients provided in Perl, Python, and Java or who would like to use Docker containers to integrate the resources into analysis pipelines and workflows. © 2024 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: Retrieving data from EMBL-EBI using Dbfetch via the web interface Alternate Protocol 1: Retrieving data from EMBL-EBI using WSDbfetch via the REST interface Alternate Protocol 2: Retrieving data from EMBL-EBI using Dbfetch via RESTful web services with Python client Support Protocol 1: Installing Python REST web services clients Basic Protocol 2: Sequence similarity search using FASTA search via the web interface Alternate Protocol 3: Sequence similarity search using FASTA via RESTful web services with Perl client Support Protocol 2: Installing Perl REST web services clients Basic Protocol 3: Sequence similarity search using NCBI BLAST+ RESTful web services with Python client Basic Protocol 4: Sequence similarity search using HMMER3 phmmer REST web services with Perl client and Docker Support Protocol 3: Installing Docker and running the EMBL-EBI client container Basic Protocol 5: Protein functional analysis using InterProScan 5 RESTful web services with the Python client and Docker Alternate Protocol 4: Protein functional analysis using InterProScan 5 RESTful web services with the Java client Support Protocol 4: Installing Java web services clients Basic Protocol 6: Multiple sequence alignment using Clustal Omega via web interface Alternate Protocol 5: Multiple sequence alignment using Clustal Omega with Perl client and Docker Support Protocol 5: Exploring the RESTful API with OpenAPI User Inferface.

Accelerated protein retention expansion microscopy using microwave radiation

Protein retention expansion microscopy (ExM) retains genetically encoded fluorescent proteins or antibody-conjugated fluorescent probes in fixed tissue and isotropically expands the tissue through a swellable polymer network to allow nanoscale (<70 nm) resolution on diffraction-limited confocal microscopes. Despite numerous advantages ExM brings to biological studies, the full protocol is time-consuming and can take multiple days to complete. Here, we adapted the ExM protocol to the vibratome-sectioned brain tissue of Xenopus laevis tadpoles and implemented a microwave-assisted protocol to reduce the workflow from days to hours. In addition to the significantly accelerated processing time, our microwave-assisted ExM (M/WExM) protocol maintains the superior resolution and signal-to-noise ratio of the original ExM protocol. Furthermore, the M/WExM protocol yields higher magnitude of expansion, suggesting that in addition to accelerating the process through increased diffusion rate of reagents, microwave radiation may also facilitate the expansion process. To demonstrate the applicability of this method to other specimens and protocols, we adapted the microwave-accelerated protocol to whole mount adult brain tissue of Drosophila melanogaster fruit flies, and successfully reduced the total processing time of a widely-used Drosophila IHC-ExM protocol from 6 days to 2 days. Our results demonstrate that with appropriate adjustment of the microwave parameters (wattage, pulse duration, interval, and number of cycles), this protocol can be readily adapted to different model organisms and tissue types to greatly increase the efficiency of ExM experiments.

Comparative interactome analysis of α-arrestin families in human and Drosophila

The α-arrestins form a large family of evolutionally conserved modulators that control diverse signaling pathways, including both G-protein-coupled receptor (GPCR)-mediated and non-GPCR-mediated pathways, across eukaryotes. However, unlike β-arrestins, only a few α-arrestin targets and functions have been characterized. Here, using affinity purification and mass spectrometry, we constructed interactomes for 6 human and 12 Drosophila α-arrestins. The resulting high-confidence interactomes comprised 307 and 467 prey proteins in human and Drosophila, respectively. A comparative analysis of these interactomes predicted not only conserved binding partners, such as motor proteins, proteases, ubiquitin ligases, RNA splicing factors, and GTPase-activating proteins, but also those specific to mammals, such as histone modifiers and the subunits of V-type ATPase. Given the manifestation of the interaction between the human α-arrestin, TXNIP, and the histone-modifying enzymes, including HDAC2, we undertook a global analysis of transcription signals and chromatin structures that were affected by TXNIP knockdown. We found that TXNIP activated targets by blocking HDAC2 recruitment to targets, a result that was validated by chromatin immunoprecipitation assays. Additionally, the interactome for an uncharacterized human α-arrestin ARRDC5 uncovered multiple components in the V-type ATPase, which plays a key role in bone resorption by osteoclasts. Our study presents conserved and species-specific protein-protein interaction maps for α-arrestins, which provide a valuable resource for interrogating their cellular functions for both basic and clinical research.

Evolutionary diversification reveals distinct somatic versus germline cytoskeletal functions of the Arp2 branched actin nucleator protein

Branched actin networks are critical in many cellular processes, including cell motility and division. Arp2, a protein within the seven-membered Arp2/3 complex, is responsible for generating branched actin. Given its essential roles, Arp2 evolves under stringent sequence conservation throughout eukaryotic evolution. We unexpectedly discovered recurrent evolutionary diversification of Arp2 in Drosophila, yielding independently arising paralogs Arp2D in obscura species and Arp2D2 in montium species. Both paralogs are unusually testis-enriched in expression relative to Arp2. We investigated whether their sequence divergence from canonical Arp2 led to functional specialization by replacing Arp2 in D. melanogaster with either Arp2D or Arp2D2. Despite their divergence, we surprisingly found that both complement Arp2's essential function in somatic tissue, suggesting they have preserved the ability to polymerize branched actin even in a non-native species. However, we found that Arp2D- and Arp2D2-expressing males display defects throughout sperm development, with Arp2D resulting in more pronounced deficiencies and subfertility, suggesting the Arp2 paralogs are cross-species incompatible in the testis. We focused on Arp2D and pinpointed two highly diverged structural regions-the D-loop and C terminus-and found that they contribute to germline defects in D. melanogaster sperm development. However, while the Arp2D C terminus is suboptimal in the D. melanogaster testis, it is essential for Arp2D somatic function. Testis cytology of the paralogs' native species revealed striking differences in germline actin structures, indicating unique cytoskeletal requirements. Our findings suggest canonical Arp2 function differs between somatic versus germline contexts, and Arp2 paralogs may have recurrently evolved for species-specialized actin branching in the testis.

Two-step CRISPR-Cas9 protocol for transposable element deletion in D. melanogaster natural populations

We present a protocol for generating a precise deletion, without altering the genetic background of the strain, of a transposable element (TE) in a natural population of Drosophila melanogaster using two steps of CRISPR-Cas9 homology-directed repair. We describe steps for replacing the TE by a fluorescent marker and for subsequent marker removal using single-guide RNAs, repair plasmids, and microinjection. We also detail steps for screening the deletion of the TE and generating a homozygous mutant strain. For complete details on the use and execution of this protocol, please refer to Merenciano and Gonzalez.1.

TnaA, a trithorax group protein, modulates wingless expression in different regions of the Drosophila wing imaginal disc

wingless expression is exquisitely regulated by different factors and enhancers in the imaginal wing discs of Drosophila melanogaster in four domains: the dorsal band, the dorso-ventral boundary, and the inner and outer ring domains. tonalli is a trithorax group gene that encodes a putative SUMO E3 ligase that binds to chromatin to regulate the expression of its targets, including the Hox genes. However, its role in modulating gene expression is barely known. Here, we show that TnaA modulates the wingless expression at two domains of the wing disc, the dorso-ventral boundary and the inner ring. At first, tonalli interacts genetically with Notch to form the wing margin. In the inner ring domain, TnaA modulates wingless transcription. When the dosage of TnaA increases in or near the inner ring since early larval stages, this domain expands with a rapid increase in wingless expression. TnaA occupies the wingless Inner Ring Enhancer at the wing disc, meanwhile it does not affect wingless expression directed by the Ventral Disc Enhancer in leg discs, suggesting that TnaA acts as a wingless enhancer-specific factor. We describe for the first time the presence of TnaA at the Inner Ring Enhancer as a specific regulator of wingless in the development of wing boundaries.

Drosophila as a Robust Model System for Assessing Autophagy: A Review

Autophagy is the process through which a body breaks down and recycles its own cellular components, primarily inside lysosomes. It is a cellular response to starvation and stress, which plays decisive roles in various biological processes such as senescence, apoptosis, carcinoma, and immune response. Autophagy, which was first discovered as a survival mechanism during starvation in yeast, is now known to serve a wide range of functions in more advanced organisms. It plays a vital role in how cells respond to stress, starvation, and infection. While research on yeast has led to the identification of many key components of the autophagy process, more research into autophagy in more complex systems is still warranted. This review article focuses on the use of the fruit fly Drosophila melanogaster as a robust testing model in further research on autophagy. Drosophila provides an ideal environment for exploring autophagy in a living organism during its development. Additionally, Drosophila is a well-suited compact tool for genetic analysis in that it serves as an intermediate between yeast and mammals because evolution conserved the molecular machinery required for autophagy in this species. Experimental tractability of host-pathogen interactions in Drosophila also affords great convenience in modeling human diseases on analogous structures and tissues.

Midgut transcriptomic responses to dengue and chikungunya viruses in the vectors Aedes albopictus and Aedes malayensis

Dengue (DENV) and chikungunya (CHIKV) viruses are among the most preponderant arboviruses. Although primarily transmitted through the bite of Aedes aegypti mosquitoes, Aedes albopictus and Aedes malayensis are competent vectors and have an impact on arbovirus epidemiology. Here, to fill the gap in our understanding of the molecular interactions between secondary vectors and arboviruses, we used transcriptomics to profile the whole-genome responses of A. albopictus to CHIKV and of A. malayensis to CHIKV and DENV at 1 and 4 days post-infection (dpi) in midguts. In A. albopictus, 1793 and 339 genes were significantly regulated by CHIKV at 1 and 4 dpi, respectively. In A. malayensis, 943 and 222 genes upon CHIKV infection, and 74 and 69 genes upon DENV infection were significantly regulated at 1 and 4 dpi, respectively. We reported 81 genes that were consistently differentially regulated in all the CHIKV-infected conditions, identifying a CHIKV-induced signature. We identified expressed immune genes in both mosquito species, using a de novo assembled midgut transcriptome for A. malayensis, and described the immune architectures. We found the JNK pathway activated in all conditions, generalizing its antiviral function to Aedines. Our comprehensive study provides insight into arbovirus transmission by multiple Aedes vectors.

Tau polarizes an aging transcriptional signature to excitatory neurons and glia

Aging is a major risk factor for Alzheimer's disease (AD), and cell-type vulnerability underlies its characteristic clinical manifestations. We have performed longitudinal, single-cell RNA-sequencing in Drosophila with pan-neuronal expression of human tau, which forms AD neurofibrillary tangle pathology. Whereas tau- and aging-induced gene expression strongly overlap (93%), they differ in the affected cell types. In contrast to the broad impact of aging, tau-triggered changes are strongly polarized to excitatory neurons and glia. Further, tau can either activate or suppress innate immune gene expression signatures in a cell-type-specific manner. Integration of cellular abundance and gene expression pinpoints nuclear factor kappa B signaling in neurons as a marker for cellular vulnerability. We also highlight the conservation of cell-type-specific transcriptional patterns between Drosophila and human postmortem brain tissue. Overall, our results create a resource for dissection of dynamic, age-dependent gene expression changes at cellular resolution in a genetically tractable model of tauopathy.

P-element invasion fuels molecular adaptation in laboratory populations of Drosophila melanogaster

Transposable elements (TEs) are mobile genetic parasites that frequently invade new host genomes through horizontal transfer. Invading TEs often exhibit a burst of transposition, followed by reduced transposition rates as repression evolves in the host. We recreated the horizontal transfer of P-element DNA transposons into a Drosophila melanogaster host and followed the expansion of TE copies and evolution of host repression in replicate laboratory populations reared at different temperatures. We observed that while populations maintained at high temperatures rapidly go extinct after TE invasion, those maintained at lower temperatures persist, allowing for TE spread and the evolution of host repression. We also surprisingly discovered that invaded populations experienced recurrent insertion of P-elements into a specific long non-coding RNA, lncRNA:CR43651, and that these insertion alleles are segregating at unusually high frequency in experimental populations, indicative of positive selection. We propose that, in addition to driving the evolution of repression, transpositional bursts of invading TEs can drive molecular adaptation.

GTPase-activating protein TBC1D5 coordinates with retromer to constrain synaptic growth by inhibiting BMP signaling

Formation and plasticity of neural circuits rely on precise regulation of synaptic growth. At Drosophila neuromuscular junction (NMJ), Bone Morphogenetic Protein (BMP) signaling is critical for many aspects of synapse formation and function. The evolutionarily conserved retromer complex and its associated GTPase-activating protein TBC1D5 are critical regulators of membrane trafficking and cellular signaling. However, their functions in regulating the formation of NMJ are less understood. Here, we report that TBC1D5 is required for inhibition of synaptic growth, and loss of TBC1D5 leads to abnormal presynaptic terminal development, including excessive satellite boutons and branch formation. Ultrastructure analysis reveals that the size of synaptic vesicles and the density of subsynaptic reticulum are increased in TBC1D5 mutant boutons. Disruption of interactions of TBC1D5 with Rab7 and retromer phenocopies the loss of TBC1D5. Unexpectedly, we find that TBC1D5 is functionally linked to Rab6, in addition to Rab7, to regulate synaptic growth. Mechanistically, we show that loss of TBC1D5 leads to upregulated BMP signaling by increasing the protein level of BMP type II receptor Wishful Thinking (Wit) at NMJ. Overall, our data establish that TBC1D5 in coordination with retromer constrains synaptic growth by regulating Rab7 activity, which negatively regulates BMP signaling through inhibiting Wit level.

Evolutionary diversification reveals distinct somatic versus germline cytoskeletal functions of the Arp2 branched actin nucleator protein

Branched actin networks are critical in many cellular processes, including cell motility and division. Arp2, a protein within the 7-membered Arp2/3 complex, is responsible for generating branched actin. Given its essential roles, Arp2 evolves under stringent sequence conservation throughout eukaryotic evolution. We unexpectedly discovered recurrent evolutionary diversification of Arp2 in Drosophila, yielding independently arising paralogs Arp2D in obscura species and Arp2D2 in montium species. Both paralogs are unusually testis-enriched in expression relative to Arp2. We investigated whether their sequence divergence from canonical Arp2 led to functional specialization by replacing Arp2 in D. melanogaster with either Arp2D or Arp2D2. Despite their divergence, we surprisingly found both complement Arp2's essential function in the soma, suggesting they have preserved the ability to polymerize branched actin even in a non-native species. However, we found that Arp2D-expressing males are subfertile and display many defects throughout sperm development. We pinpointed two highly diverged structural regions in Arp2D that contribute to these defects: subdomain 2 and the C-terminus. We expected that germline function would be rescued by replacing Arp2D's long and charged C-terminus with Arp2's short C-terminus, yet surprisingly, the essential somatic function of Arp2D was lost. Therefore, while Arp2D's structural divergence is incompatible with D. melanogaster sperm development, its unique C-terminus has evolved a critical role in actin polymerization. Our findings suggest canonical Arp2's function differs between somatic versus germline contexts, and Arp2 paralogs have recurrently evolved and specialized for actin branching in the testis.

Dietary Restriction Impacts Peripheral Circadian Clock Output Important for Longevity in Drosophila

Circadian clocks may mediate lifespan extension by caloric or dietary restriction (DR). We find that the core clock transcription factor Clock is crucial for a robust longevity and fecundity response to DR in Drosophila. To identify clock-controlled mediators, we performed RNA-sequencing from abdominal fat bodies across the 24 h day after just 5 days under control or DR diets. In contrast to more chronic DR regimens, we did not detect significant changes in the rhythmic expression of core clock genes. Yet we discovered that DR induced de novo rhythmicity or increased expression of rhythmic clock output genes. Network analysis revealed that DR increased network connectivity in one module comprised of genes encoding proteasome subunits. Adult, fat body specific RNAi knockdown demonstrated that proteasome subunits contribute to DR-mediated lifespan extension. Thus, clock control of output links DR-mediated changes in rhythmic transcription to lifespan extension.

Enrichment of Hard Sweeps on the X Chromosome in Drosophila melanogaster

The characteristic properties of the X chromosome, such as male hemizygosity and its unique inheritance pattern, expose it to natural selection in a way that can be different from the autosomes. Here, we investigate the differences in the tempo and mode of adaptation on the X chromosome and autosomes in a population of Drosophila melanogaster. Specifically, we test the hypothesis that due to hemizygosity and a lower effective population size on the X, the relative proportion of hard sweeps, which are expected when adaptation is gradual, compared with soft sweeps, which are expected when adaptation is rapid, is greater on the X than on the autosomes. We quantify the incidence of hard versus soft sweeps in North American D. melanogaster population genomic data with haplotype homozygosity statistics and find an enrichment of the proportion of hard versus soft sweeps on the X chromosome compared with the autosomes, confirming predictions we make from simulations. Understanding these differences may enable a deeper understanding of how important phenotypes arise as well as the impact of fundamental evolutionary parameters on adaptation, such as dominance, sex-specific selection, and sex-biased demography.

Genetic variation in P-element dysgenic sterility is associated with double-strand break repair and alternative splicing of TE transcripts

The germline mobilization of transposable elements (TEs) by small RNA mediated silencing pathways is conserved across eukaryotes and critical for ensuring the integrity of gamete genomes. However, genomes are recurrently invaded by novel TEs through horizontal transfer. These invading TEs are not targeted by host small RNAs, and their unregulated activity can cause DNA damage in germline cells and ultimately lead to sterility. Here we use hybrid dysgenesis-a sterility syndrome of Drosophila caused by transposition of invading P-element DNA transposons-to uncover host genetic variants that modulate dysgenic sterility. Using a panel of highly recombinant inbred lines of Drosophila melanogaster, we identified two linked quantitative trait loci (QTL) that determine the severity of dysgenic sterility in young and old females, respectively. We show that ovaries of fertile genotypes exhibit increased expression of splicing factors that suppress the production of transposase encoding transcripts, which likely reduces the transposition rate and associated DNA damage. We also show that fertile alleles are associated with decreased sensitivity to double-stranded breaks and enhanced DNA repair, explaining their ability to withstand high germline transposition rates. Together, our work reveals a diversity of mechanisms whereby host genotype modulates the cost of an invading TE, and points to genetic variants that were likely beneficial during the P-element invasion.

The importance of being the HGNC

The HUGO Gene Nomenclature Committee (HGNC) has been providing standardized symbols and names for human genes since the late 1970s. As funding agencies change their priorities, finding financial support for critical biomedical resources such as the HGNC becomes ever more challenging. In this article, we outline the key roles the HGNC currently plays in aiding communication and the need for these activities to be maintained.

The Drosophila BEAF insulator protein interacts with the polybromo subunit of the PBAP chromatin remodeling complex

The Drosophila Boundary Element-Associated Factor of 32 kDa (BEAF) binds in promoter regions of a few thousand mostly housekeeping genes. BEAF is implicated in both chromatin domain boundary activity and promoter function, although molecular mechanisms remain elusive. Here, we show that BEAF physically interacts with the polybromo subunit (Pbro) of PBAP, a SWI/SNF-class chromatin remodeling complex. BEAF also shows genetic interactions with Pbro and other PBAP subunits. We examine the effect of this interaction on gene expression and chromatin structure using precision run-on sequencing and micrococcal nuclease sequencing after RNAi-mediated knockdown in cultured S2 cells. Our results are consistent with the interaction playing a subtle role in gene activation. Fewer than 5% of BEAF-associated genes were significantly affected after BEAF knockdown. Most were downregulated, accompanied by fill-in of the promoter nucleosome-depleted region and a slight upstream shift of the +1 nucleosome. Pbro knockdown caused downregulation of several hundred genes and showed a correlation with BEAF knockdown but a better correlation with promoter-proximal GAGA factor binding. Micrococcal nuclease sequencing supports that BEAF binds near housekeeping gene promoters while Pbro is more important at regulated genes. Yet there is a similar general but slight reduction of promoter-proximal pausing by RNA polymerase II and increase in nucleosome-depleted region nucleosome occupancy after knockdown of either protein. We discuss the possibility of redundant factors keeping BEAF-associated promoters active and masking the role of interactions between BEAF and the Pbro subunit of PBAP in S2 cells. We identify Facilitates Chromatin Transcription (FACT) and Nucleosome Remodeling Factor (NURF) as candidate redundant factors.

The insect perspective on Z-disc structure and biology

Myofibrils are the intracellular structures formed by actin and myosin filaments. They are paracrystalline contractile cables with unusually well-defined dimensions. The sliding of actin past myosin filaments powers contractions, and the entire system is held in place by a structure called the Z-disc, which anchors the actin filaments. Myosin filaments, in turn, are anchored to another structure called the M-line. Most of the complex architecture of myofibrils can be reduced to studying the Z-disc, and recently, important advances regarding the arrangement and function of Z-discs in insects have been published. On a very small scale, we have detailed protein structure information. At the medium scale, we have cryo-electron microscopy maps, super-resolution microscopy and protein-protein interaction networks, while at the functional scale, phenotypic data are available from precise genetic manipulations. All these data aim to answer how the Z-disc works and how it is assembled. Here, we summarize recent data from insects and explore how it fits into our view of the Z-disc, myofibrils and, ultimately, muscles.

Divergent evolutionary trajectories shape the postmating transcriptional profiles of conspecifically and heterospecifically mated cactophilic Drosophila females

Postmating-prezygotic (PMPZ) reproductive isolation is hypothesized to result from divergent coevolutionary trajectories of sexual selection and/or sexual conflict in isolated populations. However, the genetic basis of PMPZ incompatibilities between species is poorly understood. Here, we use a comparative framework to compare global gene expression in con- and heterospecifically mated Drosophila mojavensis and D. arizonae female reproductive tracts. We find striking divergence between the species in the female postmating transcriptional response to conspecific mating, including differences in differential expression (DE), alternative splicing (AS), and intron retention (IR). As predicted, heterospecific matings produce disrupted transcriptional profiles, but the overall patterns of misregulation are different between the reciprocal crosses. Moreover, we find a positive correlation between postmating transcriptional divergence between species and levels of transcriptional disruption in heterospecific crosses. This result indicates that mating responsive genes that have diverged more in expression also have more disrupted transcriptional profiles in heterospecifically mated females. Overall, our results provide insights into the evolution of PMPZ isolation and lay the foundation for future studies aimed at identifying specific genes involved in PMPZ incompatibilities and the evolutionary forces that have contributed to their divergence in closely related species.

Comparison of global gene expression patterns in con- and heterospecifically mated Drosophila mojavensis and Drosophila arizonae suggest that mating-responsive genes with divergent expression also exhibit more disrupted transcriptional profiles in heterospecifically mated females, providing further insight into the evolution of postmating-prezygotic reproductive isolation.

Dynamics of sex-biased gene expression during development in the stick insect Timema californicum

Sexually dimorphic phenotypes are thought to arise primarily from sex-biased gene expression during development. Major changes in developmental strategies, such as the shift from hemimetabolous to holometabolous development, are therefore expected to have profound consequences for the dynamics of sex-biased gene expression. However, no studies have previously examined sex-biased gene expression during development in hemimetabolous insects, precluding comparisons between developmental strategies. Here we characterized sex-biased gene expression at three developmental stages in a hemimetabolous stick insect (Timema californicum): hatchlings, juveniles, and adults. As expected, the proportion of sex-biased genes gradually increased during development, mirroring the gradual increase of phenotypic sexual dimorphism. Sex-biased genes identified at early developmental stages were generally consistently male- or female-biased at later stages, suggesting their importance in sexual differentiation. Additionally, we compared the dynamics of sex-biased gene expression during development in T. californicum to those of the holometabolous fly Drosophila melanogaster by reanalyzing publicly available RNA-seq data from third instar larval, pupal and adult stages. In D. melanogaster, 84% of genes were sex-biased at the adult stage (compared to only 20% in T. californicum), and sex-biased gene expression increased abruptly at the adult stage when morphological sexual dimorphism is manifested. Our findings are consistent with the prediction that the dynamics of sex-biased gene expression during development differ extensively between holometabolous and hemimetabolous insect species.

DYRK1a Inhibitor Mediated Rescue of Drosophila Models of Alzheimer’s Disease-Down Syndrome Phenotypes

Alzheimer’s disease (AD) is the most common neurodegenerative disease which is becoming increasingly prevalent due to ageing populations resulting in huge social, economic, and health costs to the community. Despite the pathological processing of genes such as Amyloid Precursor Protein (APP) into Amyloid-β and Microtubule Associated Protein Tau (MAPT) gene, into hyperphosphorylated Tau tangles being known for decades, there remains no treatments to halt disease progression. One population with increased risk of AD are people with Down syndrome (DS), who have a 90% lifetime incidence of AD, due to trisomy of human chromosome 21 (HSA21) resulting in three copies of APP and other AD-associated genes, such as DYRK1A (Dual specificity tyrosine-phosphorylation-regulated kinase 1A) overexpression. This suggests that blocking DYRK1A might have therapeutic potential. However, it is still not clear to what extent DYRK1A overexpression by itself leads to AD-like phenotypes and how these compare to Tau and Amyloid-β mediated pathology. Likewise, it is still not known how effective a DYRK1A antagonist may be at preventing or improving any Tau, Amyloid-β and DYRK1a mediated phenotype. To address these outstanding questions, we characterised Drosophila models with targeted overexpression of human Tau, human Amyloid-β or the fly orthologue of DYRK1A, called minibrain (mnb). We found targeted overexpression of these AD-associated genes caused degeneration of photoreceptor neurons, shortened lifespan, as well as causing loss of locomotor performance, sleep, and memory. Treatment with the experimental DYRK1A inhibitor PST-001 decreased pathological phosphorylation of human Tau [at serine (S) 262]. PST-001 reduced degeneration caused by human Tau, Amyloid-β or mnb lengthening lifespan as well as improving locomotion, sleep and memory loss caused by expression of these AD and DS genes. This demonstrated PST-001 effectiveness as a potential new therapeutic targeting AD and DS pathology.

Pseudocleavage furrows restrict plasma membrane-associated PH domain in syncytial Drosophila embryos

Syncytial cells contain multiple nuclei and have local distribution and function of cellular components despite being synthesized in a common cytoplasm. The syncytial Drosophila blastoderm embryo shows reduced spread of organelle and plasma membrane-associated proteins between adjacent nucleo-cytoplasmic domains. Anchoring to the cytoarchitecture within a nucleo-cytoplasmic domain is likely to decrease the spread of molecules; however, its role in restricting this spread has not been assessed. In order to analyze the cellular mechanisms that regulate the rate of spread of plasma membrane-associated molecules in the syncytial Drosophila embryos, we express a pleckstrin homology (PH) domain in a localized manner at the anterior of the embryo by tagging it with the bicoid mRNA localization signal. Anteriorly expressed PH domain forms an exponential gradient in the anteroposterior axis with a longer length scale compared with Bicoid. Using a combination of experiments and theoretical modeling, we find that the characteristic distribution and length scale emerge due to plasma membrane sequestration and restriction within an energid. Loss of plasma membrane remodeling to form pseudocleavage furrows shows an enhanced spread of PH domain but not Bicoid. Modeling analysis suggests that the enhanced spread of the PH domain occurs due to the increased spread of the cytoplasmic population of the PH domain in pseudocleavage furrow mutants. Our analysis of cytoarchitecture interaction in regulating plasma membrane protein distribution and constraining its spread has implications on the mechanisms of spread of various molecules, such as morphogens in syncytial cells.

Transcriptomic and functional genetic evidence for distinct ecophysiological responses across complex life cycle stages

Organisms with complex life cycles demonstrate a remarkable ability to change their phenotypes across development, presumably as an evolutionary adaptation to developmentally variable environments. Developmental variation in environmentally sensitive performance, and thermal sensitivity in particular, has been well documented in holometabolous insects. For example, thermal performance in adults and juvenile stages exhibit little genetic correlation (genetic decoupling) and can evolve independently, resulting in divergent thermal responses. Yet, we understand very little about how this genetic decoupling occurs. We tested the hypothesis that genetic decoupling of thermal physiology is driven by fundamental differences in physiology between life stages, despite a potentially conserved Cellular Stress Response. We used RNAseq to compare transcript expression in response to a cold stressor in Drosophila melanogaster larvae and adults and used RNAi (RNA interference) to test whether knocking down nine target genes differentially affected larval and adult cold tolerance. Transcriptomic responses of whole larvae and adults during and following exposure to -5°C were largely unique both in identity of responding transcripts and in temporal dynamics. Further, we analyzed the tissue-specificity of differentially-expressed transcripts from FlyAtlas 2 data, and concluded that stage-specific differences in transcription were not simply driven by differences in tissue composition. In addition, RNAi of target genes resulted in largely stage-specific and sometimes sex-specific effects on cold tolerance. The combined evidence suggests that thermal physiology is largely stage-specific at the level of gene expression, and thus natural selection may be acting on different loci during the independent thermal adaptation of different life stages.

Population-scale long-read sequencing uncovers transposable elements associated with gene expression variation and adaptive signatures in Drosophila

High quality reference genomes are crucial to understanding genome function, structure and evolution. The availability of reference genomes has allowed us to start inferring the role of genetic variation in biology, disease, and biodiversity conservation. However, analyses across organisms demonstrate that a single reference genome is not enough to capture the global genetic diversity present in populations. In this work, we generate 32 high-quality reference genomes for the well-known model species D. melanogaster and focus on the identification and analysis of transposable element variation as they are the most common type of structural variant. We show that integrating the genetic variation across natural populations from five climatic regions increases the number of detected insertions by 58%. Moreover, 26% to 57% of the insertions identified using long-reads were missed by short-reads methods. We also identify hundreds of transposable elements associated with gene expression variation and new TE variants likely to contribute to adaptive evolution in this species. Our results highlight the importance of incorporating the genetic variation present in natural populations to genomic studies, which is essential if we are to understand how genomes function and evolve.

Metabolic Processes Are Differentially Regulated During Wild-Type and Attenuated Dengue Virus Infection in Aedes aegypti

Successful completion of the dengue virus (DENV) life cycle in its mosquito vectors is important for efficient human-mosquito-human cycle of transmission, but the virus-mosquito interactions that underpin this critical event are poorly defined. To understand the virus-host interactions that determine viral infection by Aedes aegypti, the principal DENV vector, the authors compared transcriptomic changes in the head/thorax of the mosquito after intrathoracic infection with the wild-type DENV2 16681 strain and its attenuated derivative, PDK53. Using high-throughput RNA-sequencing, the authors identified 1,629 differentially expressed genes (DEGs) during 16681 infection, compared with only 22 DEGs identified during PDK53 infection, indicating that 16681 infection triggers a more robust host transcriptomic response compared with PDK53 infection. The authors further found that 16681 infection, but not PDK53 infection, altered metabolism in these heads/thoraces. Altogether, our findings reveal differential regulation of metabolic processes during wild-type and attenuated DENV infection, and suggest the need for future work to study the role of metabolic processes in determining DENV infection and replication in its mosquito vectors.

Endoplasmic reticulum-translocation is essential for APOL1 cellular toxicity

Two variants at the APOL1 gene, encoding apolipoprotein L1, account for more than 70% of the increased risk for chronic kidney disease in individuals of African ancestry. While the initiating event for APOL1 risk variant cell injury remains to be clarified, we explored the possibility of blocking APOL1 toxicity at a more upstream level. We demonstrate that deletion of the first six amino acids of exon 4 abrogates APOL1 cytotoxicity by impairing APOL1 translocation to the lumen of ER and splicing of the signal peptide. Likewise, in orthologous systems, APOL1 lethality was partially abrogated in yeast strains and flies with reduced dosage of genes encoding ER translocon proteins. An inhibitor of ER to Golgi trafficking reduced lethality as well. We suggest that targeting the MSALFL sequence or exon 4 skipping may serve as potential therapeutic approaches to mitigate the risk of CKD caused by APOL1 renal risk variants.

A Random-Screening Approach to Identify RNAi Targets for the Control of Western Corn Rootworm (Diabrotica. virgifera virgifera Le Conte)

Identification of active target genes in bioassay screening is the first important step for application of RNA interference (RNAi) for pest control. Here, we describe the methodology for performing high-throughput RNAi target screening against important agriculture pest, Western corn rootworm in 96-well microplate. Two approaches are presented to identify active targets from random-cDNA library or testing a certain group of specific targets via in silico sequence analysis. Methods of PCR primer design, DNA template preparation, and dsRNA production described here can be applied for other pests.

Decoding gene regulation in the fly brain

The Drosophila brain is a frequently used model in neuroscience. Single-cell transcriptome analysis^1-6, three-dimensional morphological classification⁷ and electron microscopy mapping of the connectome^8,9 have revealed an immense diversity of neuronal and glial cell types that underlie an array of functional and behavioural traits in the fly. The identities of these cell types are controlled by gene regulatory networks (GRNs), involving combinations of transcription factors that bind to genomic enhancers to regulate their target genes. Here, to characterize GRNs at the cell-type level in the fly brain, we profiled the chromatin accessibility of 240,919 single cells spanning 9 developmental timepoints and integrated these data with single-cell transcriptomes. We identify more than 95,000 regulatory regions that are used in different neuronal cell types, of which 70,000 are linked to developmental trajectories involving neurogenesis, reprogramming and maturation. For 40 cell types, uniquely accessible regions were associated with their expressed transcription factors and downstream target genes through a combination of motif discovery, network inference and deep learning, creating enhancer GRNs. The enhancer architectures revealed by DeepFlyBrain lead to a better understanding of neuronal regulatory diversity and can be used to design genetic driver lines for cell types at specific timepoints, facilitating their characterization and manipulation.

Glial immune-related pathways mediate effects of closed head traumatic brain injury on behavior and lethality in Drosophila

In traumatic brain injury (TBI), the initial injury phase is followed by a secondary phase that contributes to neurodegeneration, yet the mechanisms leading to neuropathology in vivo remain to be elucidated. To address this question, we developed a Drosophila head-specific model for TBI termed Drosophila Closed Head Injury (dCHI), where well-controlled, nonpenetrating strikes are delivered to the head of unanesthetized flies. This assay recapitulates many TBI phenotypes, including increased mortality, impaired motor control, fragmented sleep, and increased neuronal cell death. TBI results in significant changes in the transcriptome, including up-regulation of genes encoding antimicrobial peptides (AMPs). To test the in vivo functional role of these changes, we examined TBI-dependent behavior and lethality in mutants of the master immune regulator NF-κB, important for AMP induction, and found that while sleep and motor function effects were reduced, lethality effects were enhanced. Similarly, loss of most AMP classes also renders flies susceptible to lethal TBI effects. These studies validate a new Drosophila TBI model and identify immune pathways as in vivo mediators of TBI effects.

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

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

A connectome is not enough - what is still needed to understand the brain of Drosophila?

Understanding the structure and operation of any nervous system has been a subject of research for well over a century. A near-term opportunity in this quest is to understand the brain of a model species, the fruit fly Drosophila melanogaster. This is an enticing target given its relatively small size (roughly 200,000 neurons), coupled with the behavioral richness that this brain supports, and the wide variety of techniques now available to study both brain and behavior. It is clear that within a few years we will possess a connectome for D. melanogaster: an electron-microscopy-level description of all neurons and their chemical synaptic connections. Given what we will soon have, what we already know and the research that is currently underway, what more do we need to know to enable us to understand the fly's brain? Here, we itemize the data we will need to obtain, collate and organize in order to build an integrated model of the brain of D. melanogaster.

Dia- and Rok-dependent enrichment of capping proteins in a cortical region

Rho signaling with its major targets the formin Dia, Rho kinase (Rok) and non-muscle myosin II (MyoII, encoded by zip in flies) control turnover, amount and contractility of actomyosin. Much less investigated has been a potential function for the distribution of F-actin plus and minus ends. In syncytial Drosophila embryos, Rho1 signaling is high between actin caps, i.e. the cortical intercap region. Capping protein binds to free plus ends of F-actin to prevent elongation of the filament. Capping protein has served as a marker to visualize the distribution of F-actin plus ends in cells and in vitro. In the present study, we probed the distribution of plus ends with capping protein in syncytial Drosophila embryos. We found that capping proteins are specifically enriched in the intercap region similar to Dia and MyoII but distinct from overall F-actin. The intercap enrichment of Capping protein was impaired in dia mutants and embryos, in which Rok and MyoII activation was inhibited. Our observations reveal that Dia and Rok-MyoII control Capping protein enrichment and support a model that Dia and Rok-MyoII control the organization of cortical actin cytoskeleton downstream of Rho1 signaling.

This article has an associated First Person interview with the first authors of the paper.

Summary: Plus ends of actin filaments are enriched at cortical regions rich in Rho signaling in syncytial Drosophila embryos depending on the actin regulator Dia and Rho kinase.

The Drosophila Baramicin polypeptide gene protects against fungal infection

The fruit fly Drosophila melanogaster combats microbial infection by producing a battery of effector peptides that are secreted into the haemolymph. Technical difficulties prevented the investigation of these short effector genes until the recent advent of the CRISPR/CAS era. As a consequence, many putative immune effectors remain to be formally described, and exactly how each of these effectors contribute to survival is not well characterized. Here we describe a novel Drosophila antifungal peptide gene that we name Baramicin A. We show that BaraA encodes a precursor protein cleaved into multiple peptides via furin cleavage sites. BaraA is strongly immune-induced in the fat body downstream of the Toll pathway, but also exhibits expression in other tissues. Importantly, we show that flies lacking BaraA are viable but susceptible to the entomopathogenic fungus Beauveria bassiana. Consistent with BaraA being directly antimicrobial, overexpression of BaraA promotes resistance to fungi and the IM10-like peptides produced by BaraA synergistically inhibit growth of fungi in vitro when combined with a membrane-disrupting antifungal. Surprisingly, BaraA mutant males but not females display an erect wing phenotype upon infection. Here, we characterize a new antifungal immune effector downstream of Toll signalling, and show it is a key contributor to the Drosophila antimicrobial response.

Finding differentially expressed sRNA-Seq regions with srnadiff

Small RNAs (sRNAs) encompass a great variety of molecules of different kinds, such as microRNAs, small interfering RNAs, Piwi-associated RNA, among others. These sRNAs have a wide range of activities, which include gene regulation, protection against virus, transposable element silencing, and have been identified as a key actor in determining the development of the cell. Small RNA sequencing is thus routinely used to assess the expression of the diversity of sRNAs, usually in the context of differentially expression, where two conditions are compared. Tools that detect differentially expressed microRNAs are numerous, because microRNAs are well documented, and the associated genes are well defined. However, tools are lacking to detect other types of sRNAs, which are less studied, and whose precursor RNA is not well characterized. We present here a new method, called srnadiff, which finds all kinds of differentially expressed sRNAs. To the extent of our knowledge, srnadiff is the first tool that detects differentially expressed sRNAs without the use of external information, such as genomic annotation or additional sequences of sRNAs.

Sexual antagonism, temporally fluctuating selection, and variable dominance affect a regulatory polymorphism in Drosophila melanogaster

Understanding how genetic variation is maintained within species is a major goal of evolutionary genetics that can shed light on the preservation of biodiversity. Here, we examined the maintenance of a regulatory single-nucleotide polymorphism (SNP) of the X-linked Drosophila melanogaster gene fezzik. The derived variant at this site is at intermediate frequency in many worldwide populations but absent in populations from the ancestral species range in sub-Saharan Africa. We collected and genotyped wild-caught individuals from a single European population biannually over a period of 5 years, which revealed an overall difference in allele frequency between the sexes and a consistent change in allele frequency across seasons in females but not in males. Modeling based on the observed allele and genotype frequencies suggested that both sexually antagonistic and temporally fluctuating selection may help maintain variation at this site. The derived variant is predicted to be female-beneficial and mostly recessive; however, there was uncertainty surrounding our dominance estimates and long-term modeling projections suggest that it is more likely to be dominant. By examining gene expression phenotypes, we found that phenotypic dominance was variable and dependent upon developmental stage and genetic background, suggesting that dominance may be variable at this locus. We further determined that fezzik expression and genotype are associated with starvation resistance in a sex-dependent manner, suggesting a potential phenotypic target of selection. By characterizing the mechanisms of selection acting on this SNP, our results improve our understanding of how selection maintains genetic and phenotypic variation in natural populations.

The H3.3K27M oncohistone antagonizes reprogramming in Drosophila

Development proceeds by the activation of genes by transcription factors and the inactivation of others by chromatin-mediated gene silencing. In certain cases development can be reversed or redirected by mis-expression of master regulator transcription factors. This must involve the activation of previously silenced genes, and such developmental aberrations are thought to underlie a variety of cancers. Here, we express the wing-specific Vestigial master regulator to reprogram the developing eye, and test the role of silencing in reprogramming using an H3.3K27M oncohistone mutation that dominantly inhibits histone H3K27 trimethylation. We find that production of the oncohistone blocks eye-to-wing reprogramming. CUT&Tag chromatin profiling of mutant tissues shows that H3K27me3 of domains is generally reduced upon oncohistone production, suggesting that a previous developmental program must be silenced for effective transformation. Strikingly, Vg and H3.3K27M synergize to stimulate overgrowth of eye tissue, a phenotype that resembles that of mutations in Polycomb silencing components. Transcriptome profiling of elongating RNA Polymerase II implicates the mis-regulation of signaling factors in overgrowth. Our results demonstrate that growth dysregulation can result from the simple combination of crippled silencing and transcription factor mis-expression, an effect that may explain the origins of oncohistone-bearing cancers.

The differentiation of cell fates in multicellular organisms requires that certain genes be activated, and genes for alternative cell fates are repressed by chromatin silencing. Specific histone mutations that cripple silencing have been found associated with brain cancers in human patients, and these cancers may originate from instability of cell fates. We tested this idea by expressing a wing specification factor in the Drosophila eye to reprogram cell fates and create winged eyes. To test if defects in chromatin silencing increased cell reprogramming, we simultaneously expressed a crippling mutant histone. Contrary to expectations, we found that wing-to-eye reprogramming no longer occurs and instead the eye overgrows, a phenotype reminiscent of the cancers where the histone mutation was first identified. We suggest that reprogramming requires chromatin silencing of the previous developmental program, and that blocking reprogramming can uncouple growth-promoting effects from developmental one of tissue specification transcription factors.

Evolutionary dynamics of sex-biased genes expressed in cricket brains and gonads

Sex-biased gene expression, particularly sex-biased expression in the gonad, has been linked to rates of protein sequence evolution (nonsynonymous to synonymous substitutions, dN/dS) in animals. However, in insects, sex-biased expression studies remain centred on a few holometabolous species. Moreover, other major tissue types such as the brain remain underexplored. Here, we studied sex-biased gene expression and protein evolution in a hemimetabolous insect, the cricket Gryllus bimaculatus. We generated novel male and female RNA-seq data for two sexual tissue types, the gonad and somatic reproductive system, and for two core components of the nervous system, the brain and ventral nerve cord. From a genome-wide analysis, we report several core findings. Firstly, testis-biased genes had accelerated evolution, as compared to ovary-biased and unbiased genes, which was associated with positive selection events. Secondly, although sex-biased brain genes were much less common than for the gonad, they exhibited a striking tendency for rapid protein sequence evolution, an effect that was stronger for the female than male brain. Further, some sex-biased brain genes were linked to sexual functions and mating behaviours, which we suggest may have accelerated their evolution via sexual selection. Thirdly, a tendency for narrow cross-tissue expression breadth, suggesting low pleiotropy, was observed for sex-biased brain genes, suggesting relaxed purifying selection, which we speculate may allow enhanced freedom to evolve adaptive protein functional changes. The findings of rapid evolution of testis-biased genes and male and female-biased brain genes are discussed with respect to pleiotropy, positive selection and the mating biology of this cricket.

Databases for Protein-Protein Interactions

Protein-protein interaction networks have a crucial role in biological processes. Proteins perform multiple functions in forming physical and functional interactions in cellular systems. Information concerning an enormous number of protein interactions in a wide range of species has accumulated and has been integrated into various resources for molecular biology and systems biology. This chapter provides a review of the representative databases and the major computational methods used for protein-protein interactions.

Clinal and seasonal changes are correlated in Drosophila melanogaster natural populations

Spatial and seasonal variations in the environment are ubiquitous. Environmental heterogeneity can affect natural populations and lead to covariation between environment and allele frequencies. Drosophila melanogaster is known to harbor polymorphisms that change both with latitude and seasons. Identifying the role of selection in driving these changes is not trivial, because nonadaptive processes can cause similar patterns. Given the environment changes in similar ways across seasons and along the latitudinal gradient, one promising approach may be to look for parallelism between clinal and seasonal changes. Here, we test whether there is a genome-wide correlation between clinal and seasonal changes, and whether the pattern is consistent with selection. Allele frequency estimates were obtained from pooled samples from seven different locations along the east coast of the United States, and across seasons within Pennsylvania. We show that there is a genome-wide correlation between clinal and seasonal variations, which cannot be explained by linked selection alone. This pattern is stronger in genomic regions with higher functional content, consistent with natural selection. We derive a way to biologically interpret these correlations and show that around 3.7% of the common, autosomal variants could be under parallel seasonal and spatial selection. Our results highlight the contribution of natural selection in driving fluctuations in allele frequencies in natural fly populations and point to a shared genomic basis to climate adaptation that happens over space and time in D. melanogaster.

The genomic basis of evolutionary differentiation among honey bees

In contrast to the western honey bee, Apis mellifera, other honey bee species have been largely neglected despite their importance and diversity. The genetic basis of the evolutionary diversification of honey bees remains largely unknown. Here, we provide a genome-wide comparison of three honey bee species, each representing one of the three subgenera of honey bees, namely the dwarf (Apis florea), giant (A. dorsata), and cavity-nesting (A. mellifera) honey bees with bumblebees as an outgroup. Our analyses resolve the phylogeny of honey bees with the dwarf honey bees diverging first. We find that evolution of increased eusocial complexity in Apis proceeds via increases in the complexity of gene regulation, which is in agreement with previous studies. However, this process seems to be related to pathways other than transcriptional control. Positive selection patterns across Apis reveal a trade-off between maintaining genome stability and generating genetic diversity, with a rapidly evolving piRNA pathway leading to genomes depleted of transposable elements, and a rapidly evolving DNA repair pathway associated with high recombination rates in all Apis species. Diversification within Apis is accompanied by positive selection in several genes whose putative functions present candidate mechanisms for lineage-specific adaptations, such as migration, immunity, and nesting behavior.

Automated generation of context-specific gene regulatory networks with a weighted approach in Drosophila melanogaster

The regulation of gene expression is a key factor in the development and maintenance of life in all organisms. Even so, little is known at whole genome scale for most genes and contexts. We propose a method, Tool for Weighted Epigenomic Networks in Drosophila melanogaster (Fly T-WEoN), to generate context-specific gene regulatory networks starting from a reference network that contains all known gene regulations in the fly. Unlikely regulations are removed by applying a series of knowledge-based filters. Each of these filters is implemented as an independent module that considers a type of experimental evidence, including DNA methylation, chromatin accessibility, histone modifications and gene expression. Fly T-WEoN is based on heuristic rules that reflect current knowledge on gene regulation in D. melanogaster obtained from the literature. Experimental data files can be generated with several standard procedures and used solely when and if available. Fly T-WEoN is available as a Cytoscape application that permits integration with other tools and facilitates downstream network analysis. In this work, we first demonstrate the reliability of our method to then provide a relevant application case of our tool: early development of D. melanogaster. Fly T-WEoN together with its step-by-step guide is available at https://weon.readthedocs.io.

Gene expression and alternative splicing dynamics are perturbed in female head transcriptomes following heterospecific copulation

BACKGROUND

Despite the growing interest in the female side of copulatory interactions, the roles played by differential expression and alternative splicing mechanisms of pre-RNA on tissues outside of the reproductive tract have remained largely unknown. Here we addressed these questions in the context of con- vs heterospecific matings between Drosophila mojavensis and its sister species, D. arizonae. We analyzed transcriptional responses in female heads using an integrated investigation of genome-wide patterns of gene expression, including differential expression (DE), alternative splicing (AS) and intron retention (IR).

RESULTS

Our results indicated that early transcriptional responses were largely congruent between con- and heterospecific matings but are substantially perturbed over time. Conspecific matings induced functional pathways related to amino acid balance previously associated with the brain's physiology and female postmating behavior. Heterospecific matings often failed to activate regulation of some of these genes and induced expression of additional genes when compared with those of conspecifically-mated females. These mechanisms showed functional specializations with DE genes mostly linked to pathways of proteolysis and nutrient homeostasis, while AS genes were more related to photoreception and muscle assembly pathways. IR seems to play a more general role in DE regulation during the female postmating response.

CONCLUSIONS

We provide evidence showing that AS genes substantially perturbed by heterospecific matings in female heads evolve at slower evolutionary rates than the genome background. However, DE genes evolve at evolutionary rates similar, or even higher, than those of male reproductive genes, which highlights their potential role in sexual selection and the evolution of reproductive barriers.

False and true positives in arthropod thermal adaptation candidate gene lists

Genome-wide studies are prone to false positives due to inherently low priors and statistical power. One approach to ameliorate this problem is to seek validation of reported candidate genes across independent studies: genes with repeatedly discovered effects are less likely to be false positives. Inversely, genes reported only as many times as expected by chance alone, while possibly representing novel discoveries, are also more likely to be false positives. We show that, across over 30 genome-wide studies that reported Drosophila and Daphnia genes with possible roles in thermal adaptation, the combined lists of candidate genes and orthologous groups are rapidly approaching the total number of genes and orthologous groups in the respective genomes. This is consistent with the expectation of high frequency of false positives. The majority of these spurious candidates have been identified by one or a few studies, as expected by chance alone. In contrast, a noticeable minority of genes have been identified by numerous studies with the probabilities of such discoveries occurring by chance alone being exceedingly small. For this subset of genes, different studies are in agreement with each other despite differences in the ecological settings, genomic tools and methodology, and reporting thresholds. We provide a reference set of presumed true positives among Drosophila candidate genes and orthologous groups involved in response to changes in temperature, suitable for cross-validation purposes. Despite this approach being prone to false negatives, this list of presumed true positives includes several hundred genes, consistent with the "omnigenic" concept of genetic architecture of complex traits.

Proteomics of protein trafficking by in vivo tissue-specific labeling

Conventional approaches to identify secreted factors that regulate homeostasis are limited in their abilities to identify the tissues/cells of origin and destination. We established a platform to identify secreted protein trafficking between organs using an engineered biotin ligase (BirA*G3) that biotinylates, promiscuously, proteins in a subcellular compartment of one tissue. Subsequently, biotinylated proteins are affinity-enriched and identified from distal organs using quantitative mass spectrometry. Applying this approach in Drosophila, we identify 51 muscle-secreted proteins from heads and 269 fat body-secreted proteins from legs/muscles, including CG2145 (human ortholog ENDOU) that binds directly to muscles and promotes activity. In addition, in mice, we identify 291 serum proteins secreted from conditional BirA*G3 embryo stem cell-derived teratomas, including low-abundance proteins with hormonal properties. Our findings indicate that the communication network of secreted proteins is vast. This approach has broad potential across different model systems to identify cell-specific secretomes and mediators of interorgan communication in health or disease.

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.