A genome-wide association study for nutritional indices in Drosophila

DGRPool, a web tool leveraging harmonized Drosophila Genetic Reference Panel phenotyping data for the study of complex traits

Genome-wide association studies have advanced our understanding of complex traits, but studying how a GWAS variant can affect a specific trait in the human population remains challenging due to environmental variability. Drosophila melanogaster is in this regard an excellent model organism for studying the relationship between genetic and phenotypic variation due to its simple handling, standardized growth conditions, low cost, and short lifespan. The Drosophila Genetic Reference Panel (DGRP) in particular has been a valuable tool for studying complex traits, but proper harmonization and indexing of DGRP phenotyping data is necessary to fully capitalize on this resource. To address this, we created a web tool called DGRPool (dgrpool.epfl.ch), which aggregates phenotyping data of 1034 phenotypes across 135 DGRP studies in a common environment. DGRPool enables users to download data and run various tools such as genome-wide (GWAS) and phenome-wide (PheWAS) association studies. As a proof-of-concept, DGRPool was used to study the longevity phenotype and uncovered both established and unexpected correlations with other phenotypes such as locomotor activity, starvation resistance, desiccation survival, and oxidative stress resistance. DGRPool has the potential to facilitate new genetic and molecular insights of complex traits in Drosophila and serve as a valuable, interactive tool for the scientific community.

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

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

Genetic Background Matters: Population-Based Studies in Model Organisms for Translational Research

We are all similar but a bit different. These differences are partially due to variations in our genomes and are related to the heterogeneity of symptoms and responses to treatments that patients exhibit. Most animal studies are performed in one single strain with one manipulation. However, due to the lack of variability, therapies are not always reproducible when treatments are translated to humans. Panels of already sequenced organisms are valuable tools for mimicking human phenotypic heterogeneities and gene mapping. This review summarizes the current knowledge of mouse, fly, and yeast panels with insightful applications for translational research.

Anesthetic Preconditioning of Traumatic Brain Injury Is Ineffective in a Drosophila Model of Obesity

We tested the hypothesis that obesity influences the pharmacodynamics of volatile general anesthetics (VGAs) by comparing effects of anesthetic exposure on mortality from traumatic brain injury (TBI) in lean and obese Drosophila melanogaster We induced TBI with a high-impact trauma device. Starvation-selection over multiple generations resulted in an obese phenotype (SS flies). Fed flies served as lean controls (FC flies). Adult (1-8-day-old) SS and FC flies were exposed to equianesthetic doses of isoflurane or sevoflurane either before or after TBI. The principal outcome was percent mortality 24 hours after injury, expressed as the Mortality Index at 24 hours (MI24). TBI resulted in a lower MI24 in FC than in SS flies [21 (2.35) and 57.8 (2.14), respectively n = 12, P = 0.0001]. Pre-exposure to isoflurane or sevoflurane preconditioned FC flies to TBI, reducing the risk of death to 0.53 (0.25 to 1.13) and 0.82 (0.43 to 1.58), respectively, but had no preconditioning effect in SS flies. Postexposure to isoflurane or sevoflurane increased the risk of death in SS flies, but only postexposure to isoflurane increased the risk in FC flies [1.39 (0.81 to 2.38)]. Thus, obesity affects the pharmacodynamics of VGAs, thwarting the preconditioning effect of isoflurane and sevoflurane in TBI. SIGNIFICANCE STATEMENT: Inadvertent preconditioning in models of traumatic brain injury (TBI) is a recognized confounder. The findings in a fruit fly (Drosophila melanogaster) model of closed-head TBI indicate that anesthetic pharmacodynamics are profoundly affected by obesity. Specifically, obesity thwarts the brain-protective effect of anesthetic preconditioning. This finding is important for experimental studies of TBI and supports the versatility of the fruit fly as a model for the exploration of anesthetic pharmacodynamics in a wide parameter space.

Genome-wide analysis in Drosophila reveals diet-by-gene interactions and uncovers diet-responsive genes

Genetic and environmental factors play a major role in metabolic health. However, they do not act in isolation, as a change in an environmental factor such as diet may exert different effects based on an individual’s genotype. Here, we sought to understand how such gene–diet interactions influenced nutrient storage and utilization, a major determinant of metabolic disease. We subjected 178 inbred strains from the Drosophila genetic reference panel (DGRP) to diets varying in sugar, fat, and protein. We assessed starvation resistance, a holistic phenotype of nutrient storage and utilization that can be robustly measured. Diet influenced the starvation resistance of most strains, but the effect varied markedly between strains such that some displayed better survival on a high carbohydrate diet (HCD) compared to a high-fat diet while others had opposing responses, illustrating a considerable gene × diet interaction. This demonstrates that genetics plays a major role in diet responses. Furthermore, heritability analysis revealed that the greatest genetic variability arose from diets either high in sugar or high in protein. To uncover the genetic variants that contribute to the heterogeneity in starvation resistance, we mapped 566 diet-responsive SNPs in 293 genes, 174 of which have human orthologs. Using whole-body knockdown, we identified two genes that were required for glucose tolerance, storage, and utilization. Strikingly, flies in which the expression of one of these genes, CG4607 a putative homolog of a mammalian glucose transporter, was reduced at the whole-body level, displayed lethality on a HCD. This study provides evidence that there is a strong interplay between diet and genetics in governing survival in response to starvation, a surrogate measure of nutrient storage efficiency and obesity. It is likely that a similar principle applies to higher organisms thus supporting the case for nutrigenomics as an important health strategy.

