The Drosophila melanogaster Genetic Reference Panel

What can insects teach us about hearing loss?

Over the last three decades, insects have been utilized to provide a deep and fundamental understanding of many human diseases and disorders. Here, we present arguments for insects as models to understand general principles underlying hearing loss. Despite ∼600 million years since the last common ancestor of vertebrates and invertebrates, we share an overwhelming degree of genetic homology particularly with respect to auditory organ development and maintenance. Despite the anatomical differences between human and insect auditory organs, both share physiological principles of operation. We explain why these observations are expected and highlight areas in hearing loss research in which insects can provide insight. We start by briefly introducing the evolutionary journey of auditory organs, the reasons for using insect auditory organs for hearing loss research, and the tools and approaches available in insects. Then, the first half of the review focuses on auditory development and auditory disorders with a genetic cause. The second half analyses the physiological and genetic consequences of ageing and short- and long-term changes as a result of noise exposure. We finish with complex age and noise interactions in auditory systems. In this review, we present some of the evidence and arguments to support the use of insects to study mechanisms and potential treatments for hearing loss in humans. Obviously, insects cannot fully substitute for all aspects of human auditory function and loss of function, although there are many important questions that can be addressed in an animal model for which there are important ethical, practical and experimental advantages.

Rapid seasonal changes in phenotypes in a wild Drosophila population

Seasonal environmental change is one of the most rapid and striking environmental variables. Although relatively rapid adaptation to environmental changes over several years or several decades has been described in many taxa, rapid responses to seasonal environments are delicate, and therefore, the detection of the evolutionary responses requires sensitive methods. In this study, we examined seasonal changes in phenotypes related to thermal tolerance and morphological traits of Drosophila lutescens collected at the spring and autumn periods from a single location. We first demonstrated that flies in the two seasonal periods were almost genetically identical using double-digest restriction site-associated DNA sequencing and analysis. Using an experimental design to eliminate the effect of possible confounding factors that influence phenotypes (i.e., maternal effects and the environmental conditions in which each phenotype was analyzed), we showed that the heat tolerance of D. lutescens was significantly higher in the autumn population than in the spring population. Furthermore, cold tolerance was slightly higher in the spring population than in the autumn one. Although wing length and thorax length did not change significantly between seasons, the ratio of wing length to thorax length changed significantly between them. These results suggest that seasonal environmental heterogeneity induces rapid phenotypic changes within a year. Finally, we discuss the possibility of rapid evolutionary responses to seasonal changes.

Between the genotype and the phenotype lies the microbiome: symbiosis and the making of 'postgenomic' knowledge

Emphatic claims of a "microbiome revolution" aside, the study of the gut microbiota and its role in organismal development and evolution is a central feature of so-called postgenomics; namely, a conceptual and/or practical turn in contemporary life sciences, which departs from genetic determinism and reductionism to explore holism, emergentism and complexity in biological knowledge-production. This paper analyses the making of postgenomic knowledge about developmental symbiosis in Drosophila melanogaster by a specific group of microbiome scientists. Drawing from both practical philosophy of science and Science and Technology Studies, the paper documents epistemological questions of artefactuality and representativeness of model organisms as they emerge in the day-to-day labour producing and being produced by the "microbiome revolution." Specifically, the paper builds on all the written and editorial exchanges involved in the troubled publication of a research paper studying the symbiotic role of the microbiota in the flies' development. These written materials permit us to delimit the network of justifications, evidence, standards of knowledge-production, trust in the tools and research designs that make up the conditions of possibility of a postgenomic fact. More than reframing the organism as a radically novel multiplicity of reactive genomes, we conclude, doing postgenomic research on the microbiota and symbiosis means producing a story that deviates from the scripts embedded into the sociotechnical experimental systems of post-Human Genome Project life sciences.

Identification and Interpretation of A-to-I RNA Editing Events in Insect Transcriptomes

Adenosine-to-inosine (A-to-I) RNA editing is the most prevalent RNA modification in the nervous systems of metazoans. To study the biological significance of RNA editing, we first have to accurately identify these editing events from the transcriptome. The genome-wide identification of RNA editing sites remains a challenging task. In this review, we will first introduce the occurrence, regulation, and importance of A-to-I RNA editing and then describe the established bioinformatic procedures and difficulties in the accurate identification of these sit esespecially in small sized non-model insects. In brief, (1) to obtain an accurate profile of RNA editing sites, a transcriptome coupled with the DNA resequencing of a matched sample is favorable; (2) the single-cell sequencing technique is ready to be applied to RNA editing studies, but there are a few limitations to overcome; (3) during mapping and variant calling steps, various issues, like mapping and base quality, soft-clipping, and the positions of mismatches on reads, should be carefully considered; (4) Sanger sequencing of both RNA and the matched DNA is the best verification of RNA editing sites, but other auxiliary evidence, like the nonsynonymous-to-synonymous ratio or the linkage information, is also helpful for judging the reliability of editing sites. We have systematically reviewed the understanding of the biological significance of RNA editing and summarized the methodology for identifying such editing events. We also raised several promising aspects and challenges in this field. With insightful perspectives on both scientific and technical issues, our review will benefit the researchers in the broader RNA editing community.

Pleiotropic fitness effects across sexes and ages in the Drosophila genome and transcriptome

Selection varies between categories of individuals, with far-reaching ramifications: Sex-specific selection can impede or accelerate adaptation, and differences in selection between young and old individuals are ultimately responsible for senescence. Here, we measure early- and late-life fitness in adults of both sexes from the Drosophila genetic reference panel and perform quantitative genetic and transcriptomic analyses. Fitness was heritable, showed positive pleiotropy across sexes and age classes, and appeared to be influenced by very large numbers of loci with small effects plus a smaller number with moderate effects. Most loci affected male and female fitness in the same direction; relatively few candidate sexually antagonistic loci were found, though these were enriched on the X chromosome as predicted by theory. The expression level of many genes showed an opposite correlation with fitness in males and females, consistent with unresolved sexual conflict over transcription. The load of deleterious mutations correlated negatively with fitness across genotypes, and we found some evidence for the mutation accumulation (but not the antagonistic pleiotropy) theory of aging.

Historic Museum Samples Provide Evidence for a Recent Replacement of Wolbachia Types in European Drosophila melanogaster

