The genetic basis of variation in Drosophila melanogaster mating behavior.

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

Pleiotropy, epistasis and the genetic architecture of quantitative traits.

Pleiotropy (whereby one genetic polymorphism affects multiple traits) and epistasis (whereby non-linear interactions between genetic polymorphisms affect the same trait) are fundamental aspects of the genetic architecture of quantitative traits. Recent advances in the ability to characterize the effects of polymorphic variants on molecular and organismal phenotypes in human and model organism populations have revealed the prevalence of pleiotropy and unexpected shared molecular genetic bases among quantitative traits, including diseases. By contrast, epistasis is common between polymorphic loci associated with quantitative traits in model organisms, such that alleles at one locus have different effects in different genetic backgrounds, but is rarely observed for human quantitative traits and common diseases. Here, we review the concepts and recent inferences about pleiotropy and epistasis, and discuss factors that contribute to similarities and differences between the genetic architecture of quantitative traits in model organisms and humans.

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.

Genetic and Genomic Analyses of Drosophila melanogaster Models of Chromatin Modification Disorders.

Switch/Sucrose Non-Fermentable (SWI/SNF)-related intellectual disability disorders (SSRIDDs) and Cornelia de Lange syndrome are rare syndromic neurodevelopmental disorders with overlapping clinical phenotypes. SSRIDDs are associated with the BAF (Brahma-Related Gene-1 Associated Factor) complex, whereas CdLS is a disorder of chromatin modification associated with the cohesin complex. Here, we used RNA interference in Drosophila melanogaster to reduce expression of six genes (brm, osa, Snr1, SMC1, SMC3, vtd) orthologous to human genes associated with SSRIDDs and CdLS. These fly models exhibit changes in sleep, activity, startle behavior (a proxy for sensorimotor integration) and brain morphology. Whole genome RNA sequencing identified 9,657 differentially expressed genes (FDR < 0.05), 156 of which are differentially expressed in both sexes in SSRIDD- and CdLS-specific analyses, including Bap60, which is orthologous to SMARCD1, a SSRIDD-associated BAF component, k-means clustering reveals genes co-regulated within and across SSRIDD and CdLS fly models. RNAi-mediated reduction of expression of six genes co-regulated with focal genes brm, osa, and/or Snr1 recapitulated changes in behavior of the focal genes. Based on the assumption that fundamental biological processes are evolutionarily conserved, Drosophila models can be used to understand underlying molecular effects of variants in chromatin-modification pathways and may aid in discovery of drugs that ameliorate deleterious phenotypic effects.

Genetic and genomic analyses of Drosophila melanogaster models of chromatin modification disorders.

Switch/sucrose nonfermentable (SWI/SNF)-related intellectual disability disorders (SSRIDDs) and Cornelia de Lange syndrome are rare syndromic neurodevelopmental disorders with overlapping clinical phenotypes. SSRIDDs are associated with the BAF (Brahma-Related Gene-1 associated factor) complex, whereas CdLS is a disorder of chromatin modification associated with the cohesin complex. Here, we used RNA interference in Drosophila melanogaster to reduce the expression of six genes (brm, osa, Snr1, SMC1, SMC3, vtd) orthologous to human genes associated with SSRIDDs and CdLS. These fly models exhibit changes in sleep, activity, startle behavior (a proxy for sensorimotor integration), and brain morphology. Whole genome RNA sequencing identified 9,657 differentially expressed genes (FDR < 0.05), 156 of which are differentially expressed in both sexes in SSRIDD- and CdLS-specific analyses, including Bap60, which is orthologous to SMARCD1, an SSRIDD-associated BAF component. k-means clustering reveals genes co-regulated within and across SSRIDD and CdLS fly models. RNAi-mediated reduction of expression of six genes co-regulated with focal genes brm, osa, and/or Snr1 recapitulated changes in the behavior of the focal genes. Based on the assumption that fundamental biological processes are evolutionarily conserved, Drosophila models can be used to understand underlying molecular effects of variants in chromatin-modification pathways and may aid in the discovery of drugs that ameliorate deleterious phenotypic effects.

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.

The genetic architecture of behavioral canalization.

