Skip to main content


Quantitative genomics of starvation stress resistance in Drosophila



A major challenge of modern biology is to understand the networks of interacting genes regulating complex traits, and the subset of these genes that affect naturally occurring quantitative genetic variation. Previously, we used P-element mutagenesis and quantitative trait locus (QTL) mapping in Drosophila to identify candidate genes affecting resistance to starvation stress, and variation in resistance to starvation stress between the Oregon-R (Ore) and 2b strains. Here, we tested the efficacy of whole-genome transcriptional profiling for identifying genes affecting starvation stress resistance.


We evaluated whole-genome transcript abundance for males and females of Ore, 2b, and four recombinant inbred lines derived from them, under control and starved conditions. There were significant differences in transcript abundance between the sexes for nearly 50% of the genome, while the transcriptional response to starvation stress involved approximately 25% of the genome. Nearly 50% of P-element insertions in 160 genes with altered transcript abundance during starvation stress had mutational effects on starvation tolerance. Approximately 5% of the genome exhibited genetic variation in transcript abundance, which was largely attributable to regulation by unlinked genes. Genes exhibiting variation in transcript abundance among lines did not cluster within starvation resistance QTLs, and none of the candidate genes affecting variation in starvation resistance between Ore and 2b exhibited significant differences in transcript abundance between lines.


Expression profiling is a powerful method for identifying networks of pleiotropic genes regulating complex traits, but the relationship between variation in transcript abundance among lines used to map QTLs and genes affecting variation in quantitative traits is complicated.


Quantitative traits affecting morphology, physiology, behavior, disease susceptibility and reproductive fitness are controlled by multiple interacting genes whose effects are conditional on the genetic, sexual and external environments [1]. Advances in medicine, agriculture, and an understanding of adaptive evolution depend on discovering the genes that regulate these complex traits, and determining the genetic and molecular properties of alleles at loci that cause segregating genetic variation in natural populations. Assessing subtle effects of induced mutations on quantitative trait phenotypes in model organisms is a straightforward approach to identify genes regulating complex traits [13]. However, the large number of potential mutations to evaluate, the necessity to induce mutations in a common inbred background, and the level of replication required to detect subtle effects [1] all limit the feasibility of systematic whole-genome mutagenesis screens for complex traits in higher eukaryotes. Mapping quantitative trait loci (QTLs) affecting variation in complex traits to broad genomic regions by linkage to polymorphic molecular markers is also straightforward. However, our ability to determine what genes in the QTL regions cause the trait variation is hampered by the large number of recombinants required for high-resolution mapping, and the small and environmentally sensitive effects of QTL alleles [1, 4].

There has been great excitement recently about the utility of whole-genome transcriptional profiling to identify candidate genes regulating complex traits, by assessing changes in gene expression in the background of single mutations affecting the trait [5, 6], between lines selected for different phenotypic values of the trait [7], and in response to environmental stress and aging [812]. Transcript abundance is also a quantitative trait for which there is considerable variation between wild-type strains [11, 1317], and for which expression QTLs (eQTLs) [18] have been mapped [1517, 19]. Thus, candidate genes affecting variation in quantitative trait phenotypes are those for which the map positions of trait QTL and eQTL coincide [16, 20].

Transcript profiling typically implicates hundreds to thousands of genes in the regulation of quantitative traits and associated with trait variation between strains; the majority of these genes are computationally predicted genes that have not been experimentally verified. To what extent do changes in transcript abundance predicate effects of induced mutations and allelic variants between strains on quantitative trait phenotypes? It is encouraging that several studies have confirmed the phenotypic effects of mutations in genes implicated by changes in expression [57]. However, limited numbers of genes were tested, and their choice was not unbiased. None of the candidate QTLs nominated by transcriptional profiling has been validated according to the rigorous standards necessary to prove that any candidate gene corresponds to a QTL [1, 4]. To begin to answer this question, we need to compare gene-expression data with genes known to affect the trait from independent mutagenesis and QTL mapping studies. This comparison has not been possible to date because there are only a few complex traits for which the genetic architecture is known at this level of detail, one of which is resistance to starvation stress in Drosophila.

Previously, we used P-element mutagenesis in an isogenic background to identify 383 candidate genes affecting starvation tolerance in D. melanogaster [21]. Further, we mapped QTLs affecting variation in starvation resistance between two isogenic Drosophila strains, Oregon-R (Ore) and 2b [21], followed by complementation tests to mutations to identify twelve candidate genes affecting variation in starvation resistance between these strains [21]. Here, we used Affymetrix Drosophila GeneChips to examine expression profiles of two starvation-resistant and two starvation-sensitive recombinant inbred (RI) lines, as well as parental lines Ore and 2b, under normal and starvation stress conditions. We used a statistically rigorous analysis to identify genes whose expression was altered between the sexes, during starvation stress treatment, between lines, and interactions between these main effects. In the comparison of expression profiling with the P-element mutagenesis performed previously, we found nearly 50% concordance between the effects of 160 P-element mutations on starvation stress resistance and changes in gene expression during starvation - 77 mutations with significant effects also had significant changes in transcript abundance, while 83 mutations did not affect the starvation resistance phenotype, yet had significant changes in transcript level. We identified 153 novel candidate genes for which there was variation in gene expression between the lines and which co-localized with starvation resistance QTLs. However, we did not detect genetic variation in expression for any of the candidate genes identified by complementation tests. Our efforts to associate genetic variation in expression with variation in quantitative trait phenotypes is confounded by the observation of widespread regulation of transcript abundance by unlinked genes, the difficulty in detecting rare transcripts that may be expressed in only a few cell types at a particular period of development, and genetic variation between QTL alleles that is not regulated at the level of transcription.


The sexually dimorphic transcriptome