Wolbachia is one of the most common bacterial endosymbionts, which is frequently found in numerous arthropods and nematode taxa. Wolbachia infections can have a strong influence on the evolutionary dynamics of their hosts since these bacteria are reproductive manipulators that affect the fitness and life history of their host species for their own benefit. Host-symbiont interactions with Wolbachia are perhaps best studied in the model organism Drosophila melanogaster, which is naturally infected with at least 5 different variants among which wMel and wMelCS are the most frequent ones. Comparisons of infection types between natural flies and long-term lab stocks have previously indicated that wMelCS represents the ancestral type, which was only very recently replaced by the nowadays dominant wMel in most natural populations. In this study, we took advantage of recently sequenced museum specimens of D. melanogaster that have been collected 90 to 200 yr ago in Northern Europe to test this hypothesis. Our comparison to contemporary Wolbachia samples provides compelling support for the replacement hypothesis. Our analyses show that sequencing data from historic museum specimens and their bycatch are an emerging and unprecedented resource to address fundamental questions about evolutionary dynamics in host-symbiont interactions. However, we also identified contamination with DNA from crickets that resulted in co-contamination with cricket-specific Wolbachia in several samples. These results underpin the need for rigorous quality assessments of museomic data sets to account for contamination as a source of error that may strongly influence biological interpretations if it remains undetected.

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.

Genetic independence between traits separated by metamorphosis is widespread but varies with biological function

Why is metamorphosis so pervasive? Does it facilitate the independent (micro)evolution of quantitative traits in distinct life stages, similarly to how it enables some limbs and organs to develop at specific life stages? We tested this hypothesis by measuring the expression of 6400 genes in 41 Drosophila melanogaster inbred lines at larval and adult stages. Only 30% of the genes showed significant genetic correlations between larval and adult expression. By contrast, 46% of the traits showed some level of genetic independence between stages. Gene ontology terms enrichment revealed that across stages correlated traits were often involved in proteins synthesis, insecticide resistance and innate immunity, while a vast number of genes expression traits associated with energy metabolism were independent between life stages. We compared our results to a similar case: genetic constraints between males and females in gonochoric species (i.e. sexual antagonism). We expected selection for the separation between males and females to be higher than between juvenile and adult functions, as gonochorism is a more common strategy in the animal kingdom than metamorphosis. Surprisingly, we found that inter-stage constraints were lower than inter-sexual genetic constraints. Overall, our results show that metamorphosis enables a large part of the transcriptome to evolve independently at different life stages.

Testing a candidate meiotic drive locus identified by pool sequencing

Meiotic drive biases the transmission of alleles in heterozygous individuals, such that Mendel's law of equal segregation is violated. Most examples of meiotic drive have been discovered over the past century based on causing sex ratio distortion or the biased transmission of easily scoreable genetic markers that were linked to drive alleles. More recently, several approaches have been developed that attempt to identify distortions of Mendelian segregation genome wide. Here, we test a candidate female meiotic drive locus in Drosophila melanogaster, identified previously as causing a ∼54:46 distortion ratio using sequencing of large pools of backcross progeny. We inserted fluorescent visible markers near the candidate locus and scored transmission in thousands of individual progeny. We observed a small but significant deviation from the Mendelian expectation; however, it was in the opposite direction to that predicted based on the original experiments. We discuss several possible causes of the discrepancy between the 2 approaches, noting that subtle viability effects are particularly challenging to disentangle from potential small-effect meiotic drive loci. We conclude that pool sequencing approaches remain a powerful method to identify candidate meiotic drive loci but that genotyping of individual progeny at early developmental stages may be required for robust confirmation.

Male-derived transcripts isolated from the mated female reproductive tract in Drosophila melanogaster

In species with internal fertilization, sperm, and seminal fluid are transferred from male to female during mating. While both sperm and seminal fluid contain various types of molecules, including RNA, the role of most of these molecules in the coordination of fertilization or in other possible functions is poorly understood. In Drosophila, exosomes from the accessory gland, which produces seminal fluid, are transferred to females, but their potential cargoes have not been described. Moreover, while the RNA composition of sperm has been described in several mammalian species, little work on this problem has occurred in Drosophila. Here we use single nucleotide polymorphism differences between males and females from a set of highly inbred lines of D. melanogaster, and transcriptome data from the female reproductive tract, sperm, testis, and accessory gland, to investigate the potential origin, male vs female, RNA molecules isolated from 3 female reproductive tract organs, the seminal receptacle and spermatheca, which store sperm, and the parovaria, which does not. We find that mated females carry male-derived transcripts from many genes, including those that are markers of the accessory gland and known seminal fluid proteins. Our observations also support the idea that intact sperm transcripts can be isolated from the female sperm storage organs.

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

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

Evolution of secondary cell number and position in the Drosophila accessory gland

In animals with internal fertilization, males transfer gametes and seminal fluid during copulation, both of which are required for successful reproduction. In Drosophila and other insects, seminal fluid is produced in the paired accessory gland (AG), the ejaculatory duct, and the ejaculatory bulb. The D. melanogaster AG has emerged as an important model system for this component of male reproductive biology. Seminal fluid proteins produced in the Drosophila AG are required for proper storage and use of sperm by the females, and are also critical for establishing and maintaining a suite of short- and long-term postcopulatory female physiological responses that promote reproductive success. The Drosophila AG is composed of two main cell types. The majority of AG cells, which are referred to as main cells, are responsible for production of many seminal fluid proteins. A minority of cells, about 4%, are referred to as secondary cells. These cells, which are restricted to the distal tip of the D. melanogaster AG, may play an especially important role in the maintenance of the long-term female post-mating response. Many studies of Drosophila AG evolution have suggested that the proteins produced in the gland evolve quickly, as does the transcriptome. Here, we investigate the evolution of secondary cell number and position in the AG in a collection of eight species spanning the entire history of the Drosophila genus. We document a heretofore underappreciated rapid evolutionary rate for both number and position of these specialized AG cells, raising several questions about the developmental, functional, and evolutionary significance of this variation.

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.

Evaluating power to detect recurrent selective sweeps under increasingly realistic evolutionary null models

The detection of selective sweeps from population genomic data often relies on the premise that the beneficial mutations in question have fixed very near the sampling time. As it has been previously shown that the power to detect a selective sweep is strongly dependent on the time since fixation as well as the strength of selection, it is naturally the case that strong, recent sweeps leave the strongest signatures. However, the biological reality is that beneficial mutations enter populations at a rate, one that partially determines the mean wait time between sweep events and hence their age distribution. An important question thus remains about the power to detect recurrent selective sweeps when they are modeled by a realistic mutation rate and as part of a realistic distribution of fitness effects, as opposed to a single, recent, isolated event on a purely neutral background as is more commonly modeled. Here we use forward-in-time simulations to study the performance of commonly used sweep statistics, within the context of more realistic evolutionary baseline models incorporating purifying and background selection, population size change, and mutation and recombination rate heterogeneity. Results demonstrate the important interplay of these processes, necessitating caution when interpreting selection scans; specifically, false-positive rates are in excess of true-positive across much of the evaluated parameter space, and selective sweeps are often undetectable unless the strength of selection is exceptionally strong.

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.

Uncovering the genetic architecture of parallel evolution

