FlyBase: updates to the Drosophila melanogaster knowledge base

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.

DMLS: an automated pipeline to extract the Drosophila modular transcription regulators and targets from massive literature articles

Transcription regulation in multicellular species is mediated by modular transcription factor (TF) binding site combinations termed cis-regulatory modules (CRMs). Such CRM-mediated transcription regulation determines the gene expression patterns during development. Biologists frequently investigate CRM transcription regulation on gene expressions. However, the knowledge of the target genes and regulatory TFs participating in the CRMs under study is mostly fragmentary throughout the literature. Researchers need to afford tremendous human resources to fully surf through the articles deposited in biomedical literature databases in order to obtain the information. Although several novel text-mining systems are now available for literature triaging, these tools do not specifically focus on CRM-related literature prescreening, failing to correctly extract the information of the CRM target genes and regulatory TFs from the literature. For this reason, we constructed a supportive auto-literature prescreener called Drosophila Modular transcription-regulation Literature Screener (DMLS) that achieves the following: (i) prescreens articles describing experiments on modular transcription regulation, (ii) identifies the described target genes and TFs of the CRMs under study for each modular transcription-regulation-describing article and (iii) features an automated and extendable pipeline to perform the task. We demonstrated that the final performance of DMLS in extracting the described target gene and regulatory TF lists of CRMs under study for given articles achieved test macro area under the ROC curve (auROC) = 89.7% and area under the precision-recall curve (auPRC) = 77.6%, outperforming the intuitive gene name-occurrence-counting method by at least 19.9% in auROC and 30.5% in auPRC. The web service and the command line versions of DMLS are available at https://cobis.bme.ncku.edu.tw/DMLS/  and  https://github.com/cobisLab/DMLS/, respectively. Database Tool URL: https://cobis.bme.ncku.edu.tw/DMLS/.

Osiris gene family defines the cuticle nanopatterns of Drosophila

Nanostructures of pores and protrusions in the insect cuticle modify molecular permeability and surface wetting and help insects sense various environmental cues. However, the cellular mechanisms that modify cuticle nanostructures are poorly understood. Here, we elucidate how insect-specific Osiris family genes are expressed in various cuticle-secreting cells in the Drosophila head during the early stages of cuticle secretion and cover nearly the entire surface of the head epidermis. Furthermore, we demonstrate how each sense organ cell with various cuticular nanostructures expressed a unique combination of Osiris genes. Osiris gene mutations cause various cuticle defects in the corneal nipples and pores of the chemosensory sensilla. Thus, our study emphasizes on the importance of Osiris genes for elucidating cuticle nanopatterning in insects.

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

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

The genetic basis of variation in Drosophila melanogaster mating behavior

Mating behavior is an essential fitness trait. We used the inbred, sequenced lines of the Drosophila Genetic Reference Panel (DGRP) to gain insights into the evolution of mating success and to evaluate the overlap in genetic architecture of mating behavior between the sexes. We found significant genetic variation for mating success when DGRP males and females from the same line were mated together, and when DGRP males and females were mated to an unrelated strain of the opposite sex. The mating success of DGRP males and females was not correlated when they were paired with the unrelated strain, suggesting independent genetic architecture of mating success in males and females that was confirmed by genome-wide association analyses. However, the mating success between pairs of the same or different DGRP lines was predicted accurately by the respective female and male mating success with the unrelated line.

The Arabidopsis Information Resource in 2024

Since 1999, The Arabidopsis Information Resource (www.arabidopsis.org) has been curating data about the Arabidopsis thaliana genome. Its primary focus is integrating experimental gene function information from the peer-reviewed literature and codifying it as controlled vocabulary annotations. Our goal is to produce a "gold standard" functional annotation set that reflects the current state of knowledge about the Arabidopsis genome. At the same time, the resource serves as a nexus for community-based collaborations aimed at improving data quality, access, and reuse. For the past decade, our work has been made possible by subscriptions from our global user base. This update covers our ongoing biocuration work, some of our modernization efforts that contribute to the first major infrastructure overhaul since 2011, the introduction of JBrowse2, and the resource's role in community activities such as organizing the structural reannotation of the genome. For gene function assessment, we used gene ontology annotations as a metric to evaluate: (1) what is currently known about Arabidopsis gene function and (2) the set of "unknown" genes. Currently, 74% of the proteome has been annotated to at least one gene ontology term. Of those loci, half have experimental support for at least one of the following aspects: molecular function, biological process, or cellular component. Our work sheds light on the genes for which we have not yet identified any published experimental data and have no functional annotation. Drawing attention to these unknown genes highlights knowledge gaps and potential sources of novel discoveries.

FlyBase: updates to the Drosophila genes and genomes database

FlyBase (flybase.org) is a model organism database and knowledge base about Drosophila melanogaster, commonly known as the fruit fly. Researchers from around the world rely on the genetic, genomic, and functional information available in FlyBase, as well as its tools to view and interrogate these data. In this article, we describe the latest developments and updates to FlyBase. These include the introduction of single-cell RNA sequencing data, improved content and display of functional information, updated orthology pipelines, new chemical reports, and enhancements to our outreach resources.

Updates to the Alliance of Genome Resources central infrastructure

The Alliance of Genome Resources (Alliance) is an extensible coalition of knowledgebases focused on the genetics and genomics of intensively studied model organisms. The Alliance is organized as individual knowledge centers with strong connections to their research communities and a centralized software infrastructure, discussed here. Model organisms currently represented in the Alliance are budding yeast, Caenorhabditis elegans, Drosophila, zebrafish, frog, laboratory mouse, laboratory rat, and the Gene Ontology Consortium. The project is in a rapid development phase to harmonize knowledge, store it, analyze it, and present it to the community through a web portal, direct downloads, and application programming interfaces (APIs). Here, we focus on developments over the last 2 years. Specifically, we added and enhanced tools for browsing the genome (JBrowse), downloading sequences, mining complex data (AllianceMine), visualizing pathways, full-text searching of the literature (Textpresso), and sequence similarity searching (SequenceServer). We enhanced existing interactive data tables and added an interactive table of paralogs to complement our representation of orthology. To support individual model organism communities, we implemented species-specific "landing pages" and will add disease-specific portals soon; in addition, we support a common community forum implemented in Discourse software. We describe our progress toward a central persistent database to support curation, the data modeling that underpins harmonization, and progress toward a state-of-the-art literature curation system with integrated artificial intelligence and machine learning (AI/ML).

