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.

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.

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.

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.

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.

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.

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.

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.

High-Throughput Method for Measuring Alcohol Sedation Time of Individual Drosophila melanogaster

Drosophila melanogaster provides an excellent model to study the genetic underpinnings of alcohol sensitivity. In contrast to studies in human populations, the Drosophila model allows strict control over genetic background, and virtually unlimited numbers of individuals of the same genotype can be reared rapidly under well-controlled environmental conditions without regulatory restrictions and at relatively low cost. Flies exposed to ethanol undergo physiological and behavioral changes that resemble human alcohol intoxication, including loss of postural control, sedation, and development of tolerance. Here, we describe a simple, low-cost, high-throughput assay for assessing alcohol sedation sensitivity in large numbers of single flies. The assay is based on video recording of single flies introduced without anesthesia in 24-well cell culture plates in a set-up that enables synchronous initiation of alcohol exposure. The system enables a single person to collect individual ethanol sedation data on as many as 2,000 flies within an 8 h work period. The assay can, in principle, be extended to assess the effects of exposure to any volatile substance and applied to measure effects of acute toxicity of volatiles on other insects, including other fly species.

Gene expression networks in the Drosophila Genetic Reference Panel

A major challenge in modern biology is to understand how naturally occurring variation in DNA sequences affects complex organismal traits through networks of intermediate molecular phenotypes. This question is best addressed in a genetic mapping population in which all molecular polymorphisms are known and for which molecular endophenotypes and complex traits are assessed on the same genotypes. Here, we performed deep RNA sequencing of 200 Drosophila Genetic Reference Panel inbred lines with complete genome sequences and for which phenotypes of many quantitative traits have been evaluated. We mapped expression quantitative trait loci for annotated genes, novel transcribed regions, transposable elements, and microbial species. We identified host variants that affect expression of transposable elements, independent of their copy number, as well as microbiome composition. We constructed sex-specific expression quantitative trait locus regulatory networks. These networks are enriched for novel transcribed regions and target genes in heterochromatin and euchromatic regions of reduced recombination, as well as genes regulating transposable element expression. This study provides new insights regarding the role of natural genetic variation in regulating gene expression and generates testable hypotheses for future functional analyses.

Context-dependent genetic architecture of Drosophila life span

Understanding the genetic basis of variation in life span is a major challenge that is difficult to address in human populations. Evolutionary theory predicts that alleles affecting natural variation in life span will have properties that enable them to persist in populations at intermediate frequencies, such as late-life-specific deleterious effects, antagonistic pleiotropic effects on early and late-age fitness components, and/or sex- and environment-specific or antagonistic effects. Here, we quantified variation in life span in males and females reared in 3 thermal environments for the sequenced, inbred lines of the Drosophila melanogaster Genetic Reference Panel (DGRP) and an advanced intercross outbred population derived from a subset of DGRP lines. Quantitative genetic analyses of life span and the micro-environmental variance of life span in the DGRP revealed significant genetic variance for both traits within each sex and environment, as well as significant genotype-by-sex interaction (GSI) and genotype-by-environment interaction (GEI). Genome-wide association (GWA) mapping in both populations implicates over 2,000 candidate genes with sex- and environment-specific or antagonistic pleiotropic allelic effects. Over 1,000 of these genes are associated with variation in life span in other D. melanogaster populations. We functionally assessed the effects of 15 candidate genes using RNA interference (RNAi): all affected life span and/or micro-environmental variance of life span in at least one sex and environment and exhibited sex-and environment-specific effects. Our results implicate novel candidate genes affecting life span and suggest that variation for life span may be maintained by variable allelic effects in heterogeneous environments.

Chemosensation and Evolution of Drosophila Host Plant Selection

The ability to respond to chemosensory cues is critical for survival of most organisms. Among insects, Drosophila melanogaster has the best characterized olfactory system, and the availability of genome sequences of 30 Drosophila species provides an ideal scenario for studies on evolution of chemosensation. Gene duplications of chemoreceptor genes allow for functional diversification of the rapidly evolving chemoreceptor repertoire. Although some species of the genus Drosophila are generalists for host plant selection, rapid evolution of olfactory receptors, gustatory receptors, odorant-binding proteins, and cytochrome P450s has enabled diverse host specializations of different members of the genus. Here, I review diversification of the chemoreceptor repertoire among members of the genus Drosophila along with co-evolution of detoxification mechanisms that may have enabled occupation of diverse host plant ecological niches.

Evolution of Reproductive Behavior

Behaviors associated with reproduction are major contributors to the evolutionary success of organisms and are subject to many evolutionary forces, including natural and sexual selection, and sexual conflict. Successful reproduction involves a range of behaviors, from finding an appropriate mate, courting, and copulation, to the successful production and (in oviparous animals) deposition of eggs following mating. As a consequence, behaviors and genes associated with reproduction are often under strong selection and evolve rapidly. Courtship rituals in flies follow a multimodal pattern, mediated through visual, chemical, tactile, and auditory signals. Premating behaviors allow males and females to assess the species identity, reproductive state, and condition of their partners. Conflicts between the "interests" of individual males, and/or between the reproductive strategies of males and females, often drive the evolution of reproductive behaviors. For example, seminal proteins transmitted by males often show evidence of rapid evolution, mediated by positive selection. Postmating behaviors, including the selection of oviposition sites, are highly variable and Drosophila species span the spectrum from generalists to obligate specialists. Chemical recognition features prominently in adaptation to host plants for feeding and oviposition. Selection acting on variation in pre-, peri-, and postmating behaviors can lead to reproductive isolation and incipient speciation. Response to selection at the genetic level can include the expansion of gene families, such as those for detecting pheromonal cues for mating, or changes in the expression of genes leading to visual cues such as wing spots that are assessed during mating. Here, we consider the evolution of reproductive behavior in Drosophila at two distinct, yet complementary, scales. Some studies take a microevolutionary approach, identifying genes and networks involved in reproduction, and then dissecting the genetics underlying complex behaviors in D. melanogaster Other studies take a macroevolutionary approach, comparing reproductive behaviors across the genus Drosophila and how these might correlate with environmental cues. A full synthesis of this field will require unification across these levels.

