Relative importance of intensity and spectrum of artificial light at night in disrupting behavior of a nocturnal rodent.

The influence of light spectral properties on circadian rhythms is of substantial interest to laboratory-based investigation of the circadian system and to field-based understanding of the effects of artificial light at night. The tradeoffs between intensity and spectrum regarding masking behaviors are largely unknown, even for well-studied organisms. We used a custom LED illumination system to document the response of wild type house mice (Mus musculus) to 1-hr nocturnal exposure of all combinations of four intensity levels (0.01, 0.5, 5, and 50 lx) and three correlated color temperatures (CCT; 1750, 1950, and 3000 K). Higher intensities of light (50 lx) suppressed cage activity substantially, and consistently more for the higher CCT light (91% for 3000 K; 53% for 1750 K). At the lower intensities (0.01 lx), mean activity was increased, with the greatest increases for the lowest CCT (12.3% increase at 1750 K; 3% increase at 3000 K). Multiple linear regression confirmed the influence of both CCT (p<0.001) and intensity (p<0.001) on changes in activity (r2=0.66, F9,171=3.33; p<0.001) with the scaled effect size of intensity 3.6 times greater than CCT. Activity suppression was significantly lower for male than female mice (p<0.0001). Assessment of light-evoked cFos expression in the suprachiasmatic nucleus at 50 lx showed no significant difference between high and low CCT exposure. The significant differences by spectral composition illustrate a need to account for light spectrum in circadian studies of behavior and confirm that spectral controls can mitigate some, but certainly not all, of the effects of light pollution on species in the wild.

Cellular age explains variation in age-related cell-to-cell transcriptome variability.

Organs and tissues age at different rates within a single individual. Such asynchrony in aging has been widely observed at multiple levels, from functional hallmarks, such as anatomical structures and physiological processes, to molecular endophenotypes, such as the transcriptome and metabolome. However, we lack a conceptual framework to understand why some components age faster than others. Just as demographic models explain why aging evolves, here we test the hypothesis that demographic differences among cell types, determined by cell-specific differences in turnover rate, can explain why the transcriptome shows signs of aging in some cell types but not others. Through analysis of mouse single-cell transcriptome data across diverse tissues and ages, we find that cellular age explains a large proportion of the variation in the age-related increase in transcriptome variance. We further show that long-lived cells are characterized by relatively high expression of genes associated with proteostasis and that the transcriptome of long-lived cells shows greater evolutionary constraint than short-lived cells. In contrast, in short-lived cell types, the transcriptome is enriched for genes associated with DNA repair. Based on these observations, we develop a novel heuristic model that explains how and why aging rates differ among cell types.

Members of the CUGBP Elav-like family of RNA-binding proteins are expressed in distinct populations of primary sensory neurons.

Primary sensory dorsal root ganglia (DRG) neurons are diverse, with distinct populations that respond to specific stimuli. Previously, we observed that functionally distinct populations of DRG neurons express mRNA transcript variants with different 3' untranslated regions (3'UTRs). 3'UTRs harbor binding sites for interaction with RNA-binding proteins (RBPs) for transporting mRNAs to subcellular domains, modulating transcript stability, and regulating the rate of translation. In the current study, analysis of publicly available single-cell RNA-sequencing data generated from adult mice revealed that 17 3'UTR-binding RBPs were enriched in specific populations of DRG neurons. This included four members of the CUG triplet repeat (CUGBP) Elav-like family (CELF): CELF2 and CELF4 were enriched in peptidergic, CELF6 in both peptidergic and nonpeptidergic, and CELF3 in tyrosine hydroxylase-expressing neurons. Immunofluorescence studies confirmed that 60% of CELF4+ neurons are small-diameter C fibers and 33% medium-diameter myelinated (likely Aδ) fibers and showed that CELF4 is distributed to peripheral termini. Coexpression analyses using transcriptomic data and immunofluorescence revealed that CELF4 is enriched in nociceptive neurons that express GFRA3, CGRP, and the capsaicin receptor TRPV1. Reanalysis of published transcriptomic data from macaque DRG revealed a highly similar distribution of CELF members, and reanalysis of single-nucleus RNA-sequencing data derived from mouse and rat DRG after sciatic injury revealed differential expression of CELFs in specific populations of sensory neurons. We propose that CELF RBPs may regulate the fate of mRNAs in populations of nociceptors, and may play a role in pain and/or neuronal regeneration following nerve injury.

