Skip to main content

Advertisement

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

BMC Genomics volume 17, Article number: 867 (2016) Cite this article

Abstract

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.

Background

Obesity is a worldwide epidemic and among the leading causes of morbidity and mortality in the US, increasing susceptibility to chronic pathologies that impact many body systems [1]. Against a backdrop of landmark advances in dissections of the basic biology of metabolism [2], geneticists continue to search for novel genes that modulate body composition and diet responses, as potential targets for drug development. Such factors have emerged from mapping studies in human populations [3] and in mice reared on a high-fat diet [4–6], though only a handful have been validated at the single-gene level [3]. As a complement to mammals, invertebrate model organisms offer a far quicker route to the discovery and validation of nutrient response genes. In the fruit fly, classic mutagenesis screens have enabled landmark discoveries of obesity and body composition genes conserved in vertebrates [7–12]. These screens have been applied to one qualitative, easily-scored trait at a time. Consequently, geneticists have had an incomplete understanding of how multiple aspects of metabolism are coordinated in response to the organism’s changing needs, even for the well-studied insulin [13, 14] and target of rapamycin (TOR) [15] pathways in the fly and other model organisms.

The advent of the Drosophila Genetic Reference Panel (DGRP) enables the mapping of natural phenotypic variation to DNA sequence in flies with an association paradigm [16, 17]. Genome-wide association studies (GWAS), by using populations in which many natural genetic variants are segregating at once, can require orders of magnitude less phenotyping effort than systematic screens perturbing each gene in turn. As such, GWAS makes possible the genetic dissection of traits that require assays too nuanced or time-consuming for mutagenesis studies, including many aspects of metabolism and dietary response. And alongside its use in gene discovery, GWAS also reveals the genetic architecture of natural variation in a trait of interest. These strengths have motivated previous studies focused on metabolic rate and body size and composition in the fly, yielding a number of intriguing candidate genes and polymorphisms [16, 18–22].

A potential strategy to boost power in GWAS is the use of cross-phenotype association tests, to search for cases in which a genetic variant associates with multiple traits at once [23]. The recent success of principal component-based methods for this purpose in the fly [24] has led to a first set of candidate cases in which master regulatory loci associate jointly with multiple facets of metabolism. We reasoned that any such candidate would serve as a high-value candidate for experimental validation. We also expected that association tests for gene-by-diet interactions, which have not as yet been implemented in the fly [19, 21, 25, 26], could be harnessed to reveal novel genetic determinants of metabolism. We thus set out to use cross-phenotype association tests in the DGRP, along with single-phenotype GWAS approaches, to generate testable hypotheses about genes underlying metabolic traits and the response to changing diet.

Results

We surveyed natural genetic variation across genetically distinct fly lines in metabolic behaviors and their dependence on dietary yeast restriction. Dietary yeast levels are known to modulate fly body composition, stress resistance, and lifespan in laboratory strains [27–30]. We reared 171–181 DGRP lines on two diets (ad libitum, AL, 5 % yeast; dietary restriction, DR, 0.5 % yeast), using non-virgin females to avoid effects of mating status on diet–responsive phenotypes [31]. In each, we measured wet body mass, resistance to acute starvation, and whole-body triglyceride and glucose levels (Additional file 1: Figure S1). As expected [32–36], on average DR treatment was associated with an increase in all traits except wet body mass, which decreased in DR flies relative to AL flies (Additional file 1: Figure S1). Phenotype values of a given strain reared on the two diets were correlated (R 2 = 0.53-0.69; Additional file 1: Table S1), indicating a sizeable genetic contribution to phenotype, independent of diet. Indeed, heritabilities ranged from modest to high for the GWAS phenotypes (47–93 %; Additional file 1: Table S2). For most traits, the ratio of the respective values measured in flies reared on the AL and DR diets was also highly heritable (Additional file 1: Table S2).

GWAS of metabolic traits

To begin to identify the genetic basis of variation across DGRP lines in metabolic behaviors, we tested each trait in turn for association to each of the 1.9 million variants segregating in the population with minor allele frequency ≥5 %. Our association paradigm used a linear regression method that included terms in the model for the effect of genotype at a given locus, its interaction with diet, and diet itself, and we corrected for multiple testing via permutation. At an experiment-wise threshold for significance corresponding to false discovery rate ≤ 10 %, we mapped diet-independent effects of one polymorphic locus, nimb3, to measurements of resistance of acute starvation, and another locus, 14-3-3 epsilon, to whole-body triglyceride content (Table 1). The top-scoring polymorphisms in these genes explained 18–20 % of the genotypic variance in their respective phenotypes (Table 1). Together, these association results provided a first line of evidence that our GWAS could identify genes that govern variation in fly metabolism.

Table 1 Variants associated with metabolic traits
Full size table

We used the results from each genome-wide scan as input into tests of groups of genes of common function for enriched association signal, using GOglm [37]. Emerging from this analysis were eight GO terms whose genes tended to associate with different aspects of metabolism (Additional file 1: Table S3). Most notable was the enrichment of association signal, for body mass and glucose levels, in sensory perception and antennal morphogenesis genes, suggesting the possibility that some variants modulating these traits act by tuning animals’ perception of food [38–40].

Cross-phenotype association mapping

We reasoned that a complementary strategy to mine our panel of metabolic trait measurements in DGRP lines would be to test a given variant for association with multiple traits at once. For this purpose, we used each pair of phenotypes in turn (from among our set of four phenotypes in animals reared on each diet; Additional file 1: Figure S1) as input into the multiple-phenotype association test (SMAT; [41]). Genome-wide scans with this paradigm were not sufficiently well-powered to map loci at reasonable false discovery rates after permutation-based multiple testing correction (see Methods and Additional file 1: Table S4). To improve power, and motivated by the clinical interest in human diseases where nutrition is a risk factor, we developed a test pipeline specialized to fly homologs of a broad range of human disease genes (as curated by Online Mendelian Inheritance in Man [42]; see Methods). Our rationale was that any genes of this set, if mapped in our association calculations, would likely be only a subset of the full complement of loci causal for the respective trait, but would be of particular interest as candidates for translational applications. Scans of human disease loci using SMAT attained power in a number of trait pairs, when considered independently (Additional file 1: Table S4) and when corrected for multiple testing across pairs (see Methods). Emerging from the results was strong signal for six polymorphisms jointly associating with starvation resistance and body mass, in DR flies and, separately, in flies reared on the AL diet (Table 1). These polymorphisms were in eIF4G2, a paralog of the nutrient-sensing factor eIF4G [43]; the ceramide synthase gene schlank; and the Drosophila homolog of huntingtin, the gene underlying human Huntington’s disease (Table 1). These associations explained 8–12 % of the genetic variance in starvation resistance or body mass (Table 1). Cross-phenotype mapping also found associations between starvation resistance and triglyceride levels and two loci (Table 1): CG7560, an FAD-linked oxidoreductase involved in methionine metabolism [44], and Cht12, a chitinase involved in cuticle formation [45]. Tests for gene sets with enriched association signal in cross-phenotype tests revealed association between starvation resistance and glucose content and genes annotated in ventral cord development (Additional file 1: Table S3), raising the possibility that these characters could be affected by variants that act in the central nervous system, plausibly via locomotor hyperactivity during starvation [38–40]. These data established our cross-phenotype association tests as a useful tool in the mapping of genotype to phenotype for metabolic traits.

A new role for schlank in nutrient response and lifespan in adult flies

Among our cross-phenotype association hits, as a first test of functional validation we focused on schlank, given its known role as a determinant of lipid accumulation in fly larvae [46] and, in yeast, as a determinant of lifespan [47]. The association between intronic variants in schlank and starvation resistance and body mass reached experiment-wide significance in DR flies (Table 1) and was also apparent in flies on the AL diet (Fig. 1a-b). Since loss of schlank function compromises viability beyond the larval stage [46], to investigate its role in metabolic traits during adulthood we crossed a stock ubiquitously expressing the GeneSwitch (GS) regulator, which is inducible upon treatment with RU486 [48], to a TRiP background stock [49] harboring a GAL4-controlled schlank RNAi transgene. In the progeny, upon induction we confirmed schlank knockdown by qRT-PCR (Additional file 1: Figure S2). We first focused on acute starvation resistance, which increases in wild-type flies reared under chronic nutrient deprivation in a manner dependent on changes in fat metabolism [36]. Strikingly, schlank knockdown during adulthood completely abrogated this advantage of the DR treatment: starvation resistance in knockdown flies reared on DR was 59 % of wild-type controls, on par with that of animals reared under AL (Fig. 1c). Further testing revealed a diet-independent, 5–13 % reduction in body mass in schlank knockdown flies (Fig. 1d), and also a loss of 10 % in triglyceride content (Fig. 1e). Negative controls ruled out effects on these traits from RU486 treatment alone (Additional file 1: Figures S3 and S4). These data make clear that schlank is of critical importance for the enhanced starvation resistance among DR flies, likely via remodeling of fat storage and catabolism in flies reared on this diet [36]. Lifespan assays suggested a potential role for schlank in DR-induced longevity, as knockdown of this gene shortened lifespan to a slightly greater degree in flies reared on the restricted diet (24 versus 19 % on AL; Additional file 1: Figure S5). Together, our discoveries of a role for schlank in nutrient-dependent changes in starvation resistance, as well as body weight, provide a first compelling validation of multi-phenotype association signal in the DGRP as a signpost for genes with true metabolic function.

