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

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

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.

Selection for long and short sleep duration in Drosophila melanogaster reveals the complex genetic network underlying natural variation in sleep.

Why do some individuals need more sleep than others? Forward mutagenesis screens in flies using engineered mutations have established a clear genetic component to sleep duration, revealing mutants that convey very long or short sleep. Whether such extreme long or short sleep could exist in natural populations was unknown. We applied artificial selection for high and low night sleep duration to an outbred population of Drosophila melanogaster for 13 generations. At the end of the selection procedure, night sleep duration diverged by 9.97 hours in the long and short sleeper populations, and 24-hour sleep was reduced to 3.3 hours in the short sleepers. Neither long nor short sleeper lifespan differed appreciably from controls, suggesting little physiological consequences to being an extreme long or short sleeper. Whole genome sequence data from seven generations of selection revealed several hundred thousand changes in allele frequencies at polymorphic loci across the genome. Combining the data from long and short sleeper populations across generations in a logistic regression implicated 126 polymorphisms in 80 candidate genes, and we confirmed three of these genes and a larger genomic region with mutant and chromosomal deficiency tests, respectively. Many of these genes could be connected in a single network based on previously known physical and genetic interactions. Candidate genes have known roles in several classic, highly conserved developmental and signaling pathways-EGFR, Wnt, Hippo, and MAPK. The involvement of highly pleiotropic pathway genes suggests that sleep duration in natural populations can be influenced by a wide variety of biological processes, which may be why the purpose of sleep has been so elusive.

Genome-Wide Association Study of Circadian Behavior in Drosophila melanogaster.

Circadian rhythms influence physiological processes from sleep-wake cycles to body temperature and are controlled by highly conserved cycling molecules. Although the mechanistic basis of the circadian clock has been known for decades, the extent to which circadian rhythms vary in nature and the underlying genetic basis for that variation is not well understood. We measured circadian period (Æ®) and rhythmicity index in the Drosophila Genetic Reference Panel (DGRP) and observed extensive genetic variation in both. Seven DGRP lines had sexually dimorphic arrhythmicity and one line had an exceptionally long Æ®. Genome-wide analyses identified 584 polymorphisms in 268 genes. We observed differences among transcripts for nine genes predicted to interact among themselves and canonical clock genes in the long period line and a control. Mutations/RNAi knockdown targeting these genes also affected circadian behavior. Our observations reveal that complex genetic interactions influence high levels of variation in circadian phenotypes.

Comparison of normalization and differential expression analyses using RNA-Seq data from 726 individual Drosophila melanogaster.

BACKGROUND: A generally accepted approach to the analysis of RNA-Seq read count data does not yet exist. We sequenced the mRNA of 726 individuals from the Drosophila Genetic Reference Panel in order to quantify differences in gene expression among single flies. One of our experimental goals was to identify the optimal analysis approach for the detection of differential gene expression among the factors we varied in the experiment: genotype, environment, sex, and their interactions. Here we evaluate three different filtering strategies, eight normalization methods, and two statistical approaches using our data set. We assessed differential gene expression among factors and performed a statistical power analysis using the eight biological replicates per genotype, environment, and sex in our data set. RESULTS: We found that the most critical considerations for the analysis of RNA-Seq read count data were the normalization method, underlying data distribution assumption, and numbers of biological replicates, an observation consistent with previous RNA-Seq and microarray analysis comparisons. Some common normalization methods, such as Total Count, Quantile, and RPKM normalization, did not align the data across samples. Furthermore, analyses using the Median, Quantile, and Trimmed Mean of M-values normalization methods were sensitive to the removal of low-expressed genes from the data set. Although it is robust in many types of analysis, the normal data distribution assumption produced results vastly different than the negative binomial distribution. In addition, at least three biological replicates per condition were required in order to have sufficient statistical power to detect expression differences among the three-way interaction of genotype, environment, and sex. CONCLUSIONS: The best analysis approach to our data was to normalize the read counts using the DESeq method and apply a generalized linear model assuming a negative binomial distribution using either edgeR or DESeq software. Genes having very low read counts were removed after normalizing the data and fitting it to the negative binomial distribution. We describe the results of this evaluation and include recommended analysis strategies for RNA-Seq read count data.