Outcomes after Genicular Artery Embolization Vary According to the Radiographic Severity of Osteoarthritis: Results from a Prospective Single-Center Study.

This prospective single-center study compared outcomes after genicular artery embolization (GAE) in people with mild radiographic knee osteoarthritis (OA) (n = 11) versus moderate to severe radiographic knee OA (n = 22). The mild OA group was older and had experienced symptoms for less time (P < .05). All participants underwent complete embolization of neovessels from genicular arteries. The primary outcome was the proportion of responders at 6 months according to predetermined improvements in pain, function, and/or global change. Results indicated that a higher proportion of participants (n = 9, 81.8%) with mild OA fulfilled responder criteria after treatment compared with people with moderate to severe OA (n = 8, 36.4%) (P = .014). Secondary outcomes for pain, quality of life, and global change were also better in the mild OA group (P < .05). No serious adverse events occurred, including the absence of osteonecrosis on magnetic resonance imaging. The study demonstrated that outcomes after GAE varied according to the severity of baseline radiographic OA.

Using Drosophila to identify naturally occurring genetic modifiers of amyloid beta 42- and tau-induced toxicity.

Alzheimer's disease is characterized by 2 pathological proteins, amyloid beta 42 and tau. The majority of Alzheimer's disease cases in the population are sporadic and late-onset Alzheimer's disease, which exhibits high levels of heritability. While several genetic risk factors for late-onset Alzheimer's disease have been identified and replicated in independent studies, including the ApoE ε4 allele, the great majority of the heritability of late-onset Alzheimer's disease remains unexplained, likely due to the aggregate effects of a very large number of genes with small effect size, as well as to biases in sample collection and statistical approaches. Here, we present an unbiased forward genetic screen in Drosophila looking for naturally occurring modifiers of amyloid beta 42- and tau-induced ommatidial degeneration. Our results identify 14 significant SNPs, which map to 12 potential genes in 8 unique genomic regions. Our hits that are significant after genome-wide correction identify genes involved in neuronal development, signal transduction, and organismal development. Looking more broadly at suggestive hits (P < 10-5), we see significant enrichment in genes associated with neurogenesis, development, and growth as well as significant enrichment in genes whose orthologs have been identified as significantly or suggestively associated with Alzheimer's disease in human GWAS studies. These latter genes include ones whose orthologs are in close proximity to regions in the human genome that are associated with Alzheimer's disease, but where a causal gene has not been identified. Together, our results illustrate the potential for complementary and convergent evidence provided through multitrait GWAS in Drosophila to supplement and inform human studies, helping to identify the remaining heritability and novel modifiers of complex diseases.

Genicular artery embolization for early-stage knee osteoarthritis: results from a triple-blind single-centre randomized controlled trial.

