FlyBase: updates to the Drosophila melanogaster knowledge base

Investigation of epilepsy-related genes in a Drosophila model

Complex genetic architecture is the major cause of heterogeneity in epilepsy, which poses challenges for accurate diagnosis and precise treatment. A large number of epilepsy candidate genes have been identified from clinical studies, particularly with the widespread use of next-generation sequencing. Validating these candidate genes is emerging as a valuable yet challenging task. Drosophila serves as an ideal animal model for validating candidate genes associated with neurogenetic disorders such as epilepsy, due to its rapid reproduction rate, powerful genetic tools, and efficient use of ethological and electrophysiological assays. Here, we systematically summarize the advantageous techniques of the Drosophila model used to investigate epilepsy genes, including genetic tools for manipulating target gene expression, ethological assays for seizure-like behaviors, electrophysiological techniques, and functional imaging for recording neural activity. We then introduce several typical strategies for identifying epilepsy genes and provide new insights into gene‒gene interactions in epilepsy with polygenic causes. We summarize well-established precision medicine strategies for epilepsy and discuss prospective treatment options, including drug therapy and gene therapy for genetic epilepsy based on the Drosophila model. Finally, we also address genetic counseling and assisted reproductive technology as potential approaches for the prevention of genetic epilepsy.

Gene model for the ortholog of Ilp3 in Drosophila ananassae

Gene model for the ortholog of Insulin-like peptide 3 ( Ilp3 ) in the May 2011 (Agencourt dana_caf1/DanaCAF1) Genome Assembly (GenBank Accession: GCA_000005115.1 ) of Drosophila ananassae . This ortholog was characterized as part of a developing dataset to study the evolution of the Insulin/insulin-like growth factor signaling pathway (IIS) across the genus Drosophila using the Genomics Education Partnership gene annotation protocol for Course-based Undergraduate Research Experiences.

Gene model for the ortholog of lin-28 in Drosophila simulans

Gene model for the ortholog of lin-28 ( lin-28 ) in the May 2017 (Princeton ASM75419v2/DsimGB2) Genome Assembly (GenBank Accession: GCA_000754195.3 ) of Drosophila simulans . This ortholog was characterized as part of a developing dataset to study the evolution of the Insulin/insulin-like growth factor signaling pathway (IIS) across the genus Drosophila using the Genomics Education Partnership gene annotation protocol for Course-based Undergraduate Research Experiences.

Gene model for the ortholog of Roc1a in Drosophila ananassae

Gene model for the ortholog of Regulator of cullins 1a ( Roc1a ) in the Drosophila ananassae May 2011 (Agencourt dana_caf1/DanaCAF1) Genome Assembly (GenBank Accession: GCA_000005115.1 ). This ortholog was characterized as part of a developing dataset to study the evolution of the Insulin/insulin-like growth factor signaling pathway (IIS) across the genus Drosophila using the Genomics Education Partnership gene annotation protocol for Course-based Undergraduate Research Experiences.

Regulatory changes in the fatty acid elongase eloF underlie the evolution of sex-specific pheromone profiles in Drosophila prolongata

BACKGROUND

Pheromones play a key role in regulating sexual behavior throughout the animal kingdom. In Drosophila and other insects, many cuticular hydrocarbons (CHCs) are sexually dimorphic, and some are known to perform pheromonal functions. However, the genetic control of sex-specific CHC production is poorly understood outside of the model species D. melanogaster. A recent evolutionary change is found in D. prolongata, which, compared to its closest relatives, shows greatly increased sexual dimorphism in both CHCs and the chemosensory system responsible for their perception. A key transition involves a male-specific increase in the proportion of long-chain CHCs.

RESULTS

Perfuming D. prolongata females with the male-biased long-chain CHCs reduces copulation success, suggesting that these compounds function as sex pheromones. The evolutionary change in CHC profiles correlates with a male-specific increase in the expression of multiple genes involved in CHC biosynthesis, including fatty acid elongases, reductases and other key enzymes. In particular, elongase F, which is responsible for producing female-specific pheromones in D. melanogaster, is strongly upregulated in D. prolongata males compared both to females and to males of the sibling species. Mutations in eloF reduce the amount of long-chain CHCs, resulting in a partial feminization of pheromone profiles in D. prolongata males. Transgenic experiments show that sex-biased expression of eloF is caused in part by a putative transposable element honghaier insertion in its regulatory region.

CONCLUSIONS

These results show that cis-regulatory changes in the eloF gene, along with other changes in the CHC synthesis pathway, contribute to the evolution of sexual communication.

ArtSymbioCyc, a metabolic network database collection dedicated to arthropod symbioses: a case study, the tripartite cooperation in Sipha maydis