Identification of genetic modifiers of lifespan on a high sugar diet in the Drosophila Genetic Reference Panel

Genome-wide association studies (GWAS) have become beneficial in identifying genetic variants underlying susceptibility to various complex diseases and conditions, including obesity. Utilizing the Drosophila Genetic Reference Panel (DGRP), we performed a GWAS of lifespan of 193 genetically distinct lines on a high sugar diet (HSD). The DGRP analysis pipeline determined the most significant lifespan associated polymorphisms were within loci of genes involved in: neural processes, behavior, development, and apoptosis, among other functions. Next, based on the relevance to obesity pathology, and the availability of transgenic RNAi lines targeting the genes we identified, whole-body in vivo knockdown of several candidate genes was performed. We utilized the GAL4-UAS binary expression system to independently validate the impacts of these loci on Drosophila lifespan during HSD. These loci were largely confirmed to affect lifespan in that HSD setting, as well as a normal diet setting. However, we also detected unexpected dietary effects of the HSD, including inconsistent diet effects on lifespan relative to a normal diet and a strong downregulation of feeding quantity.

Mitochondrial haplotypes affect metabolic phenotypes in the Drosophila Genetic Reference Panel

The nature and extent of mitochondrial DNA variation in a population and how it affects traits is poorly understood. Here we resequence the mitochondrial genomes of 169 Drosophila Genetic Reference Panel lines, identifying 231 variants that stratify along 12 mitochondrial haplotypes. We identify 1,845 cases of mitonuclear allelic imbalances, thus implying that mitochondrial haplotypes are reflected in the nuclear genome. However, no major fitness effects are associated with mitonuclear imbalance, suggesting that such imbalances reflect population structure at the mitochondrial level rather than genomic incompatibilities. Although mitochondrial haplotypes have no direct impact on mitochondrial respiration, some haplotypes are associated with stress- and metabolism-related phenotypes, including food intake in males. Finally, through reciprocal swapping of mitochondrial genomes, we demonstrate that a mitochondrial haplotype associated with high food intake can rescue a low food intake phenotype. Together, our findings provide new insight into population structure at the mitochondrial level and point to the importance of incorporating mitochondrial haplotypes in genotype-phenotype relationship studies.

Genome-Wide Association Analysis of Anoxia Tolerance in Drosophila melanogaster

As the genetic bases to variation in anoxia tolerance are poorly understood, we used the Drosophila Genetics Reference Panel (DGRP) to conduct a genome-wide association study (GWAS) of anoxia tolerance in adult and larval Drosophila melanogaster Survival ranged from 0-100% in adults exposed to 6 h of anoxia and from 20-98% for larvae exposed to 1 h of anoxia. Anoxia tolerance had a broad-sense heritability of 0.552 in adults and 0.433 in larvae. Larval and adult phenotypes were weakly correlated but the anoxia tolerance of adult males and females were strongly correlated. The GWA identified 180 SNPs in adults and 32 SNPs in larvae associated with anoxia tolerance. Gene ontology enrichment analysis indicated that many of the 119 polymorphic genes associated with adult anoxia-tolerance were associated with ionic transport or immune function. In contrast, the 22 polymorphic genes associated with larval anoxia-tolerance were mostly associated with regulation of transcription and DNA replication. RNAi of mapped genes generally supported the hypothesis that disruption of these genes reduces anoxia tolerance. For two ion transport genes, we tested predicted directional and sex-specific effects of SNP alleles on adult anoxia tolerance and found strong support in one case but not the other. Correlating our phenotype to prior DGRP studies suggests that genes affecting anoxia tolerance also influence stress-resistance, immune function and ionic balance. Overall, our results provide evidence for multiple new potential genetic influences on anoxia tolerance and provide additional support for important roles of ion balance and immune processes in determining variation in anoxia tolerance.

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.

Genetic basis of thermal plasticity variation in Drosophila melanogaster body size