Knockdown expression of Syndecan in the fat body impacts nutrient metabolism and the organismal response to environmental stresses in Drosophila melanogaster.

The heparan sulfate proteoglycan syndecans are transmembrane proteins involved in multiple physiological processes, including cell-matrix adhesion and inflammation. Recent evidence from model systems and humans suggest that syndecans have a role in energy balance and nutrient metabolism regulation. However, much remains to be learned about the mechanisms through which syndecans influence these phenotypes. Previously, we reported that Drosophila melanogaster Syndecan (Sdc) mutants had reduced metabolic activity compared to controls. Here, we knocked down endogenous Sdc expression in the fat body (the functional equivalent of mammalian adipose tissue and liver) to investigate whether the effects on metabolism originate from this tissue. We found that knocking down Sdc in the fat body leads to flies with higher levels of glycogen and fat and that survive longer during starvation, likely due to their extra energy reserves and an increase in gluconeogenesis. However, compared to control flies, they are also more sensitive to environmental stresses (e.g. bacterial infection and cold) and have reduced metabolic activity under normal feeding conditions. Under the same conditions, fat-body Sdc reduction enhances expression of genes involved in glyceroneogenesis and gluconeogenesis and induces a drastic decrease in phosphorylation levels of AKT and extracellular signal regulated kinase 1/2 (ERK1/2). Altogether, these findings strongly suggest that Drosophila fat body Sdc is involved in a mechanism that shifts resources to different physiological functions according to nutritional status.

Genome-wide association in Drosophila identifies a role for Piezo and Proc-R in sleep latency.

Sleep latency, the amount of time that it takes an individual to fall asleep, is a key indicator of sleep need. Sleep latency varies considerably both among and within species and is heritable, but lacks a comprehensive description of its underlying genetic network. Here we conduct a genome-wide association study of sleep latency. Using previously collected sleep and activity data on a wild-derived population of flies, we calculate sleep latency, confirming significant, heritable genetic variation for this complex trait. We identify 520 polymorphisms in 248 genes contributing to variability in sleep latency. Tests of mutations in 23 candidate genes and additional putative pan-neuronal knockdown of 9 of them implicated CG44153, Piezo, Proc-R and Rbp6 in sleep latency. Two large-effect mutations in the genes Proc-R and Piezo were further confirmed via genetic rescue. This work greatly enhances our understanding of the genetic factors that influence variation in sleep latency.

Genome-wide association study of sleep in Drosophila melanogaster.

BACKGROUND: Sleep is a highly conserved behavior, yet its duration and pattern vary extensively among species and between individuals within species. The genetic basis of natural variation in sleep remains unknown. RESULTS: We used the Drosophila Genetic Reference Panel (DGRP) to perform a genome-wide association (GWA) study of sleep in D. melanogaster. We identified candidate single nucleotide polymorphisms (SNPs) associated with differences in the mean as well as the environmental sensitivity of sleep traits; these SNPs typically had sex-specific or sex-biased effects, and were generally located in non-coding regions. The majority of SNPs (80.3%) affecting sleep were at low frequency and had moderately large effects. Additive models incorporating multiple SNPs explained as much as 55% of the genetic variance for sleep in males and females. Many of these loci are known to interact physically and/or genetically, enabling us to place them in candidate genetic networks. We confirmed the role of seven novel loci on sleep using insertional mutagenesis and RNA interference. CONCLUSIONS: We identified many SNPs in novel loci that are potentially associated with natural variation in sleep, as well as SNPs within genes previously known to affect Drosophila sleep. Several of the candidate genes have human homologues that were identified in studies of human sleep, suggesting that genes affecting variation in sleep are conserved across species. Our discovery of genetic variants that influence environmental sensitivity to sleep may have a wider application to all GWA studies, because individuals with highly plastic genotypes will not have consistent phenotypes.