Most arthropods live in close association with bacteria. The genomes of associated partners have co-evolved, creating situations of interdependence that are complex to decipher despite the availability of their complete sequences. We developed ArtSymbioCyc, a metabolism-oriented database collection gathering genomic resources for arthropods and their associated bacteria. ArtSymbioCyc uses the powerful tools of the BioCyc community to produce high-quality annotations and to analyze and compare metabolic networks on a genome-wide scale. We used ArtSymbioCyc to study the case of the tripartite symbiosis of the cereal aphid Sipha maydis focusing on amino acid and vitamin metabolisms, as these compounds are known to be important in this strictly phloemophagous insect. We showed that the metabolic pathways of the insect host and its two obligate bacterial associates are interdependent and specialized in the exploitation of Poaceae phloem, particularly for the biosynthesis of sulfur-containing amino acids and most vitamins. This demonstrates that ArtSymbioCyc does not only reveal the individual metabolic capacities of each partner and their respective contributions to the holobiont they constitute but also allows to predict the essential inputs that must come from host nutrition.IMPORTANCEThe evolution has driven the emergence of complex arthropod-microbe symbiotic systems, whose metabolic integration is difficult to unravel. With its user-friendly interface, ArtSymbioCyc (https://artsymbiocyc.cycadsys.org) eases and speeds up the analysis of metabolic networks by enabling precise inference of compound exchanges between associated partners and helps unveil the adaptive potential of arthropods in contexts such as conservation or agricultural control.

Gene model for the ortholog of rictor in Drosophila ananassae

Gene model for the ortholog of rapamycin-insensitive companion of Tor ( rictor ) in the May 2011 (Agencourt dana_caf1/DanaCAF1) Genome Assembly (GenBank Accession: GCA_000005115.1 ) of Drosophila ananassae . This ortholog was characterized as part of a developing dataset to study the evolution of the Insulin/insulin-like growth factor signaling pathway (IIS) across the genus Drosophila using the Genomics Education Partnership gene annotation protocol for Course-based Undergraduate Research Experiences.

Complex determinants of R-loop formation at transposable elements and major DNA satellites

Aberrant activation of transposable elements (TEs) has been a well-documented source of genomic instability and disease, stemming from their insertion into genes and their imposition of epigenetic effects on nearby loci. However, the extent to which their disruptive effects involve concomitant or subsequent formation of DNA:RNA hybrids (R-loops) remains unknown. Here, we used DNA:RNA immunoprecipitation followed by high-throughput sequencing (DRIP-seq) to map the R-loop profiles of TEs and satellites in Drosophila melanogaster ovaries in control and rhino knockout flies, where dozens of TE families are derepressed. We observe that R-loops form primarily in LTR retrotransposons that carry A/T-rich sequence motifs, which are known to favor R-loop formation at genes in Drosophila and other species. We also report evidence of R-loop formation at 11 of 14 highly abundant D. melanogaster DNA satellites. R-loop formation is positively correlated with expression level for both TEs and satellites; however, neither sequence content nor expression fully explain which repeat families form R-loops, suggesting other factors are at play. Finally, by analyzing population frequencies of R-loop-forming TEs, we present evidence that TE copies with high R-loop signal may be under stronger negative selection, which suggests that R-loop formation by TEs may be deleterious to their host. Collectively, these results provide insight into the determinants of R-loop formation at repetitive elements.

Gene model for the ortholog of Glys in Drosophila busckii

Gene model for the ortholog of glycogen synthase ( Glys ) in the Drosophia busckii Sep. 2015 (UC Berkeley ASM127793v1/DbusGB1) Genome Assembly (GenBank Accession: GCA_001277935.1 ). This ortholog was characterized as part of a developing dataset to study the evolution of the Insulin/insulin-like growth factor signaling pathway (IIS) across the genus Drosophila using the Genomics Education Partnership gene annotation protocol for Course-based Undergraduate Research Experiences.

Gene model for the ortholog of Sik3 in Drosophila mojavensis

Gene model for the ortholog of Salt-inducible kinase 3 ( Sik3 ) in the May 2011 (Agencourt dmoj_caf1/DmojCAF1) Genome Assembly (GenBank Accession: GCA_000005175.1 ) of Drosophila mojavensis . This ortholog was characterized as part of a developing dataset to study the evolution of the Insulin/insulin-like growth factor signaling pathway (IIS) across the genus Drosophila using the Genomics Education Partnership gene annotation protocol for Course-based Undergraduate Research Experiences.

Gene model for the ortholog of wrd in Drosophila ananassae

Gene model for the ortholog of well-rounded ( wrd ) in the May 2011 (Agencourt dana_caf1/DanaCAF1) Genome Assembly (GenBank Accession: GCA_000005115.1 ) of Drosophila ananassae . This ortholog was characterized as part of a developing dataset to study the evolution of the Insulin/insulin-like growth factor signaling pathway (IIS) across the genus Drosophila using the Genomics Education Partnership gene annotation protocol for Course-based Undergraduate Research Experiences.

Gene model for the ortholog of Roc1a in Drosophila eugracilis

Gene model for the ortholog of Regulator of cullins 1a ( Roc1a ) in the Drosophila eugracilis Apr. 2013 (BCM-HGSC/Deug_2.0, DeugGB2) Genome Assembly (GenBank Accession: GCA_000236325.2). This ortholog was characterized as part of a developing dataset to study the evolution of the Insulin/insulin-like growth factor signaling pathway (IIS) across the genus Drosophila using the Genomics Education Partnership gene annotation protocol for Course-based Undergraduate Research Experiences.

COME: contrastive mapping learning for spatial reconstruction of single-cell RNA sequencing data

MOTIVATION

Single-cell RNA sequencing (scRNA-seq) enables high-throughput transcriptomic profiling at single-cell resolution. The inherent spatial location is crucial for understanding how single cells orchestrate multicellular functions and drive diseases. However, spatial information is often lost during tissue dissociation. Spatial transcriptomic (ST) technologies can provide precise spatial gene expression atlas, while their practicality is constrained by the number of genes they can assay or the associated costs at a larger scale and the fine-grained cell-type annotation. By transferring knowledge between scRNA-seq and ST data through cell correspondence learning, it is possible to recover the spatial properties inherent in scRNA-seq datasets.

RESULTS

In this study, we introduce COME, a COntrastive Mapping lEarning approach that learns mapping between ST and scRNA-seq data to recover the spatial information of scRNA-seq data. Extensive experiments demonstrate that the proposed COME method effectively captures precise cell-spot relationships and outperforms previous methods in recovering spatial location for scRNA-seq data. More importantly, our method is capable of precisely identifying biologically meaningful information within the data, such as the spatial structure of missing genes, spatial hierarchical patterns, and the cell-type compositions for each spot. These results indicate that the proposed COME method can help to understand the heterogeneity and activities among cells within tissue environments.

AVAILABILITY AND IMPLEMENTATION

The COME is freely available in GitHub (https://github.com/cindyway/COME).

Gene model for the ortholog of Pdk1 in Drosophila grimshawi

Gene model for the ortholog of Phosphoinositide-dependent kinase 1 ( Pdk1 ) in the May 2011 (Agencourt dgri_caf1/DgriCAF1) Genome Assembly (GenBank Accession: GCA_000005155.1 ) of Drosophila grimshawi . This ortholog was characterized as part of a developing dataset to study the evolution of the Insulin/insulin-like growth factor signaling pathway (IIS) across the genus Drosophila using the Genomics Education Partnership gene annotation protocol for Course-based Undergraduate Research Experiences.

The hypertrophic cardiomyopathy-associated A331P actin variant enhances basal contractile activity and elicits resting muscle dysfunction

Previous studies aimed at defining the mechanistic basis of hypertrophic cardiomyopathy caused by A331P cardiac actin have reported conflicting results. The mutation is located along an actin surface strand, proximal to residues that interact with tropomyosin. These F-actin-tropomyosin associations are vital for proper contractile inhibition. To help resolve disease pathogenesis, we implemented a multidisciplinary approach. Transgenic Drosophila, expressing A331P actin, displayed skeletal muscle hypercontraction and elevated basal myocardial activity. A331P thin filaments, reconstituted using recombinant human cardiac actin, exhibited higher in vitro myosin-based sliding speeds, exclusively at low Ca2+ concentrations. Cryo-EM-based reconstructions revealed no detectable A331P-related structural perturbations in F-actin. In silico, however, the P331-containing actin surface strand was less mobile and established diminished van der Waal's attractive forces with tropomyosin, which correlated with greater variability in inhibitory tropomyosin positioning. Such mutation-induced effects potentially elevate resting contractile activity among our models and may stimulate pathology in patients.

Gene model for the ortholog of Roc1a in Drosophila mojavensis

Gene model for the ortholog of Regulator of cullins 1a ( Roc1a ) in the Drosophila mojavensis May 2011 (Agencourt dmoj_caf1/DmojCAF1) Genome Assembly (GenBank Accession: GCA_000005175.1 ). This ortholog was characterized as part of a developing dataset to study the evolution of the Insulin/insulin-like growth factor signaling pathway (IIS) across the genus Drosophila using the Genomics Education Partnership gene annotation protocol for Course-based Undergraduate Research Experiences.

LitSumm: large language models for literature summarization of noncoding RNAs

Curation of literature in life sciences is a growing challenge. The continued increase in the rate of publication, coupled with the relatively fixed number of curators worldwide, presents a major challenge to developers of biomedical knowledgebases. Very few knowledgebases have resources to scale to the whole relevant literature and all have to prioritize their efforts. In this work, we take a first step to alleviating the lack of curator time in RNA science by generating summaries of literature for noncoding RNAs using large language models (LLMs). We demonstrate that high-quality, factually accurate summaries with accurate references can be automatically generated from the literature using a commercial LLM and a chain of prompts and checks. Manual assessment was carried out for a subset of summaries, with the majority being rated extremely high quality. We apply our tool to a selection of >4600 ncRNAs and make the generated summaries available via the RNAcentral resource. We conclude that automated literature summarization is feasible with the current generation of LLMs, provided that careful prompting and automated checking are applied. Database URL: https://rnacentral.org/.

Gene model for the ortholog of eIF4E1 in Drosophila yakuba

Gene model for the ortholog of eukaryotic translation initiation factor 4E1 ( eIF4E1 ) in the Dyak_CAF1 Genome Assembly (GenBank Accession: GCA_000005975.1 ) of Drosophila yakuba . This ortholog was characterized as part of a developing dataset to study the evolution of the Insulin/insulin-like growth factor signaling pathway (IIS) across the genus Drosophila using the Genomics Education Partnership gene annotation protocol for Course-based Undergraduate Research Experiences.

Gene model for the ortholog of eIF4E1 in Drosophila ananassae

Gene model for the ortholog of eukaryotic translation initiation factor 4E1 ( eIF4E1 ) in the May 2011 (Agencourt dana_caf1/DanaCAF1) Genome Assembly (GenBank Accession: GCA_000005115.1 ) of Drosophila ananassae . This ortholog was characterized as part of a developing dataset to study the evolution of the Insulin/insulin-like growth factor signaling pathway (IIS) across the genus Drosophila using the Genomics Education Partnership gene annotation protocol for Course-based Undergraduate Research Experiences.

The anillin knockdown in the Drosophila nervous system shows locomotor and learning defects

Anillin (Ani) is an evolutionarily conserved protein with a multi-domain structure that cross-links cytoskeletal proteins and plays an essential role in the formation of the contractile ring during cytokinesis. However, Ani is highly expressed in the human central nervous system (CNS), and it scaffolds myelin in the CNS of mice and modulates neuronal migration and growth in Caenorhabditis elegans. Although Ani is also highly expressed in the Drosophila CNS, its role remains unclear. In the present study, we showed that Ani is not only highly expressed in larval neuroblasts of the CNS, but also weakly expressed in the neuromuscular junction (NMJ) and axons. In addition, the ani knockdown in the nervous system led to pupal lethality, larval locomotor defects, and learning disability, along with abnormal morphology of the NMJ and distribution patterns of the mature neuropil in the CNS. These results show that Ani plays an important role also in the Drosophila nervous system.

Characterization and automated classification of sentences in the biomedical literature: a case study for biocuration of gene expression and protein kinase activity

Biological knowledgebases are essential resources for biomedical researchers, providing ready access to gene function and genomic data. Professional, manual curation of knowledgebases, however, is labor-intensive and thus high-performing machine learning methods that improve biocuration efficiency are needed. Here we report on sentence-level classification to identify biocuration-relevant sentences in the full text of published references for two gene function data types: gene expression and protein kinase activity. We performed a detailed characterization of sentences from references in the WormBase bibliography and used this characterization to define three tasks for classifying sentences as either 1) fully curatable, 2) fully and partially curatable, or 3) all language-related. We evaluated various machine learning (ML) models applied to these tasks and found that GPT and BioBERT achieve the highest average performance, resulting in F1 performance scores ranging from 0.89 to 0.99 depending upon the task. Our findings demonstrate the feasibility of extracting biocuration-relevant sentences from full text. Integrating these models into professional biocuration workflows, such as those used by the Alliance of Genome Resources and the ACKnowledge community curation platform, might well facilitate efficient and accurate annotation of the biomedical literature.

