Analysis of wild tomato introgression lines elucidates the genetic basis of transcriptome and metabolome variation underlying fruit traits and pathogen response.

Wild tomato species represent a rich gene pool for numerous desirable traits lost during domestication. Here, we exploited an introgression population representing wild desert-adapted species and a domesticated cultivar to establish the genetic basis of gene expression and chemical variation accompanying the transfer of wild-species-associated fruit traits. Transcriptome and metabolome analysis of 580 lines coupled to pathogen sensitivity assays resulted in the identification of genomic loci associated with levels of hundreds of transcripts and metabolites. These associations occurred in hotspots representing coordinated perturbation of metabolic pathways and ripening-related processes. Here, we identify components of the Solanum alkaloid pathway, as well as genes and metabolites involved in pathogen defense and linking fungal resistance with changes in the fruit ripening regulatory network. Our results outline a framework for understanding metabolism and pathogen resistance during tomato fruit ripening and provide insights into key fruit quality traits.

Pathways to defense metabolites and evading fruit bitterness in genus Solanum evolved through 2-oxoglutarate-dependent dioxygenases.

The genus Solanum comprises three food crops (potato, tomato, and eggplant), which are consumed on daily basis worldwide and also producers of notorious anti-nutritional steroidal glycoalkaloids (SGAs). Hydroxylated SGAs (i.e. leptinines) serve as precursors for leptines that act as defenses against Colorado Potato Beetle (Leptinotarsa decemlineata Say), an important pest of potato worldwide. However, SGA hydroxylating enzymes remain unknown. Here, we discover that 2-OXOGLUTARATE-DEPENDENT-DIOXYGENASE (2-ODD) enzymes catalyze SGA-hydroxylation across various Solanum species. In contrast to cultivated potato, Solanum chacoense, a widespread wild potato species, has evolved a 2-ODD enzyme leading to the formation of leptinines. Furthermore, we find a related 2-ODD in tomato that catalyzes the hydroxylation of the bitter α-tomatine to hydroxytomatine, the first committed step in the chemical shift towards downstream ripening-associated non-bitter SGAs (e.g. esculeoside A). This 2-ODD enzyme prevents bitterness in ripe tomato fruit consumed today which otherwise would remain unpleasant in taste and more toxic.

A Solanum neorickii introgression population providing a powerful complement to the extensively characterized Solanum pennellii population.

We present a complementary resource for trait fine-mapping in tomato to those based on the intra-specific cross between cultivated tomato and the wild tomato species Solanum pennellii, which have been extensively used for quantitative genetics in tomato over the last 20 years. The current population of backcross inbred lines (BILs) is composed of 107 lines derived after three backcrosses of progeny of the wild species Solanum neorickii (LA2133) and cultivated tomato (cultivar TA209) and is freely available to the scientific community. These S. neorickii BILs were genotyped using the 10K SolCAP single nucleotide polymorphism chip, and 3111 polymorphic markers were used to map recombination break points relative to the physical map of Solanum lycopersicum. The BILs harbor on average 4.3 introgressions per line, with a mean introgression length of 34.7 Mbp, allowing partitioning of the genome into 340 bins and thereby facilitating rapid trait mapping. We demonstrate the power of using this resource in comparison with archival data from the S. pennellii resources by carrying out metabolic quantitative trait locus analysis following gas chromatography-mass spectrometry on fruits harvested from the S. neorickii BILs. The metabolic candidate genes phenylalanine ammonia-lyase and cystathionine gamma-lyase were then tested and validated in F2 populations and via agroinfiltration-based overexpression in order to exemplify the fidelity of this method in identifying the genes that drive tomato metabolic phenotypes.

Quantitative Trait Loci Analysis Identifies a Prominent Gene Involved in the Production of Fatty Acid-Derived Flavor Volatiles in Tomato.

To gain insight into the genetic regulation of lipid metabolism in tomato, we conducted metabolic trait loci (mQTL) analysis following the lipidomic profiling of fruit pericarp and leaf tissue of the Solanum pennellii introgression lines (IL). To enhance mapping resolution for selected fruit-specific mQTL, we profiled the lipids in a subset of independently derived S. pennellii backcross inbred lines, as well as in a near-isogenic sub-IL population. We identified a putative lecithin:cholesterol acyltransferase that controls the levels of several lipids, and two members of the class III lipase family, LIP1 and LIP2, that were associated with decreased levels of diacylglycerols (DAGs) and triacylglycerols (TAGs). Lipases of this class cleave fatty acids from the glycerol backbone of acylglycerols. The released fatty acids serve as precursors of flavor volatiles. We show that LIP1 expression correlates with fatty acid-derived volatile levels. We further confirm the function of LIP1 in TAG and DAG breakdown and volatile synthesis using transgenic plants. Taken together, our study extensively characterized the genetic architecture of lipophilic compounds in tomato and demonstrated at molecular level that release of free fatty acids from the glycerol backbone can have a major impact on downstream volatile synthesis.

