FlyBase 2.0: the next generation

A shift in chromatin binding of phosphorylated p38 precedes transcriptional changes upon oxidative stress

P38 mitogen-activated protein kinases are key in the regulation of the cellular response to stressors. P38 is known to regulate transcription, mRNA processing, stability, and translation. The transcriptional changes mediated by phosphorylated p38 (P-p38) in response to extracellular stimuli have been thoroughly analyzed in many tissues and organisms. However, the genomic localization of chromatin-associated P-p38 remains poorly understood. Here, we analyze the chromatin binding of activated P-p38 and its role in the response to reactive oxygen species (ROS) in Drosophila S2 cells. We found that P-p38 is already bound to chromatin in basal conditions. After ROS exposure, chromatin-associated P-p38 relocates towards genes involved in the recovery process. Our findings highlight the role of P-p38 dynamic chromatin binding in orchestrating gene expression responses to oxidative stress.

Validation of heat-inducible Ixodes scapularis HSP70 and tick-specific 3xP3 promoters in ISE6 cells

Ixodes scapularis is an important vector of many pathogens, including the causative agent of Lyme disease. The gene function studies in I. scapularis and other ticks are hampered by the lack of genetic tools, including an inducible promoter for temporal control over transgene-encoding protein or double-stranded RNA. We characterized an intergenic sequence upstream of a heat shock protein 70 (HSP70) gene that can drive Renilla luciferase and mCherry expression in the I. scapularis cell line ISE6 (IsHSP70). In another construct, we replaced the Drosophila melanogaster minimal HSP70 promoter of the 3xP3 promoter with a minimal portion of IsHSP70 promoter and generated an I. scapularis-specific 3xP3 (Is3xP3) promoter. Both IsHSP70 and Is3xP3 have a heat-inducible expression of mCherry fluorescence in ISE6 cells with an approximately 10-fold increase in the percentage of fluorescent cells upon 2 h heat shock. These promoters described will be valuable tools for gene function studies.

Actin Cytoskeleton and Integrin Components Are Interdependent for Slit Diaphragm Maintenance in Drosophila Nephrocytes

In nephrotic syndrome, the podocyte filtration structures are damaged in a process called foot process effacement. This is mediated by the actin cytoskeleton; however, which actins are involved and how they interact with other filtration components, like the basement membrane, remains poorly understood. Here, we used the well-established Drosophila pericardial nephrocyte-the equivalent of podocytes in flies-knockdown models (RNAi) to study the interplay of the actin cytoskeleton (Act5C, Act57B, Act42A, and Act87E), alpha- and beta-integrin (basement membrane), and the slit diaphragm (Sns and Pyd). Knockdown of an actin gene led to variations of formation of actin stress fibers, the internalization of Sns, and a disrupted slit diaphragm cortical pattern. Notably, deficiency of Act5C, which resulted in complete absence of nephrocytes, could be partially mitigated by overexpressing Act42A or Act87E, suggesting at least partial functional redundancy. Integrin localized near the actin cytoskeleton as well as slit diaphragm components, but when the nephrocyte cytoskeleton or slit diaphragm was disrupted, this switched to colocalization, both at the surface and internalized in aggregates. Altogether, the data show that the interdependence of the slit diaphragm, actin cytoskeleton, and integrins is key to the structure and function of the Drosophila nephrocyte.

Pervasive tissue-, genetic background-, and allele-specific gene expression effects in Drosophila melanogaster

The pervasiveness of gene expression variation and its contribution to phenotypic variation and evolution is well known. This gene expression variation is context dependent, with differences in regulatory architecture often associated with intrinsic and environmental factors, and is modulated by regulatory elements that can act in cis (linked) or in trans (unlinked) relative to the genes they affect. So far, little is known about how this genetic variation affects the evolution of regulatory architecture among closely related tissues during population divergence. To address this question, we analyzed gene expression in the midgut, hindgut, and Malpighian tubule as well as microbiome composition in the two gut tissues in four Drosophila melanogaster strains and their F1 hybrids from two divergent populations: one from the derived, European range and one from the ancestral, African range. In both the transcriptome and microbiome data, we detected extensive tissue- and genetic background-specific effects, including effects of genetic background on overall tissue specificity. Tissue-specific effects were typically stronger than genetic background-specific effects, although the two gut tissues were not more similar to each other than to the Malpighian tubules. An examination of allele specific expression revealed that, while both cis and trans effects were more tissue-specific in genes expressed differentially between populations than genes with conserved expression, trans effects were more tissue-specific than cis effects. Despite there being highly variable regulatory architecture, this observation was robust across tissues and genetic backgrounds, suggesting that the expression of trans variation can be spatially fine-tuned as well as or better than cis variation during population divergence and yielding new insights into cis and trans regulatory evolution.

Composition and abundance of midgut plasma membrane proteins in two major hemipteran vectors of plant viruses, Bemisia tabaci and Myzus persicae

Multiple species within the order Hemiptera cause severe agricultural losses on a global scale. Aphids and whiteflies are of particular importance due to their role as vectors for hundreds of plant viruses, many of which enter the insect via the gut. To facilitate the identification of novel targets for disruption of plant virus transmission, we compared the relative abundance and composition of the gut plasma membrane proteomes of adult Bemisia tabaci (Hemiptera: Aleyrodidae) and Myzus persicae (Hemiptera: Aphididae), representing the first study comparing the gut plasma membrane proteomes of two different insect species. Brush border membrane vesicles were prepared from dissected guts, and proteins extracted, identified and quantified from triplicate samples via timsTOF mass spectrometry. A total of 1699 B. tabaci and 1175 M. persicae proteins were identified. Following bioinformatics analysis and manual curation, 151 B. tabaci and 115 M. persicae proteins were predicted to localize to the plasma membrane of the gut microvilli. These proteins were further categorized based on molecular function and biological process according to Gene Ontology terms. The most abundant gut plasma membrane proteins were identified. The ten plasma membrane proteins that differed in abundance between the two insect species were associated with the terms "protein binding" and "viral processes." In addition to providing insight into the gut physiology of hemipteran insects, these gut plasma membrane proteomes provide context for appropriate identification of plant virus receptors based on a combination of bioinformatic prediction and protein localization on the surface of the insect gut.

Characterization of the chemoreceptor repertoire of a highly specialized fly with comparisons to other Drosophila species

To explore the diversity of scenarios in nature, animals have evolved tools to interact with different environmental conditions. Chemoreceptors are an important interface component and among them, olfactory receptors (ORs) and gustatory receptors (GRs) can be used to find food and detect healthy resources. Drosophila is a model organism in many scientific fields, in part due to the diversity of species and niches they occupy. The contrast between generalists and specialists Drosophila species provides an important model for studying the evolution of chemoreception. Here, we compare the repertoire of chemoreceptors of different species of Drosophila with that of D. incompta, a highly specialized species whose ecology is restricted to Cestrum flowers, after reporting the preferences of D. incompta to the odor of Cestrum flowers in olfactory tests. We found evidence that the chemoreceptor repertoire in D. incompta is smaller than that presented by species in the Sophophora subgenus. Similar patterns were found in other non-Sophophora species, suggesting the presence of underlying phylogenetic trends. Nevertheless, we also found autapomorphic gene losses and detected some genes that appear to be under positive selection in D. incompta, suggesting that the specific lifestyle of these flies may have shaped the evolution of individual genes in each of these gene families.

