FlyBase: introduction of the Drosophila melanogaster Release 6 reference genome assembly and large-scale migration of genome annotations

Chromosome-level genome assembly of the synanthropic fly Chrysomya megacephala: insights into oviposition location

The oriental latrine fly, Chrysomya megacephala (Diptera: Calliphoridae), is a medically important synanthropic blow fly species characterized by its necrophagy and coprophagy, often observed near carrion and animal feces. Notably, C. megacephala always arrives at carcass earlier than other species. To elucidate the underlying mechanisms behind the host choice in C. megacephala, we present the chromosome-scale genome assembly for this species. The genome size is 816.79 Mb, with a contig N50 of 1.60 Mb. The Hi-C data were anchored to six chromosomes, accounting for 99.93% of the draft assembled genome. Comparative genomic analysis revealed significant expansions in pathways of ligand-gated ion channel activity, passive transmembrane transporter activity, and protein methyltransferase activity, which may be closely associated with host localization and oviposition. After identifying 69 odor-binding proteins (OBPs) in the assembled genome, phylogenetic analysis showed that DmelOBP99b and CmegOBP99b exhibited high homology. Transcriptome analysis demonstrated that the relative expression of CmegOBP99b was consistently the highest during the metamorphosis, and RT-qPCR further confirmed the similar results. Additionally, CmegOBP99b exhibited a strong binding affinity to DMDS (dimethyl disulfide) as determined by molecular docking. To determine the protein expression level of CmegOBP99b in various body parts, we prepared recombinant CmegOBP99b protein and anti-CmegOBP99b polyclonal antibodies. Western blot analysis showed that CmegOBP99b was significantly expressed in the female's head compared to other parts, which is consistent with RT-qPCR results. Therefore, CmegOBP99b may be the primary odor-binding protein responsible for olfactory recognition and the behavioral coordination of C. megacephala. This study not only provides valuable insights into the molecular mechanisms of oviposition localization in C. megacephala but also facilitates further research into the genetic diversity and phylogeny of the Calliphoridae family.

Isolating selective from non-selective forces using site frequency ratios

A new method is introduced for estimating the distribution of mutation fitness effects using site frequency spectra. Unlike previous methods, which make assumptions about non-selective factors, or that try to incorporate such factors into the underlying model, this new method mostly avoids non-selective effects by working with the ratios of counts of selected sites to neutral sites. An expression for the likelihood of a set of selected/neutral ratios is found by treating the ratio of two Poisson random variables as the ratio of two gaussian random variables. This approach also avoids the need to estimate the relative mutation rates of selected and neutral sites. Simulations over a wide range of demographic models, with linked selection effects show that the new SFRatios method performs well for statistical tests of selection, and it performs well for estimating the distribution of selection effects. Performance was better with weak selection models and for expansion and structured demographic models than for bottleneck models. Applications to two populations of Drosophila melanogaster reveal clear but very weak selection on synonymous sites. For nonsynonymous sites, selection was found to be consistent with previous estimates and stronger for an African population than for one from North Carolina.

Behavioral Manipulation of Ixodes scapularis by Ehrlichia muris eauclairensis: Implications for Tick-Borne Disease Transmission

Tick-borne diseases pose significant risks to both animals and humans, with emerging pathogens like Ehrlichia muris eauclairensis (EME) underscoring the need for a deeper understanding of pathogen-vector interactions and tick fitness. This study investigates the impact of EME on Ixodes scapularis nymphs, revealing significant behavioral changes in EME-positive ticks. These ticks exhibited increased movement speed, faster bite site-seeking for attachment, and prolonged feeding durations compared to control ticks. Proteomic analyses of the tick synganglion during resting and feeding phases identified 196 differentially expressed proteins in EME-positive ticks, including multiple proteins associated with nicotinic acetylcholine signaling pathways. Our findings indicated altered neuropeptide expression related to stimulus response and activity, suggesting changes in neurophysiology. This research provides the first evidence of behavioral manipulation by an Ehrlichia species, indicating that the tick nervous system is a site of bacterial influence and a potential target for interventions. These findings offer new insights into pathogen-vector dynamics that could lead to the development of transmission-blocking therapies, significantly impacting tick fitness and disease transmission.

fos genes in mainly invertebrate model systems: A review of commonalities and some diversities

fos genes, transcription factors with a common basic region and leucine zipper domains binding to a consensus DNA sequence (TGA{}TCA), are evolutionarily conserved in eukaryotes. Homologs can be found in many different species from yeast to vertebrates. In yeast, the homologous GCN4 gene is required to mediate "emergency" situations like nutrient deprivation and the unfolded protein response. The C. elegans homolog fos-1 is required for reproduction and vulval development, as well as in adult homeostasis. In Drosophila melanogaster, there is also a sole fos homolog: the gene kayak, with five isoforms. The kayak locus has been studied in detail. It was originally described as embryonic lethal with a "dorsal open" phenotype. Since then, kayak has been shown to be required for oocyte maturation and as a source for piRNA; for early dorsoventral specification, macrophage function, dorsal closure, endoderm differentiation, and finally during metamorphosis in wing and eye-antennal development. In mammals there are multiple fos loci, each one with alternative splicing giving rise to multiple isoforms. Overall, mammalian fos genes are required for bone, cartilage and tooth formation, and in some instances for placental angiogenesis and retinal function. We review here mainly what is known about fos function in invertebrate model systems, especially during embryogenesis. We propose that fos genes, evolutionarily conserved transcription factors, evolved early during eukaryotic development, and from its inception as part of an environmental stress response machinery, were co-opted several times during development to regulate processes that may require similar cellular responses.

The chromosome-level genome provides insights into the adaptive evolution of the visual system in Oratosquilla oratoria

BACKGROUND

The marine crustacean Oratosquilla oratoria is economically significant in seafood and aquaculture industries. However, the lack of high-quality genome assembly has hindered our understanding of O. oratoria, particularly the mechanisms underlying its developed visual system.

RESULTS

We generated a chromosome-level genome assembly for O. oratoria (2.97 Gb, 44 pseudo-chromosomes) using combination sequencing strategies. Our analysis revealed that more than half of the genome was covered by repeat sequences, and LINE elements showed significant expansion. Furthermore, phylogenetic analysis revealed a close relationship of O. oratoria with Dendrobranchiata and Pleocyemata. In addition, the evolutionary rate of O. oratoria was slightly faster than that of Dendrobranchiata and slower than that of Pleocyemata. Interestingly, we observed the significant expansion of middle-wavelength-sensitive (MWS) opsins in the O. oratoria genome by tandem duplication, partially contributing to their unique visual capabilities. Compared with other crustaceans, O. oratoria has evolved a thicker cornea that was possibly driven by visual adaptations and ecological requirements. Employing comparative transcriptome analysis, we identified a tandemly duplicated cuticle protein (CP) cluster that was specifically expanded and expressed in the ocular tissues of O. oratoria, potentially contributing to the thick cornea of O. oratoria.

CONCLUSIONS

Our study established the first chromosome-level genome for Stomatopoda species, providing a valuable genomic resource for studying the molecular mechanisms underlying the developed visual system of O. oratoria.

Female germline expression of OVO transcription factor bridges Drosophila generations

OVO is required for female germ cell viability but has no known function in the male germline in Drosophila. ovo is autoregulated by 2 antagonistic isoforms, OVO-A and OVO-B. All ovo- alleles were created as partial revertants of the antimorphic ovoD1 allele. Creation of new targeted alleles in an ovo+ background indicated that disrupting the germline-specific exon extension of ovo-B leads to an arrested egg chamber phenotype, rather than germ cell death. RNA sequencing analysis, including >1 K full-length cDNAs, indicates that ovo has several unannotated splice variations in the extended exon and a minor population of ovo-B transcripts has an alternative splice. This indicates that classical ovo alleles, such as ovoD1rv23, are not truly null for ovo and are likely to be weak antimorphs. To generate bonafide nulls, we deleted the ovo-A and ovo-B promoters showing that only ovo-B is required for female germ cell viability, and there is an early and continual developmental requirement for ovo-B in the female germline. To visualize OVO expression and localization, we endogenously tagged ovo and found nuclear OVO in all differentiating female germ cells throughout oogenesis in adults. We also found that OVO is maternally deposited into the embryo, where it showed nuclear localization in newly formed pole cells. Maternal OVO persisted in embryonic germ cells until zygotic OVO expression was detectable, suggesting that there is continuous nuclear OVO expression in the female germline in the transition from one generation to the next.

The S6 kinase gene in the fruit fly, Drosophila melanogaster, is essential for metabolic regulation