Body size is a quantitative trait that is closely associated to fitness and under the control of both genetic and environmental factors. While developmental plasticity for this and other traits is heritable and under selection, little is known about the genetic basis for variation in plasticity that can provide the raw material for its evolution. We quantified genetic variation for body size plasticity in Drosophila melanogaster by measuring thorax and abdomen length of females reared at two temperatures from a panel representing naturally segregating alleles, the Drosophila Genetic Reference Panel (DGRP). We found variation between genotypes for the levels and direction of thermal plasticity in size of both body parts. We then used a Genome-Wide Association Study (GWAS) approach to unravel the genetic basis of inter-genotype variation in body size plasticity, and used different approaches to validate selected QTLs and to explore potential pleiotropic effects. We found mostly “private QTLs”, with little overlap between the candidate loci underlying variation in plasticity for thorax versus abdomen size, for different properties of the plastic response, and for size versus size plasticity. We also found that the putative functions of plasticity QTLs were diverse and that alleles for higher plasticity were found at lower frequencies in the target population. Importantly, a number of our plasticity QTLs have been targets of selection in other populations. Our data sheds light onto the genetic basis of inter-genotype variation in size plasticity that is necessary for its evolution.

Environmental conditions can influence development and lead to the production of phenotypes adjusted to the conditions adults will live in. This developmental plasticity, which can help organisms cope with environmental heterogeneity, is heritable and under selection. Its evolution will depend on available genetic variation. Using a panel of D. melanogaster flies representing naturally segregating alleles, we identified DNA sequence variants associated to variation in thermal plasticity for body size. We found that these variants correspond to a diverse set of functions and that their effects differ between body parts and properties of the thermal response. Our results shed new light onto the long discussed genes for plasticity.

The Sleep Inbred Panel, a Collection of Inbred Drosophila melanogaster with Extreme Long and Short Sleep Duration

Understanding how genomic variation causes differences in observable phenotypes remains a major challenge in biology. It is difficult to trace the sequence of events originating from genomic variants to changes in transcriptional responses or protein modifications. Ideally, one would conduct experiments with individuals that are at either extreme of the trait of interest, but such resources are often not available. Further, advances in genome editing will enable testing of candidate polymorphisms individually and in combination. Here we have created a resource for the study of sleep with 39 inbred lines of Drosophila-the Sleep Inbred Panel (SIP). SIP lines have stable long- and short-sleeping phenotypes developed from naturally occurring polymorphisms. These lines are fully sequenced, enabling more accurate targeting for genome editing and transgenic constructs. This panel facilitates the study of intermediate transcriptional and proteomic correlates of sleep, and supports genome editing studies to verify polymorphisms associated with sleep duration.

Obesity and Aging in the Drosophila Model

Being overweight increases the risk of many metabolic disorders, but how it affects lifespan is not completely clear. Not all obese people become ill, and the exact mechanism that turns excessive fat storage into a health-threatening state remains unknown. Drosophila melanogaster has served as an excellent model for many diseases, including obesity, diabetes, and hyperglycemia-associated disorders, such as cardiomyopathy or nephropathy. Here, we review the connections between fat storage and aging in different types of fly obesity. Whereas obesity induced by high-fat or high-sugar diet is associated with hyperglycemia, cardiomyopathy, and in some cases, shortening of lifespan, there are also examples in which obesity correlates with longevity. Transgenic lines with downregulations of the insulin/insulin-like growth factor (IIS) and target of rapamycin (TOR) signaling pathways, flies reared under dietary restriction, and even certain longevity selection lines are obese, yet long-lived. The mechanisms that underlie the differential lifespans in distinct types of obesity remain to be elucidated, but fat turnover, inflammatory pathways, and dysregulations of glucose metabolism may play key roles. Altogether, Drosophila is an excellent model to study the physiology of adiposity in both health and disease.

The impact of genome variation and diet on the metabolic phenotype and microbiome composition of Drosophila melanogaster

The metabolic phenotype of an organism depends on a complex regulatory network, which integrates the plethora of intrinsic and external information and prioritizes the flow of nutrients accordingly. Given the rise of metabolic disorders including obesity, a detailed understanding of this regulatory network is in urgent need. Yet, our level of understanding is far from completeness and complicated by the discovery of additional layers in metabolic regulation, such as the impact of the microbial community present in the gut on the hosts’ energy storage levels. Here, we investigate the interplay between genome variation, diet and the gut microbiome in the shaping of a metabolic phenotype. For this purpose, we reared a set of fully sequenced wild type Drosophila melanogaster flies under basal and nutritionally challenged conditions and performed metabolic and microbiome profiling experiments. Our results introduce the fly as a model system to investigate the impact of genome variation on the metabolic response to diet alterations and reveal candidate single nucleotide polymorphisms associated with different metabolic traits, as well as metabolite-metabolite and metabolite-microbe correlations. Intriguingly, the dietary changes affected the microbiome composition less than anticipated. These results challenge the current view of a rapidly changing microbiome in response to environmental fluctuations.

A Genomic Reference Panel for Drosophila serrata