Fig. 1
figure1

The ceramide synthase schlank governs body mass, starvation, and lifespan. a and b Trait association to schlank genotype in the DGRP. In a given panel, each point reports phenotype measurements (y-axis) from one DGRP strain homozygous for the schlank allele indicated on the x-axis. Each black horizontal bar reports the mean phenotype across all strains harboring the indicated allele. c-e Validation of schlank function. c Each trace reports survival of a strain harboring a schlank RNAi construct and a GAL4-GeneSwitch (GS) activator under the daughterless promoter, pre-treated with RU486 (+ drug, dotted lines) or vehicle (no drug, solid lines) for 10 days on the indicated diet and then switched to acute starvation media. d Each column reports body mass per fly of the indicated strain, with genotypes as in c (n = 2 populations of 15 flies per condition). e Each column reports whole body triglyceride levels per fly of the indicated strain, with genotypes as in c (n = 1 population of 15 flies per condition). For d and e top, middle, and bottom horizontal bars of a given vertical box denote the respective quartiles over batches and/or technical replicates, and the top and bottom short horizontal bars report minimum and maximum, respectively. AL, ad libitum diet; DR, dietary restriction. *, p < 0.05, **, p < 0.01, ***, p < 10-3, ****, p <10-7, *****, p < 10-15. Red asterisks denote significance of the effect of the genetic perturbation in animals on AL food, blue denotes significance in animals on DR food, and purple denotes significance of the interaction between diet and genetic perturbation. For strain details see Additional file 1: Table S6

Full size image
Fig. 2
figure2

A metabolic function for fly huntingtin (htt). a and b Trait association to htt genotype in the DGRP. Data are as in Fig. 1a-b, except that an intronic indel in htt is analyzed. c and d Validation of htt function in laboratory strains. c Each trace reports survival of a strain expressing a GAL4 activator under the daughterless promoter and harboring an htt RNAi construct (RNAi) or no RNAi transgene (control), pre-treated for 10 days on the indicated diet and then switched to acute starvation media. d Each column reports body mass per fly of the indicated strain, with genotypes as in c. (n = 5 populations of 15 flies per condition). Symbols and abbreviations are as in Fig. 1. For strain details see Additional file 1: Table S6

Full size image

A metabolic role for fly htt

As a second case study for functional confirmation of our cross-phenotype mapping results, we focused on htt, the Drosophila huntingtin gene. In our survey of DGRP lines, a single-base insertion in an intron of fly htt reached significance for joint association to body mass and starvation resistance, in animals on the AL diet (Table 1), and a similar effect was discernable upon manual inspection in DR flies (Fig. 2a-b). A metabolic role for fly htt would be consistent with the symptoms of human Huntington’s patients, who, even on prescribed diets of thousands of excess calories a day, often fail to maintain body mass [50, 51]. To verify the dependence of fly nutrient responsive traits on htt, we bred a line with constitutive, ubiquitous expression of an htt RNAi construct, and confirmed knockdown in this strain by qRT-PCR (Additional file 1: Figure S2). This htt knockdown line exhibited marked improvements in starvation resistance relative to isogenic controls, in DR flies and, to a lesser degree, in animals reared on the AL diet (Fig. 2c). Knockdown of htt also conferred robust increases in body mass, with a slightly stronger effect in animals on AL food than on DR (24 versus 22 %, respectively; Fig. 2d), but had no significant effect on triglyceride content (Additional file 1: Figure S4). To test htt’s role in longevity while avoiding the potential for knockdown effects during development, we bred flies expressing drug-inducible, ubiquitous htt RNAi, and induced the transgene only after eclosion, confirming knockdown by qRT-PCR (Additional file 1: Figure S2). The resulting htt knockdown in adult animals conferred a 9 % lifespan extension on the AL diet and decreased lifespan in DR animals by 5 % relative to isogenic controls, effects that could not be attributed to inducer treatment alone (Additional file 1: Figure S5). These data make clear that reduced htt expression can alter longevity and promote starvation resistance and increase body mass, likely through mechanisms that do not involve changes to steady-state triglyceride levels. As nutrient restriction enhanced some of these effects, mitigated others, and had no detectable effect on the natural variant originally emerging in our association scan, htt likely participates in a complex relationship with other dietary response regulators.

Common variants associated with metabolic traits

As an additional strategy to improve mapping power in the DGRP, we developed a scheme for analysis of polymorphisms at which the minor allele was very common. We expected that common variants mapped to a given trait might only represent a fraction of its full genetic architecture, but we considered this restriction a justifiable trade-off in light of the expected high power in association tests of common variants where the number of observations among individuals bearing the minor allele is maximized. For a first evaluation of this strategy, we applied our linear model association test framework for a given trait, including genotype and gene-by-diet interaction effects, to the 734,670 markers in the DGRP with minor allele frequency ≥25 %. This genome-scale scan uncovered significant, diet-independent association between genotype at nimb3 and acute starvation resistance, as we had observed in scans of variants at lower allele frequencies; we also detected association to body mass at the uncharacterized gene CG43921 (Table 1). Gene-by-diet association was significant between whole-body glucose measurements and a locus encompassing the cleavage and polyadenylation factor pic and the uncharacterized gene CG7966 (Table 1). These loci explained up to 18 % of the variance in their respective traits (Table 1). Functional genomic tests revealed enrichment in association to starvation resistance among genes categories annotated in actin filament production (Additional file 1: Table S5). In the latter, top-scoring genes were Arp2, Arpc2, and SCAR (association p = 0.0008, 0.001, and 0.003 respectively), which have roles in the development of neural processes and axonal defasciculation [52, 53]. Also emerging from functional genomic tests was an enrichment in association to starvation resistance among heterophilic cell-cell adhesion genes (Additional file 1: Table S5); these included Notch (association p = 0.006), among whose many biological roles is that of axon guidance [54], and “beaten path” genes such as beat-Ib, beat-IIIb, and beat-Va (association p = 0.00001, p = 0.00008, p = 0.0002 respectively), which are involved in motor neuron defasciculation [55]. These results suggest that variants acting to tune development of the nervous system, if they modulate feeding or physical activity, could have downstream influences on metabolism [56].

CG43921/heavyweight is a novel determinant of body mass

For a validation of a top-scoring locus from our association tests of common variants, we focused on the uncharacterized gene CG43921, in which an intronic polymorphism associated with body mass across DGRP lines in a diet-independent manner (Table 1 and Fig. 3a). To verify the dependence of fly body mass on CG43921, we bred a line with constitutive expression of a CG43921 RNAi construct under a ubiquitous, constitutive driver, and we confirmed knockdown in this strain by qRT-PCR (Additional file 1: Figure S2). Knockdown of CG43921 by this strategy increased body mass by 9 % relative to isogenic controls in flies reared on the AL diet, with no significant effect in DR animals (Fig. 3b). We did not find significant changes in triglyceride levels (Additional file 1: Figure S4) or lifespan (Additional file 1: Figure S5) and we observed only minor extension of starvation resistance (Additional file 1: Figure S6) in CG43921 knockdown flies. These data implicate CG43921 primarily in the control of body mass, the trait that showed the strongest association signal, and as such we give this gene the common name heavyweight (hwt).

Fig. 3
figure3

CG43921/heavyweight is a novel determinant of body mass. a DGRP trait association with heavyweight (hwt) genotype. Data are as in Fig. 1a-b except that an intronic variant in hwt is analyzed. b Validation of hwt function in laboratory strains. Each column reports body mass per fly of a strain expressing a GAL4 activator under the daughterless promoter and harboring an hwt RNAi construct (RNAi) or no RNAi transgene (control) (n = 2 populations of 15 flies per condition). Symbols and abbreviations are as in Fig. 1. For strain details see Additional file 1: Table S6