Nearly one-half of the genome (6,569 probe sets) exhibited significantly different transcript levels between the sexes (P(Sex) < 0.001), with 3,965 probe sets upregulated in females and 2,604 probe sets upregulated in males (the complete list is given in Additional data file 1). The greatest differences in transcript abundance between the sexes were for probe sets implicated in sex-specific functions: chorion, vitelline membrane, and yolk proteins involved in egg production were upregulated in females; and accessory gland peptides, male-specific RNAs, and protein ejaculatory bulb components were upregulated in males. However, the probe sets exhibiting sex dimorphism in expression fell into 28 biological process and 41 molecular function Gene Ontology (GO) categories; for most of these categories, differences in expression between the sexes was unexpected. We determined which GO categories contained significantly different numbers of upregulated probe sets in males and females (Table 1). Genes involved in the biological process categories of cell communication, cell growth and/or maintenance, development, and cell death were upregulated more often in females than in males. Genes involved in the molecular function categories of binding, most enzymes, signal transduction, structural molecules, and regulation of transcription and translation were upregulated in females more often than in males; however, genes encoding oxidoreductase enzymes, carrier transporters and ion transporters were upregulated in males more often than in females (Table 1).

Table 1 Gene Ontology categories with sex-biased gene expression

The genomic distribution of sex-biased genes was not random (Figure 1). There was a paucity of male-biased genes on the X and fourth chromosomes, and an excess on chromosome 2R25 = 100.77; P < 0.0001). There was a deficit of female-biased genes on chromosome 4, and an excess on chromosome 2R25 = 29.18; P < 0.0001).

Figure 1

Chromosome locations of genes differentially expressed by sex. (a) Observed (magenta) and expected (blue) number of probe sets upregulated in males. (b) Observed (magenta) and expected (blue) numbers of probe sets upregulated in females.

Transcriptional response to starvation stress

We found 3,451 probe sets with significantly different mean transcript levels between the control and starved conditions (P(treatment) < 0.001): 1,736 were downregulated (some by as much as 40-fold) and 1,715 were upregulated (at most by 7.2-fold) during starvation (the complete list is available as Additional data file 2). These probe sets fell into 24 biological process and 25 molecular function GO categories. We determined which GO categories had a significantly different number of up- and downregulated probe sets in response to starvation stress. Genes affecting the biological processes of protein and nucleic-acid metabolism (protein biosynthesis; protein catabolism, folding, localization, modification, and repair; biosynthesis of nucleic acid macromolecules and lipids) were upregulated during starvation (Table 2). The expression of genes in three molecular function categories (nucleotide binding, hydrolases binding to acid anhydrides, and ribosome structure) increased during starvation; while defense/immunity proteins, peptidases, cuticle structural proteins, and carrier transport proteins were downregulated (Table 2).

Table 2 Gene Ontology categories with increased or decreased gene expression during starvation

The treatment × sex interaction term was significant (P < 0.001) for 817 probe sets, of which 715 had significant treatment effects for one or both sexes in the separate sex analyses (Additional data file 3). We categorized these 715 probe sets as sex-specific if significant expression changes in response to starvation occurred in one sex only; as sex-biased if expression levels changed in the same direction in both sexes, but were of different magnitude; or as sex-antagonistic if expression levels significantly changed in both sexes, but in opposite directions (Figure 2a-c). Most probe sets exhibited sex-specific or sex-biased expression, with only two genes, CG14095 and Rpd3, meeting the sex-antagonistic criterion. More probe sets exhibiting sex-specific or sex-biased expression were downregulated (454) than upregulated (263) during starvation. Starvation stress was accompanied by reduced expression of genes involved in the developmental processes of gametogenesis and sex determination as well as signal transduction in females, and of genes involved in mechanosensory and reproductive behavior in males (Table 2).

Figure 2

Genetic architecture of transcription. (a-c) Sex × treatment interaction for females (magenta)and males (blue): (a) Chorion protein 38; (b) Alkaline phosphatase 4; (c) Phosphogluconate dehydrogenase. (d-k) Interactions with line. Ore (black), 2b (red), RI 14 (green), RI 21 (dark blue), RI 35(magenta), RI 42 (light blue). (d, e) Sex × line interaction, averaged over treatments: (d) modulo; (e) l(2) giant larvae. (f-i) line × treatment interaction, averaged over sex: (f) CG11089; (g) Nervana 1; (h) Cyp9b2; (i) Peroxiredoxin 2540. (j, k) Sex × line × treatment interaction. The difference in expression between the starved and control treatments is plotted for females (magenta) and males (blue): (j) sallimus; (k) Esterase 6. (l-o) Regulation of transcript abundance. The same letters denote expression levels that are not significantly different. Magenta indicates 2b and blue indicates Ore genome. (l, m) Linked regulation of variation in transcript abundance: (l) UDP-glycosyltransferase 35b; (m) Signal recognition particle receptor b. (n, o) Unlinked regulation of variation in transcript abundance: (n) Arrestin 2; (o) Klarsicht.

Transcript abundance versus mutations

The genes represented by probe sets with significant treatment and/or treatment × sex effects are candidate genes for starvation resistance. Previously, we screened 933 co-isogenic single P-element insertion lines for their effect on starvation resistance [21]. Of these insertions, 383 had significant effects on starvation resistance, while the remaining 550 did not [21]. Of the 933 lines, we know the locations of the 385 of the inserts and that genes tagged by these inserts are represented on the array. Thus, we can directly compare the extent to which effects of P-element mutations on the starvation phenotype correspond to changes in transcript abundance in response to starvation. This comparison allows us to assess the hypothesis that changes in transcript abundance can be used to identify candidate genes with effects on phenotype, an hypothesis implicit in previous microarray studies [57]. Overall, there was no statistical association between the phenotypic and transcript data (χ21 = 0.0006, P = 1). For 194 genes, there was agreement between the phenotype and the expression level. Seventy-seven genes had significant differences in both transcript profile and mutant phenotypes, and 117 genes affected neither phenotype nor expression level (Additional data file 4). There was disagreement between the expression and phenotypic analyses for 191 genes (49.6%): 108 of the genes tagged by P-elements affected starvation resistance, but did not display differences in transcript level in response to starvation stress, and P-element insertions in 83 genes that exhibited significant differences in transcription in response to starvation did not have significant phenotypic effects on starvation tolerance (Additional data file 4).

The genetic architecture of transcription