The S6 kinase (S6K) enzyme phosphorylates the S6 ribosomal protein, promoting protein translation and growth. Here we investigated in flies whether hypomorphic conditions in S6K affect intermediate metabolism and oxidative homeostasis, besides carbohydrates and growth. We also employed partial activation of the nuclear factor 2 erythroid related factor 2 (Nrf2) in a S6K hypomorphic background and controls. S6K is activated by the target of rapamycin (TOR) kinase, a key kinase regulating metabolism, downstream of the insulin receptor in flies. The insulin pathway is a general anabolic pathway, and key regulator of glucose homeostasis. The Nrf2 counters pro-oxidative conditions, also involved in inflammatory responses and metabolism. The Nrf2 fly homolog is Cap'n'collar C (CncC). We quantified glucose, glycogen, and total lipids in control and different pro-oxidative conditions. We corroborated an accumulation of lipids and carbohydrates in the mutants, and document sexual differences. We document also metabolic and survival differences in the responses to "mild" pro-oxidative conditions in young flies (seven days old), with females being most affected. We compare 10 mM paraquat survival of virgin flies to mated mixed-sex flies housed together. We used females to study transcriptomic differences between wild type and S6k hypomorphs. Results highlight dysregulation of lipid and antioxidant enzymes and genes, in agreement with lipid and oxidative metabolism data. Our results are consistent with the insulin/TOR pathway acting as an integrator of intermediate metabolism and oxidative homeostasis (this last together with the CncC pathway).

Superior target genes and pathways for RNAi-mediated pest control revealed by genome-wide analysis in the beetle Tribolium castaneum

BACKGROUND

An increasing human population, the emergence of resistances against pesticides and their potential impact on the environment call for the development of new eco-friendly pest control strategies. RNA interference (RNAi)-based pesticides have emerged as a new option with the first products entering the market. Essentially, double-stranded RNAs targeting essential genes of pests are either expressed in the plants or sprayed on their surface. Upon feeding, pests mount an RNAi response and die. However, it has remained unclear whether RNAi-based insecticides should target the same pathways as classic pesticides or whether the different mode-of-action would favor other processes. Moreover, there is no consensus on the best genes to be targeted.

RESULTS

We performed a genome-wide screen in the red flour beetle to identify 905 RNAi target genes. Based on a validation screen and clustering, we identified the 192 most effective target genes in that species. The transfer to oral application in other beetle pests revealed a list of 34 superior target genes, which are an excellent starting point for application in other pests. Gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) analyses of our genome-wide dataset revealed that genes with high efficacy belonged mainly to basic cellular processes such as gene expression and protein homeostasis - processes not targeted by classic insecticides.

CONCLUSION

Our work revealed the best target genes and target processes for RNAi-based pest control and we propose a procedure to transfer our short list of superior target genes to other pests. © 2024 The Author(s). Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

TimeFlies: an snRNA-seq aging clock for the fruit fly head sheds light on sex-biased aging

Although multiple high-performing epigenetic aging clocks exist, few are based directly on gene expression. Such transcriptomic aging clocks allow us to extract age-associated genes directly. However, most existing transcriptomic clocks model a subset of genes and are limited in their ability to predict novel biomarkers. With the growing popularity of single-cell sequencing, there is a need for robust single-cell transcriptomic aging clocks. Moreover, clocks have yet to be applied to investigate the elusive phenomenon of sex differences in aging. We introduce TimeFlies, a pan-cell-type scRNA-seq aging clock for the Drosophila melanogaster head. TimeFlies uses deep learning to classify the donor age of cells based on genome-wide gene expression profiles. Using explainability methods, we identified key marker genes contributing to the classification, with lncRNAs showing up as highly enriched among predicted biomarkers. The top biomarker gene across cell types is lncRNA:roX1, a regulator of X chromosome dosage compensation, a pathway previously identified as a top biomarker of aging in the mouse brain. We validated this finding experimentally, showing a decrease in survival probability in the absence of roX1 in vivo. Furthermore, we trained sex-specific TimeFlies clocks and noted significant differences in model predictions and explanations between male and female clocks, suggesting that different pathways drive aging in males and females.

Double trouble: two retrotransposons triggered a cascade of invasions in Drosophila species within the last 50 years

Horizontal transfer of genetic material in eukaryotes has rarely been documented over short evolutionary timescales. Here, we show that two retrotransposons, Shellder and Spoink, invaded the genomes of multiple species of the melanogaster subgroup within the last 50 years. Through horizontal transfer, Spoink spread in D. melanogaster during the 1980s, while both Shellder and Spoink invaded D. simulans in the 1990s. Possibly following hybridization, D. simulans infected the island endemic species D. mauritiana (Mauritius) and D. sechellia (Seychelles) with both TEs after 1995. In the same approximate time-frame, Shellder also invaded D. teissieri, a species confined to sub-Saharan Africa. We find that the donors of Shellder and Spoink are likely American Drosophila species from the willistoni, cardini, and repleta groups. Thus, the described cascade of TE invasions could only become feasible after D. melanogaster and D. simulans extended their distributions into the Americas 200 years ago, likely aided by human activity. Our work reveals that cascades of TE invasions, likely initiated by human-mediated range expansions, could have an impact on the genomic and phenotypic evolution of geographically dispersed species. Within a few decades, TEs could invade many species, including island endemics, with distributions very distant from the donor of the TE.

Cis-regulatory Variation in Relation to Sex and Sexual Dimorphism in Drosophila melanogaster

Much of sexual dimorphism is likely due to sex-biased gene expression, which results from differential regulation of a genome that is largely shared between males and females. Here, we use allele-specific expression to explore cis-regulatory variation in Drosophila melanogaster in relation to sex. We develop a Bayesian framework to infer the transcriptome-wide joint distribution of cis-regulatory effects across the sexes. We also examine patterns of cis-regulatory variation with respect to two other levels of variation in sexual dimorphism: (i) across genes that vary in their degree of sex-biased expression and (ii) among tissues that vary in their degree of dimorphism (e.g. relatively low dimorphism in heads vs. high dimorphism in gonads). We uncover evidence of widespread cis-regulatory variation in all tissues examined, with female-biased genes being especially enriched for this variation. A sizeable proportion of cis-regulatory variation is inferred to have sex-specific effects, with sex-dependent cis effects being much more frequent in gonads than in heads. Finally, we find some genes where 1 allele contributes to more than 50% of a gene's expression in heterozygous males but <50% of its expression in heterozygous females. Such variants could provide a mechanism for sex-specific dominance reversals, a phenomenon important for sexually antagonistic balancing selection. However, tissue differences in allelic imbalance are approximately as frequent as sex differences, perhaps suggesting that sexual conflict may not be particularly unique in shaping patterns of expression variation.

OVO positively regulates essential maternal pathways by binding near the transcriptional start sites in the Drosophila female germline

Differentiation of female germline stem cells into a mature oocyte includes the expression of RNAs and proteins that drive early embryonic development in Drosophila. We have little insight into what activates the expression of these maternal factors. One candidate is the zinc-finger protein OVO. OVO is required for female germline viability and has been shown to positively regulate its own expression, as well as a downstream target, ovarian tumor, by binding to the transcriptional start site (TSS). To find additional OVO targets in the female germline and further elucidate OVO's role in oocyte development, we performed ChIP-seq to determine genome-wide OVO occupancy, as well as RNA-seq comparing hypomorphic and wild type rescue ovo alleles. OVO preferentially binds in close proximity to target TSSs genome-wide, is associated with open chromatin, transcriptionally active histone marks, and OVO-dependent expression. Motif enrichment analysis on OVO ChIP peaks identified a 5'-TAACNGT-3' OVO DNA binding motif spatially enriched near TSSs. However, the OVO DNA binding motif does not exhibit precise motif spacing relative to the TSS characteristic of RNA polymerase II complex binding core promoter elements. Integrated genomics analysis showed that 525 genes that are bound and increase in expression downstream of OVO are known to be essential maternally expressed genes. These include genes involved in anterior/posterior/germ plasm specification (bcd, exu, swa, osk, nos, aub, pgc, gcl), egg activation (png, plu, gnu, wisp, C(3)g, mtrm), translational regulation (cup, orb, bru1, me31B), and vitelline membrane formation (fs(1)N, fs(1)M3, clos). This suggests that OVO is a master transcriptional regulator of oocyte development and is responsible for the expression of structural components of the egg as well as maternally provided RNAs that are required for early embryonic development.

Transcriptome analysis of Drosophila suzukii reveals molecular mechanisms conferring pyrethroid and spinosad resistance