Identifying the genetic architecture underlying adaptive traits is exceptionally challenging in natural populations. This is because associations between traits not only mask the targets of selection but also create correlated patterns of genomic divergence that hinder our ability to isolate causal genetic effects. Here, we examine the repeated evolution of components of the auxin pathway that have contributed to the replicated loss of gravitropism (i.e. the ability of a plant to bend in response to gravity) in multiple populations of the Senecio lautus species complex in Australia. We use a powerful approach which combines parallel population genomics with association mapping in a Multiparent Advanced Generation Inter-Cross (MAGIC) population to break down genetic and trait correlations to reveal how adaptive traits evolve during replicated evolution. We sequenced auxin and shoot gravitropism-related gene regions in 80 individuals from six natural populations (three parallel divergence events) and 133 individuals from a MAGIC population derived from two of the recently diverged natural populations. We show that artificial tail selection on gravitropism in the MAGIC population recreates patterns of parallel divergence in the auxin pathway in the natural populations. We reveal a set of 55 auxin gene regions that have evolved repeatedly during the evolution of the species, of which 50 are directly associated with gravitropism divergence in the MAGIC population. Our work creates a strong link between patterns of genomic divergence and trait variation contributing to replicated evolution by natural selection, paving the way to understand the origin and maintenance of adaptations in natural populations.

Senolytic and senomorphic secondary metabolites as therapeutic agents in Drosophila melanogaster models of Parkinson's disease

Drosophila melanogaster is a valuable model organism for a wide range of biological exploration. The well-known advantages of D. melanogaster include its relatively simple biology, the ease with which it is genetically modified, the relatively low financial and time costs associated with their short gestation and life cycles, and the large number of offspring they produce per generation. D. melanogaster has facilitated the discovery of many significant insights into the pathology of Parkinson's disease (PD) and has served as an excellent preclinical model of PD-related therapeutic discovery. In this review, we provide an overview of the major D. melanogaster models of PD, each of which provide unique insights into PD-relevant pathology and therapeutic targets. These models are discussed in the context of their past, current, and future potential use for studying the utility of secondary metabolites as therapeutic agents in PD. Over the last decade, senolytics have garnered an exponential interest in their ability to mitigate a broad spectrum of diseases, including PD. Therefore, an emphasis is placed on the senolytic and senomorphic properties of secondary metabolites. It is expected that D. melanogaster will continue to be critical in the effort to understand and improve treatment of PD, including their involvement in translational studies focused on secondary metabolites.

A novel panel of Drosophila TAFAZZIN mutants in distinct genetic backgrounds as a resource for therapeutic testing

Barth Syndrome is a rare, X-linked disorder caused by mutation of the gene TAFAZZIN (TAZ). The corresponding Tafazzin protein is involved in the remodeling of cardiolipin, a phospholipid with critical roles in mitochondrial function. While recent clinical trials have been promising, there is still no cure for Barth Syndrome. Because TAZ is highly conserved, multiple animal and cell culture models exist for pre-clinical testing of therapeutics. However, since the same mutation in different patients can lead to different symptoms and responses to treatment, isogenized experimental models can't fully account for human disease conditions. On the other hand, isogenized animal models allow for sufficient numbers to thoroughly establish efficacy for a given genetic background. Therefore, a combined method for testing treatments in a panel of isogenized cohorts that are genetically distinct from each other would be transformative for testing emerging pre-clinical therapies. To aid in this effort, we've created a novel panel of 10 Drosophila lines, each with the same TAZ mutation in highly diverse genetic backgrounds, to serve as a helpful resource to represent natural variation in background genetics in pre-clinical studies. As a proof of principle, we test our panel here using nicotinamide riboside (NR), a treatment with established therapeutic value, to evaluate how robust this therapy is across the 10 genetic backgrounds in this novel reference panel. We find substantial variation in the response to NR across backgrounds. We expect this resource will be valuable in pre-clinical testing of emerging therapies for Barth Syndrome.

A predominant role of genotypic variation in both expression of sperm competition genes and paternity success in Drosophila melanogaster

Sperm competition is a crucial aspect of male reproductive success in many species, including Drosophila melanogaster, and seminal fluid proteins (Sfps) can influence sperm competitiveness. However, the combined effect of environmental and genotypic variation on sperm competition gene expression remains poorly understood. Here, we used Drosophila Genetic Reference Panel (DGRP) inbred lines and manipulated developmental population density (i.e. larval density) to test the effects of genotype, environment and genotype-by-environment interactions (GEI) on the expression of the known sperm competition genes Sex Peptide, Acp36DE and CG9997. High larval density resulted in reduced adult body size, but expression of sperm competition genes remained unaffected. Furthermore, we found no significant GEI but genotypic effects in the expression of SP and Acp36DE. Our results also revealed GEI for relative competitive paternity success (second male paternity; P2), with genes' expression positively correlated with P2. Given the effect of genotype on the expression of genes, we conducted a genome-wide association study (GWAS) and identified polymorphisms in putative cis-regulatory elements as predominant factors regulating the expression of SP and Acp36DE. The association of genotypic variation with sperm competition outcomes, and the resilience of sperm competition genes' expression against environmental challenges, demonstrates the importance of genome variation background in reproductive fitness.

Nacα protects the larval fat body from cell death by maintaining cellular proteostasis in Drosophila

Protein homeostasis (proteostasis) is crucial for the maintenance of cellular homeostasis. Impairment of proteostasis activates proteotoxic and unfolded protein response pathways to resolve cellular stress or induce apoptosis in damaged cells. However, the responses of individual tissues to proteotoxic stress and evoking cell death program have not been extensively explored in vivo. Here, we show that a reduction in Nascent polypeptide-associated complex protein alpha subunit (Nacα) specifically and progressively induces cell death in Drosophila fat body cells. Nacα mutants disrupt both ER integrity and the proteasomal degradation system, resulting in caspase activation through JNK and p53. Although forced activation of the JNK and p53 pathways was insufficient to induce cell death in the fat body, the reduction of Nacα sensitized fat body cells to intrinsic and environmental stresses. Reducing overall protein synthesis by mTor inhibition or Minute mutants alleviated the cell death phenotype in Nacα mutant fat body cells. Our work revealed that Nacα is crucial for protecting the fat body from cell death by maintaining cellular proteostasis, thus demonstrating the coexistence of a unique vulnerability and cell death resistance in the fat body.

Genetic Variation in Sexual Size Dimorphism Is Associated with Variation in Sex-Specific Plasticity in Drosophila