A total of 706 probe sets exhibited variation in expression among the six lines; 640 probe sets were significant (P < 0.001) for the main effect of line, 190 for the line × sex interaction, 200 for the line × treatment interaction, and 85 for the three-way interaction of line × sex × treatment (Additional data file 5, and Figure 2d-k). Thus, transcript abundance exhibits both genotype by sex and genotype by environment interaction.

We used post-hoc Tukey tests to group lines with similar levels of gene expression, and compared the expression clusters with the Ore and 2b genotype of the six lines. There are three possible scenarios by which genetic variation in transcript abundance could arise. First, genetic variation in regulatory regions of gene A causes variation in the expression of gene A (cis-acting regulatory variation). Second, genetic variation in regulation of gene B causes variation in expression of A, which is itself not genetically variable (trans-acting regulatory variation). Third, genetic variation in both gene A and gene B affect the transcript abundance of gene A (cis- and trans-acting regulatory variation). These two-locus interactions could be additive or epistatic. We observe whether or not expression of gene A co-segregates with markers differentiating the two parental strains. Co-segregation will always be observed in case 1. It could also be observed in cases 2 and 3 if gene B is tightly linked to gene A, such that it is not separated by recombination from A in the genotypes tested. However, co-segregation will not be observed if gene A and gene B are unlinked. The most prevalent observation was regulation of expression by unlinked genes. For example, there were unambiguous interpretations for 246 probe sets that were significant for the main effect of line only: 65 (26.4%) were regulated by linked genes and 181 (73.5%) were regulated by unlinked genes (Additional data file 6, and Figure 2l-o). We also inferred linkage of genes regulating expression levels under control and starved conditions separately. There were unambiguous Tukey interpretations for 277 probe sets under control conditions, of which 32 exhibited linked regulatory variation (11.6%) and 245 were regulated by variation at unlinked genes (88.4%). For 244 probe sets under starved conditions, 46 were regulated by polymorphism at linked loci, (18.9%) and 198 were regulated by variation at unlinked genes (81.1%) (Additional data file 7).

Association of genetic variance in transcription with QTLs

Probe sets from the three-way ANOVA that are significant for the main effect of line and/or line × sex (P < 0.001), but not significant for the line × treatment interaction terms, exhibit genetic variation in transcription among the six lines that is independent of the starvation treatment. A total of 489 probe sets met these criteria, and we know the cytological locations of 475 of the corresponding genes. Previously, RI lines derived from Ore and 2b have been used to map QTL affecting variation in life span [2225], sensory bristle numbers [26], ovariole number [27], courtship signal [28], olfactory behavior [29], metabolism and flight [30], as well as starvation resistance [21]. Genes that exhibit significant differences for the main effect of line and/or line × sex which are located within QTL regions are putative candidate genes corresponding to the QTL [16, 20]. We identified several novel putative candidate genes affecting these traits (Additional data file 5). We examined whether probe sets with significant line and/or line × sex effects tended to cluster within regions containing QTL mapped under standard culture conditions, as would be the case if QTL regions were enriched for genes exhibiting transcriptional variation between the parental lines. We found no evidence for such clustering; indeed, the only trait showing a non-random association of probe sets with QTL that survived a Bonferroni correction for multiple tests was in the direction of a deficiency of probe sets in the QTL region (Table 3).

Table 3 Association of genetic variation in transcription with genetic variation in quantitative traits

The 217 probe sets with significant line × treatment and/or line × treatment × sex terms (Additional data file 5) represent genetic differences among the lines in response to the starvation treatment. Are these probe sets enriched in regions to which starvation resistance QTL map? We found that 47 of the probe sets meeting these criteria, representing 45 unique genes, fell within starvation resistance QTL regions; and the remaining 170 probe sets, representing 169 unique genes, fell outside the QTL intervals. These probe sets were not over-represented within starvation resistance QTL (χ21 = 0.26, P > 0.05).

There is significant variation in starvation half-life among the six lines (P < 0.0001; Additional data file 8). For those probe sets previously identified as having significant differences in transcript level among the lines, we assessed the extent to which variation in transcript abundance was associated with variation in starvation half-life. We found 281 probe sets with significant correlations (P < 0.05) between starvation phenotype and transcript level, for 273 of which the cytological location was known (Additional data file 5). However, 66 of the probe sets associated with starvation half-life mapped to starvation resistance QTL, and 207 did not. Again, these probe sets were not over-represented within starvation resistance QTL (χ21 = 0.45, P > 0.05).

Although there is no tendency for genes exhibiting variation in transcript abundance among lines to cluster within starvation resistance QTLs, those that do co-localize with the QTLs are candidate genes affecting variation in starvation tolerance between Ore and 2b. We found 155 probe sets, corresponding to 153 candidate genes, which met one or more of the above criteria (Additional data file 5). Most (114, 75%) were predicted genes. The remaining genes (Table 4) are reasonable candidates for starvation resistance QTLs, affecting the processes of protein metabolism, defense/immune response, proteolysis and peptidolysis, and transport.

Table 4 Candidate QTLs for starvation resistance

Complementation tests to mutations have implicated several candidate genes affecting variation between Ore and 2b in olfactory behavior [29] (Vanaso), longevity [31, 32] (Dopa decarboxylase, shuttle craft and ms(2)35Ci) and starvation resistance [21] (spalt major, Ryanodine receptor 44F, crooked legs, NaCP60E, Phosphoglucose isomerase, bellwether, numb, Punch, l(2)rG270, l(2)k17002, l(2)k00611, and l(2)k03205). None of these genes exhibited significant differences in transcript abundance between lines.


The sexually dimorphic transcriptome

Consistent with previous reports [5, 11, 33, 34], we observed highly significant differences in transcript abundance between males and females for nearly half the genome. These differences in transcriptional profiles were not confined to stereotypical sex-specific biological processes. Female transcript levels were upregulated for genes involved in protein biosynthesis, metabolism, and transcription regulation, while male transcript levels were higher for probe sets involved in ion and carrier transporters, as in a previous study of sex differences in transcription in Drosophila heads [5]. Differences in transcript abundance between the sexes may be an underlying mechanism for commonly observed sex-specific effects of QTLs associated with a variety of complex traits in Drosophila [2126, 32, 35, 36] and other organisms [37]. Males and females are effectively different environments in which genes act. The chromosomal locations of genes with sex-dependent expression were non-random. We confirmed the apparently general phenomenon that the Drosophila X chromosome is depauperate for genes that are upregulated in males [33, 34]; X-chromosome demasculinization is perhaps attributable to selection against genes that are advantageous in males but deleterious to females [33]. In contrast to previous studies, we observed that chromosome 2R harbored an excess, and chromosome 4 a deficiency, of genes that were upregulated in both males and females.