Drosophila suzukii lay eggs in soft-skinned, ripening fruits, making this insect a serious threat to berry production. Since its 2008 introduction into North America, growers have used insecticides, such as pyrethroids and spinosads, as the primary approach for D. suzukii management, resulting in development of insecticide resistance in this pest. This study sought to identify the molecular mechanisms conferring insecticide resistance in these populations. We sequenced the transcriptomes of two pyrethroid- and two spinosad-resistant isofemale lines. In both pyrethroid-resistant lines and one spinosad-resistant line, we identified overexpression of metabolic genes that are implicated in resistance in other insect pests. In the other spinosad-resistant line, we observed an overexpression of cuticular genes that have been linked to resistance. Our findings enabled the development of molecular diagnostics that we used to confirm persistence of insecticide resistance in California, U.S.A. To validate these findings, we leveraged D. melanogaster mutants with reduced expression of metabolic or cuticular genes that were found to be upregulated in resistant D. suzukii to demonstrate that these genes are involved in promoting resistance. This study is the first to characterize the molecular mechanisms of insecticide resistance in D. suzukii and provides insights into how current management practices can be optimized.

Drosophila melanogaster pigmentation demonstrates adaptive phenotypic parallelism but genomic unpredictability over multiple timescales

Populations are capable of responding to environmental change over ecological timescales via adaptive tracking. However, the translation from patterns of allele frequency change to rapid adaptation of complex traits remains unresolved. We used abdominal pigmentation in Drosophila melanogaster as a model phenotype to address the nature, genetic architecture, and repeatability of rapid adaptation in the field. We show that D. melanogaster pigmentation evolves as a highly parallel and deterministic response to shared environmental gradients across latitude and season in natural North American populations. We then experimentally evolved replicate, genetically diverse fly populations in field mesocosms to remove any confounding effects of demography and/or cryptic structure that may drive patterns in wild populations; we show that pigmentation rapidly responds, in parallel, in fewer than ten generations. Thus, pigmentation evolves concordantly in response to spatial and temporal climatic gradients. We next examined whether phenotypic differentiation was associated with allele frequency change at loci with established links to genetic variance in pigmentation in natural populations. We found that across all spatial and temporal scales, phenotypic patterns were associated with variation at pigmentation-related loci, and the sets of genes we identified in each context were largely nonoverlapping. Therefore, our findings suggest that parallel phenotypic evolution is associated with an unpredictable genomic response, with distinct components of the polygenic architecture shifting across each environmental gradient to produce redundant adaptive patterns.

Significance Statement

Shifts in global climate conditions have heightened our need to understand the dynamics and pace of adaptation in natural populations. In order to anticipate the population-level response to rapidly changing environmental conditions, we need to understand whether trait evolution is predictable over short timescales, and whether the genetic basis of adaptation is shared or distinct across multiple timescales. Here, we explored parallelism in the adaptive response of a complex phenotype, D. melanogaster pigmentation, to shared conditions that varied over multiple spatiotemporal scales. Our results demonstrate that while phenotypic adaptation proceeds as a predictable response to environmental gradients, even over short timescales, the genetic basis of the adaptive response is variable and nuanced across spatial and temporal contexts.

Chromosome structure in Drosophila is determined by boundary pairing not loop extrusion

Two different models have been proposed to explain how the endpoints of chromatin looped domains ('TADs') in eukaryotic chromosomes are determined. In the first, a cohesin complex extrudes a loop until it encounters a boundary element roadblock, generating a stem-loop. In this model, boundaries are functionally autonomous: they have an intrinsic ability to halt the movement of incoming cohesin complexes that is independent of the properties of neighboring boundaries. In the second, loops are generated by boundary:boundary pairing. In this model, boundaries are functionally non-autonomous, and their ability to form a loop depends upon how well they match with their neighbors. Moreover, unlike the loop-extrusion model, pairing interactions can generate both stem-loops and circle-loops. We have used a combination of MicroC to analyze how TADs are organized, and experimental manipulations of the even skipped TAD boundary, homie, to test the predictions of the 'loop-extrusion' and the 'boundary-pairing' models. Our findings are incompatible with the loop-extrusion model, and instead suggest that the endpoints of TADs in flies are determined by a mechanism in which boundary elements physically pair with their partners, either head-to-head or head-to-tail, with varying degrees of specificity. Although our experiments do not address how partners find each other, the mechanism is unlikely to require loop extrusion.

Stem-loop and circle-loop TADs generated by directional pairing of boundary elements have distinct physical and regulatory properties

The chromosomes in multicellular eukaryotes are organized into a series of topologically independent loops called TADs. In flies, TADs are formed by physical interactions between neighboring boundaries. Fly boundaries exhibit distinct partner preferences, and pairing interactions between boundaries are typically orientation-dependent. Pairing can be head-to-tail or head-to-head. The former generates a stem-loop TAD, while the latter gives a circle-loop TAD. The TAD that encompasses the Drosophila even skipped (eve) gene is formed by the head-to-tail pairing of the nhomie and homie boundaries. To explore the relationship between loop topology and the physical and regulatory landscape, we flanked the nhomie boundary region with two attP sites. The attP sites were then used to generate four boundary replacements: λ DNA, nhomie forward (WT orientation), nhomie reverse (opposite of WT orientation), and homie forward (same orientation as WT homie). The nhomie forward replacement restores the WT physical and regulatory landscape: in MicroC experiments, the eve TAD is a 'volcano' triangle topped by a plume, and the eve gene and its regulatory elements are sequestered from interactions with neighbors. The λ DNA replacement lacks boundary function: the endpoint of the 'new' eve TAD on the nhomie side is ill-defined, and eve stripe enhancers activate a nearby gene, eIF3j. While nhomie reverse and homie forward restore the eve TAD, the topology is a circle-loop, and this changes the local physical and regulatory landscape. In MicroC experiments, the eve TAD interacts with its neighbors, and the plume at the top of the eve triangle peak is converted to a pair of 'clouds' of contacts with the next-door TADs. Consistent with the loss of isolation afforded by the stem-loop topology, the eve enhancers weakly activate genes in the neighboring TADs. Conversely, eve function is partially disrupted.

The Drosophila histone methyltransferase SET1 coordinates multiple signaling pathways in regulating male germline stem cell maintenance and differentiation

Many tissue-specific adult stem cell lineages maintain a balance between proliferation and differentiation. Here, we study how the H3K4me3 methyltransferase Set1 regulates early-stage male germ cells in Drosophila. Early-stage germline-specific knockdown of Set1 results in temporally progressive defects, arising as germ cell loss and developing into overpopulated early-stage germ cells. These germline defects also impact the niche architecture and cyst stem cell lineage non-cell-autonomously. Additionally, wild-type Set1, but not the catalytically inactive Set1, rescues the Set1 knockdown phenotypes, highlighting the functional importance of the methyltransferase activity of Set1. Further, RNA-sequencing experiments reveal key signaling pathway components, such as the JAK-STAT pathway gene Stat92E and the BMP pathway gene Mad, which are upregulated upon Set1 knockdown. Genetic interaction assays support the functional relationships between Set1 and JAK-STAT or BMP pathways, as both Stat92E and Mad mutations suppress the Set1 knockdown phenotypes. These findings enhance our understanding of the balance between proliferation and differentiation in an adult stem cell lineage. The phenotype of germ cell loss followed by over-proliferation when inhibiting a histone methyltransferase also raises concerns about using their inhibitors in cancer therapy.

Biases in ARG-Based Inference of Historical Population Size in Populations Experiencing Selection

Inferring the demographic history of populations provides fundamental insights into species dynamics and is essential for developing a null model to accurately study selective processes. However, background selection and selective sweeps can produce genomic signatures at linked sites that mimic or mask signals associated with historical population size change. While the theoretical biases introduced by the linked effects of selection have been well established, it is unclear whether ancestral recombination graph (ARG)-based approaches to demographic inference in typical empirical analyses are susceptible to misinference due to these effects. To address this, we developed highly realistic forward simulations of human and Drosophila melanogaster populations, including empirically estimated variability of gene density, mutation rates, recombination rates, purifying, and positive selection, across different historical demographic scenarios, to broadly assess the impact of selection on demographic inference using a genealogy-based approach. Our results indicate that the linked effects of selection minimally impact demographic inference for human populations, although it could cause misinference in populations with similar genome architecture and population parameters experiencing more frequent recurrent sweeps. We found that accurate demographic inference of D. melanogaster populations by ARG-based methods is compromised by the presence of pervasive background selection alone, leading to spurious inferences of recent population expansion, which may be further worsened by recurrent sweeps, depending on the proportion and strength of beneficial mutations. Caution and additional testing with species-specific simulations are needed when inferring population history with non-human populations using ARG-based approaches to avoid misinference due to the linked effects of selection.

Single-fly genome assemblies fill major phylogenomic gaps across the Drosophilidae Tree of Life