Evolution of Epistatic Networks and the Genetic Basis of Innate Behaviors

Instinctive behaviors are genetically programmed behaviors that occur independent of experience. How genetic programs that give rise to the manifestation of such behaviors evolve remains an unresolved question. I propose that evolution of species-specific innate behaviors is accomplished through progressive modifications of pre-existing genetic networks composed of allelic variants. I hypothesize that changes in frequencies of one or more constituent allelic variants within the network leads to changes in gene network connectivity and the emergence of a reorganized network that can support the emergence of a novel behavioral phenotype and becomes stabilized when key allelic variants are driven to fixation.

Systems genetics of the Drosophila metabolome

How effects of DNA sequence variants are transmitted through intermediate endophenotypes to modulate organismal traits remains a central question in quantitative genetics. This problem can be addressed through a systems approach in a population in which genetic polymorphisms, gene expression traits, metabolites, and complex phenotypes can be evaluated on the same genotypes. Here, we focused on the metabolome, which represents the most proximal link between genetic variation and organismal phenotype, and quantified metabolite levels in 40 lines of the Drosophila melanogaster Genetic Reference Panel. We identified sex-specific modules of genetically correlated metabolites and constructed networks that integrate DNA sequence variation and variation in gene expression with variation in metabolites and organismal traits, including starvation stress resistance and male aggression. Finally, we asked to what extent SNPs and metabolites can predict trait phenotypes and generated trait- and sex-specific prediction models that provide novel insights about the metabolomic underpinnings of complex phenotypes.

Genetics of cocaine and methamphetamine consumption and preference in Drosophila melanogaster

Illicit use of psychostimulants, such as cocaine and methamphetamine, constitutes a significant public health problem. Whereas neural mechanisms that mediate the effects of these drugs are well-characterized, genetic factors that account for individual variation in susceptibility to substance abuse and addiction remain largely unknown. Drosophila melanogaster can serve as a translational model for studies on substance abuse, since flies have a dopamine transporter that can bind cocaine and methamphetamine, and exposure to these compounds elicits effects similar to those observed in people, suggesting conserved evolutionary mechanisms underlying drug responses. Here, we used the D. melanogaster Genetic Reference Panel to investigate the genetic basis for variation in psychostimulant drug consumption, to determine whether similar or distinct genetic networks underlie variation in consumption of cocaine and methamphetamine, and to assess the extent of sexual dimorphism and effect of genetic context on variation in voluntary drug consumption. Quantification of natural genetic variation in voluntary consumption, preference, and change in consumption and preference over time for cocaine and methamphetamine uncovered significant genetic variation for all traits, including sex-, exposure- and drug-specific genetic variation. Genome wide association analyses identified both shared and drug-specific candidate genes, which could be integrated in genetic interaction networks. We assessed the effects of ubiquitous RNA interference (RNAi) on consumption behaviors for 34 candidate genes: all affected at least one behavior. Finally, we utilized RNAi knockdown in the nervous system to implicate dopaminergic neurons and the mushroom bodies as part of the neural circuitry underlying experience-dependent development of drug preference.

Illicit use of cocaine and methamphetamine is a major public health problem. Whereas the neurological effects of these drugs are well characterized, it remains challenging to determine genetic risk factors for substance abuse in human populations. The fruit fly, Drosophila melanogaster, presents an excellent model for identifying evolutionarily conserved genes that affect drug consumption, since genetic background and exposure can be controlled precisely. We took advantage of natural variation in a panel of inbred wild derived fly lines with complete genome sequences to assess the extent of genetic variation among these lines for voluntary consumption of cocaine and methamphetamine and to explore whether some genetic backgrounds might show experience-dependent development of drug preference. The drug consumption traits were highly variable among the lines with strong sex-, drug- and exposure time-specific components. We identified candidate genes and gene networks associated with variation in consumption of cocaine and methamphetamine and development of drug preference. Using tissue-specific suppression of gene expression, we were able to functionally implicate candidate genes that affected at least one consumption trait in at least one drug and sex. In humans, the mesolimbic dopaminergic projection plays a role in drug addiction. We asked whether in Drosophila the mushroom bodies could play an analogous role, as they are integrative brain centers associated with experience-dependent learning. Indeed, our results suggest that variation in consumption and development of preference for both cocaine and methamphetamine is mediated, at least in part, through a neural network that comprises dopaminergic projections to the mushroom bodies.

Genetics of alcohol consumption in Drosophila melanogaster

Individual variation in alcohol consumption in human populations is determined by genetic, environmental, social and cultural factors. In contrast to humans, genetic contributions to complex behavioral phenotypes can be readily dissected in Drosophila, where both the genetic background and environment can be controlled and behaviors quantified through simple high-throughput assays. Here, we measured voluntary consumption of ethanol in ∼3000 individuals of each sex from an advanced intercross population derived from 37 lines of the Drosophila melanogaster Genetic Reference Panel. Extreme quantitative trait loci mapping identified 385 differentially segregating allelic variants located in or near 291 genes at P < 10^-8 . The effects of single nucleotide polymorphisms associated with voluntary ethanol consumption are sex-specific, as found for other alcohol-related phenotypes. To assess causality, we used RNA interference knockdown or P{MiET1} mutants and their corresponding controls and functionally validated 86% of candidate genes in at least one sex. We constructed a genetic network comprised of 23 genes along with a separate trio and a pair of connected genes. Gene ontology analyses showed enrichment of developmental genes, including development of the nervous system. Furthermore, a network of human orthologs showed enrichment for signal transduction processes, protein metabolism and developmental processes, including nervous system development. Our results show that the genetic architecture that underlies variation in voluntary ethanol consumption is sexually dimorphic and partially overlaps with genetic factors that control variation in feeding behavior and alcohol sensitivity. This integrative genetic architecture is rooted in evolutionarily conserved features that can be extrapolated to human genetic interaction networks.