Here we describe a collection of re-sequenced inbred lines of Drosophila serrata, sampled from a natural population situated deep within the species endemic distribution in Brisbane, Australia. D. serrata is a member of the speciose montium group whose members inhabit much of south east Asia and has been well studied for aspects of climatic adaptation, sexual selection, sexual dimorphism, and mate recognition. We sequenced 110 lines that were inbred via 17-20 generations of full-sib mating at an average coverage of 23.5x with paired-end Illumina reads. 15,228,692 biallelic SNPs passed quality control after being called using the Joint Genotyper for Inbred Lines (JGIL). Inbreeding was highly effective and the average levels of residual heterozygosity (0.86%) were well below theoretical expectations. As expected, linkage disequilibrium decayed rapidly, with r² dropping below 0.1 within 100 base pairs. With the exception of four closely related pairs of lines which may have been due to technical errors, there was no statistical support for population substructure. Consistent with other endemic populations of other Drosophila species, preliminary population genetic analyses revealed high nucleotide diversity and, on average, negative Tajima’s D values. A preliminary GWAS was performed on a cuticular hydrocarbon trait, 2-Me-C₂₈ revealing 4 SNPs passing Bonferroni significance residing in or near genes. One gene Cht9 may be involved in the transport of CHCs from the site of production (oenocytes) to the cuticle. Our panel will facilitate broader population genomic and quantitative genetic studies of this species and serve as an important complement to existing D. melanogaster panels that can be used to test for the conservation of genetic architectures across the Drosophila genus.

Genotype Influences Day-to-Day Variability in Sleep in Drosophila melanogaster

Patterns of sleep often vary among individuals. But sleep and activity may also vary within an individual, fluctuating in pattern across time. One possibility is that these daily fluctuations in sleep are caused by the underlying genotype of the individual. However, differences attributable to genetic causes are difficult to distinguish from environmental factors in outbred populations such as humans. We therefore employed Drosophila as a model of intra-individual variability in sleep using previously collected sleep and activity data from the Drosophila Genetic Reference Panel, a collection of wild-derived inbred lines. Individual flies had significant daily fluctuations in their sleep patterns, and these fluctuations were heritable. Using the standard deviation of sleep parameters as a metric, we conducted a genome-wide association study. We found 663 polymorphisms in 104 genes associated with daily fluctuations in sleep. We confirmed the effects of 12 candidate genes on the standard deviation of sleep parameters. Our results suggest that daily fluctuations in sleep patterns are due in part to gene activity.

Population genetic diversity in zebrafish lines

Toxicological and pharmacological researchers have seized upon the many benefits of zebrafish, including the short generation time, well-characterized development, and early maturation as clear embryos. A major difference from many model organisms is that standard husbandry practices in zebrafish are designed to maintain population diversity. While this diversity is attractive for translational applications in human and ecological health, it raises critical questions on how interindividual genetic variation might contribute to chemical exposure or disease susceptibility differences. Findings from pooled samples of zebrafish support this supposition of diversity yet cannot directly measure allele frequencies for reference versus alternate alleles. Using the Tanguay lab Tropical 5D zebrafish line (T5D), we performed whole genome sequencing on a large group (n = 276) of individual zebrafish embryos. Paired-end reads were collected on an Illumina 3000HT, then aligned to the most recent zebrafish reference genome (GRCz10). These data were used to compare observed population genetic variation across species (humans, mice, zebrafish), then across lines within zebrafish. We found more single nucleotide polymorphisms (SNPs) in T5D than have been reported in SNP databases for any of the WIK, TU, TL, or AB lines. We theorize that some subset of the novel SNPs may be shared with other zebrafish lines but have not been identified in other studies due to the limitations of capturing population diversity in pooled sequencing strategies. We establish T5D as a model that is representative of diversity levels within laboratory zebrafish lines and demonstrate that experimental design and analysis can exert major effects when characterizing genetic diversity in heterogeneous populations.

The genetics of egg retention and fertilization success in Drosophila: One step closer to understanding the transition from facultative to obligate viviparity

Oviparous, facultative egg retention enables Drosophila females to withhold fertilized eggs in their reproductive tracts until circumstances favor oviposition. The propensity to retain fertilized eggs varies greatly between species, and is correlated with other reproductive traits, such as egg size and ovariole number. While previous studies have described the phenomenon, no study to date has characterized within-species variation or the genetic basis of the trait. Here, we develop a novel microscope-based method for measuring egg retention in Drosophila females and determine the range of phenotypic variation in mated female egg retention in a subset of 91 Drosophila Genetic Reference Panel (DGRP) lines. We inferred the genetic basis of egg retention using a genome-wide association study (GWAS). Further, the scoring of more than 95,000 stained, staged eggs enabled estimates of fertilization success for each line. We found evidence that ovary- and spermathecae-related genes as well as genes affecting olfactory behavior, male mating behavior, male-female attraction and sperm motility may play a crucial role in post-mating physiology. Based on our findings we also propose potential evolutionary routes toward obligate viviparity. In particular, we propose that the loss of fecundity incurred by viviparity could be offset by benefits arising from enhanced mate discrimination, resource specialization, or modified egg morphology.

Genetic Architecture of Natural Variation Underlying Adult Foraging Behavior That Is Essential for Survival of Drosophila melanogaster