This study investigated the effects of transcatheter arterial embolization (TAE) on pain, function, and quality of life in people with early-stage symptomatic knee osteoarthritis (OA) compared to a sham procedure. A total of 59 participants with symptomatic Kellgren-Lawrence grade 2 knee OA were randomly allocated to TAE or a sham procedure. The intervention group underwent TAE of one or more genicular arteries. The control group received a blinded sham procedure. The primary outcome was knee pain at 12 months according to the Knee injury and Osteoarthritis Outcome Score (KOOS) pain scale. Secondary outcomes included self-reported function and quality of life (KOOS, EuroQol five-dimension five-level questionnaire (EQ-5D-5L)), self-reported Global Change, six-minute walk test, 30-second chair stand test, and adverse events. Subgroup analyses compared participants who received complete embolization of all genicular arteries (as distinct from embolization of some arteries) (n = 17) with the control group (n = 29) for KOOS and Global Change scores at 12 months. Continuous variables were analyzed with quantile regression, adjusting for baseline scores. Dichotomized variables were analyzed with chi-squared tests. Overall, 58 participants provided questionnaire data at 12 months. No significant differences were found for the primary and secondary outcomes, with both groups improving following the procedure. At 12 months, KOOS pain scores improved by 41.3% and 29.4% in the intervention and control groups, respectively. No adverse events occurred. Subgroup analysis indicated that the complete embolization group had significantly better KOOS Sports and Recreation, KOOS Quality of Life, and Global Change scores than the control group; 76.5% of participants who received complete embolization reporting being moderately or much better compared to 37.9% of the control group. TAE might produce benefits above placebo, but only when complete embolization of all genicular arteries is performed. Further comparative studies are required before definitive conclusions regarding the effectiveness of TAE can be made. Level of evidence: I.

Longitudinal RNA Sequencing of Skin and DRG Neurons in Mice with Paclitaxel-Induced Peripheral Neuropathy.

Paclitaxel-induced peripheral neuropathy is a condition of nerve degeneration induced by chemotherapy, which afflicts up to 70% of treated patients. Therapeutic interventions are unavailable due to an incomplete understanding of the underlying mechanisms. We previously discovered that major physiological changes in the skin underlie paclitaxel-induced peripheral neuropathy in zebrafish and rodents. The precise molecular mechanisms are only incompletely understood. For instance, paclitaxel induces the upregulation of MMP-13, which, when inhibited, prevents axon degeneration. To better understand other gene regulatory changes induced by paclitaxel, we induced peripheral neuropathy in mice following intraperitoneal injection either with vehicle or paclitaxel every other day four times total. Skin and dorsal root ganglion neurons were collected based on distinct behavioural responses categorised as "pain onset" (d4), "maximal pain" (d7), "beginning of pain resolution" (d11), and "recovery phase" (d23) for comparative longitudinal RNA sequencing. The generated datasets validate previous discoveries and reveal additional gene expression changes that warrant further validation with the goal to aid in the development of drugs that prevent or reverse paclitaxel-induced peripheral neuropathy.

In search of a Drosophila core cellular network with single-cell transcriptome data.

Along with specialized functions, cells of multicellular organisms also perform essential functions common to most if not all cells. Whether diverse cells do this by using the same set of genes, interacting in a fixed coordinated fashion to execute essential functions, or a subset of genes specific to certain cells, remains a central question in biology. Here, we focus on gene coexpression to search for a core cellular network across a whole organism. Single-cell RNA-sequencing measures gene expression of individual cells, enabling researchers to discover gene expression patterns that contribute to the diversity of cell functions. Current efforts to study cellular functions focus primarily on identifying differentially expressed genes across cells. However, patterns of coexpression between genes are probably more indicative of biological processes than are the expression of individual genes. We constructed cell-type-specific gene coexpression networks using single-cell transcriptome datasets covering diverse cell types from the fruit fly, Drosophila melanogaster. We detected a set of highly coordinated genes preserved across cell types and present this as the best estimate of a core cellular network. This core is very small compared with cell-type-specific gene coexpression networks and shows dense connectivity. Gene members of this core tend to be ancient genes and are enriched for those encoding ribosomal proteins. Overall, we find evidence for a core cellular network in diverse cell types of the fruit fly. The topological, structural, functional, and evolutionary properties of this core indicate that it accounts for only a minority of essential functions.

Thirty second chair stand test: Test-retest reliability, agreement and minimum detectable change in people with early-stage knee osteoarthritis.