A Drosophila model for toxicogenomics: Genetic variation in susceptibility to heavy metal exposure

The genetic factors that give rise to variation in susceptibility to environmental toxins remain largely unexplored. Studies on genetic variation in susceptibility to environmental toxins are challenging in human populations, due to the variety of clinical symptoms and difficulty in determining which symptoms causally result from toxic exposure; uncontrolled environments, often with exposure to multiple toxicants; and difficulty in relating phenotypic effect size to toxic dose, especially when symptoms become manifest with a substantial time lag. Drosophila melanogaster is a powerful model that enables genome-wide studies for the identification of allelic variants that contribute to variation in susceptibility to environmental toxins, since the genetic background, environmental rearing conditions and toxic exposure can be precisely controlled. Here, we used extreme QTL mapping in an outbred population derived from the D. melanogaster Genetic Reference Panel to identify alleles associated with resistance to lead and/or cadmium, two ubiquitous environmental toxins that present serious health risks. We identified single nucleotide polymorphisms (SNPs) associated with variation in resistance to both heavy metals as well as SNPs associated with resistance specific to each of them. The effects of these SNPs were largely sex-specific. We applied mutational and RNAi analyses to 33 candidate genes and functionally validated 28 of them. We constructed networks of candidate genes as blueprints for orthologous networks of human genes. The latter not only provided functional contexts for known human targets of heavy metal toxicity, but also implicated novel candidate susceptibility genes. These studies validate Drosophila as a translational toxicogenomics gene discovery system.

Although physiological effects of environmental toxins are well documented, we know little about the genetic factors that determine individual variation in susceptibility to toxins. Such information is difficult to obtain in human populations due to heterogeneity in genetic background and environmental exposure, and the diversity of symptoms and time lag with which they appear after toxic exposure. Here, we show that the fruit fly, Drosophila, can serve as a powerful genetic model system to elucidate the genetic underpinnings that contribute to individual variation in resistance to toxicity, using lead and cadmium exposure as an experimental paradigm. We identified genes that harbor genetic variants that contribute to individual variation in resistance to heavy metal exposure. Furthermore, we constructed genetic networks on which we could superimpose human counterparts of Drosophila genes. We were able to place human genes previously implicated in heavy metal toxicity in biological context and identify novel targets for heavy metal toxicity. Thus, we demonstrate that based on evolutionary conservation of fundamental biological processes, we can use Drosophila as a powerful translational model for toxicogenomics studies.

Regulation of Drosophila Lifespan by bellwether Promoter Alleles

Longevity varies among individuals, but how natural genetic variation contributes to variation in lifespan is poorly understood. Drosophila melanogaster presents an advantageous model system to explore the genetic underpinnings of longevity, since its generation time is brief and both the genetic background and rearing environment can be precisely controlled. The bellwether (blw) gene encodes the α subunit of mitochondrial ATP synthase. Since metabolic rate may influence lifespan, we investigated whether alternative haplotypes in the blw promoter affect lifespan when expressed in a co-isogenic background. We amplified 521 bp upstream promoter sequences containing alternative haplotypes and assessed promoter activity both in vitro and in vivo using a luciferase reporter system. The AG haplotype showed significantly greater expression of luciferase than the GT haplotype. We then overexpressed a blw cDNA construct driven by either the AG or GT haplotype promoter in transgenic flies and showed that the AG haplotype also results in greater blw cDNA expression and a significant decrease in lifespan relative to the GT promoter haplotype, in male flies only. Thus, our results show that naturally occurring regulatory variants of blw affect lifespan in a sex-specific manner.

Obp56h Modulates Mating Behavior in Drosophila melanogaster

Social interactions in insects are driven by conspecific chemical signals that are detected via olfactory and gustatory neurons. Odorant binding proteins (Obps) transport volatile odorants to chemosensory receptors, but their effects on behaviors remain poorly characterized. Here, we report that RNAi knockdown of Obp56h gene expression in Drosophila melanogaster enhances mating behavior by reducing courtship latency. The change in mating behavior that results from inhibition of Obp56h expression is accompanied by significant alterations in cuticular hydrocarbon (CHC) composition, including reduction in 5-tricosene (5-T), an inhibitory sex pheromone produced by males that increases copulation latency during courtship. Whole genome RNA sequencing confirms that expression of Obp56h is virtually abolished in Drosophila heads. Inhibition of Obp56h expression also affects expression of other chemoreception genes, including upregulation of lush in both sexes and Obp83ef in females, and reduction in expression of Obp19b and Or19b in males. In addition, several genes associated with lipid metabolism, which underlies the production of cuticular hydrocarbons, show altered transcript abundances. Our data show that modulation of mating behavior through reduction of Obp56h is accompanied by altered cuticular hydrocarbon profiles and implicate 5-T as a possible ligand for Obp56h.

Genetic architecture of natural variation in cuticular hydrocarbon composition in Drosophila melanogaster

Insect cuticular hydrocarbons (CHCs) prevent desiccation and serve as chemical signals that mediate social interactions. Drosophila melanogaster CHCs have been studied extensively, but the genetic basis for individual variation in CHC composition is largely unknown. We quantified variation in CHC profiles in the D. melanogaster Genetic Reference Panel (DGRP) and identified novel CHCs. We used principal component (PC) analysis to extract PCs that explain the majority of CHC variation and identified polymorphisms in or near 305 and 173 genes in females and males, respectively, associated with variation in these PCs. In addition, 17 DGRP lines contain the functional Desat2 allele characteristic of African and Caribbean D. melanogaster females (more 5,9-C27:2 and less 7,11-C27:2, female sex pheromone isomers). Disruption of expression of 24 candidate genes affected CHC composition in at least one sex. These genes are associated with fatty acid metabolism and represent mechanistic targets for individual variation in CHC composition.

DOI:http://dx.doi.org/10.7554/eLife.09861.001