AbstractThe difference in body size between females and males, or sexual size dimorphism (SSD), is ubiquitous, yet we have a poor understanding of the developmental genetic mechanisms that generate it and how these mechanisms may vary within and among species. Such an understanding of the genetic architecture of SSD is important if we are to evaluate alternative models of SSD evolution, but the genetic architecture is difficult to describe because SSD is a characteristic of populations, not individuals. Here, we overcome this challenge by using isogenic lineages of Drosophila to measure SSD for 196 genotypes. We demonstrate extensive genetic variation for SSD, primarily driven by higher levels of genetic variation for body size among females than among males. While we observe a general increase in SSD with sex-averaged body size (pooling for sex) among lineages, most of the variation in SSD is independent of sex-averaged body size and shows a strong genetic correlation with sex-specific plasticity, such that increased female-biased SSD is associated with increased body size plasticity in females. Our data are consistent with the condition dependence hypothesis of sexual dimorphism and suggest that SSD in Drosophila is a consequence of selection on the developmental genetic mechanisms that regulate the plasticity of body size.

Using Drosophila to identify naturally occurring genetic modifiers of amyloid beta 42- and tau-induced toxicity

Alzheimer's disease is characterized by 2 pathological proteins, amyloid beta 42 and tau. The majority of Alzheimer's disease cases in the population are sporadic and late-onset Alzheimer's disease, which exhibits high levels of heritability. While several genetic risk factors for late-onset Alzheimer's disease have been identified and replicated in independent studies, including the ApoE ε4 allele, the great majority of the heritability of late-onset Alzheimer's disease remains unexplained, likely due to the aggregate effects of a very large number of genes with small effect size, as well as to biases in sample collection and statistical approaches. Here, we present an unbiased forward genetic screen in Drosophila looking for naturally occurring modifiers of amyloid beta 42- and tau-induced ommatidial degeneration. Our results identify 14 significant SNPs, which map to 12 potential genes in 8 unique genomic regions. Our hits that are significant after genome-wide correction identify genes involved in neuronal development, signal transduction, and organismal development. Looking more broadly at suggestive hits (P < 10-5), we see significant enrichment in genes associated with neurogenesis, development, and growth as well as significant enrichment in genes whose orthologs have been identified as significantly or suggestively associated with Alzheimer's disease in human GWAS studies. These latter genes include ones whose orthologs are in close proximity to regions in the human genome that are associated with Alzheimer's disease, but where a causal gene has not been identified. Together, our results illustrate the potential for complementary and convergent evidence provided through multitrait GWAS in Drosophila to supplement and inform human studies, helping to identify the remaining heritability and novel modifiers of complex diseases.

Basidiomycota species in Drosophila gut are associated with host fat metabolism

The importance of bacterial microbiota on host metabolism and obesity risk is well documented. However, the role of fungal microbiota on host storage metabolite pools is largely unexplored. We aimed to investigate the role of microbiota on D. melanogaster fat metabolism, and examine interrelatedness between fungal and bacterial microbiota, and major metabolic pools. Fungal and bacterial microbiota profiles, fat, glycogen, and trehalose metabolic pools are measured in a context of genetic variation represented by whole genome sequenced inbred Drosophila Genetic Reference Panel (DGRP) samples. Increasing Basidiomycota, Acetobacter persici, Acetobacter pomorum, and Lactobacillus brevis levels correlated with decreasing triglyceride levels. Host genes and biological pathways, identified via genome-wide scans, associated with Basidiomycota and triglyceride levels were different suggesting the effect of Basidiomycota on fat metabolism is independent of host biological pathways that control fungal microbiota or host fat metabolism. Although triglyceride, glycogen and trehalose levels were highly correlated, microorganisms' effect on triglyceride pool were independent of glycogen and trehalose levels. Multivariate analyses suggested positive interactions between Basidiomycota, A. persici, and L. brevis that collectively correlated negatively with fat and glycogen pools. In conclusion, fungal microbiota can be a major player in host fat metabolism. Interactions between fungal and bacterial microbiota may exert substantial control over host storage metabolite pools and influence obesity risk.

Natural selection has driven the recurrent loss of an immunity gene that protects Drosophila against a major natural parasite

Polymorphisms in immunity genes can have large effects on susceptibility to infection. To understand the origins of this variation, we have investigated the genetic basis of resistance to the parasitoid wasp Leptopilina boulardi in Drosophila melanogaster. We found that increased expression of the gene lectin-24A after infection by parasitic wasps was associated with a faster cellular immune response and greatly increased rates of killing the parasite. lectin-24A encodes a protein that is strongly up-regulated in the fat body after infection and localizes to the surface of the parasite egg. In certain susceptible lines, a deletion upstream of the lectin-24A has largely abolished expression. Other mutations predicted to abolish the function of this gene have arisen recurrently in this gene, with multiple loss-of-expression alleles and premature stop codons segregating in natural populations. The frequency of these alleles varies greatly geographically, and in some southern African populations, natural selection has driven them near to fixation. We conclude that natural selection has favored the repeated loss of an important component of the immune system, suggesting that in some populations, a pleiotropic cost to lectin-24A expression outweighs the benefits of resistance.

GenomeMUSter mouse genetic variation service enables multi-trait, multi-population data integration and analyses