Foraging behavior is critical for the fitness of individuals. However, the genetic basis of variation in foraging behavior and the evolutionary forces underlying such natural variation have rarely been investigated. We developed a systematic approach to assay the variation in survival rate in a foraging environment for adult flies derived from a wild Drosophila melanogaster population. Despite being such an essential trait, there is substantial variation of foraging behavior among D. melanogaster strains. Importantly, we provided the first evaluation of the potential caveats of using inbred Drosophila strains to perform genome-wide association studies on life-history traits, and concluded that inbreeding depression is unlikely a major contributor for the observed large variation in adult foraging behavior. We found that adult foraging behavior has a strong genetic component and, unlike larval foraging behavior, depends on multiple loci. Identified candidate genes are enriched in those with high expression in adult heads and, demonstrated by expression knock down assay, are involved in maintaining normal functions of the nervous system. Our study not only identified candidate genes for foraging behavior that is relevant to individual fitness, but also shed light on the initial stage underlying the evolution of the behavior.

Charting the genotype-phenotype map: lessons from the Drosophila melanogaster Genetic Reference Panel

Understanding the genetic architecture (causal molecular variants, their effects, and frequencies) of quantitative traits is important for precision agriculture and medicine and predicting adaptive evolution, but is challenging in most species. The Drosophila melanogaster Genetic Reference Panel (DGRP) is a collection of 205 inbred strains with whole genome sequences derived from a single wild population in Raleigh, NC, USA. The large amount of quantitative genetic variation, lack of population structure, and rapid local decay of linkage disequilibrium in the DGRP and outbred populations derived from DGRP lines present a favorable scenario for performing genome-wide association (GWA) mapping analyses to identify candidate causal genes, polymorphisms, and pathways affecting quantitative traits. The many GWA studies utilizing the DGRP have revealed substantial natural genetic variation for all reported traits, little evidence for variants with large effects but enrichment for variants with low P-values, and a tendency for lower frequency variants to have larger effects than more common variants. The variants detected in the GWA analyses rarely overlap those discovered using mutagenesis, and often are the first functional annotations of computationally predicted genes. Variants implicated in GWA analyses typically have sex-specific and genetic background-specific (epistatic) effects, as well as pleiotropic effects on other quantitative traits. Studies in the DGRP reveal substantial genetic control of environmental variation. Taking account of genetic architecture can greatly improve genomic prediction in the DGRP. These features of the genetic architecture of quantitative traits are likely to apply to other species, including humans. WIREs Dev Biol 2018, 7:e289. doi: 10.1002/wdev.289 This article is categorized under: Invertebrate Organogenesis > Flies.

Charting the genotype-phenotype map: lessons from the Drosophila melanogaster Genetic Reference Panel

Understanding the genetic architecture (causal molecular variants, their effects, and frequencies) of quantitative traits is important for precision agriculture and medicine and predicting adaptive evolution, but is challenging in most species. The Drosophila melanogaster Genetic Reference Panel (DGRP) is a collection of 205 inbred strains with whole genome sequences derived from a single wild population in Raleigh, NC, USA. The large amount of quantitative genetic variation, lack of population structure, and rapid local decay of linkage disequilibrium in the DGRP and outbred populations derived from DGRP lines present a favorable scenario for performing genome-wide association (GWA) mapping analyses to identify candidate causal genes, polymorphisms, and pathways affecting quantitative traits. The many GWA studies utilizing the DGRP have revealed substantial natural genetic variation for all reported traits, little evidence for variants with large effects but enrichment for variants with low P-values, and a tendency for lower frequency variants to have larger effects than more common variants. The variants detected in the GWA analyses rarely overlap those discovered using mutagenesis, and often are the first functional annotations of computationally predicted genes. Variants implicated in GWA analyses typically have sex-specific and genetic background-specific (epistatic) effects, as well as pleiotropic effects on other quantitative traits. Studies in the DGRP reveal substantial genetic control of environmental variation. Taking account of genetic architecture can greatly improve genomic prediction in the DGRP. These features of the genetic architecture of quantitative traits are likely to apply to other species, including humans. WIREs Dev Biol 2018, 7:e289. doi: 10.1002/wdev.289 This article is categorized under: Invertebrate Organogenesis > Flies.

Charting the genotype-phenotype map: lessons from the Drosophila melanogaster Genetic Reference Panel

