- Research article
- Open Access
Invariance and plasticity in the Drosophila melanogastermetabolomic network in response to temperature
© Hariharan et al.; licensee BioMed Central. 2015
- Received: 19 August 2014
- Accepted: 11 December 2014
- Published: 24 December 2014
Metabolomic responses to extreme thermal stress have recently been investigated in Drosophila melanogaster. However, a network level understanding of metabolomic responses to longer and less drastic temperature changes, which more closely reflect variation in natural ambient temperatures experienced during development and adulthood, is currently lacking. Here we use high-resolution, non-targeted metabolomics to dissect metabolomic changes in D. melanogaster elicited by moderately cool (18°C) or warm (27°C) developmental and adult temperature exposures.
We find that temperature at which larvae are reared has a dramatic effect on metabolomic network structure measured in adults. Using network analysis, we are able to identify modules that are highly differentially expressed in response to changing developmental temperature, as well as modules whose correlation structure is strongly preserved across temperature.
Our results suggest that the effect of temperature on the metabolome provides an easily studied and powerful model for understanding the forces that influence invariance and plasticity in biological networks.
- Drosophila melanogaster
- Differential coexpression
Temperature has profound effects on cellular and organismal biochemistry , with enzymes functioning best at specific temperatures . The effects of temperature on physiology and fitness have been particularly well-studied in the fruit fly, Drosophila melanogaster, which can survive transient exposures to a wide range of temperatures, from −10°C to 40°C . Early studies found a strong negative correlation between survival rate and chronic exposure to high temperatures, while showing that transient exposure to extreme temperatures could increase survival (so-called ‘heat hardening’) -. Molecular studies have shown that this effect is due, in part, to both heat- and cold-shock inducing the expression of heat shock proteins (HSPs), which help cells to counteract the deleterious effect of thermal shock through multiple mechanisms -.
Over the past several years, researchers have begun to dissect the underlying molecular mechanisms of temperature responses beyond HSPs, with numerous studies in D. melanogaster looking at the effect of temperature on the transcriptome, the proteome - and the metabolome ,,.
The value of studying metabolomic changes in response to temperature shifts in D. melanogaster has been demonstrated by numerous studies. For example, it has recently been shown that long-term cold acclimation and heat hardening alters the metabolomic profile of D. melanogaster larvae ,. Other studies have sought to identify temperature-associated metabolomic signatures of inbreeding in D. melanogaster ,. These metabolomic studies have led to a better understanding of the global effects of the temperature stress response in D. melanogaster, moving beyond earlier studies that focused on just a very small number of metabolites in the D. melanogaster temperature stress-response .
Here we extend these analyses in three ways. First, we take advantage of recent advances in high-resolution mass spectrometry -. This approach allows us to greatly increase the number of metabolites and metabolic pathways that we can assay for effects of temperature.
Second, we incorporate recent developments in differential network analysis ,. Numerous gene expression studies have shown that in some cases, perturbations can have limited effects on the magnitude of transcripts while causing substantial changes in the correlations between transcripts (e.g., ). More recently, researchers have brought this approach to the study of metabolomic networks -. In light of these results, here we explore not only the structure of the metabolomic network , but more specifically, how temperature changes the structure of specific modules within the larger metabolomic network. In particular, we are able to identify modules whose structure is highly constant in response to temperature, as well as modules that change their structure dramatically in response to temperature. Our focus on modules that maintain correlation structure, or alter that structure, in response to temperature, can pave the way to a better understanding of the mechanisms underlying robustness and plasticity of networks in response to environmental changes . To our knowledge, this study is the first to identify temperature-dependent constant and plastic modules in the metabolome, and suggests novel approaches to better understand how poikilothermic organisms might adapt to a changing environment. Moreover, the approaches we use here could lead to important insights into network evolution and the role that network structure plays in the ability of organisms to cope with stressors in general .
Finally, by studying the effect of developmental and adult temperature on high-resolution, untargeted metabolomic profiles, we identify novel associations between rearing conditions and metabolites and metabolic pathways.
Our dataset consisted of data from technical duplicates for 95 samples across four Drosophila Genome Reference Panel (DGRP) genotypes . These included 48 male and 47 female samples, with 24 samples per genotype except for DGRP 25189 which contributed only 23 samples to our analysis. After applying quality control procedures (see Methods), our dataset included 4359 features from a C18 column, and 2961 features from an AE column. We were able to assign putative matches to 1141 metabolites from the C18 column and 926 metabolites from the AE column. We thus had data for a total of 7,320 features with 2027 putatively identified features from both columns, noting that some features overlap between the two columns.
Pathway enrichment analysis using mummichog
For this analysis, we investigated the effects of developmental temperature and adult temperature on male and female flies separately.
Developmental temperature effects
Metabolic pathways altered by developmental temperature in males
Overlap size (b)
Pathway size (c)
Overlap features (id) (d)
Glycogen degradation I/ Trehalose biosynthesis
Salvage pathways of adenine, hypoxanthine, and their nucleosides
AMP, Deoxyadenosine, Xanthine, Adenine
Acyl carrier protein metabolism
Selenocysteine biosynthesis II (archaea and eukaryotes)
Metabolic pathways altered by developmental temperature in females
Overlap size (b)
Pathway size (c)
Overlap features (id) (d)
Glycogen degradation I
GLC-1-P, Alpha-Glucose, Maltotetraose, GLC, Maltotriose, GLC-6-P, Maltose, Alpha-GLC-6-P
Lactose/Melibiose degradation III
Lactose, Melibiose, Galactose, GLC
Trehalose biosynthesis I
Trehalose-6P, Trehalose, Alpha-GLC-6-P
tRNA charging pathway
Val, Pro, Thr, Phe, Arg, Trp
Salvage pathways of adenine, hypoxanthine, and their nucleosides
Xanthine, Adenine, Inosine, Deoxyadenosine, Deoxyinosine, AMP, Hypoxanthine
Sphingosine and sphingosine-1-phosphate metabolism
CPD3DJ-11366, NADP, Sphingosine, Palmitaldehyde
Adult temperature effects
Metabolic pathways altered by adult temperature in males and females
Overlap size (b)
Pathway size (c)
Overlap features (id) (d)
*Degradation of purine ribonucleosides
Adenosine, Guanosine, Adenine
4-hydroxyproline degradation I
GLT, L-4-hydroxy-proline, NAD
Arginine degradation I (arginase pathway)
L- glutamate_gamma-semialdehyde, GLT, NAD
NAD biosynthesis from 2-amino-3-carboxymuconate semialdehyde
Deamido-NAD, GLT, NAD
Warm and cool temperature effects
Metabolic pathways altered by developmental and adult temperature in males
Overlap size (b)
Pathway size (c)
Overlap features (id) (d)
Serotonin and melatonin biosynthesis
Trp, N-acetyl-serotonin, S-adenosylmethionine, N-acetyl-5-Methoxy-tryptamine
Metabolic pathways altered by developmental and adult temperature in females
Overlap size (b)
Pathway size (c)
Overlap features (id) (d)
Arginine degradation VI (arginase 2 pathway)
L-glutamate_gamma-semialdehyde, Pro, Arg
Salvage pathways of guanine, xanthine, and their nucleosides
Deoxyguanosine, GMP, Guanine
We were also interested in identifying metabolites whose intensities were affected significantly by genotype, i.e., metabolites whose levels differ among genotypes. In males, we identified 354 (C18) and 297 (AE) metabolites, while in females we detected 355 (C18) and 364 (AE) metabolites that exhibited significant differential abundance. Significantly enriched metabolic pathways in these sets of metabolites were identified using mummichog (Additional file 1). Seven and five metabolomic pathways were detected as being significantly affected by genotype in female and male flies respectively.
By contrast, we found very few metabolites whose temperature response was genotype-dependent. Specifically, we detected less than six metabolites each from the interaction analysis between different factors (genotype, developmental temperature, and adult temperature).
For illustrative purposes, we visualized representative DiffCoEx and Similarity matrix modules in Cytoscape , choosing the smallest of the detected modules for each sex. Comparison of the magenta module (Figure 3a and Additional file 4) from the male fly metabolome reveals that this module becomes much more tightly correlated when developmental temperature changes from 18°C to 27°C. In contrast, a visual comparison of the red module (Figure 3b and Additional file 5) for female samples suggests that this module loses a large proportion of its correlations when developmental temperature becomes warmer. Cytoscape based network visualizations for the most striking of the Similarity matrix modules for the male and female metabolomes are shown in Figure 4.
For each DiffCoEx and Similarity matrix modules, we used mummichog to test for metabolic pathway enrichment. We found significant enrichment in six out of 12 DiffCoEx modules in the male fly data, and in three out of 11 DiffCoEx modules in the female fly data (Additional file 6 and Additional file 7). Several novel pathways not detected by differential expression analysis were detected as significant in the network analysis, including glycolysis, lipoate biosynthesis and N-acetyl-glucosamine degradation. We did not detect any enriched metabolic pathways among the preserved models.
In this study, we have carried out two levels of analysis. The first level of analysis allowed us to identify main effects of temperature on metabolite levels. The second level of analysis, focused on the metabolomic network, enabled us to ask questions about network invariance and plasticity.
We investigated metabolomic thermal responses to moderate temperatures (27°C and 18°C). In contrast to earlier studies that have focused on the effects of stressful or extreme temperatures on the fly metabolome, here we focus on longer term exposure to temperatures that are considered to be relatively benign for this organism . Thermal responses have been extensively investigated in D. melanogaster, mostly looking at the effects of stressful temperatures during both development and adulthood on fly physiology ,-. The metabolomic effects of heat-shock, rapid cold hardening, long-term cold acclimation, and of benign and stressful temperatures on inbreeding, have all been previously studied in D. melanogaster ,,. In addition to these metabolomic studies, molecular mechanisms of many thermal responses in the fly have been investigated at the level of the transcriptome, and more recently, the proteome -,. One benefit of focusing on less extreme temperatures is that we can gain insight into ways that the metabolome responds to natural variation in temperature specifically, rather than response to stress.
Main effects of temperature on the metabolome
Several metabolic pathways are known to be affected by temperature in D. melanogaster, and our analysis looking at main effects of temperature recapitulated many of these earlier findings while extending them in new directions.
From previous studies we know that brief exposures to stressful temperatures, including cold stress, can cause accumulations of polysaccharides in D. melanogaster ,. Mechanisms by which maltose and trehalose protect cells during stressful conditions have been described previously . Stressful temperatures can also cause accumulations of proline in flies .
Cold acclimation in D. melanogaster is also known to be associated with significant metabolic changes involving sugars (sucrose, fructose, and trehalose), polyamines and a few metabolic intermediates . While that study used targeted metabolomics, results from our untargeted metabolomics approach also point to altered trehalose biosynthesis in both male and female flies reared at a lower developmental temperature. Our results suggests that these pathways play a role not only in the response to acute, potentially lethal stressors, but also to the range of temperatures a fly is likely to encounter from day to day.
A recent study implicated dopamine in modulating behavioral response to temperature changes in D. melanogaster . We also observed changes in dopamine metabolism with thermal responses in D. melanogaster. The agreement between our study and the earlier study supports a role for dopamine in thermal responses in the fly, which warrants deeper investigation in future studies. Specifically, future work should examine the effect of up- or down-regulation of dopamine  on temperature-dependent fitness.
Earlier studies have observed that cold acclimation can alter the glycerophospholipid composition of biomembranes ,. In line with these earlier findings, our results suggest that membrane sphingolipid metabolism might be important in how flies respond to changes in developmental temperature. Taken together with earlier studies, it seems likely that temperature has a major effect on cell membrane biochemistry. Interestingly, the role of sphingosine might go beyond membrane structure, as it is also involved in diverse signaling processes ,.
We also observed a significant effect of developmental temperature on zymosterol metabolism. While zymosterol is a key intermediate in cholesterol , it is likely that zymosterol in flies was bioaccumulated from yeast in fly media . This suggests that in D. melanogaster, developmental temperature modulates bioaccumulation of zymosterol from fly media, a hypothesis one could test by adding zymosterol to yeast-free holidic fly medium .
In contrast to developmental temperature, the effect of adult temperature on the metabolome was more modest (Table 3). This might be because adult flies were not exposed to the different temperatures for the same duration nor over different developmental stages, and developing fly larvae might be biochemically more malleable than adults. In our comparison of lifetime temperature exposure (combining samples with the same developmental and adult temperature), we observed changes in several metabolic pathways (Tables 4 and 5), at least one of which is of potentially interesting biological relevance. Spermine is a polyamine formed from spermidine and has been shown previously to mediate stress resistance in Drosophila .
These findings come with an important caveat. As with all metazoan species, the fruit fly metabolome has yet to be fully curated. We do not know how many metabolites are in the fly, and of those we have measured, most have unknown identities. Mummichog uses a metabolic pathways and network-based approach to find enriched pathways and annotate metabolites . For the majority of metabolites in the enriched pathways calculated by mummichog, we were able to find independent confirmation of their identities by mining the METLIN database using their m/z values (Additional file 8) . While validation of all putative matches is an important goal, it is beyond the scope of this current study, but highlights the need for a curated fly metabolome.
The fact that we confirm previous findings allows us to place greater confidence in our non-targeted, high-throughput analysis, as well as in the novel biological insights gained from our study. We report developmental temperature effects on several pathways not previously associated with temperature response. These include salvage pathways of adenine, hypoxanthine and their nucleosides, and acyl carrier protein metabolism. These findings suggest novel hypotheses for temperature adaptation that one could test through knock-down of fly genes coding for enzymes used in these pathways.
Effect of genotype
We were primarily interested in investigating the effects of developmental and adult temperature in the fly metabolome. While we found strong effects of temperature on the metabolome, genotype proved to affect an even greater number of metabolites than temperature. Interestingly, most of the pathways affected by genotype do not overlap with those modified by temperature. To our surprise, very few metabolites showed significant genotype-by-temperature interactions. Thus, our data suggest that the effect of temperature on the metabolome is consistent across genotypes. The lack of genotype-by-temperature effects might be due, in part, to the fact that we only included four genotypes in this study. Studies on a broader array of genotypes are needed to rule out the effect of genotype on temperature effects both on absolute metabolite levels and metabolite network structure.
Effect of temperature on network structure
To go beyond the identification of main effects, we performed a second level of analysis, examining the effect of temperature on the structure of metabolite co-expression networks and carrying out detailed analysis on such networks in response to temperature.
Differential co-expression analysis of metabolomic data is relatively novel ,, and until now has not been used to identify pathways associated with environmental perturbations. Our analysis has allowed us to develop novel hypotheses for metabolite modules that respond to temperature. Our network analysis identified both highly invariant as well as highly plastic metabolite modules. Our discovery of metabolite modules that were significantly preserved across the two temperature conditions suggests that there are not just individual metabolites, but large modules whose coordinate function is invariant in the face of at least mild temperature changes.
Importantly, we do not know whether the invariant correlation structure among modules is an adaptation that facilitates homeostasis, if the correlations are neutral, or if, in fact, the failure of those correlations to change in the face of different temperatures is actually maladaptive . To distinguish between these hypotheses, one would need to compare temperature-dependent fitness in groups in which the correlation structure was broken, either through dietary or genetic manipulation of two or more metabolites simultaneously, or by comparing genotypes (such as the DGRP ) in which levels of two or more metabolites deviate to a greater or lesser extent from that predicted from these correlations.
While the preserved modules found in this study cannot directly answer questions about network robustness, they can serve as a foundation for carrying out further experimental and computational studies to glean deeper insights into network robustness in response to temperature changes.
In summary, our results demonstrate that high resolution metabolomic profiles coupled with differential co-expression network analysis provides a powerful tool to understand the molecular response to shifts in temperature. Transcriptomic studies have used different methods to detect and quantify differentially co-expressed genes in a wide variety of settings, from cancer tissues ,, to Alzheimer’s disease , to cardiovascular disease . Here we show that network correlation analysis developed for gene expression data  can successfully be extended to analyze metabolomic data to address important biological questions.
Others have already noted that the metabolome offers a powerful intermediate step to link genetic and environmental variation to downstream phenotypes . This study illustrates how network approaches can further increase that explanatory power, suggesting novel hypotheses to be tested in future. Of course, we still have far to go—only a fraction of all metabolites are known, even in humans. As metabolomics catches up with genomic and proteomic technology, network analysis is likely to tell us much about organismal function.
All experiments were carried out using four genetically distinct and well characterized inbred strains (FlyBase 25180, 25184, 25189, and 25198) that are part of the Drosophila Genetic Reference Panel (DGRP) . Fly stocks were maintained in glass vials at 24°C on a 12/12 light–dark cycle at approximately 50% humidity. Flies were cultured on standard yeast-molasses-agar-cornmeal medium with propionic acid added as an anti-fungal agent.
For each strain, eight fresh vials were created, with ten males and ten females in each. After 48 hours of copulation and egg laying at 24°C, adults were discarded. The vials containing eggs were placed at either 18°C or 27°C until eclosion. Virgin collections for adult flies at developmental temperatures 27°C and 18°C were initiated at 10 and 19 days, respectively, following incubation at that temperature. Virgins were collected from each developmental temperature, and then distributed equally between the 18°C and 27°C incubators. This resulted in four temperature treatment groups: 18°C → 18°C, 18°C → 27°C, 27°C → 18°C, 27°C → 27°C. Adult flies were kept at this temperature for four to five days before being frozen in liquid nitrogen, and stored at −80°C.
Sample preparation and analysis were carried out as described previously . Briefly, acetonitrile extracts from frozen samples of three adult whole body D. melanogaster were analyzed in a dual column chromatography-mass spectrometry (LC-MS) platform. In this DC/LC-MS platform, the C18 (or ‘reverse phase chromatography’) column retains and separates chemicals with partial hydrophobic character, and the Anion Exchange (AE) column retains and separates negatively charged analytes. Fractionated samples, after electrospray ionization, were detected using a Fourier-transform mass spectrometer. Data were extracted as non-annotated mass/charge (m/z) features, column retention time, and ion intensity. Fly standards were also run alongside investigational samples. Specifically, we made a pooled reference sample and stored a large number of aliquots at −80°C to allow comparison of analytic behavior over long periods of time. The pooled reference samples were for the purposes of quality control (i.e., to ensure relative consistency among identical samples within days) and for quality assurance (i.e., to ensure consistent results between days). They did not contribute data to downstream statistical analysis.
Since the column chemistries are different, we analyzed data from the C18 and AE columns, separately. First, we applied various quality control procedures to the data, following . We included only those metabolites with a signal-to-noise ratio of at least 14. We then log-transformed the data and removed any analytes missing from more than 5% of male or more than 5% of female samples. After this, between 0.28% (AE) and 0.31% (C18) of all cells consisted of missing values (i.e., missing fewer than 1 in 300 m/z values in each case). To impute these missing values, we used the ‘EMimpute array’ method in the LSimpute algorithm .
Statistical analyses were carried out using the statistical package R (R Development Core team 2008) and Microsoft Excel 2013.
where y i is the intensity of metabolite i, μ is the mean, and ∊ is the error.
To identify differentially expressed metabolites for each factor in equation 1, we used a conservative false discovery rate (FDR) threshold of 0.01 . To investigate the effect of genotype on temperature treatment and its interactions, we carried out a likelihood ratio test .
Metabolic pathway enrichment and metabolite annotation
We sought to identify metabolic pathways that were enriched within the sets of differentially expressed metabolites. To do this, we used the program mummichog, which provides putative annotation for metabolites based on mass-charge ratios and carries out statistical tests for enrichment . Mummichog provides adjusted P-values, correcting for the fact that different m/z features can map to a single metabolite and single m/z values can map to different metabolites . Mummichog uses BioCyc (biocyc.org) as a source for metabolite ontologies where some of the fly metabolic pathways in BioCyc are based on experimental evidence while others are inferred. We chose an adjusted P-value cutoff of 0.01 to identify enriched metabolic pathways from Mummichog output.
We used the DiffCoEx algorithm to identify sets of metabolites or modules whose correlation structure changed significantly across experimental conditions . DiffCoEx runs on the Weighted Gene Correlation Network Analysis (WGCNA) platform, which is an R-based software package designed for network analysis . To calculate the significance of the co-expression differences within and between modules, we used permutation testing. Each dataset was permuted 1000 times for each temperature condition. The dispersion value (d) for each module and for each pair of modules were calculated. P-values were calculated as the proportion of permuted dispersion values greater than or equal to the original d . As before, we used mummichog to identify metabolic pathways that were enriched within specific modules. We also used a slight modification of DiffCoEx to detect significantly preserved modules in both male and female fly metabolomes. To do this, we modified the DiffCoEx algorithm to calculate a Similarity Matrix S instead of a dissimilarity matrix. Modules whose correlation structure changes between conditions were defined by hierarchical clustering on a difference matrix, D , where D = |(A − A ')/2|β/2, where A and A ’ are the adjacency matrixes under two conditions and β is the soft-thresholding parameter (from WGCNA ), set to β = 5 in all of our analyses. β is a positive integer and serves to transform the correlations so that weights of large correlation differences are given more emphasis compared to smaller differences which are less meaningful. Higher values of β increase statistical stringency by placing lower emphasis on smaller correlation differences. To identify modules with constant structure across treatments, we calculated a similarity matrix, S = (1 − D ) * |A | * |A '|. By this definition, pairs of metabolites whose correlations are similar and high (i.e., r ij,18 ≈ r ij,27 → 1, where r ij,t is the correlation coefficient between two metabolites i and j at temperature t), are weighted more heavily than pairs of metabolites whose correlations are similar but non-significant (i.e., r ij,18 ≈ r ij,27 → 0). At the limit where r ij,18 = r ij,27 , S ij = r 2.
To visualize representative modules, we used Cytoscape . Metabolites were defined as nodes, and the correlation between any two metabolites was treated as a weighted edge in Cytoscape.
For a detailed description of DiffCoEx parameters used in this study, see Additional file 9.
Ramkumar Hariharan and Jessica M Hoffman are co-first authors.
We thank all members of the Promislow lab for their helpful suggestions and critical comments on this manuscript.
- Gillooly JF, Brown JH, West GB, Savage VM, Charnov EL: Effects of size and temperature on metabolic rate. Science. 2001, 293 (5538): 2248-2251. 10.1126/science.1061967.View ArticlePubMedGoogle Scholar
- Pichaud N, Chatelain EH, Ballard JW, Tanguay R, Morrow G, Blier PU: Thermal sensitivity of mitochondrial metabolism in two distinct mitotypes of Drosophila simulans: evaluation of mitochondrial plasticity. J Exp Biol. 2010, 213 (10): 1665-1675. 10.1242/jeb.040261.View ArticlePubMedGoogle Scholar
- Koštál V, Šimek P, Zahradníčková H, Cimlová J, Štětina T: Conversion of the chill susceptible fruit fly larva (Drosophila melanogaster) to a freeze tolerant organism. Proc Natl Acad Sci U S A. 2012, 109 (9): 3270-3274. 10.1073/pnas.1119986109.PubMed CentralView ArticlePubMedGoogle Scholar
- Malmendal A, Overgaard J, Bundy JG, Sørensen JG, Nielsen NC, Loeschcke V, Holmstrup M: Metabolomic profiling of heat stress: hardening and recovery of homeostasis in Drosophila . Am J Physiol Regul Integr Comp Physiol. 2006, 291 (1): R205-R212. 10.1152/ajpregu.00867.2005.View ArticlePubMedGoogle Scholar
- Sejerkilde M, Sorensen JG, Loeschcke V: Effects of cold- and heat hardening on thermal resistance in Drosophila melanogaster . J Insect Physiol. 2003, 49 (8): 719-726. 10.1016/S0022-1910(03)00095-7.View ArticlePubMedGoogle Scholar
- Maynard Smith J: Acclimatization to high temperatures in inbred and outbred Drosophila subobscura. J. Genet. 1956, 54, 497–505. J Genet. 1956, 84 (1): 37-45. 10.1007/BF02715887.View ArticleGoogle Scholar
- Feder ME, Hofmann GE: Heat-shock proteins, molecular chaperones, and the stress response: evolutionary and ecological physiology. Annu Rev Physiol. 1999, 61: 243-282. 10.1146/annurev.physiol.61.1.243.View ArticlePubMedGoogle Scholar
- Arya R, Mallik M, Lakhotia SC: Heat shock genes - integrating cell survival and death. J Biosci. 2007, 32 (3): 595-610. 10.1007/s12038-007-0059-3.View ArticlePubMedGoogle Scholar
- Koštál V, Korbelová J, Rozsypal J, Zahradníčková H, Cimlová J, Tomčala A, Šimek P: Long-term cold acclimation extends survival time at 0 degrees C and modifies the metabolomic profiles of the larvae of the fruit fly Drosophila melanogaster . PLoS One. 2011, 6 (9): e25025-10.1371/journal.pone.0025025.PubMed CentralView ArticlePubMedGoogle Scholar
- Vesala L, Salminen TS, Koštál V, Zahradníčková H, Hoikkala A: Myo-inositol as a main metabolite in overwintering flies: seasonal metabolomic profiles and cold stress tolerance in a northern drosophilid fly. J Exp Biol. 2012, 215 (16): 2891-2897. 10.1242/jeb.069948.View ArticlePubMedGoogle Scholar
- Zhang J, Marshall KE, Westwood JT, Clark MS, Sinclair BJ: Divergent transcriptomic responses to repeated and single cold exposures in Drosophila melanogaster . J Exp Biol. 2011, 214 (Pt 23): 4021-4029. 10.1242/jeb.059535.View ArticlePubMedGoogle Scholar
- Sinclair BJ, Gibbs AG, Roberts SP: Gene transcription during exposure to, and recovery from, cold and desiccation stress in Drosophila melanogaster . Insect Mol Biol. 2007, 16 (4): 435-443. 10.1111/j.1365-2583.2007.00739.x.View ArticlePubMedGoogle Scholar
- Qin W, Neal SJ, Robertson RM, Westwood JT, Walker VK: Cold hardening and transcriptional change in Drosophila melanogaster . Insect Mol Biol. 2005, 14 (6): 607-613. 10.1111/j.1365-2583.2005.00589.x.View ArticlePubMedGoogle Scholar
- Colinet H, Overgaard J, Com E, Sørensen JG: Proteomic profiling of thermal acclimation in Drosophila melanogaster . Insect Biochem Mol Biol. 2013, 43 (4): 352-365. 10.1016/j.ibmb.2013.01.006.View ArticlePubMedGoogle Scholar
- Vermeulen CJ, Pedersen KS, Beck HC, Petersen J, Gagalova KK, Loeschcke V: Proteomic Characterization of Inbreeding-Related Cold Sensitivity in Drosophila melanogaster. Plos One 2013, 8(5):e62680.,Google Scholar
- Colinet H, Hoffmann AA: Comparing phenotypic effects and molecular correlates of developmental, gradual and rapid cold acclimation responses in Drosophila melanogaster . Funct Ecol. 2012, 26 (1): 84-93. 10.1111/j.1365-2435.2011.01898.x.View ArticleGoogle Scholar
- Pedersen KS, Kristensen TN, Loeschcke V, Petersen BO, Duus JØ, Nielsen NC, Malmendal A: Metabolomic signatures of inbreeding at benign and stressful temperatures in Drosophila melanogaster . Genetics. 2008, 180 (2): 1233-1243. 10.1534/genetics.108.089144.PubMed CentralView ArticlePubMedGoogle Scholar
- Sanders MM, Kon C: Glutamine is a powerful effector of heat-shock protein expression in Drosophila Kc cells. J Cell Physiol. 1991, 146 (1): 180-190. 10.1002/jcp.1041460123.View ArticlePubMedGoogle Scholar
- Wishart DS: Computational approaches to metabolomics. Methods Mol Biol. 2010, 593: 283-313. 10.1007/978-1-60327-194-3_14.View ArticlePubMedGoogle Scholar
- Jones DP, Park Y, Ziegler TR: Nutritional metabolomics: progress in addressing complexity in diet and health. Annu Rev Nutr. 2012, 32: 183-202. 10.1146/annurev-nutr-072610-145159.PubMed CentralView ArticlePubMedGoogle Scholar
- Park YH, Lee K, Soltow QA, Strobel FH, Brigham KL, Parker RE, Wilson ME, Sutliff RL, Mansfield KG, Wachtman LM, Ziegler TR, Jones DP: High-performance metabolic profiling of plasma from seven mammalian species for simultaneous environmental chemical surveillance and bioeffect monitoring. Toxicology. 2012, 295 (1–3): 47-55. 10.1016/j.tox.2012.02.007.PubMed CentralView ArticlePubMedGoogle Scholar
- Tesson BM, Breitling R, Jansen RC: DiffCoEx: a simple and sensitive method to find differentially coexpressed gene modules. BMC Bioinformatics 2010, 11:497.,Google Scholar
- Langfelder P, Horvath S: WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics. 2008, 9: 559-10.1186/1471-2105-9-559.PubMed CentralView ArticlePubMedGoogle Scholar
- Amar D, Safer H, Shamir R: Dissection of regulatory networks that are altered in disease via differential co-expression. PLoS Comput Biol. 2013, 9 (3): e1002955-10.1371/journal.pcbi.1002955.PubMed CentralView ArticlePubMedGoogle Scholar
- Fukushima A, Kusano M, Redestig H, Arita M, Saito K: Metabolomic correlation-network modules in Arabidopsis based on a graph-clustering approach. BMC Syst Biol. 2011, 5: 1-10.1186/1752-0509-5-1.PubMed CentralView ArticlePubMedGoogle Scholar
- Kotze HL, Armitage EG, Sharkey KJ, Allwood JW, Dunn WB, Williams KJ, Goodacre R: A novel untargeted metabolomics correlation-based network analysis incorporating human metabolic reconstructions. BMC Syst Biol. 2013, 7: 107-10.1186/1752-0509-7-107.PubMed CentralView ArticlePubMedGoogle Scholar
- Muller-Linow M, Weckwerth W, Hutt MT: Consistency analysis of metabolic correlation networks. BMC Syst Biol. 2007, 1: 44-10.1186/1752-0509-1-44.PubMed CentralView ArticlePubMedGoogle Scholar
- Oms-Oliu G, Hertog MLATM, Van de Poel B, Ampofo-Asiama J, Geeraerd AH, Nicolai BM: Metabolic characterization of tomato fruit during preharvest development, ripening, and postharvest shelf-life. Postharvest Biology and Technology. 2011, 62 (1): 7-16. 10.1016/j.postharvbio.2011.04.010.View ArticleGoogle Scholar
- Barkai N, Leibler S: Robustness in simple biochemical networks. Nature. 1997, 387 (6636): 913-917. 10.1038/43199.View ArticlePubMedGoogle Scholar
- Gibson G: Robustness and evolvability in living systems. Science. 2005, 310 (5746): 237-237. 10.1126/science.1119501.View ArticleGoogle Scholar
- Mackay TF, Richards S, Stone EA, Barbadilla A, Ayroles JF, Zhu D, Casillas S, Han Y, Magwire MM, Cridland JM, Richardson MF, Anholt RR, Barrón M, Bess C, Blankenburg KP, Carbone MA, Castellano D, Chaboub L, Duncan L, Harris Z, Javaid M, Jayaseelan JC, Jhangiani SN, Jordan KW, Lara F, Lawrence F, Lee SL, Librado P, Linheiro RS, Lyman RF, et al: The Drosophila melanogaster genetic reference panel. Nature. 2012, 482 (7384): 173-178. 10.1038/nature10811.PubMed CentralView ArticlePubMedGoogle Scholar
- Li S, Park Y, Duraisingham S, Strobel FH, Khan N, Soltow QA, Jones DP, Pulendran B: Predicting network activity from high throughput metabolomics. PLoS Comput Biol. 2013, 9 (7): e1003123-10.1371/journal.pcbi.1003123.PubMed CentralView ArticlePubMedGoogle Scholar
- Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T: Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003, 13 (11): 2498-2504. 10.1101/gr.1239303.PubMed CentralView ArticlePubMedGoogle Scholar
- Hoffmann AA: Physiological climatic limits in Drosophila: patterns and implications. J Exp Biol. 2010, 213 (6): 870-880. 10.1242/jeb.037630.View ArticlePubMedGoogle Scholar
- Clark MS, Worland MR: How insects survive the cold: molecular mechanisms-a review. J Comp Physiol B. 2008, 178 (8): 917-933. 10.1007/s00360-008-0286-4.View ArticlePubMedGoogle Scholar
- Feder ME, Blair N, Figueras H: Natural thermal stress and heat-shock protein expression in Drosophila larvae and pupae. Funct Ecol. 1997, 11 (1): 90-100. 10.1046/j.1365-2435.1997.00060.x.View ArticleGoogle Scholar
- Dillon ME, Wang G, Garrity PA, Huey RB: Thermal preference in Drosophila . J Therm Biol. 2009, 34 (3): 109-119. 10.1016/j.jtherbio.2008.11.007.PubMed CentralView ArticlePubMedGoogle Scholar
- Czarnoleski M, Cooper BS, Kierat J, Angilletta MJ Jr: Flies developed small bodies and small cells in warm and in thermally fluctuating environments. J Exp Biol. 2013, 216 (Pt 15): 2896-2901. 10.1242/jeb.083535.PubMed CentralView ArticlePubMedGoogle Scholar
- Ueno T, Tomita J, Kume S, Kume K: Dopamine Modulates Metabolic Rate and Temperature Sensitivity in Drosophila melanogaster. Plos One 2012, 7(2):e31513.,Google Scholar
- Hirsh J, Riemensperger T, Coulom H, Iché M, Coupar J, Birman S: Roles of dopamine in circadian rhythmicity and extreme light sensitivity of circadian entrainment. Curr Biol. 2010, 20 (3): 209-214. 10.1016/j.cub.2009.11.037.View ArticlePubMedGoogle Scholar
- Alexander GJ, Schwenk E: Studies on biosynthesis of cholesterol. IX. Zymosterol as a precursor of cholesterol. Arch Biochem Biophys. 1957, 66 (2): 381-387. 10.1016/S0003-9861(57)80013-7.View ArticlePubMedGoogle Scholar
- Gault CR, Obeid LM, Hannun YA: An overview of sphingolipid metabolism: from synthesis to breakdown. Adv Exp Med Biol. 2010, 688: 1-23. 10.1007/978-1-4419-6741-1_1.PubMed CentralView ArticlePubMedGoogle Scholar
- Carvalho M, Sampaio JL, Palm W, Brankatschk M, Eaton S, Shevchenko A: Effects of diet and development on the Drosophila lipidome. Mol Syst Biol. 2012, 8: 600-10.1038/msb.2012.29.PubMed CentralView ArticlePubMedGoogle Scholar
- Piper MD, Blanc E, Leitão-Gonçalves R, Yang M, He X, Linford NJ, Hoddinott MP, Hopfen C, Soultoukis GA, Niemeyer C, Kerr F, Pletcher SD, Ribeiro C, Partridge L: A holidic medium for Drosophila melanogaster . Nat Methods. 2014, 11 (1): 100-105. 10.1038/nmeth.2731.View ArticlePubMedGoogle Scholar
- Sigrist SJ, Carmona-Gutierrez D, Gupta VK, Bhukel A, Mertel S, Eisenberg T, Madeo F: Spermidine-triggered autophagy ameliorates memory during aging. Autophagy. 2014, 10 (1): 178-179. 10.4161/auto.26918.View ArticlePubMedGoogle Scholar
- Tautenhahn R, Cho K, Uritboonthai W, Zhu Z, Patti GJ, Siuzdak G: An accelerated workflow for untargeted metabolomics using the METLIN database. Nat Biotechnol. 2012, 30 (9): 826-828. 10.1038/nbt.2348.PubMed CentralView ArticlePubMedGoogle Scholar
- Castro C, Krumsiek J, Lehrbach NJ, Murfitt SA, Miska EA, Griffin JL: A study of Caenorhabditis elegans DAF-2 mutants by metabolomics and differential correlation networks. Mol Biosyst. 2013, 9 (7): 1632-1642. 10.1039/c3mb25539e.View ArticlePubMedGoogle Scholar
- Ghalambor CK, McKay JK, Carroll SP, Reznick DN: Adaptive versus non-adaptive phenotypic plasticity and the potential for contemporary adaptation in new environments. Funct Ecol. 2007, 21 (3): 394-407. 10.1111/j.1365-2435.2007.01283.x.View ArticleGoogle Scholar
- Bockmayr M, Klauschen F, Györffy B, Denkert C, Budczies J: New network topology approaches reveal differential correlation patterns in breast cancer. Systems Biology 2013, 7:78.,Google Scholar
- Bhattacharyya M, Bandyopadhyay S: Studying the differential co-expression of microRNAs reveals significant role of white matter in early Alzheimer’s progression. Mol Biosyst. 2013, 9 (3): 457-466. 10.1039/c2mb25434d.View ArticlePubMedGoogle Scholar
- Rotival M, Petretto E: Leveraging gene co-expression networks to pinpoint the regulation of complex traits and disease, with a focus on cardiovascular traits. Brief Funct Genomics. 2014, 13 (1): 66-78. 10.1093/bfgp/elt030.View ArticlePubMedGoogle Scholar
- Dumas ME: Metabolome 2.0: quantitative genetics and network biology of metabolic phenotypes. Mol Biosyst. 2012, 8 (10): 2494-2502. 10.1039/c2mb25167a.View ArticlePubMedGoogle Scholar
- Hoffman JM, Soltow QA, Li S, Sidik A, Jones DP, Promislow DE: Effects of age, sex, and genotype on high-sensitivity metabolomic profiles in the fruit fly, Drosophila melanogaster. Aging Cell 2014, 13(4):596–604.,Google Scholar
- Bo TH, Dysvik B, Jonassen I: LSimpute: accurate estimation of missing values in microarray data with least squares methods. Nucleic Acids Res. 2004, 32 (3): e34-10.1093/nar/gnh026.PubMed CentralView ArticlePubMedGoogle Scholar
- Benjamini Y, Hochberg Y: Controlling the false discovery rate - a practical and powerful approach to multiple testing. Journal of the Royal Statistical Society Series B-Methodological. 1995, 57 (1): 289-300.Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. 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.