Shifted levels of sleep and activity under darkness as mechanisms underlying ectoparasite resistance

Parasites harm host fitness and are pervasive agents of natural selection to evolve host defense strategies Host defensive traits in natural populations typically show genetic variation, which may be maintained when parasite resistance imposes fitness costs on the host in the absence of parasites. Previously we demonstrated significant evolutionary responses to artificial selection for increasing behavioral immunity to Gamasodes queenslandicus mites in replicate lines of Drosophila melanogaster. Here, we report transcriptional shifts in metabolic processes between selected and control fly lines based on RNA-seq analyses. We also show decreased starvation resistance and increased use of nutrient reserves in flies from mite-resistant lines. Additionally, mite-resistant lines exhibited increased behavioral activity, such as, reduced sleep and elevated oxygen consumption under conditions of darkness. The link between resistance and sleep was confirmed in an independent panel of D. melanogaster genetic lines exhibiting variable sleep durations, showing a positive correlation between mite resistance and reduced sleep. Experimentally restraining the activity of artificially selected mite-resistant flies during exposure to parasites under dark conditions reduced their resistance advantage relative to control flies. The results suggest that ectoparasite resistance in this system involves increased dark-condition activity and metabolic gene expression at the expense of nutrient reserves and starvation resistance.

An important role for triglyceride in regulating spermatogenesis

Drosophila is a powerful model to study how lipids affect spermatogenesis. Yet, the contribution of neutral lipids, a major lipid group which resides in organelles called lipid droplets (LD), to sperm development is largely unknown. Emerging evidence suggests LD are present in the testis and that loss of neutral lipid- and LD-associated genes causes subfertility; however, key regulators of testis neutral lipids and LD remain unclear. Here, we show LD are present in early-stage somatic and germline cells within the Drosophila testis. We identified a role for triglyceride lipase brummer (bmm) in regulating testis LD, and found that whole-body loss of bmm leads to defects in sperm development. Importantly, these represent cell-autonomous roles for bmm in regulating testis LD and spermatogenesis. Because lipidomic analysis of bmm mutants revealed excess triglyceride accumulation, and spermatogenic defects in bmm mutants were rescued by genetically blocking triglyceride synthesis, our data suggest that bmm-mediated regulation of triglyceride influences sperm development. This identifies triglyceride as an important neutral lipid that contributes to Drosophila sperm development, and reveals a key role for bmm in regulating testis triglyceride levels during spermatogenesis.

Large-language models facilitate discovery of the molecular signatures regulating sleep and activity

Sleep, locomotor and social activities are essential animal behaviors, but their reciprocal relationships and underlying mechanisms remain poorly understood. Here, we elicit information from a cutting-edge large-language model (LLM), generative pre-trained transformer (GPT) 3.5, which interprets 10.2-13.8% of Drosophila genes known to regulate the 3 behaviors. We develop an instrument for simultaneous video tracking of multiple moving objects, and conduct a genome-wide screen. We have identified 758 fly genes that regulate sleep and activities, including mre11 which regulates sleep only in the presence of conspecifics, and NELF-B which regulates sleep regardless of whether conspecifics are present. Based on LLM-reasoning, an educated signal web is modeled for understanding of potential relationships between its components, presenting comprehensive molecular signatures that control sleep, locomotor and social activities. This LLM-aided strategy may also be helpful for addressing other complex scientific questions.

The alpha-ketoacid dehydrogenase complexes of Drosophila melanogaster

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

Recognition of nonself is necessary to activate Drosophila's immune response against an insect parasite

BACKGROUND

Innate immune responses can be activated by pathogen-associated molecular patterns (PAMPs), danger signals released by damaged tissues, or the absence of self-molecules that inhibit immunity. As PAMPs are typically conserved across broad groups of pathogens but absent from the host, it is unclear whether they allow hosts to recognize parasites that are phylogenetically similar to themselves, such as parasitoid wasps infecting insects.

RESULTS

Parasitoids must penetrate the cuticle of Drosophila larvae to inject their eggs. In line with previous results, we found that the danger signal of wounding triggers the differentiation of specialized immune cells called lamellocytes. However, using oil droplets to mimic infection by a parasitoid wasp egg, we found that this does not activate the melanization response. This aspect of the immune response also requires exposure to parasite molecules. The unidentified factor enhances the transcriptional response in hemocytes and induces a specific response in the fat body.

CONCLUSIONS

We conclude that a combination of danger signals and the recognition of nonself molecules is required to activate Drosophila's immune response against parasitic insects.

TRPA5 encodes a thermosensitive ankyrin ion channel receptor in a triatomine insect

As ectotherms, insects need heat-sensitive receptors to monitor environmental temperatures and facilitate thermoregulation. We show that TRPA5, a class of ankyrin transient receptor potential (TRP) channels absent in dipteran genomes, may function as insect heat receptors. In the triatomine bug Rhodnius prolixus (order: Hemiptera), a vector of Chagas disease, the channel RpTRPA5B displays a uniquely high thermosensitivity, with biophysical determinants including a large channel activation enthalpy change (72 kcal/mol), a high temperature coefficient (Q10 = 25), and in vitro temperature-induced currents from 53°C to 68°C (T0.5 = 58.6°C), similar to noxious TRPV receptors in mammals. Monomeric and tetrameric ion channel structure predictions show reliable parallels with fruit fly dTRPA1, with structural uniqueness in ankyrin repeat domains, the channel selectivity filter, and potential TRP functional modulator regions. Overall, the finding of a member of TRPA5 as a temperature-activated receptor illustrates the diversity of insect molecular heat detectors.

Genomes of historical specimens reveal multiple invasions of LTR retrotransposons in Drosophila melanogaster during the 19th century

Transposable element invasions have a profound impact on the evolution of genomes and phenotypes. It is thus an important open question how often such TE invasions occur. To address this question, we utilize the genomes of historical specimens, sampled about 200 y ago. We found that the LTR retrotransposons Blood, Opus, and 412 spread in Drosophila melanogaster in the 19th century. These invasions constitute second waves, as degraded fragments were found for all three TEs. The composition of Opus and 412, but not of Blood, shows a pronounced geographic heterogeneity, likely due to founder effects during the invasions. Finally, we identified species from the Drosophila simulans complex as the likely origin of the TEs. We show that in total, seven TE families invaded D. melanogaster during the last 200y, thereby increasing the genome size by up to 1.2Mbp. We suggest that this high rate of TE invasions was likely triggered by human activity. Based on the analysis of strains and specimens sampled at different times, we provide a detailed timeline of TE invasions, making D. melanogaster the first organism where the invasion history of TEs during the last two centuries could be inferred.

Palmitoylation of proteolipid protein M6 promotes tricellular junction assembly in epithelia of Drosophila