Long-read sequencing is driving rapid progress in genome assembly across all major groups of life, including species of the family Drosophilidae, a longtime model system for genetics, genomics, and evolution. We previously developed a cost-effective hybrid Oxford Nanopore (ONT) long-read and Illumina short-read sequencing approach and used it to assemble 101 drosophilid genomes from laboratory cultures, greatly increasing the number of genome assemblies for this taxonomic group. The next major challenge is to address the laboratory culture bias in taxon sampling by sequencing genomes of species that cannot easily be reared in the lab. Here, we build upon our previous methods to perform amplification-free ONT sequencing of single wild flies obtained either directly from the field or from ethanol-preserved specimens in museum collections, greatly improving the representation of lesser studied drosophilid taxa in whole-genome data. Using Illumina Novaseq X Plus and ONT P2 sequencers with R10.4.1 chemistry, we set a new benchmark for inexpensive hybrid genome assembly at US $150 per genome while assembling genomes from as little as 35 ng of genomic DNA from a single fly. We present 183 new genome assemblies for 179 species as a resource for drosophilid systematics, phylogenetics, and comparative genomics. Of these genomes, 62 are from pooled lab strains and 121 from single adult flies. Despite the sample limitations of working with small insects, most single-fly diploid assemblies are comparable in contiguity (>1 Mb contig N50), completeness (>98% complete dipteran BUSCOs), and accuracy (>QV40 genome-wide with ONT R10.4.1) to assemblies from inbred lines. We present a well-resolved multi-locus phylogeny for 360 drosophilid and 4 outgroup species encompassing all publicly available (as of August 2023) genomes for this group. Finally, we present a Progressive Cactus whole-genome, reference-free alignment built from a subset of 298 suitably high-quality drosophilid genomes. The new assemblies and alignment, along with updated laboratory protocols and computational pipelines, are released as an open resource and as a tool for studying evolution at the scale of an entire insect family.

Odorant receptor co-receptors affect expression of tuning receptors in Drosophila

Insects detect odorants using two large families of heteromeric receptors, the Odorant Receptors (ORs) and Ionotropic Receptors (IRs). Most OR and IR genes encode odorant-binding "tuning" subunits, whereas four (Orco, Ir8a, Ir25a, and Ir76b) encode co-receptor subunits required for receptor function. Olfactory neurons are thought to degenerate in the absence of Orco in ants and bees, and limited data suggest this may happen to some olfactory neurons in Drosophila fruit flies as well. Here, we thoroughly examined the role of co-receptors on olfactory neuron survival in Drosophila. Leveraging knowledge that olfactory neuron classes are defined by the expression of different tuning receptors, we used tuning receptor expression in antennal transcriptomes as a proxy for the survival of distinct olfactory neuron classes. Consistent with olfactory neuron degeneration, expression of many OR-family tuning receptors is decreased in Orco mutants relative to controls, and transcript loss is progressive with age. The effects of Orco are highly receptor-dependent, with expression of some receptor transcripts nearly eliminated and others unaffected. Surprisingly, further studies revealed that olfactory neuron classes with reduced tuning receptor expression generally survive in Orco mutant flies. Furthermore, there is little apoptosis or neuronal loss in the antenna of these flies. We went on to investigate the effects of IR family co-receptor mutants using similar approaches and found that expression of IR tuning receptors is decreased in the absence of Ir8a and Ir25a, but not Ir76b. As in Orco mutants, Ir8a-dependent olfactory neurons mostly endure despite near-absent expression of associated tuning receptors. Finally, we used differential expression analysis to identify other antennal genes whose expression is changed in IR and OR co-receptor mutants. Taken together, our data indicate that odorant co-receptors are necessary for maintaining expression of many tuning receptors at the mRNA level. Further, most Drosophila olfactory neurons persist in OR and IR co-receptor mutants, suggesting that the impact of co-receptors on neuronal survival may vary across insect species.

OVO Positively Regulates Essential Maternal Pathways by Binding Near the Transcriptional Start Sites in the Drosophila Female Germline

Differentiation of female germline stem cells into a mature oocyte includes the expression of RNAs and proteins that drive early embryonic development in Drosophila. We have little insight into what activates the expression of these maternal factors. One candidate is the zinc-finger protein OVO. OVO is required for female germline viability and has been shown to positively regulate its own expression, as well as a downstream target, ovarian tumor, by binding to the transcriptional start site (TSS). To find additional OVO targets in the female germline and further elucidate OVO's role in oocyte development, we performed ChIP-seq to determine genome-wide OVO occupancy, as well as RNA-seq comparing hypomorphic and wild type rescue ovo alleles. OVO preferentially binds in close proximity to target TSSs genome-wide, is associated with open chromatin, transcriptionally active histone marks, and OVO-dependent expression. Motif enrichment analysis on OVO ChIP peaks identified a 5'-TAACNGT-3' OVO DNA binding motif spatially enriched near TSSs. However, the OVO DNA binding motif does not exhibit precise motif spacing relative to the TSS characteristic of RNA Polymerase II complex binding core promoter elements. Integrated genomics analysis showed that 525 genes that are bound and increase in expression downstream of OVO are known to be essential maternally expressed genes. These include genes involved in anterior/posterior/germ plasm specification (bcd, exu, swa, osk, nos, aub, pgc, gcl), egg activation (png, plu, gnu, wisp, C(3)g, mtrm), translational regulation (cup, orb, bru1, me31B), and vitelline membrane formation (fs(1)N, fs(1)M3, clos). This suggests that OVO is a master transcriptional regulator of oocyte development and is responsible for the expression of structural components of the egg as well as maternally provided RNAs that are required for early embryonic development.

Earl Grey: A Fully Automated User-Friendly Transposable Element Annotation and Analysis Pipeline

Transposable elements (TEs) are major components of eukaryotic genomes and are implicated in a range of evolutionary processes. Yet, TE annotation and characterization remain challenging, particularly for nonspecialists, since existing pipelines are typically complicated to install, run, and extract data from. Current methods of automated TE annotation are also subject to issues that reduce overall quality, particularly (i) fragmented and overlapping TE annotations, leading to erroneous estimates of TE count and coverage, and (ii) repeat models represented by short sections of total TE length, with poor capture of 5' and 3' ends. To address these issues, we present Earl Grey, a fully automated TE annotation pipeline designed for user-friendly curation and annotation of TEs in eukaryotic genome assemblies. Using nine simulated genomes and an annotation of Drosophila melanogaster, we show that Earl Grey outperforms current widely used TE annotation methodologies in ameliorating the issues mentioned above while scoring highly in benchmarking for TE annotation and classification and being robust across genomic contexts. Earl Grey provides a comprehensive and fully automated TE annotation toolkit that provides researchers with paper-ready summary figures and outputs in standard formats compatible with other bioinformatics tools. Earl Grey has a modular format, with great scope for the inclusion of additional modules focused on further quality control and tailored analyses in future releases.

The evolution of sexual dimorphism in gene expression in response to a manipulation of mate competition

Many genes are differentially expressed between males and females and patterns of sex-biased gene expression (SBGE) vary among species. Some of this variation is thought to have evolved in response to differences in mate competition among species that cause varying patterns of sex-specific selection. We used experimental evolution to test this by quantifying SBGE and sex-specific splicing in 15 Drosophila melanogaster populations that evolved for 104 generations in mating treatments that removed mate competition via enforced monogamy, or allowed mate competition in either small, simple, or larger, structurally more complex mating environments. Consistent with sex-specific selection affecting SBGE, initially sex-biased genes diverged in expression more among treatments than unbiased genes, and there was greater expression divergence for male- than female-biased genes. It has been suggested the transcriptome should be "feminized" under monogamy because of the removal of sexual selection on males; we did not observe this, likely because selection differs in additional ways between monogamy vs. polygamy. Significant divergence in average expression dimorphism between treatments was observed and, in some treatment comparisons, the direction of the divergence differed across different sex-bias categories. There was not a generalized reduction in expression dimorphism under enforced monogamy.

scRNA-seq Reveals Novel Genetic Pathways and Sex Chromosome Regulation in Tribolium Spermatogenesis

Spermatogenesis is critical to sexual reproduction yet evolves rapidly in many organisms. High-throughput single-cell transcriptomics promises unparalleled insight into this important process but understanding can be impeded in nonmodel systems by a lack of known genes that can reliably demarcate biologically meaningful cell populations. Tribolium castaneum, the red flour beetle, lacks known markers for spermatogenesis found in insect species like Drosophila melanogaster. Using single-cell sequencing data collected from adult beetle testes, we implement a strategy for elucidating biologically meaningful cell populations by using transient expression stage identification markers, weighted principal component clustering, and SNP-based haploid/diploid phasing. We identify populations that correspond to observable points in sperm differentiation and find species specific markers for each stage. Our results indicate that molecular pathways underlying spermatogenesis in Coleoptera are substantially diverged from those in Diptera. We also show that most genes on the X chromosome experience meiotic sex chromosome inactivation. Temporal expression of Drosophila MSL complex homologs coupled with spatial analysis of potential chromatin entry sites further suggests that the dosage compensation machinery may mediate escape from meiotic sex chromosome inactivation and postmeiotic reactivation of the X chromosome.