Behaviors are components of fitness and contribute to adaptive evolution. Behaviors represent the interactions of an organism with its environment, yet innate behaviors display robustness in the face of environmental change, which we refer to as 'behavioral canalization'. We hypothesize that positive selection of hub genes of genetic networks stabilizes the genetic architecture for innate behaviors by reducing variation in the expression of interconnected network genes. Robustness of these stabilized networks would be protected from deleterious mutations by purifying selection or suppressing epistasis. We propose that, together with newly emerging favorable mutations, epistatically suppressed mutations can generate a reservoir of cryptic genetic variation that could give rise to decanalization when genetic backgrounds or environmental conditions change to allow behavioral adaptation.

Pleiotropic fitness effects of a Drosophila odorant-binding protein.

Insect odorant-binding proteins (OBPs) are members of a rapidly evolving multigene family traditionally thought to facilitate chemosensation. However, studies on Drosophila have shown that members of this family have evolved functions beyond chemosensation, as evident from their expression in reproductive tissues and the brain. Previous studies implicated diverse functions of Obp56h, a member of the largest gene cluster of the D. melanogaster Obp repertoire. Here, we examined the effect of CRISPR/Cas9-mediated deletion of Obp56h on 2 fitness phenotypes, on resistance to starvation stress and heat stress, and on locomotion and sleep phenotypes. Obp56h-/- mutants show a strong sexually dimorphic effect on starvation stress survival, with females being more resistant to starvation stress than the control. In contrast, Obp56h-/- females, but not males, are highly sensitive to heat stress. Both sexes show changes in locomotion and sleep patterns. Transcriptional profiling of RNA from heads of Obp56h-/- flies and the wildtype control reveals differentially expressed genes, including gene products associated with antimicrobial immune responses and members of the Turandot family of stress-induced secreted peptides. In addition, differentially expressed genes of unknown function were identified in both sexes. Genes encoding components of the mitochondrial electron transport chain, cuticular proteins, gene products associated with regulation of feeding behavior (Lst and CCHa2), ribosomal proteins, lncRNAs, snoRNAs, tRNAs, and snRNAs show changes in transcript abundances in Obp56h-/- females. These differentially expressed genes are likely to contribute to Obp56h-mediated effects on the diverse phenotypes that arise upon deletion of this OBP.

Pleiotropic fitness effects of the lncRNA Uhg4 in Drosophila melanogaster.

BACKGROUND: Long noncoding RNAs (lncRNAs) are a diverse class of RNAs that are critical for gene regulation, DNA repair, and splicing, and have been implicated in development, stress response, and cancer. However, the functions of many lncRNAs remain unknown. In Drosophila melanogaster, U snoRNA host gene 4 (Uhg4) encodes an antisense long noncoding RNA that is host to seven small nucleolar RNAs (snoRNAs). Uhg4 is expressed ubiquitously during development and in all adult tissues, with maximal expression in ovaries; however, it has no annotated function(s). RESULTS: We used CRISPR-Cas9 germline gene editing to generate multiple deletions spanning the promoter region and first exon of Uhg4. Females showed arrested egg development and both males and females were sterile. In addition, Uhg4 deletion mutants showed delayed development and decreased viability, and changes in sleep and responses to stress. Whole-genome RNA sequencing of Uhg4 deletion flies and their controls identified co-regulated genes and genetic interaction networks associated with Uhg4. Gene ontology analyses highlighted a broad spectrum of biological processes, including regulation of transcription and translation, morphogenesis, and stress response. CONCLUSION: Uhg4 is a lncRNA essential for reproduction with pleiotropic effects on multiple fitness traits.

Gregor Mendel's legacy in quantitative genetics.

Gregor Mendel's discovery of the laws of segregation and independent assortment and his inference of the existence of non-mendelian interactions between loci remain at the heart of today's explorations of the genetic architecture of quantitative traits.

Modulation of the Drosophila transcriptome by developmental exposure to alcohol.