Understanding the genetic architecture (causal molecular variants, their effects, and frequencies) of quantitative traits is important for precision agriculture and medicine and predicting adaptive evolution, but is challenging in most species. The Drosophila melanogaster Genetic Reference Panel (DGRP) is a collection of 205 inbred strains with whole genome sequences derived from a single wild population in Raleigh, NC, USA. The large amount of quantitative genetic variation, lack of population structure, and rapid local decay of linkage disequilibrium in the DGRP and outbred populations derived from DGRP lines present a favorable scenario for performing genome-wide association (GWA) mapping analyses to identify candidate causal genes, polymorphisms, and pathways affecting quantitative traits. The many GWA studies utilizing the DGRP have revealed substantial natural genetic variation for all reported traits, little evidence for variants with large effects but enrichment for variants with low P-values, and a tendency for lower frequency variants to have larger effects than more common variants. The variants detected in the GWA analyses rarely overlap those discovered using mutagenesis, and often are the first functional annotations of computationally predicted genes. Variants implicated in GWA analyses typically have sex-specific and genetic background-specific (epistatic) effects, as well as pleiotropic effects on other quantitative traits. Studies in the DGRP reveal substantial genetic control of environmental variation. Taking account of genetic architecture can greatly improve genomic prediction in the DGRP. These features of the genetic architecture of quantitative traits are likely to apply to other species, including humans. WIREs Dev Biol 2018, 7:e289. doi: 10.1002/wdev.289 This article is categorized under: Invertebrate Organogenesis > Flies.

Variation in Position Effect Variegation Within a Natural Population

Changes in chromatin state may drive changes in gene expression, and it is of growing interest to understand the population genetic forces that drive differences in chromatin state. Here, we use the phenomenon of position effect variegation (PEV), a well-studied proxy for chromatin state, to survey variation in PEV among a naturally derived population. Further, we explore the genetic architecture of natural variation in factors that modify PEV. While previous mutation screens have identified over 150 suppressors and enhancers of PEV, it remains unknown to what extent allelic variation in these modifiers mediate interindividual variation in PEV. Is natural variation in PEV mediated by segregating genetic variation in known Su(var) and E(var) genes, or is the trait polygenic, with many variants mapping elsewhere in the genome? We designed a dominant mapping study that directly answers this question and suggests that the bulk of the variance in PEV does not map to genes with prior annotated impact to PEV. Instead, we find enrichment of top P-value ranked associations that suggest impact to active promoter and transcription start site proximal regions. This work highlights extensive variation in PEV within a population, and provides a quantitative view of the role naturally segregating autosomal variants play in modifying PEV-a phenomenon that continues to shape our understanding of chromatin state and epigenetics.

Leger RJ. The genetic basis for variation in resistance to infection in the Drosophila melanogaster genetic reference panel

Individuals vary extensively in the way they respond to disease but the genetic basis of this variation is not fully understood. We found substantial individual variation in resistance and tolerance to the fungal pathogen Metarhizium anisopliae Ma549 using the Drosophila melanogaster Genetic Reference Panel (DGRP). In addition, we found that host defense to Ma549 was correlated with defense to the bacterium Pseudomonas aeruginosa Pa14, and several previously published DGRP phenotypes including oxidative stress sensitivity, starvation stress resistance, hemolymph glucose levels, and sleep indices. We identified polymorphisms associated with differences between lines in both their mean survival times and microenvironmental plasticity, suggesting that lines differ in their ability to adapt to variable pathogen exposures. The majority of polymorphisms increasing resistance to Ma549 were sex biased, located in non-coding regions, had moderately large effect and were rare, suggesting that there is a general cost to defense. Nevertheless, host defense was not negatively correlated with overall longevity and fecundity. In contrast to Ma549, minor alleles were concentrated in the most Pa14-susceptible as well as the most Pa14-resistant lines. A pathway based analysis revealed a network of Pa14 and Ma549-resistance genes that are functionally connected through processes that encompass phagocytosis and engulfment, cell mobility, intermediary metabolism, protein phosphorylation, axon guidance, response to DNA damage, and drug metabolism. Functional testing with insertional mutagenesis lines indicates that 12/13 candidate genes tested influence susceptibility to Ma549. Many candidate genes have homologs identified in studies of human disease, suggesting that genes affecting variation in susceptibility are conserved across species.

Drosophila Genotype Influences Commensal Bacterial Levels

Host genotype can influence the composition of the commensal bacterial community in some organisms. Composition, however, is only one parameter describing a microbial community. Here, we test whether a second parameter—abundance of bacteria—is a heritable trait by quantifying the presence of four commensal bacterial strains within 36 gnotobiotic inbred lines of Drosophila melanogaster. We find that D. melanogaster genotype exerts a significant effect on microbial levels within the fly. When introduced as monocultures into axenic flies, three of the four bacterial strains were reliably detected within the fly. The amounts of these different strains are strongly correlated, suggesting that the host regulates commensal bacteria through general, not bacteria-specific, means. While the correlation does not appear to be driven by simple variation in overall gut dimensions, a genetic association study suggests that variation in commensal bacterial load may largely be attributed to physical aspects of host cell growth and development.

The Genetic Basis of Natural Variation in Drosophila (Diptera: Drosophilidae) Virgin Egg Retention