The haplolethal gene wupA of Drosophila exhibits potential as a target for an X-poisoning gene drive

A synthetic gene drive that targets haplolethal genes on the X chromosome can skew the sex ratio toward males. Like an "X-shredder," it does not involve "homing," and that has advantages including the reduction of gene drive resistance allele formation. We examine this "X-poisoning" strategy by targeting 4 of the 11 known X-linked haplolethal/haplosterile genes of Drosophila melanogaster with CRISPR/Cas9. We find that targeting the wupA gene during spermatogenesis skews the sex ratio so fewer than 14% of progeny are daughters. That is unless we cross the mutagenic males to X^XY female flies that bear attached-X chromosomes, which reverses the inheritance of the poisoned X chromosome so that sons inherit it from their father, in which case only 2% of the progeny are sons. These sex ratio biases suggest that most of the CRISPR/Cas9 mutants we induced in the wupA gene are haplolethal but some are recessive lethal. The males generating wupA mutants do not suffer from reduced fertility; rather, the haplolethal mutants arrest development in the late stages of embryogenesis well after fertilized eggs have been laid. This provides a distinct advantage over genetic manipulation strategies involving sterility which can be countered by the remating of females. We also find that wupA mutants that destroy the nuclear localization signal of shorter isoforms are not haplolethal as long as the open reading frame remains intact. Like D. melanogaster, wupA orthologs of Drosophila suzukii and Anopheles mosquitos are found on X chromosomes making wupA a viable X-poisoning target in multiple species.

Enhancer-promoter interactions become more instructive in the transition from cell-fate specification to tissue differentiation

To regulate expression, enhancers must come in proximity to their target gene. However, the relationship between the timing of enhancer-promoter (E-P) proximity and activity remains unclear, with examples of uncoupled, anticorrelated and correlated interactions. To assess this, we selected 600 characterized enhancers or promoters with tissue-specific activity in Drosophila embryos and performed Capture-C in FACS-purified myogenic or neurogenic cells during specification and tissue differentiation. This enabled direct comparison between E-P proximity and activity transitioning from OFF-to-ON and ON-to-OFF states across developmental conditions. This showed remarkably similar E-P topologies between specified muscle and neuronal cells, which are uncoupled from activity. During tissue differentiation, many new distal interactions emerge where changes in E-P proximity reflect changes in activity. The mode of E-P regulation therefore appears to change as embryogenesis proceeds, from largely permissive topologies during cell-fate specification to more instructive regulation during terminal tissue differentiation, when E-P proximity is coupled to activation.

Oxidative stress-induced changes in the transcriptomic profile of extracellular vesicles

Extracellular vesicles (EVs) have been proposed to play dual roles in cellular homeostasis, functioning both to remove unwanted intracellular molecules, and to enable communication between cells as a means of modulating cellular responses in different physiological and pathological scenarios. EVs contain a broad range of cargoes, including multiple biotypes of RNA, which can vary depending on the cell status, and may function as signalling molecules. In this study, we carried out comparative transcriptomic analysis of Drosophila EVs and cells, demonstrating that the RNA profile of EVs is distinct from cells and shows dose-dependent changes in response to oxidative stress. We identified a high abundance of snoRNAs in EVs, alongside an enrichment of intronic and untranslated regions (UTRs) of mRNAs under stress. We also observed an increase in the relative abundance of either aberrant or modified mRNAs under stress. These findings suggest that EVs may function both for the elimination of specific cellular RNAs, and for the incorporation of RNAs that may hold signalling potential.

Chigno/CG11180 and SUMO are Chinmo-interacting proteins with a role in Drosophila testes somatic support cells

Stem cells are critical for replenishment of cells lost to death, damage or differentiation. Drosophila testes are a key model system for elucidating mechanisms regulating stem cell maintenance and differentiation. An intriguing gene identified through such studies is the transcription factor, chronologically inappropriate morphogenesis (Chinmo). Chinmo is a downstream effector of the Jak-STAT signaling pathway that acts in testis somatic stem cells to ensure maintenance of male stem cell fate and sexual identity. Defects in these processes can lead to infertility and the formation of germ cell tumors. While Chinmo's effect on testis stem cell behavior has been investigated in detail, there is still much to be learned about its structure, function, and interactions with other proteins. Using a two-hybrid screen, we find that Chinmo interacts with itself, the small ubiquitin-like modifier SUMO, the novel protein CG11180, and four other proteins (CG4318, Ova (ovaries absent), Taf3 (TBP-associated factor 3), and CG18269). Since both Chinmo and CG11180 contain sumoylation sites and SUMO-interacting motifs (SIMs), we analyzed their interaction in more detail. Using site-directed mutagenesis of a unique SIM in CG11180, we demonstrate that Chinmo's interaction with CG11180 is SUMO-dependent. Furthermore, to assess the functional relevance of both SUMO and CG11180, we performed RNAi-mediated knockdown of both proteins in somatic cells of the Drosophila testis. Using this approach, we find that CG11180 and SUMO are required in somatic cells of adult testes, and that reduction of either protein causes formation of germ cell tumors. Overall, our work suggests that SUMO may be involved in the interaction of Chinmo and CG11180 and that these genes are required in somatic cells of the adult Drosophila testis. Consistent with the CG11180 knockdown phenotype in male testes, and to underscore its connection to Chinmo, we propose the name Chigno (Childless Gambino) for CG11180.