BACKGROUND AND PURPOSE: To determine intra-session test-retest reliability, agreement and minimum detectable change (MDC) of the 30 CST across three tests in people with knee osteoarthritis (OA). METHODS: A test-retest reliability study was performed with 93 people with mild radiological knee OA. Participants were asked to complete three attempts of the 30 CST 1-2 min apart according to a standardised protocol. Participants completed three attempts on two occasions: baseline and 6 months later. Change between tests within each session was assessed with ANOVA's and post-hoc t-tests. Reliability was assessed with intra-class correlation coefficients (ICC[2,1] ). Measurement error was expressed as MDC for an individual (MDCind ) and a group (MDCgroup ). Floor effects were considered present if more than 15% of participants scored zero for a test. RESULTS: Scores increased by 0.5 and 0.8 stands between the first and second test (p < 0.05) at the baseline and 6-month assessments respectively, and then stabilised between the second the third tests at the baseline assessment (p > 0.05) or decreased (0.3 stands) at the 6-month assessment (p < 0.05). Scores demonstrated excellent reliability (ICCs >0.9). MDCind was approximately 2.5 stands and MDCgroup was 0.3-0.4 stands. No floor effects were apparent. DISCUSSION: The 30CST demonstrated a practice effect between the first and second tests, which was no longer apparent by the third test. Despite this, scores demonstrated excellent intra-session reliability. MDC estimates provide clinicians and researchers with the smallest change that can be detected by the instrument beyond measurement error for individuals and groups in community-dwelling adults with knee OA.

Correction: Genetic and metabolomic architecture of variation in diet restriction-mediated lifespan extension in Drosophila.

[This corrects the article DOI: 10.1371/journal.pgen.1008835.].

The metabolome as a biomarker of aging in Drosophila melanogaster.

Many biomarkers have been shown to be associated not only with chronological age but also with functional measures of biological age. In human populations, it is difficult to show whether variation in biological age is truly predictive of life expectancy, as such research would require longitudinal studies over many years, or even decades. We followed adult cohorts of 20 Drosophila Genetic Reference Panel (DGRP) strains chosen to represent the breadth of lifespan variation, obtain estimates of lifespan, baseline mortality, and rate of aging, and associate these parameters with age-specific functional traits including fecundity and climbing activity and with age-specific targeted metabolomic profiles. We show that activity levels and metabolome-wide profiles are strongly associated with age, that numerous individual metabolites show a strong association with lifespan, and that the metabolome provides a biological clock that predicts not only sample age but also future mortality rates and lifespan. This study with 20 genotypes and 87 metabolites, while relatively small in scope, establishes strong proof of principle for the fly as a powerful experimental model to test hypotheses about biomarkers and aging and provides further evidence for the potential value of metabolomic profiles as biomarkers of aging.

Modular Evolution of the Drosophila Metabolome.

Comparative phylogenetic studies offer a powerful approach to study the evolution of complex traits. Although much effort has been devoted to the evolution of the genome and to organismal phenotypes, until now relatively little work has been done on the evolution of the metabolome, despite the fact that it is composed of the basic structural and functional building blocks of all organisms. Here we explore variation in metabolite levels across 50 My of evolution in the genus Drosophila, employing a common garden design to measure the metabolome within and among 11 species of Drosophila. We find that both sex and age have dramatic and evolutionarily conserved effects on the metabolome. We also find substantial evidence that many metabolite pairs covary after phylogenetic correction, and that such metabolome coevolution is modular. Some of these modules are enriched for specific biochemical pathways and show different evolutionary trajectories, with some showing signs of stabilizing selection. Both observations suggest that functional relationships may ultimately cause such modularity. These coevolutionary patterns also differ between sexes and are affected by age. We explore the relevance of modular evolution to fitness by associating modules with lifespan variation measured in the same common garden. We find several modules associated with lifespan, particularly in the metabolome of older flies. Oxaloacetate levels in older females appear to coevolve with lifespan, and a lifespan-associated module in older females suggests that metabolic associations could underlie 50 My of lifespan evolution.