Transcriptional response to starvation stress

The transcriptional response to starvation stress involved approximately 25% of the genome. The stress profile indicates upregulation of genes involved in growth and maintenance processes and protein biosynthesis, with increased transcription of genes encoding translation initiation and elongation factors, mitochondrial and cytosolic ribosomal structural proteins, and hydrolases involving acid anhydrides. This increase in protein biosynthesis and hydrolase activity can be interpreted as an attempt to use available proteins for nourishment. A similar phenomenon has been observed in the response of yeast [38] and mammalian cells [39] to starvation, where substantial protein and organelle degradation provides substrate to starving cells [40]. Our observation that peptidases, which catalyze the hydrolysis of peptide bonds, were significantly downregulated in response to starvation, is consistent with the preservation of nascent protein chains. The downregulation of carrier activity and defense/immunity proteins indicates that transport across cell membranes slows and the immune response is compromised in starving flies.

We compared our results to those of a previous microarray study investigating gene-expression changes in starved larvae [41]. We found 21 probe sets that were significantly altered in both studies during starvation. Many of these genes have predicted functions that have not been verified experimentally; however, a few of the genes have known functions. Insulin-like Receptor, Serine pyruvate aminotransferase, Amylase distal, and mitochondrial carnitine palmitoyltransferase I, genes known to be involved in metabolism, were common to the two studies. Interestingly, Peroxidasin, a gene involved in oxygen and reactive oxygen species metabolism was upregulated fourfold in larvae, while it was downregulated 1.61-fold in our study.

Starvation stress was accompanied by reduced expression of genes affecting gametogenesis, by as much as 66-fold in starved female flies. Egg components such as chorion, yolk, and vitelline membrane proteins were among the most severely restricted transcripts, implicating suppression of female reproductive function during starvation. This depressed reproductive function is not unique to flies, as female mice on a calorically restricted diet experience a cessation in estrous cycle [42] and amenorrhea is one of the hallmarks of anorexia nervosa in human females [43]. Several male accessory gland proteins were also downregulated by as much as 6.5-fold during starvation stress. Oddly, six genes affecting spermatogenesis had significantly different levels of transcript abundance between the control and starved flies in both males and females; we found no male-specific differences in transcript abundance for genes involved in spermatogenesis (Additional data files 1 and 2).

Transcription of Rpd3 and CG14095 was upregulated in females and downregulated in males during starvation. Rpd3 is a transcriptional co-repressor, while the function of CG14095 is unknown. Sex-antagonistic patterns of expression have been observed in liver tissue studies of ethanol-fed rats [44], suggesting that these expression patterns may not be unique to flies.

The large number of transcripts altered during starvation implies massive pleiotropy; even more so when our conservative significance threshold is taken into account. This is consistent with our previous observation that 383 of 933 single P-element insertion lines tested (41%) had direct effects on starvation tolerance [21]. Further, candidate genes identified from the P-element screen and from complementation tests of QTL alleles to mutations at positional candidate genes are pleiotropic, and affect cell fate specification, cell proliferation, oogenesis, metabolism, and feeding behaviors [21].

Transcript abundance versus mutations

To what extent do candidate genes affecting response to starvation stress identified from changes in transcript abundance coincide with those implicated by assessing quantitative effects of P-element insertions on starvation tolerance? The resounding lack of an overall statistical association between the two methods is somewhat deceptive. While there was no association overall, if we had only tested the 160 P-element mutations corresponding to genes with altered transcript abundance during starvation, we would have found that 77 (48%) actually had phenotypic effects on starvation resistance. The lack of association was caused by 108 genes tagged by P-elements that affected starvation resistance, but did not display differences in transcript level in response to starvation stress, and P-element insertions in 83 genes that exhibited significant differences in transcription in response to starvation but did not have significant phenotypic effects on starvation tolerance. Genes affecting starvation that are regulated post-transcriptionally, or for which differences in transcript abundance that are undetectable on the array have large phenotypic consequences, contribute to the first source of discordance between the two methods. The second source of discordance could arise if the genes exhibiting expression changes during starvation are truly candidate genes affecting starvation resistance, but the particular P-element insertional mutation tested was not in a region affecting the starvation phenotype; a P-element insertion or point mutation in another location might produce a significant effect on starvation tolerance [45]. Another possibility is that the gene is downregulated during starvation; thus, a P-element mutation in the gene might not have an effect on the starvation resistance phenotype. Alternatively, a fraction of these probe sets could be false positives. Therefore, we conclude that assessing the effects of mutations at genes exhibiting changes in transcript abundance in response to an environmental (or genetic [5, 7]) perturbation is a highly efficient strategy for identifying networks of pleiotropic genes regulating complex traits.

Genetic variation in transcript abundance and quantitative trait phenotypes

The prospects for easily identifying genes corresponding to QTLs using microarray profiling seem less rosy at present. It has been proposed that candidate genes corresponding to QTLs are those for which expression differs between the parental strains used to construct the QTL mapping population, and which are located in the regions to which the QTLs map [20]. However, differences in expression between lines could be due to polymorphisms between the tested strains and the strain used to construct the probe sets on the array. Further, the lines differ for many traits, and QTLs affecting them overlap; unless the QTLs are mapped with very high resolution, candidate genes chosen by this criterion alone could affect another trait. The issue of polymorphism can be circumvented for traits with environmentally conditional expression by considering probe sets exhibiting a line × treatment environment interaction, and trait specificity can be addressed by correlating expression levels with the trait phenotype. None of these criteria led to an enrichment of candidate genes with variation in expression within QTL regions.