BACKGROUND: Prenatal exposure to ethanol can cause fetal alcohol spectrum disorder (FASD), a prevalent, preventable pediatric disorder. Identifying genetic risk alleles for FASD is challenging since time, dose, and frequency of exposure are often unknown, and manifestations of FASD are diverse and evident long after exposure. Drosophila melanogaster is an excellent model to study the genetic basis of the effects of developmental alcohol exposure since many individuals of the same genotype can be reared under controlled environmental conditions. RESULTS: We used 96 sequenced, wild-derived inbred lines from the Drosophila melanogaster Genetic Reference Panel (DGRP) to profile genome-wide transcript abundances in young adult flies that developed on ethanol-supplemented medium or standard culture medium. We found substantial genetic variation in gene expression in response to ethanol with extensive sexual dimorphism. We constructed sex-specific genetic networks associated with alcohol-dependent modulation of gene expression that include protein-coding genes, Novel Transcribed Regions (NTRs, postulated to encode long non-coding RNAs) and female-specific coordinated regulation of snoRNAs that regulate pseudouridylation of ribosomal RNA. We reared DGRP lines which showed extreme upregulation or downregulation of snoRNA expression during developmental alcohol exposure on standard or ethanol supplemented medium and demonstrated that developmental exposure to ethanol has genotype-specific effects on adult locomotor activity and sleep. CONCLUSIONS: There is significant and sex-specific natural genetic variation in the transcriptional response to developmental exposure to ethanol in Drosophila that comprises networks of genes affecting nervous system development and ethanol metabolism as well as networks of regulatory non-coding RNAs.

Ibrutinib as a potential therapeutic for cocaine use disorder.

Cocaine use presents a worldwide public health problem with high socioeconomic cost. No current pharmacologic treatments are available for cocaine use disorder (CUD) or cocaine toxicity. To explore pharmaceutical treatments for tthis disorder and its sequelae we analyzed gene expression data from post-mortem brain tissue of individuals with CUD who died from cocaine-related causes with matched cocaine-free controls (n = 71, Mage = 39.9, 100% male, 49% with CUD, 3 samples/brain regions). To match molecular signatures from brain pathology with potential therapeutics, we leveraged the L1000 database honing in on neuronal mRNA profiles of 825 repurposable compounds (e.g., FDA approved). We identified 16 compounds that were negatively associated with CUD gene expression patterns across all brain regions (padj < 0.05), all of which outperformed current targets undergoing clinical trials for CUD (all padj > 0.05). An additional 43 compounds were positively associated with CUD expression. We performed an in silico follow-up potential therapeutics using independent transcriptome-wide in vitro (neuronal cocaine exposure; n = 18) and in vivo (mouse cocaine self-administration; n = 12-15) datasets to prioritize candidates for experimental validation. Among these medications, ibrutinib was consistently linked with the molecular profiles of both neuronal cocaine exposure and mouse cocaine self-administration. We assessed the therapeutic efficacy of ibrutinib using the Drosophila melanogaster model. Ibrutinib reduced cocaine-induced startle response and cocaine-induced seizures (n = 61-142 per group; sex: 51% female), despite increasing cocaine consumption. Our results suggest that ibrutinib could be used for the treatment of cocaine use disorder.

Epistasis for head morphology in Drosophila melanogaster.

Epistasis-gene-gene interaction-is common for mutations with large phenotypic effects in humans and model organisms. Epistasis impacts quantitative genetic models of speciation, response to natural and artificial selection, genetic mapping, and personalized medicine. However, the existence and magnitude of epistasis between alleles with small quantitative phenotypic effects are controversial and difficult to assess. Here, we use the Drosophila melanogaster Genetic Reference Panel of sequenced inbred lines to evaluate the magnitude of naturally occurring epistasis modifying the effects of mutations in jing and inv, two transcription factors that have subtle quantitative effects on head morphology as homozygotes. We find significant epistasis for both mutations and performed single marker genome-wide association analyses to map candidate modifier variants and loci affecting head morphology. A subset of these loci was significantly enriched for a known genetic interaction network, and mutations of the candidate epistatic modifier loci also affect head morphology.

Physiological and metabolomic consequences of reduced expression of the Drosophila brummer triglyceride Lipase.

Disruption of lipolysis has widespread effects on intermediary metabolism and organismal phenotypes. Defects in lipolysis can be modeled in Drosophila melanogaster through genetic manipulations of brummer (bmm), which encodes a triglyceride lipase orthologous to mammalian Adipose Triglyceride Lipase. RNAi-mediated knock-down of bmm in all tissues or metabolic specific tissues results in reduced locomotor activity, altered sleep patterns and reduced lifespan. Metabolomic analysis on flies in which bmm is downregulated reveals a marked reduction in medium chain fatty acids, long chain saturated fatty acids and long chain monounsaturated and polyunsaturated fatty acids, and an increase in diacylglycerol levels. Elevated carbohydrate metabolites and tricarboxylic acid intermediates indicate that impairment of fatty acid mobilization as an energy source may result in upregulation of compensatory carbohydrate catabolism. bmm downregulation also results in elevated levels of serotonin and dopamine neurotransmitters, possibly accounting for the impairment of locomotor activity and sleep patterns. Physiological phenotypes and metabolomic changes upon reduction of bmm expression show extensive sexual dimorphism. Altered metabolic states in the Drosophila model are relevant for understanding human metabolic disorders, since pathways of intermediary metabolism are conserved across phyla.