The outermost layer of an insect’s body is called the epicuticle and is made of a blend of fat molecules. “Cuticular hydrocarbons” (or CHCs) are the most common fat molecules in the epicuticle, and play an important role in protecting the insect’s body from harsh, dry habitats. CHCs also have other roles in insect behavior. For example, these molecules act as chemical cues when insects search for mates (i.e. pheromones), and they can even contribute to camouflage.

Insects are amongst the most diverse groups of animals on Earth, and different species have different blends of CHC molecules in their epicuticles. Fruit flies are a useful model to understand the genetics of CHC production, including CHCs that act as sex pheromones. Previous research has analyzed the CHCs made by both sexes in several fruit fly strains. However this work was unable to uncover which genes influence how much of a given CHC an individual fly will make.

Dembeck et al. have now looked into CHC production in a collection of 205 different fly strains, all of which have already had their total genetic material sequenced and studied. Comparing these known sequences and looking for associations between genetic differences and particular CHCs uncovered 24 genes that may be involved in CHC manufacture. Only six of the genes had been identified previously. Dembeck et al. found that interfering with the activity of any of the 24 genes had a knock-on effect on many other CHCs present in the flies’ epicuticle.

These 24 genes could to be pieced together in a network that is needed to make and recycle CHCs. The complexity and flexibility of this network can explain in part how insects have been able to build epicuticles for almost every environment. These data set the stage for future work directed towards understanding the evolutionary significance of variation in CHC composition in many fruit fly populations.

DOI:http://dx.doi.org/10.7554/eLife.09861.002

Polymorphisms in early neurodevelopmental genes affect natural variation in alcohol sensitivity in adult drosophila

BACKGROUND

Alcohol abuse and alcoholism are significant public health problems, but the genetic basis for individual variation in alcohol sensitivity remains poorly understood. Drosophila melanogaster presents a powerful model system for dissecting the genetic underpinnings that determine individual variation in alcohol-related phenotypes. We performed genome wide association analyses for alcohol sensitivity using the sequenced, inbred lines of the D. melanogaster Genetic Reference Panel (DGRP) together with extreme QTL mapping in an advanced intercross population derived from sensitive and resistant DGRP lines.

RESULTS

The DGRP harbors substantial genetic variation for alcohol sensitivity and tolerance. We identified 247 candidate genes affecting alcohol sensitivity in the DGRP or the DGRP-derived advanced intercross population, some of which met a Bonferroni-corrected significance threshold, while others occurred among the top candidate genes associated with variation in alcohol sensitivity in multiple analyses. Among these were candidate genes associated with development and function of the nervous system, including several genes in the Dopamine decarboxylase (Ddc) cluster involved in catecholamine synthesis. We found that 58 of these genes formed a genetic interaction network. We verified candidate genes using mutational analysis, targeted gene disruption through RNAi knock-down and transcriptional profiling. Two-thirds of the candidate genes have been implicated in previous Drosophila, mouse and human studies of alcohol-related phenotypes.

CONCLUSIONS

Individual variation in alcohol sensitivity in Drosophila is highly polygenic and in part determined by variation in evolutionarily conserved signaling pathways that are associated with catecholamine neurotransmitter biosynthesis and early development of the nervous system.

The genetic basis for variation in olfactory behavior in Drosophila melanogaster

The genetic underpinnings that contribute to variation in olfactory perception are not fully understood. To explore the genetic basis of variation in olfactory perception, we measured behavioral responses to 14 chemically diverse naturally occurring odorants in 260400 flies from 186 lines of the Drosophila melanogaster Genetic Reference Panel, a population of inbred wild-derived lines with sequenced genomes. We observed variation in olfactory behavior for all odorants. Low to moderate broad-sense heritabilities and the large number of tests for genotype–olfactory phenotype association performed precluded any individual variant from reaching formal significance. However, the top variants (nominal P < 5×10⁻⁵) were highly enriched for genes involved in nervous system development and function, as expected for a behavioral trait. Further, pathway enrichment analyses showed that genes tagged by the top variants included components of networks centered on cyclic guanosine monophosphate and inositol triphosphate signaling, growth factor signaling, Rho signaling, axon guidance, and regulation of neural connectivity. Functional validation with RNAi and mutations showed that 15 out of 17 genes tested indeed affect olfactory behavior. Our results show that in addition to chemoreceptors, variation in olfactory perception depends on polymorphisms that can result in subtle variations in synaptic connectivity within the nervous system.

Genome-wide association analysis of tolerance to methylmercury toxicity in Drosophila implicates myogenic and neuromuscular developmental pathways

Methylmercury (MeHg) is a persistent environmental toxin present in seafood that can compromise the developing nervous system in humans. The effects of MeHg toxicity varies among individuals, despite similar levels of exposure, indicating that genetic differences contribute to MeHg susceptibility. To examine how genetic variation impacts MeHg tolerance, we assessed developmental tolerance to MeHg using the sequenced, inbred lines of the Drosophila melanogaster Genetic Reference Panel (DGRP). We found significant genetic variation in the effects of MeHg on development, measured by eclosion rate, giving a broad sense heritability of 0.86. To investigate the influence of dietary factors, we measured MeHg toxicity with caffeine supplementation in the DGRP lines. We found that caffeine counteracts the deleterious effects of MeHg in the majority of lines, and there is significant genetic variance in the magnitude of this effect, with a broad sense heritability of 0.80. We performed genome-wide association (GWA) analysis for both traits, and identified candidate genes that fall into several gene ontology categories, with enrichment for genes involved in muscle and neuromuscular development. Overexpression of glutamate-cysteine ligase, a MeHg protective enzyme, in a muscle-specific manner leads to a robust rescue of eclosion of flies reared on MeHg food. Conversely, mutations in kirre, a pivotal myogenic gene identified in our GWA analyses, modulate tolerance to MeHg during development in accordance with kirre expression levels. Finally, we observe disruptions of indirect flight muscle morphogenesis in MeHg-exposed pupae. Since the pathways for muscle development are evolutionarily conserved, it is likely that the effects of MeHg observed in Drosophila can be generalized across phyla, implicating muscle as an additional hitherto unrecognized target for MeHg toxicity. Furthermore, our observations that caffeine can ameliorate the toxic effects of MeHg show that nutritional factors and dietary manipulations may offer protection against the deleterious effects of MeHg exposure.

