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

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.

Dietary stress remodels the genetic architecture of lifespan variation in outbred Drosophila

Evolutionary theory suggests that lifespan-reducing alleles should be purged from the gene pool, and yet decades of genome-wide association and model organism studies have shown that they persist. One potential explanation is that alleles that regulate lifespan do so only in certain environmental contexts. We exposed outbred Drosophila to control and high-sugar diets and genotyped more than 10,000 adult flies to track allele frequency changes over the course of a single adult lifespan. We identified thousands of lifespan-associated alleles associated with early versus late-life trade-offs, late-onset effects and genotype-by-environment interactions. Remarkably, a third of lifespan-associated genetic variation had environmentally dependent effects on lifespan. We find that lifespan-reducing alleles are often recently derived, have stronger effects on a high-sugar diet and show signatures of selection in wild Drosophila populations, consistent with the evolutionary mismatch hypothesis. Our results provide insight into the highly polygenic and context-dependent genetic architecture of lifespan variation and the evolutionary processes that shape this key trait.

The path to immortalization of cells starts by managing stress through gene duplications

The genomes of immortalized cell lines (and cancer cells) are characterized by multiple types of aberrations, ranging from single nucleotide polymorphisms (SNPs) to structural rearrangements that have accumulated over time. Consequently, it is difficult to estimate the relative impact of different aberrations, the order of events, and which gene functions were under selective pressure at the early stage towards cellular immortalization. Here, we have established novel cell cultures derived from Drosophila melanogaster embryos that were sampled at multiple time points over a one-year period. Using short-read DNA sequencing, we show that copy-number gain in preferentially stress-related genes were acquired in a dominant fraction of cells in 300-days old cultures. Furthermore, transposable elements were active in cells of all cultures. Only a few (<1%) SNPs could be followed over time, and these showed no trend to increase or decrease. We conclude that the early cellular responses of a novel culture comprise sequence duplication and transposable element activity. During immortalization, positive selection first occurs on genes that are related to stress response before shifting to genes that are related to growth.

MED13 and glycolysis are conserved modifiers of α-synuclein-associated neurodegeneration

α-Synuclein (α-syn) is important in synucleinopathies such as Parkinson's disease (PD). While genome-wide association studies (GWASs) of synucleinopathies have identified many risk loci, the underlying genes have not been shown for most loci. Using Drosophila, we screened 3,471 mutant chromosomes for genetic modifiers of α-synuclein and identified 12 genes. Eleven modifiers have human orthologs associated with diseases, including MED13 and CDC27, which lie within PD GWAS loci. Drosophila Skd/Med13 and glycolytic enzymes are co-upregulated by α-syn-associated neurodegeneration. While elevated α-syn compromises mitochondrial function, co-expressing skd/Med13 RNAi and α-syn synergistically increase the ratio of oxidized-to-reduced glutathione. The resulting neurodegeneration can be suppressed by overexpressing a glycolytic enzyme or treatment with deferoxamine, suggesting that compensatory glycolysis is neuroprotective. In addition, the functional relationship between α-synuclein, MED13, and glycolytic enzymes is conserved between flies and mice. We propose that hypoxia-inducible factor and MED13 are part of a druggable pathway for PD.

Drosophila SUMM4 complex couples insulator function and DNA replication control

Asynchronous replication of chromosome domains during S phase is essential for eukaryotic genome function, but the mechanisms establishing which domains replicate early versus late in different cell types remain incompletely understood. Intercalary heterochromatin domains replicate very late in both diploid chromosomes of dividing cells and in endoreplicating polytene chromosomes where they are also underreplicated. Drosophila SNF2-related factor SUUR imparts locus-specific underreplication of polytene chromosomes. SUUR negatively regulates DNA replication fork progression; however, its mechanism of action remains obscure. Here, we developed a novel method termed MS-Enabled Rapid protein Complex Identification (MERCI) to isolate a stable stoichiometric native complex SUMM4 that comprises SUUR and a chromatin boundary protein Mod(Mdg4)-67.2. Mod(Mdg4) stimulates SUUR ATPase activity and is required for a normal spatiotemporal distribution of SUUR in vivo. SUUR and Mod(Mdg4)-67.2 together mediate the activities of gypsy insulator that prevent certain enhancer-promoter interactions and establish euchromatin-heterochromatin barriers in the genome. Furthermore, SuUR or mod(mdg4) mutations reverse underreplication of intercalary heterochromatin. Thus, SUMM4 can impart late replication of intercalary heterochromatin by attenuating the progression of replication forks through euchromatin/heterochromatin boundaries. Our findings implicate a SNF2 family ATP-dependent motor protein SUUR in the insulator function, reveal that DNA replication can be delayed by a chromatin barrier, and uncover a critical role for architectural proteins in replication control. They suggest a mechanism for the establishment of late replication that does not depend on an asynchronous firing of late replication origins.