Canalization of Tomato Fruit Metabolism.

To explore the genetic robustness (canalization) of metabolism, we examined the levels of fruit metabolites in multiple harvests of a tomato introgression line (IL) population. The IL partitions the whole genome of the wild species Solanum pennellii in the background of the cultivated tomato (Solanum lycopersicum). We identified several metabolite quantitative trait loci that reduce variability for both primary and secondary metabolites, which we named canalization metabolite quantitative trait loci (cmQTL). We validated nine cmQTL using an independent population of backcross inbred lines, derived from the same parents, which allows increased resolution in mapping the QTL previously identified in the ILs. These cmQTL showed little overlap with QTL for the metabolite levels themselves. Moreover, the intervals they mapped to harbored few metabolism-associated genes, suggesting that the canalization of metabolism is largely controlled by regulatory genes.

De Novo Assembly of a New Solanum pennellii Accession Using Nanopore Sequencing.

Updates in nanopore technology have made it possible to obtain gigabases of sequence data. Prior to this, nanopore sequencing technology was mainly used to analyze microbial samples. Here, we describe the generation of a comprehensive nanopore sequencing data set with a median read length of 11,979 bp for a self-compatible accession of the wild tomato species Solanum pennellii We describe the assembly of its genome to a contig N50 of 2.5 MB. The assembly pipeline comprised initial read correction with Canu and assembly with SMARTdenovo. The resulting raw nanopore-based de novo genome is structurally highly similar to that of the reference S. pennellii LA716 accession but has a high error rate and was rich in homopolymer deletions. After polishing the assembly with Illumina reads, we obtained an error rate of <0.02% when assessed versus the same Illumina data. We obtained a gene completeness of 96.53%, slightly surpassing that of the reference S. pennellii Taken together, our data indicate that such long read sequencing data can be used to affordably sequence and assemble gigabase-sized plant genomes.

Solanum pennellii backcross inbred lines (BILs) link small genomic bins with tomato traits.

We present a resource for fine mapping of traits derived from the wild tomato species Solanum pennellii (LA0716). The population of backcross inbred lines (BILs) is composed of 446 lines derived after a few generations of backcrosses of the wild species with cultivated tomato (cultivar M82; LA3475), followed by more than seven generations of self-pollination. The BILs were genotyped using the 10K SOL-CAP single nucleotide polymorphism (SNP) -Chip, and 3700 polymorphic markers were used to map recombination break points relative to the physical map of Solanum lycopersicum. The BILs carry, on average, 2.7 introgressions per line, with a mean introgression length of 11.7 Mbp. Whereas the classic 76 introgression lines (ILs) partitioned the genome into 106 mapping bins, the BILs generated 633 bins, thereby enhancing the mapping resolution of traits derived from the wild species. We demonstrate the power of the BILs for rapid fine mapping of simple and complex traits derived from the wild tomato species.

A Feedback-Insensitive Isopropylmalate Synthase Affects Acylsugar Composition in Cultivated and Wild Tomato.

Acylsugars are insecticidal specialized metabolites produced in the glandular trichomes of plants in the Solanaceae family. In the tomato clade of the Solanum genus, acylsugars consist of aliphatic acids of different chain lengths esterified to sucrose, or less frequently to glucose. Through liquid chromatography-mass spectrometry screening of introgression lines, we previously identified a region of chromosome 8 in the Solanum pennellii LA0716 genome (IL8-1/8-1-1) that causes the cultivated tomato Solanum lycopersicum to shift from producing acylsucroses with abundant 3-methylbutanoic acid acyl chains derived from leucine metabolism to 2-methylpropanoic acid acyl chains derived from valine metabolism. We describe multiple lines of evidence implicating a trichome-expressed gene from this region as playing a role in this shift. S. lycopersicum M82 SlIPMS3 (Solyc08g014230) encodes a functional end product inhibition-insensitive version of the committing enzyme of leucine biosynthesis, isopropylmalate synthase, missing the carboxyl-terminal 160 amino acids. In contrast, the S. pennellii LA0716 IPMS3 allele found in IL8-1/8-1-1 encodes a nonfunctional truncated IPMS protein. M82 transformed with an SlIPMS3 RNA interference construct exhibited an acylsugar profile similar to that of IL8-1-1, whereas the expression of SlIPMS3 in IL8-1-1 partially restored the M82 acylsugar phenotype. These IPMS3 alleles are polymorphic in 14 S. pennellii accessions spread throughout the geographical range of occurrence for this species and are associated with acylsugars containing varying amounts of 2-methylpropanoic acid and 3-methylbutanoic acid acyl chains.