Full size image

Diet-dependent metabolic function of rdgA

Paralleling our analyses of low-frequency variants, we used the set of common variants in the DGRP as input into tests for cross-phenotype associations, analyzing each pair of our phenotypes on each diet. Two significant single-locus mapping results emerged from analysis of flies reared on DR. Genotype at CG43462, a suspected determinant of fly lifespan [57], was associated with starvation resistance and triglyceride levels; the diacylglycerol kinase gene rdgA (Table 1), whose mammalian homologs have been implicated in diabetes and obesity risk [58, 59] and which also modulates fly lifespan [60], associated with body mass and starvation resistance. Functional genomic tests detected enriched association, in our joint analysis of starvation resistance and body mass, in two sets of genes annotated in alcohol catabolism (Additional file 1: Table S5). Top-scoring genes in these categories were three sphingomyelin phosphodiesterases (CG15533, CG15534, and CG3376; association p = 0.0099, 0.0012, and 0.0013, respectively) and two phosphatidylinositol phosphatases (Plip and Pten, association p = 0.0022 and 0.0056, respectively), highlighting the potential for variants throughout phospholipid pathways that could perturb metabolism in the DGRP.

We chose rdgA for an additional validation case study, on the strength of previous evidence that it interacts with the TOR pathway and modulates lifespan [60]. The natural variant segregating in the DGRP in rdgA was associated with stronger effects in animals reared on the DR diet than on AL (Fig. 4a and b), though this diet dependence did not reach significance in our starvation resistance GWAS. We reasoned that disruption of rdgA might affect metabolic characters in a diet-specific manner. To test this, we first phenotyped flies bearing a transposon insertion in rdgA [61]. Relative to background-matched controls, these mutants survived acute starvation for 51 % longer after rearing on the DR diet, with a more modest 16 % increase after pre-treatment on the AL diet (Fig. 4c). rdgA mutants also exhibited an increase in mass, with an effect significantly stronger in DR flies than those on AL food (24 and 3 %, respectively; Fig. 4d). Consistent with the gain in body mass on DR, there was a significant 12 % increase in triglyceride levels in rdgA mutant flies relative to controls on DR but not AL (Fig. 4e). To investigate the role of rdgA in lifespan, we used a drug-inducible RNAi construct, triggering knockdown only in adulthood. Mirroring results of a previous study in flies reared on a single diet [60], we found that neuron-specific knockdown of rdgA was sufficient for lifespan extension; however, our data revealed this effect to be most pronounced during DR (9 % extension), and almost undetectable on the AL diet (Additional file 1: Figure S5). We conclude that laboratory-induced perturbations to rdgA impact starvation resistance, body mass, and lifespan. The striking diet dependence of effects of these mutations, as with natural polymorphisms at the locus, reflect a likely role for rdgA in nutrient response, as expected given its interactions with TOR signaling [60].

Fig. 4
figure4

The diacylglycerol kinase rdgA governs body mass, starvation resistance, and lifespan. a and b Trait association to rdgA genotype in the DGRP. Data are as in Fig. 1a-b, except that an intronic variant in rdgA is analyzed. c-e Validation of rdgA function in laboratory strains. c Each trace reports survival of a strain harboring a transposon insertion in rdgA (mutant) or isogenic wild-type control, pre-treated for 10 days on the indicated diet and then switched to acute starvation media. d Each column reports body mass per fly of the indicated strain, with genotypes as in c (n = 2 populations of 15 flies per condition). e Each column reports whole body triglyceride levels normalized to body mass in the indicated strain, with genotypes as in c (n = 1 population of 15 flies per condition). Symbols and abbreviations are as in Fig. 1. For strain details see Additional file 1: Table S6

Full size image

Discussion

A direct consequence of obesity is the increased probability of age-related pathologies, including diabetes, cardiovascular disease, neurodegenerative disorders and cancer [1]. Conversely, fasting and dietary restriction prolong lifespan and delay the onset of chronic diseases [29, 62, 63]. The genes that underlie these effects are prime targets for pro-healthspan drugs and personalized diets, and natural variation mapping strategies are well suited for expedient screens to discover these factors. In this work, we have implemented several strategies for highly-powered GWAS in the fly, including a focus on homologs of genes involved in human disease and common variants, and tests of a given variant for joint relationships to multiple traits in parallel. Our work also provides a first window onto gene-by-diet interactions in the DGRP. In total, our association scans of common alleles and disease gene homologs have implicated 17 variants, in 12 genes, in nutrient response traits. And we identified groups of genes of common function with enriched association to metabolic traits, although formally, nuanced statistical biases would need to be ruled out as the potential source for any such pathway signal [64].

Among the single-locus GWAS hits validated in our study, three have confirmed or hypothesized roles in phospholipid signaling and lipid metabolism. rdgA encodes a diacylglycerol kinase that catalyzes the formation of phosphatidic acid, a phospholipid with pleiotropic signaling effects [60, 65]. Our discovery of rdgA as a regulator of body mass, starvation resistance, and lifespan dovetails with the remarkable pleiotropy of this gene, including roles in taste [66] and smell [67] as well as phototransduction [68], sound [69], temperature perception [70], and lifespan [60]. Under one plausible model, rdgA’s main link to body mass and nutrient response could be in the control of taste preference and/or feeding behavior. Alternatively, rdgA’s metabolic role could involve its known link to TOR signaling [60, 65, 71], plausibly via regulation of macroautophagy in the fat body [72] or global regulation of translation in response to stress [73]. The latter would be consistent with our finding that rdgA mutants are only long-lived when on a restricted diet. Another of our validated GWAS hit loci, the ceramide synthase schlank (homologous to the yeast lifespan gene LAG1 [47]), is known to control larval body size and fat accumulation [46]. Compromised lipid storage is therefore a potential proximal cause for the drop in body mass, starvation resistance, and lifespan we observe upon schlank knockdown in adult flies. The latter effects could involve the redirection of biosynthetic resources to other lipid pathways as ceramide and sphingolipid biosynthesis is blocked, since schlank loss of function is known to perturb levels of a number of lipid species [46]. However, it is tempting to speculate that schlank could act as a master metabolic regulator, in light of its ubiquitous expression and the evidence for ceramide as a mediator of apoptosis [74, 75] and MAPKKK/JNK mediated stress-response signaling [74, 76, 77] in the fly. A third gene on which we focus here, CG43921 (which we have named heavyweight), is of unknown biochemical function but is annotated with a FYVE domain which may bind phosphatidylinositol 3-phosphate (PI3P) [78], as well as a src homology 2 domain [79]. Future work will be necessary to establish whether rdgA, schlank, and hwt, and the cohort of additional phospholipid metabolic genes in which we detected enriched association signal (Additional file 1: Table S5), indeed exert effects on body mass and starvation resistance via phospholipid signaling. A particularly compelling model would invoke phospholipid regulation of the insulin pathway, to modulate apoptosis and cell growth [80] and/or glucose homeostasis [81–84].

Our results also reveal a novel metabolic function for fly huntingtin, reminiscent of the weight loss seen in human Huntington’s patients [50, 51] and in some mouse models of the disease [85–88]. The molecular-genetic origin of these phenotypes is a subject of ongoing debate, since in principle it could manifest from either a toxic gain-of-function in the mutant protein or a loss of wildtype huntingtin function [51, 86, 89]. As we observed dramatic increases in body size and starvation resistance in htt knockdown flies, our study provides a first compelling rationale for refuting the hypothesis that weight loss results from huntingtin loss of function. Our data leaves open the question of how native htt acts as a negative regulator of body mass and starvation resistance. Conceivably, its metabolic function could be mediated by circulating leptin and adiponectin [85, 89] and/or IGF-1 [89], which are perturbed by Huntington’s disease mutations in mammals. In the face of the early lethality of mouse huntingtin nulls [90] and the developmental defects in huntingtin knockdown zebrafish [91], expedient invertebrate models are sorely needed in the field. Our analyses of fly htt lay the groundwork for tests of mechanistic hypotheses centered around this key human disease gene.

Conclusions

We have conducted genome-wide scans for variants associated with starvation resistance, body mass and composition and their response to changing diets. We have validated the effects of schlank, rdgA, and htt on diet-dependent traits, and we have uncovered a panel of additional candidates for future analyses. Alongside these novel genetic discoveries, our experimental design sets a precedent for the use of gene-by-diet and cross-phenotype association tests in fruit fly GWAS.