Bradysia (Sciara) coprophila larvae up-regulate DNA repair pathways and down-regulate developmental regulators in response to ionizing radiation

The level of resistance to radiation and the developmental and molecular responses can vary between species, and even between developmental stages of one species. For flies (order: Diptera), prior studies concluded that the fungus gnat Bradysia (Sciara) coprophila (sub-order: Nematocera) is more resistant to irradiation-induced mutations that cause visible phenotypes than the fruit fly Drosophila melanogaster (sub-order: Brachycera). Therefore, we characterized the effects of and level of resistance to ionizing radiation on B. coprophila throughout its life cycle. Our data show that B. coprophila embryos are highly sensitive to even low doses of gamma-irradiation, whereas late-stage larvae can tolerate up to 80 Gy (compared to 40 Gy for D. melanogaster) and still retain their ability to develop to adulthood, though with a developmental delay. To survey the genes involved in the early transcriptional response to irradiation of B. coprophila larvae, we compared larval RNA-seq profiles with and without radiation treatment. The up-regulated genes were enriched for DNA damage response genes, including those involved in DNA repair, cell cycle arrest, and apoptosis, whereas the down-regulated genes were enriched for developmental regulators, consistent with the developmental delay of irradiated larvae. Interestingly, members of the PARP and AGO families were highly up-regulated in the B. coprophila radiation response. We compared the transcriptome responses in B. coprophila to the transcriptome responses in D. melanogaster from 3 previous studies: whereas pathway responses are highly conserved, specific gene responses are less so. Our study lays the groundwork for future work on the radiation responses in Diptera.

Deciphering the comprehensive knowledgebase landscape featuring infertility with IDDB Xtra

Infertility affects ∼15% of couples globally and half of cases are related to genetic disorders. Despite growing data and unprecedented improvements in high-throughput sequencing technologies, accumulated fertility-related issues concerning genetic diagnosis and potential treatment are urgent to be solved. However, there is a lack of comprehensive platforms that characterise various infertility-related records to provide research applications for exploring infertility in-depth and genetic counselling of infertility couple. To solve this problem, we provide IDDB Xtra by further integrating phenotypic manifestations, genomic datasets, epigenetics, modulators in collaboration with numerous interactive tools into our previous infertility database, IDDB. IDDB Xtra houses manually-curated 2369 genes of human and nine model organisms, 273 chromosomal abnormalities, 884 phenotypes, 60 genomic datasets, 464 epigenetic records, 1144 modulators relevant to infertility diagnosis and treatment. Additionally, IDDB Xtra incorporated customized graphical applications for researchers and clinicians to decipher in-depth disease mechanisms from the perspectives of developmental atlas, mutation effects, and clinical manifestations. Users can browse genes across developmental stages of human and mouse, filter candidate genes, mine potential variants and retrieve infertility biomedical network in an intuitive web interface. In summary, IDDB Xtra not only captures valuable research and data, but also provides useful applications to facilitate the genetic counselling and drug discovery of infertility. IDDB Xtra is freely available at https://mdl.shsmu.edu.cn/IDDB/and http://www.allostery.net/IDDB.

Analysis of SMAD1/5 target genes in a sea anemone reveals ZSWIM4-6 as a novel BMP signaling modulator

BMP signaling has a conserved function in patterning the dorsal-ventral body axis in Bilateria and the directive axis in anthozoan cnidarians. So far, cnidarian studies have focused on the role of different BMP signaling network components in regulating pSMAD1/5 gradient formation. Much less is known about the target genes downstream of BMP signaling. To address this, we generated a genome-wide list of direct pSMAD1/5 target genes in the anthozoan Nematostella vectensis, several of which were conserved in Drosophila and Xenopus. Our ChIP-seq analysis revealed that many of the regulatory molecules with documented bilaterally symmetric expression in Nematostella are directly controlled by BMP signaling. We identified several so far uncharacterized BMP-dependent transcription factors and signaling molecules, whose bilaterally symmetric expression may be indicative of their involvement in secondary axis patterning. One of these molecules is zswim4-6, which encodes a novel nuclear protein that can modulate the pSMAD1/5 gradient and potentially promote BMP-dependent gene repression.

Strategies for dissecting the complexity of neurodevelopmental disorders

Neurodevelopmental disorders (NDDs) are associated with a wide range of clinical features, affecting multiple pathways involved in brain development and function. Recent advances in high-throughput sequencing have unveiled numerous genetic variants associated with NDDs, which further contribute to disease complexity and make it challenging to infer disease causation and underlying mechanisms. Herein, we review current strategies for dissecting the complexity of NDDs using model organisms, induced pluripotent stem cells, single-cell sequencing technologies, and massively parallel reporter assays. We further highlight single-cell CRISPR-based screening techniques that allow genomic investigation of cellular transcriptomes with high efficiency, accuracy, and throughput. Overall, we provide an integrated review of experimental approaches that can be applicable for investigating a broad range of complex disorders.

Orai-mediated calcium entry determines activity of central dopaminergic neurons by regulation of gene expression