Gene model for the ortholog of Glys in Drosophila simulans

Gene model for the ortholog of glycogen synthase ( Glys ) in the Drosophila simulans May 2017 (Princeton ASM75419v2/DsimGB2) Genome Assembly (GenBank Accession: GCA_000754195.3 ). This ortholog was characterized as part of a developing dataset to study the evolution of the Insulin/insulin-like growth factor signaling pathway (IIS) across the genus Drosophila using the Genomics Education Partnership gene annotation protocol for Course-based Undergraduate Research Experiences.

FlyRNAi.org 2025 update-expanded resources for new technologies and species

The design, analysis and mining of large-scale 'omics studies with the goal of advancing biological and biomedical understanding require use of a range of bioinformatics tools, including approaches tailored to needs specific to a given species and/or technology. The FlyRNAi database at the Drosophila RNAi Screening Center and Transgenic RNAi Project (DRSC/TRiP) Functional Genomics Resources (https://fgr.hms.harvard.edu/tools) supports an increasingly broad group of technologies and species. Recently, for example, we expanded the database to include additional new data-centric resources that facilitate mining and analysis of single-cell transcriptomics. In addition, we have applied our approaches to CRISPR reagent and gene-centric bioinformatics approaches in Drosophila to arthropod vectors of infectious diseases. Building on our previous comprehensive reports on the FlyRNAi database, here we focus on new and updated resources with a primary focus on data-centric tools. Altogether, our suite of online resources supports various stages of functional genomics studies for Drosophila and other arthropods, and facilitate a wide range of reagent design, analysis, data mining and analysis approaches by biologists and biomedical experts studying Drosophila, other common genetic model species, arthropod vectors and/or human biology.