Developmental Alcohol Exposure in Drosophila: Effects on Adult Phenotypes and Gene Expression in the Brain.

Fetal alcohol exposure can lead to developmental abnormalities, intellectual disability, and behavioral changes, collectively termed fetal alcohol spectrum disorder (FASD). In 2015, the Centers for Disease Control found that 1 in 10 pregnant women report alcohol use and more than 3 million women in the USA are at risk of exposing their developing fetus to alcohol. Drosophila melanogaster is an excellent genetic model to study developmental effects of alcohol exposure because many individuals of the same genotype can be reared rapidly and economically under controlled environmental conditions. Flies exposed to alcohol undergo physiological and behavioral changes that resemble human alcohol-related phenotypes. Here, we show that adult flies that developed on ethanol-supplemented medium have decreased viability, reduced sensitivity to ethanol, and disrupted sleep and activity patterns. To assess the effects of exposure to alcohol during development on brain gene expression, we performed single cell RNA sequencing and resolved cell clusters with differentially expressed genes which represent distinct neuronal and glial populations. Differential gene expression showed extensive sexual dimorphism with little overlap between males and females. Gene expression differences following developmental alcohol exposure were similar to previously reported differential gene expression following cocaine consumption, suggesting that common neural substrates respond to both drugs. Genes associated with glutathione metabolism, lipid transport, glutamate and GABA metabolism, and vision feature in sexually dimorphic global multi-cluster interaction networks. Our results provide a blueprint for translational studies on alcohol-induced effects on gene expression in the brain that may contribute to or result from FASD in human populations.

A High Throughput Microplate Feeder Assay for Quantification of Consumption in Drosophila.

Quantifying food intake in Drosophila is used to study the genetic and physiological underpinnings of consumption-associated traits, their environmental factors, and the toxicological and pharmacological effects of numerous substances. Few methods currently implemented are amenable to high throughput measurement. The Microplate Feeder Assay (MFA) was developed for quantifying the consumption of liquid food for individual flies using absorbance. In this assay, flies consume liquid food medium from select wells of a 1536-well microplate. By incorporating a dilute tracer dye into the liquid food medium and loading a known volume into each well, absorbance measurements of the well acquired before and after consumption reflect the resulting change in volume (i.e., volume consumed). To enable high throughput analysis with this method, a 3D-printed coupler was designed that allows flies to be sorted individually into 96-well microplates. This device precisely orients 96- and 1536-well microplates to give each fly access to up to 4 wells for consumption, thus enabling food preference quantification in addition to regular consumption. Furthermore, the device has barrier strips that toggle between open and closed positions to allow for controlled containment and release of a column of samples at a time. This method enables high throughput measurements of consumption of aqueous solutions by many flies simultaneously. It also has the potential to be adapted to other insects and to screen consumption of nutrients, toxins, or pharmaceuticals.

Genetic basis of variation in cocaine and methamphetamine consumption in outbred populations of Drosophila melanogaster.

We used Drosophila melanogaster to map the genetic basis of naturally occurring variation in voluntary consumption of cocaine and methamphetamine. We derived an outbred advanced intercross population (AIP) from 37 sequenced inbred wild-derived lines of the Drosophila melanogaster Genetic Reference Panel (DGRP), which are maximally genetically divergent, have minimal residual heterozygosity, are not segregating for common inversions, and are not infected with Wolbachia pipientis We assessed consumption of sucrose, methamphetamine-supplemented sucrose, and cocaine-supplemented sucrose and found considerable phenotypic variation for consumption of both drugs, in both sexes. We performed whole-genome sequencing and extreme quantitative trait locus (QTL) mapping on the top 10% of consumers for each replicate, sex, and condition and an equal number of randomly selected flies. We evaluated changes in allele frequencies among high consumers and control flies and identified 3,033 variants significantly (P < 1.9 × 10-8) associated with increased consumption, located in or near 1,962 genes. Many of these genes are associated with nervous system development and function, and 77 belong to a known gene-gene interaction subnetwork. We assessed the effects of RNA interference (RNAi) on drug consumption for 22 candidate genes; 17 had a significant effect in at least one sex. We constructed allele-specific AIPs that were homozygous for alternative candidate alleles for 10 single-nucleotide polymorphisms (SNPs) and measured average consumption for each population; 9 SNPs had significant effects in at least one sex. The genetic basis of voluntary drug consumption in Drosophila is polygenic and implicates genes with human orthologs and associated variants with sex- and drug-specific effects.