Transposable element finder (TEF): finding active transposable elements from next generation sequencing data

BACKGROUND

Detection of newly transposed events by transposable elements (TEs) from next generation sequence (NGS) data is difficult, due to their multiple distribution sites over the genome containing older TEs. The previously reported Transposon Insertion Finder (TIF) detects TE transpositions on the reference genome from NGS short reads using end sequences of target TE. TIF requires the sequence of target TE and is not able to detect transpositions for TEs with an unknown sequence.

RESULT

The new algorithm Transposable Element Finder (TEF) enables the detection of TE transpositions, even for TEs with an unknown sequence. TEF is a finding tool of transposed TEs, in contrast to TIF as a detection tool of transposed sites for TEs with a known sequence. The transposition event is often accompanied with a target site duplication (TSD). Focusing on TSD, two algorithms to detect both ends of TE, TSDs and target sites are reported here. One is based on the grouping with TSDs and direct comparison of k-mers from NGS without similarity search. The other is based on the junction mapping of TE end sequence candidates. Both methods succeed to detect both ends and TSDs of known active TEs in several tests with rice, Arabidopsis and Drosophila data and discover several new TEs in new locations. PCR confirmed the detected transpositions of TEs in several test cases in rice.

CONCLUSIONS

TEF detects transposed TEs with TSDs as a result of TE transposition, sequences of both ends and their inserted positions of transposed TEs by direct comparison of NGS data between two samples. Genotypes of transpositions are verified by counting of junctions of head and tail, and non-insertion sequences in NGS reads. TEF is easy to run and independent of any TE library, which makes it useful to detect insertions from unknown TEs bypassed by common TE annotation pipelines.

m6A modification of U6 snRNA modulates usage of two major classes of pre-mRNA 5' splice site

Alternative splicing of messenger RNAs is associated with the evolution of developmentally complex eukaryotes. Splicing is mediated by the spliceosome, and docking of the pre-mRNA 5' splice site into the spliceosome active site depends upon pairing with the conserved ACAGA sequence of U6 snRNA. In some species, including humans, the central adenosine of the ACAGA box is modified by N6 methylation, but the role of this m6A modification is poorly understood. Here, we show that m6A modified U6 snRNA determines the accuracy and efficiency of splicing. We reveal that the conserved methyltransferase, FIONA1, is required for Arabidopsis U6 snRNA m6A modification. Arabidopsis fio1 mutants show disrupted patterns of splicing that can be explained by the sequence composition of 5' splice sites and cooperative roles for U5 and U6 snRNA in splice site selection. U6 snRNA m6A influences 3' splice site usage. We generalise these findings to reveal two major classes of 5' splice site in diverse eukaryotes, which display anti-correlated interaction potential with U5 snRNA loop 1 and the U6 snRNA ACAGA box. We conclude that U6 snRNA m6A modification contributes to the selection of degenerate 5' splice sites crucial to alternative splicing.

Chronic sleep loss disrupts rhythmic gene expression in Drosophila

Genome-wide profiling of rhythmic gene expression has offered new avenues for studying the contribution of circadian clock to diverse biological processes. Sleep has been considered one of the most important physiological processes that are regulated by the circadian clock, however, the effects of chronic sleep loss on rhythmic gene expression remain poorly understood. In the present study, we exploited Drosophila sleep mutants insomniac 1 (inc 1 ) and wide awake D2 (wake D2 ) as models for chronic sleep loss. We profiled the transcriptomes of head tissues collected from 4-week-old wild type flies, inc 1 and wake D2 at timepoints around the clock. Analysis of gene oscillation revealed a substantial loss of rhythmicity in inc 1 and wake D2 compared to wild type flies, with most of the affected genes common to both mutants. The disruption of gene oscillation was not due to changes in average gene expression levels. We also identified a subset of genes whose loss of rhythmicity was shared among animals with chronic sleep loss and old flies, suggesting a contribution of aging to chronic, sleep-loss-induced disruption of gene oscillation.