The filamins of Drosophila

The actin cytoskeleton is a dynamic mesh of filaments that provide structural support for cells and respond to external deformation forces. Active sensing of these forces is crucial for the function of the actin cytoskeleton, and some actin crosslinkers accomplish it. One such crosslinker is filamin, a highly conserved actin crosslinker dimeric protein with an elastic region capable of responding to mechanical changes in the actin cytoskeleton. Filamins are required across various cells and tissues. In Drosophila early and recent studies have provided many details about filamin functions. This review centers on the two Drosophila filamins encoded by the cheerio and jitterbu g genes. We examine the structural and evolutionary aspects of filamin genes in flies, contrasting them with those of other model organisms. Then, we synthesize phenotypic data across diverse cell types. Additionally, we outline the genetic tools available for both genes. We also propose to divide filamins into typical and atypical based on the number of actin-binding domains and their relationship with other filamins.

The three-dimensional genome drives the evolution of asymmetric gene duplicates via enhancer capture-divergence

Previous evolutionary models of duplicate gene evolution have overlooked the pivotal role of genome architecture. Here, we show that proximity-based regulatory recruitment by distally duplicated genes is an efficient mechanism for modulating tissue-specific production of preexisting proteins. By leveraging genomic asymmetries, we performed a coexpression analysis on Drosophila melanogaster tissue data to show the generality of enhancer capture-divergence (ECD) as a significant evolutionary driver of asymmetric, distally duplicated genes. We use the recently evolved gene HP6/Umbrea as an example of the ECD process. By assaying genome-wide chromosomal conformations in multiple Drosophila species, we show that HP6/Umbrea was inserted near a preexisting, long-distance three-dimensional genomic interaction. We then use this data to identify a newly found enhancer (FLEE1), buried within the coding region of the highly conserved, essential gene MFS18, that likely neofunctionalized HP6/Umbrea. Last, we demonstrate ancestral transcriptional coregulation of HP6/Umbrea's future insertion site, illustrating how enhancer capture provides a highly evolvable, one-step solution to Ohno's dilemma.

Gene model for the ortholog of Glys in Drosophila yakuba

Gene model for the ortholog of glycogen synthase ( Glys ) in the May 2011 (WUGSC dyak_caf1/DyakCAF1) Genome Assembly (GenBank Accession: GCA_000005975.1 ) of Drosophila yakuba . This ortholog was characterized as part of a developing dataset to study the evolution of the Insulin/insulin-like growth factor signaling pathway (IIS) across the genus Drosophila using the Genomics Education Partnership gene annotation protocol for Course-based Undergraduate Research Experiences.

Gene model for the ortholog of S6k in Drosophila yakuba