Maturation and fine-tuning of neural circuits frequently require neuromodulatory signals that set the excitability threshold, neuronal connectivity, and synaptic strength. Here, we present a mechanistic study of how neuromodulator-stimulated intracellular Ca2+ signals, through the store-operated Ca2+ channel Orai, regulate intrinsic neuronal properties by control of developmental gene expression in flight-promoting central dopaminergic neurons (fpDANs). The fpDANs receive cholinergic inputs for release of dopamine at a central brain tripartite synapse that sustains flight (Sharma and Hasan, 2020). Cholinergic inputs act on the muscarinic acetylcholine receptor to stimulate intracellular Ca2+ release through the endoplasmic reticulum (ER) localised inositol 1,4,5-trisphosphate receptor followed by ER-store depletion and Orai-mediated store-operated Ca2+ entry (SOCE). Analysis of gene expression in fpDANs followed by genetic, cellular, and molecular studies identified Orai-mediated Ca2+ entry as a key regulator of excitability in fpDANs during circuit maturation. SOCE activates the transcription factor trithorax-like (Trl), which in turn drives expression of a set of genes, including Set2, that encodes a histone 3 lysine 36 methyltransferase (H3K36me3). Set2 function establishes a positive feedback loop, essential for receiving neuromodulatory cholinergic inputs and sustaining SOCE. Chromatin-modifying activity of Set2 changes the epigenetic status of fpDANs and drives expression of key ion channel and signalling genes that determine fpDAN activity. Loss of activity reduces the axonal arborisation of fpDANs within the MB lobe and prevents dopamine release required for the maintenance of long flight.

Intrinsic protein disorder is insufficient to drive subnuclear clustering in embryonic transcription factors

Modern microscopy has revealed that core nuclear functions, including transcription, replication, and heterochromatin formation, occur in spatially restricted clusters. Previous work from our lab has shown that subnuclear high-concentration clusters of transcription factors may play a role in regulating RNA synthesis in the early Drosophila embryo. A nearly ubiquitous feature of eukaryotic transcription factors is that they contain intrinsically disordered regions (IDRs) that often arise from low complexity amino acid sequences within the protein. It has been proposed that IDRs within transcription factors drive co-localization of transcriptional machinery and target genes into high-concentration clusters within nuclei. Here, we test that hypothesis directly, by conducting a broad survey of the subnuclear localization of IDRs derived from transcription factors. Using a novel algorithm to identify IDRs in the Drosophila proteome, we generated a library of IDRs from transcription factors expressed in the early Drosophila embryo. We used this library to perform a high-throughput imaging screen in Drosophila Schneider-2 (S2) cells. We found that while subnuclear clustering does not occur when the majority of IDRs are expressed alone, it is frequently seen in full-length transcription factors. These results are consistent in live Drosophila embryos, suggesting that IDRs are insufficient to drive the subnuclear clustering behavior of transcription factors. Furthermore, the clustering of transcription factors in living embryos was unaffected by the deletion of IDR sequences. Our results demonstrate that IDRs are unlikely to be the primary molecular drivers of the clustering observed during transcription, suggesting a more complex and nuanced role for these disordered protein sequences.

DrosOMA: the Drosophila Orthologous Matrix browser

Background

Comparative genomic analyses to delineate gene evolutionary histories inform the understanding of organismal biology by characterising gene and gene family origins, trajectories, and dynamics, as well as enabling the tracing of speciation, duplication, and loss events, and facilitating the transfer of gene functional information across species. Genomic data are available for an increasing number of species from the genus Drosophila, however, a dedicated resource exploiting these data to provide the research community with browsable results from genus-wide orthology delineation has been lacking.

Methods

Using the OMA Orthologous Matrix orthology inference approach and browser deployment framework, we catalogued orthologues across a selected set of Drosophila species with high-quality annotated genomes. We developed and deployed a dedicated instance of the OMA browser to facilitate intuitive exploration, visualisation, and downloading of the genus-wide orthology delineation results.

Results

DrosOMA - the Drosophila Orthologous Matrix browser, accessible from https://drosoma.dcsr.unil.ch/ - presents the results of orthology delineation for 36 drosophilids from across the genus and four outgroup dipterans. It enables querying and browsing of the orthology data through a feature-rich web interface, with gene-view, orthologous group-view, and genome-view pages, including comprehensive gene name and identifier cross-references together with available functional annotations and protein domain architectures, as well as tools to visualise local and global synteny conservation.

Conclusions

The DrosOMA browser demonstrates the deployability of the OMA browser framework for building user-friendly orthology databases with dense sampling of a selected taxonomic group. It provides the Drosophila research community with a tailored resource of browsable results from genus-wide orthology delineation.

Integrating non-mammalian model organisms in the diagnosis of rare genetic diseases in humans

Next-generation sequencing technology has rapidly accelerated the discovery of genetic variants of interest in individuals with rare diseases. However, showing that these variants are causative of the disease in question is complex and may require functional studies. Use of non-mammalian model organisms - mainly fruitflies (Drosophila melanogaster), nematode worms (Caenorhabditis elegans) and zebrafish (Danio rerio) - enables the rapid and cost-effective assessment of the effects of gene variants, which can then be validated in mammalian model organisms such as mice and in human cells. By probing mechanisms of gene action and identifying interacting genes and proteins in vivo, recent studies in these non-mammalian model organisms have facilitated the diagnosis of numerous genetic diseases and have enabled the screening and identification of therapeutic options for patients. Studies in non-mammalian model organisms have also shown that the biological processes underlying rare diseases can provide insight into more common mechanisms of disease and the biological functions of genes. Here, we discuss the opportunities afforded by non-mammalian model organisms, focusing on flies, worms and fish, and provide examples of their use in the diagnosis of rare genetic diseases.

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.

Comparison of insect and human cytochrome b561 proteins: Insights into candidate ferric reductases in insects