Deep learning of cross-species single-cell landscapes identifies conserved regulatory programs underlying cell types

Despite extensive efforts to generate and analyze reference genomes, genetic models to predict gene regulation and cell fate decisions are lacking for most species. Here, we generated whole-body single-cell transcriptomic landscapes of zebrafish, Drosophila and earthworm. We then integrated cell landscapes from eight representative metazoan species to study gene regulation across evolution. Using these uniformly constructed cross-species landscapes, we developed a deep-learning-based strategy, Nvwa, to predict gene expression and identify regulatory sequences at the single-cell level. We systematically compared cell-type-specific transcription factors to reveal conserved genetic regulation in vertebrates and invertebrates. Our work provides a valuable resource and offers a new strategy for studying regulatory grammar in diverse biological systems.

Vitelline Membrane Protein 26 Mutagenesis, Using CRISPR/Cas9, Results in Egg Collapse in Plutella xylostella

Vitelline membrane proteins (VMPs) are the main proteins that form the inner shell (vitelline membrane layer) of insect eggs and are an integral part of egg formation and embryo development. Here, we characterized the molecular structure and expression patterns of the VMP26 gene and analyzed its reproductive functions in diamondback moth, Plutella xylostella (L.), a worldwide migratory pest of cruciferous plants. The PxVMP26 gene was shown to be a single exon gene that contained an open reading frame of 852 base pairs (bp) encoding 283 amino acids. Both qPCR and western blot analyses showed that PxVMP26 was specifically expressed in female adults and was significantly highly expressed in the ovary. Further anatomical analysis indicated that the expression level of PxVMP26 in the ovarian tube with an incomplete yolk was significantly higher than that in the ovarian tube with a complete yolk. CRISPR/Cas9-induced PxVMP26 knockout successfully created two homozygous strains with 8- and 46-bp frameshift mutations. The expression deficiency of the PxVMP26 protein was detected in the mutant strains using immunofluorescence and western blot. No significant difference was found in the number of eggs laid within three days between wild and mutant individuals, but there was a lower egg hatchability. The loss of the PxVMP26 gene changed the mean egg size, damaged the structure of the vitelline membrane, and increased the proportion of abnormal eggs due to water loss, resulting in egg collapse. This first analysis of the roles of the VMP gene in the oocyte formation and embryonic development of P. xylostella, using CRISPR/Cas9 technology, provides a basis for screening new genetic control targets of P. xylostella.

Altering enhancer-promoter linear distance impacts promoter competition in cis and in trans

In Drosophila, pairing of maternal and paternal homologs can permit trans-interactions between enhancers on one homolog and promoters on another, an example of a phenomenon called transvection. When chromosomes are paired, promoters in cis and in trans to an enhancer can compete for the enhancer's activity, but the parameters that govern this competition are as yet poorly understood. To assess how the linear spacing between an enhancer and promoter can influence promoter competition in Drosophila, we employed transgenic constructs wherein the eye-specific enhancer GMR is placed at varying distances from a heterologous hsp70 promoter driving a fluorescent reporter. While GMR activates the reporter to a high degree when the enhancer and promoter are spaced by a few hundred base pairs, activation is strongly attenuated when the enhancer is moved 3 kilobases away. By examining transcription of endogenous genes near the point of transgene insertion, we show that linear spacing of 3 kb between GMR and the hsp70 promoter results in elevated transcription of neighboring promoters, suggesting a loss of specificity between the enhancer and its intended transgenic target promoter. Furthermore, increasing spacing between GMR and hsp70 by just 100 bp can enhance transvection, resulting in increased activation of a promoter on a paired homolog at the expense of a promoter in cis to the enhancer. Finally, cis-/trans-promoter competition assays in which one promoter carries mutations to key core promoter elements show that GMR will skew its activity toward a wild type promoter, suggesting that an enhancer is in a balanced competition between its potential target promoters in cis and in trans.

ZipHiC: a novel Bayesian framework to identify enriched interactions and experimental biases in Hi-C data

Motivation

