To date, variation in nectar chemistry of flowering plants has not been studied in detail. Such variation exerts considerable influence on pollinator-plant interactions, as well as on flower traits that play important roles in the selection of a plant for visitation by specific pollinators. Over the past 60 years the Aquilegia genus has been used as a key model for speciation studies. In this study, we defined the metabolomic profiles of flower samples of two Aquilegia species, A. Canadensis and A. pubescens. We identified a total of 75 metabolites that were classified into six main categories: organic acids, fatty acids, amino acids, esters, sugars, and unknowns. The mean abundances of 25 of these metabolites were significantly different between the two species, providing insights into interspecies variation in floral chemistry. Using the PlantSEED biochemistry database, we found that the majority of these metabolites are involved in biosynthetic pathways. Finally, we explored the annotated genome of A. coerulea, using the PlantSEED pipeline and reconstructed the metabolic network of Aquilegia. This network, which contains the metabolic pathways involved in generating the observed chemical variation, is now publicly available from the DOE Systems Biology Knowledge Base (KBase; http://kbase.us).