Triple Targeting of Breast Tumors Driven by Hormonal Receptors and HER2.

Breast cancers that express hormonal receptors (HR) and HER2 display resistance to targeted therapy. Tumor-promotional signaling from the HER2 and estrogen receptor (ER) pathways converges at the cyclin D1 and cyclin-dependent kinases (CDK) 4 and 6 complex, which drives cell-cycle progression and development of therapeutic resistance. Therefore, we hypothesized that co-targeting of ER, HER2, and CDK4/6 may result in improved tumoricidal activity and suppress drug-resistant subclones that arise on therapy. We tested the activity of the triple targeted combination therapy with tucatinib (HER2 small-molecule inhibitor), palbociclib (CKD4/6 inhibitor), and fulvestrant (selective ER degrader) in HR+/HER2+ human breast tumor cell lines and xenograft models. In addition, we evaluated whether triple targeted combination prevents growth of tucatinib or palbociclib-resistant subclones in vitro and in vivo Triple targeted combination significantly reduced HR+/HER2+ tumor cell viability, clonogenic survival, and in vivo growth. Moreover, survival of HR+/HER2+ cells that were resistant to the third drug in the regimen was reduced by the other two drugs in combination. We propose that a targeted triple combination approach will be clinically effective in the treatment of otherwise drug-resistant tumors, inducing robust responses in patients.

Evolution of natural lifespan variation and molecular strategies of extended lifespan in yeast.

To understand the genetic basis and selective forces acting on longevity, it is useful to examine lifespan variation among closely related species, or ecologically diverse isolates of the same species, within a controlled environment. In particular, this approach may lead to understanding mechanisms underlying natural variation in lifespan. Here, we analyzed 76 ecologically diverse wild yeast isolates and discovered a wide diversity of replicative lifespan (RLS). Phylogenetic analyses pointed to genes and environmental factors that strongly interact to modulate the observed aging patterns. We then identified genetic networks causally associated with natural variation in RLS across wild yeast isolates, as well as genes, metabolites, and pathways, many of which have never been associated with yeast lifespan in laboratory settings. In addition, a combined analysis of lifespan-associated metabolic and transcriptomic changes revealed unique adaptations to interconnected amino acid biosynthesis, glutamate metabolism, and mitochondrial function in long-lived strains. Overall, our multiomic and lifespan analyses across diverse isolates of the same species shows how gene-environment interactions shape cellular processes involved in phenotypic variation such as lifespan.

Genetic Analysis of Methyl Anthranilate, Mesifurane, Linalool, and Other Flavor Compounds in Cultivated Strawberry (Fragaria × ananassa).

The cultivated strawberry (Fragaria × ananassa) is an economically important fruit crop that is intensively bred for improved sensory qualities. The diversity of fruit flavors and aromas in strawberry results mainly from the interactions of sugars, acids, and volatile organic compounds (VOCs) that are derived from diverse biochemical pathways influenced by the expression of many genes. This study integrates multiomic analyses to identify QTL and candidate genes for multiple aroma compounds in a complex strawberry breeding population. Novel fruit volatile QTL was discovered for methyl anthranilate, methyl 2-hexenoate, methyl 2-methylbutyrate, mesifurane, and a shared QTL on Chr 3 was found for nine monoterpene and sesquiterpene compounds, including linalool, 3-carene, ß-phellandrene, α-limonene, linalool oxide, nerolidol, α-caryophellene, α-farnesene, and ß-farnesene. Fruit transcriptomes from a subset of 64 individuals were used to support candidate gene identification. For methyl esters including the grape-like methyl anthranilate, a novel ANTHANILIC ACID METHYL TRANSFERASE-like gene was identified. Two mesifurane QTL correspond with the known biosynthesis gene O-METHYL TRANSFERASE 1 and a novel FURANEOL GLUCOSYLTRANSFERASE. The shared terpene QTL contains multiple fruit-expressed terpenoid pathway-related genes including NEROLIDOL SYNTHASE 1 (FanNES1). The abundance of linalool and other monoterpenes is partially governed by a co-segregating expression-QTL (eQTL) for FanNES1 transcript variation, and there is additional evidence for quantitative effects from other terpenoid-pathway genes in this narrow genomic region. These QTLs present new opportunities in breeding for improved flavor in commercial strawberry.