A major difficulty in using changes in gene expression between two strains to identify candidate genes corresponding to QTLs arises because variation in transcript abundance for positional candidate genes could arise from several causes. First, variation in transcript abundance is attributable to regulatory polymorphism in the candidate gene itself. Second, the candidate gene is itself not genetically variable, but regulatory variation in a second gene affects variation in its expression. Third, variation in transcript abundance at the candidate gene is attributable to interacting regulatory polymorphisms in both the candidate gene and a second gene. These interactions could be additive or epistatic. Positional candidate genes with variation in transcript abundance arising from the first or third cause could potentially correspond to genetically variable QTLs. However, it is becoming clear that genetic variation in transcript abundance is largely attributable to regulation by unlinked genes (see [1517, 19] and this paper). Indeed, single P-element insertions can alter the transcript expression of as many as 161 genes compared to a co-isogenic control line [5]. This low signal-to-noise ratio means that choosing positional candidate genes for further study based only on differences in transcript abundance between parental lines does not have a high likelihood of success.

In the future, the falling cost of whole-genome expression analysis will facilitate assessing transcriptional variation and variation in trait phenotypes in the same large QTL mapping populations. Co-localization of QTLs with main effects jointly affecting variation in transcription and trait phenotype will help winnow out monomorphic genes that are regulated by unlinked loci, and such data would enable direct tests for epistasis at the level of transcription and the trait. It is unlikely that this approach will completely supplant high-resolution QTL mapping and complementation tests to mutations for elucidating the genetic architecture of complex traits in Drosophila. None of the 12 candidate genes affecting variation in starvation resistance between Ore and 2b [21] exhibited variation in transcript abundance in this study. Possibly any transcriptional differences between Ore and 2b alleles at these loci are rare messages below the threshold of detection, or that are expressed in only a few cell types or at a particular period of development. In addition, not all allelic differences between QTL alleles are necessarily regulated at the level of transcription. Nevertheless, incorporation of knowledge about variation in transcript abundance will greatly inform our choice of candidate genes for confirmation by mutant complementation tests and association studies, which is currently biased by our poor understanding of the pleiotropic and epistatic consequences of variation in positional candidate genes on variation in trait phenotypes.

Materials and methods


We used the isogenic lines 2b [22, 46] and Oregon-R [47] (Ore) to establish 98 RI lines for mapping QTLs affecting starvation resistance [21]. Survival times for Oregon-R flies were 36.0 and 51.6 h for males and females, respectively. For 2b, survival times were 29.2 h for males and 40.4 h for females. Here, we assessed transcriptional profiles under control conditions and during starvation for 2b, Ore, two starvation resistant (RI.14, RI.21) and two starvation sensitive (RI.35, RI.42) RI lines. Recombination breakpoints for the RI lines have been determined previously [23] and are resolved to the nearest cytological lettered subdivision. We maintained control flies on cornmeal-agar-molasses medium, and starved flies on non-nutritive (1.5% agar and water) medium, under standard culture conditions (25°C, 70% humidity, and a 12-h light: 12-h dark cycle).

Starvation half-life

We assessed survival of all six lines under starvation conditions by placing two replicates of ten flies each per sex on starvation medium, and recording the number of flies alive at 8-h intervals until all were dead. We used these survival curves to infer the starvation half-life for each line/sex combination. We used an analysis of variance (ANOVA) model Y = μ + L + S + L × S + R(L × S) + E, to partition variance in survival times into sources attributable to the cross-classified main effects of lines (L), sex (S), variance between replicate vials (R), and within-vial environmental variance (E).

Transcriptional profiling

For each of two independent replicates, we collected 300 male and 300 female virgins from all lines, aged 2-5 days post-eclosion. The control treatment consisted of 100 non-starved flies/line/sex. We placed the remaining 200 flies/line/sex on starvation medium, and collected approximately 100 flies/line/sex at the predetermined starvation half-life. Starved flies from all lines should therefore be in roughly the same physiological condition. We extracted whole-body RNA from each of the 48 independent samples (6 lines × 2 treatments × 2 sexes × 2 replicates) with Triazol reagent (Gibco BRL), followed by DNase digestion (RQ1 DNase, Promega,) and a 1:1 phenol (Sigma-Aldrich)-chloroform (Fisher Scientific) extraction. We hybridized biotinylated cRNA probes to single-color whole-genome Affymetrix Drosophila GeneChip arrays as described in the Affymetrix GeneChip Expression Analysis 2000 manual.

Data analysis

We normalized the expression data by scaling overall probe set intensity to 100 on each chip using standard reference probe sets on each chip for the normalization procedure. Each probe set on the array consists of 14 perfect match (PM) and single nucleotide mismatch (MM) pairs. We used the average difference (AD) in normalized RNA expression between the 14 perfect match (PM) and mismatch (MM) probe pairs per probe set (Affymetrix Microarray Suite, Version 4.0) as the analysis variable. We calculated the minimum AD threshold value [5] as AD = 30. If the mean AD of a probe set was less than 30, and the maximum AD value was also less than 30, we eliminated the probe set from further consideration. We set all remaining AD scores < 30, to AD = 30. We performed a three-way factorial ANOVA of AD for each probe set, according to the model: Y = μ + S + T + L + S × T + S × L + T ×L + S × T × L + E, where S, T, and L represent, respectively, the fixed cross-classified effects of sex, treatment (control versus starved), and line, and E is the replicate variance between arrays. We determined F-ratio tests of significance for each term in the ANOVA, and considered probe sets with P values ≤ 0.001 for any term to be significant. (There are approximately 14,000 probe sets on the array; thus 14 false positives would be expected at this significance threshold.)

We computed the female:male ratio of AD values, averaged over all lines and treatments, for probe sets for which the main effect of S was significant. Similarly, we computed the starved:control ratio of AD values, averaged over lines and sex, for probe sets with significant T terms. We categorized these probe sets according to their gene ontology (GO) for biological process and molecular function [48]. We assessed significant differences in GO categories between up- and downregulated probe sets using G tests [49], under the null hypothesis of equal numbers of up- and down- regulated probe sets in each category, and using Bonferroni corrections to account for multiple tests. For probe sets with significant T × S terms, we ran two-way ANOVAs separately by sex using the reduced model Y = μ + L + T + L × T + E.