Genetic basis of variation in thermal developmental plasticity for Drosophila melanogaster body pigmentation

Seasonal differences in insect pigmentation are attributed to the influence of ambient temperature on pigmentation development. This thermal plasticity is adaptive and heritable, and thereby capable of evolving. However, the specific genes contributing to the variation in plasticity that can drive its evolution remain largely unknown. To address this, we analysed pigmentation and pigmentation plasticity in Drosophila melanogaster. We measured two components of pigmentation in the thorax and abdomen: overall darkness and the proportion of length covered by darker pattern elements (a trident in the thorax and bands in the abdomen) in females from two developmental temperatures (17 or 28°C) and 191 genotypes. Using a GWAS approach to identify the genetic basis of variation in pigmentation and its response to temperature, we identified numerous dispersed QTLs, including some mapping to melanogenesis genes (yellow, ebony, and tan). Remarkably, we observed limited overlap between QTLs for variation within specific temperatures and those influencing thermal plasticity, as well as minimal overlap between plasticity QTLs across pigmentation components and across body parts. For most traits, consistent with selection favouring the retention of plasticity, we found that lower plasticity alleles were often at lower frequencies. The functional analysis of selected candidate QTLs and pigmentation genes largely confirmed their contributions to variation in pigmentation and/or pigmentation plasticity. Overall, our study reveals the existence and underlying basis of extensive and trait-specific genetic variation for pigmentation and pigmentation plasticity, offering a rich reservoir of raw material for natural selection to shape the evolution of these traits independently.

The Drosophila histone methyl-transferase SET1 coordinates multiple signaling pathways in regulating male germline stem cell maintenance and differentiation

Many cell types come from tissue-specific adult stem cells that maintain the balance between proliferation and differentiation. Here, we study how the H3K4me3 methyltransferase, Set1, regulates early-stage male germ cell proliferation and differentiation in Drosophila. Early-stage germline-specific knockdown of set1 results in a temporally progressed defects, arising as germ cell loss and developing to overpopulated early-stage germ cells. These germline defects also impact the niche architecture and cyst stem cell lineage in a non-cell-autonomous manner. Additionally, wild-type Set1, but not the catalytically inactive Set1, could rescue the set1 knockdown phenotypes, highlighting the functional importance of the methyl-transferase activity of the Set1 enzyme. Further, RNA-seq experiments reveal key signaling pathway components, such as the JAK-STAT pathway gene stat92E and the BMP pathway gene mad, that are upregulated upon set1 knockdown. Genetic interaction assays support the functional relationships between set1 and JAK-STAT or BMP pathways, as mutations of both the stat92E and mad genes suppress the set1 knockdown phenotypes. These findings enhance our understanding of the balance between proliferation and differentiation in an adult stem cell lineage. The germ cell loss followed by over-proliferation phenotypes when inhibiting a histone methyl-transferase raise concerns about using their inhibitors in cancer therapy.

Hippocampome.org 2.0 is a knowledge base enabling data-driven spiking neural network simulations of rodent hippocampal circuits

Hippocampome.org is a mature open-access knowledge base of the rodent hippocampal formation focusing on neuron types and their properties. Previously, Hippocampome.org v1.0 established a foundational classification system identifying 122 hippocampal neuron types based on their axonal and dendritic morphologies, main neurotransmitter, membrane biophysics, and molecular expression (Wheeler et al., 2015). Releases v1.1 through v1.12 furthered the aggregation of literature-mined data, including among others neuron counts, spiking patterns, synaptic physiology, in vivo firing phases, and connection probabilities. Those additional properties increased the online information content of this public resource over 100-fold, enabling numerous independent discoveries by the scientific community. Hippocampome.org v2.0, introduced here, besides incorporating over 50 new neuron types, now recenters its focus on extending the functionality to build real-scale, biologically detailed, data-driven computational simulations. In all cases, the freely downloadable model parameters are directly linked to the specific peer-reviewed empirical evidence from which they were derived. Possible research applications include quantitative, multiscale analyses of circuit connectivity and spiking neural network simulations of activity dynamics. These advances can help generate precise, experimentally testable hypotheses and shed light on the neural mechanisms underlying associative memory and spatial navigation.

Transvection between nonallelic genomic positions in Drosophila

In Drosophila, pairing of maternal and paternal homologous chromosomes can permit trans-interactions between enhancers on one homolog and promoters on another, an example of transvection. Although trans-interactions have been observed at many loci in the Drosophila genome and in other organisms, the parameters that govern enhancer action in trans remain poorly understood. Using a transgenic reporter system, we asked whether enhancers and promoters at nonallelic, but nearby, genomic positions can communication in trans. Using one transgenic insertion carrying the synthetic enhancer GMR and another nearby insertion carrying the hsp70 promoter driving a fluorescent reporter, we show that transgenes separated by 2.6 kb of linear distance can support enhancer action in trans at the 53F8 locus. Furthermore, transvection between the nonallelic insertions can be augmented by a small deletion flanking one insert, likely via changes to the paired configuration of the homologs. Subsequent analyses of other insertions in 53F8 that carry different transgenic sequences demonstrate that the capacity to support transvection between nonallelic sites varies greatly, suggesting that factors beyond the linear distance between insertion sites play an important role. Finally, analysis of transvection between nearby nonallelic sites at other genomic locations shows evidence of position effects, where one locus supported GMR action in trans over a linear distance of over 10 kb, whereas another locus showed no evidence of transvection over a span <200 bp. Overall, our data demonstrate that transvection between nonallelic sites represents a complex interplay between genomic context, interallelic distance, and promoter identity.

Essential genes identification model based on sequence feature map and graph convolutional neural network

BACKGROUND

Essential genes encode functions that play a vital role in the life activities of organisms, encompassing growth, development, immune system functioning, and cell structure maintenance. Conventional experimental techniques for identifying essential genes are resource-intensive and time-consuming, and the accuracy of current machine learning models needs further enhancement. Therefore, it is crucial to develop a robust computational model to accurately predict essential genes.

RESULTS

In this study, we introduce GCNN-SFM, a computational model for identifying essential genes in organisms, based on graph convolutional neural networks (GCNN). GCNN-SFM integrates a graph convolutional layer, a convolutional layer, and a fully connected layer to model and extract features from gene sequences of essential genes. Initially, the gene sequence is transformed into a feature map using coding techniques. Subsequently, a multi-layer GCN is employed to perform graph convolution operations, effectively capturing both local and global features of the gene sequence. Further feature extraction is performed, followed by integrating convolution and fully-connected layers to generate prediction results for essential genes. The gradient descent algorithm is utilized to iteratively update the cross-entropy loss function, thereby enhancing the accuracy of the prediction results. Meanwhile, model parameters are tuned to determine the optimal parameter combination that yields the best prediction performance during training.

CONCLUSIONS

Experimental evaluation demonstrates that GCNN-SFM surpasses various advanced essential gene prediction models and achieves an average accuracy of 94.53%. This study presents a novel and effective approach for identifying essential genes, which has significant implications for biology and genomics research.

The ftz upstream element drives late ftz stripes but is not required for regulation of Ftz target genes

The regulation of gene expression in precise, rapidly changing spatial patterns is essential for embryonic development. Multiple enhancers have been identified for the evolving expression patterns of the cascade of Drosophila segmentation genes that establish the basic body plan of the fly. Classic reporter transgene experiments identified multiple cis-regulatory elements (CREs) that are sufficient to direct various aspects of the evolving expression pattern of the pair-rule gene fushi tarazu (ftz). These include enhancers that coordinately activate expression in all seven stripes and stripe-specific elements that activate expression in one or more ftz stripes. Of the two 7-stripe enhancers, analysis of reporter transgenes demonstrated that the upstream element (UPS) is autoregulatory, requiring direct binding of Ftz protein to direct striped expression. Here, we asked about the endogenous role of the UPS by precisely deleting this 7-stripe enhancer. In ftzΔUPS7S homozygotes, ftz stripes appear in the same order as wildtype, and all but stripe 4 are expressed at wildtype levels by the end of the cellular blastoderm stage. This suggests that the zebra element and UPS harbor information to direct stripe 4 expression, although previous deletion analyses failed to identify a stripe-specific CRE within these two 7-stripe enhancers. However, the UPS is necessary for late ftz stripe expression, with all 7 stripes decaying earlier than wildtype in ftzΔUPS7S homozygotes. Despite this premature loss of ftz expression, downstream target gene regulation proceeds as in wildtype, and segmentation is unperturbed in the overwhelming majority of animals. We propose that this late-acting enhancer provides a buffer against perturbations in gene expression but is not required for establishment of Ftz cell fates. Overall, our results demonstrate that multiple enhancers, each directing distinct aspects of an overall gene expression pattern, contribute to fine-tuning the complex patterns necessary for embryonic development.