Drosophila melanogaster is able to thrive in harsh northern climates through adaptations in life-history traits and physiological mechanisms that allow for survival through the winter. We examined the genetic basis of natural variation in one such trait, female virgin egg retention, which was previously shown to vary clinally and seasonally. To further our understanding of the genetic basis and evolution of virgin egg retention, we performed a genome-wide association study (GWAS) using the previously sequenced Drosophila Genetic Reference Panel (DGRP) mapping population. We found 29 single nucleotide polymorphisms (SNPs) associated with virgin egg retention and assayed 6 available mutant lines, each harboring a mutation in a candidate gene, for effects on egg retention time. We found that four out of the six mutant lines had defects in egg retention time as compared with the respective controls: mun, T48, Mes-4, and Klp67A Surprisingly, none of these genes has a recognized role in ovulation control, but three of the four genes have known effects on fertility or have high expression in the ovaries. We also found that the SNP set associated with egg retention time was enriched for clinal SNPs. The majority of clinal SNPs had alleles associated with longer egg retention present at higher frequencies in higher latitudes. Our results support previous studies that show higher frequency of long retention times at higher latitude, providing evidence for the adaptive value of virgin egg-retention.

Cross-phenotype association tests uncover genes mediating nutrient response in Drosophila

Background

Obesity-related diseases are major contributors to morbidity and mortality in the developed world. Molecular diagnostics and targets of therapies to combat nutritional imbalance are urgently needed in the clinic. Invertebrate animals have been a cornerstone of basic research efforts to dissect the genetics of metabolism and nutrient response. We set out to use fruit flies reared on restricted and nutrient-rich diets to identify genes associated with starvation resistance, body mass and composition, in a survey of genetic variation across the Drosophila Genetic Reference Panel (DGRP).

Results

We measured starvation resistance, body weight and composition in DGRP lines on each of two diets and used several association mapping strategies to harness this panel of phenotypes for molecular insights. We tested DNA sequence variants for a relationship with single metabolic traits and with multiple traits at once, using a scheme for cross-phenotype association mapping; we focused our association tests on homologs of human disease genes and common polymorphisms; and we tested for gene-by-diet interactions. The results revealed gene and gene-by-diet associations between 17 variants and body mass, whole-body triglyceride and glucose content, or starvation resistance. Focused molecular experiments validated the role in body mass of an uncharacterized gene, CG43921 (which we rename heavyweight), and previously unknown functions for the diacylglycerol kinase rdgA, the huntingtin homolog htt, and the ceramide synthase schlank in nutrient-dependent body mass, starvation resistance, and lifespan.

Conclusions

Our findings implicate a wealth of gene candidates in fly metabolism and nutrient response, and ascribe novel functions to htt, rdgA, hwt and schlank.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3137-9) contains supplementary material, which is available to authorized users.

Effects of larval crowding on quantitative variation for development time and viability in Drosophila melanogaster

Competition between individuals belonging to the same species is a universal feature of natural populations and is the process underpinning organismal adaptation. Despite its importance, still comparatively little is known about the genetic variation responsible for competitive traits. Here, we measured the phenotypic variation and quantitative genetics parameters for two fitness-related traits-egg-to-adult viability and development time-across a panel of Drosophila strains under varying larval densities. Both traits exhibited substantial genetic variation at all larval densities, as well as significant genotype-by-environment interactions (GEIs). GEI was attributable to changes in the rank order of reaction norms for both traits, and additionally to differences in the between-line variance for development time. The coefficient of genetic variation increased under stress conditions for development time, while it was higher at both high and low densities for viability. While development time also correlated negatively with fitness at high larval densities-meaning that fast developers have high fitness-there was no correlation with fitness at low density. This result suggests that GEI may be a common feature of fitness-related genetic variation and, further, that trait values under noncompetitive conditions could be poor indicators of individual fitness. The latter point could have significant implications for animal and plant breeding programs, as well as for conservation genetics.

Genome-Wide Analysis Reveals Novel Regulators of Growth in Drosophila melanogaster

Organismal size depends on the interplay between genetic and environmental factors. Genome-wide association (GWA) analyses in humans have implied many genes in the control of height but suffer from the inability to control the environment. Genetic analyses in Drosophila have identified conserved signaling pathways controlling size; however, how these pathways control phenotypic diversity is unclear. We performed GWA of size traits using the Drosophila Genetic Reference Panel of inbred, sequenced lines. We find that the top associated variants differ between traits and sexes; do not map to canonical growth pathway genes, but can be linked to these by epistasis analysis; and are enriched for genes and putative enhancers. Performing GWA on well-studied developmental traits under controlled conditions expands our understanding of developmental processes underlying phenotypic diversity.

Genetic studies in Drosophila have elucidated conserved signaling pathways and environmental factors that together control organismal size. In humans, hundreds of genes are associated with height variation, but these associations have not been performed in a controlled environment. As a result we are still lacking an understanding of the mechanisms creating size variability within a species. Here, under carefully controlled environmental conditions, we identify naturally occurring genetic variants that are associated with size diversity in Drosophila. We identify a cluster of associations close to the kek1 locus, a well-characterized growth regulator, but otherwise find that most variants are located in or close to genes that do not belong to the conserved pathways but may interact with these in a biological network. We validate 33 novel growth regulatory genes that participate in diverse cellular processes, most notably cellular metabolism and cell polarity. This study is the first genome-wide association analysis of natural variants underlying size in Drosophila and our results complement the knowledge we have accumulated on this trait from mutational studies of single genes.