Probe sets with significant L, L × S, L × T or L × T × S terms are candidate QTLs for traits that vary among the lines. We performed post-hoc Tukey tests for all probe sets for which these terms were significant to determine in which lines transcription was up-or downregulated in response to starvation stress. For probe sets that were significant for the main effect of L, but not any of the interaction terms, we conducted Tukey tests using the expression values pooled across control and starved conditions and both sexes. We computed Tukey tests separately for males and females, averaged over both treatments, for probe sets that were significant for the L × S interaction; and separately by treatment, for probe sets significant for the L × T interaction. The Tukey analyses separated the lines into groups within which AD values were not significantly different. Since the genotype for each recombinant inbred line at any given location is known, we used the Tukey analyses to classify probe sets as exhibiting linked or unlinked regulation of transcript abundance. We considered linked factors to regulate transcript abundance if Ore and 2b differ in transcript abundance, and this difference is reflected in the RI lines according to their Ore and 2b genotype in the region to which the gene maps. Conversely, we inferred that unlinked factors regulate transcript abundance in cases where there is not a 1:1 correspondence between parental line genotype and Tukey grouping. We determined the fold-change between Tukey groupings by calculating the ratio of the deviant line(s) expression level to the mean expression level of the parental or common group. Most Tukey analyses were unambiguous; where multiple interpretations were possible, we calculated the fold-change for all possibilities.

Statistical analyses

We used SAS procedures for all statistical analyses [50].

Additional data files