Cytochrome b561 (cytb561) proteins comprise a family of transmembrane oxidoreductases that transfer single electrons across a membrane. Most eukaryotic species, including insects, possess multiple cytb561 homologs. To learn more about this protein family in insects, we carried out a bioinformatics-based investigation of cytb561 family members from nine species representing eight insect orders. We performed a phylogenetic analysis to classify insect cytb561 ortholog groups. We then conducted sequence analyses and analyzed protein models to predict structural elements that may impact the biological functions and localization of these proteins, with a focus on possible ferric reductase activity. Our study revealed three ortholog groups, designated CG1275, Nemy, and CG8399, and a fourth group of less-conserved genes. We found that CG1275 and Nemy proteins are similar to a human ferric reductase, duodenal cytochrome b561 (Dcytb), and have many conserved amino acid residues that function in substrate binding in Dcytb. Notably, CG1275 and Nemy proteins contain a conserved histidine and other residues that play a role in ferric ion reduction by Dcytb. Nemy proteins were distinguished by a novel cysteine-rich cytoplasmic loop sequence. CG8399 orthologs are similar to a putative ferric reductase in humans, stromal cell-derived receptor 2. Like other members of the CYBDOM class of cytb561 proteins, these proteins contain reeler, DOMON, and cytb561 domains. Drosophila melanogaster CG8399 is the only insect cytb561 with known ferric reductase activity. Our investigation of the DOMON domain in CG8399 proteins revealed a probable heme-binding site and a possible site for ferric reduction. The fourth group includes a subgroup of proteins with a conserved "KXXXXKXH" non-cytoplasmic loop motif that may be a substrate binding site and is present in a potential ferric reductase, human tumor suppressor cytochrome b561. This study provides a foundation for future investigations of the biological functions of cytb561 genes in insects.

Systematic assessment of transcriptomic and metabolic reprogramming by blue light exposure coupled with aging

The prevalent use of light-emitting diodes (LEDs) has caused revolutionary changes in modern life, but the potential hazards to health of blue light are poorly understood. N6-methyladenosine (m6A) is the most prevalent posttranscriptional modification in eukaryotes and can modulate diverse physiological processes by regulating mRNA fate. Here, to understand the effects and molecular mechanisms of daily low-intensity blue light exposure (BLE) and ascertain whether m6A methylation plays a role in BLE-induced phenotypes, we constructed a series of Drosophila models under different durations of daily low-intensity BLE and obtained multiomics profiles. Our results revealed that BLE could induce transcriptomic, m6A epitranscriptomic, and metabolomic reprogramming in Drosophila along with aging process. Importantly, the m6A methylation sites enriched in the 5' untranslated regions (UTRs) of Drosophila transcripts showed strong age specificity and could be altered by BLE. We experimentally validated that aging-related gene Tor and circadian rhythm-related gene per were regulated by 5' UTR-enriched m6A methylation. Overall, our study provides a systematic assessment of m6A RNA methylome reprogramming by BLE and aging in Drosophila model.

Artificial intelligence for neurodegenerative experimental models

INTRODUCTION

Experimental models are essential tools in neurodegenerative disease research. However, the translation of insights and drugs discovered in model systems has proven immensely challenging, marred by high failure rates in human clinical trials.

METHODS

Here we review the application of artificial intelligence (AI) and machine learning (ML) in experimental medicine for dementia research.

RESULTS

Considering the specific challenges of reproducibility and translation between other species or model systems and human biology in preclinical dementia research, we highlight best practices and resources that can be leveraged to quantify and evaluate translatability. We then evaluate how AI and ML approaches could be applied to enhance both cross-model reproducibility and translation to human biology, while sustaining biological interpretability.

DISCUSSION

AI and ML approaches in experimental medicine remain in their infancy. However, they have great potential to strengthen preclinical research and translation if based upon adequate, robust, and reproducible experimental data.

HIGHLIGHTS

There are increasing applications of AI in experimental medicine. We identified issues in reproducibility, cross-species translation, and data curation in the field. Our review highlights data resources and AI approaches as solutions. Multi-omics analysis with AI offers exciting future possibilities in drug discovery.

Updates to the Alliance of Genome Resources Central Infrastructure Alliance of Genome Resources Consortium

The Alliance of Genome Resources (Alliance) is an extensible coalition of knowledgebases focused on the genetics and genomics of intensively-studied model organisms. The Alliance is organized as individual knowledge centers with strong connections to their research communities and a centralized software infrastructure, discussed here. Model organisms currently represented in the Alliance are budding yeast, C. elegans, Drosophila, zebrafish, frog, laboratory mouse, laboratory rat, and the Gene Ontology Consortium. The project is in a rapid development phase to harmonize knowledge, store it, analyze it, and present it to the community through a web portal, direct downloads, and APIs. Here we focus on developments over the last two years. Specifically, we added and enhanced tools for browsing the genome (JBrowse), downloading sequences, mining complex data (AllianceMine), visualizing pathways, full-text searching of the literature (Textpresso), and sequence similarity searching (SequenceServer). We enhanced existing interactive data tables and added an interactive table of paralogs to complement our representation of orthology. To support individual model organism communities, we implemented species-specific "landing pages" and will add disease-specific portals soon; in addition, we support a common community forum implemented in Discourse. We describe our progress towards a central persistent database to support curation, the data modeling that underpins harmonization, and progress towards a state-of-the art literature curation system with integrated Artificial Intelligence and Machine Learning (AI/ML).

Comparative Population Transcriptomics Provide New Insight into the Evolutionary History and Adaptive Potential of World Ocean Krill

Genetic variation is instrumental for adaptation to changing environments but it is unclear how it is structured and contributes to adaptation in pelagic species lacking clear barriers to gene flow. Here, we applied comparative genomics to extensive transcriptome datasets from 20 krill species collected across the Atlantic, Indian, Pacific, and Southern Oceans. We compared genetic variation both within and between species to elucidate their evolutionary history and genomic bases of adaptation. We resolved phylogenetic interrelationships and uncovered genomic evidence to elevate the cryptic Euphausia similis var. armata into species. Levels of genetic variation and rates of adaptive protein evolution vary widely. Species endemic to the cold Southern Ocean, such as the Antarctic krill Euphausia superba, showed less genetic variation and lower evolutionary rates than other species. This could suggest a low adaptive potential to rapid climate change. We uncovered hundreds of candidate genes with signatures of adaptive evolution among Antarctic Euphausia but did not observe strong evidence of adaptive convergence with the predominantly Arctic Thysanoessa. We instead identified candidates for cold-adaptation that have also been detected in Antarctic fish, including genes that govern thermal reception such as TrpA1. Our results suggest parallel genetic responses to similar selection pressures across Antarctic taxa and provide new insights into the adaptive potential of important zooplankton already affected by climate change.

