Solanum pennellii is a wild tomato species endemic to Andean regions in South America, where it has evolved to thrive in arid habitats. Because of its extreme stress tolerance and unusual morphology, it is an important donor of germplasm for the cultivated tomato Solanum lycopersicum 1 . Introgression lines (ILs) in which large genomic regions of S. lycopersicum are replaced with the corresponding segments from S. pennellii can show remarkably superior agronomic performance 2 . Here we describe a high-quality genome assembly of the parents of the IL population. By anchoring the S. pennellii genome to the genetic map, we define candidate genes for stress tolerance and provide evidence that transposable elements had a role in the evolution of these traits. Our work paves a path toward further tomato improvement and for deciphering the mechanisms underlying the myriad other agronomic traits that can be improved with S. pennellii germplasm.Crosses between distantly related plants can lead to substantial improvements in performance. Notably, S. pennellii × S. lycopersicum ILs have been used to define numerous quantitative trait loci (QTLs) for superior yield, chemical composition, morphology, abiotic stress tolerance and extreme heterosis 3,4 . Although genetic studies have proven informative, few genes underlying specific QTLs have been cloned, largely because of the lack of a S. pennellii genome sequence. To support QTL analyses, we sequenced the genome of S. pennellii using Illumina sequencing with ~190-fold coverage ( Fig. 1 and Supplementary Tables 1-5). The initial assembly size was 942 Mb, with a scaffold N50 value of 1.7 Mb and N90 value of 0.43 Mb (Table 1 and Supplementary Tables 6 and 7). We estimated the total genome size to be about 1.2 Gb using a k-mer-based analysis ( Supplementary Fig. 1 and Supplementary Table 8), in accordance with previous estimations 3,4 . We anchored 97.1% of the genome assembly to chromosomes using genetic maps and restriction site-associated DNA sequencing (RAD-seq)-based markers from the IL population 5 (Supplementary Note). Comparison of the assembly to publicly available BAC sequences indicated an accuracy of >99.9%, and a satisfactory accuracy of gap-filled regions was shown by realigning reads (Supplementary Fig. 2 and Supplementary Table 9). Of the 307,350 S. lycopersicum and 7,812 S. pennellii publicly available ESTs, 93% and >96% could be aligned to the genome, respectively (Supplementary Table 10), indicating comprehensive coverage of the gene-rich regions. We predicted 32,273 high-confidence genes and a potential set of 44,966 protein-coding genes and checked these

Metabolomics, the large-scale study of the metabolic complement of the cell [1][2][3] , is a mature science that has been practiced for over 20 years 4 . Indeed, it is now a commonly used experimental systems biology tool with demonstrated utility in both fundamental and applied aspects of plant, microbial and mammalian research [5][6][7][8][9][10][11][12][13][14][15] . Among the many thousands of studies published in this area over the last 20 years, notable highlights [5][6][7][8]10,11,16 are briefly described in Supplementary Note 1.Despite the insight afforded by such studies, the nature of metabolites, particularly their diversity (in both chemical structure and dynamic range of abundance 9,12 ), remains a major challenge with regard to the ability to provide adequate coverage of the metabolome that can complement that achieved for the genome, transcriptome and proteome. Despite these comparative limitations, enormous advances have been made with regard to the number of analytes about which accurate quantitative information can be acquired, and a vast number of studies have yielded important biological information and biologically active metabolites across the kingdoms of life 14 . We have previously estimated that upwards of 1 million different metabolites occur across the tree of life, with between 1,000 and 40,000 estimated to occur in a single species 4 .

Phenylpropanoids comprise an important class of plant secondary metabolites. A number of transcription factors have been used to upregulate-specific branches of phenylpropanoid metabolism, but by far the most effective has been the fruit-specific expression of AtMYB12 in tomato, which resulted in as much as 10% of fruit dry weight accumulating as flavonols and hydroxycinnamates. We show that AtMYB12 not only increases the demand of flavonoid biosynthesis but also increases the supply of carbon from primary metabolism, energy and reducing power, which may fuel the shikimate and phenylalanine biosynthetic pathways to supply more aromatic amino acids for secondary metabolism. AtMYB12 directly binds promoters of genes encoding enzymes of primary metabolism. The enhanced supply of precursors, energy and reducing power achieved by AtMYB12 expression can be harnessed to engineer high levels of novel phenylpropanoids in tomato fruit, offering an effective production system for bioactives and other high value ingredients.