Identification and mode of inheritance of quantitative trait loci for secondary metabolite abundance in tomato.

A large-scale metabolic quantitative trait loci (mQTL) analysis was performed on the well-characterized Solanum pennellii introgression lines to investigate the genomic regions associated with secondary metabolism in tomato fruit pericarp. In total, 679 mQTLs were detected across the 76 introgression lines. Heritability analyses revealed that mQTLs of secondary metabolism were less affected by environment than mQTLs of primary metabolism. Network analysis allowed us to assess the interconnectivity of primary and secondary metabolism as well as to compare and contrast their respective associations with morphological traits. Additionally, we applied a recently established real-time quantitative PCR platform to gain insight into transcriptional control mechanisms of a subset of the mQTLs, including those for hydroxycinnamates, acyl-sugar, naringenin chalcone, and a range of glycoalkaloids. Intriguingly, many of these compounds displayed a dominant-negative mode of inheritance, which is contrary to the conventional wisdom that secondary metabolite contents decreased on domestication. We additionally performed an exemplary evaluation of two candidate genes for glycolalkaloid mQTLs via the use of virus-induced gene silencing. The combined data of this study were compared with previous results on primary metabolism obtained from the same material and to other studies of natural variance of secondary metabolism.

The genome of the stress-tolerant wild tomato species Solanum pennellii.

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. 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. 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.

Resolution by recombination: breaking up Solanum pennellii introgressions.

Quantitative trait locus (QTL) genetics retains an important role in the study of biological and agronomic processes; however, its genetic resolution is often comparatively low. Community-based strategies are thus required to address this issue. Here we detail such a strategy wherein the widely used Solanum pennellii introgression lines (ILs) in the genetic background of the cultivated tomato (Solanum lycopersicum) are broken up into molecular marker-defined sublines as a community resource for map-based cloning.

Identification of enzyme activity quantitative trait loci in a Solanum lycopersicum x Solanum pennellii introgression line population.

Activities of 28 enzymes from central carbon metabolism were measured in pericarp tissue of ripe tomato fruits from field trials with an introgression line (IL) population generated by introgressing segments of the genome of the wild relative Solanum pennellii (LA0716) into the modern tomato cultivar Solanum lycopersicum M82. Enzyme activities were determined using a robotized platform in optimized conditions, where the activities largely reflect the level of the corresponding proteins. Two experiments were analyzed from years with markedly different climate conditions. A total of 27 quantitative trait loci were shared in both experiments. Most resulted in increased enzyme activity when a portion of the S. lycopersicum genome was substituted with the corresponding portion of the genome of S. pennellii. This reflects the change in activity between the two parental genotypes. The mode of inheritance was studied in a heterozygote IL population. A similar proportion of quantitative trait loci (approximately 30%) showed additive, recessive, and dominant modes of inheritance, with only 5% showing overdominance. Comparison with the location of putative genes for the corresponding proteins indicates a large role of trans-regulatory mechanisms. These results point to the genetic control of individual enzyme activities being under the control of a complex program that is dominated by a network of trans-acting genes.

The making of a compound inflorescence in tomato and related nightshades.

Variation in the branching of plant inflorescences determines flower number and, consequently, reproductive success and crop yield. Nightshade (Solanaceae) species are models for a widespread, yet poorly understood, program of eudicot growth, where short side branches are initiated upon floral termination. This "sympodial" program produces the few-flowered tomato inflorescence, but the classical mutants compound inflorescence (s) and anantha (an) are highly branched, and s bears hundreds of flowers. Here we show that S and AN, which encode a homeobox transcription factor and an F-box protein, respectively, control inflorescence architecture by promoting successive stages in the progression of an inflorescence meristem to floral specification. S and AN are sequentially expressed during this gradual phase transition, and the loss of either gene delays flower formation, resulting in additional branching. Independently arisen alleles of s account for inflorescence variation among domesticated tomatoes, and an stimulates branching in pepper plants that normally have solitary flowers. Our results suggest that variation of Solanaceae inflorescences is modulated through temporal changes in the acquisition of floral fate, providing a flexible evolutionary mechanism to elaborate sympodial inflorescence shoots.

The SOL Genomics Network: a comparative resource for Solanaceae biology and beyond.