The following additional data files are available with the online version of this article. Additional data file 1 contains a list of all probe sets with significantly different expression in females and males. Additional data file 2 contains a list of all probe sets with significantly different expression under control and starved conditions. Additional data file 3 lists the probe sets for which the sex by treatment interaction term is significant. Additional data file 4 shows the correspondence between the results of a screen for the effects on resistance to starvation stress for single P-element inserts, in a co-isogenic background [21], and changes in transcript abundance between control and starved treatments. Additional data file 5 summarizes probe sets for which there is significant genetic variation in transcript abundance. Additional data file 6 shows the probe sets for which the only significant genetic term was the main effect of line. Additional data file 7 gives the same information as Additional data file 6, but separately for the control and starved treatments, and with the results of the analyses pooled over sexes, and for males and females separately. Additional data file 8 is the ANOVA of starvation half-life for Ore, 2b, RI.14, RI.21, RI.35 and RI.42. Additional data files 9 and 10 give the raw expression data and presence/absence calls for the control and starved treatments, respectively.


  1. 1.

    Mackay TFC: The genetic architecture of quantitative traits. Annu Rev Genet. 2001, 35: 303-339. 10.1146/annurev.genet.35.102401.090633.

  2. 2.

    Lee SS, Lee RYN, Fraser AG, Kamath RS, Ahringer J, Ruvkun G: A systematic RNAi screen identifies a critical role for mitochondria in C. elegans longevity. Nat Genet. 2003, 33: 40-48. 10.1038/ng1056.

  3. 3.

    Norga KK, Gurganus MC, Dilda CL, Yamamoto A, Lyman RF, Patel PH, Rubin GM, Hoskins RA, Mackay TFC, Bellen HJ: Quantitative analysis of bristle number in Drosophila mutants identifies genes involved in neural development. Curr Biol. 2003, 13: 1388-1397. 10.1016/S0960-9822(03)00546-3.

  4. 4.

    Glazier AM, Nadeau JH, Aitman TJ: Finding genes that underlie complex traits. Science. 2002, 298: 2345-2349. 10.1126/science.1076641.

  5. 5.

    Anholt RRH, Dilda CL, Chang S, Fanara JJ, Kulkarni NH, Ganguly I, Rollmann SM, Kamdar KP, Mackay TFC: The genetic architecture of odor-guided behavior in Drosophila: epistasis and the transcriptome. Nat Genet. 2003, 35: 180-184. 10.1038/ng1240.

  6. 6.

    Murphy CT, McCarroll SA, Bargmann CI, Fraser A, Kamath RS, Ahringer J, Li H, Kenyon C: Genes that act downstream of DAF-16 to influence the lifespan of Caenorhabditis elegans. Nature. 2003, 424: 277-283. 10.1038/nature01789.

  7. 7.

    Toma DP, White KP, Hirsch J, Greenspan RJ: Identification of genes involved in Drosophila melanogaster geotaxis, a complex behavioral trait. Nat Genet. 2002, 31: 349-353.

  8. 8.

    Lee CK, Klopp RG, Weindruch R, Prolla TA: Gene expression profile of aging and its retardation by caloric restriction. Science. 1999, 285: 1390-1393. 10.1126/science.285.5432.1390.

  9. 9.

    Lee CK, Weindruch R, Prolla TA: Gene-expression profile of the aging brain in mice. Nat Genet. 2000, 25: 294-297. 10.1038/77046.

  10. 10.

    Zou S, Meadows S, Sharp L, Jan LY, Jan JN: Genome-wide study of aging and oxidative stress response in Drosophila melanogaster. Proc Natl Acad Sci USA. 2000, 97: 13726-13731. 10.1073/pnas.260496697.

  11. 11.

    Jin W, Riley RM, Wolfinger RD, White KP, Passador-Gurgel G, Gibson G: The contributions of sex, genotype, and age to transcriptional variance in Drosophila melanogaster. Nat Genet. 2001, 29: 389-395. 10.1038/ng766.

  12. 12.

    Pletcher SD, Macdonald SJ, Marguerie R, Certa U, Stearns SC, Goldstein DB, Partridge L: Genome-wide transcript profiles in aging and calorically restricted Drosophila melanogaster. Curr Biol. 2002, 12: 712-723. 10.1016/S0960-9822(02)00808-4.

  13. 13.

    Cowles CR, Hirschhorn JN, Altshuler D, Lander ES: Detection of regulatory variation in mouse genes. Nat Genet. 2002, 32: 432-437. 10.1038/ng992.

  14. 14.

    Oleksiak MF, Churchill GA, Crawford DL: Variation in gene expression within and among natural populations. Nat Genet. 2002, 32: 261-266. 10.1038/ng983.

  15. 15.

    Brem RB, Yvert G, Clinton R, Kruglyak L: Genetic dissection of transcriptional regulation in budding yeast. Science. 2002, 296: 752-755. 10.1126/science.1069516.

  16. 16.

    Schadt EE, Monks SA, Drake TA, Lusis AJ, Che N, Colinayo V, Ruff TG, Milligan SB, Lamb JR, Cavet G, et al: Genetics of gene expression surveyed in maize, mouse and man. Nature. 2003, 422: 297-302. 10.1038/nature01434.

  17. 17.

    Yvert G, Brem RB, Whittle J, Akey JM, Foss E, Smith EN, Mackelprang R, Kruglyak L: Trans-acting regulatory variation in Saccharomyces cerevisiae and the role of transcription factors. Nat Genet. 2003, 35: 57-64. 10.1038/ng1222.

  18. 18.

    Jansen RC, Nap JP: Genetical genomics: the added value from segregation. Trends Genet. 2001, 17: 388-391. 10.1016/S0168-9525(01)02310-1.

  19. 19.

    Morley M, Molony CM, Weber TM, Devlin JL, Ewens KG, Spielman RS, Cheung VG: Genetic analysis of genome-wide variation in human gene expression. Nature. 2004, 430: 743-747. 10.1038/nature02797.

  20. 20.

    Wayne ML, McIntyre LM: Combining mapping and arraying: an approach to candidate gene identification. Proc Natl Acad Sci USA. 2002, 99: 14903-14906. 10.1073/pnas.222549199.

  21. 21.

    Harbison ST, Yamamoto AH, Fanara JJ, Norga KK, Mackay TFC: Quantitative trait loci affecting starvation resistance in Drosophila melanogaster. Genetics. 2004, 166: 1807-1823. 10.1534/genetics.166.4.1807.

  22. 22.

    Vieira C, Pasyukova EG, Zeng ZB, Hackett JB, Lyman RF, Mackay TFC: Genotype-environment interaction for quantitative trait loci affecting life span in Drosophila melanogaster. Genetics. 2000, 154: 213-227.

  23. 23.

    Nuzhdin SV, Pasyukova EG, Dilda CL, Zeng ZB, Mackay TFC: Sex-specific quantitative trait loci affecting longevity in Drosophila melanogaster. Proc Natl Acad Sci USA. 1997, 94: 9734-9739. 10.1073/pnas.94.18.9734.

  24. 24.

    Leips J, Mackay TFC: Quantitative trait loci for life span in Drosophila melanogaster : Interactions with genetic background and larval density. Genetics. 2000, 155: 1773-1788.

  25. 25.

    Leips J, Mackay TFC: The complex genetic architecture of Drosophila life span. Exp Aging Res. 2002, 28: 361-390. 10.1080/03610730290080399.

  26. 26.

    Gurganus MC, Fry JD, Nuzhdin SV, Pasyukova EG, Lyman RF, Mackay TFC: Genotype-environment interaction at quantitative trait loci affecting sensory bristle number in Drosophila melanogaster. Genetics. 1998, 149: 1883-1898.

  27. 27.

    Wayne ML, Hackett JB, Dilda CL, Nuzhdin SV, Pasyukova EG, Mackay TFC: Quantitative trait loci mapping of fitness-related traits in Drosophila melanogaster. Genet Res. 2001, 77: 107-116. 10.1017/S0016672300004894.

  28. 28.

    Gleason JM, Nuzhdin SV, Ritchie MG: Quantitative trait loci affecting a courtship signal in Drosophila melanogaster. Heredity. 2002, 89: 1-6. 10.1038/sj.hdy.6800099.

  29. 29.

    Fanara JJ, Robinson KO, Rollmann SM, Anholt RRH, Mackay TFC: Vanaso is a candidate quantitative trait gene for Drosophila olfactory behavior. Genetics. 2002, 162: 1321-1328.

  30. 30.

    Montooth KL, Marden JH, Clark AG: Mapping determinants of variation in energy metabolism, respiration and flight in Drosophila. Genetics. 2003, 165: 623-635.

  31. 31.

    De Luca M, Roshina NV, Geiger-Thornsberry GL, Lyman RF, Pasyukova EG, Mackay TFC: Dopa decarboxylase (Ddc) affects variation in Drosophila longevity. Nat Genet. 2003, 34: 429-433. 10.1038/ng1218.

  32. 32.

    Pasyukova EG, Roshina NV, Mackay TFC: shuttle craft : A candidate quantitative trait gene for Drosophila life span. The Ageing Cell. 2004, 3: 297-307. 10.1111/j.1474-9728.2004.00114.x.

  33. 33.

    Parisi M, Nuttall R, Naiman D, Bouffard G, Malley J, Andrews J, Eastman S, Oliver B: Paucity of genes on the Drosophila X chromosome showing male-biased expression. Science. 2003, 299: 697-700. 10.1126/science.1079190.

  34. 34.

    Rantz JM, Castillo-Davis CI, Meiklejohn CD, Hartl DL: Sex-dependent gene expression and evolution of the Drosophila transcriptome. Science. 2003, 300: 1742-1745. 10.1126/science.1085881.

  35. 35.

    Reiwitch SG, Nuzhdin SV: Quantitative trait loci for lifespan of mated Drosophila melanogaster affect both sexes. Genet Res. 2002, 80: 225-230. 10.1017/S0016672302005943.

  36. 36.

    Dilda CL, Mackay TFC: The genetic architecture of Drosophila sensory bristle number. Genetics. 2002, 162: 1655-1674.

  37. 37.

    Butterfield RJ, Roper RJ, Rhein DM, Melvold RW, Haynes L, Ma RZ, Doerge RW, Teuscher C: Sex-specific quantitative trait loci govern susceptibility to Theiler's murine encephalomyelitis virus-induced demyelination. Genetics. 2003, 163: 1041-1046.

  38. 38.

    Teichert U, Mechler B, Muller H, Wolf DH: Lysosomal (vacuolar) proteinases of yeast are essential catalysts for protein degradation, differentiation, and cell survival. J Biol Chem. 1989, 264: 16037-16045.

  39. 39.

    Seglen PO, Bohley P: Autophagy and other vacuolar protein degradation mechanims. Experientia. 1992, 48: 158-172.

  40. 40.

    Kim J, Klionsky DJ: Autophagy, cytoplasm-to-vacuole targeting pathway, and pexophagy in yeast and mammalian cells. Annu Rev Biochem. 2000, 69: 303-342. 10.1146/annurev.biochem.69.1.303.

  41. 41.

    Zinke I, Schutz CS, Katzenberger JD, Bauer M, Pankratz MJ: Nutrient control of gene expression in Drosophila: microarray analysis of starvation and sugar-dependent response. EMBO J. 2002, 21: 6162-6173. 10.1093/emboj/cdf600.

  42. 42.

    Nelson JF, Karelus K, Bergman MD, Felicio LS: Neuroendocrine involvement in aging: Evidence from studies of reproductive aging and caloric restriction. Neurobiol Aging. 1995, 16: 837-843. 10.1016/0197-4580(95)00072-M.

  43. 43.

    Reijonen JH, Pratt HD, Patel DR, Greydanus DE: Eating disorders in the adolescent population: an overview. J Adolescent Res. 2003, 18: 209-222. 10.1177/0743558403018003002.

  44. 44.

    Tadic SD, Elm MS, Li HS, Van Londen GJ, Subbotin VM, Whitcomb DC, Eagon PK: Sex differences in hepatic gene expression in a rat model of ethanol induced liver injury. J Appl Physiol. 2002, 93: 1057-1068.

  45. 45.

    Sokolowski MB: Drosophila: genetics meets behavior. Nat Rev Genet. 2001, 2: 879-890. 10.1038/35098592.

  46. 46.

    Pasyukova EG, Nuzhdin SV: Doc and copia instability in an isogenic Drosophila melanogaster stock. Mol Gen Genet. 1993, 240: 302-306. 10.1007/BF00277071.

  47. 47.

    Lindsley DL, Zimm GG: The Genome of Drosophila melanogaster. 1992, San Diego, CA: Academic Press

  48. 48.

    The FlyBase Consortium: The FlyBase database of the Drosophila genome projects and community literature. Nucleic Acids Res. 2003, 31: 172-175. 10.1093/nar/gkg094. []

  49. 49.

    Sokal RR, Rohlf FJ: Biometry. 1995, New York: WH Freeman, 697-698. Third

  50. 50.

    SAS Institute: SAS/STAT User's Guide, Release 6.12 Edition. 1988, Cary, NC: SAS Institute

