Plant metabolites are important to world food security due to their roles in crop yield and nutritional quality. Here we report the metabolic profile of 300 tomato accessions (Solanum lycopersicum and related wild species) by quantifying 60 primary and secondary metabolites, including volatile organic compounds, over a period of 2 yr. Metabolite content and genetic inheritance of metabolites varied broadly, both within and between different genetic groups. Using genotype information gained from 10 000 single nucleotide polymorphism markers, we performed a metabolite genome-wide association mapping (GWAS) study. We identified 79 associations influencing 13 primary and 19 secondary metabolites with large effects at high resolution. Four genome regions were detected, highlighting clusters of associations controlling the variation of several metabolites. Local linkage disequilibrium analysis and allele mining identified possible candidate genes which may modulate the content of metabolites that are of significant importance for human diet and fruit consumption. We precisely characterized two associations involved in fruit acidity and phenylpropanoid volatile production. Taken together, this study reveals complex and distinct metabolite regulation in tomato subspecies and demonstrates that GWAS is a powerful tool for gene-metabolite annotation and identification, pathways elucidation, and further crop improvement.

HighlightTomato quality could be improved under deficit irrigation while maintaining yield. The underlying genetic architecture is polygenic and varies with water availability. Candidate genes related to primary metabolism were identified.

As a result of climate change, drought will increasingly limit crop production in the future. Studying genotype by watering regime interactions is necessary to improve plant adaptation to low water availability. In cultivated tomato (Solanum lycopersicum L.), extensively grown in dry areas, well-mastered water deficits can stimulate metabolite production, increasing plant defenses and concentration of compounds involved in fruit quality, at the same time. However, few tomato Quantitative Trait Loci (QTLs) and genes involved in response to drought are identified or only in wild species. In this study, we phenotyped a population of 119 recombinant inbred lines derived from a cross between a cherry tomato and a large fruit tomato, grown in greenhouse under two watering regimes, in two locations. A large genetic variability was measured for 19 plant and fruit traits, under the two watering treatments. Highly significant genotype by watering regime interactions were detected and resulted from re-ranking more than scale changes. The population was genotyped for 679 SNP markers to develop a genetic map. In total, 56 QTLs were identified among which 11 were interactive between watering regimes. These later mainly exhibited antagonist effects according to watering treatment. Variation in gene expression in leaves of parental accessions revealed 2259 differentially expressed genes, among which candidate genes presenting sequence polymorphisms were identified under two main interactive QTLs. Our results provide knowledge about the genetic control of genotype by watering regime interactions in cultivated tomato and the possible use of deficit irrigation to improve tomato quality.

Characterizing the natural diversity of gene expression across environments is an important step in understanding how genotype-by-environment interactions shape phenotypes. Here, we analyzed the impact of water deficit onto gene expression levels in tomato at the genome-wide scale. We sequenced the transcriptome of growing leaves and fruit pericarps at cell expansion stage in a cherry and a large fruited accession and their F hybrid grown under two watering regimes. Gene expression levels were steadily affected by the genotype and the watering regime. Whereas phenotypes showed mostly additive inheritance, ~80% of the genes displayed non-additive inheritance. By comparing allele-specific expression (ASE) in the F hybrid to the allelic expression in both parental lines, respectively, 3005 genes in leaf and 2857 genes in fruit deviated from 1:1 ratio independently of the watering regime. Among these genes, ~55% were controlled by cis factors, ~25% by trans factors and ~20% by a combination of both types of factors. A total of 328 genes in leaf and 113 in fruit exhibited significant ASE-by-watering regime interaction, among which ~80% presented trans-by-watering regime interaction, suggesting a response to water deficit mediated through a majority of trans-acting loci in tomato. We cross-validated the expression levels of 274 transcripts in fruit and leaves of 124 recombinant inbred lines (RILs) and identified 163 expression quantitative trait loci (eQTLs) mostly confirming the divergences identified by ASE. Combining phenotypic and expression data, we observed a complex network of variation between genes encoding enzymes involved in the sugar metabolism.

Water deficit (WD) leads to significant phenotypic changes in crops resulting from complex stress regulation mechanisms involving responses at the physiological, biochemical and molecular levels. Tomato growth and fruit quality have been shown to be significantly affected by WD stress. Understanding the molecular mechanism underlying response to WD is crucial to develop tomato cultivars with relatively high performance under low watering conditions. Transcriptome response to WD was investigated through the RNA sequencing of fruit and leaves in eight accessions grown under two irrigation conditions, in order to get insight into the complex genetic regulation of WD response in tomato. Significant differences in genotype WD response were first observed at the phenotypic level for fruit composition and plant development traits. At the transcriptome level, a total of 14,065 differentially expressed genes (DEGs) in response to WD were detected, among which 7393 (53%) and 11,059 (79%) were genotype- and organ-specific, respectively. Water deficit induced transcriptome variations much stronger in leaves than in fruit. A significant effect of the genetic background on expression variation was observed compared to the WD effect, along with the presence of a set of genes showing a significant genotype x watering regime interaction. Integrating the DEGs with previously identified WD response quantitative trait loci (QTLs) mapped in a multi-parental population derived from the crossing of the eight genotypes narrowed the candidate gene lists to within the confidence intervals surrounding the QTLs. The results present valuable resources for further study to decipher the genetic determinants of tomato response to WD.

No abstract