Transcriptional and epigenetic responses to mating and aging in Drosophila melanogaster

BACKGROUND

Phenotypic plasticity allows organisms to respond rapidly to changing environmental circumstances, and understanding its genomic basis can yield insights regarding the underlying genes and genetic networks affecting complex phenotypes. Female Drosophila melanogaster undergo dramatic physiological changes mediated by seminal fluid components transferred upon mating, including decreased longevity. Their physiological and behavioral effects have been well characterized, but little is known about resulting changes in regulation of gene expression or the extent to which mating-induced changes in gene expression are the same as those occurring during aging.

RESULTS

We assessed genome-wide mRNA, microRNA, and three common histone modifications implicated in gene activation for young and aged virgin and mated female D. melanogaster in a factorial design. We identified phenotypically plastic transcripts and epigenetic modifications associated with mating and aging. We used these data to derive phenotypically plastic regulatory networks associated with mating of young flies, and aging of virgin and mated flies. Many of the mRNAs, microRNAs and epigenetic modifications associated with mating of young flies also occur with age in virgin flies, which may reflect mating-induced accelerated aging. We functionally tested the plastic regulatory networks by overexpressing environmentally sensitive microRNAs. Overexpression resulted in altered expression of ~70% of candidate target genes, and in all cases affected oviposition.

CONCLUSIONS

Our results implicate microRNAs as mediators of phenotypic plasticity associated with mating and provide a comprehensive documentation of the genomic and epigenomic changes that accompany mating- and aging-induced physiological changes in female D. melanogaster.

Natural variation in genome architecture among 205 Drosophila melanogaster Genetic Reference Panel lines

The Drosophila melanogaster Genetic Reference Panel (DGRP) is a community resource of 205 sequenced inbred lines, derived to improve our understanding of the effects of naturally occurring genetic variation on molecular and organismal phenotypes. We used an integrated genotyping strategy to identify 4,853,802 single nucleotide polymorphisms (SNPs) and 1,296,080 non-SNP variants. Our molecular population genomic analyses show higher deletion than insertion mutation rates and stronger purifying selection on deletions. Weaker selection on insertions than deletions is consistent with our observed distribution of genome size determined by flow cytometry, which is skewed toward larger genomes. Insertion/deletion and single nucleotide polymorphisms are positively correlated with each other and with local recombination, suggesting that their nonrandom distributions are due to hitchhiking and background selection. Our cytogenetic analysis identified 16 polymorphic inversions in the DGRP. Common inverted and standard karyotypes are genetically divergent and account for most of the variation in relatedness among the DGRP lines. Intriguingly, variation in genome size and many quantitative traits are significantly associated with inversions. Approximately 50% of the DGRP lines are infected with Wolbachia, and four lines have germline insertions of Wolbachia sequences, but effects of Wolbachia infection on quantitative traits are rarely significant. The DGRP complements ongoing efforts to functionally annotate the Drosophila genome. Indeed, 15% of all D. melanogaster genes segregate for potentially damaged proteins in the DGRP, and genome-wide analyses of quantitative traits identify novel candidate genes. The DGRP lines, sequence data, genotypes, quality scores, phenotypes, and analysis and visualization tools are publicly available.

Olfactomedin proteins: central players in development and disease

Olfactomedin proteins are characterized by a conserved domain of \texorpdfstring~\textasciitilde250 amino acids corresponding to the olfactomedin archetype first discovered in olfactory neuroepithelium. They arose early in evolution and occur throughout the animal kingdom. In mice and humans olfactomedin proteins comprise a diverse array of glycoproteins, many of which are critical for early development and functional organization of the nervous system as well as hematopoiesis. Olfactomedin domains appear to facilitate protein-protein interactions, intercellular interactions, and cell adhesion. Several members of the family have been implicated in various common diseases, notably myocilin in glaucoma and OLFM4 in cancer. This review highlights this important, hitherto understudied family of proteins.

Modulation of feeding behavior by odorant-binding proteins in Drosophila melanogaster

Nutrient intake and avoidance of toxins are essential for survival and controlled by attractive and aversive feeding responses. Drosophila melanogaster presents one of the best characterized systems for studies on chemosensation, which is mediated by multigene families of chemoreceptors, including olfactory receptors, gustatory receptors, and odorant-binding proteins (OBPs). Although the response profiles of gustatory receptors have been well studied, the contribution of OBPs to food intake is largely unknown. As most aversive (“bitter”) tastants are hydrophobic, we hypothesized that OBPs may fulfill an essential function in transporting bitter tastants to gustatory receptors to modulate feeding behavior. Here, we used 16 RNAi lines that inhibit expression of individual target Obp genes and show that OBPs modulate sucrose intake in response to a panel of nine bitter compounds. Similar to their function in olfaction, OBPs appear to interact with bitter compounds in a combinatorial and sex-dependent manner. RNAi-mediated reduction in expression of individual Obp genes resulted either in enhanced or reduced intake of sucrose in the presence of bitter compounds, consistent with roles for OBPs in transporting tastants to bitter taste receptors, sequestering them to limit their access to these receptors, or interacting directly with gustatory neurons that respond to sucrose.

A molecular mechanism for glaucoma: endoplasmic reticulum stress and the unfolded protein response

Primary open angle glaucoma (POAG) is a common late-onset neurodegenerative disease. Ocular hypertension represents a major risk factor, but POAG etiology remains poorly understood. Some cases of early-onset congenital glaucoma and adult POAG are linked to mutations in myocilin, a secreted protein of poorly defined function. Transgenic overexpression of myocilin in Drosophila and experiments in mice and human populations implicate the unfolded protein response (UPR) in the pathogenesis of glaucoma. We postulate that compromised ability of the UPR to eliminate misfolded mutant or damaged proteins, including myocilin, causes endoplasmic reticulum stress, resulting in functional impairment of trabecular meshwork cells that regulate intraocular pressure. This mechanism of POAG is reminiscent of other age-dependent neurodegenerative diseases that involve accumulation of protein aggregates.