A conserved role for frizzled in sleep architecture.

Previous studies of natural variants in Drosophila melanogaster implicated the Wnt signaling receptor frizzled in sleep. Given that the Wnt signaling pathway is highly conserved across species, we hypothesized that frizzled class receptor 1 (Fzd1), the murine homolog of frizzled, would also have a role in sleep. Using a CRISPR transgenic approach, we removed most of the Fzd1 coding region from C57BL/6N mice. We used a video assay to measure sleep characteristics in Fzd1-deficient mice. As Wnt signaling is known to affect visuospatial memory, we also examined the impact of the deletion on learning and memory using the novel object recognition (NOR) paradigm. Fzd1-deficient mice had altered sleep compared to littermate controls. The mice did not respond differently to the NOR paradigm compared to controls but did display anxiety-like behavior. Our strategy demonstrates that the study of natural variation in Drosophila sleep translates into candidate genes for sleep in vertebrate species such as the mouse.

Whole Animal Feeding FLat (WAFFL): a complete and comprehensive validation of a novel, high-throughput fly experimentation system.

Non-mammalian model organisms have been essential for our understanding of the mechanisms that control development, disease, and physiology, but they are underutilized in pharmacological and toxicological phenotypic screening assays due to their low throughput in comparison with cell-based screens. To increase the utility of using Drosophila melanogaster in screening, we designed the Whole Animal Feeding FLat (WAFFL), a novel, flexible, and complete system for feeding, monitoring, and assaying flies in a high-throughput format. Our 3D printed system is compatible with inexpensive and readily available, commercial 96-well plate consumables and equipment. Experimenters can change the diet at will during the experiment and video record for behavior analysis, enabling precise dosing, measurement of feeding, and analysis of behavior in a 96-well plate format.

Nuclear genomic control of naturally occurring variation in mitochondrial function in Drosophila melanogaster.

BACKGROUND: Mitochondria are organelles found in nearly all eukaryotic cells that play a crucial role in cellular survival and function. Mitochondrial function is under the control of nuclear and mitochondrial genomes. While the latter has been the focus of most genetic research, we remain largely ignorant about the nuclear-encoded genomic control of inter-individual variability in mitochondrial function. Here, we used Drosophila melanogaster as our model organism to address this question. RESULTS: We quantified mitochondrial state 3 and state 4 respiration rates and P:O ratio in mitochondria isolated from the thoraces of 40 sequenced inbred lines of the Drosophila Genetic Reference Panel. We found significant within-population genetic variability for all mitochondrial traits. Hence, we performed genome-wide association mapping and identified 141 single nucleotide polymorphisms (SNPs) associated with differences in mitochondrial respiration and efficiency (P ≤1 × 10-5). Gene-centered regression models showed that 2-3 SNPs can explain 31, 13, and 18% of the phenotypic variation in state 3, state 4, and P:O ratio, respectively. Most of the genes tagged by the SNPs are involved in organ development, second messenger-mediated signaling pathways, and cytoskeleton remodeling. One of these genes, sallimus (sls), encodes a component of the muscle sarcomere. We confirmed the direct effect of sls on mitochondrial respiration using two viable mutants and their coisogenic wild-type strain. Furthermore, correlation network analysis revealed that sls functions as a transcriptional hub in a co-regulated module associated with mitochondrial respiration and is connected to CG7834, which is predicted to encode a protein with mitochondrial electron transfer flavoprotein activity. This latter finding was also verified in the sls mutants. CONCLUSIONS: Our results provide novel insights into the genetic factors regulating natural variation in mitochondrial function in D. melanogaster. The integrative genomic approach used in our study allowed us to identify sls as a novel hub gene responsible for the regulation of mitochondrial respiration in muscle sarcomere and to provide evidence that sls might act via the electron transfer flavoprotein/ubiquinone oxidoreductase complex.