Gene-environment interaction in molar-incisor hypomineralization.

Molar incisor hypomineralization (MIH) is an enamel condition characterized by lesions ranging in color from white to brown which present rapid caries progression, and mainly affects permanent first molars and incisors. These enamel defects usually occur when there are disturbances during the mineralization or maturation stage of amelogenesis. Both genetic and environmental factors have been suggested to play roles in MIH's development, but no conclusive risk factors have shown the source of the disease. During head and neck development, the interferon regulatory factor 6 (IRF6) gene is involved in the structure formation of the oral and maxillofacial regions, and the transforming growth factor alpha (TGFA) is an essential cell regulator, acting during proliferation, differentiation, migration and apoptosis. In this present study, it was hypothesized that these genes interact and contribute to predisposition of MIH. Environmental factors affecting children that were 3 years of age or older were also hypothesized to play a role in the disease etiology. Those factors included respiratory issues, malnutrition, food intolerance, infection of any sort and medication intake. A total of 1,065 salivary samples from four different cohorts were obtained, and DNA was extracted from each sample and genotyped for nine different single nucleotide polymorphisms. Association tests and logistic regression implemented in PLINK were used for analyses. A potential interaction between TGFA rs930655 with all markers tested in the cohort from Turkey was identified. These interactions were not identified in the remaining cohorts. Associations (p<0.05) between the use of medication after three years of age and MIH were also found, suggesting that conditions acquired at the age children start to socialize might contribute to the development of MIH.

Genicular Artery Embolization to Improve Pain and Function in Early-Stage Knee Osteoarthritis-24-Month Pilot Study Results.

Ten patients with painful early-stage knee osteoarthritis underwent transcatheter arterial embolization of neoangiogenic branches from genicular arteries. At 12 months, 6 patients had responded to treatment according to study criteria. Median pain, function, and quality-of-life scores had improved by 15.4%, 21.3%, and 100%, respectively, for the 10 patients. Median 30-second chair stand test and 6-minute walk test scores had improved by 43% and 26%, respectively. No major adverse events occurred. Improvements following transcatheter arterial embolization provide evidence that embolization might be an effective treatment technique for people with early-stage knee osteoarthritis, although more rigorous evaluation is required.

Genetic and metabolomic architecture of variation in diet restriction-mediated lifespan extension in Drosophila.

In most organisms, dietary restriction (DR) increases lifespan. However, several studies have found that genotypes within the same species vary widely in how they respond to DR. To explore the mechanisms underlying this variation, we exposed 178 inbred Drosophila melanogaster lines to a DR or ad libitum (AL) diet, and measured a panel of 105 metabolites under both diets. Twenty four out of 105 metabolites were associated with the magnitude of the lifespan response. These included proteinogenic amino acids and metabolites involved in α-ketoglutarate (α-KG)/glutamine metabolism. We confirm the role of α-KG/glutamine synthesis pathways in the DR response through genetic manipulations. We used covariance network analysis to investigate diet-dependent interactions between metabolites, identifying the essential amino acids threonine and arginine as "hub" metabolites in the DR response. Finally, we employ a novel metabolic and genetic bipartite network analysis to reveal multiple genes that influence DR lifespan response, some of which have not previously been implicated in DR regulation. One of these is CCHa2R, a gene that encodes a neuropeptide receptor that influences satiety response and insulin signaling. Across the lines, variation in an intronic single nucleotide variant of CCHa2R correlated with variation in levels of five metabolites, all of which in turn were correlated with DR lifespan response. Inhibition of adult CCHa2R expression extended DR lifespan of flies, confirming the role of CCHa2R in lifespan response. These results provide support for the power of combined genomic and metabolomic analysis to identify key pathways underlying variation in this complex quantitative trait.