Epistasis dominates the genetic architecture of Drosophila quantitative traits

Epistasis-nonlinear genetic interactions between polymorphic loci-is the genetic basis of canalization and speciation, and epistatic interactions can be used to infer genetic networks affecting quantitative traits. However, the role that epistasis plays in the genetic architecture of quantitative traits is controversial. Here, we compared the genetic architecture of three Drosophila life history traits in the sequenced inbred lines of the Drosophila melanogaster Genetic Reference Panel (DGRP) and a large outbred, advanced intercross population derived from 40 DGRP lines (Flyland). We assessed allele frequency changes between pools of individuals at the extremes of the distribution for each trait in the Flyland population by deep DNA sequencing. The genetic architecture of all traits was highly polygenic in both analyses. Surprisingly, none of the SNPs associated with the traits in Flyland replicated in the DGRP and vice versa. However, the majority of these SNPs participated in at least one epistatic interaction in the DGRP. Despite apparent additive effects at largely distinct loci in the two populations, the epistatic interactions perturbed common, biologically plausible, and highly connected genetic networks. Our analysis underscores the importance of epistasis as a principal factor that determines variation for quantitative traits and provides a means to uncover genetic networks affecting these traits. Knowledge of epistatic networks will contribute to our understanding of the genetic basis of evolutionarily and clinically important traits and enhance predictive ability at an individualized level in medicine and agriculture.

Genetics of aggression

Aggression mediates competition for food, mating partners, and habitats and, among social animals, establishes stable dominance hierarchies. In humans, abnormal aggression is a hallmark of neuropsychiatric disorders and can be elicited by environmental factors acting on an underlying genetic susceptibility. Identifying the genetic architecture that predisposes to aggressive behavior in people is challenging because of difficulties in quantifying the phenotype, genetic heterogeneity, and uncontrolled environmental conditions. Studies on mice have identified single-gene mutations that result in hyperaggression, contingent on genetic background. These studies can be complemented by systems genetics approaches in Drosophila melanogaster, in which mutational analyses together with genome-wide transcript analyses, artificial selection studies, and genome-wide analysis of epistasis have revealed that a large segment of the genome contributes to the manifestation of aggressive behavior with widespread epistatic interactions. Comparative genomic analyses based on the principle of evolutionary conservation are needed to enable a complete dissection of the neurogenetic underpinnings of this universal fitness trait.

Genome-wide association for sensitivity to chronic oxidative stress in Drosophila melanogaster

Reactive oxygen species (ROS) are a common byproduct of mitochondrial energy metabolism, and can also be induced by exogenous sources, including UV light, radiation, and environmental toxins. ROS generation is essential for maintaining homeostasis by triggering cellular signaling pathways and host defense mechanisms. However, an imbalance of ROS induces oxidative stress and cellular death and is associated with human disease, including age-related locomotor impairment. To identify genes affecting sensitivity and resistance to ROS-induced locomotor decline, we assessed locomotion of aged flies of the sequenced, wild-derived lines from the Drosophila melanogaster Genetics Reference Panel on standard medium and following chronic exposure to medium supplemented with 3 mM menadione sodium bisulfite (MSB). We found substantial genetic variation in sensitivity to oxidative stress with respect to locomotor phenotypes. We performed genome-wide association analyses to identify candidate genes associated with variation in sensitivity to ROS-induced decline in locomotor performance, and confirmed the effects for 13 of 16 mutations tested in these candidate genes. Candidate genes associated with variation in sensitivity to MSB-induced oxidative stress form networks of genes involved in neural development, immunity, and signal transduction. Many of these genes have human orthologs, highlighting the utility of genome-wide association in Drosophila for studying complex human disease.

Genome-wide association analysis of oxidative stress resistance in Drosophila melanogaster

BACKGROUND

Aerobic organisms are susceptible to damage by reactive oxygen species. Oxidative stress resistance is a quantitative trait with population variation attributable to the interplay between genetic and environmental factors. Drosophila melanogaster provides an ideal system to study the genetics of variation for resistance to oxidative stress.

METHODS AND FINDINGS

We used 167 wild-derived inbred lines of the Drosophila Genetic Reference Panel for a genome-wide association study of acute oxidative stress resistance to two oxidizing agents, paraquat and menadione sodium bisulfite. We found significant genetic variation for both stressors. Single nucleotide polymorphisms (SNPs) associated with variation in oxidative stress resistance were often sex-specific and agent-dependent, with a small subset common for both sexes or treatments. Associated SNPs had moderately large effects, with an inverse relationship between effect size and allele frequency. Linear models with up to 12 SNPs explained 67-79% and 56-66% of the phenotypic variance for resistance to paraquat and menadione sodium bisulfite, respectively. Many genes implicated were novel with no known role in oxidative stress resistance. Bioinformatics analyses revealed a cellular network comprising DNA metabolism and neuronal development, consistent with targets of oxidative stress-inducing agents. We confirmed associations of seven candidate genes associated with natural variation in oxidative stress resistance through mutational analysis.

CONCLUSIONS

We identified novel candidate genes associated with variation in resistance to oxidative stress that have context-dependent effects. These results form the basis for future translational studies to identify oxidative stress susceptibility/resistance genes that are evolutionary conserved and might play a role in human disease.

Phenotypic plasticity of the Drosophila transcriptome

Phenotypic plasticity is the ability of a single genotype to produce different phenotypes in response to changing environments. We assessed variation in genome-wide gene expression and four fitness-related phenotypes of an outbred Drosophila melanogaster population under 20 different physiological, social, nutritional, chemical, and physical environments; and we compared the phenotypically plastic transcripts to genetically variable transcripts in a single environment. The environmentally sensitive transcriptome consists of two transcript categories, which comprise ∼15% of expressed transcripts. Class I transcripts are genetically variable and associated with detoxification, metabolism, proteolysis, heat shock proteins, and transcriptional regulation. Class II transcripts have low genetic variance and show sexually dimorphic expression enriched for reproductive functions. Clustering analysis of Class I transcripts reveals a fragmented modular organization and distinct environmentally responsive transcriptional signatures for the four fitness-related traits. Our analysis suggests that a restricted environmentally responsive segment of the transcriptome preserves the balance between phenotypic plasticity and environmental canalization.