Gene model for the ortholog of Ribosomal protein S6 kinase ( S6k ) in the Dyak_CAF1 Genome Assembly (GenBank Accession: GCA_000005975.1) of Drosophila yakuba . This ortholog was characterized as part of a developing dataset to study the evolution of the Insulin/insulin-like growth factor signaling pathway (IIS) across the genus Drosophila using the Genomics Education Partnership gene annotation protocol for Course-based Undergraduate Research Experiences.

Gene model for the ortholog of Ilp5 in Drosophila ananassae

Gene model for the ortholog of Insulin-like peptide 5 ( Ilp5 ) in the D. ananassae May 2011 (Agencourt dana_caf1/DanaCAF1) Genome Assembly (GenBank Accession: GCA_000005115.1 ) of Drosophila ananassae . This ortholog was characterized as part of a developing dataset to study the evolution of the Insulin/insulin-like growth factor signaling pathway (IIS) across the genus Drosophila using the Genomics Education Partnership gene annotation protocol for Course-based Undergraduate Research Experiences.

Gene model for the ortholog of Dsor1 in Drosophila persimilis

Gene model for the ortholog of Downstream of raf1 ( Dsor1 ) in the May 2011 (Broad dper_caf1/DperCAF1) Genome Assembly (GenBank Accession: GCA_000005195.1 ) of Drosophila persimilis . This ortholog was characterized as part of a developing dataset to study the evolution of the Insulin/insulin-like growth factor signaling pathway (IIS) across the genus Drosophila using the Genomics Education Partnership gene annotation protocol for Course-based Undergraduate Research Experiences.

Gene model for the ortholog of ImpL2 in Drosophila simulans

Gene model for the ortholog of Ecdysone-inducible gene L2 ( ImpL2 ) in the May 2017 (Princeton ASM75419v2/DsimGB2) Genome Assembly (GenBank Accession: GCA_000754195.3 ) of Drosophila simulans . This ortholog was characterized as part of a developing dataset to study the evolution of the Insulin/insulin-like growth factor signaling pathway (IIS) across the genus Drosophila using the Genomics Education Partnership gene annotation protocol for Course-based Undergraduate Research Experiences.

Gene model for the ortholog Myc in Drosophila ananassae

Gene model for the ortholog of Myc ( Myc ) in the May 2011 (Agencourt dana_caf1/DanaCAF1) Genome Assembly (GenBank Accession: GCA_000005115.1 ) of Drosophila ananassae . This ortholog was characterized as part of a developing dataset to study the evolution of the Insulin/insulin-like growth factor signaling pathway (IIS) across the genus Drosophila using the Genomics Education Partnership gene annotation protocol for Course-based Undergraduate Research Experiences.

Steering veridical large language model analyses by correcting and enriching generated database queries: first steps toward ChatGPT bioinformatics

Large language models (LLMs) leverage factual knowledge from pretraining. Yet this knowledge remains incomplete and sometimes challenging to retrieve-especially in scientific domains not extensively covered in pretraining datasets and where information is still evolving. Here, we focus on genomics and bioinformatics. We confirm and expand upon issues with plain ChatGPT functioning as a bioinformatics assistant. Poor data retrieval and hallucination lead ChatGPT to err, as do incorrect sequence manipulations. To address this, we propose a system basing LLM outputs on up-to-date, authoritative facts and facilitating LLM-guided data analysis. Specifically, we introduce NagGPT, a middleware tool to insert between LLMs and databases, designed to bridge gaps in LLM knowledge and usage of database application programming interfaces. NagGPT proxies LLM-generated database queries, with special handling of incorrect queries. It acts as a gatekeeper between query responses and the LLM prompt, redirecting large responses to files but providing a synthesized snippet and injecting comments to steer the LLM. A companion OpenAI custom GPT, Genomics Fetcher-Analyzer, connects ChatGPT with NagGPT. It steers ChatGPT to generate and run Python code, performing bioinformatics tasks on data dynamically retrieved from a dozen common genomics databases (e.g. NCBI, Ensembl, UniProt, WormBase, and FlyBase). We implement partial mitigations for encountered challenges: detrimental interactions between code generation style and data analysis, confusion between database identifiers, and hallucination of both data and actions taken. Our results identify avenues to augment ChatGPT as a bioinformatics assistant and, more broadly, to improve factual accuracy and instruction following of unmodified LLMs.

Comparative Proteomics Highlights that GenX Exposure Leads to Metabolic Defects and Inflammation in Astrocytes

Exposure to PFAS such as GenX (HFPO dimer acid) has become increasingly common due to the replacement of older generation PFAS in manufacturing processes. While neurodegenerative and developmental effects of legacy PFAS exposure have been studied in depth, there is a limited understanding specific to the effects of GenX exposure. To investigate the effects of GenX exposure, we exposed Drosophila melanogaster to GenX and assessed the motor behavior and performed quantitative proteomics of fly brains to identify molecular changes in the brain. Additionally, metabolic network-based analysis using the iDrosophila1 model unveiled a potential link between GenX exposure and neurodegeneration. Since legacy PFAS exposure has been linked to Parkinson's disease (PD), we compared the proteome data sets between GenX-exposed flies and a fly model of PD expressing human α-synuclein. Considering the proteomic data- and network-based analyses that revealed GenX may be regulating GABA-associated pathways and the immune system, we next explored the effects of GenX on astrocytes, as astrocytes in the brain can regulate GABA. An array of assays demonstrated GenX exposure may lead to mitochondrial dysfunction and neuroinflammatory response in astrocytes, possibly linking non-cell autonomous neurodegeneration to the motor deficits associated with GenX exposure.