Resistance to state transitions in responsiveness is differentially modulated by different volatile anaesthetics in male mice.

BACKGROUND: Recent studies point to a fundamental distinction between population-based and individual-based anaesthetic pharmacology. At the population level, anaesthetic potency is defined as the relationship between drug concentration and the likelihood of response to a stimulus. At the individual level, even when the anaesthetic concentration is held constant, fluctuations between the responsive and unresponsive states are observed. Notably, these spontaneous fluctuations exhibit resistance to state transitions Rst. Therefore, the response probability in each individual depends not just upon the drug concentration, but also upon responses to previous stimuli. Here, we hypothesise that Rst is distinct from drug potency and is differentially modulated by different anaesthetics. METHODS: Adult (14-24 weeks old) C57BL/6J male mice (n=60) were subjected to repeated righting reflex (RR) assays at equipotent steady-state concentrations of isoflurane (0.6 vol%), sevoflurane (1.0 vol%), and halothane (0.4 vol%). RESULTS: Fluctuations in RR were observed for all tested anaesthetics. Analysis of these fluctuations revealed that Rst was differentially modulated by different anaesthetics (F[2, 56.01]=49.59; P<0.0001). Fluctuations in RR were modelled using a stochastic dynamical system. This analysis confirmed that the amount of noise that drives behavioural state transitions depends on the anaesthetic agent (F[2, 42.86]=16.72; P<0.0001). CONCLUSIONS: Whilst equipotent doses of distinct anaesthetics produce comparable population response probabilities, they engage dramatically different dynamics in each individual animal. This manifests as a differential aggregate propensity to exhibit state transitions. Thus, resistance to state transitions is a fundamentally distinct, novel measure of individualised anaesthetic pharmacology.

The metabolome as a link in the genotype-phenotype map for peroxide resistance in the fruit fly, Drosophila melanogaster.

BACKGROUND: Genetic association studies that seek to explain the inheritance of complex traits typically fail to explain a majority of the heritability of the trait under study. Thus, we are left with a gap in the map from genotype to phenotype. Several approaches have been used to fill this gap, including those that attempt to map endophenotype such as the transcriptome, proteome or metabolome, that underlie complex traits. Here we used metabolomics to explore the nature of genetic variation for hydrogen peroxide (H2O2) resistance in the sequenced inbred Drosophila Genetic Reference Panel (DGRP). RESULTS: We first studied genetic variation for H2O2 resistance in 179 DGRP lines and along with identifying the insulin signaling modulator u-shaped and several regulators of feeding behavior, we estimate that a substantial amount of phenotypic variation can be explained by a polygenic model of genetic variation. We then profiled a portion of the aqueous metabolome in subsets of eight 'high resistance' lines and eight 'low resistance' lines. We used these lines to represent collections of genotypes that were either resistant or sensitive to the stressor, effectively modeling a discrete trait. Across the range of genotypes in both populations, flies exhibited surprising consistency in their metabolomic signature of resistance. Importantly, the resistance phenotype of these flies was more easily distinguished by their metabolome profiles than by their genotypes. Furthermore, we found a metabolic response to H2O2 in sensitive, but not in resistant genotypes. Metabolomic data further implicated at least two pathways, glycogen and folate metabolism, as determinants of sensitivity to H2O2. We also discovered a confounding effect of feeding behavior on assays involving supplemented food. CONCLUSIONS: This work suggests that the metabolome can be a point of convergence for genetic variation influencing complex traits, and can efficiently elucidate mechanisms underlying trait variation.