Tricellular junctions (TCJs) seal epithelial cell vertices and are essential for tissue integrity and physiology, but how TCJs are assembled and maintained is poorly understood. In Drosophila, the transmembrane proteins Anakonda (Aka, also known as Bark), Gliotactin (Gli) and M6 organize occluding TCJs. Aka and M6 localize in an interdependent manner to vertices and act jointly to localize Gli, but how these proteins interact to assemble TCJs was not previously known. Here, we show that the proteolipid protein M6 physically interacts with Aka and with itself, and that M6 is palmitoylated on conserved juxta-membrane cysteine residues. This modification promotes vertex localization of M6 and binding to Aka, but not to itself, and becomes essential when TCJ protein levels are reduced. Abolishing M6 palmitoylation leads to delayed localization of M6 and Aka but does not affect the rate of TCJ growth or mobility of M6 or Aka. Our findings suggest that palmitoylation-dependent recruitment of Aka by M6 promotes initiation of TCJ assembly, whereas subsequent TCJ growth relies on different mechanisms that are independent of M6 palmitoylation.

The genome of the blind bee louse fly reveals deep convergences with its social host and illuminates Drosophila origins

Social insects' nests harbor intruders known as inquilines,1 which are usually related to their hosts.2,3 However, distant non-social inquilines may also show convergences with their hosts,4,5 although the underlying genomic changes remain unclear. We analyzed the genome of the wingless and blind bee louse fly Braula coeca, an inquiline kleptoparasite of the western honey bee, Apis mellifera.6,7 Using large phylogenomic data, we confirmed recent accounts that the bee louse fly is a drosophilid8,9 and showed that it had likely evolved from a sap-breeder ancestor associated with honeydew and scale insects' wax. Unlike many parasites, the bee louse fly genome did not show significant erosion or strict reliance on an endosymbiont, likely due to a relatively recent age of inquilinism. However, we observed a horizontal transfer of a transposon and a striking parallel evolution in a set of gene families between the honey bee and the bee louse fly. Convergences included genes potentially involved in metabolism and immunity and the loss of nearly all bitter-tasting gustatory receptors, in agreement with life in a protective nest and a diet of honey, pollen, and beeswax. Vision and odorant receptor genes also exhibited rapid losses. Only genes whose orthologs in the closely related Drosophila melanogaster respond to honey bee pheromone components or floral aroma were retained, whereas the losses included orthologous receptors responsive to the anti-ovarian honey bee queen pheromones. Hence, deep genomic convergences can underlie major phenotypic transitions during the evolution of inquilinism between non-social parasites and their social hosts.

Gene expression variation underlying tissue-specific responses to copper stress in Drosophila melanogaster

Copper is one of a handful of biologically necessary heavy metals that is also a common environmental pollutant. Under normal conditions, copper ions are required for many key physiological processes. However, in excess, copper results in cell and tissue damage ranging in severity from temporary injury to permanent neurological damage. Because of its biological relevance, and because many conserved copper-responsive genes respond to nonessential heavy metal pollutants, copper resistance in Drosophila melanogaster is a useful model system with which to investigate the genetic control of the heavy metal stress response. Because heavy metal toxicity has the potential to differently impact specific tissues, we genetically characterized the control of the gene expression response to copper stress in a tissue-specific manner in this study. We assessed the copper stress response in head and gut tissue of 96 inbred strains from the Drosophila Synthetic Population Resource using a combination of differential expression analysis and expression quantitative trait locus mapping. Differential expression analysis revealed clear patterns of tissue-specific expression. Tissue and treatment specific responses to copper stress were also detected using expression quantitative trait locus mapping. Expression quantitative trait locus associated with MtnA, Mdr49, Mdr50, and Sod3 exhibited both genotype-by-tissue and genotype-by-treatment effects on gene expression under copper stress, illuminating tissue- and treatment-specific patterns of gene expression control. Together, our data build a nuanced description of the roles and interactions between allelic and expression variation in copper-responsive genes, provide valuable insight into the genomic architecture of susceptibility to metal toxicity, and highlight candidate genes for future functional characterization.

The Complex Landscape of Structural Divergence Between the Drosophila pseudoobscura and D. persimilis Genomes

Structural genomic variants are key drivers of phenotypic evolution. They can span hundreds to millions of base pairs and can thus affect large numbers of genetic elements. Although structural variation is quite common within and between species, its characterization depends upon the quality of genome assemblies and the proportion of repetitive elements. Using new high-quality genome assemblies, we report a complex and previously hidden landscape of structural divergence between the genomes of Drosophila persimilis and D. pseudoobscura, two classic species in speciation research, and study the relationships among structural variants, transposable elements, and gene expression divergence. The new assemblies confirm the already known fixed inversion differences between these species. Consistent with previous studies showing higher levels of nucleotide divergence between fixed inversions relative to collinear regions of the genome, we also find a significant overrepresentation of INDELs inside the inversions. We find that transposable elements accumulate in regions with low levels of recombination, and spatial correlation analyses reveal a strong association between transposable elements and structural variants. We also report a strong association between differentially expressed (DE) genes and structural variants and an overrepresentation of DE genes inside the fixed chromosomal inversions that separate this species pair. Interestingly, species-specific structural variants are overrepresented in DE genes involved in neural development, spermatogenesis, and oocyte-to-embryo transition. Overall, our results highlight the association of transposable elements with structural variants and their importance in driving evolutionary divergence.

Spoink, a LTR retrotransposon, invaded D. melanogaster populations in the 1990s

During the last few centuries D. melanogaster populations were invaded by several transposable elements, the most recent of which was thought to be the P-element between 1950 and 1980. Here we describe a novel TE, which we named Spoink, that has invaded D. melanogaster. It is a 5216nt LTR retrotransposon of the Ty3/gypsy superfamily. Relying on strains sampled at different times during the last century we show that Spoink invaded worldwide D. melanogaster populations after the P-element between 1983 and 1993. This invasion was likely triggered by a horizontal transfer from the D. willistoni group, much as the P-element. Spoink is probably silenced by the piRNA pathway in natural populations and about 1/3 of the examined strains have an insertion into a canonical piRNA cluster such as 42AB. Given the degree of genetic investigation of D. melanogaster it is perhaps surprising that Spoink was able to invade unnoticed.

Characterisation of the in-vivo miRNA landscape in Drosophila ribonuclease mutants reveals Pacman-mediated regulation of the highly conserved let-7 cluster during apoptotic processes

The control of gene expression is a fundamental process essential for correct development and to maintain homeostasis. Many post-transcriptional mechanisms exist to maintain the correct levels of each RNA transcript within the cell. Controlled and targeted cytoplasmic RNA degradation is one such mechanism with the 5'-3' exoribonuclease Pacman (XRN1) and the 3'-5' exoribonuclease Dis3L2 playing crucial roles. Loss of function mutations in either Pacman or Dis3L2 have been demonstrated to result in distinct phenotypes, and both have been implicated in human disease. One mechanism by which gene expression is controlled is through the function of miRNAs which have been shown to be crucial for the control of almost all cellular processes. Although the biogenesis and mechanisms of action of miRNAs have been comprehensively studied, the mechanisms regulating their own turnover are not well understood. Here we characterise the miRNA landscape in a natural developing tissue, the Drosophila melanogaster wing imaginal disc, and assess the importance of Pacman and Dis3L2 on the abundance of miRNAs. We reveal a complex landscape of miRNA expression and show that whilst a null mutation in dis3L2 has a minimal effect on the miRNA expression profile, loss of Pacman has a profound effect with a third of all detected miRNAs demonstrating Pacman sensitivity. We also reveal a role for Pacman in regulating the highly conserved let-7 cluster (containing miR-100, let-7 and miR-125) and present a genetic model outlining a positive feedback loop regulated by Pacman which enhances our understanding of the apoptotic phenotype observed in Pacman mutants.