The SOL Genomics Network (SGN; http://sgn.cornell.edu) is a rapidly evolving comparative resource for the plants of the Solanaceae family, which includes important crop and model plants such as potato (Solanum tuberosum), eggplant (Solanum melongena), pepper (Capsicum annuum), and tomato (Solanum lycopersicum). The aim of SGN is to relate these species to one another using a comparative genomics approach and to tie them to the other dicots through the fully sequenced genome of Arabidopsis (Arabidopsis thaliana). SGN currently houses map and marker data for Solanaceae species, a large expressed sequence tag collection with computationally derived unigene sets, an extensive database of phenotypic information for a mutagenized tomato population, and associated tools such as real-time quantitative trait loci. Recently, the International Solanaceae Project (SOL) was formed as an umbrella organization for Solanaceae research in over 30 countries to address important questions in plant biology. The first cornerstone of the SOL project is the sequencing of the entire euchromatic portion of the tomato genome. SGN is collaborating with other bioinformatics centers in building the bioinformatics infrastructure for the tomato sequencing project and implementing the bioinformatics strategy of the larger SOL project. The overarching goal of SGN is to make information available in an intuitive comparative format, thereby facilitating a systems approach to investigations into the basis of adaptation and phenotypic diversity in the Solanaceae family, other species in the Asterid clade such as coffee (Coffea arabica), Rubiaciae, and beyond.

Zooming in on a quantitative trait for tomato yield using interspecific introgressions.

To explore natural biodiversity we developed and examined introgression lines (ILs) containing chromosome segments of wild species (Solanum pennellii) in the background of the cultivated tomato (S. lycopersicum). We identified Brix9-2-5, which is a S. pennellii quantitative trait locus (QTL) that increases sugar yield of tomatoes and was mapped within a flower- and fruit-specific invertase (LIN5). QTL analysis representing five different tomato species delimited the functional polymorphism of Brix9-2-5 to an amino acid near the catalytic site of the invertase crystal, affecting enzyme kinetics and fruit sink strength. These results underline the power of diverse ILs for high-resolution perspectives on complex phenotypes.

Functional divergence of a syntenic invertase gene family in tomato, potato, and Arabidopsis.

Comparative analysis of complex developmental pathways depends on our ability to resolve the function of members of gene families across taxonomic groups. LIN5, which belongs to a small gene family of apoplastic invertases in tomato (Lycopersicon esculentum), is a quantitative trait locus that modifies fruit sugar composition. We have compared the genomic organization and expression of this gene family in the two distantly related species: tomato and Arabidopsis. Invertase family members reside on segmental duplications in the near-colinear genomes of tomato and potato (Solanum tuberosum). These chromosomal segments are syntenically duplicated in the model plant Arabidopsis. On the basis of phylogenetic analysis of genes in the microsyntenic region, we conclude that these segmental duplications arose independently after the separation of the tomato/potato clade from Arabidopsis. Rapid regulatory divergence is characteristic of the invertase family. Interestingly, although the processes of gene duplication and specialization of expression occurred separately in the two species, synteny-based orthologs from both clades acquired similar organ-specific expression. This similar expression pattern of the genes is evidence of comparable evolutionary constraints (parallel evolution) rather than of functional orthology. The observation that functional orthology cannot be identified through analysis of expression similarity highlights the caution that needs to be exercised in extrapolating developmental networks from a model organism.

O-methyltransferases involved in the biosynthesis of volatile phenolic derivatives in rose petals.

Rose (Rosa hybrida) flowers produce and emit a diverse array of volatiles, characteristic to their unique scent. One of the most prominent compounds in the floral volatiles of many rose varieties is the methoxylated phenolic derivative 3,5-dimethoxytoluene (orcinol dimethyl ether). Cell-free extracts derived from developing rose petals displayed O-methyltransferase (OMT) activities toward several phenolic substrates, including 3,5-dihydroxytoluene (orcinol), 3-methoxy,5-hydroxytoluene (orcinol monomethyl ether), 1-methoxy, 2-hydroxy benezene (guaiacol), and eugenol. The activity was most prominent in rose cv Golden Gate, a variety that produces relatively high levels of orcinol dimethyl ether, as compared with rose cv Fragrant Cloud, an otherwise scented variety but which emits almost no orcinol dimethyl ether. Using a functional genomics approach, we have identified and characterized two closely related cDNAs from a rose petal library that each encode a protein capable of methylating the penultimate and immediate precursors (orcinol and orcinol monomethyl ether, respectively) to give the final orcinol dimethyl ether product. The enzymes, designated orcinol OMTs (OOMT1 and OOMT2), are closely related to other plant methyltransferases whose substrates range from isoflavones to phenylpropenes. The peak in the levels of OOMT1 and OOMT2 transcripts in the flowers coincides with peak OMT activity and with the emission of orcinol dimethyl ether.