Understanding the neurogenetics of sleep: progress from Drosophila.

Most behaviors manifest themselves through interactions with environments. Sleep, however, is characterized by immobility and reduced responsiveness. Although nearly all animals sleep, the purpose of sleep remains an enduring puzzle. Drosophila melanogaster exhibits all the behavioral characteristics of mammalian sleep, enabling the use of powerful genetic approaches to dissect conserved fundamental neurogenetic aspects of sleep. Drosophila studies over the past four years have identified novel genes and pathways modulating sleep, such as Shaker and sleepless, and candidate brain regions known to function in circadian regulation and learning and memory. Advances in systems genetics coupled with the ability to target specific brain regions enable the characterization of transcriptional networks and neural circuits contributing to phenotypic variation in sleep.

Extensive epistasis for olfactory behaviour, sleep and waking activity in Drosophila melanogaster.

Epistasis is an important feature of the genetic architecture of quantitative traits, but the dynamics of epistatic interactions in natural populations and the relationship between epistasis and pleiotropy remain poorly understood. Here, we studied the effects of epistatic modifiers that segregate in a wild-derived Drosophila melanogaster population on the mutational effects of P-element insertions in Semaphorin-5C (Sema-5c) and Calreticulin (Crc), pleiotropic genes that affect olfactory behaviour and startle behaviour and, in the case of Crc, sleep phenotypes. We introduced Canton-S B (CSB) third chromosomes with or without a P-element insertion at the Crc or Sema-5c locus in multiple wild-derived inbred lines of the Drosophila melanogaster Genetic Reference Panel (DGRP) and assessed the effects of epistasis on the olfactory response to benzaldehyde and, for Crc, also on sleep. In each case, we found substantial epistasis and significant variation in the magnitude of epistasis. The predominant direction of epistatic effects was to suppress the mutant phenotype. These observations support a previous study on startle behaviour using the same D. melanogaster chromosome substitution lines, which concluded that suppressing epistasis may buffer the effects of new mutations. However, epistatic effects are not correlated among the different phenotypes. Thus, suppressing epistasis appears to be a pervasive general feature of natural populations to protect against the effects of new mutations, but different epistatic interactions modulate different phenotypes affected by mutations at the same pleiotropic gene.

What have we learned about sleep from selective breeding strategies?

Selective breeding is a classic technique that enables an experimenter to modify a heritable target trait as desired. Direct selective breeding for extreme sleep and circadian phenotypes in flies successfully alters these behaviors, and sleep and circadian perturbations emerge as correlated responses to selection for other traits in mice, rats, and dogs. The application of sequencing technologies to the process of selective breeding identifies the genetic network impacting the selected trait in a holistic way. Breeding techniques preserve the extreme phenotypes generated during selective breeding, generating community resources for further functional testing. Selective breeding is thus a unique strategy that can explore the phenotypic limits of sleep and circadian behavior, discover correlated responses of traits having shared genetic architecture with the target trait, identify naturally-occurring genomic variants and gene expression changes that affect trait variability, and pinpoint genes with conserved roles.

Identification of Genes Contributing to a Long Circadian Period in Drosophila Melanogaster.

The endogenous circadian period of animals and humans is typically very close to 24 h. Individuals with much longer circadian periods have been observed, however, and in the case of humans, these deviations have health implications. Previously, we observed a line of Drosophila with a very long average period of 31.3 h for locomotor activity behavior. Preliminary mapping indicated that the long period did not map to known canonical clock genes but instead mapped to multiple chromosomes. Using RNA-Seq, we surveyed the whole transcriptome of fly heads from this line across time and compared it with a wild-type control. A three-way generalized linear model revealed that approximately two-thirds of the genes were expressed differentially among the two genotypes, while only one quarter of the genes varied across time. Using these results, we applied algorithms to search for genes that oscillated over 24 h, identifying genes not previously known to cycle. We identified 166 differentially expressed genes that overlapped with a previous Genome-wide Association Study (GWAS) of circadian behavior, strongly implicating them in the long-period phenotype. We tested mutations in 45 of these genes for their effect on the circadian period. Mutations in Alk, alph, CG10089, CG42540, CG6034, Kairos (CG6123), CG8768, klg, Lar, sick, and tinc had significant effects on the circadian period, with seven of these mutations increasing the circadian period of locomotor activity behavior. Genetic rescue of mutant Kairos restored the circadian period to wild-type levels, suggesting it has a critical role in determining period length in constant darkness.