Candidate genetic modifiers of retinitis pigmentosa identified by exploiting natural variation in Drosophila

Individuals carrying the same pathogenic mutation can present with a broad range of disease outcomes. While some of this variation arises from environmental factors, it is increasingly recognized that the background genetic variation of each individual can have a profound effect on the expressivity of a pathogenic mutation. In order to understand this background effect on disease-causing mutations, studies need to be performed across a wide range of backgrounds. Recent advancements in model organism biology allow us to test mutations across genetically diverse backgrounds and identify the genes that influence the expressivity of a mutation. In this study, we used the Drosophila Genetic Reference Panel, a collection of ∼200 wild-derived strains, to test the variability of the retinal phenotype of the Rh1(G69D) Drosophila model of retinitis pigmentosa (RP). We found that the Rh1(G69D) retinal phenotype is quite a variable quantitative phenotype. To identify the genes driving this extensive phenotypic variation, we performed a genome-wide association study. We identified 106 candidate genes, including 14 high-priority candidates. Functional testing by RNAi indicates that 10/13 top candidates tested influence the expressivity of Rh1(G69D). The human orthologs of the candidate genes have not previously been implicated as RP modifiers and their functions are diverse, including roles in endoplasmic reticulum stress, apoptosis and retinal degeneration and development. This study demonstrates the utility of studying a pathogenic mutation across a wide range of genetic backgrounds. These candidate modifiers provide new avenues of inquiry that may reveal new RP disease mechanisms and therapies.

Quantitative Genetics of Food Intake in Drosophila melanogaster

Food intake is an essential animal activity, regulated by neural circuits that motivate food localization, evaluate nutritional content and acceptance or rejection responses through the gustatory system, and regulate neuroendocrine feedback loops that maintain energy homeostasis. Excess food consumption in people is associated with obesity and metabolic and cardiovascular disorders. However, little is known about the genetic basis of natural variation in food consumption. To gain insights in evolutionarily conserved genetic principles that regulate food intake, we took advantage of a model system, Drosophila melanogaster, in which food intake, environmental conditions and genetic background can be controlled precisely. We quantified variation in food intake among 182 inbred, sequenced lines of the Drosophila melanogaster Genetic Reference Panel (DGRP). We found significant genetic variation in the mean and within-line environmental variance of food consumption and observed sexual dimorphism and genetic variation in sexual dimorphism for both food intake traits (mean and variance). We performed genome wide association (GWA) analyses for mean food intake and environmental variance of food intake (using the coefficient of environmental variation, CV_E, as the metric for environmental variance) and identified molecular polymorphisms associated with both traits. Validation experiments using RNAi-knockdown confirmed 24 of 31 (77%) candidate genes affecting food intake and/or variance of food intake, and a test cross between selected DGRP lines confirmed a SNP affecting mean food intake identified in the GWA analysis. The majority of the validated candidate genes were novel with respect to feeding behavior, and many had mammalian orthologs implicated in metabolic diseases.

The complex contributions of genetics and nutrition to immunity in Drosophila melanogaster

Both malnutrition and undernutrition can lead to compromised immune defense in a diversity of animals, and "nutritional immunology" has been suggested as a means of understanding immunity and determining strategies for fighting infection. The genetic basis for the effects of diet on immunity, however, has been largely unknown. In the present study, we have conducted genome-wide association mapping in Drosophila melanogaster to identify the genetic basis for individual variation in resistance, and for variation in immunological sensitivity to diet (genotype-by-environment interaction, or GxE). D. melanogaster were reared for several generations on either high-glucose or low-glucose diets and then infected with Providencia rettgeri, a natural bacterial pathogen of D. melanogaster. Systemic pathogen load was measured at the peak of infection intensity, and several indicators of nutritional status were taken from uninfected flies reared on each diet. We find that dietary glucose level significantly alters the quality of immune defense, with elevated dietary glucose resulting in higher pathogen loads. The quality of immune defense is genetically variable within the sampled population, and we find genetic variation for immunological sensitivity to dietary glucose (genotype-by-diet interaction). Immune defense was genetically correlated with indicators of metabolic status in flies reared on the high-glucose diet, and we identified multiple genes that explain variation in immune defense, including several that have not been previously implicated in immune response but which are confirmed to alter pathogen load after RNAi knockdown. Our findings emphasize the importance of dietary composition to immune defense and reveal genes outside the conventional "immune system" that can be important in determining susceptibility to infection. Functional variation in these genes is segregating in a natural population, providing the substrate for evolutionary response to pathogen pressure in the context of nutritional environment.