Decoupled evolution of the Sex Peptide gene family and Sex Peptide Receptor in Drosophilidae

Across internally fertilising species, males transfer ejaculate proteins that trigger wide-ranging changes in female behaviour and physiology. Much theory has been developed to explore the drivers of ejaculate protein evolution. The accelerating availability of high-quality genomes now allows us to test how these proteins are evolving at fine taxonomic scales. Here, we use genomes from 264 species to chart the evolutionary history of Sex Peptide (SP), a potent regulator of female post-mating responses in Drosophila melanogaster. We infer that SP first evolved in the Drosophilinae subfamily and has followed markedly different evolutionary trajectories in different lineages. Outside of the Sophophora-Lordiphosa, SP exists largely as a single-copy gene with independent losses in several lineages. Within the Sophophora-Lordiphosa, the SP gene family has repeatedly and independently expanded. Up to seven copies, collectively displaying extensive sequence variation, are present in some species. Despite these changes, SP expression remains restricted to the male reproductive tract. Alongside, we document considerable interspecific variation in the presence and morphology of seminal microcarriers that, despite the critical role SP plays in microcarrier assembly in D. melanogaster, appear to be independent of changes in the presence/absence or sequence of SP. We end by providing evidence that SP's evolution is decoupled from that of its receptor, SPR, in which we detect no evidence of correlated diversifying selection. Collectively, our work describes the divergent evolutionary trajectories that a novel gene has taken following its origin and finds a surprisingly weak coevolutionary signal between a supposedly sexually antagonistic protein and its receptor.

Paralog transcriptional differentiation in the D. melanogaster-specific gene family Sdic across populations and spermatogenesis stages

How recently originated gene copies become stable genomic components remains uncertain as high sequence similarity of young duplicates precludes their functional characterization. The tandem multigene family Sdic is specific to Drosophila melanogaster and has been annotated across multiple reference-quality genome assemblies. Here we show the existence of a positive correlation between Sdic copy number and total expression, plus vast intrastrain differences in mRNA abundance among paralogs, using RNA-sequencing from testis of four strains with variable paralog composition. Single cell and nucleus RNA-sequencing data expose paralog expression differentiation in meiotic cell types within testis from third instar larva and adults. Additional RNA-sequencing across synthetic strains only differing in their Y chromosomes reveal a tissue-dependent trans-regulatory effect on Sdic: upregulation in testis and downregulation in male accessory gland. By leveraging paralog-specific expression information from tissue- and cell-specific data, our results elucidate the intraspecific functional diversification of a recently expanded tandem gene family.

The composition of piRNA clusters in Drosophila melanogaster deviates from expectations under the trap model

BACKGROUND

It is widely assumed that the invasion of a transposable element (TE) in mammals and invertebrates is stopped when a copy of the TE jumps into a piRNA cluster (i.e., the trap model). However, recent works, which for example showed that deletion of three major piRNA clusters has no effect on TE activity, cast doubt on the trap model.

RESULTS

Here, we test the trap model from a population genetics perspective. Our simulations show that the composition of regions that act as transposon traps (i.e., potentially piRNA clusters) ought to deviate from regions that have no effect on TE activity. We investigated TEs in five Drosophila melanogaster strains using three complementary approaches to test whether the composition of piRNA clusters matches these expectations. We found that the abundance of TE families inside and outside of piRNA clusters is highly correlated, although this is not expected under the trap model. Furthermore, the distribution of the number of TE insertions in piRNA clusters is also much broader than expected.

CONCLUSIONS

We found that the observed composition of piRNA clusters is not in agreement with expectations under the simple trap model. Dispersed piRNA producing TE insertions and temporal as well as spatial heterogeneity of piRNA clusters may account for these deviations.

ChimeraTE: a pipeline to detect chimeric transcripts derived from genes and transposable elements

Transposable elements (TEs) produce structural variants and are considered an important source of genetic diversity. Notably, TE-gene fusion transcripts, i.e. chimeric transcripts, have been associated with adaptation in several species. However, the identification of these chimeras remains hindered due to the lack of detection tools at a transcriptome-wide scale, and to the reliance on a reference genome, even though different individuals/cells/strains have different TE insertions. Therefore, we developed ChimeraTE, a pipeline that uses paired-end RNA-seq reads to identify chimeric transcripts through two different modes. Mode 1 is the reference-guided approach that employs canonical genome alignment, and Mode 2 identifies chimeras derived from fixed or insertionally polymorphic TEs without any reference genome. We have validated both modes using RNA-seq data from four Drosophila melanogaster wild-type strains. We found ∼1.12% of all genes generating chimeric transcripts, most of them from TE-exonized sequences. Approximately ∼23% of all detected chimeras were absent from the reference genome, indicating that TEs belonging to chimeric transcripts may be recent, polymorphic insertions. ChimeraTE is the first pipeline able to automatically uncover chimeric transcripts without a reference genome, consisting of two running Modes that can be used as a tool to investigate the contribution of TEs to transcriptome plasticity.

Single-fly assemblies fill major phylogenomic gaps across the Drosophilidae Tree of Life

Long-read sequencing is driving rapid progress in genome assembly across all major groups of life, including species of the family Drosophilidae, a longtime model system for genetics, genomics, and evolution. We previously developed a cost-effective hybrid Oxford Nanopore (ONT) long-read and Illumina short-read sequencing approach and used it to assemble 101 drosophilid genomes from laboratory cultures, greatly increasing the number of genome assemblies for this taxonomic group. The next major challenge is to address the laboratory culture bias in taxon sampling by sequencing genomes of species that cannot easily be reared in the lab. Here, we build upon our previous methods to perform amplification-free ONT sequencing of single wild flies obtained either directly from the field or from ethanol-preserved specimens in museum collections, greatly improving the representation of lesser studied drosophilid taxa in whole-genome data. Using Illumina Novaseq X Plus and ONT P2 sequencers with R10.4.1 chemistry, we set a new benchmark for inexpensive hybrid genome assembly at US $150 per genome while assembling genomes from as little as 35 ng of genomic DNA from a single fly. We present 183 new genome assemblies for 179 species as a resource for drosophilid systematics, phylogenetics, and comparative genomics. Of these genomes, 62 are from pooled lab strains and 121 from single adult flies. Despite the sample limitations of working with small insects, most single-fly diploid assemblies are comparable in contiguity (>1Mb contig N50), completeness (>98% complete dipteran BUSCOs), and accuracy (>QV40 genome-wide with ONT R10.4.1) to assemblies from inbred lines. We present a well-resolved multi-locus phylogeny for 360 drosophilid and 4 outgroup species encompassing all publicly available (as of August 2023) genomes for this group. Finally, we present a Progressive Cactus whole-genome, reference-free alignment built from a subset of 298 suitably high-quality drosophilid genomes. The new assemblies and alignment, along with updated laboratory protocols and computational pipelines, are released as an open resource and as a tool for studying evolution at the scale of an entire insect family.

Analysis of asymptomatic Drosophila models for ALS and SMA reveals convergent impact on functional protein complexes linked to neuro-muscular degeneration

BACKGROUND

Spinal Muscular Atrophy (SMA) and Amyotrophic Lateral Sclerosis (ALS) share phenotypic and molecular commonalities, including the fact that they can be caused by mutations in ubiquitous proteins involved in RNA metabolism, namely SMN, TDP-43 and FUS. Although this suggests the existence of common disease mechanisms, there is currently no model to explain the resulting motor neuron dysfunction. In this work we generated a parallel set of Drosophila models for adult-onset RNAi and tagged neuronal expression of the fly orthologues of the three human proteins, named Smn, TBPH and Caz, respectively. We profiled nuclear and cytoplasmic bound mRNAs using a RIP-seq approach and characterized the transcriptome of the RNAi models by RNA-seq. To unravel the mechanisms underlying the common functional impact of these proteins on neuronal cells, we devised a computational approach based on the construction of a tissue-specific library of protein functional modules, selected by an overall impact score measuring the estimated extent of perturbation caused by each gene knockdown.

RESULTS

Transcriptome analysis revealed that the three proteins do not bind to the same RNA molecules and that only a limited set of functionally unrelated transcripts is commonly affected by their knock-down. However, through our integrative approach we were able to identify a concerted effect on protein functional modules, albeit acting through distinct targets. Most strikingly, functional annotation revealed that these modules are involved in critical cellular pathways for motor neurons, including neuromuscular junction function. Furthermore, selected modules were found to be significantly enriched in orthologues of human neuronal disease genes.