Gene model for the ortholog of DENR in Drosophila yakuba

Gene model for the ortholog of Density regulated protein ( DENR ) in the May 2011 (WUGSC dyak_caf1/DyakCAF1) Genome Assembly (GenBank Accession: GCA_000005975.1 ) of Drosophila yakuba . This ortholog was characterized as part of a developing dataset to study the evolution of the Insulin/insulin-like growth factor signaling pathway (IIS) across the genus Drosophila using the Genomics Education Partnership gene annotation protocol for Course-based Undergraduate Research Experiences.

Gene model for the ortholog of Thor in Drosophila yakuba

Gene model for the ortholog of Thor ( Thor ) in the D. yakuba May 2011 (WUGSC dyak_caf1/DyakCAF1) Genome Assembly (GenBank Accession: GCA_000005975.1 ) of Drosophila yakuba . This ortholog was characterized as part of a developing dataset to study the evolution of the Insulin/insulin-like growth factor signaling pathway (IIS) across the genus Drosophila using the Genomics Education Partnership gene annotation protocol for Course-based Undergraduate Research Experiences.

Cephalic ganglia transcriptomics of the American cockroach Periplaneta americana (Blattodea: Blattidae)

The American cockroach Periplaneta americana (L.) (Blattodea, Blattidae) has been a model organism for biochemical and physiological study for almost a century, however, its use does not benefit from the genetic tools found in key model species such as Drosophila melanogaster. To facilitate the use of the cockroach as a model system in neuroscience and to serve as a foundation for functional and translational experimentation, a transcriptome of the cephalic ganglia was assembled and annotated, and differential expression profiles between these ganglia were assessed. The transcriptome assembly yielded >400 k transcripts, with >40 k putative coding sequences. Gene ontology and protein domain searches indicate the cerebral and gnathal ganglia (GNG) have distinct genetic expression profiles. The developmental Toll signaling pathway appears to be active in the adult central nervous system (CNS), which may suggest a separate role for this pathway besides innate immune activation or embryonic development. The catabolic glycolytic and citric acid cycle enzymes are well represented in both ganglia, but key enzymes are more highly expressed in the GNG. Both ganglia express gluconeogenic and trehaloneogenic enzymes, suggesting a larger role of the CNS in regulating hemolymph sugar homeostasis than previously appreciated. The annotation and quantification of the cephalic ganglia transcriptome reveal both canonical and novel pathways in signaling and metabolism in an adult insect and lay a foundation for future functional and genetic analysis.

Regulatory Changes in the Fatty Acid Elongase eloF Underlie the Evolution of Sex-specific Pheromone Profiles in Drosophila prolongata

Pheromones play a key role in regulating sexual behavior throughout the animal kingdom. In Drosophila and other insects, many cuticular hydrocarbons (CHCs) are sexually dimorphic, and some are known to perform pheromonal functions. However, the genetic control of sex-specific CHC production is not understood outside of the model species D. melanogaster. A recent evolutionary change is found in D. prolongata, which, compared to its closest relatives, shows greatly increased sexual dimorphism in both CHCs and the chemosensory system responsible for their perception. A key transition involves a male-specific increase in the proportion of long-chain CHCs. Perfuming D. prolongata females with the male-biased CHCs reduces copulation success, suggesting that these compounds function as sex pheromones. The evolutionary change in CHC profiles correlates with a male-specific increase in the expression of multiple genes involved in CHC biosynthesis, including fatty acid elongases and reductases and other key enzymes. In particular, elongase F, which is responsible for producing female-specific pheromones in D. melanogaster, is strongly upregulated in D. prolongata males compared both to females and to males of the sibling species. Induced mutations in eloF reduce the amount of long-chain CHCs, resulting in a partial feminization of pheromone profiles in D. prolongata males while having minimal effect in females. Transgenic experiments show that sex-biased expression of eloF is caused in part by a putative transposable element insertion in its regulatory region. These results reveal one of the genetic mechanisms responsible for a recent evolutionary change in sexual communication.

The conserved genetic program of male germ cells uncovers ancient regulators of human spermatogenesis

Male germ cells share a common origin across animal species, therefore they likely retain a conserved genetic program that defines their cellular identity. However, the unique evolutionary dynamics of male germ cells coupled with their widespread leaky transcription pose significant obstacles to the identification of the core spermatogenic program. Through network analysis of the spermatocyte transcriptome of vertebrate and invertebrate species, we describe the conserved evolutionary origin of metazoan male germ cells at the molecular level. We estimate the average functional requirement of a metazoan male germ cell to correspond to the expression of approximately 10,000 protein-coding genes, a third of which defines a genetic scaffold of deeply conserved genes that has been retained throughout evolution. Such scaffold contains a set of 79 functional associations between 104 gene expression regulators that represent a core component of the conserved genetic program of metazoan spermatogenesis. By genetically interfering with the acquisition and maintenance of male germ cell identity, we uncover 161 previously unknown spermatogenesis genes and three new potential genetic causes of human infertility. These findings emphasize the importance of evolutionary history on human reproductive disease and establish a cross-species analytical pipeline that can be repurposed to other cell types and pathologies.

Gene model for the ortholog of Pten in Drosophila miranda

Gene model for the ortholog of Phosphatase and tensin homolog ( Pten ) in the D. miranda Apr. 2013 (UC Berkeley DroMir_2.2/DmirGB2) Genome Assembly (GenBank Accession: GCA_000269505.2 ) of Drosophila miranda . This ortholog was characterized as part of a developing dataset to study the evolution of the Insulin/insulin-like growth factor signaling pathway (IIS) across the genus Drosophila using the Genomics Education Partnership gene annotation protocol for Course-based Undergraduate Research Experiences.