Crop-GPA: an integrated platform of crop gene-phenotype associations

With the increasing availability of large-scale biology data in crop plants, there is an urgent demand for a versatile platform that fully mines and utilizes the data for modern molecular breeding. We present Crop-GPA ( https://crop-gpa.aielab.net ), a comprehensive and functional open-source platform for crop gene-phenotype association data. The current Crop-GPA provides well-curated information on genes, phenotypes, and their associations (GPAs) to researchers through an intuitive interface, dynamic graphical visualizations, and efficient online tools. Two computational tools, GPA-BERT and GPA-GCN, are specifically developed and integrated into Crop-GPA, facilitating the automatic extraction of gene-phenotype associations from bio-crop literature and predicting unknown relations based on known associations. Through usage examples, we demonstrate how our platform enables the exploration of complex correlations between genes and phenotypes in crop plants. In summary, Crop-GPA serves as a valuable multi-functional resource, empowering the crop research community to gain deeper insights into the biological mechanisms of interest.

Pri peptides temporally coordinate transcriptional programs during epidermal differentiation

To achieve a highly differentiated state, cells undergo multiple transcriptional processes whose coordination and timing are not well understood. In Drosophila embryonic epidermal cells, polished-rice (Pri) smORF peptides act as temporal mediators of ecdysone to activate a transcriptional program leading to cell shape remodeling. Here, we show that the ecdysone/Pri axis concomitantly represses the transcription of a large subset of cuticle genes to ensure proper differentiation of the insect exoskeleton. The repression relies on the transcription factor Ken and persists for several days throughout early larval stages, during which a soft cuticle allows larval crawling. The onset of these cuticle genes normally awaits the end of larval stages when the rigid pupal case assembles, and their premature expression triggers abnormal sclerotization of the larval cuticle. These results uncovered a temporal switch to set up distinct structures of cuticles adapted to the animal lifestyle and which might be involved in the evolutionary history of insects.

Bone morphogenetic protein signaling: the pathway and its regulation

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

Unraveling the G protein-coupled receptor superfamily in aphids: Contractions and duplications linked to phloem feeding

The G Protein-Coupled Receptor (GPCR) superfamily is the largest and most diverse transmembrane receptor family, playing crucial roles in regulating various physiological processes. As one of the most destructive pests, aphids have been subject to previous studies, which revealed fewer GPCR superfamily members in Acyrthosiphon pisum and Aphis gossypii and the loss of multiple neuropeptide GPCRs. To elucidate the contraction patterns and evolutionary features of the aphid GPCR superfamily, we identified 97, 105, and 95 GPCR genes in Rhopalosiphum maidis, A. pisum, and A. gossypii, respectively. Comparative analysis and phylogenetic investigations with other hemipteran insects revealed a contracted GPCR superfamily in aphids. This contraction mainly occurred in biogenic amine receptors, GABA-B-R, and fz families, and several neuropeptide receptors such as ACPR, CrzR, and PTHR were completely lost. This phenomenon may be related to the parasitic nature of aphids. Additionally, several GPCRs associated with aphid feeding and water balance underwent duplication, including Lkr, NPFR, CCHa1-R, and DH-R, Type A LGRs, but the SK/CCKLR that inhibits feeding was completely lost, indicating changes in feeding genes that underpin the aphid's prolonged phloem feeding behavior. Furthermore, we observed fine-tuning in opsins, with reduced long-wavelength opsins and additional duplications of short-wavelength opsin, likely associated with daytime activity. Lastly, we found variations in the number of mthl genes in aphids. In conclusion, our investigation sheds light on the GPCR superfamily in aphids, revealing its association with diet lifestyle and laying the foundation for understanding and developing control strategies for the aphid GPCR superfamily.

IntroUNET: Identifying introgressed alleles via semantic segmentation

A growing body of evidence suggests that gene flow between closely related species is a widespread phenomenon. Alleles that introgress from one species into a close relative are typically neutral or deleterious, but sometimes confer a significant fitness advantage. Given the potential relevance to speciation and adaptation, numerous methods have therefore been devised to identify regions of the genome that have experienced introgression. Recently, supervised machine learning approaches have been shown to be highly effective for detecting introgression. One especially promising approach is to treat population genetic inference as an image classification problem, and feed an image representation of a population genetic alignment as input to a deep neural network that distinguishes among evolutionary models (i.e. introgression or no introgression). However, if we wish to investigate the full extent and fitness effects of introgression, merely identifying genomic regions in a population genetic alignment that harbor introgressed loci is insufficient-ideally we would be able to infer precisely which individuals have introgressed material and at which positions in the genome. Here we adapt a deep learning algorithm for semantic segmentation, the task of correctly identifying the type of object to which each individual pixel in an image belongs, to the task of identifying introgressed alleles. Our trained neural network is thus able to infer, for each individual in a two-population alignment, which of those individual's alleles were introgressed from the other population. We use simulated data to show that this approach is highly accurate, and that it can be readily extended to identify alleles that are introgressed from an unsampled "ghost" population, performing comparably to a supervised learning method tailored specifically to that task. Finally, we apply this method to data from Drosophila, showing that it is able to accurately recover introgressed haplotypes from real data. This analysis reveals that introgressed alleles are typically confined to lower frequencies within genic regions, suggestive of purifying selection, but are found at much higher frequencies in a region previously shown to be affected by adaptive introgression. Our method's success in recovering introgressed haplotypes in challenging real-world scenarios underscores the utility of deep learning approaches for making richer evolutionary inferences from genomic data.

Comparative genomics reveals evolutionary drivers of the dietary shift in Hemiptera

Hemiptera insects exhibit a close relationship to plants and demonstrate a diverse range of dietary preferences, encompassing phytophagy as the predominant feeding habit while a minority engages in carnivorous or haematophagous behaviour. To counteract the challenges posed by phytophagous insects, plants have developed an array of toxic compounds, causing significant evolutionary selection pressure on these insects. In this study, we employed a comparative genomics approach to analyse the expansion and contraction of gene families specific to phytophagous insect lineages, along with their adaptive evolutionary traits, utilising representative species from the Hemiptera order. Our investigation revealed substantial expansions of gene families within the phytophagous lineages, especially in the Pentatomomorpha branch represented by Oncopeltus fasciatus and Riptortus pedestris. Notably, these expansions of gene families encoding enzymes are potentially involved in hemipteran-plant interactions. Moreover, the adaptive evolutionary analysis of these lineages revealed a higher prevalence of adaptively evolved genes in the Pentatomomorpha branch. The observed branch-specific gene expansions and adaptive evolution likely contribute significantly to the diversification of species within Hemiptera. These results help enhance our understanding of the genomic characteristics of the evolution of different feeding habits in hemipteran insects.

Regulatory evolution tuning pigmentation intensity quantitatively in Drosophila

Quantitative variation in attributes such as color, texture, or stiffness dominates morphological diversification. It results from combinations of alleles at many Mendelian loci. Here, we identify an additional source of quantitative variation among species, continuous evolution in a gene regulatory region. Specifically, we examined the modulation of wing pigmentation in a group of fly species and showed that inter-species variation correlated with the quantitative expression of the pigmentation gene yellow. This variation results from an enhancer of yellow determining darkness through species-specific activity. We mapped the divergent activities between two sister species and found the changes to be broadly distributed along the enhancer. Our results demonstrate that enhancers can act as dials fueling quantitative morphological diversification by modulating trait properties.

Recent insights into crosstalk between genetic parasites and their host genome

The bulk of higher order organismal genomes is comprised of transposable element (TE) copies, i.e. genetic parasites. The host-parasite relation is multi-faceted, varying across genomic region (genic versus intergenic), life-cycle stages, tissue-type and of course in health versus pathological state. The reach of functional genomics though, in investigating genotype-to-phenotype relations, has been limited when TEs are involved. The aim of this review is to highlight recent progress made in understanding how TE origin biochemical activity interacts with the central dogma stages of the host genome. Such interaction can also bring about modulation of the immune context and this could have important repercussions in disease state where immunity has a role to play. Thus, the review is to instigate ideas and action points around identifying evolutionary adaptations that the host genome and the genetic parasite have evolved and why they could be relevant.

Role of BicDR in bristle shaft construction and support of BicD functions

Cell polarization requires asymmetric localization of numerous mRNAs, proteins and organelles. The movement of cargo towards the minus end of microtubules mostly depends on cytoplasmic dynein motors. In the dynein-dynactin-Bicaudal-D transport machinery, Bicaudal-D (BicD) links the cargo to the motor. Here, we focus on the role of Drosophila BicD-related (BicDR, CG32137) in the development of the long bristles. Together with BicD, it contributes to the organization and stability of the actin cytoskeleton in the not-yet-chitinized bristle shaft. BicD and BicDR also support the stable expression and distribution of Rab6 and Spn-F in the bristle shaft, including the distal tip localization of Spn-F, pointing to the role of microtubule-dependent vesicle trafficking for bristle construction. BicDR supports the function of BicD, and we discuss the hypothesis whereby BicDR might transport cargo more locally, with BicD transporting cargo over long distances, such as to the distal tip. We also identified embryonic proteins that interact with BicDR and appear to be BicDR cargo. For one of them, EF1γ (also known as eEF1γ), we show that the encoding gene EF1γ interacts with BicD and BicDR in the construction of the bristles.

BioKA: a curated and integrated biomarker knowledgebase for animals

Biomarkers play an important role in various area such as personalized medicine, drug development, clinical care, and molecule breeding. However, existing animals' biomarker resources predominantly focus on human diseases, leaving a significant gap in non-human animal disease understanding and breeding research. To address this limitation, we present BioKA (Biomarker Knowledgebase for Animals, https://ngdc.cncb.ac.cn/bioka), a curated and integrated knowledgebase encompassing multiple animal species, diseases/traits, and annotated resources. Currently, BioKA houses 16 296 biomarkers associated with 951 mapped diseases/traits across 31 species from 4747 references, including 11 925 gene/protein biomarkers, 1784 miRNA biomarkers, 1043 mutation biomarkers, 773 metabolic biomarkers, 357 circRNA biomarkers and 127 lncRNA biomarkers. Furthermore, BioKA integrates various annotations such as GOs, protein structures, protein-protein interaction networks, miRNA targets and so on, and constructs an interactive knowledge network of biomarkers including circRNA-miRNA-mRNA associations, lncRNA-miRNA associations and protein-protein associations, which is convenient for efficient data exploration. Moreover, BioKA provides detailed information on 308 breeds/strains of 13 species, and homologous annotations for 8784 biomarkers across 16 species, and offers three online application tools. The comprehensive knowledge provided by BioKA not only advances human disease research but also contributes to a deeper understanding of animal diseases and supports livestock breeding.

Kinesin-1 patterns Par-1 and Rho signaling at the cortex of syncytial embryos of Drosophila

The cell cortex of syncytial Drosophila embryos is patterned into cap and intercap regions by centrosomes, specific sets of proteins that are restricted to their respective regions by unknown mechanisms. Here, we found that Kinesin-1 is required for the restriction of plus- and minus-ends of centrosomal and non-centrosomal microtubules to the cap region, marked by EB1 and Patronin/Shot, respectively. Kinesin-1 also directly or indirectly restricts proteins and Rho signaling to the intercap, including the RhoGEF Pebble, Dia, Myosin II, Capping protein-α, and the polarity protein Par-1. Furthermore, we found that Par-1 is required for cap restriction of Patronin/Shot, and vice versa Patronin, for Par-1 enrichment at the intercap. In summary, our data support a model that Kinesin-1 would mediate the restriction of centrosomal and non-centrosomal microtubules to a region close to the centrosomes and exclude Rho signaling and Par-1. In addition, mutual antagonistic interactions would refine and maintain the boundary between cap and intercap and thus generate a distinct cortical pattern.

The genetic and molecular basis of haploinsufficiency in flowering plants

In diploid organisms, haploinsufficiency can be defined as the requirement for more than one fully functional copy of a gene. In contrast to most genes, whose loss-of-function alleles are recessive, loss-of-function alleles of haploinsufficient genes are dominant. However, forward and reverse genetic screens are biased toward obtaining recessive, loss-of-function mutations, and therefore, dominant mutations of all types are underrepresented in mutant collections. Despite this underrepresentation, haploinsufficient loci have intriguing implications for studies of genome evolution, gene dosage, stability of protein complexes, genetic redundancy, and gene expression. Here we review examples of haploinsufficiency in flowering plants and describe the underlying molecular mechanisms and evolutionary forces driving haploinsufficiency. Finally, we discuss the masking of haploinsufficiency by genetic redundancy, a widespread phenomenon among angiosperms.

Drosophila melanogaster Transcriptome Response to Different Wolbachia Strains

Wolbachia is a maternally inherited, intercellular bacterial symbiont of insects and some other invertebrates. Here, we investigated the effect of two different Wolbachia strains, differing in a large chromosomal inversion, on the differential expression of genes in D. melanogaster females. We revealed significant changes in the transcriptome of the infected flies compared to the uninfected ones, as well as in the transcriptome of flies infected with the Wolbachia strain, wMelPlus, compared to flies infected with the wMelCS112 strain. We linked differentially expressed genes (DEGs) from two pairwise comparisons, "uninfected-wMelPlus-infected" and "uninfected-wMelCS112-infected", into two gene networks, in which the following functional groups were designated: "Proteolysis", "Carbohydrate transport and metabolism", "Oxidation-reduction process", "Embryogenesis", "Transmembrane transport", "Response to stress" and "Alkaline phosphatases". Our data emphasized similarities and differences between infections by different strains under study: a wMelPlus infection results in more than double the number of upregulated DEGs and half the number of downregulated DEGs compared to a wMelCS112 infection. Thus, we demonstrated that Wolbachia made a significant contribution to differential expression of host genes and that the bacterial genotype plays a vital role in establishing the character of this contribution.

Finding information about uncharacterized Drosophila melanogaster genes

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

Evolution of the expression and regulation of the nuclear hormone receptor ERR gene family in the chordate lineage

The Estrogen Related Receptor (ERR) nuclear hormone receptor genes have a wide diversity of roles in vertebrate development. In embryos, ERR genes are expressed in several tissues, including the central and peripheral nervous systems. Here we seek to establish the evolutionary history of chordate ERR genes, their expression and their regulation. We examine ERR expression in mollusc, amphioxus and sea squirt embryos, finding the single ERR orthologue is expressed in the nervous system in all three, with muscle expression also found in the two chordates. We show that most jawed vertebrates and lampreys have four ERR paralogues, and that vertebrate ERR genes were ancestrally linked to Estrogen Receptor genes. One of the lamprey paralogues shares conserved expression domains with jawed vertebrate ERRγ in the embryonic vestibuloacoustic ganglion, eye, brain and spinal cord. Hypothesising that conserved expression derives from conserved regulation, we identify a suite of pan-vertebrate conserved non-coding sequences in ERR introns. We use transgenesis in lamprey and chicken embryos to show that these sequences are regulatory and drive reporter gene expression in the nervous system. Our data suggest an ancient association between ERR and the nervous system, including expression in cells associated with photosensation and mechanosensation. This includes the origin in the vertebrate common ancestor of a suite of regulatory elements in the 3' introns that drove nervous system expression and have been conserved from this point onwards.

Validation of a heat-inducible Ixodes scapularis HSP70 promoter and developing a tick-specific 3xP3 promoter sequence in ISE6 cells

Ixodes scapularis is an important vector of many pathogens, including the causative agent of Lyme disease, tick-borne encephalitis, and anaplasmosis. The study of gene function in I. scapularis and other ticks has been hampered by the lack of genetic tools, such as an inducible promoter to permit temporal control over transgenes encoding protein or double-stranded RNA expression. Studies of vector-pathogen relationships would also benefit from the capability to activate anti-pathogen genes at different times during pathogen infection and dissemination. We have characterized an intergenic sequence upstream of the heat shock protein 70 (HSP70) gene that can drive Renilla luciferase expression and mCherry fluorescence in the I. scapularis cell line ISE6. In another construct, we replaced the Drosophila melanogaster minimal HSP70 promoter in the synthetic 3xP3 promoter with a minimal portion of the I. scapularis HSP70 promoter and generated an I. scapularis specific 3xP3 (Is3xP3) promoter. Both promoter constructs, IsHSP70 and Is3xP3, allow for heat-inducible expression of mCherry fluorescence in ISE6 cells with an approximately 10-fold increase in the percentage of fluorescent positive cells upon exposure to a 2 h heat shock. These promoters described here will be valuable tools for gene function studies and temporal control of gene expression, including anti-pathogen genes.

Poly(ADP-ribosyl)ating enzymes coordinate changes in the expression of metabolic genes with developmental progression

Metabolism, known to be temporally regulated to meet evolving energy demands, plays a crucial role in shaping developmental pace. Recent studies have demonstrated that two key proteins PARP1 and PARG play a regulatory role in the transcription of both morphogenic and metabolic genes. Intriguingly, in Drosophila, the depletion of PARP1 or PARG proteins causes a developmental arrest before pupation, resulting in individuals unable to complete their development. This phenotype highlights the critical involvement of poly(ADP-ribosyl)ating enzymes in regulating the metamorphic process. In this study, we provide compelling evidence that these enzymes intricately coordinate transcriptional changes in both developmental and metabolic pathways during metamorphosis. Specifically, they promote the expression of genes crucial for pupation, while simultaneously negatively regulating the expression of metabolic genes before the transition to the pupal stage. Additionally, these enzymes suppress the expression of genes that are no longer required during this transformative period. Our findings shed light on the intricate interplay between poly(ADP-ribosyl)ating enzymes, developmental processes, and metabolic regulation before metamorphosis and highlight a new role of poly(ADP-ribosyl)ating enzymes in the global regulation of transcription.

Revised fission yeast gene and allele nomenclature guidelines for machine readability

Standardized nomenclature for genes, gene products, and isoforms is crucial to prevent ambiguity and enable clear communication of scientific data, facilitating efficient biocuration and data sharing. Standardized genotype nomenclature, which describes alleles present in a specific strain that differ from those in the wild-type reference strain, is equally essential to maximize research impact and ensure that results linking genotypes to phenotypes are Findable, Accessible, Interoperable, and Reusable (FAIR). In this publication, we extend the fission yeast clade gene nomenclature guidelines to support the curation efforts at PomBase (www.pombase.org), the Schizosaccharomyces pombe Model Organism Database. This update introduces nomenclature guidelines for noncoding RNA genes, following those set forth by the Human Genome Organisation Gene Nomenclature Committee. Additionally, we provide a significant update to the allele and genotype nomenclature guidelines originally published in 1987, to standardize the diverse range of genetic modifications enabled by the fission yeast genetic toolbox. These updated guidelines reflect a community consensus between numerous fission yeast researchers. Adoption of these rules will improve consistency in gene and genotype nomenclature, and facilitate machine-readability and automated entity recognition of fission yeast genes and alleles in publications or datasets. In conclusion, our updated guidelines provide a valuable resource for the fission yeast research community, promoting consistency, clarity, and FAIRness in genetic data sharing and interpretation.

Use of Farnesyl Transferase Inhibitors in an Ageing Model in Drosophila

The presence of farnesylated proteins at the inner nuclear membrane (INM), such as the Lamins or Kugelkern in Drosophila, leads to specific changes in the nuclear morphology and accelerated ageing on the organismal level reminiscent of the Hutchinson-Gilford progeria syndrome (HGPS). Farnesyl transferase inhibitors (FTIs) can suppress the phenotypes of the nuclear morphology in cultured fibroblasts from HGPS patients and cultured cells overexpressing farnesylated INM proteins. Similarly, FTIs have been reported to suppress the shortened lifespan in model organisms. Here, we report an experimental system combining cell culture and Drosophila flies for testing the activity of substances on the HGPS-like nuclear morphology and lifespan, with FTIs as an experimental example. Consistent with previous reports, we show that FTIs were able to ameliorate the nuclear phenotypes induced by the farnesylated nuclear proteins Progerin, Kugelkern, or truncated Lamin B in cultured cells. The subsequent validation in Drosophila lifespan assays demonstrated the applicability of the experimental system: treating adult Drosophila with the FTI ABT-100 reversed the nuclear phenotypes and extended the lifespan of experimentally induced short-lived flies. Since kugelkern-expressing flies have a significantly shorter average lifespan, half the time is needed for testing substances in the lifespan assay.

Sexually dimorphic mechanisms of VGLUT-mediated protection from dopaminergic neurodegeneration

Parkinson's disease (PD) targets some dopamine (DA) neurons more than others. Sex differences offer insights, with females more protected from DA neurodegeneration. The mammalian vesicular glutamate transporter VGLUT2 and Drosophila ortholog dVGLUT have been implicated as modulators of DA neuron resilience. However, the mechanisms by which VGLUT2/dVGLUT protects DA neurons remain unknown. We discovered DA neuron dVGLUT knockdown increased mitochondrial reactive oxygen species in a sexually dimorphic manner in response to depolarization or paraquat-induced stress, males being especially affected. DA neuron dVGLUT also reduced ATP biosynthetic burden during depolarization. RNA sequencing of VGLUT+ DA neurons in mice and flies identified candidate genes that we functionally screened to further dissect VGLUT-mediated DA neuron resilience across PD models. We discovered transcription factors modulating dVGLUT-dependent DA neuroprotection and identified dj-1β as a regulator of sex-specific DA neuron dVGLUT expression. Overall, VGLUT protects DA neurons from PD-associated degeneration by maintaining mitochondrial health.

Unanticipated domain requirements for Drosophila Wnk kinase in vivo

WNK (With no Lysine [K]) kinases have critical roles in the maintenance of ion homeostasis and the regulation of cell volume. Their overactivation leads to pseudohypoaldosteronism type II (Gordon syndrome) characterized by hyperkalemia and high blood pressure. More recently, WNK family members have been shown to be required for the development of the nervous system in mice, zebrafish, and flies, and the cardiovascular system of mice and fish. Furthermore, human WNK2 and Drosophila Wnk modulate canonical Wnt signaling. In addition to a well-conserved kinase domain, animal WNKs have a large, poorly conserved C-terminal domain whose function has been largely mysterious. In most but not all cases, WNKs bind and activate downstream kinases OSR1/SPAK, which in turn regulate the activity of various ion transporters and channels. Here, we show that Drosophila Wnk regulates Wnt signaling and cell size during the development of the wing in a manner dependent on Fray, the fly homolog of OSR1/SPAK. We show that the only canonical RF(X)V/I motif of Wnk, thought to be essential for WNK interactions with OSR1/SPAK, is required to interact with Fray in vitro. However, this motif is unexpectedly dispensable for Fray-dependent Wnk functions in vivo during fly development and fluid secretion in the Malpighian (renal) tubules. In contrast, a structure function analysis of Wnk revealed that the less-conserved C-terminus of Wnk, that recently has been shown to promote phase transitions in cell culture, is required for viability in vivo. Our data thus provide novel insights into unexpected in vivo roles of specific WNK domains.

Genotype-by-environment interactions influence the composition of the Drosophila seminal proteome

Ejaculate proteins are key mediators of post-mating sexual selection and sexual conflict, as they can influence both male fertilization success and female reproductive physiology. However, the extent and sources of genetic variation and condition dependence of the ejaculate proteome are largely unknown. Such knowledge could reveal the targets and mechanisms of post-mating selection and inform about the relative costs and allocation of different ejaculate components, each with its own potential fitness consequences. Here, we used liquid chromatography coupled with tandem mass spectrometry to characterize the whole-ejaculate protein composition across 12 isogenic lines of Drosophila melanogaster that were reared on a high- or low-quality diet. We discovered new proteins in the transferred ejaculate and inferred their origin in the male reproductive system. We further found that the ejaculate composition was mainly determined by genotype identity and genotype-specific responses to larval diet, with no clear overall diet effect. Nutrient restriction increased proteolytic protein activity and shifted the balance between reproductive function and RNA metabolism. Our results open new avenues for exploring the intricate role of genotypes and their environment in shaping ejaculate composition, or for studying the functional dynamics and evolutionary potential of the ejaculate in its multivariate complexity.

Identifying candidate de novo genes expressed in the somatic female reproductive tract of Drosophila melanogaster

Most eukaryotic genes have been vertically transmitted to the present from distant ancestors. However, variable gene number across species indicates that gene gain and loss also occurs. While new genes typically originate as products of duplications and rearrangements of preexisting genes, putative de novo genes-genes born out of ancestrally nongenic sequence-have been identified. Previous studies of de novo genes in Drosophila have provided evidence that expression in male reproductive tissues is common. However, no studies have focused on female reproductive tissues. Here we begin addressing this gap in the literature by analyzing the transcriptomes of 3 female reproductive tract organs (spermatheca, seminal receptacle, and parovaria) in 3 species-our focal species, Drosophila melanogaster-and 2 closely related species, Drosophila simulans and Drosophila yakuba, with the goal of identifying putative D. melanogaster-specific de novo genes expressed in these tissues. We discovered several candidate genes, located in sequence annotated as intergenic. Consistent with the literature, these genes tend to be short, single exon, and lowly expressed. We also find evidence that some of these genes are expressed in other D. melanogaster tissues and both sexes. The relatively small number of intergenic candidate genes discovered here is similar to that observed in the accessory gland, but substantially fewer than that observed in the testis.

Heuristic-enabled active machine learning: A case study of predicting essential developmental stage and immune response genes in Drosophila melanogaster

Computational prediction of absolute essential genes using machine learning has gained wide attention in recent years. However, essential genes are mostly conditional and not absolute. Experimental techniques provide a reliable approach of identifying conditionally essential genes; however, experimental methods are laborious, time and resource consuming, hence computational techniques have been used to complement the experimental methods. Computational techniques such as supervised machine learning, or flux balance analysis are grossly limited due to the unavailability of required data for training the model or simulating the conditions for gene essentiality. This study developed a heuristic-enabled active machine learning method based on a light gradient boosting model to predict essential immune response and embryonic developmental genes in Drosophila melanogaster. We proposed a new sampling selection technique and introduced a heuristic function which replaces the human component in traditional active learning models. The heuristic function dynamically selects the unlabelled samples to improve the performance of the classifier in the next iteration. Testing the proposed model with four benchmark datasets, the proposed model showed superior performance when compared to traditional active learning models (random sampling and uncertainty sampling). Applying the model to identify conditionally essential genes, four novel essential immune response genes and a list of 48 novel genes that are essential in embryonic developmental condition were identified. We performed functional enrichment analysis of the predicted genes to elucidate their biological processes and the result evidence our predictions. Immune response and embryonic development related processes were significantly enriched in the essential immune response and embryonic developmental genes, respectively. Finally, we propose the predicted essential genes for future experimental studies and use of the developed tool accessible at http://heal.covenantuniversity.edu.ng for conditional essentiality predictions.

Gene expression variation underlying tissue-specific responses to copper stress in Drosophila melanogaster

Copper is one of a handful of biologically necessary heavy metals that is also a common environmental pollutant. Under normal conditions, copper ions are required for many key physiological processes. However, in excess, copper quickly results in cell and tissue damage that can range in severity from temporary injury to permanent neurological damage. Because of its biological relevance, and because many conserved copper-responsive genes also respond to other non-essential heavy metal pollutants, copper resistance in Drosophila melanogaster is a useful model system with which to investigate the genetic control of the response to heavy metal stress. Because heavy metal toxicity has the potential to differently impact specific tissues, we genetically characterized the control of the gene expression response to copper stress in a tissue-specific manner in this study. We assessed the copper stress response in head and gut tissue of 96 inbred strains from the Drosophila Synthetic Population Resource (DSPR) using a combination of differential expression analysis and expression quantitative trait locus (eQTL) mapping. Differential expression analysis revealed clear patterns of tissue-specific expression, primarily driven by a more pronounced gene expression response in gut tissue. eQTL mapping of gene expression under control and copper conditions as well as for the change in gene expression following copper exposure (copper response eQTL) revealed hundreds of genes with tissue-specific local cis-eQTL and many distant trans-eQTL. eQTL associated with MtnA, Mdr49, Mdr50, and Sod3 exhibited genotype by environment effects on gene expression under copper stress, illuminating several tissue- and treatment-specific patterns of gene expression control. Together, our data build a nuanced description of the roles and interactions between allelic and expression variation in copper-responsive genes, provide valuable insight into the genomic architecture of susceptibility to metal toxicity, and highlight many candidate genes for future functional characterization.

The Drosophila Nab2 RNA binding protein inhibits m6A methylation and male-specific splicing of Sex lethal transcript in female neuronal tissue

The Drosophila polyadenosine RNA binding protein Nab2, which is orthologous to a human protein lost in a form of inherited intellectual disability, controls adult locomotion, axon projection, dendritic arborization, and memory through a largely undefined set of target RNAs. Here, we show a specific role for Nab2 in regulating splicing of ~150 exons/introns in the head transcriptome and focus on retention of a male-specific exon in the sex determination factor Sex-lethal (Sxl) that is enriched in female neurons. Previous studies have revealed that this splicing event is regulated in females by N6-methyladenosine (m6A) modification by the Mettl3 complex. At a molecular level, Nab2 associates with Sxl pre-mRNA in neurons and limits Sxl m6A methylation at specific sites. In parallel, reducing expression of the Mettl3, Mettl3 complex components, or the m6A reader Ythdc1 rescues mutant phenotypes in Nab2 flies. Overall, these data identify Nab2 as an inhibitor of m6A methylation and imply significant overlap between Nab2 and Mettl3 regulated RNAs in neuronal tissue.

The Drosophila gene encoding JIG protein (CG14850) is critical for CrebA nuclear trafficking during development

Coordination of mitochondrial and nuclear processes is key to the cellular health; however, very little is known about the molecular mechanisms regulating nuclear-mitochondrial crosstalk. Here, we report a novel molecular mechanism controlling the shuttling of CREB (cAMP response element-binding protein) protein complex between mitochondria and nucleoplasm. We show that a previously unknown protein, herein termed as Jig, functions as a tissue-specific and developmental timing-specific coregulator in the CREB pathway. Our results demonstrate that Jig shuttles between mitochondria and nucleoplasm, interacts with CrebA protein and controls its delivery to the nucleus, thus triggering CREB-dependent transcription in nuclear chromatin and mitochondria. Ablating the expression of Jig prevents CrebA from localizing to the nucleoplasm, affecting mitochondrial functioning and morphology and leads to Drosophila developmental arrest at the early third instar larval stage. Together, these results implicate Jig as an essential mediator of nuclear and mitochondrial processes. We also found that Jig belongs to a family of nine similar proteins, each of which has its own tissue- and time-specific expression profile. Thus, our results are the first to describe the molecular mechanism regulating nuclear and mitochondrial processes in a tissue- and time-specific manner.

Characterizing the Protein Isoforms of foraging (for), the PKGI Ortholog in Drosophila melanogaster

The foraging (for) gene of Drosophila melanogaster encodes a cGMP-dependent protein kinase (PKG), which is a major effector of the cGMP signaling pathway involved in the regulation of behaviour and metabolic traits. Despite being well studied at the transcript level, little is known about the for gene at the protein level. Here, we provide a detailed characterization of the for gene protein (FOR) products and present new tools for their study, including five isoform-specific antibodies and a transgenic strain that carries an HA-labelled for allele (for^BAC::HA). Our results showed that multiple FOR isoforms were expressed in the larval and adult stages of D. melanogaster and that the majority of whole-body FOR expression arises from three (P1, P1α, and P3) of eight putative protein isoforms. We found that FOR expression differed between the larval and adult stages and between the dissected larval organs we analyzed, which included the central nervous system (CNS), fat body, carcass, and intestine. Moreover, we showed that the FOR expression differed between two allelic variants of the for gene, namely, for^s (sitter) and for^R (rover), that are known to differ in many food-related traits. Together, our in vivo identification of FOR isoforms and the existence of temporal, spatial, and genetic differences in their expression lay the groundwork for determining their functional significance.

Role of BicDR in bristle shaft construction, tracheal development, and support of BicD functions

Cell polarization requires asymmetric localization of numerous mRNAs, proteins, and organelles. The movement of cargo towards the minus end of microtubules mostly depends on cytoplasmic dynein motors, which function as multiprotein complexes. In the dynein/dynactin/Bicaudal-D (DDB) transport machinery, Bicaudal-D (BicD) links the cargo to the motor. Here we focus on the role of BicD-related (BicDR) and its contribution to microtubule-dependent transport processes. Drosophila BicDR is required for the normal development of bristles and dorsal trunk tracheae. Together with BicD, it contributes to the organization and stability of the actin cytoskeleton in the not-yet-chitinized bristle shaft and the localization of Spn-F and Rab6 at the distal tip. We show that BicDR supports the function of BicD in bristle development and our results suggest that BicDR transports cargo more locally whereas BicD is more responsible for delivering functional cargo over the long distance to the distal tip. We identified the proteins that interact with BicDR and appear to be BicDR cargo in embryonic tissues. For one of them, EF1γ, we showed that EF1γ genetically interacts with BicD and BicDR in the construction of the bristles.

Comparative study of neuropeptide signaling systems in Hemiptera

Numerous physiological processes in insects are tightly regulated by neuropeptides and their receptors. Although they form an ancient signaling system, there is still a great deal of variety in neuropeptides and their receptors among different species within the same order. Neuropeptides and their receptors have been documented in many hemipteran insects, but the differences among them have been poorly characterized. Commercial grapevines worldwide are plagued by the bug Daktulosphaira vitifoliae (Hemiptera: Sternorrhyncha). Here, 33 neuropeptide precursors and 48 putative neuropeptide G protein-coupled receptor (GPCR) genes were identified in D. vitifoliae. Their expression profiles at the probe and feeding stages reflected potential regulatory roles in probe behavior. By comparison, we found that the Releasing Hormone-Related Peptides (GnRHs) system of Sternorrhyncha was differentiated from those of the other 2 suborders in Hemiptera. Independent secondary losses of the adipokinetic hormone/corazonin-related peptide receptor (ACP) and corazonin (CRZ) occurred during the evolution of Sternorrhyncha. Additionally, we discovered that the neuropeptide signaling systems of Sternorrhyncha were very different from those of Heteroptera and Auchenorrhyncha, which was consistent with Sternorrhyncha's phylogenetic position at the base of the order. This research provides more knowledge on neuropeptide systems and sets the groundwork for the creation of novel D. vitifoliae management strategies that specifically target these signaling pathways.

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

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