Inferring multi-locus selection in admixed populations

Admixture, the exchange of genetic information between distinct source populations, is thought to be a major source of adaptive genetic variation. Unlike mutation events, which periodically generate single alleles, admixture can introduce many selected alleles simultaneously. As such, the effects of linkage between selected alleles may be especially pronounced in admixed populations. However, existing tools for identifying selected mutations within admixed populations only account for selection at a single site, overlooking phenomena such as linkage among proximal selected alleles. Here, we develop and extensively validate a method for identifying and quantifying the individual effects of multiple linked selected sites on a chromosome in admixed populations. Our approach numerically calculates the expected local ancestry landscape in an admixed population for a given multi-locus selection model, and then maximizes the likelihood of the model. After applying this method to admixed populations of Drosophila melanogaster and Passer italiae, we found that the impacts between linked sites may be an important contributor to natural selection in admixed populations. Furthermore, for the situations we considered, the selection coefficients and number of selected sites are overestimated in analyses that do not consider the effects of linkage among selected sites. Our results imply that linkage among selected sites may be an important evolutionary force in admixed populations. This tool provides a powerful generalized method to investigate these crucial phenomena in diverse populations.

Smurfness-based two-phase model of ageing helps deconvolve the ageing transcriptional signature

Ageing is characterised at the molecular level by six transcriptional 'hallmarks of ageing', that are commonly described as progressively affected as time passes. By contrast, the 'Smurf' assay separates high-and-constant-mortality risk individuals from healthy, zero-mortality risk individuals, based on increased intestinal permeability. Performing whole body total RNA sequencing, we found that Smurfness distinguishes transcriptional changes associated with chronological age from those associated with biological age. We show that transcriptional heterogeneity increases with chronological age in non-Smurf individuals preceding the other five hallmarks of ageing that are specifically associated with the Smurf state. Using this approach, we also devise targeted pro-longevity genetic interventions delaying entry in the Smurf state. We anticipate that increased attention to the evolutionary conserved Smurf phenotype will bring about significant advances in our understanding of the mechanisms of ageing.

RNAinsecta: A tool for prediction of precursor microRNA in insects and search for their target in the model organism Drosophila melanogaster

INTRODUCTION AND BACKGROUND

Pre-MicroRNAs are the hairpin loops from which microRNAs are produced that have been found to negatively regulate gene expression in several organisms. In insects, microRNAs participate in several biological processes including metamorphosis, reproduction, immune response, etc. Numerous tools have been designed in recent years to predict novel pre-microRNA using binary machine learning classifiers where prediction models are trained with true and pseudo pre-microRNA hairpin loops. Currently, there are no existing tool that is exclusively designed for insect pre-microRNA detection.

AIM

Application of machine learning algorithms to develop an open source tool for prediction of novel precursor microRNA in insects and search for their miRNA targets in the model insect organism, Drosophila melanogaster.

METHODS

Machine learning algorithms such as Random Forest, Support Vector Machine, Logistic Regression and K-Nearest Neighbours were used to train insect true and false pre-microRNA features with 10-fold Cross Validation on SMOTE and Near-Miss datasets. miRNA targets IDs were collected from miRTarbase and their corresponding transcripts were collected from FlyBase. We used miRanda algorithm for the target searching.

RESULTS

In our experiment, SMOTE performed significantly better than Near-Miss for which it was used for modelling. We kept the best performing parameters after obtaining initial mean accuracy scores >90% of Cross Validation. The trained models on Support Vector Machine achieved accuracy of 92.19% while the Random Forest attained an accuracy of 80.28% on our validation dataset. These models are hosted online as web application called RNAinsecta. Further, searching target for the predicted pre-microRNA in Drosophila melanogaster has been provided in RNAinsecta.

Genetic mapping of the p47 L.3.2 mutation in Drosophilamelanogaster

An EMS-based forward genetic screen was conducted in an apoptotic null background to identify genetic aberrations that contribute to regulation of cell growth in Drosophila melanogaster . The current work maps the genomic location of one of the identified mutants, L.3.2 . Genetic crosses conducted through the Fly-CURE consortium determined that the gene locus for the L.3.2 mutation is p47 on chromosome 2R.

Enthralling genetic regulatory mechanisms meddling insecticide resistance development in insects: role of transcriptional and post-transcriptional events

Insecticide resistance in insects severely threatens both human health and agriculture, making insecticides less compelling and valuable, leading to frequent pest management failures, rising input costs, lowering crop yields, and disastrous public health. Insecticide resistance results from multiple factors, mainly indiscriminate insecticide usage and mounted selection pressure on insect populations. Insects respond to insecticide stress at the cellular level by modest yet significant genetic propagations. Transcriptional, co-transcriptional, and post-transcriptional regulatory signals of cells in organisms regulate the intricate processes in gene expressions churning the genetic information in transcriptional units into proteins and non-coding transcripts. Upregulation of detoxification enzymes, notably cytochrome P450s (CYPs), glutathione S-transferases (GSTs), esterases [carboxyl choline esterase (CCE), carboxyl esterase (CarE)] and ATP Binding Cassettes (ABC) at the transcriptional level, modification of target sites, decreased penetration, or higher excretion of insecticides are the noted insect physiological responses. The transcriptional regulatory pathways such as AhR/ARNT, Nuclear receptors, CncC/Keap1, MAPK/CREB, and GPCR/cAMP/PKA were found to regulate the detoxification genes at the transcriptional level. Post-transcriptional changes of non-coding RNAs (ncRNAs) such as microRNAs (miRNA), long non-coding RNAs (lncRNA), and epitranscriptomics, including RNA methylation, are reported in resistant insects. Additionally, genetic modifications such as mutations in the target sites and copy number variations (CNV) are also influencing insecticide resistance. Therefore, these cellular intricacies may decrease insecticide sensitivity, altering the concentrations or activities of proteins involved in insecticide interactions or detoxification. The cellular episodes at the transcriptional and post-transcriptional levels pertinent to insecticide resistance responses in insects are extensively covered in this review. An overview of molecular mechanisms underlying these biological rhythms allows for developing alternative pest control methods to focus on insect vulnerabilities, employing reverse genetics approaches like RNA interference (RNAi) technology to silence particular resistance-related genes for sustained insect management.