Hundreds of inbred laboratory mouse strains and intercross populations have been used to functionalize genetic variants that contribute to disease. Thousands of disease relevant traits have been characterized in mice and made publicly available. New strains and populations including the Collaborative Cross, expanded BXD and inbred wild-derived strains add to set of complex disease mouse models, genetic mapping resources and sensitized backgrounds against which to evaluate engineered mutations. The genome sequences of many inbred strains, along with dense genotypes from others could allow integrated analysis of trait - variant associations across populations, but these analyses are not feasible due to the sparsity of genotypes available. Moreover, the data are not readily interoperable with other resources. To address these limitations, we created a uniformly dense data resource by harmonizing multiple variant datasets. Missing genotypes were imputed using the Viterbi algorithm with a data-driven technique that incorporates local phylogenetic information, an approach that is extensible to other model organism species. The result is a web- and programmatically-accessible data service called GenomeMUSter ( https://muster.jax.org ), comprising allelic data covering 657 strains at 106.8M segregating sites. Interoperation with phenotype databases, analytic tools and other resources enable a wealth of applications including multi-trait, multi-population meta-analysis. We demonstrate this in a cross-species comparison of the meta-analysis of Type 2 Diabetes and of substance use disorders, resulting in the more specific characterization of the role of human variant effects in light of mouse phenotype data. Other applications include refinement of mapped loci and prioritization of strain backgrounds for disease modeling to further unlock extant mouse diversity for genetic and genomic studies in health and disease.

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

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

Evolution of chemosensory and detoxification gene families across herbivorous Drosophilidae

Herbivorous insects are exceptionally diverse, accounting for a quarter of all known eukaryotic species, but the genomic basis of adaptations that enabled this dietary transition remains poorly understood. Many studies have suggested that expansions and contractions of chemosensory and detoxification gene families-genes directly mediating interactions with plant chemical defenses-underlie successful plant colonization. However, this hypothesis has been challenging to test because the origins of herbivory in many insect lineages are ancient (>150 million years ago (mya)), obscuring genomic evolutionary patterns. Here, we characterized chemosensory and detoxification gene family evolution across Scaptomyza, a genus nested within Drosophila that includes a recently derived (<15 mya) herbivore lineage of mustard (Brassicales) specialists and carnation (Caryophyllaceae) specialists, and several nonherbivorous species. Comparative genomic analyses revealed that herbivorous Scaptomyza has among the smallest chemosensory and detoxification gene repertoires across 12 drosophilid species surveyed. Rates of gene turnover averaged across the herbivore clade were significantly higher than background rates in over half of the surveyed gene families. However, gene turnover was more limited along the ancestral herbivore branch, with only gustatory receptors and odorant-binding proteins experiencing strong losses. The genes most significantly impacted by gene loss, duplication, or changes in selective constraint were those involved in detecting compounds associated with feeding on living plants (bitter or electrophilic phytotoxins) or their ancestral diet (fermenting plant volatiles). These results provide insight into the molecular and evolutionary mechanisms of plant-feeding adaptations and highlight gene candidates that have also been linked to other dietary transitions in Drosophila.

Sex-specific estimation of cis and trans regulation of gene expression in heads and gonads of Drosophila melanogaster

The regulatory architecture of gene expression is known to differ substantially between sexes in Drosophila, but most studies performed so far used whole-body data and only single crosses, which may have limited their scope to detect patterns that are robust across tissues and biological replicates. Here, we use allele-specific gene expression of parental and reciprocal hybrid crosses between 6 Drosophila melanogaster inbred lines to quantify cis- and trans-regulatory variation in heads and gonads of both sexes separately across 3 replicate crosses. Our results suggest that female and male heads, as well as ovaries, have a similar regulatory architecture. On the other hand, testes display more and substantially different cis-regulatory effects, suggesting that sex differences in the regulatory architecture that have been previously observed may largely derive from testis-specific effects. We also examine the difference in cis-regulatory variation of genes across different levels of sex bias in gonads and heads. Consistent with the idea that intersex correlations constrain expression and can lead to sexual antagonism, we find more cis variation in unbiased and moderately biased genes in heads. In ovaries, reduced cis variation is observed for male-biased genes, suggesting that cis variants acting on these genes in males do not lead to changes in ovary expression. Finally, we examine the dominance patterns of gene expression and find that sex- and tissue-specific patterns of inheritance as well as trans-regulatory variation are highly variable across biological crosses, although these were performed in highly controlled experimental conditions. This highlights the importance of using various genetic backgrounds to infer generalizable patterns.

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

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

Nonlinear expression patterns and multiple shifts in gene network interactions underlie robust phenotypic change in Drosophila melanogaster selected for night sleep duration

All but the simplest phenotypes are believed to result from interactions between two or more genes forming complex networks of gene regulation. Sleep is a complex trait known to depend on the system of feedback loops of the circadian clock, and on many other genes; however, the main components regulating the phenotype and how they interact remain an unsolved puzzle. Genomic and transcriptomic data may well provide part of the answer, but a full account requires a suitable quantitative framework. Here we conducted an artificial selection experiment for sleep duration with RNA-seq data acquired each generation. The phenotypic results are robust across replicates and previous experiments, and the transcription data provides a high-resolution, time-course data set for the evolution of sleep-related gene expression. In addition to a Hierarchical Generalized Linear Model analysis of differential expression that accounts for experimental replicates we develop a flexible Gaussian Process model that estimates interactions between genes. 145 gene pairs are found to have interactions that are different from controls. Our method appears to be not only more specific than standard correlation metrics but also more sensitive, finding correlations not significant by other methods. Statistical predictions were compared to experimental data from public databases on gene interactions. Mutations of candidate genes implicated by our results affected night sleep, and gene expression profiles largely met predicted gene-gene interactions.

Perspectives on the Drosophila melanogaster Model for Advances in Toxicological Science

The use of Drosophila melanogaster for studies of toxicology has grown considerably in the last decade. The Drosophila model has long been appreciated as a versatile and powerful model for developmental biology and genetics because of its ease of handling, short life cycle, low cost of maintenance, molecular genetic accessibility, and availability of a wide range of publicly available strains and data resources. These features, together with recent unique developments in genomics and metabolomics, make the fly model especially relevant and timely for the development of new approach methodologies and movements toward precision toxicology. Here, we offer a perspective on how flies can be leveraged to identify risk factors relevant to environmental exposures and human health. First, we review and discuss fundamental toxicologic principles for experimental design with Drosophila. Next, we describe quantitative and systems genetics approaches to resolve the genetic architecture and candidate pathways controlling susceptibility to toxicants. Finally, we summarize the current state and future promise of the emerging field of Drosophila metabolomics for elaborating toxic mechanisms. © 2023 The Authors. Current Protocols published by Wiley Periodicals LLC.

Deep data mining reveals variable abundance and distribution of microbial reproductive manipulators within and among diverse host species

Bacterial symbionts that manipulate the reproduction of their hosts are important factors in invertebrate ecology and evolution, and are being leveraged for host biological control. Infection prevalence restricts which biological control strategies are possible and is thought to be strongly influenced by the density of symbiont infection within hosts, termed titer. Current methods to estimate infection prevalence and symbiont titers are low-throughput, biased towards sampling infected species, and rarely measure titer. Here we develop a data mining approach to estimate symbiont infection frequencies within host species and titers within host tissues. We applied this approach to screen ~32,000 publicly available sequence samples from the most common symbiont host taxa, discovering 2,083 arthropod and 119 nematode infected samples. From these data, we estimated that Wolbachia infects approximately 44% of all arthropod and 34% of all nematode species, while other reproductive manipulators only infect 1-8% of arthropod and nematode species. Although relative titers within hosts were highly variable within and between arthropod species, a combination of arthropod host species and Wolbachia strain explained approximately 36% of variation in Wolbachia titer across the dataset. To explore potential mechanisms for host control of symbiont titer, we leveraged population genomic data from the model system Drosophila melanogaster. In this host, we found a number of SNPs associated with titer in candidate genes potentially relevant to host interactions with Wolbachia. Our study demonstrates that data mining is a powerful tool to detect bacterial infections and quantify infection intensities, thus opening an array of previously inaccessible data for further analysis in host-symbiont evolution.

Evolutionary genetics of the mitochondrial genome: insights from Drosophila

Mitochondria are key to energy conversion in virtually all eukaryotes. Intriguingly, despite billions of years of evolution inside the eukaryote, mitochondria have retained their own small set of genes involved in the regulation of oxidative phosphorylation (OXPHOS) and protein translation. Although there was a long-standing assumption that the genetic variation found within the mitochondria would be selectively neutral, research over the past 3 decades has challenged this assumption. This research has provided novel insight into the genetic and evolutionary forces that shape mitochondrial evolution and broader implications for evolutionary ecological processes. Many of the seminal studies in this field, from the inception of the research field to current studies, have been conducted using Drosophila flies, thus establishing the species as a model system for studies in mitochondrial evolutionary biology. In this review, we comprehensively review these studies, from those focusing on genetic processes shaping evolution within the mitochondrial genome, to those examining the evolutionary implications of interactions between genes spanning mitochondrial and nuclear genomes, and to those investigating the dynamics of mitochondrial heteroplasmy. We synthesize the contribution of these studies to shaping our understanding of the evolutionary and ecological implications of mitochondrial genetic variation.

Hybrid breakdown in male reproduction between recently diverged Drosophila melanogaster populations has a complex and variable genetic architecture

Secondary contact between formerly isolated populations may result in hybrid breakdown, in which untested allelic combinations in hybrids are maladaptive and limit genetic exchange. Studying early-stage reproductive isolation may yield key insights into the genetic architectures and evolutionary forces underlying the first steps toward speciation. Here, we leverage the recent worldwide expansion of Drosophila melanogaster to test for hybrid breakdown between populations that diverged within the last 13,000 years. We found clear evidence for hybrid breakdown in male reproduction, but not female reproduction or viability, supporting the prediction that hybrid breakdown affects the heterogametic sex first. The frequency of non-reproducing F2 males varied among different crosses involving the same southern African and European populations, as did the qualitative effect of cross direction, implying a genetically variable basis of hybrid breakdown and a role for uniparentally inherited factors. The levels of breakdown observed in F2 males were not recapitulated in backcrossed individuals, consistent with the existence of incompatibilities with at least three partners. Thus, some of the very first steps toward reproductive isolation could involve incompatibilities with complex and variable genetic architectures. Collectively, our findings emphasize this system's potential for future studies on the genetic and organismal basis of early-stage reproductive isolation.

A single-cell atlas of the sexually dimorphic Drosophila foreleg and its sensory organs during development

To respond to the world around them, animals rely on the input of a network of sensory organs distributed throughout the body. Distinct classes of sensory organs are specialized for the detection of specific stimuli such as strain, pressure, or taste. The features that underlie this specialization relate both to the neurons that innervate sensory organs and the accessory cells they comprise. To understand the genetic basis of this diversity of cell types, both within and between sensory organs, we performed single-cell RNA sequencing on the first tarsal segment of the male Drosophila melanogaster foreleg during pupal development. This tissue displays a wide variety of functionally and structurally distinct sensory organs, including campaniform sensilla, mechanosensory bristles, and chemosensory taste bristles, as well as the sex comb, a recently evolved male-specific structure. In this study, we characterize the cellular landscape in which the sensory organs reside, identify a novel cell type that contributes to the construction of the neural lamella, and resolve the transcriptomic differences among support cells within and between sensory organs. We identify the genes that distinguish between mechanosensory and chemosensory neurons, resolve a combinatorial transcription factor code that defines 4 distinct classes of gustatory neurons and several types of mechanosensory neurons, and match the expression of sensory receptor genes to specific neuron classes. Collectively, our work identifies core genetic features of a variety of sensory organs and provides a rich, annotated resource for studying their development and function.

Social modulation of oogenesis and egg laying in Drosophila melanogaster

Being part of a group facilitates cooperation between group members but also creates competition for resources. This is a conundrum for gravid females, whose future offspring benefit from being in a group only if there are enough resources relative to group size. Females may therefore be expected to modulate reproductive output depending on social context. In the fruit fly Drosophila melanogaster, females actively attract conspecifics to lay eggs on the same resources, generating groups in which individuals may cooperate or compete. The genetic tractability of this species allows dissecting the mechanisms underlying physiological adaptation to social context. Here, we show that females produce eggs increasingly faster as group size increases. By laying eggs faster when grouped than when isolated, females reduce competition between offspring and increase offspring survival. In addition, grouped females lay eggs during the day, while isolated females lay them at night. We show that responses to the presence of others requires visual input and that flies from any sex, mating status, or species can trigger these responses. The mechanisms of this modulation of egg laying by group is connected to a lifting of the inhibition of light on oogenesis and egg laying, possibly mediated in part by an increase in juvenile hormone activity. Because modulation of reproduction by social context is a hallmark of animals with higher levels of sociality, our findings in a species considered solitary question the validity of this nomenclature and suggest a widespread and profound influence of social context on reproduction.

Evaluating power to detect recurrent selective sweeps under increasingly realistic evolutionary null models

The detection of selective sweeps from population genomic data often relies on the premise that the beneficial mutations in question have fixed very near the sampling time. As it has been previously shown that the power to detect a selective sweep is strongly dependent on the time since fixation as well as the strength of selection, it is naturally the case that strong, recent sweeps leave the strongest signatures. However, the biological reality is that beneficial mutations enter populations at a rate, one that partially determines the mean wait time between sweep events and hence their age distribution. An important question thus remains about the power to detect recurrent selective sweeps when they are modelled by a realistic mutation rate and as part of a realistic distribution of fitness effects (DFE), as opposed to a single, recent, isolated event on a purely neutral background as is more commonly modelled. Here we use forward-in-time simulations to study the performance of commonly used sweep statistics, within the context of more realistic evolutionary baseline models incorporating purifying and background selection, population size change, and mutation and recombination rate heterogeneity. Results demonstrate the important interplay of these processes, necessitating caution when interpreting selection scans; specifically, false positive rates are in excess of true positive across much of the evaluated parameter space, and selective sweeps are often undetectable unless the strength of selection is exceptionally strong.

Teaser Text

Outlier-based genomic scans have proven a popular approach for identifying loci that have potentially experienced recent positive selection. However, it has previously been shown that an evolutionarily appropriate baseline model that incorporates non-equilibrium population histories, purifying and background selection, and variation in mutation and recombination rates is necessary to reduce often extreme false positive rates when performing genomic scans. Here we evaluate the power to detect recurrent selective sweeps using common SFS-based and haplotype-based methods under these increasingly realistic models. We find that while these appropriate evolutionary baselines are essential to reduce false positive rates, the power to accurately detect recurrent selective sweeps is generally low across much of the biologically relevant parameter space.

Unveiling Subtle Geographical Clines: Phenotypic Effects and Dynamics of Circadian Clock Gene Polymorphisms

Our understanding of the gene regulatory network that constitutes the circadian clock has greatly increased in recent decades, notably due to the use of Drosophila as a model system. In contrast, the analysis of natural genetic variation that enables the robust function of the clock under a broad range of environments has developed more slowly. In the current study, we analyzed comprehensive genome sequencing data from wild European populations of Drosophila, which were densely sampled through time and space. We identified hundreds of single nucleotide polymorphisms (SNPs) in nine genes associated with the clock, 276 of which exhibited a latitudinal cline in their allele frequencies. While the effect sizes of these clinal patterns were small, indicating subtle adaptations driven by natural selection, they provided important insights into the genetic dynamics of circadian rhythms in natural populations. We selected nine SNPs in different genes and assessed their impact on circadian and seasonal phenotypes by reconstructing outbred populations fixed for either of the SNP alleles, from inbred DGRP strains. The circadian free-running period of the locomotor activity rhythm was affected by an SNP in doubletime (dbt) and eyes absent (Eya). The SNPs in Clock (Clk), Shaggy (Sgg), period (per), and timeless (tim) affected the acrophase. The alleles of the SNP in Eya conferred different levels of diapause and the chill coma recovery response.

Belly roll, a GPI-anchored Ly6 protein, regulates Drosophila melanogaster escape behaviors by modulating the excitability of nociceptive peptidergic interneurons

Appropriate modulation of escape behaviors in response to potentially damaging stimuli is essential for survival. Although nociceptive circuitry has been studied, it is poorly understood how genetic contexts affect relevant escape responses. Using an unbiased genome-wide association analysis, we identified an Ly6/α-neurotoxin family protein, Belly roll (Bero), which negatively regulates Drosophila nociceptive escape behavior. We show that Bero is expressed in abdominal leucokinin-producing neurons (ABLK neurons) and bero knockdown in ABLK neurons resulted in enhanced escape behavior. Furthermore, we demonstrated that ABLK neurons responded to activation of nociceptors and initiated the behavior. Notably, bero knockdown reduced persistent neuronal activity and increased evoked nociceptive responses in ABLK neurons. Our findings reveal that Bero modulates an escape response by regulating distinct neuronal activities in ABLK neurons.

Population genomics reveals mechanisms and dynamics of de novo expressed open reading frame emergence in Drosophila melanogaster

Novel genes are essential for evolutionary innovations and differ substantially even between closely related species. Recently, multiple studies across many taxa showed that some novel genes arise de novo, that is, from previously noncoding DNA. To characterize the underlying mutations that allowed de novo gene emergence and their order of occurrence, homologous regions must be detected within noncoding sequences in closely related sister genomes. So far, most studies do not detect noncoding homologs of de novo genes because of incomplete assemblies and annotations, and long evolutionary distances separating genomes. Here, we overcome these issues by searching for de novo expressed open reading frames (neORFs), the not-yet fixed precursors of de novo genes that emerged within a single species. We sequenced and assembled genomes with long-read technology and the corresponding transcriptomes from inbred lines of Drosophila melanogaster, derived from seven geographically diverse populations. We found line-specific neORFs in abundance but few neORFs shared by lines, suggesting a rapid turnover. Gain and loss of transcription is more frequent than the creation of ORFs, for example, by forming new start and stop codons. Consequently, the gain of ORFs becomes rate limiting and is frequently the initial step in neORFs emergence. Furthermore, transposable elements (TEs) are major drivers for intragenomic duplications of neORFs, yet TE insertions are less important for the emergence of neORFs. However, highly mutable genomic regions around TEs provide new features that enable gene birth. In conclusion, neORFs have a high birth-death rate, are rapidly purged, but surviving neORFs spread neutrally through populations and within genomes.

Investigating the outcomes of virus coinfection within and across host species

Interactions between coinfecting pathogens have the potential to alter the course of infection and can act as a source of phenotypic variation in susceptibility between hosts. This phenotypic variation may influence the evolution of host-pathogen interactions within host species and interfere with patterns in the outcomes of infection across host species. Here, we examine experimental coinfections of two Cripaviruses-Cricket Paralysis Virus (CrPV), and Drosophila C Virus (DCV)-across a panel of 25 Drosophila melanogaster inbred lines and 47 Drosophilidae host species. We find that interactions between these viruses alter viral loads across D. melanogaster genotypes, with a ~3 fold increase in the viral load of DCV and a ~2.5 fold decrease in CrPV in coinfection compared to single infection, but we find little evidence of a host genetic basis for these effects. Across host species, we find no evidence of systematic changes in susceptibility during coinfection, with no interaction between DCV and CrPV detected in the majority of host species. These results suggest that phenotypic variation in coinfection interactions within host species can occur independently of natural host genetic variation in susceptibility, and that patterns of susceptibility across host species to single infections can be robust to the added complexity of coinfection.

Indel driven rapid evolution of core nuclear pore protein gene promoters

Nuclear pore proteins (Nups) prominently are among the few genes linked to speciation from hybrid incompatibility in Drosophila. These studies have focused on coding sequence evolution of Nup96 and Nup160 and shown evidence of positive selection driving nucleoporin evolution. Intriguingly, channel Nup54 functionality is required for neuronal wiring underlying the female post-mating response induced by male-derived sex-peptide. A region of rapid evolution in the core promoter of Nup54 suggests a critical role for general transcriptional regulatory elements at the onset of speciation, but whether this is a general feature of Nup genes has not been determined. Consistent with findings for Nup54, additional channel Nup58 and Nup62 promoters also rapidly accumulate insertions/deletions (indels). Comprehensive examination of Nup upstream regions reveals that core Nup complex gene promoters accumulate indels rapidly. Since changes in promoters can drive changes in expression, these results indicate an evolutionary mechanism driven by indel accumulation in core Nup promoters. Compensation of such gene expression changes could lead to altered neuronal wiring, rapid fixation of traits caused by promoter changes and subsequently the rise of new species. Hence, the nuclear pore complex may act as a nexus for species-specific changes via nucleo-cytoplasmic transport regulated gene expression.

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

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

Ovarian transcriptional response to Wolbachia infection in D. melanogaster in the context of between-genotype variation in gene expression

Wolbachia is a maternally transmitted endosymbiotic bacteria that infects a wide variety of arthropod and nematode hosts. The effects of Wolbachia on host biology are far-reaching and include changes in host gene expression. However, previous work on the host transcriptional response has generally been investigated in the context of a single host genotype. Thus, the relative effect of Wolbachia infection versus vs. host genotype on gene expression is unknown. Here, we explicitly test the relative roles of Wolbachia infection and host genotype on host gene expression by comparing the ovarian transcriptomes of 4 strains of Drosophila melanogaster (D. melanogaster) infected and uninfected with Wolbachia. Our data suggest that infection explains a small amount of transcriptional variation, particularly in comparison to variation in gene expression among strains. However, infection specifically affects genes related to cell cycle, translation, and metabolism. We also find enrichment of cell division and recombination processes among genes with infection-associated differential expression. Broadly, the transcriptomic changes identified in this study provide novel understanding of the relative magnitude of the effect of Wolbachia infection on gene expression in the context of host genetic variation and also point to genes that are consistently differentially expressed in response to infection among multiple genotypes.

Sex-Biased Expression Is Associated With Chromatin State in Drosophila melanogaster and Drosophila simulans

In Drosophila melanogaster and D. simulans head tissue, 60% of orthologous genes show evidence of sex-biased expression in at least one species. Of these, ∼39% (2,192) are conserved in direction. We hypothesize enrichment of open chromatin in the sex where we see expression bias and closed chromatin in the opposite sex. Male-biased orthologs are significantly enriched for H3K4me3 marks in males of both species (∼89% of male-biased orthologs vs. ∼76% of unbiased orthologs). Similarly, female-biased orthologs are significantly enriched for H3K4me3 marks in females of both species (∼90% of female-biased orthologs vs. ∼73% of unbiased orthologs). The sex-bias ratio in female-biased orthologs was similar in magnitude between the two species, regardless of the closed chromatin (H3K27me2me3) marks in males. However, in male-biased orthologs, the presence of H3K27me2me3 in both species significantly reduced the correlation between D. melanogaster sex-bias ratio and the D. simulans sex-bias ratio. Male-biased orthologs are enriched for evidence of positive selection in the D. melanogaster group. There are more male-biased genes than female-biased genes in both species. For orthologs with gains/losses of sex-bias between the two species, there is an excess of male-bias compared to female-bias, but there is no consistent pattern in the relationship between H3K4me3 or H3K27me2me3 chromatin marks and expression. These data suggest chromatin state is a component of the maintenance of sex-biased expression and divergence of sex-bias between species is reflected in the complexity of the chromatin status.

Genetic variation of morphological scaling in Drosophila melanogaster

Morphological scaling relationships between the sizes of individual traits and the body captures the characteristic shape of a species, and their evolution is the primary mechanism of morphological diversification. However, we have almost no knowledge of the genetic variation of scaling, which is critical if we are to understand how scaling evolves. Here we explore the genetics of population scaling relationships (scaling relationships fit to multiple genetically-distinct individuals in a population) by describing the distribution of individual scaling relationships (genotype-specific scaling relationships that are unseen or cryptic). These individual scaling relationships harbor the genetic variation in the developmental mechanisms that regulate trait growth relative to body growth, and theoretical studies suggest that their distribution dictates how the population scaling relationship will respond to selection. Using variation in nutrition to generate size variation within 197 isogenic lineages of Drosophila melanogaster, we reveal extensive variation in the slopes of the wing-body and leg-body individual scaling relationships among genotypes. This variation reflects variation in the nutritionally-induced size plasticity of the wing, leg, and body. Surprisingly, we find that variation in the slope of individual scaling relationships primarily results from variation in nutritionally-induced plasticity of body size, not leg or wing size. These data allow us to predict how different selection regimes affect scaling in Drosophila, and is the first step in identifying the genetic targets of such selection. More generally, our approach provides a framework for understanding the genetic variation of scaling, an important prerequisite to explaining how selection changes scaling and morphology.

Genetic variation in chromatin state across multiple tissues in Drosophila melanogaster

We use ATAC-seq to examine chromatin accessibility for four different tissues in Drosophila melanogaster: adult female brain, ovaries, and both wing and eye-antennal imaginal discs from males. Each tissue is assayed in eight different inbred strain genetic backgrounds, seven associated with a reference quality genome assembly. We develop a method for the quantile normalization of ATAC-seq fragments and test for differences in coverage among genotypes, tissues, and their interaction at 44099 peaks throughout the euchromatic genome. For the strains with reference quality genome assemblies, we correct ATAC-seq profiles for read mis-mapping due to nearby polymorphic structural variants (SVs). Comparing coverage among genotypes without accounting for SVs results in a highly elevated rate (55%) of identifying false positive differences in chromatin state between genotypes. After SV correction, we identify 1050, 30383, and 4508 regions whose peak heights are polymorphic among genotypes, among tissues, or exhibit genotype-by-tissue interactions, respectively. Finally, we identify 3988 candidate causative variants that explain at least 80% of the variance in chromatin state at nearby ATAC-seq peaks.

Natural genetic variation in a dopamine receptor is associated with variation in female fertility in Drosophila melanogaster

Fertility is a major component of fitness but its genetic architecture remains poorly understood. Using a full diallel cross of 50 Drosophila Genetic Reference Panel inbred lines with whole genome sequences, we found substantial genetic variation in fertility largely attributable to females. We mapped genes associated with variation in female fertility by genome-wide association analysis of common variants in the fly genome. Validation of candidate genes by RNAi knockdown confirmed the role of the dopamine 2-like receptor (Dop2R) in promoting egg laying. We replicated the Dop2R effect in an independently collected productivity dataset and showed that the effect of the Dop2R variant was mediated in part by regulatory gene expression variation. This study demonstrates the strong potential of genome-wide association analysis in this diverse panel of inbred strains and subsequent functional analyses for understanding the genetic architecture of fitness traits.

Indirect genetic effects for social network structure in Drosophila melanogaster

The position an individual holds in a social network is dependent on both its direct and indirect social interactions. Because social network position is dependent on the actions and interactions of conspecifics, it is likely that the genotypic composition of individuals within a social group impacts individuals' network positions. However, we know very little about whether social network positions have a genetic basis, and even less about how the genotypic makeup of a social group impacts network positions and structure. With ample evidence indicating that network positions influence various fitness metrics, studying how direct and indirect genetic effects shape network positions is crucial for furthering our understanding of how the social environment can respond to selection and evolve. Using replicate genotypes of Drosophila melanogaster fruit flies, we created social groups that varied in their genotypic makeup. Social groups were videoed, and networks were generated using motion-tracking software. We found that both an individual's own genotype and the genotypes of conspecifics in its social group affect its position within a social network. These findings provide an early example of how indirect genetic effects and social network theory can be linked, and shed new light on how quantitative genetic variation shapes the structure of social groups. This article is part of a discussion meeting issue 'Collective behaviour through time'.

The importance of reaction norms in dietary restriction and ageing research

Ageing research has progressed rapidly through our ability to modulate the ageing process. Pharmacological and dietary treatments can increase lifespan and have been instrumental in our understanding of the mechanisms of ageing. Recently, several studies have reported genetic variance in response to these anti-ageing interventions, questioning their universal application and making a case for personalised medicine in our field. As an extension of these findings the response to dietary restriction was found to not be repeatable when the same genetic mouse lines were retested. We show here that this effect is more widespread with the response to dietary restriction also showing low repeatability across genetic lines in the fly (Drosophila melanogaster). We further argue that variation in reaction norms, the relationship between dose and response, can explain such conflicting findings in our field. We simulate genetic variance in reaction norms and show that such variation can: 1) lead to over- or under-estimation of treatment responses, 2) dampen the response measured if a genetically heterogeneous population is studied, and 3) illustrate that genotype-by-dose-by-environment interactions can lead to low repeatability of DR and potentially other anti-ageing interventions. We suggest that putting experimental biology and personalised geroscience in a reaction norm framework will aid progress in ageing research.