Understanding the genotype -phenotype map and how variation at different levels of biological organization are associated are central topics in modern biology. Fast developments in sequencing technologies and other molecular omic tools enable researchers to obtain detailed information on variation at DNA level and on intermediate endophenotypes, such as RNA, proteins, metabolites. This facilitates our understanding of the link between genotypes and molecular and functional organismal phenotypes. Here, we use the Drosophila Genetic Reference Panel and nuclear magnetic resonance (NMR) metabolomics to investigate the ability of the metabolome to predict organismal phenotypes. We do this by performing NMR metabolomics on four replicate pools of male flies from each of 170 isogenic lines. We show first that metabolite profiles are variable among the investigated lines and that this variation is highly heritable. Second, we identify genes associated with metabolome variation. Third, we show that for four of five traits related to behaviour and stress resistance the metabolome predicts phenotypes more accurately than genomic data, and that the metabolites driving the increased accuracy was trait specific. Our comprehensive characterization of population-scale diversity of metabolomes and its genetic basis illustrates that metabolites have large potential as predictors of organismal phenotypes. This finding has strong applied importance e.g. in human medicine and animal and plant breeding.1 is a design matrix linking the genomic ( 1 ) and metabolomic ( 1 ) effects to the phenotypes. The random genomic effects are 1~( , [1,1] ! " ), and the metabolomic effects are 1~( , [1,1] 5 " ), where is the additive genomic relationship matrix as specified previously, and is the metabolomic relationship matrix.The metabolomic relationship matrix was computed as = ′ NMR ⁄ , where is a × NMR matrix of *77 18* T 6 , N 6 ), and was compared (using paired t-test corrected for multiple testing by a false discovery rate (FDR) of <0.05) within and across 9/ clusters to identify the NMR features resulting in the largest prediction accuracy. We only considered Drosophila melanogaster Genetic Reference Panel lines.