Transcriptional pausing factor M1BP regulates cellular homeostasis by suppressing autophagy and apoptosis in Drosophila eye

During organogenesis cellular homeostasis plays a crucial role in patterning and growth. The role of promoter proximal pausing of RNA polymerase II, which regulates transcription of several developmental genes by GAGA factor or Motif 1 Binding Protein (M1BP), has not been fully understood in cellular homeostasis. Earlier, we reported that M1BP, a functional homolog of ZKSCAN3, regulates wingless and caspase-dependent cell death (apoptosis) in the Drosophila eye. Further, blocking apoptosis does not fully rescue the M1BPRNAi phenotype of reduced eye. Therefore, we looked for other possible mechanism(s). In a forward genetic screen, members of the Jun-amino-terminal-(NH2)-Kinase (JNK) pathway were identified. Downregulation of M1BP ectopically induces JNK, a pro-death pathway known to activate both apoptosis and caspase-independent (autophagy) cell death. Activation of JNK pathway components can enhance M1BPRNAi phenotype and vice-versa. Downregulation of M1BP ectopically induced JNK signaling, which leads to apoptosis and autophagy. Apoptosis and autophagy are regulated independently by their genetic circuitry. Here, we found that blocking either apoptosis or autophagy alone rescues the reduced eye phenotype of M1BP downregulation; whereas, blocking both apoptosis and autophagy together significantly rescues the M1BP reduced eye phenotype to near wild-type in nearly 85% progeny. This data suggests that the cellular homeostasis response demonstrated by two independent cell death mechanisms, apoptosis and autophagy, can be regulated by a common transcriptional pausing mechanism orchestrated by M1BP. Since these fundamental processes are conserved in higher organisms, this novel functional link between M1BP and regulation of both apoptosis and autophagy can be extrapolated to humans.

Chromosomal-level Genome Assembly of the Coffee Bee Hawk Moth Reveals the Evolution of Chromosomes and the Molecular Basis of Distinct Phenotypes

Cephonodes hylas, the coffee bee hawk moth is a hawk moth species with unique characteristics, such as larvae feeding on gardenia, overcoming the toxicity of its iridoid glycosides, diurnal adults, and transparent wings. Although C. hylas is a fascinating model for molecular biological research, genome sequence analysis-based genetic approaches to elucidate these peculiarities have not yet been undertaken. We successfully achieved de novo genome assembly at the chromosome level of C. hylas comparable to the Lepidoptera model organism, silkworm. Additionally, 16,854 protein-coding genes were annotated, and the constructed genome sequence and annotated genes were of the highest quality BUSCO completion compared to closely related species. Comparative genome analysis revealed the process of chromosomal evolution from the Bombycoidea ancestral (n = 31) genome and changes in turnover at the chromosome level associated with chromosomal fusion events, such as the rate of repetitive sequence insertion. These analyses were only possible because the genome was constructed at the chromosome level. Additionally, increased the nonsynonymous/synonymous rate (dN/dS) ratios were observed in multiple photoreceptor-related genes that were strongly associated with the acquisition of diurnal activity. Furthermore, tandemly duplicated expanded genes containing many digestive and other enzymes and larval midgut-specific expression were also confirmed. These genes may be involved in the metabolism of genipin, a toxin found in gardenias. Using the genome sequence of C. hylas determined at the chromosome level, we have successfully identified new insights into the chromosomal evolution of Bombycoidea, as well as the relationship between the genome sequence and its characteristic traits.

The Genomic Basis of Adaptation to High Elevations in Africanized Honey Bees

A range of different genetic architectures underpin local adaptation in nature. Honey bees (Apis mellifera) in the Eastern African Mountains harbor high frequencies of two chromosomal inversions that likely govern adaptation to this high-elevation habitat. In the Americas, honey bees are hybrids of European and African ancestries and adaptation to latitudinal variation in climate correlates with the proportion of these ancestries across the genome. It is unknown which, if either, of these forms of genetic variation governs adaptation in honey bees living at high elevations in the Americas. Here, we performed whole-genome sequencing of 29 honey bees from both high- and low-elevation populations in Colombia. Analysis of genetic ancestry indicated that both populations were predominantly of African ancestry, but the East African inversions were not detected. However, individuals in the higher elevation population had significantly higher proportions of European ancestry, likely reflecting local adaptation. Several genomic regions exhibited particularly high differentiation between highland and lowland bees, containing candidate loci for local adaptation. Genes that were highly differentiated between highland and lowland populations were enriched for functions related to reproduction and sperm competition. Furthermore, variation in levels of European ancestry across the genome was correlated between populations of honey bees in the highland population and populations at higher latitudes in South America. The results are consistent with the hypothesis that adaptation to both latitude and elevation in these hybrid honey bees are mediated by variation in ancestry at many loci across the genome.

Tailored machine learning models for functional RNA detection in genome-wide screens

The in silico prediction of non-coding and protein-coding genetic loci has received considerable attention in comparative genomics aiming in particular at the identification of properties of nucleotide sequences that are informative of their biological role in the cell. We present here a software framework for the alignment-based training, evaluation and application of machine learning models with user-defined parameters. Instead of focusing on the one-size-fits-all approach of pervasive in silico annotation pipelines, we offer a framework for the structured generation and evaluation of models based on arbitrary features and input data, focusing on stable and explainable results. Furthermore, we showcase the usage of our software package in a full-genome screen of Drosophila melanogaster and evaluate our results against the well-known but much less flexible program RNAz.