CONCLUSIONS

The results presented here show that SMA and ALS disease-associated genes linked to RNA metabolism functionally converge on neuronal protein complexes, providing a new hypothesis to explain the common motor neuron phenotype. The functional modules identified represent promising biomarkers and therapeutic targets, namely given their alteration in asymptomatic settings.

The genomic distribution of transposable elements is driven by spatially variable purifying selection

It is widely accepted that the genomic distribution of transposable elements (TEs) mainly reflects the outcome of purifying selection and insertion bias (1). Nevertheless, the relative importance of these two evolutionary forces could not be tested thoroughly. Here, we introduce an experimental system, which allows separating purifying selection from TE insertion bias. We used experimental evolution to study the TE insertion patterns in Drosophila simulans founder populations harboring 1040 insertions of an active P-element. After 10 generations at a large population size, we detected strong selection against P-element insertions. The exception were P-element insertions in genomic regions for which a strong insertion bias has been proposed (2-4). Because recurrent P-element insertions cannot explain this pattern, we conclude that purifying selection, with variable strength along the chromosomes, is the major determinant of the genomic distribution of P-elements. Genomic regions with relaxed purifying selection against P-element insertions exhibit normal levels of purifying selection against base substitutions. This suggests that different types of purifying selection operate on base substitutions and P-element insertions. Our results highlight the power of experimental evolution to understand basic evolutionary processes, which are difficult to infer from patterns of natural variation alone.

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

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

Female germline expression of OVO transcription factor bridges Drosophila generations

OVO is required for karyotypically female germ cell viability but has no known function in the male germline in Drosophila. ovo is autoregulated by two antagonistic isoforms, OVO-A and OVO-B. All ovo- alleles were created as partial revertants of the antimorphic ovoD1 allele. Creation of new targeted alleles in an ovo+ background indicated that disrupting the germline-specific exon extension of ovo-B leads to an arrested egg chamber phenotype, rather than germ cell death. RNA-seq analysis, including >1K full length cDNAs, indicates that ovo utilizes a number of unannotated splice variations in the extended exon and a minor population of ovo-B transcripts utilizes an alternative splice. This indicates that classical ovo alleles such as ovoD1rv23, are not truly null for ovo, and are likely to be weak antimorphs. To generate bonafide nulls, we deleted the ovo-A and ovo-B promoters showing that only ovo-B is required for female germ cell viability and there is an early and polyphasic developmental requirement for ovo-B in the female germline. To visualize OVO expression and localization, we endogenously tagged ovo and found nuclear OVO in all differentiating female germ cells throughout oogenesis in adults. We also found that OVO is maternally deposited into the embryo, where it showed nuclear localization in newly formed pole cells. Maternal OVO persisted in embryonic germ cells until zygotic OVO expression was detectable, suggesting that there is continuous nuclear OVO expression in the female germline in the transition from one generation to the next.

lncRNA read-through regulates the BX-C insulator Fub-1

Though long non-coding RNAs (lncRNAs) represent a substantial fraction of the Pol II transcripts in multicellular animals, only a few have known functions. Here we report that the blocking activity of the Bithorax complex (BX-C) Fub-1 boundary is segmentally regulated by its own lncRNA. The Fub-1 boundary is located between the Ultrabithorax (Ubx) gene and the bxd/pbx regulatory domain, which is responsible for regulating Ubx expression in parasegment PS6/segment A1. Fub-1 consists of two hypersensitive sites, HS1 and HS2. HS1 is an insulator while HS2 functions primarily as an lncRNA promoter. To activate Ubx expression in PS6/A1, enhancers in the bxd/pbx domain must be able to bypass Fub-1 blocking activity. We show that the expression of the Fub-1 lncRNAs in PS6/A1 from the HS2 promoter inactivates Fub-1 insulating activity. Inactivation is due to read-through as the HS2 promoter must be directed toward HS1 to disrupt blocking.

The Transmission Patterns of the Endosymbiont Wolbachia within the Hawaiian Drosophilidae Adaptive Radiation

The evolution of endosymbionts and their hosts can lead to highly dynamic interactions with varying fitness effects for both the endosymbiont and host species. Wolbachia, a ubiquitous endosymbiont of arthropods and nematodes, can have both beneficial and detrimental effects on host fitness. We documented the occurrence and patterns of transmission of Wolbachia within the Hawaiian Drosophilidae and examined the potential contributions of Wolbachia to the rapid diversification of their hosts. Screens for Wolbachia infections across a minimum of 140 species of Hawaiian Drosophila and Scaptomyza revealed species-level infections of 20.0%, and across all 399 samples, a general infection rate of 10.3%. Among the 44 Wolbachia strains we identified using a modified Wolbachia multi-locus strain typing scheme, 30 (68.18%) belonged to supergroup B, five (11.36%) belonged to supergroup A, and nine (20.45%) had alleles with conflicting supergroup assignments. Co-phylogenetic reconciliation analysis indicated that Wolbachia strain diversity within their endemic Hawaiian Drosophilidae hosts can be explained by vertical (e.g., co-speciation) and horizontal (e.g., host switch) modes of transmission. Results from stochastic character trait mapping suggest that horizontal transmission is associated with the preferred oviposition substrate of the host, but not the host's plant family or island of occurrence. For Hawaiian Drosophilid species of conservation concern, with 13 species listed as endangered and 1 listed as threatened, knowledge of Wolbachia strain types, infection status, and potential for superinfection could assist with conservation breeding programs designed to bolster population sizes, especially when wild populations are supplemented with laboratory-reared, translocated individuals. Future research aimed at improving the understanding of the mechanisms of Wolbachia transmission in nature, their impact on the host, and their role in host species formation may shed light on the influence of Wolbachia as an evolutionary driver, especially in Hawaiian ecosystems.

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

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

Integrating lipid metabolism, pheromone production and perception by Fruitless and Hepatocyte Nuclear Factor 4

Sexual attraction and perception are crucial for mating and reproductive success. In Drosophila melanogaster, the male-specific isoform of Fruitless (Fru), FruM, is a known master neuro-regulator of innate courtship behavior to control the perception of sex pheromones in sensory neurons. Here, we show that the non-sex-specific Fru isoform (FruCOM) is necessary for pheromone biosynthesis in hepatocyte-like oenocytes for sexual attraction. Loss of FruCOM in oenocytes resulted in adults with reduced levels of cuticular hydrocarbons (CHCs), including sex pheromones, and show altered sexual attraction and reduced cuticular hydrophobicity. We further identify Hepatocyte nuclear factor 4 (Hnf4) as a key target of FruCOM in directing fatty acid conversion to hydrocarbons. Fru or Hnf4 depletion in oenocytes disrupts lipid homeostasis, resulting in a sex-dimorphic CHC profile that differs from doublesex- and transformer-dependent CHC dimorphism. Thus, Fru couples pheromone perception and production in separate organs to regulate chemosensory communications and ensure efficient mating behavior.

Genomic benchmarks: a collection of datasets for genomic sequence classification

BACKGROUND

Recently, deep neural networks have been successfully applied in many biological fields. In 2020, a deep learning model AlphaFold won the protein folding competition with predicted structures within the error tolerance of experimental methods. However, this solution to the most prominent bioinformatic challenge of the past 50 years has been possible only thanks to a carefully curated benchmark of experimentally predicted protein structures. In Genomics, we have similar challenges (annotation of genomes and identification of functional elements) but currently, we lack benchmarks similar to protein folding competition.

RESULTS

Here we present a collection of curated and easily accessible sequence classification datasets in the field of genomics. The proposed collection is based on a combination of novel datasets constructed from the mining of publicly available databases and existing datasets obtained from published articles. The collection currently contains nine datasets that focus on regulatory elements (promoters, enhancers, open chromatin region) from three model organisms: human, mouse, and roundworm. A simple convolution neural network is also included in a repository and can be used as a baseline model. Benchmarks and the baseline model are distributed as the Python package 'genomic-benchmarks', and the code is available at https://github.com/ML-Bioinfo-CEITEC/genomic_benchmarks .

CONCLUSIONS

Deep learning techniques revolutionized many biological fields but mainly thanks to the carefully curated benchmarks. For the field of Genomics, we propose a collection of benchmark datasets for the classification of genomic sequences with an interface for the most commonly used deep learning libraries, implementation of the simple neural network and a training framework that can be used as a starting point for future research. The main aim of this effort is to create a repository for shared datasets that will make machine learning for genomics more comparable and reproducible while reducing the overhead of researchers who want to enter the field, leading to healthy competition and new discoveries.