Unlike Mendelian traits, where the genotype allows a direct prediction of the phenotype, predicting phenotypic values is not straightforward for complex traits, which arise from multiple segregating genes and their interactions with the environment. Here, a single genotype can often express different phenotypes in different environments. Such phenotypic plasticity is the counterpoint to “environmental canalization,” whereby genotypes produce the same phenotype in different environments. Whereas phenotypic plasticity allows organisms to respond rapidly to changing environments, environmental canalization buffers phenotypes against environmental perturbations. The balance between plasticity and robustness is crucial for optimal fitness, but the genetic basis for phenotypic plasticity is poorly defined. Here, we present the most comprehensive analysis to date of variation in genome-wide gene expression of an outbred Drosophila melanogaster population under 20 different environments. We find that a restricted environmentally responsive segment of the transcriptome (∼15%) preserves the balance between phenotypic plasticity and environmental canalization. Environmentally plastic transcripts can be grouped into two categories. Class I transcripts are genetically variable and associated with detoxification, metabolism, proteolysis, heat shock proteins, and transcriptional regulation. Class II transcripts have low genetic variance and show sexually dimorphic expression enriched for reproductive functions. Despite low genetic variance these transcripts evolve rapidly.

Functional dissection of Odorant binding protein genes in Drosophila melanogaster

Most organisms rely on olfaction for survival and reproduction. The olfactory system of Drosophila melanogaster is one of the best characterized chemosensory systems and serves as a prototype for understanding insect olfaction. Olfaction in Drosophila is mediated by multigene families of odorant receptors and odorant binding proteins (OBPs). Although molecular response profiles of odorant receptors have been well documented, the contributions of OBPs to olfactory behavior remain largely unknown. Here, we used RNAi-mediated suppression of Obp gene expression and measurements of behavioral responses to 16 ecologically relevant odorants to systematically dissect the functions of 17 OBPs. We quantified the effectiveness of RNAi-mediated suppression by quantitative real-time polymerase chain reaction and used a proteomic liquid chromatography and tandem mass spectrometry procedure to show target-specific suppression of OBPs expressed in the antennae. Flies in which expression of a specific OBP is suppressed often show altered behavioral responses to more than one, but not all, odorants, in a sex-dependent manner. Similarly, responses to a specific odorant are frequently affected by suppression of expression of multiple, but not all, OBPs. These results show that OBPs are essential for mediating olfactory behavioral responses and suggest that OBP-dependent odorant recognition is combinatorial.

Tuning the chemosensory window: a fly's perspective

Accurate perception of chemical signals from the environment is critical for the fitness of most animals. Drosophila melanogaster experiences its chemical environment through families of chemoreceptors that include olfactory receptors, gustatory receptors and odorant binding proteins. Its chemical environment, however, changes during its life cycle and the interpretation of chemical signals is dependent on dynamic social and physical surroundings. Phenotypic plasticity of gene expression of the chemoreceptor repertoire allows flies to adjust the chemosensory window through which they "view" their world and to modify the ensemble of expressed chemoreceptor proteins in line with their developmental and physiological state and according to their needs to locate food and oviposition sites under different social and physical environmental conditions. Furthermore, males and females differ in their expression profiles of chemoreceptor genes. Thus, each sex experiences its chemical environment via combinatorial activation of distinct chemoreceptor ensembles. The remarkable phenotypic plasticity of the chemoreceptor repertoire raises several fundamental questions. What are the mechanisms that translate environmental cues into regulation of chemoreceptor gene expression? How are gustatory and olfactory cues integrated perceptually? What is the relationship between ensembles of odorant binding proteins and odorant receptors? And, what is the significance of co-regulated chemoreceptor transcriptional networks?

The soluble proteome of the Drosophila antenna

The olfactory system of Drosophila melanogaster is one of the best characterized chemosensory systems. Identification of proteins contained in the third antennal segment, the main olfactory organ, has previously relied primarily on immunohistochemistry, and although such studies and in situ hybridization studies are informative, they focus generally on one or few gene products at a time, and quantification is difficult. In addition, purification of native proteins from the antenna is challenging because it is small and encased in a hard cuticle. Here, we describe a simple method for the large-scale detection of soluble proteins from the Drosophila antenna by chromatographic separation of tryptic peptides followed by tandem mass spectrometry with femtomole detection sensitivities. Examination of the identities of these proteins indicates that they originate both from the extracellular perilymph and from the cytoplasm of disrupted cells. We identified enzymes involved with intermediary metabolism, proteins associated with regulation of gene expression, nucleic acid metabolism and protein metabolism, proteins associated with microtubular transport, 8 odorant-binding proteins, protective enzymes associated with antibacterial defense and defense against oxidative damage, cuticular proteins, and proteins of unknown function, which represented about one-third of all soluble proteins. The procedure described here opens the way for precise quantification of any target protein in the Drosophila antenna and should be readily applicable to antennae from other insects.

Plasticity of the chemoreceptor repertoire in Drosophila melanogaster