In-silico identification and comparison of transcription factor binding sites cluster in anterior-posterior patterning genes in Drosophila melanogaster and Tribolium castaneum

The cis-regulatory data that help in transcriptional regulation is arranged into modular pieces of a few hundred base pairs called CRMs (cis-regulatory modules) and numerous binding sites for multiple transcription factors are prominent characteristics of these cis-regulatory modules. The present study was designed to localize transcription factor binding site (TFBS) clusters on twelve Anterior-posterior (A-P) genes in Tribolium castaneum and compare them to their orthologous gene enhancers in Drosophila melanogaster. Out of the twelve A-P patterning genes, six were gap genes (Kruppel, Knirps, Tailless, Hunchback, Giant, and Caudal) and six were pair rule genes (Hairy, Runt, Even-skipped, Fushi-tarazu, Paired, and Odd-skipped). The genes along with 20 kb upstream and downstream regions were scanned for TFBS clusters using the Motif Cluster Alignment Search Tool (MCAST), a bioinformatics tool that looks for set of nucleotide sequences for statistically significant clusters of non-overlapping occurrence of a given set of motifs. The motifs used in the current study were Hunchback, Caudal, Giant, Kruppel, Knirps, and Even-skipped. The results of the MCAST analysis revealed the maximum number of TFBS for Hunchback, Knirps, Caudal, and Kruppel in both D. melanogaster and T. castaneum, while Bicoid TFBS clusters were found only in D. melanogaster. The size of all the predicted TFBS clusters was less than 1kb in both insect species. These sequences revealed more transversional sites (Tv) than transitional sites (Ti) and the average Ti/Tv ratio was 0.75.

cyjShiny: A cytoscape.js R Shiny Widget for network visualization and analysis

cyjShiny is an open-source R package that allows users to embed network visualization into Shiny apps and R Markdown documents. cyjShiny (https://github.com/cytoscape/cyjShiny) builds on the cytoscape.js Javascript graph library. Additionally, the package provides helper functions to convert common R data representations (e.g., data.frame) into forms compatible with cytoscape.js.

Single-cell profiling of the developing embryonic heart in Drosophila

Drosophila is an important model for studying heart development and disease. Yet, single-cell transcriptomic data of its developing heart have not been performed. Here, we report single-cell profiling of the entire fly heart using ∼3000 Hand-GFP embryos collected at five consecutive developmental stages, ranging from bilateral migrating rows of cardiac progenitors to a fused heart tube. The data revealed six distinct cardiac cell types in the embryonic fly heart: cardioblasts, both Svp+ and Tin+ subtypes; and five types of pericardial cell (PC) that can be distinguished by four key transcription factors (Eve, Odd, Ct and Tin) and include the newly described end of the line PC. Notably, the embryonic fly heart combines transcriptional signatures of the mammalian first and second heart fields. Using unique markers for each heart cell type, we defined their number and location during heart development to build a comprehensive 3D cell map. These data provide a resource to track the expression of any gene in the developing fly heart, which can serve as a reference to study genetic perturbations and cardiac diseases.

Y-Linked Copy Number Polymorphism of Target of Rapamycin Is Associated with Sexual Size Dimorphism in Seed Beetles

The Y chromosome is theorized to facilitate evolution of sexual dimorphism by accumulating sexually antagonistic loci, but empirical support is scarce. Due to the lack of recombination, Y chromosomes are prone to degenerative processes, which poses a constraint on their adaptive potential. Yet, in the seed beetle, Callosobruchus maculatus segregating Y linked variation affects male body size and thereby sexual size dimorphism (SSD). Here, we assemble C. maculatus sex chromosome sequences and identify molecular differences associated with Y-linked SSD variation. The assembled Y chromosome is largely euchromatic and contains over 400 genes, many of which are ampliconic with a mixed autosomal and X chromosome ancestry. Functional annotation suggests that the Y chromosome plays important roles in males beyond primary reproductive functions. Crucially, we find that, besides an autosomal copy of the gene target of rapamycin (TOR), males carry an additional TOR copy on the Y chromosome. TOR is a conserved regulator of growth across taxa, and our results suggest that a Y-linked TOR provides a male specific opportunity to alter body size. A comparison of Y haplotypes associated with male size difference uncovers a copy number variation for TOR, where the haplotype associated with decreased male size, and thereby increased sexual dimorphism, has two additional TOR copies. This suggests that sexual conflict over growth has been mitigated by autosome to Y translocation of TOR followed by gene duplications. Our results reveal that despite of suppressed recombination, the Y chromosome can harbor adaptive potential as a male-limited supergene.

The genetic basis of variation in immune defense against Lysinibacillus fusiformis infection in Drosophila melanogaster

The genetic causes of phenotypic variation often differ depending on the population examined, particularly if the populations were founded by relatively small numbers of genotypes. Similarly, the genetic causes of phenotypic variation among similar traits (resistance to different xenobiotic compounds or pathogens) may also be completely different or only partially overlapping. Differences in genetic causes for variation in the same trait among populations suggests context dependence for how selection acts on those traits. Similarities in the genetic causes of variation for different traits, on the other hand, suggests pleiotropy which would also influence how natural selection shapes variation in a trait. We characterized immune defense against a natural Drosophila pathogen, the Gram-positive bacterium Lysinibacillus fusiformis, in three different populations and found almost no overlap in the genetic architecture of variation in survival post infection. However, when comparing our results to a similar experiment with the fungal pathogen, B. bassiana, we found a convincing shared QTL peak for both pathogens. This peak contains the Bomanin cluster of Drosophila immune effectors. Loss of function mutants and RNAi knockdown experiments confirms a role of some of these genes in immune defense against both pathogens. This suggests that natural selection may act on the entire cluster of Bomanin genes (and the linked region under the QTL) or specific peptides for specific pathogens.