Long non-coding RNAs involved in Drosophila development and regeneration

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

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.

Genome-Wide Association Study and transcriptome analysis reveals a complex gene network that regulates opsin gene expression and cell fate determination in Drosophila R7 photoreceptor cells

Background

An animal's ability to discriminate between differing wavelengths of light (i.e., color vision) is mediated, in part, by a subset of photoreceptor cells that express opsins with distinct absorption spectra. In Drosophila R7 photoreceptors, expression of the rhodopsin molecules, Rh3 or Rh4, is determined by a stochastic process mediated by the transcription factor spineless. The goal of this study was to identify additional factors that regulate R7 cell fate and opsin choice using a Genome Wide Association Study (GWAS) paired with transcriptome analysis via RNA-Seq.

Results

We examined Rh3 and Rh4 expression in a subset of fully-sequenced inbred strains from the Drosophila Genetic Reference Panel and performed a GWAS to identify 42 naturally-occurring polymorphisms-in proximity to 28 candidate genes-that significantly influence R7 opsin expression. Network analysis revealed multiple potential interactions between the associated candidate genes, spineless and its partners. GWAS candidates were further validated in a secondary RNAi screen which identified 12 lines that significantly reduce the proportion of Rh3 expressing R7 photoreceptors. Finally, using RNA-Seq, we demonstrated that all but four of the GWAS candidates are expressed in the pupal retina at a critical developmental time point and that five are among the 917 differentially expressed genes in sevenless mutants, which lack R7 cells.

Conclusions

Collectively, these results suggest that the relatively simple, binary cell fate decision underlying R7 opsin expression is modulated by a larger, more complex network of regulatory factors. Of particular interest are a subset of candidate genes with previously characterized neuronal functions including neurogenesis, neurodegeneration, photoreceptor development, axon growth and guidance, synaptogenesis, and synaptic function.

Gene model for the ortholog of Myc in Drosophila eugracilis

Gene model for the ortholog of Myc ( Myc ) in the D. eugracilis Apr. 2013 (BCM-HGSC/Deug_2.0) (DeugGB2) Genome Assembly (GenBank Accession: GCA_000236325.2) of Drosophila eugracilis . This ortholog was characterized as part of a developing dataset to study the evolution of the Insulin/insulin-like growth factor signaling pathway (IIS) across the genus Drosophila using the Genomics Education Partnership gene annotation protocol for Course-based Undergraduate Research Experiences.

Adenine nucleotide translocase 2 (Ant2) is required for individualization of spermatogenesis of Drosophila melanogaster

Successful completion of spermatogenesis is crucial for the perpetuation of the species. In Drosophila, spermatid individualization, a process involving changes in mitochondrial structure and function is critical to produce functional mature sperm. Ant2, encoding a mitochondrial adenine nucleotide translocase, is highly expressed in male testes and plays a role in energy metabolism in the mitochondria. However, its molecular function remains unclear. Here, we identified an important role of Ant2 in spermatid individualization. In Ant2 knockdown testes, spermatid individualization complexes composed of F-actin cones exhibited a diffuse distribution, and mature sperms were absent in the seminal vesicle, thus leading to male sterility. The most striking effects in Ant2-knockdown spermatids were decrease in tubulin polyglycylation and disruption of proper mitochondria derivatives function. Excessive apoptotic cells were also observed in Ant2-knockdown testes. To further investigate the phenotype of Ant2 knockdown in testes at the molecular level, complementary transcriptome and proteome analyses were performed. At the mRNA level, 868 differentially expressed genes were identified, of which 229 genes were upregulated and 639 were downregulated induced via Ant2 knockdown. iTRAQ-labeling proteome analysis revealed 350 differentially expressed proteins, of which 117 proteins were upregulated and 233 were downregulated. The expression of glutathione transferase (GstD5, GstE5, GstE8, and GstD3), proteins involved in reproduction were significantly regulated at both the mRNA and protein levels. These results indicate that Ant2 is crucial for spermatid maturation by affecting mitochondrial morphogenesis.

The genetic basis of incipient sexual isolation in Drosophila melanogaster

Speciation is a fundamental evolutionary process but the genetic changes accompanying speciation are difficult to determine since true species do not produce viable and fertile offspring. Partially reproductively isolated incipient species are useful for assessing genetic changes that occur prior to speciation. Drosophila melanogaster from Zimbabwe, Africa are partially sexually isolated from other D. melanogaster populations whose males have poor mating success with Zimbabwe females. We used the North American D. melanogaster Genetic Reference Panel (DGRP) to show that there is significant genetic variation in mating success of DGRP males with Zimbabwe females, to map genetic variants and genes associated with variation in mating success and to determine whether mating success to Zimbabwe females is associated with other quantitative traits previously measured in the DGRP. Incipient sexual isolation is highly polygenic and associated with the common African inversion In(3R)K and the amount of the sex pheromone 5,9-heptacosadiene in DGRP females. We functionally validated the effect of eight candidate genes using RNA interference to provide testable hypotheses for future studies investigating the molecular genetic basis of incipient sexual isolation in D. melanogaster.

Ecological genomics in the Northern krill uncovers loci for local adaptation across ocean basins

