Abstract
Semi-polar metabolites such as flavonoids, phenolic acids, and alkaloids are very important health-related compounds in tomato. As a first step to identify genes responsible for the synthesis of semi-polar metabolites, quantitative trait loci (QTLs) that influence the semi-polar metabolite content in red-ripe tomato fruit were identified, by characterizing fruits of a population of introgression lines (ILs) derived from a cross between the cultivated tomato Solanum lycopersicum and the wild species Solanum chmielewskii. By analyzing fruits of plants grown at two different locations, we were able to identify robust metabolite QTLs for changes in phenylpropanoid glycoconjugation on chromosome 9, for accumulation of flavonol glycosides on chromosome 5, and for alkaloids on chromosome 7. To further characterize the QTLs we used a combination of genome sequencing, transcriptomics and targeted metabolomics to identify candidate key genes underlying the observed metabolic variation.
Introduction
Tomato (Solanum lycopersicum) is one of the most important vegetable crops worldwide with more than 160 million tons produced in 2013 ()1. This crop has served as a model organism for fleshy fruit plants and the complete genome sequence of one reference genome and up to 500 re-sequenced accessions is now available through the Solanaceae Genome Network (SGN)2 (). As with many other crop plants, tomato has been subjected to intensive domestication and breeding activities, which reduced the genetic variability in commercial materials. Domestication has been focused on yield, disease resistance, color and shape, while taste and nutritional value have long been neglected (). Currently, there is a growing demand to introduce novel genetic variation in commercial tomato in order to improve quality traits such as flavor and nutritional value. This genetic variation can be found in mutagenized populations, in core collections and in wild species and in introgression lines (ILs) derived from those. The potential of wild species as sources for genetic improvement of crops is increasingly recognized. A major goal of modern tomato breeding is to screen crossable wild Solanum species, such as Solanum lycopersicoides, Solanum pennelli, Solanum habrochaites, Solanum chmielewskii, Solanum pimpinellifolium, Solanum neorickii, Solanum peruvianum, and Solanum cheesmanii for valuable traits, such as resistance against various biotic and abiotic stresses (; ), primary metabolites () and secondary metabolites (). Wild species have been used as a source to develop ILs in S. lycopersicum, resulting in a set of lines each carrying a single or a few well-defined chromosome segments from the exotic germplasm source. These populations can be used to identify quantitative trait loci (QTLs) that improve crop quality once introgressed into an elite genetic background (). In addition to desirable traits, wild species also carry many agriculturally undesirable traits. Molecular genetic studies can identify the genetic and physical position of the underlying QTLs and introgression breeding can transfer the desirable traits into commercial varieties, while selecting against the undesirable ones.
The quality of tomato, in terms of nutritional value, taste, fragrance and appearance is essentially determined by its biochemical composition. To improve the quality of the crop, currently much research is devoted to the elucidation of the pathways and mechanisms that lead to the synthesis and accumulation of quality-related metabolites. The identification of QTLs that influence the chemical composition of ripe fruit, by screening IL populations, is an effective first step toward the identification of the underlying key genes that influence the nutritional quality of tomatoes. One of the best examples of this approach is the use of the founder tomato IL population, derived from a cross between the cultivated S. lycopersicum cv M82 and the green fruited wild species S. pennellii LA0716 (). This population has been used to identify QTLs for primary metabolites, volatile compounds, as well as semi-polar secondary metabolites, such as flavonoids and alkaloids (; ; ; ; ; ). These analyses also led to the identification of candidate genes involved in specific QTLs (; ; ; ), some of which have been shown to be the key gene underlying a specific QTL by reverse genetics studies (), while for others this still remains to be demonstrated. Despite the large number of studies related to primary metabolites and yield-associated traits, far less is known about QTLs determining secondary metabolites, such as flavonoids and alkaloids.
Flavonoids represent a large family of low molecular weight polyphenolic secondary metabolites that are widespread over the plant kingdom. To date, more than 6000 different flavonoids have been described and the number is still growing (). Based on their aglycone structure they can be grouped into several classes, such as chalcones, flavanones, flavonols, anthocyanins, and others. Flavonoids are involved in a diverse range of biological processes, such as pigmentation to attract pollinators and seed dispersers, protection against damage from ultraviolet light and pathogen resistance. In addition, they are associated with human health-promoting properties (; ). In tomato fruits, accumulation of flavonoids is restricted to the peel (, ; ). The main flavonoids present in tomato fruit peel are the chalcone naringenin-chalcone and various sugar conjugates of the flavonols quercetin and kaempferol. The structural information available about flavonoids and other semi-polar metabolites present in tomato increased substantially in the past decade, thanks to advances made in metabolomics tools, such as liquid chromatography and mass spectrometry (, ; ; ). However, our understanding of the genetic network regulating the accumulation of these compounds in tomato fruit is still incomplete. As indicated above, QTL analyses in interspecific IL populations can be used as a tool to identify key genes of this network in two ways: (i) qualitative and quantitative variation within and between metabolites, established by metabolic profiling of the complete set of ILs, can be used to determine the functional nature of the underlying key genes and (ii) precise knowledge of map positions of introgressions and the tomato genome sequence could facilitate the molecular cloning of these candidate genes. Previously, we demonstrated the success of this approach, by using an IL population derived from a cross between the commercial tomato cultivar S. lycopersicum cv. Moneyberg and the wild species S. chmielewskii (accession LA1840) to unravel the molecular and biochemical basis underlying the y mutation in tomato, which leads to pink-colored tomato fruits ().
Alkaloids are generally considered as anti-nutritional factors in our diet. Their biological effects in humans range from highly toxic, such as α-solanine and α-chaconine in potato tubers, to bitter tasting, such as α-tomatine in tomato. Domestication and breeding efforts have focused on reducing the levels of these anti-nutrients, but the success has been limited and some of these substances still remain in our daily diet (, ). In recent years, significant progress has been made in the elucidation of the steroidal glycoalkaloid pathway in Solanaceous species (, ; ; , ; ). In fruit of the cultivated tomato, the bitter tasting α-tomatine is present at high levels in early developmental stages and its levels decrease upon ripening due to its conversion into the acetyl glucosylated forms lycoperoside G, F or esculeoside A, which are not bitter. Putative intermediates in this conversion are hydroxytomatine (also called lycoperoside H), lycoperoside A, B, or C and hydroxylycoperoside A, B, or C, resulting from subsequent hydroxylation, acetylation and a second hydroxylation reactions. Fruits of many wild tomato species accumulate mostly the early, bitter, type of alkaloids ().
In the current study, we used the S. chmielewskii IL population to identify genomic regions controlling the production of semi-polar secondary metabolites, such as alkaloids, flavonoids and other phenylpropanoids, in tomato fruit. By combining biochemical pathway knowledge and genomic information, several candidate genes were identified. Further analysis of a major QTL on chromosome 5 for flavonols revealed the flavonoid pathway gene chalcone isomerase 1 (CHI1) as the key gene underlying the variation in quercetin- and kaempferol glycosides.
Materials and Methods
Plant Material and Growth Conditions
The IL population is composed of 34 indeterminate lines containing single or multiple introgressions from the wild species Solanum chmielewskii (LA1840) in the background of the commercial tomato variety Solanum lycopersicum cv. Moneyberg. The 34 ILs were grown in two greenhouses located in Avignon (Southern France) and Wageningen (The Netherlands) during spring and summer of 2007. From them, 25 were grown in both locations, five only in Avignon and four only in Wageningen. The day/night temperature set points were 25/15°C and 21/19°C in Avignon and Wageningen, respectively. At least nine plants were grown per each IL and each biological replicate consisted of at least six ripe fruit obtained from three different plants. Whole fruit was sampled from the plants grown in Wageningen, while fruit pericarp from plants grown in Avignon. After harvesting and sampling, the fruit material was immediately frozen in liquid nitrogen, ground to a fine frozen powder using an analytical electric mill and stored at -80°C until used for further analyses.
Selected ILs were grown again during the spring and summer of 2008 in Wageningen (The Netherlands). Samples were harvested at four stages of ripening (mature green (G), breaker (B), turning (T), and red (R)), which were judged by the fruit appearance and firmness. Subsequently, the fruit peel was carefully separated from the rest of the fruit (the flesh tissue) using a scalpel. Both flesh and peel tissues were immediately frozen in liquid nitrogen and stored at -80°C until used. Each biological replicate consisted of at least six fruit of the same ripening stage obtained from two different plants.
IL Genotyping
Illumina® infinium bead array was used for a high resolution mapping of the IL population. This analysis was performed as described in according to the Illumina® Infinium® HD Assay protocol: (Illumina® Infinium® HD Assay Ultra Protocol Guide. California, USA: ©Illumina, Inc; 2009. pp. 1–224. Catalog #WG-901-4007). The complete information of the SNP markers used is available at http://www.plantbreeding.wur.nl/Publications/SNP/4072SNP-Sequences.xlsx.
In addition, a set of PCR-based markers consisting of 130 COSII markers () and three simple sequence repeats (), previously mapped in the tomato genome and covering all 12 tomato chromosomes, were used to genotype the S. chmielewskii IL population (; ). Sequences of the primers are available on the Solanaceae Genomics Network Web site.
For the COSII markers, amplicon size differences between the two parents were detected in 12% of the cases and were used to genotype the IL population directly; in the other cases, the amplicons were digested with different restriction enzymes (TaqI, HinfI, AluI, DraI, RsaI, and MseI) to identify polymorphisms. Where no polymorphisms were detected, single-band amplicons were purified and sequenced. Amplicon sequences were aligned and examined for polymorphisms using the program CAPSdesigner3. Thereafter, the IL population was genotyped via cleaved-amplified polymorphic sequence assays ().
Analysis of Semi-polar Metabolites by LC-PDA-QTOF-MS
Semi-polar metabolites were extracted according to . Briefly, 500 mg of fresh weight tissue were extracted with 1.5 mL pure methanol (final methanol concentration in the extract approximately 75%). The samples were sonicated for 15 min, filtered through 0.2 μm inorganic membrane filter and 5 μL were used for the analysis.
Liquid chromatography quadrupole time of flight-mass spectrometry analyses (LC-PDA-QTOF-MS) were carried out according to . The detected flavonoid compounds were identified using authentic standards and accurate mass liquid chromatography mass spectrometry analysis using public databases (; ).
RNA Isolation and qRT-PCR Gene Expression Analysis
Total RNA was isolated from 150 mg of tomato fruit tissue using 1.50 mL of Trizol reagent (Invitrogen) according to the manufacturer’s instructions. Before cDNA synthesis, total RNA was treated with DNase-I Amplification Grade (Invitrogen) and purified with an RNeasy Mini Kit (Qiagen). And aliquot of 1 μg of total RNA was used for cDNA synthesis using the iScript cDNA synthesis kit (Bio-Rad Laboratories) in a 20-μL final volume according to the manufacturer. Expression levels of each gene were measured in duplicate reactions, performed with the same cDNA pool, in the presence of fluorescent dye (iQ SYBR Green Supermix) using an iCycler iQ instrument (Bio-Rad Laboratories) with specific primer pairs (Supplementary Table 5) (). The constitutively expressed mRNA encoding ubiquitine was used as internal reference. Expression levels were determined relative to the internal reference and multiplied by a factor 10. Calculations of each sample were carried out according to the comparative Ct method.
Microarray Analysis
The transcript profiling analysis was done using whole fruit tissue. Three biological replicates – pools of at least six fruit per plant were analyzed. Total RNA was extracted as described for real-time quantitative PCR. The 100 ng of total RNA was used to synthesize cDNA using Ambion WT expression kit (Applied Biosystems/Life Technologies, Nieuwekerk a/d IJssel, The Netherlands), which was subsequently labeled with biotin using the Affymetrix GeneChip WT Terminal Labeling Kit (Affymetrix, Santa Clara, CA, USA) and hybridized to Affymetrix EUTOM3 tomato exon arrays (Affymetrix). The microarray signals were determined using MadMax microarray analysis software4. The raw data can be found in Supplementary Data Sheet 1. Further analysis was performed using Genemaths XT microarray data analysis software (Applied Maths)5. Prior to analysis, the data were normalized using 2log transformation and subsequently scaled by subtraction of the mean (for each compound over the samples).
Student’s t-test was performed in Genemath XT, the Pearson correlation coefficients were calculated using the corresponding function of Microsoft Office Excel 2010.
Cloning, Sequencing, and Mapping of CHI1 Gene
Full-length cDNA sequences were amplified from ripe whole fruit of cv. Moneyberg, IL5b and IL7d using the SMART RACE cDNA amplification kit (Clontech Laboratories). Genomic DNA was amplified from leaves of the parental lines S. lycopersicum cv. Moneyberg and S. chmielewskii LA1840 using the GenomeWalker kit (Clontech Laboratories). The amplified sequences of both full-length cDNAs and genomic DNAs were cloned into pGEM-T-Easy vector (Promega) and sequenced.
Results
Physical Mapping of the S. lycopersicum × S. chmielewskii IL Population
The objective of this study was to discover genetic and genomic regions of S. chmielewskii LA1840 that affect accumulation of secondary metabolites in fruits of the commercial tomato S. lycopersicum cv. Moneyberg, as a first step toward discovering genes acting in the related metabolic pathways. For this purpose we analyzed a population of 25 S. chmielewskii LA1840 ILs. For 20 of them the first linkage maps, based on COSII and SSR markers, were presented in and . According to these data, 14 ILs each carried a single wild chromosomal introgression. In our study a high resolution genome wide SNP array, consisting of 5,528 SNP markers () was used to determine physical boundaries of the ILs. 1,660 markers were found to be polymorphic between S. chmielewskii LA1840 and S. lycopersicum cv. Moneyberg. As a result, out of the 25 ILs analyzed, 15 ILs were found to carry a single, homozygous S. chmielewskii introgression and 10 ILs carried two or more introgressions in one or in multiple chromosomes (Figure 1; Supplementary Table 1). Four major heterozygous introgressions were found on chromosomes 1, 3, 7, and 12.
FIGURE 1
Metabolic Profiling of the S. chmielewskii ILs Using Liquid Chromatography Coupled to Mass Spectrometry (LC-MS)
The S. chmielewskii IL population was grown at two different locations, Wageningen (The Netherlands) and Avignon (France), and ripe fruits were harvested. Three biological replicates were created by pooling fruit material of three independent plants per replicate. Semi-polar metabolites of ripe fruits were profiled using LC-MS. A total of 126 compounds were putatively identified in tomato fruit based on public mass spectral databases (; ). The mass spectra and the retention times were compared with authentic chemical standards when available (Supplementary Table 2). Different biochemical families of secondary metabolites were identified, including alkaloids, flavonoids (flavanones, flavones, and flavonols) and other phenylpropanoids.
Analysis of variance (ANOVA) showed that the content of 56 compounds was significantly affected (p < 0.05) in fruits of the ILs compared to fruits of cv. Moneyberg in both growing locations (Figure 2; Supplementary Table 3). Glycosylated volatile organic compounds (VOCs), alkaloids, and flavonoids were the most representative – 18, 17, and 14 compounds, respectively. Introgression in chromosomes 9 (IL9d) and 7 (IL7d) appeared to have the largest effects on the accumulation of different types of glycosylated VOCs and alkaloids, respectively. Two introgressions, on chromosomes 4 (IL4d) and 5 (IL5b), had a major effect on the accumulation of tomato fruit flavonols.
FIGURE 2
Introgression 5b in Chromosome 5 Increases the Accumulation of Kaempferol and Quercetin Glycosides
The result from the IL screening showed that the largest quantitative changes in levels of flavonols in ripe fruits were due to the presence of the IL5b introgression (Figure 2; Supplementary Table 3). To further investigate the accumulation of kaempferol and quercetin glycosides in tomato fruit, metabolic profiling was performed in different ripening stages (mature green (G), breaker (B), turning (T), and ripe (R)) of IL5b fruits and the control cv. Moneyberg. Since flavonoids normally accumulate in the fruit peel only (; ; ), we decided to focus the LC-MS analyses on methanolic extracts of the peel of IL5b and cv. Moneyberg fruits to increase the sensitivity of the measurements. An increase in the flavonols quercetin- and kaempferol-3-O-rutinoside (denoted as Q3R and K3R, respectively), quercetin- and kaempferol-3-O-rutinoside-7-O-glucoside (Q3R7G and K3R7G) and quercetin/kaempferol-3-O-glucose (Q3G and K3G) was observed in the fruit peel of IL5b compared to cv. Moneyberg, although the extent of the differences appeared to be compound-dependent (Figure 3). In order to unravel the genetic factors associated with the different patterns of accumulation of quercetin and kaempferol glycosides caused by the introgression in IL5b, a transcriptomics analysis of the ripening fruits (G, B/T, and R) was performed using the EU-TOM3 Affimetrix microarray. A total of 511 genes predicted by the International Tomato Annotation Group (ITAG) () were located in the chromosomal region corresponding to the IL5b introgression. Expression levels of 17 genes in the IL5b introgression region were found to be upregulated threefold or higher in turning fruits of IL5b compared to turning fruits of cv. Moneyberg (Supplementary Table 4). Among them, one gene – CHALCONE ISOMERASE 1 (CHI1) (Solyc05g010320) is directly involved in the flavonoid biosynthesis pathway (). To corroborate the results from the microarray experiments, the expression of a set of known fruit-expressed biosynthetic genes involved in the phenylpropanoid/flavonoid pathway () (Supplementary Table 5) was analyzed also in ripening fruit peel samples of IL5b and cv. Moneyberg, using qRT-PCR (Figure 4). Furthermore, we tested the expression of three additional putative CHI genes, which we denoted as CHI2 (Solyc05g052240), CHI3 (Solyc02g067870, BQ505699), and CHI4 (Solyc05g010310). Most of the genes tested showed a similar ripening-correlated pattern of expression in both IL5b and cv. Moneyberg: expression levels of phenylalanine ammonia-lyase (PAL), coumaroyl-4-hydroxylase (C4H), 4-coumarate ligase (4CL), chalcone synthases (CHS1 and CHS2), chalcone isomerases 2 and 3 (CHI2 and CHI3), flavonoid-3-hydroxylase (F3H), flavonoid-3′-hydroxylase (F3′H), flavonol synthase (FLS), flavonoid-3-O-glucosyltransferase (3GT) and flavonoid 3-O-glucoside-rhamnosyltransferase (RT) increased during ripening, peaked at B/T stage and decreased in ripe fruits. In contrast, the ripening-regulated pattern of CHI1 expression in cv. Moneyberg was opposite to the other genes in the phenylpropanoid/flavonoid pathway. This gene showed a low expression at G stage, which decreased even more at the later stages of ripening. This confirmed earlier observations (; ) that low expression of CHI1 is a major bottleneck in the biosynthesis of flavonols in fruits of cultivated tomatoes, such as cv. Moneyberg. In line with the microarray results, CHI1 expression in IL5b was significantly increased compared to cv. Moneyberg at G, B and T stages, suggesting that CHI1 expression relieves the block of the pathway in fruits of IL5b, which makes it the primary candidate gene for the flavonoid QTL mapped on IL5b.
FIGURE 3
FIGURE 4
Although all other flavonoid pathway genes did not show such a dramatic difference in expression as CHI1, in general, they tended to be expressed at higher levels in breaker and/or turning fruits of IL5b compared to cv. Moneyberg. This suggests that differences in flavonoid content between IL5b and cv. Moneyberg might also be due to coordinate control of flavonoid gene expression during ripening. The MYB12 transcription factor has previously been shown to regulate flavonol biosynthesis in tomato fruit (; ). Another MYB family transcription factor (Solyc05g009720) was found among the genes up-regulated in fruits of IL5b and physically located in the introgression region – at 3.93 Mb on chromosome 5. However, no significant correlation of expression was observed between this MYB gene and the 14 biosynthetic genes involved in the phenylpropanoid/flavonoid pathway present on the microarray (Supplementary Table 6). Analysis of two near-isogenic tomato lines only differing for a S. chmielewskii introgression in chromosome 5 that starts downstream of the MYB gene, but covers the CHI1 gene, confirmed the high-flavonoid fruit phenotype caused by the presence of the S. chmielewskii introgression (results not shown). This supports our conclusion that CHI1 is the primary candidate gene underlying the flavonoid QTL on chromosome 5.
Full length cDNAs of CHI1 (Solyc05g010320) were isolated from ripe fruit of both cv. Moneyberg and IL5b using Rapid Amplification of cDNA Ends (RACE). The derived protein sequences differ at only 2 amino acid positions (N35S and D137N in cv. Moneyberg→IL5b, Figure 5A). We cannot exclude that these two amino acid differences affect the function of the protein, but consider it unlikely that they account for the changes seen in CHI1 gene expression.
FIGURE 5
A genome walking approach was used to analyze and compare the genomic structure of CHI1 in the parental lines cv. Moneyberg and S. chmielewskii. The genomic sequence between both cv. Moneyberg and S. chmielewskii showed the presence of three introns, with the highest sequence variation observed in the first intron (Figure 5B). Deletions of 16, 38, 14, and 17 bp and an insertion of 17 bp were observed in cv. Moneyberg compared to S. chmielewskii. Comparison of the promoter regions, using the TSSP/Prediction of PLANT Promoters tool at the RegSite Plant DB (Softberry Inc.), revealed an insertion in the promoter region of S. chmielewskii at – 670 bp of the transcription start site, with a size of at least 2,063 bp. The insertion sequence was searched for homology to transposon-like sequences in RepBase (). In this insertion several fragments were found with 68 to 91% similarity to (from 5′to 3′) (i) twice a Copia-38_ST Long Terminal Repeat fragment from potato, (ii) a 34 nt sequence with high similarity (91%) to the polypurine tract containing region of ToRTL1, (iii) a Ty1/Copia long terminal repeat (LTR) retro-element (), and (iv) to 3′and 5′ sequences, respectively, of a hAT-like DNA transposon from potato ().
IL7d Affects Accumulation of Alkaloids in Tomato Fruit
According to the marker data IL7d carried a S. chmielewskii introgression in the top of chromosome 7 (0–2.86 Mb) (Figure 1; Supplementary Table 1). This introgression affected accumulation of two groups of alkaloids in fruits of this IL (Figure 2; Supplementary Table 3). α-tomatine, hydroxytomatine (lycoperoside H) and lycoperoside A/B/C had higher levels in fruits of IL7d compared to the control cv. Moneyberg (Figure 6), whereas the amounts of esculeoside A and lycoperoside F/G in fruits of this introgression line were reduced by up to 45-fold (Figure 6). According to the proposed tomato alkaloid biosynthetic pathway (; ) α-tomatine undergoes a number of ripening-induced hydroxylation and glycosylation modifications to produce the esculeoside type glycoalkaloids. Therefore, the accumulation of the green fruit-type alkaloids in fruits of IL7d suggests that this genomic region harbors a genetic factor which prevents or blocks the ripening-dependent glycoalkaloid modification (Figure 6). Accumulation of putative intermediates in the proposed pathway, such as hydroxy-lycoperoside A, B, or C, which after glycosylation produce the esculeoside type alkaloids, was not observed in fruits of IL7d. This suggests that the pathway is most likely interrupted at the step of hydroxylation of acetoxytomatine. Hydroxylation of alkaloids and other secondary metabolites is often mediated by enzymes of the cytochrome P450 family.
FIGURE 6
Three P450s were found to be located in the IL7d introgression region: Solyc07g006140, Solyc07g006890, and Solyc07g007460. Of these three genes, Solyc07g006890 was the most highly expressed in fruits and its transcript level increased during ripening. In ripe fruits of IL7d this gene showed a moderate threefold decrease in expression compared to its average expression observed in fruits of cv. Moneyberg and of IL5b (Supplementary Table 7).
IL9d Affects the Accumulation of Volatile Glycosides
The introgression 9d in chromosome 9 affected the accumulation of different glycoconjugate forms of volatile compounds, such as guaiacol, methyl salicylate (MeSA), and eugenol (Figure 2; Supplementary Table 3). IL9d led to the conversion of xylosyl-glucopyranoside forms of these volatiles into the corresponding xylosyl-diglucopyranoside forms. This conversion was shown to be mediated by the Non-Smoky Glycosyltransferase 1 gene (NSGT1), located within the IL9d introgression (Figure 1) () and suggests that S. chmielewskii carries a functional version of NSGT1. Indeed, re-sequencing of the S. chmielewskii LA1840 genome using a next generation sequencing approach revealed a gene with 96% homology to NSGT1 (Supplementary Image 1). In addition to glycosides of guaiacol, MeSA and eugenol, our data revealed that diglycosides of the aroma volatiles benzyl alcohol, 2-phenylethanol, and 2- or 3-methylbutanol were modified in the same manner (Figure 2; Supplementary Table 3).
Discussion
Many metabolite QTLs (mQTLs) have been described in ILs derived from the wild tomato relatives S. pennellii or S. habrochaites (; ; ; ; ; ). S. chmielewskii is another wild species crossable with cultivated tomato, with smaller, green fruits. The influence of the fruit load on the accumulation of dry matter and sugars in tomato fruit and on primary metabolites have been described recently in tomatoes derived from the S. chmielewskii IL population (; ). In this paper, we used a combination of genomic and metabolomics approaches to identify mQTLs and candidate genes controlling the synthesis of semi-polar compounds in tomato fruit. For this, ripe fruits harvested from the S. chmielewskii IL population were analyzed for variation in semi-polar secondary metabolites, using LC-PDA-QTOF-MS. By growing the plants at two different locations we could select for robust mQTLs. The screening revealed quantitative and qualitative changes in metabolites accumulating in specific ILs. The major mQTLs were found in ILs 4d, 5b, 7d, and 9d (Figure 2).
Accumulation of Specific Flavonol Glycosides in Tomato Fruit Related to an Increase in CHI Gene Expression
Our metabolic analyses revealed a major ripening-dependent increase of several flavonol glycosides in peel of IL5b compared to cv. Moneyberg (Figure 3). Absolute quantification of the main flavonols revealed an increase in quercetin-3-O-rutinoside from 30 to 260 mg/kg FW and of kaempferol-3-O-rutinoside from 3 to 35 mg/kg FW in peel of ripe IL5b compared to cv. Moneyberg fruits (Supplementary Image 2). Compared to the levels of these compounds found among a collection of 94 cultivated tomato hybrids (), the IL9d introgression upgrades Moneyberg tomatoes from a low-flavonol to a high-flavonol round tomato type, with levels comparable to those in cherry tomatoes, which are generally regarded as a much better source of flavonols than round/beef tomatoes.
Chalcone synthase (CHS) is the first enzyme involved in the phenylpropanoid/flavonoid pathway leading to the formation of these semi-polar compounds, most of which are present in a glycosylated form. Most of the biosynthetic genes involved in the flavonoid pathway and also transcription factors involved in the regulation of the biosynthetic genes have been identified (; ; ) and, due to the availability of the tomato genome sequence (), their physical position and chromosomal location is known (Figure 1; Supplementary Table 8). However, due to the complexity of flavonoid modification and the presence of more than 500 different forms of flavonoids in tomato (; ; ), further analysis is still needed to understand the flavonoid pathway to its full extent.
Most plants do not accumulate chalcones, the first class of flavonoids at the top of the biosynthetic pathway. After its formation, naringenin chalcone is usually rapidly isomerized by chalcone isomerase (CHI) to form the flavanone naringenin, a process that may also occur spontaneously in the absence of active CHI. However, in tomato fruit, low expression of CHI is rate-limiting and naringenin chalcone is the predominant yellow pigment that accumulates in the peel (; ). Four different putative CHI genes have been annotated in tomato, which share at most 75% identity at the amino acid level. The first one (Solyc05g010320), CHI1, has been described by and its expression was low in ripe tomato fruit, explaining the accumulation of the CHI substrate naringenin chalcone. CHI2 (Solyc05g052240) expression and CHI activity was increased in Del/Ros1 transgenic plants accumulating anthocyanins (GenBank acc. no. ES893795) (), and the expression of the third one, here called CHI3 (Solyc02g067870, BQ505699), was up-regulated in transgenic plants overexpressing the MYB transcription factor ANT1 leading to anthocyanin pigmentation in the fruit (). The expression of the fourth one (Solyc05g010310, here called CHI4), is very low compared to the other CHI genes (Figure 4). Based on the genome annotation, CHI1, CHI2 and CHI4 are located on chromosome 5, while CHI3 is located on chromosome 2 (Figure 1).
There are two major arguments supporting the conclusion that CHI1 is the key gene underlying the flavonoid QTL on chromosome 5. Firstly, after analyzing the expression of the biosynthetic flavonoid genes in tomato, CHI1 showed the highest expression increase in IL5b compared to cv. Moneyberg and was the only gene whose expression was significantly increased in IL5b compared to cv. Moneyberg at all stages of fruit ripening (Figure 4; Supplementary Table 4). The increase of the expression of this gene might redirect the flux of the pathway toward the formation of flavonol glycosides, as shown in the results of semi-polar metabolites detected in the fruits of IL5b and in line with results found in transgenic plants overexpressing the petunia CHI1 gene (). Secondly, CHI1 is located within the IL5b introgression and is among the 17 genes (out of 511) mapping in this region with an expression level at least threefold higher in turning IL5b fruits relative to cv. Moneyberg fruits (Supplementary Table 4).
IL5b is not a pure line in the sense that a small (<0.5 Mb) additional introgression region in chromosome 4 was also detected in this IL (Figure 1). However, within the subset of the S. chmielewskii population analyzed in this study, there are several other ILs with introgressions overlapping with the region in this chromosome. None of these lines showed an increase of flavonol glycosides compared to cv. Moneyberg and therefore we consider it unlikely that genes of this introgressed fragment might be responsible for the IL5b flavonoid QTL. In addition, NILs only differing in a chromosome 5 introgression showed a contrasting flavonoid accumulation pattern (results not shown), supporting that the flavonoid QTL is indeed due to the chromosome 5 introgression.
The lack of CHI1 gene expression in cultivated tomato might be due to (i) a mutation in a promoter regulatory sequence (cis-effect) and/or (ii) a mutation in a transcription factor responsible for expression of CHI () (trans-effect). We cannot completely exclude the latter possibility and, in this respect, we found a possible candidate MYB TF gene in the IL5b QTL region, which was upregulated in T stage fruit of IL5b. However, the expression pattern of this MYB gene during ripening was not correlated with the expression of CHI1, nor with the expression of the other flavonoid genes tested, which argues against a causal role for this candidate gene. After analyzing the genomic structure of the CHI1 gene, our results revealed the presence of several repetitive sequences related to transposons and retrotransposons of tomato and potato. Many examples exist of (retro-)transposon insertions influencing expression of down-stream genes such as, in tomato, retrotransposon ToRTL1 driving high expression of the 3-hydroxy-3-methylglutaryl coenzyme a reductase gene 2 (HMG2) (). A possible influence of the upstream transposon-like sequences of CHI1 on its expression, however, remains to be demonstrated.
Accumulation of α-Tomatine and Lycoperosides in Tomato Fruit
In our study, fruits of IL7d showed an accumulation of α-tomatine, hydroxytomatine and lycoperoside A, B, or C, while levels of lycoperoside G, F and esculeoside A were low compared with cv. Moneyberg fruit. None of the above-mentioned intermediates could be detected in fruits of IL5b or cv. Moneyberg.
A cluster of alkaloid biosynthesis genes has been previously discovered at the bottom of tomato and potato chromosome 7 (). Our genetic and metabolic data showed that there might be another genetic factor(s) at the top of this chromosome. The accumulation of the compounds from the first steps of the putative alkaloid pathway could be due to (i) a mutation of a gene involved in the hydroxylation of lycoperoside A, B, or C, or (ii) a mutation in regulatory element, such as a transcription factor (TF), responsible for the expression of the hydroxylation gene(s). An example of TF-mediated regulation of alkaloid biosynthesis has been recently shown by . Cytochrome P450s can catalyze aromatic hydroxylations, aliphatic hydroxylations and skeleton formation in secondary metabolite pathways in plants (), and therefore CYP’s would be good candidates for further studies. The evaluation of the genomic region of the IL7d introgression region revealed the presence of three cytochrome P450 genes: Solyc07g006140 (SL2.40ch07: 984510-988395 bp), Solyc07g006890 (SL2.40ch07: 1747021-1748529 bp), and Solyc07g007460 (SL2.40ch07: 2165949-2167526 bp). Further functional analysis of these candidate genes is currently underway.
Glycosylation of Tomato Fruit Volatiles
Our results indicate that the IL9d introgression carries a functional version of the NSGT1 gene, which mediates the conversion of xylosyl-glucopyranosides of the phenylpropanoid volatiles guaiacol, methyl salicylate and eugenol into the corresponding xylosyl-diglucopyranosides. This conversion affected the release of the corresponding volatiles and subsequently the fruit aroma (, ). In addition to glycosides of guaiacol, MeSA and eugenol, the present data showed that glycosides of other volatile compounds which play a role in tomato fruit aroma, such as benzyl alcohol, 2-phenylethanol, and 2- or 3-methylbutanol were modified in the same manner. This suggests that these diglycosides may be used as a substrate of the NSGT1 enzyme as well. This hypothesis could indeed be confirmed by metabolic analysis of transgenic NSGT1fruits () (Supplementary Image 3). In contrast to phenylpropanoid volatiles, the changes in glycosylation pattern of the other volatiles did not affect their release, neither in IL9d fruits (Supplementary Image 4), nor in transgenic NSGT1 fruits (). This suggests that, in addition to NSGT1-mediated glycosylation of the third sugar, the identity of the first two sugar conjugates and the interaction with putative glycoside hydrolases are important determinants for the release of aroma volatiles. We are currently aiming to get a better understanding of the various aspects of volatile “logistics” in tomato fruit and explore the opportunities to influence tomato fruit aroma by manipulating the volatile glycosylation status and thereby the release of these volatiles.
Conclusion
We identified a number of mQTLs involved in the production of semi-polar metabolites, by examining an IL population derived from a cross between S. lycopersicum cv. Moneyberg and S. chmielewskii LA1840. The use of specific S. chmielewskii ILs in combination with the knowledge gained on these mQTLs and the underlying candidate genes can be used to breed for tomatoes with improved quality. Reverse genetics is required to further elucidate the function of specific candidate genes, in order to gain a better understanding of the biosynthetic pathways leading to the synthesis and accumulation of health-related compounds.
Statements
Author contributions
A-RB carried out the research and wrote the manuscript. YT was responsible for the metabolomics analyses. JM was involved in the research, SG was involved in writing and correcting the manuscript, MV-Z and SH were involved in the marker analyses, RM was involved in manuscript preparation and data analysis, AB was involved in supervising the project, data analysis, and writing.
Funding
This work was carried out with support of the FP7 project EU-SOL (FOOD-CT-2006-016214), the Dutch genomics initiative Centre for Biosystems Genomics (CBSG) and the COST action Quality fruit (FA1106).
Acknowledgments
We thank Keygene Netherlands for kindly providing the S. chmielewskii IL population as well as Mr. Paul Dijkhuis, and Mrs. Fien Meijer-Dekens for excellent greenhouse management and plant cultivation. We thank Bert Schipper for assistance with the LC-PDA-QTOF-MS.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fpls.2016.01428
Footnotes
3.^https://solgenomics.net/tools/caps_designer/caps_input.pl
References
1
AdatoA.MandelT.Mintz-OronS.VengerI.LevyD.YativM.et al (2009). Fruit-surface flavonoid accumulation in tomato is controlled by a SlMYB12-regulated transcriptional network.PLoS Genet.5:e1000777. 10.1371/journal.pgen.1000777
2
AlseekhS.TohgeT.WendenbergR.ScossaF.OmranianN.LiJ.et al (2015). Identification and mode of inheritance of quantitative trait loci for secondary metabolite abundance in tomato.Plant Cell27485–512. 10.1105/tpc.114.132266
3
AyabeS-iAkashiT. (2006). Cytochrome P450s in flavonoid metabolism.Phytochem. Rev.5271–282. 10.1007/s11101-006-9007-3
4
BallesterA. R.MolthoffJ.de VosR.HekkertB. T.OrzaezD.Fernandez-MorenoJ.-P.et al (2010). Biochemical and molecular analysis of pink tomatoes: deregulated expression of the gene encoding transcription factor SlMYB12 leads to pink tomato fruit color.Plant Physiol.15271–84. 10.1104/pp.109.147322
5
BermudezL.UriasU.MilsteinD.KamenetzkyL.AsisR.FernieA. R.et al (2008). A candidate gene survey of quantitative trait loci affecting chemical composition in tomato fruit.J. Exp. Bot.592875–2890. 10.1093/jxb/ern146
6
BinoR. J.De VosC. H. R.LiebermanM.HallR. D.BovyA.JonkerH. H.et al (2005). The light-hyperresponsive high pigment-2dg mutation of tomato: alterations in the fruit metabolome.New Phytol.166427–438. 10.1111/j.1469-8137.2005.01362.x
7
BovyA.De VosR.KemperM.SchijlenE.Almenar PertejoM.MuirS.et al (2002). High-flavonol tomatoes resulting from the heterologous expression of the maize transcription factor genes LC and C1.Plant Cell142509–2526. 10.1105/tpc.004218
8
BovyA. G.Gomez-RoldanV.HallR. D. (2010). “Strategies to optimize the flavonoid content of tomato fruit,” in Recent Advances in Polyphenol Research, Vol II, edsSantosC.BuelgaM. T.Escribano-BailonLattanzioV. (Oxford: Wiley-Blackwell Publishing).
9
ButelliE.TittaL.GiorgioM.MockH.-P.MatrosA.PeterekS.et al (2008). Enrichment of tomato fruit with health-promoting anthocyanins by expression of select transcription factors.Nat. Biotechnol.261301–1308. 10.1038/nbt.1506
10
CárdenasP. D.SonawaneP. D.PollierJ.Vanden BosscheR.DewanganV.WeithornE.et al (2016). GAME9 regulates the biosynthesis of steroidal alkaloids and upstream isoprenoids in the plant mevalonate pathway.Na.t Commun.7:10654. 10.1038/ncomms10654
11
CausseM.DuffeP.GomezM. C.BuretM.DamidauxR.ZamirD.et al (2004). A genetic map of candidate genes and QTLs involved in tomato fruit size and composition.J. Exp. Bot.551671–1685. 10.1093/jxb/erh207
12
ColliverS.BovyA.CollinsG.MuirS.RobinsonS.De VosC. H. R.et al (2002). Improving the nutritional content of tomatoes through reprogramming their flavonoid biosynthetic pathway.Phytochem. Rev.1113–123. 10.1023/A:1015848724102
13
DaraseliaN. D.TarchevskayaS.NaritaJ. O. (1996). The promoter for tomato 3-hydroxy-3-methylglutaryl coenzyme A reductase gene 2 has unusual regulatory elements that direct high-level expression.Plant Physiol.112727–733. 10.1104/pp.112.2.727
14
DoP. T.PrudentM.SulpiceR.CausseM.FernieA. R. (2010). The influence of fruit load on the tomato pericarp metabolome in a Solanum chmielewskii introgression line population.Plant Physiol.1541128–1142. 10.1104/pp.110.163030
15
EshedY.ZamirD. (1995). An introgression line population of Lycopersicon pennellii in the cultivated tomato enables the identification and fine mapping of yield- associated QTL.Genetics1411147–1162.
16
FAOSTAT (2015). Food and Agriculture Organization of the United Nations, FAOSTAT Database. Available at: http://faostat3.fao.org/
17
FrankelN.HassonE.IusemN. D.RossiM. S. (2003). Adaptive evolution of the water stress-induced gene Asr2 in Lycopersicon species dwelling in arid habitats.Mol. Biol. Evol.201955–1962. 10.1093/molbev/msg214
18
FraryA.XuY.LiuJ.MitchellS.TedeschiE.TanksleyS. (2005). Development of a set of PCR-based anchor markers encompassing the tomato genome and evaluation of their usefulness for genetics and breeding experiments.Theor. Appl. Genet.111291–312. 10.1007/s00122-005-2023-7
19
FridmanE.LiuY. S.Carmel-GorenL.GurA.ShoreshM.PlebanT.et al (2001). Two tightly linked QTLs modify tomato sugar content via different physiological pathways.Mol. Genet. Genomics266821–826.
20
FriedmanM. (2002). Tomato glycoalkaloids: role in the plant and in the diet.J. Agric. Food Chem.505751–5780. 10.1021/jf020560c
21
FriedmanM. (2006). Potato glycoalkaloids and metabolites: roles in the plant and in the diet.J. Agric. Food Chem.548655–8681. 10.1021/jf061471t
22
GrennanA. K. (2009). MoTo DB: a metabolic database for tomato.Plant Physiol.1511701–1702. 10.1104/pp.109.900308
23
HarborneJ. B.WilliamsC. A. (2000). Advances in flavonoid research since 1992.Phytochemistry55481–504. 10.1016/S0031-9422(00)00235-1
24
IijimaY.NakamuraY.OgataY.TanakaK.SakuraiN.SudaK.et al (2008). Metabolite annotations based on the integration of mass spectral information.Plant J.54949–962. 10.1111/j.1365-313X.2008.03434.x
25
IijimaY.WatanabeB.SasakiR.TakenakaM.OnoH.SakuraiN.et al (2013). Steroidal glycoalkaloid profiling and structures of glycoalkaloids in wild tomato fruit.Phytochemistry95145–157. 10.1016/j.phytochem.2013.07.016
26
ItkinM.HeinigU.TzfadiaO.BhideA. J.ShindeB.CardenasP. D.et al (2013). Biosynthesis of antinutritional alkaloids in solanaceous crops is mediated by clustered genes.Science341175–179. 10.1126/science.1240230
27
ItkinM.RogachevI.AlkanN.RosenbergT.MalitskyS.MasiniL.et al (2011). GLYCOALKALOID METABOLISM1 is required for steroidal alkaloid glycosylation and prevention of phytotoxicity in tomato.Plant Cell234507–4525. 10.1105/tpc.111.088732
28
JurkaJ.KapitonovV. V.PavlicekA.KlonowskiP.KohanyO.WalichiewiczJ. (2005). Repbase update, a database of eukaryotic repetitive elements.Cytogenet. Genome Res.110462–467. 10.1159/000084979
29
JurkaJ.KohanyO. (2006). hAT-2_SD hAT transposon from Solanum demissum.Repbase Rep.6492–492.
30
KoesR. E.QuattrocchioF.MolJ. N. M. (1994). The flavonoid biosynthetic pathway in plants: function and evolution.Bioessays16123–132. 10.1002/bies.950160209
31
KoniecznyA.AusubelF. M. (1993). A procedure for mapping Arabidopsis mutations using co-dominant ecotype-specific PCR-based markers.Plant J.4403–410. 10.1046/j.1365-313X.1993.04020403.x
32
LégnaniR.GognalonsP.Gébré SélassiéK.MarchouxG.MorettiA.LaterrotH. (1996). Identification and characterization of resistance to tobacco etch virus in Lycopersicon species.Plant Dis.80306–309. 10.1094/PD-80-0306
33
LinT.ZhuG.ZhangJ.XuX.YuQ.ZhengZ.et al (2014). Genomic analyses provide insights into the history of tomato breeding.Nat. Genet.461220–1226. 10.1038/ng.3117
34
MathewsH.ClendennenS. K.CaldwellC. G.LiuX. L.ConnorsK.MatheisN.et al (2003). Activation tagging in tomato identifies a transcriptional regulator of anthocyanin biosynthesis, modification, and transport.Plant Cell151689–1703. 10.1105/tpc.012963
35
MathieuS.CinV. D.FeiZ.LiH.BlissP.TaylorG.et al (2009). Flavour compounds in tomato fruits: identification of loci and potential pathways affecting volatile composition.J. Exp. Bot.60325–337. 10.1093/jxb/ern294
36
Mintz-OronS.MandelT.RogachevI.FeldbergL.LotanO.YativM.et al (2008). Gene expression and metabolism in tomato fruit surface tissues.Plant Physiol.147823–851. 10.1104/pp.108.116004
37
MocoS.BinoR. J.VorstO.VerhoevenH. A.de GrootJ.van BeekT. A.et al (2006). A liquid chromatography-mass spectrometry-based metabolome database for tomato.Plant Physiol.1411205–1218. 10.1104/pp.106.078428
38
MocoS.CapanogluE.TikunovY.BinoR. J.BoyaciogluD.HallR. D.et al (2007). Tissue specialization at the metabolite level is perceived during the development of tomato fruit.J. Exp. Bot.584131–4146. 10.1093/jxb/erm271
39
MuirS. R.CollinsG. J.RobinsonS.HughesS.BovyA.Ric de VosC. H.et al (2001). Overexpression of petunia chalcone isomerase in tomato results in fruit containing increased levels of flavonols.Nat. Biotechnol.19470–474. 10.1038/88150
40
PrudentM.CausseM.GénardM.TripodiP.GrandilloS.BertinN. (2009). Genetic and physiological analysis of tomato fruit weight and composition: Influence of carbon availability on QTL detection.J. Exp. Bot.60923–937. 10.1093/jxb/ern338
41
SchauerN.SemelY.RoessnerU.GurA.BalboI.CarrariF.et al (2006). Comprehensive metabolic profiling and phenotyping of interspecific introgression lines for tomato improvement.Nat. Biotechnol.24447–454. 10.1038/nbt1192
42
SchauerN.ZamirD.FernieA. R. (2005). Metabolic profiling of leaves and fruit of wild species tomato: a survey of the Solanum lycopersicum complex.J. Exp. Bot.56297–307. 10.1093/jxb/eri057
43
SchijlenE. G. W. M.BeekwilderJ.HallR. D.MeerI. M. (2008). Boosting beneficial phytochemicals in vegetable crop plants.CAB Rev.31–21. 10.1079/PAVSNNR20083025
44
SchijlenE. G. W. M.De VosC. H. R.Van TunenA. J.BovyA. G. (2004). Modification of flavonoid biosynthesis in crop plants.Phytochemistry652631–2648. 10.1016/j.phytochem.2004.07.028
45
SemelY.NissenbaumJ.MendaN.ZinderM.KriegerU.IssmanN.et al (2006). Overdominant quantitative trait loci for yield and fitness in tomato.Proc. Natl. Acad. Sci. U.S.A.10312981–12986. 10.1073/pnas.0604635103
46
SteinhauserM.-C.SteinhauserD.GibonY.BolgerM.ArrivaultS.UsadelB.et al (2011). Identification of enzyme activity quantitative trait loci in a Solanum lycopersicum x Solanum pennellii introgression line population.Plant Physiol.157998–1014. 10.1104/pp.111.181594
47
StevensR.BuretM.DuffeP.GarcheryC.BaldetP.RothanC.et al (2007). Candidate genes and quantitative trait loci affecting fruit ascorbic acid content in three tomato populations.Plant Physiol.1431943–1953. 10.1104/pp.106.091413
48
TapasA. R.SakarkarD. M.KakdeR. B. (2008). Flavonoids as nutraceuticals: a review.Trop. J. Pharm. Res.71089–1099. 10.4314/tjpr.v7i3.14693
49
TiemanD. M.ZeiglerM.SchmelzE. A.TaylorG.BlissP.KirstM.et al (2006). Identification of loci affecting flavour volatile emissions in tomato fruits.J. Exp. Bot.57887–896. 10.1093/jxb/erj074
50
TikunovY. M.de VosR. C. H.Gonzalez ParamasA. M.HallR. D.BovyA. G. (2010). A role for differential glycoconjugation in the emission of phenylpropanoid volatiles from tomato fruit discovered using a metabolic data fusion approach.Plant Physiol.15255–70. 10.1104/pp.109.146670
51
TikunovY. M.MolthoffJ.de VosR. C. H.BeekwilderJ.van HouwelingenA.van der HooftJ. J. J.et al (2013). Non-smoky GLYCOSYLTRANSFERASE1 prevents the release of smoky aroma from tomato fruit.Plant Cell253067–3078. 10.1105/tpc.113.114231
52
Tomato Genome Consortium (2012). The tomato genome sequence provides insights into fleshy fruit evolution.Nature485635–641. 10.1038/nature11119
53
ToubianaD.SemelY.TohgeT.BeleggiaR.CattivelliL.RosentalL.et al (2012). Metabolic profiling of a mapping population exposes new insights in the regulation of seed metabolism and seed, fruit, and plant relations.PLoS Genet.8:e1002612. 10.1371/journal.pgen.1002612
54
Viquez-ZamoraM.VosmanB.van de GeestH.BovyA. G.VisserR. G. F.FinkersR.et al (2013). Tomato breeding in the genomic era: insights from a SNP array.BMC Genet.14:354. 10.1186/1471-2164-14-354
55
WillitsG.KramerC. M.PrataR. T. N.De LucaV.PotterB. G.StephensJ. C.et al (2005). Utilization of the genetic resources of wild species to create a nontransgenic high flavonoid tomato.J. Agric. Food Chem.531231–1236. 10.1021/jf049355i
56
WuF.MuellerL. A.CrouzillatD.PetiardV.TanksleyS. D. (2006). Combining bioinformatics and phylogenetics to identify large sets of single copy, orthologous genes (COSII) for comparative, evolutionary and systematic studies: a test case in the euasterid plant clade.Genetics1741407–1420. 10.1534/genetics.106.062455
57
ZamirD. (2001). Improving plant breeding with exotic genetic libraries.Nat. Rev. Genet.2983–989. 10.1038/ni1101-983
58
ZanorM. I.OsorioS.Nunes-NesiA.CarrariF.LohseM.UsadelB.et al (2009). RNA interference of LIN5 in tomato confirms its role in controlling brix content, uncovers the influence of sugars on the levels of fruit hormones, and demonstrates the importance of sucrose cleavage for normal fruit development and fertility.Plant Physiol.1501204–1218. 10.1104/pp.109.136598
Summary
Keywords
tomato (Solanum lycopersicum), QTL analysis, flavonoids, alkaloids, introgression lines
Citation
Ballester A-R, Tikunov Y, Molthoff J, Grandillo S, Viquez-Zamora M, de Vos R, de Maagd RA, van Heusden S and Bovy AG (2016) Identification of Loci Affecting Accumulation of Secondary Metabolites in Tomato Fruit of a Solanum lycopersicum × Solanum chmielewskii Introgression Line Population. Front. Plant Sci. 7:1428. doi: 10.3389/fpls.2016.01428
Received
22 April 2016
Accepted
07 September 2016
Published
28 September 2016
Volume
7 - 2016
Edited by
Ana Margarida Fortes, University of Lisbon, Portugal
Reviewed by
Roy Navarre, United Stated Department of Agriculture, USA; Juan Capel, University of Almería, Spain
Updates
Copyright
© 2016 Ballester, Tikunov, Molthoff, Grandillo, Viquez-Zamora, de Vos, de Maagd, van Heusden and Bovy.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Arnaud G. Bovy, arnaud.bovy@wur.nl
†Present address: Ana-Rosa Ballester, Instituto de Agroquímica y Tecnología de Alimentos, Consejo Superior de Investigaciones Científicas, Valencia, Spain
This article was submitted to Plant Physiology, a section of the journal Frontiers in Plant Science
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.