Activation of the Cap'n'collar C pathway (Nrf2 pathway in vertebrates) signaling in insulin pathway compromised Drosophila melanogaster flies ameliorates the diabetic state upon pro-oxidant conditions

The insulin pathway is a crucial central system for metabolism and growth. The Nrf2 signaling pathway functions to counteract oxidative stress. Here we sought to study the consequences of an oxidative stress challenge to insulin compromised and control adult flies of different ages, varying the activation state of the Nrf2 pathway in flies, the Cap'n'collar C pathway. For this, we employed two different pro-oxidative conditions: 3 % hydrogen peroxide or 20 mM paraquat laced in the food. In both cases, wild type (control) flies die within a few days, yet there are significant differences between males and females, and also within flies of different ages (seven versus thirty days old flies). We repeated the same conditions with young (seven days old) flies that were heterozygous for a loss-of-function mutation in Keap1. There were no significant differences. We then tested two hypomorphic viable conditions of the insulin pathway (heteroallelic combination for the insulin receptor and the S6 Kinase), challenged in the same way: Whereas they also die in the pro-oxidant conditions, they fare significantly better when heterozygous for Keap1, in contrast to controls. We also monitored locomotion in all of these conditions, and, in general, found significant differences between flies without and with a mutant allele (heterozygous) for Keap1. Our results point to altered oxidative stress conditions in diabetic flies. These findings suggest that modest activation of the Cap'n'collar C pathway may be a treatment for diabetic symptoms.

Heterozygous inversion breakpoints suppress meiotic crossovers by altering recombination repair outcomes

Heterozygous chromosome inversions suppress meiotic crossover (CO) formation within an inversion, potentially because they lead to gross chromosome rearrangements that produce inviable gametes. Interestingly, COs are also severely reduced in regions nearby but outside of inversion breakpoints even though COs in these regions do not result in rearrangements. Our mechanistic understanding of why COs are suppressed outside of inversion breakpoints is limited by a lack of data on the frequency of noncrossover gene conversions (NCOGCs) in these regions. To address this critical gap, we mapped the location and frequency of rare CO and NCOGC events that occurred outside of the dl-49 chrX inversion in D. melanogaster. We created full-sibling wildtype and inversion stocks and recovered COs and NCOGCs in the syntenic regions of both stocks, allowing us to directly compare rates and distributions of recombination events. We show that COs outside of the proximal inversion breakpoint are distributed in a distance-dependent manner, with strongest suppression near the inversion breakpoint. We find that NCOGCs occur evenly throughout the chromosome and, importantly, are not suppressed near inversion breakpoints. We propose a model in which COs are suppressed by inversion breakpoints in a distance-dependent manner through mechanisms that influence DNA double-strand break repair outcome but not double-strand break formation. We suggest that subtle changes in the synaptonemal complex and chromosome pairing might lead to unstable interhomolog interactions during recombination that permits NCOGC formation but not CO formation.

Hippo pathway and Bonus control developmental cell fate decisions in the Drosophila eye

The canonical function of the Hippo signaling pathway is the regulation of organ growth. How this pathway controls cell-fate determination is less well understood. Here, we identify a function of the Hippo pathway in cell-fate decisions in the developing Drosophila eye, exerted through the interaction of Yorkie (Yki) with the transcriptional regulator Bonus (Bon), an ortholog of mammalian transcriptional intermediary factor 1/tripartite motif (TIF1/TRIM) family proteins. Instead of controlling tissue growth, Yki and Bon promote epidermal and antennal fates at the expense of the eye fate. Proteomic, transcriptomic, and genetic analyses reveal that Yki and Bon control these cell-fate decisions by recruiting transcriptional and post-transcriptional co-regulators and by repressing Notch target genes and activating epidermal differentiation genes. Our work expands the range of functions and regulatory mechanisms under Hippo pathway control.

A genome-wide computational approach to define microRNA-Polycomb/trithorax gene regulatory circuits in Drosophila

Characterization of gene regulatory networks is fundamental to understanding homeostatic development. This process can be simplified by analyzing relatively simple genomes such as the genome of Drosophila melanogaster. In this work we have developed a computational framework in Drosophila to explore for the presence of gene regulatory circuits between two large groups of transcriptional regulators: the epigenetic group of the Polycomb/trithorax (PcG/trxG) proteins and the microRNAs (miRNAs). We have searched genome-wide for miRNA targets in PcG/trxG transcripts as well as for Polycomb Response Elements (PREs) in miRNA genes. Our results show that 10% of the analyzed miRNAs could be controlling PcG/trxG gene expression, while 40% of those miRNAs are putatively controlled by the selected set of PcG/trxG proteins. The integration of these analyses has resulted in the predicted existence of 3 classes of miRNA-PcG/trxG crosstalk interactions that define potential regulatory circuits. In the first class, miRNA-PcG circuits are defined by miRNAs that reciprocally crosstalk with PcG. In the second, miRNA-trxG circuits are defined by miRNAs that reciprocally crosstalk with trxG. In the third class, miRNA-PcG/trxG shared circuits are defined by miRNAs that crosstalk with both PcG and trxG regulators. These putative regulatory circuits may uncover a novel mechanism in Drosophila for the control of PcG/trxG and miRNAs levels of expression. The computational framework developed here for Drosophila melanogaster can serve as a model case for similar analyses in other species. Moreover, our work provides, for the first time, a new and useful resource for the Drosophila community to consult prior to experimental studies investigating the epigenetic regulatory networks of miRNA-PcG/trxG mediated gene expression.

Functionally Validating Evolutionary Conserved Risk Genes for Parkinson's Disease in Drosophila melanogaster

Parkinson's disease (PD) is a heterogeneous and complex neurodegenerative disorder and large-scale genetic studies have identified >130 genes associated with PD. Although genomic studies have been decisive for our understanding of the genetic contributions underlying PD, these associations remain as statistical associations. Lack of functional validation limits the biological interpretation; however, it is labour extensive, expensive, and time consuming. Therefore, the ideal biological system for functionally validating genetic findings must be simple. The study aim was to assess systematically evolutionary conserved PD-associated genes using Drosophila melanogaster. From a literature review, a total of 136 genes have found to be associated with PD in GWAS studies, of which 11 are strongly evolutionary conserved between Homo sapiens and D. melanogaster. By ubiquitous gene expression knockdown of the PD-genes in D. melanogaster, the flies' escape response was investigated by assessing their negative geotaxis response, a phenotype that has previously been used to investigate PD in D. melanogaster. Gene expression knockdown was successful in 9/11 lines, and phenotypic consequences were observed in 8/9 lines. The results provide evidence that genetically modifying expression levels of PD genes in D. melanogaster caused reduced climbing ability of the flies, potentially supporting their role in dysfunctional locomotion, a hallmark of PD.

Apocrine secretion in the salivary glands of Drosophilidae and other dipterans is evolutionarily conserved

Apocrine secretion is a transport and secretory mechanism that remains only partially characterized, even though it is evolutionarily conserved among all metazoans, including humans. The excellent genetic model organism Drosophila melanogaster holds promise for elucidating the molecular mechanisms regulating this fundamental metazoan process. Two prerequisites for such investigations are to clearly define an experimental system to investigate apocrine secretion and to understand the evolutionarily and functional contexts in which apocrine secretion arose in that system. To this end, we recently demonstrated that, in D. melanogaster, the prepupal salivary glands utilize apocrine secretion prior to pupation to deliver innate immune and defense components to the exuvial fluid that lies between the metamorphosing pupae and its chitinous case. This finding provided a unique opportunity to appraise how this novel non-canonical and non-vesicular transport and secretory mechanism is employed in different developmental and evolutionary contexts. Here we demonstrate that this apocrine secretion, which is mechanistically and temporarily separated from the exocytotic mechanism used to produce the massive salivary glue secretion (Sgs), is shared across Drosophilidae and two unrelated dipteran species. Screening more than 30 species of Drosophila from divergent habitats across the globe revealed that apocrine secretion is a widespread and evolutionarily conserved cellular mechanism used to produce exuvial fluid. Species with longer larval and prepupal development than D. melanogaster activate apocrine secretion later, while smaller and more rapidly developing species activate it earlier. In some species, apocrine secretion occurs after the secretory material is first concentrated in cytoplasmic structures of unknown origin that we name "collectors." Strikingly, in contrast to the widespread use of apocrine secretion to provide exuvial fluid, not all species use exocytosis to produce the viscid salivary glue secretion that is seen in D. melanogaster. Thus, apocrine secretion is the conserved mechanism used to realize the major function of the salivary gland in fruitflies and related species: it produces the pupal exuvial fluid that provides an active defense against microbial invasion during pupal metamorphosis.