Krill are vital as food for many marine animals but also impacted by global warming. To learn how they and other zooplankton may adapt to a warmer world we studied local adaptation in the widespread Northern krill (Meganyctiphanes norvegica). We assemble and characterize its large genome and compare genome-scale variation among 74 specimens from the colder Atlantic Ocean and warmer Mediterranean Sea. The 19 Gb genome likely evolved through proliferation of retrotransposons, now targeted for inactivation by extensive DNA methylation, and contains many duplicated genes associated with molting and vision. Analysis of 760 million SNPs indicates extensive homogenizing gene-flow among populations. Nevertheless, we detect signatures of adaptive divergence across hundreds of genes, implicated in photoreception, circadian regulation, reproduction and thermal tolerance, indicating polygenic adaptation to light and temperature. The top gene candidate for ecological adaptation was nrf-6, a lipid transporter with a Mediterranean variant that may contribute to early spring reproduction. Such variation could become increasingly important for fitness in Atlantic stocks. Our study underscores the widespread but uneven distribution of adaptive variation, necessitating characterization of genetic variation among natural zooplankton populations to understand their adaptive potential, predict risks and support ocean conservation in the face of climate change.

Airway specific deregulation of asthma-related serpins impairs tracheal architecture and oxygenation in D. melanogaster

Serine proteases are important regulators of airway epithelial homeostasis. Altered serum or cellular levels of two serpins, Scca1 and Spink5, have been described for airway diseases but their function beyond antiproteolytic activity is insufficiently understood. To close this gap, we generated fly lines with overexpression or knockdown for each gene in the airways. Overexpression of both fly homologues of Scca1 and Spink5 induced the growth of additional airway branches, with more variable results for the respective knockdowns. Dysregulation of Scca1 resulted in a general delay in fruit fly development, with increases in larval and pupal mortality following overexpression of this gene. In addition, the morphological changes in the airways were concomitant with lower tolerance to hypoxia. In conclusion, the observed structural changes of the airways evidently had a strong impact on the airway function in our model as they manifested in a lower physical fitness of the animals. We assume that this is due to insufficient tissue oxygenation. Future work will be directed at the identification of key molecular regulators following the airway-specific dysregulation of Scca1 and Spink5 expression.

Infection and chronic disease activate a systemic brain-muscle signaling axis

Infections and neurodegenerative diseases induce neuroinflammation, but affected individuals often show nonneural symptoms including muscle pain and muscle fatigue. The molecular pathways by which neuroinflammation causes pathologies outside the central nervous system (CNS) are poorly understood. We developed multiple models to investigate the impact of CNS stressors on motor function and found that Escherichia coli infections and SARS-CoV-2 protein expression caused reactive oxygen species (ROS) to accumulate in the brain. ROS induced expression of the cytokine Unpaired 3 (Upd3) in Drosophila and its ortholog, IL-6, in mice. CNS-derived Upd3/IL-6 activated the JAK-STAT pathway in skeletal muscle, which caused muscle mitochondrial dysfunction and impaired motor function. We observed similar phenotypes after expressing toxic amyloid-β (Aβ42) in the CNS. Infection and chronic disease therefore activate a systemic brain-muscle signaling axis in which CNS-derived cytokines bypass the connectome and directly regulate muscle physiology, highlighting IL-6 as a therapeutic target to treat disease-associated muscle dysfunction.

NF-κB factors cooperate with Su(Hw)/E4F1 to balance Drosophila/human immune responses via modulating dynamic expression of miR-210

MicroRNAs (miRNAs) play crucial regulatory roles in controlling immune responses, but their dynamic expression mechanisms are poorly understood. Here, we firstly confirm that the conserved miRNA miR-210 negatively regulates innate immune responses of Drosophila and human via targeting Toll and TLR6, respectively. Secondly, our findings demonstrate that the expression of miR-210 is dynamically regulated by NF-κB factor Dorsal in immune response of Drosophila Toll pathway. Thirdly, we find that Dorsal-mediated transcriptional inhibition of miR-210 is dependent on the transcriptional repressor Su(Hw). Mechanistically, Dorsal interacts with Su(Hw) to modulate cooperatively the dynamic expression of miR-210 in a time- and dose-dependent manner, thereby controlling the strength of Drosophila Toll immune response and maintaining immune homeostasis. Fourthly, we reveal a similar mechanism in human cells, where NF-κB/RelA cooperates with E4F1 to regulate the dynamic expression of hsa-miR-210 in the TLR immune response. Overall, our study reveals a conservative regulatory mechanism that maintains animal innate immune homeostasis and provides new insights into the dynamic regulation of miRNA expression in immune response.

Rapid evolution of piRNA clusters in the Drosophila melanogaster ovary

The piRNA pathway is a highly conserved mechanism to repress transposable element (TE) activity in the animal germline via a specialized class of small RNAs called piwi-interacting RNAs (piRNAs). piRNAs are produced from discrete genomic regions called piRNA clusters (piCs). Although the molecular processes by which piCs function are relatively well understood in Drosophila melanogaster, much less is known about the origin and evolution of piCs in this or any other species. To investigate piC origin and evolution, we use a population genomic approach to compare piC activity and sequence composition across eight geographically distant strains of D. melanogaster with high-quality long-read genome assemblies. We perform annotations of ovary piCs and genome-wide TE content in each strain. Our analysis uncovers extensive variation in piC activity across strains and signatures of rapid birth and death of piCs. Most TEs inferred to be recently active show an enrichment of insertions into old and large piCs, consistent with the previously proposed "trap" model of piC evolution. In contrast, a small subset of active LTR families is enriched for the formation of new piCs, suggesting that these TEs have higher proclivity to form piCs. Thus, our findings uncover processes leading to the origin of piCs. We propose that piC evolution begins with the emergence of piRNAs from individual insertions of a few select TE families prone to seed new piCs that subsequently expand by accretion of insertions from most other TE families during evolution to form larger "trap" clusters. Our study shows that TEs themselves are the major force driving the rapid evolution of piCs.