Several computational and statistical methods have been developed to analyze data generated through the 3C-based methods, especially the Hi-C. Most of the existing methods do not account for dependency in Hi-C data.

Results

Here, we present ZipHiC, a novel statistical method to explore Hi-C data focusing on the detection of enriched contacts. ZipHiC implements a Bayesian method based on a hidden Markov random field (HMRF) model and the Approximate Bayesian Computation (ABC) to detect interactions in two-dimensional space based on a Hi-C contact frequency matrix. ZipHiC uses data on the sources of biases related to the contact frequency matrix, allows borrowing information from neighbours using the Potts model and improves computation speed using the ABC model. In addition to outperforming existing tools on both simulated and real data, our model also provides insights into different sources of biases that affects Hi-C data. We show that some datasets display higher biases from DNA accessibility or Transposable Elements content. Furthermore, our analysis in Drosophila melanogaster showed that approximately half of the detected significant interactions connect promoters with other parts of the genome indicating a functional biological role. Finally, we found that the micro-C datasets display higher biases from DNA accessibility compared to a similar Hi-C experiment, but this can be corrected by ZipHiC.

Availability and implementation

The R scripts are available at https://github.com/igosungithub/HMRFHiC.git.

Supplementary information

Supplementary data are available at Bioinformatics online.

Internal oligo(dT) priming introduces systematic bias in bulk and single-cell RNA sequencing count data

Significant advances in RNA sequencing have been recently made possible by using oligo(dT) primers for simultaneous mRNA enrichment and reverse transcription priming. The associated increase in efficiency has enabled more economical bulk RNA sequencing methods and the advent of high-throughput single-cell RNA sequencing, already one of the most widely adopted methods in transcriptomics. However, the effects of off-target oligo(dT) priming on gene expression quantification have not been appreciated. In the present study, we describe the extent, the possible causes, and the consequences of internal oligo(dT) priming across multiple public datasets obtained from various bulk and single-cell RNA sequencing platforms. To explore and address this issue, we developed a computational algorithm for RNA counting methods, which identifies the sequencing read alignments that likely resulted from internal oligo(dT) priming and removes them from the data. Directly comparing filtered datasets to those obtained by an alternative method reveals significant improvements in gene expression measurement. Finally, we infer a list of human genes whose expression quantification is most likely to be affected by internal oligo(dT) priming and predict that when measured using these methods, the expression of most genes may be inflated by at least 10% whereby some genes are affected more than others.

Comparative transcriptional analysis uncovers molecular processes in early and mature somatic cyst cells of Drosophila testes

The tight interaction between somatic and germline cells is conserved in animal spermatogenesis. The testes of Drosophila melanogaster are the model of choice to identify processes responsible for mature gamete production. However, processes of differentiation and soma-germline interactions occurring in somatic cyst cells are currently understudied. Here we focused on the comparison of transcriptome expression patterns of early and mature somatic cyst cells to find out the developmental changes taking place in them. We employed a FACS-based approach for the isolation of early and mature somatic cyst cells from fly testes, subsequent preparation of RNA-Seq libraries, and analysis of gene differential expression in the sorted cells. We found increased expression of genes involved in cell cycle-related processes in early cyst cells, which is necessary for the proliferation and self-renewal of a crucial population of early cyst cells, cyst stem cells. Genes proposedly required for lamellipodium-like projection organization for proper cyst formation were also detected among the upregulated ones in early cyst cells. Gene Ontology and interactome analyses of upregulated genes in mature cyst cells revealed a striking over-representation of gene categories responsible for metabolic and catabolic cellular processes, as well as genes supporting the energetic state of the cells provided by oxidative phosphorylation that is carried out in mitochondria. Our comparative analyses of differentially expressed genes revealed major peculiarities in early and mature cyst cells and provide novel insight into their regulation, which is important for male fertility.

Nonsense Mutations in Eukaryotes

Nonsense mutations are a type of mutations which results in a premature termination codon occurrence. In general, these mutations have been considered to be among the most harmful ones which lead to premature protein translation termination and result in shortened nonfunctional polypeptide. However, there is evidence that not all nonsense mutations are harmful as well as some molecular mechanisms exist which allow to avoid pathogenic effects of these mutations. This review addresses relevant information on nonsense mutations in eukaryotic genomes, characteristics of these mutations, and different molecular mechanisms preventing or mitigating harmful effects thereof.