Methods

DGRP phenotyping

DGRP lines were obtained from the Bloomington Drosophila Stock Center [92]. Fly husbandry and food were as previously described [36].

To measure starvation resistance for a given line, an isogenic population of 100 non-virgin female flies was reared on either DR (0.5 % yeast) or AL (5 % yeast) food from the day of eclosion, with 25 flies per vial. After 10 days, flies were switched to starvation media (agar providing water but containing no nutrients), and deaths of flies were counted three times daily until all flies died; we then tabulated the average time to death across this isogenic population for use in statistical genetics analyses below. One biological replicate (one batch of 100 animals) was assayed for 159 DGRP strains and two biological replicates (two batches of 100 animals assayed starting on two different days) for another 12 strains.

For body mass, glucose and triglyceride measurements for a given line, an isogenic population of 15 seven day old non-virgin female flies was weighed in groups of 3 flies and then homogenized in 200 μl PBS and frozen. In contrast to the constitutive RNAi experiments, for adult-induced schlank RNAi we sampled schlank flies at 14 days after eclosion to allow enough time to observe a change in body weight and body composition. This whole-fly homogenate was subjected to two technical replicates of both the StanBio Triglycerides LiquiColor® Test (Cat# 2100-225) and the Glucose Liqui-UV® Test (Cat# 1060-500). The measurements were normalized to body mass [36]. One biological replicate (one batch of 15 animals) was assayed for 88 DGRP strains and two biological replicates (two batches of 15 animals each assayed on each of 2 days) for another 93 strains.

We filtered out cases in which variation between replicates for a given DGRP strain was >20-fold greater than the median intrastrain variation (among all our measurements, only triglyceride assays on nine DGRP strains failed to meet this filter) and subjected all remaining measurements to genetic analysis as detailed below.

Heritability

Broad-sense heritabilities in Additional file 1: Table S2 were calculated for each trait as (σ2 total - σ2 intrastrain)/ σ2 total, using only the DGRP lines for which replicate measurements were carried out. For σ2 total for a given trait, we calculated the variance across all phenotype measurements for all such lines. Separately, for a given DGRP line we calculated the variance in a given trait across biological replicates; we then calculated the average of these variances across lines as our measure of σ2 intrastrain.

Genotype filtering

We used DGRP release 2 genotypes and FlyBase R5 coordinates for gene annotations [17, 93]. For the genome-wide scan of low-frequency variants, we used all DGRP polymorphic positions with minor allele frequency ≥5 % (1,932,643 total). For analysis of human disease orthologs, we subjected the set of FlyBase genes with human homologs [44] to a homolog search via DIOPT-DIST [94, 95] to find best-match human-Drosophila homologs, and then filtered for human genes with entries in the Online Mendelian Inheritance in Man database [42]. All DGRP polymorphisms with minor allele frequency ≥5 % (172,637 total) within 1 kb of a homologous gene model were used in these OMIM homolog association tests. We did not filter variants based on the coding status of the polymorphic locus. For the genome-wide scan of common variants, we used all DGRP markers with minor allele frequency ≥25 % (734,670 total).

Single-trait association tests

For a given trait and marker, we used the complete cohort of average phenotype measurements across all replicates, from flies of all DGRP lines reared on both diets, as input into an association test via ordinary least squares regression using the statsmodels module in Python [96]. The linear model was phenotype = β1 x genotype + β2 x diet + β3 x genotype x diet + intercept. Nominal p-values denoted as “genotype” in Table 1 report the probability that β1 = 0, and those denoted as “interaction” report the probability that β3 = 0. To avoid the potential for false positives at a given nominal cutoff owing to p-value inflation, we calculated false discovery rates via permutation as follows: for a given permutation i, we randomized phenotype values across DGRP lines, retaining the true diet assignment, and on this permuted data set we carried out single-trait association tests for each marker in turn as above. We counted the number of markers n i that scored above a given p-value threshold t. We tabulated the false discovery rate at t as the ratio between the average n i across ten permutations and the number of markers called at t in the real data.

Cross-phenotype association tests

We carried out cross-phenotype association tests for starvation resistance and body mass; starvation resistance and triglyceride levels; and starvation resistance and glucose levels as follows. We used the phenotype measurements for a given trait pair from flies of all DGRP lines reared on a given diet, and genotype at a given marker, as input into the R package SMAT [41]. False discovery rate estimates were calculated via permutation as above. Because triglyceride and glucose level measurements were each normalized by body mass for a given DGRP line, these phenotypes were not independent and therefore we did not subject any pair among them to cross-phenotype association testing.

Because our cross-phenotype association analysis involved testing over multiple trait pairs (Additional file 1: Table S4), we assessed the empirical false discovery rate (FDR) across the entire set of tests at a given nominal SMAT p-value threshold t as follows. We tabulated the number of loci n true with SMAT significance at or exceeding t in at least one trait pair. Next, a permuted assignment of genotype to strain was used as input into a SMAT calculation for each trait pair in turn, yielding the number of loci n perm with SMAT significance at or exceeding t in at least one trait pair for this permutation. We took the ratio between the average of n perm , across 10 permutations, and n true as the multiple testing-corrected FDR across trait pairs. At t = 1.27E-08, corresponding to the SMAT p-value attained in the real data for our least significant SMAT hit, rdgA (Table 1), this procedure resulted in an estimated FDR of 12 %. Therefore our SMAT hits have in the worst case a 12 % multiple-testing corrected FDR. In the best case, the top scoring polymorphism in schlank returned p = 4.44E-16, a level that has a 2 % multiple-testing corrected FDR.

Percent variance explained by mapped loci

For a given trait, we calculated the proportion of the variance across DGRP lines explained by a given mapped locus as the R2 from the model phenotype = β1 x genotype + intercept, fit via regression as above using measurements from all lines on both diets. The results serve as an upper bound on the proportion of variance explained, due to the Beavis effect [97, 98]. For a given locus mapped in cross-phenotype association tests, we carried out regression and calculated R2 as above for each of the two traits, and in Table 1 the greater of the two is reported.

Functional genomic analysis

We tested for enrichment of association signal among gene ontology (GO) biological processes as follows. Given the results from a single-trait or cross-phenotype association scan across markers as above, we first tabulated the single best-scoring marker for every gene tested. We then used the ranked list of these genes as input into GOglm, which tests for GO term enrichment in a ranked list without arbitrary thresholding into significant and insignificant genes, and corrects for gene length effects [37]. Only GO terms with >10 members with association results were considered for analysis. The resulting enrichment significance estimates were corrected for multiple testing with the Benjamini-Hochberg method. Annotation was from GOseq [99].

Gene knockdown by RNAi

For a given candidate gene, we used functional validation strategies that enable the study of a line in which the gene is perturbed as a comparison with an isogenic control of the same background. For this purpose, for each gene we acquired a strain harboring a UAS-driven RNAi construct [100] from VDRC [101], or from the TRiP resource at the Bloomington Drosophila Stock Center [92], and we also acquired lines expressing a constitutive GAL4 driver [102] or an RU486-inducible GeneSwitch (GS) GAL4-human progesterone receptor fusion protein driver [103]. Strains used are listed in Additional file 1: Table S6. Drivers were checked for activity by mating to UAS-GFP transgene and assaying for GFP fluorescence in the progeny. For each knockdown experiment, five males of the UAS-RNAi line of interest were crossed to 20 virgin females of a GAL4 line. As negative controls for the latter, five males lacking the RNAi construct (background-matched to the RNAi lines) were crossed to 20 females of the same GAL4 line. In the case of rdgA we compared a previously characterized transposon insertion mutant [61] to a background-matched control. For each experiment, non-virgin female progeny were collected as day-old adults and phenotyped as described below.

qRT-PCR analysis of RNAi target expression

To confirm RNAi-mediated knockdown in a given line in Additional file 1: Figure S2, total RNA was extracted from three frozen whole flies per condition via Zymo Quick RNA MiniPrep kit (R1054). For qRT-PCR, we used Superscript III Platinum SYBR Green One-Step qRT-PCT kit from Invitrogen (11736-051) and followed the manufacturer’s instructions with a Stratagene Mx3000P qPCR machine. We used a six point standard curve for each target gene to calibrate the relative abundance levels of the amplified mRNA in three technical replicate samples per condition. ROX dye was used to normalize for reaction volume loading.

Phenotyping of knockdown and mutant strains

We used knockdown and mutant strains to test candidate genes’ roles in body mass, starvation resistance and lifespan. Genotypes of control and experimental strains are listed in Additional file 1: Table S6. For body mass and starvation experiments using GS drivers, at day 1 of adulthood, flies were sorted onto food containing either the GS inducer RU486, or ethanol as a control, and body mass or starvation was assayed 10 days later as for the DGRP strains (see above). To measure body mass, cohorts of 15 flies reared on a given diet were frozen and weighed at day 7 after eclosion for htt and CG43921, and day 14 for schlank and rdgA. For body mass, significance in a given diet was tested by t-test, and interaction between the genetic perturbation and diet was tested by ANOVA. Starvation resistance assays were performed as on DGRP strains (see above); each trace in main and supplementary figures reports a representative experiment from among 2–5 biological replicates of 100 non-virgin female flies per condition. Lifespan assays used 200 non-virgin female flies per condition, 25 flies per vial, changing media every 2 days. For lifespan assays all knockdowns were induced only in adult flies. Each trace in Additional file 1: Figure S5 reports a representative experiment from among 1–5 biological replicates of 100–200 flies per condition. Survival curves for starvation resistance and lifespan were analyzed with the Cox proportional hazards method implemented in the R package “survival”.

Abbreviations

AL:

Ad libitum diet. 5 % yeast media

DGRP:

Drosophila Genetic Reference Panel

DR:

Dietary restriction. 0.5 % yeast media

FDR:

False discovery rate

GO:

Gene ontology

GS:

GeneSwitch RU486-inducible GAL4 transactivator

GWAS:

Genome-wide association study

RNAi:

RNA interference

UAS:

Upstream activation sequence

References

  1. 1.

    Jemal A. Trends in the leading causes of death in the United States, 1970-2002. JAMA. 2005;294:1255.

  2. 2.

    Williams KW, Elmquist JK. From neuroanatomy to behavior: central integration of peripheral signals regulating feeding behavior. Nat Neurosci. 2012;15:1350–5.

  3. 3.

    Fall T, Ingelsson E. Genome-wide association studies of obesity and metabolic syndrome. Mol Cell Endocrinol. 2014;382:740–57.

  4. 4.

    Drake TA, Schadt EE, Davis RC, Lusis AJ. Integrating genetic and gene expression data to study the metabolic syndrome and diabetes in mice. Am J Ther. 2005;12:503–11.

  5. 5.

    Parks BW, Nam E, Org E, Kostem E, Norheim F, Hui ST, et al. Genetic control of obesity and gut microbiota composition in response to high-fat, high-sucrose diet in mice. Cell Metab. 2013;17:141–52.

  6. 6.

    Lin C, Theodorides ML, McDaniel AH, Tordoff MG, Zhang Q, Li X, et al. QTL analysis of dietary obesity in C57BL/6byj X 129P3/J F2 mice: diet- and sex-dependent effects. PLoS One. 2013;8:e68776. Meyerhof W, editor.

  7. 7.

    Musselman LP, Fink JL, Narzinski K, Ramachandran PV, Hathiramani SS, Cagan RL, et al. A high-sugar diet produces obesity and insulin resistance in wild-type Drosophila. Model Mech. 2011;4:842–9.

  8. 8.

    Na J, Musselman LP, Pendse J, Baranski TJ, Bodmer R, Ocorr K, et al. A Drosophila model of high sugar diet-induced cardiomyopathy. PLoS Genet. 2013;9:e1003175.

  9. 9.

    Pendse J, Ramachandran PV, Na J, Narisu N, Fink JL, Cagan RL, et al. A Drosophila functional evaluation of candidates from human genome-wide association studies of type 2 diabetes and related metabolic traits identifies tissue-specific roles for dHHEX. BMC Genomics. 2013;14:136.

  10. 10.

    Teague BD, Clark AG, Doane WW. Developmental analysis of lipids from wild-type and adipose60 mutants of Drosophila melanogaster. J Exp Zool. 1986;240:95–104.

  11. 11.

    Hader T, Muller S, Aguilera M, Eulenberg KG, Steuernagel A, Ciossek T, et al. Control of triglyceride storage by a WD40/TPR-domain protein. EMBO Rep. 2003;4:511–6.

  12. 12.

    Grönke S, Beller M, Fellert S, Ramakrishnan H, Jäckle H, Kühnlein RP. Control of fat storage by a Drosophila PAT domain protein. Curr Biol CB. 2003;13:603–6.

  13. 13.

    Xu X, Gopalacharyulu P, Seppänen-Laakso T, Ruskeepää A-L, Aye CC, Carson BP, et al. Insulin signaling regulates fatty acid catabolism at the level of CoA activation. PLoS Genet. 2012;8:e1002478. Rulifson E, editor.

  14. 14.

    Karpac J, Biteau B, Jasper H. Misregulation of an adaptive metabolic response contributes to the age-related disruption of lipid homeostasis in Drosophila. Cell Rep. 2013;4:1250–61.

  15. 15.

    Birse RT, Choi J, Reardon K, Rodriguez J, Graham S, Diop S, et al. High-fat-diet-induced obesity and heart dysfunction are regulated by the TOR pathway in Drosophila. Cell Metab. 2010;12:533–44.

  16. 16.

    Mackay TF, Richards S, Stone EA, Barbadilla A, Ayroles JF, Zhu D, et al. The drosophila melanogaster genetic reference panel. Nature. 2012;482:173–8.

  17. 17.

    Huang W, Massouras A, Inoue Y, Peiffer J, Ramia M, Tarone AM, et al. Natural variation in genome architecture among 205 Drosophila melanogaster Genetic Reference Panel lines. Genome Res. 2014;24:1193–208.

  18. 18.

    Jumbo-Lucioni P, Ayroles JF, Chambers M, Jordan KW, Leips J, Mackay TF, et al. Systems genetics analysis of body weight and energy metabolism traits in Drosophila melanogaster. BMC Genomics. 2010;11:297.

  19. 19.

    Reed LK, Lee K, Zhang Z, Rashid L, Poe A, Hsieh B, et al. Systems genomics of metabolic phenotypes in wild-type Drosophila melanogaster. Genetics. 2014;197:781–93.

  20. 20.

    Dobson AJ, Chaston JM, Newell PD, Donahue L, Hermann SL, Sannino DR, et al. Host genetic determinants of microbiota-dependent nutrition revealed by genome-wide analysis of Drosophila melanogaster. Nat Commun. 2015;6:6312.

  21. 21.

    Unckless RL, Rottschaefer SM, Lazzaro BP. The complex contributions of genetics and nutrition to immunity in Drosophila melanogaster. PLoS Genet. 2015;11:e1005030. Cherry S, editor.

  22. 22.

    Vonesch SC, Lamparter D, Mackay TFC, Bergmann S, Hafen E. Genome-wide analysis reveals novel regulators of growth in Drosophila melanogaster. PLoS Genet. 2016;12:e1005616. Barsh GS, editor.

  23. 23.

    Galesloot TE, van Steen K, Kiemeney LALM, Janss LL, Vermeulen SH. A comparison of multivariate genome-wide association methods. PLoS One. 2014;9:e95923.

  24. 24.

    Unckless RL, Rottschaefer SM, Lazzaro BP. A genome-wide association study for nutritional indices in Drosophila. G3 Bethesda Md. 2015;5:417–25.

  25. 25.

    Harbison ST, Yamamoto AH, Fanara JJ, Norga KK, Mackay TFC. Quantitative trait loci affecting starvation resistance in Drosophila melanogaster. Genetics. 2004;166:1807–23.

  26. 26.

    Reed LK, Williams S, Springston M, Brown J, Freeman K, DesRoches CE, et al. Genotype-by-diet interactions drive metabolic phenotype variation in Drosophila melanogaster. Genetics. 2010;185:1009–19.

  27. 27.

    Zid BM, Rogers AN, Katewa SD, Vargas MA, Kolipinski MC, Lu TA, et al. 4E-BP extends lifespan upon dietary restriction by enhancing mitochondrial activity in Drosophila. Cell. 2009;139:149–60.

  28. 28.

    Kapahi P, Zid BM, Harper T, Koslover D, Sapin V, Benzer S. Regulation of lifespan in Drosophila by modulation of genes in the TOR signaling pathway. Curr Biol CB. 2004;14:885–90.

  29. 29.

    Mair W, Piper MDW, Partridge L. Calories do not explain extension of life span by dietary restriction in Drosophila. PLoS Biol. 2005;3:e223.

  30. 30.

    Vargas MA, Luo N, Yamaguchi A, Kapahi P. A role for S6 kinase and serotonin in postmating dietary switch and balance of nutrients in D. melanogaster. Curr Biol. 2010;20:1006–11.

  31. 31.

    Gendron CM, Kuo T-H, Harvanek ZM, Chung BY, Yew JY, Dierick HA, et al. Drosophila life span and physiology are modulated by sexual perception and reward. Science. 2014;343:544–8.

  32. 32.

    Burger JMS, Hwangbo DS, Corby-Harris V, Promislow DEL. The functional costs and benefits of dietary restriction in Drosophila. Aging Cell. 2007;6:63–71.

  33. 33.

    Morris SNS, Coogan C, Chamseddin K, Fernandez-Kim SO, Kolli S, Keller JN, et al. Development of diet-induced insulin resistance in adult Drosophila melanogaster. Biochim Biophys Acta BBA Mol Basis Dis. 1822;2012:1230–7.

  34. 34.

    Bross TG, Rogina B, Helfand SL. Behavioral, physical, and demographic changes in Drosophila populations through dietary restriction. Aging Cell. 2005;4:309–17.

  35. 35.

    Piper MDW, Skorupa D, Partridge L. Diet, metabolism and lifespan in Drosophila. Exp Gerontol. 2005;40:857–62.

  36. 36.

    Katewa SD, Demontis F, Kolipinski M, Hubbard A, Gill MS, Perrimon N, et al. Intramyocellular fatty-acid metabolism plays a critical role in mediating responses to dietary restriction in Drosophila melanogaster. Cell Metab. 2012;16:97–103.

  37. 37.

    Mi G, Di Y, Emerson S, Cumbie JS, Chang JH. Length bias correction in gene ontology enrichment analysis using logistic regression. PLoS One. 2012;7:e46128. Steinke D, editor.

  38. 38.

    Lee G, Park JH. Hemolymph sugar homeostasis and starvation-induced hyperactivity affected by genetic manipulations of the adipokinetic hormone-encoding gene in Drosophila melanogaster. Genetics. 2004;167:311–23.

  39. 39.

    Meunier N, Belgacem YH, Martin J-R. Regulation of feeding behaviour and locomotor activity by takeout in Drosophila. J Exp Biol. 2007;210:1424–34.

  40. 40.

    Johnson EC, Kazgan N, Bretz CA, Forsberg LJ, Hector CE, Worthen RJ, et al. Altered metabolism and persistent starvation behaviors caused by reduced AMPK function in Drosophila. PLoS One. 2010;5:e12799. Hassan BA, editor.

  41. 41.

    Schifano ED, Li L, Christiani DC, Lin X. Genome-wide association analysis for multiple continuous secondary phenotypes. Am J Hum Genet. 2013;92:744–59.

  42. 42.

    OMIM - Online Mendelian Inheritance in Man. McKusick-Nathans Institute of Genetic Medicine, Johns Hopkins University, Baltimore. 2015. http://omim.org/. Accessed 14 Sept 2015.

  43. 43.

    Howard A, Rogers AN. Role of translation initiation factor 4G in lifespan regulation and age-related health. Ageing Res Rev. 2014;13:115–24.

  44. 44.

    dos Santos G, Schroeder AJ, Goodman JL, Strelets VB, Crosby MA, Thurmond J, et al. FlyBase: introduction of the Drosophila melanogaster Release 6 reference genome assembly and large-scale migration of genome annotations. Nucleic Acids Res. 2015;43:D690–7.

  45. 45.

    Pesch Y-Y, Riedel D, Patil KR, Loch G, Behr M. Chitinases and Imaginal disc growth factors organize the extracellular matrix formation at barrier tissues in insects. Sci Rep. 2016;6:18340.

  46. 46.

    Bauer R, Voelzmann A, Breiden B, Schepers U, Farwanah H, Hahn I, et al. Schlank, a member of the ceramide synthase family controls growth and body fat in Drosophila. EMBO J. 2009;28:3706–16.

  47. 47.

    D’mello NP, Childress AM, Franklin DS, Kale SP, Pinswasdi C, Jazwinski SM. Cloning and characterization of LAG1, a longevity-assurance gene in yeast. J Biol Chem. 1994;269:15451–9.

  48. 48.

    Tricoire H, Battisti V, Trannoy S, Lasbleiz C, Pret A-M, Monnier V. The steroid hormone receptor EcR finely modulates Drosophila lifespan during adulthood in a sex-specific manner. Mech Ageing Dev. 2009;130:547–52.

  49. 49.

    Perkins LA, Holderbaum L, Tao R, Hu Y, Sopko R, McCall K, et al. The transgenic RNAi project at Harvard Medical School: resources and validation. Genetics. 2015;201:843–52.

  50. 50.

    Goodman AOG, Murgatroyd PR, Medina-Gomez G, Wood NI, Finer N, Vidal-Puig AJ, et al. The metabolic profile of early Huntington’s disease- a combined human and transgenic mouse study. Exp Neurol. 2008;210:691–8.

  51. 51.

    Aziz NA, van der Marck MA, Pijl H, Olde Rikkert MGM, Bloem BR, Roos RAC. Weight loss in neurodegenerative disorders. J Neurol. 2008;255:1872–80.

  52. 52.

    Zallen JA, Cohen Y, Hudson AM, Cooley L, Wieschaus E, Schejter ED. SCAR is a primary regulator of Arp2/3-dependent morphological events in Drosophila. J Cell Biol. 2002;156:689–701.

  53. 53.

    Koch N, Kobler O, Thomas U, Qualmann B, Kessels MM. Terminal axonal arborization and synaptic bouton formation critically rely on abp1 and the arp2/3 complex. PLoS One. 2014;9:e97692.

  54. 54.

    Crowner D, Le Gall M, Gates MA, Giniger E. Notch steers Drosophila ISNb motor axons by regulating the Abl signaling pathway. Curr Biol CB. 2003;13:967–72.

  55. 55.

    Pipes GC, Lin Q, Riley SE, Goodman CS. The Beat generation: a multigene family encoding IgSF proteins related to the Beat axon guidance molecule in Drosophila. Dev Camb Engl. 2001;128:4545–52.

  56. 56.

    Chen S, Spletter M, Ni X, White KP, Luo L, Long M. Frequent recent origination of brain genes shaped the evolution of foraging behavior in Drosophila. Cell Rep. 2012;1:118–32.

  57. 57.

    Seong KH, Ogashiwa T, Matsuo T, Fuyama Y, Aigaki T. Application of the gene search system to screen for longevity genes in Drosophila. Biogerontology. 2001;2:209–17.

  58. 58.

    Chibalin AV, Leng Y, Vieira E, Krook A, Björnholm M, Long YC, et al. Downregulation of diacylglycerol kinase delta contributes to hyperglycemia-induced insulin resistance. Cell. 2008;132:375–86.

  59. 59.

    Samuel VT, Petersen KF, Shulman GI. Lipid-induced insulin resistance: unravelling the mechanism. Lancet. 2010;375:2267–77.

  60. 60.

    Lin Y-H, Chen Y-C, Kao T-Y, Lin Y-C, Hsu T-E, Wu Y-C, et al. Diacylglycerol lipase regulates lifespan and oxidative stress response by inversely modulating TOR signaling in Drosophila and C. elegans. Aging Cell. 2014;13:755–64.

  61. 61.

    Bellen HJ, Levis RW, Liao G, He Y, Carlson JW, Tsang G, et al. The BDGP gene disruption project: single transposon insertions associated with 40 % of Drosophila genes. Genetics. 2004;167:761–81.

  62. 62.

    Solon-Biet SM, McMahon AC, Ballard JWO, Ruohonen K, Wu LE, Cogger VC, et al. The ratio of macronutrients, not caloric intake, dictates cardiometabolic health, aging, and longevity in Ad Libitum-Fed Mice. Cell Metab. 2014;19:418–30.

  63. 63.

    Levine ME, Suarez JA, Brandhorst S, Balasubramanian P, Cheng C-W, Madia F, et al. Low protein intake is associated with a major reduction in IGF-1, cancer, and overall mortality in the 65 and younger but not older population. Cell Metab. 2014;19:407–17.

  64. 64.

    Wang K, Li M, Hakonarson H. Analysing biological pathways in genome-wide association studies. Nat Rev Genet. 2010;11:843–54.

  65. 65.

    Yoon M-S, Sun Y, Arauz E, Jiang Y, Chen J. Phosphatidic acid activates mammalian target of rapamycin complex 1 (mTORC1) kinase by displacing FK506 binding protein 38 (FKBP38) and exerting an allosteric effect. J Biol Chem. 2011;286:29568–74.

  66. 66.

    Kain P, Badsha F, Hussain SM, Nair A, Hasan G, Rodrigues V. Mutants in phospholipid signaling attenuate the behavioral response of adult Drosophila to trehalose. Chem Senses. 2010;35:663–73.

  67. 67.

    Kain P, Chakraborty TS, Sundaram S, Siddiqi O, Rodrigues V, Hasan G. Reduced odor responses from antennal neurons of Gq, phospholipase C, and rdgA mutants in Drosophila support a role for a phospholipid intermediate in insect olfactory transduction. J Neurosci. 2008;28:4745–55.

  68. 68.

    Raghu P, Usher K, Jonas S, Chyb S, Polyanovsky A, Hardie RC. Constitutive activity of the light-sensitive channels TRP and TRPL in the Drosophila diacylglycerol kinase mutant, rdgA. Neuron. 2000;26:169–79.

  69. 69.

    Senthilan PR, Piepenbrock D, Ovezmyradov G, Nadrowski B, Bechstedt S, Pauls S, et al. Drosophila auditory organ genes and genetic hearing defects. Cell. 2012;150:1042–54.

  70. 70.

    Kwon Y, Shim H-S, Wang X, Montell C. Control of thermotactic behavior via coupling of a TRP channel to a phospholipase C signaling cascade. Nat Neurosci. 2008;11:871–3.

  71. 71.

    Yoon M-S, Rosenberger CL, Wu C, Truong N, Sweedler JV, Chen J. Rapid mitogenic regulation of the mTORC1 inhibitor, DEPTOR, by phosphatidic acid. Mol Cell. 2015;58:549–56.

  72. 72.

    Scott RC, Schuldiner O, Neufeld TP. Role and regulation of starvation-induced autophagy in the Drosophila fat body. Dev Cell. 2004;7:167–78.

  73. 73.

    Teleman AA. 4E-BP functions as a metabolic brake used under stress conditions but not during normal growth. Genes Dev. 2005;19:1844–8.

  74. 74.

    Sathyanarayana P, Barthwal MK, Kundu CN, Lane ME, Bergmann A, Tzivion G, et al. Activation of the Drosophila MLK by ceramide reveals TNF-alpha and ceramide as agonists of mammalian MLK3. Mol Cell. 2002;10:1527–33.

  75. 75.

    Pronk GJ, Ramer K, Amiri P, Williams LT. Requirement of an ICE-like protease for induction of apoptosis and ceramide generation by REAPER. Science. 1996;271:808–10.

  76. 76.

    Kraut R. Roles of sphingolipids in Drosophila development and disease: Roles of sphingolipids in Drosophila. J Neurochem. 2011;116:764–78.

  77. 77.

    Botella JA, Baines IA, Williams DD, Goberdhan DC, Proud CG, Wilson C. The Drosophila cell shape regulator c-Jun N-terminal kinase also functions as a stress-activated protein kinase. Insect Biochem Mol Biol. 2001;31:839–47.

  78. 78.

    Gaullier JM, Simonsen A, D’Arrigo A, Bremnes B, Stenmark H, Aasland R. FYVE fingers bind PtdIns(3)P. Nature. 1998;394:432–3.

  79. 79.

    Huang H, Li L, Wu C, Schibli D, Colwill K, Ma S, et al. Defining the specificity space of the human Src homology 2 domain. Mol Cell Proteomics. 2007;7:768–84.

  80. 80.

    Hietakangas V, Cohen SM. Re-evaluating AKT regulation: role of TOR complex 2 in tissue growth. Genes Dev. 2007;21:632–7.

  81. 81.

    Alb JG, Cortese JD, Phillips SE, Albin RL, Nagy TR, Hamilton BA, et al. Mice lacking phosphatidylinositol transfer protein- exhibit spinocerebellar degeneration, intestinal and hepatic steatosis, and hypoglycemia. J Biol Chem. 2003;278:33501–18.

  82. 82.

    Yea K, Kim J, Yoon JH, Kwon T, Kim JH, Lee BD, et al. Lysophosphatidylcholine activates adipocyte glucose uptake and lowers blood glucose levels in Murine models of diabetes. J Biol Chem. 2009;284:33833–40.

  83. 83.

    Rancoule C, Attané C, Grès S, Fournel A, Dusaulcy R, Bertrand C, et al. Lysophosphatidic acid impairs glucose homeostasis and inhibits insulin secretion in high-fat diet obese mice. Diabetologia. 2013;56:1394–402.

  84. 84.

    Ersoy BA, Tarun A, D’Aquino K, Hancer NJ, Ukomadu C, White MF, et al. Phosphatidylcholine transfer protein interacts with Thioesterase Superfamily Member 2 to attenuate insulin signaling. Sci Signal. 2013;6:ra64.

  85. 85.

    Phan J, Hickey MA, Zhang P, Chesselet M-F, Reue K. Adipose tissue dysfunction tracks disease progression in two Huntington’s disease mouse models. Hum Mol Genet. 2009;18:1006–16.

  86. 86.

    Aziz NA, van der Burg JMM, Landwehrmeyer GB, Brundin P, Stijnen T, EHDI Study Group, et al. Weight loss in Huntington disease increases with higher CAG repeat number. Neurology. 2008;71:1506–13.

  87. 87.

    Petersén Å, Björkqvist M. Hypothalamic-endocrine aspects in Huntington’s disease. Eur J Neurosci. 2006;24:961–7.

  88. 88.

    Mochel F, Haller RG. Energy deficit in Huntington disease: why it matters. J Clin Invest. 2011;121:493–9.

  89. 89.

    Pouladi MA, Xie Y, Skotte NH, Ehrnhoefer DE, Graham RK, Kim JE, et al. Full-length huntingtin levels modulate body weight by influencing insulin-like growth factor 1 expression. Hum Mol Genet. 2010;19:1528–38.

  90. 90.

    Nasir J, Floresco SB, O’Kusky JR, Diewert VM, Richman JM, Zeisler J, et al. Targeted disruption of the Huntington’s disease gene results in embryonic lethality and behavioral and morphological changes in heterozygotes. Cell. 1995;81:811–23.

  91. 91.

    Lumsden AL, Henshall TL, Dayan S, Lardelli MT, Richards RI. Huntingtin-deficient zebrafish exhibit defects in iron utilization and development. Hum Mol Genet. 2007;16:1905–20.

  92. 92.

    BDSC. Indiana University, Bloomington. 2015. http://flystocks.bio.indiana.edu/. Accessed 14 Sept 2015.

  93. 93.

    MacKay Lab. Genotype Files. 2014. http://dgrp2.gnets.ncsu.edu/data.html. Accessed 3 Sept 2014.

  94. 94.

    Hu Y, Flockhart I, Vinayagam A, Bergwitz C, Berger B, Perrimon N, et al. An integrative approach to ortholog prediction for disease-focused and other functional studies. BMC Bioinformatics. 2011;12:357.

  95. 95.

    Hu et al. Disease Gene Query. 2015. http://www.flyrnai.org/cgi-bin/DRSC_DG_query.pl. Accessed 14 Sep 2015.

  96. 96.

    Perktold et al. StatsModels: Statistics in Python — statsmodels 0.6.1. 2015. http://statsmodels.sourceforge.net/stable/. Accessed 15 Sep 2015. 

  97. 97.

    Beavis W. Molecular Dissection of Complex Traits. QTL Anal. Power Precis. Accuracy. New York: CRC Press; 1998. p. 145–62.

  98. 98.

    Garlapow ME, Huang W, Yarboro MT, Peterson KR, Mackay TFC. Quantitative Genetics of Food Intake in Drosophila melanogaster. PLoS One. 2015;10:e0138129. Ko DC, editor.

  99. 99.

    Young MD, Wakefield MJ, Smyth GK, Oshlack A. Gene ontology analysis for RNA-seq: accounting for selection bias. Genome Biol. 2010;11:R14.

  100. 100.

    Dietzl G, Chen D, Schnorrer F, Su K-C, Barinova Y, Fellner M, et al. A genome-wide transgenic RNAi library for conditional gene inactivation in Drosophila. Nature. 2007;448:151–6.

  101. 101.

    Vienna Biocenter Core Facilities GmbH. VDRC Stock Center: Main Page. 2015. http://stockcenter.vdrc.at/control/maincited. Accessed 14 Sep 2015.

  102. 102.

    Brand AH, Perrimon N. Targeted gene expression as a means of altering cell fates and generating dominant phenotypes. Dev Camb Engl. 1993;118:401–15.

  103. 103.

    Osterwalder T, Yoon KS, White BH, Keshishian H. A conditional tissue-specific transgene expression system using inducible GAL4. Proc. Natl. Acad. Sci. U. S. A. 2001;98:12596-601.

Download references

Acknowledgements

This work was supported by NIH/NIA grants AG049494, AG045835 and AG038688 (to PK) and a Google Research Cloud Credits Award. CN was supported by NIH/NIA F32AG047024. Stocks obtained from the Bloomington Drosophila Stock Center (NIH P40OD018537) were used in this study. The Brem lab used a LightMiner 36-core server for computations. This machine was a generous gift of Sysorex Inc.

Availability of data and materials

Phenotype data used in GWAS (shown in Additional file 1: Figure S1) is available upon request.

Authors’ contributions

Conceived and designed experiments: CSN, PK, RBB. Performed the experiments: CSN, JNB, KAW. Analyzed the data: CSN, JN, KAW, ERP. Wrote the manuscript: CSN, PK, RBB. All authors have read and approved the manuscript.

Competing interests

The authors declare that they have no competing interests.

Ethics approval and consent to participate

No ethics review committee was consulted for these studies, as there is no local review for non-mammalian studies.

Author information

Author notes

    Affiliations

    1. Buck Institute for Research on Aging, 8001 Redwood Blvd., Novato, CA, 94947, USA
      • Christopher S. Nelson
      • , Jennifer N. Beck
      • , Kenneth A. Wilson
      • , Elijah R. Pilcher
      • , Pankaj Kapahi
      •  & Rachel B. Brem
    2. Davis School of Gerontology, University of Southern California, Los Angeles, CA, USA
      • Kenneth A. Wilson
      • , Pankaj Kapahi
      •  & Rachel B. Brem
    3. Department of Urology, University of California, San Francisco, CA, USA
      • Jennifer N. Beck
      •  & Pankaj Kapahi
    4. Department of Plant and Microbial Biology, University of California, Berkeley, Berkeley, CA, USA
      • Rachel B. Brem
    Authors
    1. Search for Christopher S. Nelson in:
    2. Search for Jennifer N. Beck in:
    3. Search for Kenneth A. Wilson in:
    4. Search for Elijah R. Pilcher in:
    5. Search for Pankaj Kapahi in:
    6. Search for Rachel B. Brem in:

    Corresponding authors

    Correspondence to Pankaj Kapahi or Rachel B. Brem.

    Additional file

    Additional file 1: Table S1.

    Correlation between phenotypes in DGRP lines raised on two different diets. Table S2. Broad-sense heritabilities of metabolic traits. Table S3. GO terms enriched for association signal in scans of variants with allele frequency ≥5 %. Table S4. Phenotype pairs tested via SMAT. Table S5. GO terms enriched for association signal in scans of variants with allele frequency ≥25 %. Table S6. Strains and crosses used. Figure S1. Variation in metabolic traits and nutrient response across the DGRP. In a given bar chart, each bar reports the mean of the indicated phenotype in an isogenic population of one DGRP line reared on the indicated diet. Strains are ordered along the x-axis by AL diet phenotype. The right hand panels report the same data formatted as a scatterplot. (A), Resistance to acute starvation after rearing adults for 10 days on the indicated diet. (B), Body mass. (C), Whole body triglyceride (TG) content, normalized by body mass. (D), Whole body glucose content, normalized by body mass. a.u., dimensionless arbitrary units. AL, ad libitum diet; DR, dietary restriction. Figure S2. Decreased expression of schlank, htt, hwt, and rdgA upon RNAi. In the first, second, fourth, and sixth panels, each bar reports qRT-PCR measurements of expression of the indicated gene, in a line reared on the indicated diet and harboring both the GeneSwitch (GS) activator and a construct for RNAi of the indicated gene under the control of the indicated driver (da, daughterless; elav, embryonic lethal, abnormal vision), treated with the GS inducer RU486 (+) or a vehicle control (-); in a given diet, each expression measurement is normalized to that from the respective control-treated animals. In the third and fifth panels, each bar reports qRT-PCR measurements of expression of the indicated gene in a strain expressing an RNAi construct for the indicated gene regulated by GAL4 under the daughterless promoter (RNAi) or a background-matched control (Ctrl) with no RNAi construct; in a given diet, each expression measurement is normalized to that from the respective control animals. Figure S3. The GeneSwitch inducer RU486 has no effect on starvation resistance or body mass. (A), Each trace reports survival of flies expressing daughterless promoter (da)-driven GAL4-Gene Switch (GS) that do not harbor an RNAi construct, with drug treatment (dotted lines) or ethanol vehicle treatment (solid lines) starting from eclosion (n = 2 biological replicates of 100 flies per condition). (B), Each column reports the distribution of body mass per fly of the no-RNAi control strain harboring the da-GS driver, with genotypes and treatments as in (A) (n = 2 populations of 15 flies per condition). The checkered boxes indicate addition of the RU486 inducer drug. Top, middle, and bottom horizontal bars of a given vertical box denote the respective quartiles across batches and technical replicates, and the top and bottom short horizontal bars report minimum and maximum, respectively. No changes here were significant by t-test at alpha = 0.05. AL, ad libitum diet; DR, dietary restriction. For strain details see Additional file 1: Table S6. Figure S4. The GeneSwitch inducer RU486 has no effect on triglyceride levels in the absence of an RNAi construct, and htt and hwt knockdown does not detectably impact triglycerides. (A), Each bar reports whole-body triglyceride levels in a strain expressing a ubiquitous da-GeneSwitch driver and no RNAi transgene, treated with the GeneSwitch inducer RU486 (“+ drug”) or with ethanol alone (“- drug”). (B) and (C), Each bar reports whole-body triglyceride levels in a strain expressing constitutive da-GAL4 either with an RNAi construct for the indicated gene (“RNAi”) or without any RNAi construct for the indicated gene (“control”) (n = 1 population of 15 flies per condition). For strain details see Additional file 1: Table S6. AL, animals reared on the ad libitum diet; DR, dietary restriction. Horizontal lines in each boxplot indicate quartiles across technical replicates, and the top and bottom short horizontal bars report minimum and maximum, respectively. No changes here were significant by t-test at alpha = 0.05. Figure S5. Lifespan effects of schlank, htt, hwt and rdgA knockdown. In each panel, the right-hand plot reports survival of flies harboring an RNAi construct for one gene regulated by Act5C-driven GAL4-GeneSwitch treated with the GeneSwitch inducer RU486 (“+ drug”) or with ethanol alone (“no drug”), and the left-hand plot is for a background-matched control without an RNAi transgene. Drug or vehicle treatment was started at eclosion. (A), schlank RNAi. (B), htt RNAi. (C), hwt RNAi. (D), rdgA RNAi. AL, ad libitum diet; DR, dietary restriction. *, p < 0.05, **, p < 0.01, ***, p < 10-3, ****, p <10-7, *****, p < 10-15. Red asterisks denote significance of the effect of the genetic perturbation in animals reared on AL food, blue denotes significance in animals on DR food, and purple denotes significance of the interaction between diet and genetic perturbation. For strain details see Additional file 1: Table S6. Figure S6. hwt knockdown modestly affects resistance to acute starvation. Each trace reports survival, under acute starvation, of flies harboring a GAL4 regulator under a daughterless promoter and either a GAL4-regulated RNAi construct for hwt (dotted lines) or a matched control without a transgene (solid line). AL, animals reared before starvation treatment on the ad libitum diet; DR, dietary restriction. For strain details see Additional file 1: Table S6. *, p < 0.05, **, p < 0.01. (DOCX 2 mb)

    Rights and permissions

    Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

    Reprints and Permissions

    About this article

    Verify currency and authenticity via CrossMark

    Cite this article

    Nelson, C.S., Beck, J.N., Wilson, K.A. et al. Cross-phenotype association tests uncover genes mediating nutrient response in Drosophila . BMC Genomics 17, 867 (2016) doi:10.1186/s12864-016-3137-9

    Download citation

    Keywords

    • GWAS
    • Drosophila
    • Nutrition
    • Obesity
    • Dietary restriction

    BMC Genomics

    ISSN: 1471-2164

    Contact us