Download references


We thank K. Norga for annotating the P-element insertion lines. S. T. H. was the recipient of a W. M. Keck pre-doctoral fellowship. This work was funded by grants from the National Institutes of Health to T. F. C. M. This is a publication of the W. M. Keck Center for Behavioral Biology.

Author information

Correspondence to Trudy FC Mackay.

Electronic supplementary material

Additional File 1: The file includes all P-values from ANOVA of expression values; the mean expression in males and females, averaged across treatments and lines; the FlyBase ID, gene name, symbol and synonyms; cytological location; and molecular function, biological process and cellular component gene ontologies (XLS 3 MB)

Additional File 2: The file includes the mean expression levels under control and starved conditions, averaged over sexes and lines, as well as the information given for each probe set in Additional data file 1 (XLS 2 MB)

Additional File 3: The file includes the P-values for the treatment term in the reduced analyses for males and females separately; the mean expression values under control and starved conditions for males and females, averaged over lines; whether expression is sex-biased (SB), sex-specific (SS) or sex-antagonistic (SA); plus the information given for each probe set in additional data file 1 (XLS 474 KB)

Additional File 4: The results are grouped into four categories: (1) Transcript abundance is altered between control and starved treatments, and there is a significant effect of the P-element insertion on starvation tolerance; (2) Transcript abundance is altered between control and starved treatments, but the P-element insertion does not significantly affect starvation tolerance; (3) There is no alteration in transcript abundance between control and starved treatments, but the P-element insertion significantly affects starvation tolerance; and (4) There is no alteration in transcript abundance between control and starved treatments, and no significant effect of the P-element insertion on starvation tolerance (XLS 100 KB)

Additional File 5: The file includes all P-values from ANOVA of expression values; the gene name, symbol, synonyms, and cytological location; the molecular function, biological process and cellular location gene ontologies; co-localization of the gene with QTLs affecting life span, sensory bristle numbers, starvation resistance, ovariole number, courtship signal, flight, metabolic rate and triglycerides; and the statistical association of variation among lines in transcript levels and starvation half-life (XLS 464 KB)

Additional File 6: The cells are color-coded: purple for Ore genotype, green for 2b genotype, gray if the genotype is unknown because the gene is between an Ore and a 2b flanking marker, and the exact recombination breakpoint is not determined, and gold if the genotype is unknown because of residual heterozygosity in the RI line. The results of Tukey tests separating the lines into groups within which expression values are not significantly different are given. If more than one interpretation of the Tukey groups are possible, all are given. Linked (L) regulation of variation in transcript abundance is inferred if the Tukey groupings match the genotype. Regulation of variation in transcript abundance is inferred to be linked (IL) if the unknown genotypes could match this interpretation. Unlinked (U) regulation of variation in transcript abundance is inferred if the Tukey groupings do not match the line genotypes (XLS 322 KB)

Additional File 7: A table giving the same information as in Additional data file 6, but separately for the control and starved treatments, and with the results of the analyses pooled over sexes, and for males and females separately (XLS 282 KB)

Additional File 8: A table showing the ANOVA of starvation half-life for Ore, 2b, RI.14, RI.21, RI.35 and RI.42 (DOC 26 KB)

Additional File 9: A table showing the raw expression data and presence/absence calls for the control treatments (XLS 9 MB)

Additional File 10: A table showing the raw expression data and presence/absence calls for the starved treatments (XLS 9 MB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Rights and permissions

Reprints and Permissions

About this article


  • Quantitative Trait Locus
  • Transcript Abundance
  • Additional Data File
  • Quantitative Trait Locus Region
  • Recombinant Inbred