For most organisms, chemosensation is critical for survival and is mediated by large families of chemoreceptor proteins, whose expression must be tuned appropriately to changes in the chemical environment. We asked whether expression of chemoreceptor genes that are clustered in the genome would be regulated independently; whether expression of certain chemoreceptor genes would be especially sensitive to environmental changes; whether groups of chemoreceptor genes undergo coordinated rexpression; and how plastic the expression of chemoreceptor genes is with regard to sex, development, reproductive state, and social context. To answer these questions we used Drosophila melanogaster, because its chemosensory systems are well characterized and both the genotype and environment can be controlled precisely. Using customized cDNA microarrays, we showed that chemoreceptor genes that are clustered in the genome undergo independent transcriptional regulation at different developmental stages and between sexes. Expression of distinct subgroups of chemoreceptor genes is sensitive to reproductive state and social interactions. Furthermore, exposure of flies only to odor of the opposite sex results in altered transcript abundance of chemoreceptor genes. These genes are distinct from those that show transcriptional plasticity when flies are allowed physical contact with same or opposite sex members. We analyzed covariance in transcript abundance of chemosensory genes across all environmental conditions and found that they segregated into 20 relatively small, biologically relevant modules of highly correlated transcripts. This finely pixilated modular organization of the chemosensory subgenome enables fine tuning of the expression of the chemoreceptor repertoire in response to ecologically relevant environmental and physiological conditions.

Rapid adaptation and phenotypic plasticity to the chemical environment are essential prerequisites for survival; and, consequently, large families of genes that mediate the recognition of olfactory and gustatory cues have evolved. We asked how flexible the expression of these genes is in the face of rapidly changing conditions encountered during an individual's lifetime. We used the fruit fly, Drosophila melanogaster, to address this question, since both the genetic composition and environmental rearing conditions can be controlled precisely in this experimentally amenable model organism. By measuring expression levels of all chemosensory genes simultaneously, we identified genes that show altered expression at different developmental stages, during aging, in males and females, following mating, and in different social conditions. We asked whether chemosensory genes are regulated independently or whether their regulation is structured. We found that chemosensory genes that are located in close proximity to one another on the chromosome are often regulated independently. However, statistical analysis showed that groups of chemosensory genes are coordinately expressed in response to a range of environmental conditions, revealing an underlying modular organization of the phenotypic plasticity of the chemosensory receptor repertoire.

Systems genetics of complex traits in Drosophila melanogaster

Determining the genetic architecture of complex traits is challenging because phenotypic variation arises from interactions between multiple, environmentally sensitive alleles. We quantified genome-wide transcript abundance and phenotypes for six ecologically relevant traits in D. melanogaster wild-derived inbred lines. We observed 10,096 genetically variable transcripts and high heritabilities for all organismal phenotypes. The transcriptome is highly genetically inter-correlated, forming 241 transcriptional modules. Modules are enriched for transcripts in common pathways, gene ontology categories, tissue-specific expression, and transcription factor binding sites. The high transcriptional connectivity allows us to infer genetic networks and the function of predicted genes based on annotations of other genes in the network. Regressions of organismal phenotypes on transcript abundance implicate several hundred candidate genes that form modules of biologically meaningful correlated transcripts affecting each phenotype. Overlapping transcripts in modules associated with different traits provides insight into the molecular basis of pleiotropy between complex traits.

Overexpression of myocilin in the Drosophila eye activates the unfolded protein response: implications for glaucoma

Background

Glaucoma is the world's second leading cause of bilateral blindness with progressive loss of vision due to retinal ganglion cell death. Myocilin has been associated with congenital glaucoma and 2–4% of primary open angle glaucoma (POAG) cases, but the pathogenic mechanisms remain largely unknown. Among several hypotheses, activation of the unfolded protein response (UPR) has emerged as a possible disease mechanism.

Methodology / Principal Findings

We used a transgenic Drosophila model to analyze whole-genome transcriptional profiles in flies that express human wild-type or mutant MYOC in their eyes. The transgenic flies display ocular fluid discharge, reflecting ocular hypertension, and a progressive decline in their behavioral responses to light. Transcriptional analysis shows that genes associated with the UPR, ubiquitination, and proteolysis, as well as metabolism of reactive oxygen species and photoreceptor activity undergo altered transcriptional regulation. Following up on the results from these transcriptional analyses, we used immunoblots to demonstrate the formation of MYOC aggregates and showed that the formation of such aggregates leads to induction of the UPR, as evident from activation of the fluorescent UPR marker, xbp1-EGFP.

Conclusions / Significance

Our results show that aggregation of MYOC in the endoplasmic reticulum activates the UPR, an evolutionarily conserved stress pathway that culminates in apoptosis. We infer from the Drosophila model that MYOC-associated ocular hypertension in the human eye may result from aggregation of MYOC and induction of the UPR in trabecular meshwork cells. This process could occur at a late age with wild-type MYOC, but might be accelerated by MYOC mutants to account for juvenile onset glaucoma.

Phenotypic and transcriptional response to selection for alcohol sensitivity in Drosophila melanogaster

Gene-expression profiling combined with selection for genetically divergent Drosophila lines either highly sensitive or resistant to ethanol exposure has been used to identify candidate genes that affect alcohol sensitivity, including 23 novel genes that have human orthologs.

Background

Alcoholism is a complex disorder determined by interactions between genetic and environmental risk factors. Drosophila represents a powerful model system to dissect the genetic architecture of alcohol sensitivity, as large numbers of flies can readily be reared in defined genetic backgrounds and under controlled environmental conditions. Furthermore, flies exposed to ethanol undergo physiological and behavioral changes that resemble human alcohol intoxication, including loss of postural control, sedation, and development of tolerance.

Results

We performed artificial selection for alcohol sensitivity for 35 generations and created duplicate selection lines that are either highly sensitive or resistant to ethanol exposure along with unselected control lines. We used whole genome expression analysis to identify 1,678 probe sets with different expression levels between the divergent lines, pooled across replicates, at a false discovery rate of q < 0.001. We assessed to what extent genes with altered transcriptional regulation might be causally associated with ethanol sensitivity by measuring alcohol sensitivity of 37 co-isogenic P-element insertional mutations in 35 candidate genes, and found that 32 of these mutants differed in sensitivity to ethanol exposure from their co-isogenic controls. Furthermore, 23 of these novel genes have human orthologues.

Conclusion

Combining whole genome expression profiling with selection for genetically divergent lines is an effective approach for identifying candidate genes that affect complex traits, such as alcohol sensitivity. Because of evolutionary conservation of function, it is likely that human orthologues of genes affecting alcohol sensitivity in Drosophila may contribute to alcohol-associated phenotypes in humans.