Natural selection on sleep duration in Drosophila melanogaster.

Sleep is ubiquitous across animal species, but why it persists is not well understood. Here we observe natural selection act on Drosophila sleep by relaxing bi-directional artificial selection for extreme sleep duration for 62 generations. When artificial selection was suspended, sleep increased in populations previously selected for short sleep. Likewise, sleep decreased in populations previously selected for long sleep when artificial selection was relaxed. We measured the corresponding changes in the allele frequencies of genomic variants responding to artificial selection. The allele frequencies of these variants reversed course in response to relaxed selection, and for short sleepers, the changes exceeded allele frequency changes that would be expected under random genetic drift. These observations suggest that the variants are causal polymorphisms for sleep duration responding to natural selection pressure. These polymorphisms may therefore pinpoint the most important regions of the genome maintaining variation in sleep duration.

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.

Phenotypic variation and natural selection at catsup, a pleiotropic quantitative trait gene in Drosophila.

Quantitative traits are shaped by networks of pleiotropic genes . To understand the mechanisms that maintain genetic variation for quantitative traits in natural populations and to predict responses to artificial and natural selection, we must evaluate pleiotropic effects of underlying quantitative trait genes and define functional allelic variation at the level of quantitative trait nucleotides (QTNs). Catecholamines up (Catsup), which encodes a negative regulator of tyrosine hydroxylase , the rate-limiting step in the synthesis of the neurotransmitter dopamine, is a pleiotropic quantitative trait gene in Drosophila melanogaster. We used association mapping to determine whether the same or different QTNs at Catsup are associated with naturally occurring variation in multiple quantitative traits. We sequenced 169 Catsup alleles from a single population and detected 33 polymorphisms with little linkage disequilibrium (LD). Different molecular polymorphisms in Catsup are independently associated with variation in longevity, locomotor behavior, and sensory bristle number. Most of these polymorphisms are potentially functional variants in protein coding regions, have large effects, and are not common. Thus, Catsup is a pleiotropic quantitative trait gene, but individual QTNs do not have pleiotropic effects. Molecular population genetic analyses of Catsup sequences are consistent with balancing selection maintaining multiple functional polymorphisms.

Quantitative genetic analysis of sleep in Drosophila melanogaster.

Although intensively studied, the biological purpose of sleep is not known. To identify candidate genes affecting sleep, we assayed 136 isogenic P-element insertion lines of Drosophila melanogaster. Since sleep has been negatively correlated with energy reserves across taxa, we measured energy stores (whole-body protein, glycogen, and triglycerides) in these lines as well. Twenty-one insertions with known effects on physiology, development, and behavior affect 24-hr sleep time. Thirty-two candidate insertions significantly impact energy stores. Mutational genetic correlations among sleep parameters revealed that the genetic basis of the transition between sleep and waking states in males and females may be different. Furthermore, sleep bout number can be decoupled from waking activity in males, but not in females. Significant genetic correlations are present between sleep phenotypes and glycogen stores in males, while sleep phenotypes are correlated with triglycerides in females. Differences observed in male and female sleep behavior in flies may therefore be related to sex-specific differences in metabolic needs. Sleep thus emerges as a complex trait that exhibits extensive pleiotropy and sex specificity. The large mutational target that we observed implicates genes functioning in a variety of biological processes, suggesting that sleep may serve a number of different functions rather than a single purpose.