The Drosophila brain on cocaine at single-cell resolution.

Whereas the neurological effects of cocaine have been well documented, effects of acute cocaine consumption on genome-wide gene expression across the brain remain largely unexplored. This question cannot be readily addressed in humans but can be approached using the Drosophila melanogaster model, where gene expression in the entire brain can be surveyed at once. Flies exposed to cocaine show impaired locomotor activity, including climbing behavior and startle response (a measure of sensorimotor integration), and increased incidence of seizures and compulsive grooming. To identify specific cell populations that respond to acute cocaine exposure, we analyzed single-cell transcriptional responses in duplicate samples of flies that consumed fixed amounts of sucrose or sucrose supplemented with cocaine, in both sexes. Unsupervised clustering of the transcriptional profiles of a total of 86,224 cells yielded 36 distinct clusters. Annotation of clusters based on gene markers revealed that all major cell types (neuronal and glial) as well as neurotransmitter types from most brain regions were represented. The brain transcriptional responses to cocaine showed profound sexual dimorphism and were considerably more pronounced in males than females. Differential expression analysis within individual clusters indicated cluster-specific responses to cocaine. Clusters corresponding to Kenyon cells of the mushroom bodies and glia showed especially large transcriptional responses following cocaine exposure. Cluster specific coexpression networks and global interaction networks revealed a diverse array of cellular processes affected by acute cocaine exposure. These results provide an atlas of sexually dimorphic cocaine-modulated gene expression in a model brain.

Functional Diversification, Redundancy, and Epistasis among Paralogs of the Drosophila melanogaster Obp50a-d Gene Cluster.

Large multigene families, such as the insect odorant-binding proteins (OBPs), are thought to arise through functional diversification after repeated gene duplications. Whereas many OBPs function in chemoreception, members of this family are also expressed in tissues outside chemosensory organs. Paralogs of the Obp50 gene cluster are expressed in metabolic and male reproductive tissues, but their functions and interrelationships remain unknown. Here, we report the genetic dissection of four members of the Obp50 cluster, which are in close physical proximity without intervening genes. We used CRISPR technology to excise the entire cluster while introducing a PhiC31 reintegration site to reinsert constructs in which different combinations of the constituent Obp genes were either intact or rendered inactive. We performed whole transcriptome sequencing and assessed sexually dimorphic changes in transcript abundances (transcriptional niches) associated with each gene-edited genotype. Using this approach, we were able to estimate redundancy, additivity, diversification, and epistasis among Obp50 paralogs. We analyzed the effects of gene editing of this cluster on organismal phenotypes and found a significant skewing of sex ratios attributable to Obp50a, and sex-specific effects on starvation stress resistance attributable to Obp50d. Thus, there is functional diversification within the Obp50 cluster with Obp50a contributing to development and Obp50d to stress resistance. The deletion-reinsertion approach we applied to the Obp50 cluster provides a general paradigm for the genetic dissection of paralogs of multigene families.

Genotype by environment interaction for gene expression in Drosophila melanogaster.

The genetics of phenotypic responses to changing environments remains elusive. Using whole-genome quantitative gene expression as a model, here we study how the genetic architecture of regulatory variation in gene expression changed in a population of fully sequenced inbred Drosophila melanogaster strains when flies developed in different environments (25 °C and 18 °C). We find a substantial fraction of the transcriptome exhibited genotype by environment interaction, implicating environmentally plastic genetic architecture of gene expression. Genetic variance in expression increases at 18 °C relative to 25 °C for most genes that have a change in genetic variance. Although the majority of expression quantitative trait loci (eQTLs) for the gene expression traits in the two environments are shared and have similar effects, analysis of the environment-specific eQTLs reveals enrichment of binding sites for two transcription factors. Finally, although genotype by environment interaction in gene expression could potentially disrupt genetic networks, the co-expression networks are highly conserved across environments. Genes with higher network connectivity are under stronger stabilizing selection, suggesting that stabilizing selection on expression plays an important role in promoting network robustness.