Genomic analysis of wild tomato introgressions determining metabolism- and yield-associated traits

Ectopic expression of BrIQD35 promotes drought stress tolerance in Nicotiana benthamiana

The plant IQD gene family is responsive to a variety of stresses. In this study, we studied the structural features and functions of the gene BrIQD35 in Chinese cabbage, a member of the IQD gene family. BrIQD35 was cloned and shown to contain an IQ motif. Transient expression of BrIQD35 indicated that it was localized on the plasma membrane and was significantly upregulated under drought and salt stress in Chinese cabbage. To further identify the function of BrIQD35, it was heterologously overexpressed in Nicotiana benthamiana. Although there was no significant difference between BrIQD35-overexpressed and wild-type (WT) plants under salt stress, WT N. benthamiana showed more wilting than the BrIQD35-overexpressed plants under drought stress. Since the IQ motif has been annotated as a CaM binding site, yeast two-hybrid assays were used to explore the interaction between BrIQD35 and CaM. The results indicated that BrIQD35 interacts weakly with CaMb, but not with CaMa, suggesting that BrIQD35 may function through the Ca2+ -CaMb pathway. The findings reveal a novel gene involved in drought tolerance, which is important for plant breeding and quality improvement for Chinese cabbage.

Chemical and structural quality traits during postharvest ripening regulated by chromosome segments from a wild relative of tomato Solanum pennellii IL4-2 and IL5-1

Tomato is usually harvested at an early ripening stage with high firmness suitable for storage and transportation but lacks many quality parameters such as sugars, organic acids, and phenolics. In a recent study, we have selected introgression lines (ILs) IL4-2 and IL5-1, developed from a cross between the Solanum pennellii and the Solanum lycopersicum M82, that exhibit differentiated postharvest shelf-life characteristics in the fruit compared to M82 and the rest of the ILs. Here, we first structurally and biochemically characterized IL4-2, IL5-1, and their parent M82 to decipher the cell wall mechanistic difference between soft (IL4-2) and firm (IL5-1) lines at two postharvest ripening periods. Generally, IL4-2 had more active cell wall modifications in terms of ripening-related gene expression, water-soluble pectin, and cell wall structure under the microscope, which probably makes this line softer than IL5-1. We also evaluated these lines based on commercial quality parameters, sugars, phenolics, organic, and amino acids to gain insight into their commercial and functional quality and reveal noticeable differences. In summary, the contribution of the S. pennellii IL5-1 and IL4-2 to the shelf life of the tomato was structurally characterized, and the component differences meeting the quality criteria were revealed.

The cytosolic invertase NI6 affects vegetative growth, flowering, fruit set, and yield in tomato

Sucrose metabolism is important for most plants, both as the main source of carbon and via signaling mechanisms that have been proposed for this molecule. A cleaving enzyme, invertase (INV) channels sucrose into sink metabolism. Although acid soluble and insoluble invertases have been largely investigated, studies on the role of neutral invertases (A/N-INV) have lagged behind. Here, we identified a tomato A/N-INV encoding gene (NI6) co-localizing with a previously reported quantitative trait locus (QTL) largely affecting primary carbon metabolism in tomato. Of the eight A/N-INV genes identified in the tomato genome, NI6 mRNA is present in all organs, but its expression was higher in sink tissues (mainly roots and fruits). A NI6-GFP fusion protein localized to the cytosol of mesophyll cells. Tomato NI6-silenced plants showed impaired growth phenotype, delayed flowering and a dramatic reduction in fruit set. Global gene expression and metabolite profile analyses of these plants revealed that NI6 is not only essential for sugar metabolism, but also plays a signaling role in stress adaptation. We also identified major hubs, whose expression patterns were greatly affected by NI6 silencing; these hubs were within the signaling cascade that coordinates carbohydrate metabolism with growth and development in tomato.

Genomic Dissection of a Wild Region in a Superior Solanum pennellii Introgression Sub-Line with High Ascorbic Acid Accumulation in Tomato Fruit

The Solanum pennellii introgression lines (ILs) have been exploited to map quantitative trait loci (QTLs) and identify favorable alleles that could improve fruit quality traits in tomato varieties. Over the past few years, ILs exhibiting increased content of ascorbic acid in the fruit have been selected, among which the sub-line R182. The aims of this work were to identify the genes of the wild donor S. pennellii harbored by the sub-line and to detect genes controlling ascorbic acid accumulation by using genomics tools. A Genotyping-By-Sequencing (GBS) approach confirmed that no wild introgressions were present in the sub-line besides one region on chromosome 7. By using a dense single nucleotide polymorphism (SNP) map obtained by RNA sequencing (RNA-Seq), the wild region of the sub-line was finely identified; thus, defining 39 wild genes that replaced 33 genes of the ILs genetic background (cv. M82). The differentially expressed genes mapping in the region and the variants detected among the cultivated and the wild alleles evidenced the potential role of the novel genes present in the wild region. Interestingly, one upregulated gene, annotated as a major facilitator superfamily protein, showed a novel structure in R182, with respect to the parental lines. These genes will be further investigated using gene editing strategies.

Single Primer Enrichment Technology (SPET) for High-Throughput Genotyping in Tomato and Eggplant Germplasm

Single primer enrichment technology (SPET) is a new, robust, and customizable solution for targeted genotyping. Unlike genotyping by sequencing (GBS), and like DNA chips, SPET is a targeted genotyping technology, relying on the sequencing of a region flanking a primer. Its reliance on single primers, rather than on primer pairs, greatly simplifies panel design, and allows higher levels of multiplexing than PCR-based genotyping. Thanks to the sequencing of the regions surrounding the target SNP, SPET allows the discovery of thousands of closely linked, novel SNPs. In order to assess the potential of SPET for high-throughput genotyping in plants, a panel comprising 5k target SNPs, designed both on coding regions and introns/UTRs, was developed for tomato and eggplant. Genotyping of two panels composed of 400 tomato and 422 eggplant accessions, comprising both domesticated material and wild relatives, generated a total of 12,002 and 30,731 high confidence SNPs, respectively, which comprised both target and novel SNPs in an approximate ratio of 1:1.6, and 1:5.5 in tomato and eggplant, respectively. The vast majority of the markers was transferrable to related species that diverged up to 3.4 million years ago (Solanum pennellii for tomato and S. macrocarpon for eggplant). Maximum Likelihood phylogenetic trees and PCA outputs obtained from the whole dataset highlighted genetic relationships among accessions and species which were congruent with what was previously reported in literature. Better discrimination among domesticated accessions was achieved by using the target SNPs, while better discrimination among wild species was achieved using the whole SNP dataset. Our results reveal that SPET genotyping is a robust, high-throughput technology for genetic fingerprinting, with a high degree of cross-transferability between crops and their cultivated and wild relatives, and allows identification of duplicates and mislabeled accessions in genebanks.

Phenolic composition and antioxidant properties of ex-situ conserved tomato (Solanum lycopersicum L.) germplasm

Tomato (Solanum lycopersicum L.) local varieties represent a reservoir of genetic diversity for desirable quality traits. In this study, a representative collection of table tomato germplasm conserved ex-situ in the Portuguese Gene Bank was characterized for its polyphenols composition and antioxidant capacity. Phenolic acids, such as caffeic and p-coumaric acids bounded to a hexose and 5-O-caffeoylquinic acid, corresponded to 71-98% of the identified phenolic compounds; while the remaining fraction consisted of quercetin and kaempferol glycoside derivatives. Among the studied tomato accessions, it was possible to identify those that stand out for the analysed bioactive traits. These findings highlighted the interest of using Portuguese tomato germplasm in breeding programs or of reintroducing into cultivation these local varieties used for fresh consumption.

Duplicated flavonoid 3'-hydroxylase and flavonoid 3', 5'-hydroxylase genes in barley genome

Background

Anthocyanin compounds playing multiple biological functions can be synthesized in different parts of barley (Hordeum vulgare L.) plant. The diversity of anthocyanin molecules is related with branching the pathway to alternative ways in which dihydroflavonols may be modified either with the help of flavonoid 3'-hydroxylase (F3'H) or flavonoid 3',5'-hydroxylase (F3'5'H)-the cytochrome P450-dependent monooxygenases. The F3'H and F3'5'H gene families are among the least studied anthocyanin biosynthesis structural genes in barley. The aim of this study was to identify and characterise duplicated copies of the F3'H and F3'5'H genes in the barley genome.

Results

Four copies of the F3'5'H gene (on chromosomes 4HL, 6HL, 6HS and 7HS) and two copies of the F3'H gene (on chromosomes 1HL and 6HS) were identified in barley genome. These copies have either one or two introns. Amino acid sequences analysis demonstrated the presence of the flavonoid hydroxylase-featured conserved motifs in all copies of the F3'H and F3'5'H genes with the exception of F3'5'H-3 carrying a loss-of-function mutation in a conservative cytochrome P450 domain. It was shown that the divergence between F3'H and F3'5'H genes occurred 129 million years ago (MYA) before the emergence of monocot and dicot plant species. The F3'H copy approximately occurred 80 MYA; the appearance of F3'5'H copies occurred 8, 36 and 91 MYA. qRT-PCR analysis revealed the tissue-specific activity for some copies of the studied genes. The F3'H-1 gene was transcribed in aleurone layer, lemma and pericarp (with an increased level in the coloured pericarp), whereas the F3'H-2 gene was expressed in stems only. The F3'5'H-1 gene was expressed only in the aleurone layer, and in a coloured aleurone its expression was 30-fold higher. The transcriptional activity of F3'5'H-2 was detected in different tissues with significantly higher level in uncoloured genotype in contrast to coloured ones. The F3'5'H-3 gene expressed neither in stems nor in aleurone layer, lemma and pericarp. The F3'5'H-4 gene copy was weakly expressed in all tissues analysed.

Conclusion

F3'H and F3'5'H-coding genes involved in anthocyanin synthesis in H. vulgare were identified and characterised, from which the copies designated F3'H-1, F3'H-2, F3'5'H-1 and F3'5'H-2 demonstrated tissue-specific expression patterns. Information on these modulators of the anthocyanin biosynthesis pathway can be used in future for manipulation with synthesis of diverse anthocyanin compounds in different parts of barley plant. Finding both the copies with tissue-specific expression and a copy undergoing pseudogenization demonstrated rapid evolutionary events tightly related with functional specialization of the duplicated members of the cytochrome P450-dependent monooxygenases gene families.

New insights into tomato microRNAs

Cultivated tomato, Solanum lycopersicum, is one of the most common fruits in the global food industry. Together with the wild tomato Solanum pennellii, it is widely used for developing better cultivars. MicroRNAs affect mRNA regulation, inhibiting its translation and/or promoting its degradation. Important proteins involved in these processes are ARGONAUTE and DICER. This study aimed to identify and characterize the genes involved in the miRNA processing pathway, miRNA molecules and target genes in both species. We validated the presence of pathway genes and miRNA in different NGS libraries and 6 miRNA families using quantitative RT-PCR. We identified 71 putative proteins in S. lycopersicum and 108 in S. pennellii likely involved in small RNAs processing. Of these, 29 and 32 participate in miRNA processing pathways, respectively. We identified 343 mature miRNAs, 226 pre-miRNAs in 87 families, including 192 miRNAs, which were not previously identified, belonging to 38 new families in S. lycopersicum. In S. pennellii, we found 388 mature miRNAs and 234 pre-miRNAs contained in 85 families. All miRNAs found in S. pennellii were unpublished, being identified for the first time in our study. Furthermore, we identified 2471 and 3462 different miRNA target in S. lycopersicum and S. pennellii, respectively.

Evolution by duplication: paleopolyploidy events in plants reconstructed by deciphering the evolutionary history of VOZ transcription factors

BACKGROUND

Facilitated by the rapid progress of sequencing technology, comparative genomic studies in plants have unveiled recurrent whole genome duplication (i.e. polyploidization) events throughout plant evolution. The evolutionary past of plant genes should be analyzed in a background of recurrent polyploidy events in distinctive plant lineages. The Vascular Plant One Zinc-finger (VOZ) gene family encode transcription factors associated with a number of important traits including control of flowering time and photoperiodic pathways, but the evolutionary trajectory of this gene family remains uncharacterized.

RESULTS

In this study, we deciphered the evolutionary history of the VOZ gene family by analyses of 107 VOZ genes in 46 plant genomes using integrated methods: phylogenic reconstruction, Ks-based age estimation and genomic synteny comparisons. By scrutinizing the VOZ gene family phylogeny the core eudicot γ event was well circumscribed, and relics of the precommelinid τ duplication event were detected by incorporating genes from oil palm and banana. The more recent T and ρ polyploidy events, closely coincident with the species diversification in Solanaceae and Poaceae, respectively, were also identified. Other important polyploidy events captured included the "salicoid" event in poplar and willow, the "early legume" and "soybean specific" events in soybean, as well as the recent polyploidy event in Physcomitrella patens. Although a small transcription factor gene family, the evolutionary history of VOZ genes provided an outstanding record of polyploidy events in plants. The evolutionary past of VOZ gene family demonstrated a close correlation with critical plant polyploidy events which generated species diversification and provided answer to Darwin's "abominable mystery".

CONCLUSIONS

We deciphered the evolutionary history of VOZ transcription factor family in plants and ancestral polyploidy events in plants were recapitulated simultaneously. This analysis allowed for the generation of an idealized plant gene tree demonstrating distinctive retention and fractionation patterns following polyploidy events.

Regulatory Divergence in Wound-Responsive Gene Expression between Domesticated and Wild Tomato

The evolution of transcriptional regulatory mechanisms is central to how stress response and tolerance differ between species. However, it remains largely unknown how divergence in cis-regulatory sites and, subsequently, transcription factor (TF) binding specificity contribute to stress-responsive expression divergence, particularly between wild and domesticated species. By profiling wound-responsive gene transcriptomes in wild Solanum pennellii and domesticated S. lycopersicum, we found extensive wound response divergence and identified 493 S. lycopersicum and 278 S. pennellii putative cis-regulatory elements (pCREs) that were predictive of wound-responsive gene expression. Only 24-52% of these wound response pCREs (depending on wound response patterns) were consistently enriched in the putative promoter regions of wound-responsive genes across species. In addition, between these two species, their differences in pCRE site sequences were significantly and positively correlated with differences in wound-responsive gene expression. Furthermore, ∼11-39% of pCREs were specific to only one of the species and likely bound by TFs from different families. These findings indicate substantial regulatory divergence in these two plant species that diverged ∼3-7 million years ago. Our study provides insights into the mechanistic basis of how the transcriptional response to wounding is regulated and, importantly, the contribution of cis-regulatory components to variation in wound-responsive gene expression between a wild and a domesticated plant species.

The genetic architecture of photosynthesis and plant growth-related traits in tomato

To identify genomic regions involved in the regulation of fundamental physiological processes such as photosynthesis and respiration, a population of Solanum pennellii introgression lines was analyzed. We determined phenotypes for physiological, metabolic, and growth related traits, including gas exchange and chlorophyll fluorescence parameters. Data analysis allowed the identification of 208 physiological and metabolic quantitative trait loci with 33 of these being associated to smaller intervals of the genomic regions, termed BINs. Eight BINs were identified that were associated with higher assimilation rates than the recurrent parent M82. Two and 10 genomic regions were related to shoot and root dry matter accumulation, respectively. Nine genomic regions were associated with starch levels, whereas 12 BINs were associated with the levels of other metabolites. Additionally, a comprehensive and detailed annotation of the genomic regions spanning these quantitative trait loci allowed us to identify 87 candidate genes that putatively control the investigated traits. We confirmed 8 of these at the level of variance in gene expression. Taken together, our results allowed the identification of candidate genes that most likely regulate photosynthesis, primary metabolism, and plant growth and as such provide new avenues for crop improvement.

Back into the wild-Apply untapped genetic diversity of wild relatives for crop improvement

Deleterious effects of climate change and human activities, as well as diverse environmental stresses, present critical challenges to food production and the maintenance of natural diversity. These challenges may be met by the development of novel crop varieties with increased biotic or abiotic resistance that enables them to thrive in marginal lands. However, considering the diverse interactions between crops and environmental factors, it is surprising that evolutionary principles have been underexploited in addressing these food and environmental challenges. Compared with domesticated cultivars, crop wild relatives (CWRs) have been challenged in natural environments for thousands of years and maintain a much higher level of genetic diversity. In this review, we highlight the significance of CWRs for crop improvement by providing examples of CWRs that have been used to increase biotic and abiotic stress resistance/tolerance and overall yield in various crop species. We also discuss the surge of advanced biotechnologies, such as next-generation sequencing technologies and omics, with particular emphasis on how they have facilitated gene discovery in CWRs. We end the review by discussing the available resources and conservation of CWRs, including the urgent need for CWR prioritization and collection to ensure continuous crop improvement for food sustainability.

Phytochrome Interacting Factors (PIFs) in Solanum lycopersicum: Diversity, Evolutionary History and Expression Profiling during Different Developmental Processes

Although the importance of light for tomato plant yield and edible fruit quality is well known, the PHYTOCHROME INTERACTING FACTORS (PIFs), main components of phytochrome-mediated light signal transduction, have been studied almost exclusively in Arabidopsis thaliana. Here, the diversity, evolution and expression profile of PIF gene subfamily in Solanum lycopersicum was characterized. Eight tomato PIF loci were identified, named SlPIF1a, SlPIF1b, SlPIF3, SlPIF4, SlPIF7a, SlPIF7b, SlPIF8a and SlPIF8b. The duplication of SlPIF1, SlPIF7 and SlPIF8 genes were dated and temporally coincided with the whole-genome triplication event that preceded tomato and potato divergence. Different patterns of mRNA accumulation in response to light treatments were observed during seedling deetiolation, dark-induced senescence, diel cycle and fruit ripening. SlPIF4 showed similar expression profile as that reported for A. thaliana homologs, indicating an evolutionary conserved function of PIF4 clade. A comprehensive analysis of the evolutionary and transcriptional data allowed proposing that duplicated SlPIFs have undergone sub- and neofunctionalization at mRNA level, pinpointing the importance of transcriptional regulation for the maintenance of duplicated genes. Altogether, the results indicate that genome polyploidization and functional divergence have played a major role in diversification of the Solanum PIF gene subfamily.

Genome-wide identification and expression analysis of the IQD gene family in moso bamboo (Phyllostachys edulis)

Members of the plant-specific IQ67-domain (IQD) protein family are involved in various aspects of normal plant growth and developmental processes as well as basal defence response. Although hundreds of IQD proteins have been identified, only a small number of IQDs have been functionally characterized. Moreover, no systematic study has been performed on moso bamboo. In this study, we performed for the first time a genome-wide identification and expression analysis of the IQD gene family in moso bamboo. We identified 29 non-redundant PeIQD encoding genes. Analysis of the evolutionary patterns and divergence revealed that the IQD genes underwent a large-scale event around 12 million years ago and the division times of IQD family genes between moso bamboo and rice, and, between moso bamboo and Brachypodium, were found to be 20-35 MYA and 25-40 MYA, respectively. We surveyed the putative promoter regions of the PeIQD genes, which showed that largely stress-related cis-elements existed in these genes. The expression profiles of the IQD genes shed light on their functional divergence. Additionally, a yeast two-hybrid assay proved that PeIQD8 can interact with PeCaM2 and that IQ or I in the IQ motif is required for PeIQD8 to combine with CaM2.

Hydrophilic antioxidants from Andean tomato landraces assessed by their bioactivities in vitro and in vivo

Potential nutraceutical properties of hydrophilic antioxidants in fruits of tomato landraces collected in Andean valleys were characterised. Antioxidant metabolites were measured by HPLC-DAD-MS/MS in mature fruits and their biological activities were assessed by in vitro and in vivo methods. In vitro antioxidant capacities were established by TEAC and FRAP methods. For in vivo biological activities we used a procedure based on Caenorhabditis elegans subjected to thermal stress. In addition, Saccharomyces cerevisiae was also used as a rapid screening system to evaluate tomato antioxidant capacity. All tomato accessions displayed significant differences regarding metabolic composition, biological activity and antioxidant capacity. Metabolite composition was associated with geographical origin and fruit size. Antioxidant activities showed significant association with phenolic compounds, such as caffeoylquinic acids, ferulic acid-O-hexosides and rutin. Combination of in vitro and in vivo methods applied here allowed evaluation of the variability in nutraceutical properties of tomato landraces, which could be applied to other fruits or food products.

Constraint around Quarter-Power Allometric Scaling in Wild Tomatoes (Solanum sect. Lycopersicon; Solanaceae)

The West-Brown-Enquist (WBE) metabolic scaling theory posits that many organismal features scale predictably with body size because of selection to minimize transport costs in resource distribution networks. Many scaling exponents are quarter-powers, as predicted by WBE, but there are also biologically significant deviations that could reflect adaptation to different environments. A central but untested prediction of the WBE model is that wide deviation from optimal scaling is penalized, leading to a pattern of constraint on scaling exponents. Here, we demonstrate, using phylogenetic comparative methods, that variation in allometric scaling between mass and leaf area across 17 wild tomato taxa is constrained around a value indistinguishable from that predicted by WBE but significantly greater than 2/3 (geometric-similarity model). The allometric-scaling exponent was highly correlated with fecundity, water use, and drought response, suggesting that it is functionally significant and therefore could be under selective constraints. However, scaling was not strictly log-log linear but rather declined during ontogeny in all species, as has been observed in many plant species. We caution that although our results supported one prediction of the WBE model, it did not strongly test the model in other important respects. Nevertheless, phylogenetic comparative methods such as those used here are powerful but underutilized tools for metabolic ecology that complement existing methods to adjudicate between models.

Molecular marker applications in plants

Individuals within a population of a sexually reproducing species will have some degree of heritable genomic variation caused by mutations, insertion/deletions (INDELS), inversions, duplications, and translocations. Such variation can be detected and screened using molecular, or genetic, markers. By definition, molecular markers are genetic loci that can be easily tracked and quantified in a population and may be associated with a particular gene or trait of interest. This chapter will review the current major applications of molecular markers in plants.

The IQD gene family in soybean: structure, phylogeny, evolution and expression

Members of the plant-specific IQ67-domain (IQD) protein family are involved in plant development and the basal defense response. Although systematic characterization of this family has been carried out in Arabidopsis, tomato (Solanum lycopersicum), Brachypodium distachyon and rice (Oryza sativa), systematic analysis and expression profiling of this gene family in soybean (Glycine max) have not previously been reported. In this study, we identified and structurally characterized IQD genes in the soybean genome. A complete set of 67 soybean IQD genes (GmIQD1-67) was identified using Blast search tools, and the genes were clustered into four subfamilies (IQD I-IV) based on phylogeny. These soybean IQD genes are distributed unevenly across all 20 chromosomes, with 30 segmental duplication events, suggesting that segmental duplication has played a major role in the expansion of the soybean IQD gene family. Analysis of the Ka/Ks ratios showed that the duplicated genes of the GmIQD family primarily underwent purifying selection. Microsynteny was detected in most pairs: genes in clade 1-3 might be present in genome regions that were inverted, expanded or contracted after the divergence; most gene pairs in clade 4 showed high conservation with little rearrangement among these gene-residing regions. Of the soybean IQD genes examined, six were most highly expressed in young leaves, six in flowers, one in roots and two in nodules. Our qRT-PCR analysis of 24 soybean IQD III genes confirmed that these genes are regulated by MeJA stress. Our findings present a comprehensive overview of the soybean IQD gene family and provide insights into the evolution of this family. In addition, this work lays a solid foundation for further experiments aimed at determining the biological functions of soybean IQD genes in growth and development.

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.

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

Björn Usadel and colleagues report the genome sequence of the wild tomato species Solanum pennellii. The authors identify genes important for stress tolerance, metabolism and fruit maturation and suggest that transposable elements have had an important role in the evolution of the S. penellii stress response.

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.

Natural occurring epialleles determine vitamin E accumulation in tomato fruits

Vitamin E (VTE) content is a low heritability nutritional trait for which the genetic determinants are poorly understood. Here, we focus on a previously detected major tomato VTE quantitative trait loci (QTL; mQTL(9-2-6)) and identify the causal gene as one encoding a 2-methyl-6-phytylquinol methyltransferase (namely VTE3(1)) that catalyses one of the final steps in the biosynthesis of γ- and α-tocopherols, which are the main forms of VTE. By reverse genetic approaches, expression analyses, siRNA profiling and DNA methylation assays, we demonstrate that mQTL(9-2-6) is an expression QTL associated with differential methylation of a SINE retrotransposon located in the promoter region of VTE3(1). Promoter DNA methylation can be spontaneously reverted leading to different epialleles affecting VTE3(1) expression and VTE content in fruits. These findings indicate therefore that naturally occurring epialleles are responsible for regulation of a nutritionally important metabolic QTL and provide direct evidence of a role for epigenetics in the determination of agronomic traits.

Morphological and anatomical determinants of mesophyll conductance in wild relatives of tomato (Solanum sect. Lycopersicon, sect. Lycopersicoides; Solanaceae)

Natural selection on photosynthetic performance is a primary factor determining leaf phenotypes. The complex CO2 diffusion path from substomatal cavities to the chloroplasts - the mesophyll conductance (g(m)) - limits photosynthetic rate in many species and hence shapes variation in leaf morphology and anatomy. Among sclerophyllous and succulent taxa, structural investment in leaves, measured as the leaf dry mass per area (LMA), has been implicated in decreased gm . However, in herbaceous taxa with high g(m), it is less certain how LMA impacts CO2 diffusion and whether it significantly affects photosynthetic performance. We addressed these questions in the context of understanding the ecophysiological significance of leaf trait variation in wild tomatoes, a closely related group of herbaceous perennials. Although g(m) was high in wild tomatoes, variation in g(m) significantly affected photosynthesis. Even in these tender-leaved herbaceous species, greater LMA led to reduced g(m). This relationship between g(m) and LMA is partially mediated by cell packing and leaf thickness, although amphistomy (equal distribution of stomata on both sides of the leaf) mitigates the effect of leaf thickness. Understanding the costs of increased LMA will inform future work on the adaptive significance of leaf trait variation across ecological gradients in wild tomatoes and other systems.

Silencing of the tomato sugar partitioning affecting protein (SPA) modifies sink strength through a shift in leaf sugar metabolism

Limitations in our understanding about the mechanisms that underlie source-sink assimilate partitioning are increasingly becoming a major hurdle for crop yield enhancement via metabolic engineering. By means of a comprehensive approach, this work reports the functional characterization of a DnaJ chaperone related-protein (named as SPA; sugar partition-affecting) that is involved in assimilate partitioning in tomato plants. SPA protein was found to be targeted to the chloroplast thylakoid membranes. SPA-RNAi tomato plants produced more and heavier fruits compared with controls, thus resulting in a considerable increment in harvest index. The transgenic plants also displayed increased pigment levels and reduced sucrose, glucose and fructose contents in leaves. Detailed metabolic and enzymatic activities analyses showed that sugar phosphate intermediates were increased while the activity of phosphoglucomutase, sugar kinases and invertases was reduced in the photosynthetic organs of the silenced plants. These changes would be anticipated to promote carbon export from foliar tissues. The combined results suggested that the tomato SPA protein plays an important role in plastid metabolism and mediates the source-sink relationships by affecting the rate of carbon translocation to fruits.

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.

Using a periclinal chimera to unravel layer-specific gene expression in plants

Plant organs are made from multiple cell types, and defining the expression level of a gene in any one cell or group of cells from a complex mixture is difficult. Dicotyledonous plants normally have three distinct layers of cells, L1, L2 and L3. Layer L1 is the single layer of cells making up the epidermis, layer L2 the single cell sub-epidermal layer and layer L3 constitutes the rest of the internal cells. Here we show how it is possible to harvest an organ and characterise the level of layer-specific expression by using a periclinal chimera that has its L1 layer from Solanum pennellii and its L2 and L3 layers from Solanum lycopersicum. This is possible by measuring the level of the frequency of species-specific transcripts. RNA-seq analysis enabled the genome-wide assessment of whether a gene is expressed in the L1 or L2/L3 layers. From 13 277 genes that are expressed in both the chimera and the parental lines and with at least one polymorphism between the parental alleles, we identified 382 genes that are preferentially expressed in L1 in contrast to 1159 genes in L2/L3. Gene ontology analysis shows that many genes preferentially expressed in L1 are involved in cutin and wax biosynthesis, whereas numerous genes that are preferentially expressed in L2/L3 tissue are associated with chloroplastic processes. These data indicate the use of such chimeras and provide detailed information on the level of layer-specific expression of genes.

TARE1, a mutated Copia-like LTR retrotransposon followed by recent massive amplification in tomato

Long terminal repeat retrotransposons (LTR-RTs) are the major DNA components in flowering plants. Most LTR-RTs contain dinucleotides ‘TG’ and ‘CA’ at the ends of the two LTRs. Here we report the structure, evolution, and propensity of a tomato atypical retrotransposon element (TARE1) with both LTRs starting as ‘TA’. This family is also characterized by high copy numbers (354 copies), short LTR size (194 bp), extremely low ratio of solo LTRs to intact elements (0.05∶1), recent insertion (most within 0.75∼1.75 million years, Mys), and enrichment in pericentromeric region. The majority (83%) of the TARE1 elements are shared between S. lycopersicum and its wild relative S. pimpinellifolium, but none of them are found in potato. In the present study, we used shared LTR-RTs as molecular markers and estimated the divergence time between S. lycopersicum and S. pimpinellifolium to be <0.5 Mys. Phylogenetic analysis showed that the TARE1 elements, together with two closely related families, TARE2 and TGRE1, have formed a sub-lineage belonging to a Copia-like Ale lineage. Although TARE1 and TARE2 shared similar structural characteristics, the timing, scale, and activity of their amplification were found to be substantially different. We further propose a model wherein a single mutation from ‘G’ to ‘A’ in 3′ LTR followed by amplification is responsible for the origin of TARE1, thus providing evidence that the proliferation of a spontaneous mutation can be mediated by the amplification of LTR-RTs at the level of RNA.

Linking the potato genome to the conserved ortholog set (COS) markers

BACKGROUND

Conserved ortholog set (COS) markers are an important functional genomics resource that has greatly improved orthology detection in Asterid species. A comprehensive list of these markers is available at Sol Genomics Network (http://solgenomics.net/) and many of these have been placed on the genetic maps of a number of solanaceous species.

RESULTS

We amplified over 300 COS markers from eight potato accessions involving two diploid landraces of Solanum tuberosum Andigenum group (formerly classified as S. goniocalyx, S. phureja), and a dihaploid clone derived from a modern tetraploid cultivar of S. tuberosum and the wild species S. berthaultii, S. chomatophilum, and S. paucissectum. By BLASTn (Basic Local Alignment Search Tool of the NCBI, National Center for Biotechnology Information) algorithm we mapped the DNA sequences of these markers into the potato genome sequence. Additionally, we mapped a subset of these markers genetically in potato and present a comparison between the physical and genetic locations of these markers in potato and in comparison with the genetic location in tomato. We found that most of the COS markers are single-copy in the reference genome of potato and that the genetic location in tomato and physical location in potato sequence are mostly in agreement. However, we did find some COS markers that are present in multiple copies and those that map in unexpected locations. Sequence comparisons between species show that some of these markers may be paralogs.

CONCLUSIONS

The sequence-based physical map becomes helpful in identification of markers for traits of interest thereby reducing the number of markers to be tested for applications like marker assisted selection, diversity, and phylogenetic studies.

Metabolic engineering of tomato fruit organic acid content guided by biochemical analysis of an introgression line

Organic acid content is regarded as one of the most important quality traits of fresh tomato (Solanum lycopersicum). However, the complexity of carboxylic acid metabolism and storage means that it is difficult to predict the best way to engineer altered carboxylic acid levels. Here, we used a biochemical analysis of a tomato introgression line with increased levels of fruit citrate and malate at breaker stage to identify a metabolic engineering target that was subsequently tested in transgenic plants. Increased carboxylic acid levels in introgression line 2-5 were not accompanied by changes in the pattern of carbohydrate oxidation by pericarp discs or the catalytic capacity of tricarboxylic acid cycle enzymes measured in isolated mitochondria. However, there was a significant decrease in the maximum catalytic activity of aconitase in total tissue extracts, suggesting that a cytosolic isoform of aconitase was affected. To test the role of cytosolic aconitase in controlling fruit citrate levels, we analyzed fruit of transgenic lines expressing an antisense construct against SlAco3b, one of the two tomato genes encoding aconitase. A green fluorescent protein fusion of SlAco3b was dual targeted to cytosol and mitochondria, while the other aconitase, SlAco3a, was exclusively mitochondrial when transiently expressed in tobacco (Nicotiana tabacum) leaves. Both aconitase transcripts were decreased in fruit from transgenic lines, and aconitase activity was reduced by about 30% in the transgenic lines. Other measured enzymes of carboxylic acid metabolism were not significantly altered. Both citrate and malate levels were increased in ripe fruit of the transgenic plants, and as a consequence, total carboxylic acid content was increased by 50% at maturity.

Role of an esterase in flavor volatile variation within the tomato clade

Tomato flavor is dependent upon a complex mixture of volatiles including multiple acetate esters. Red-fruited species of the tomato clade accumulate a relatively low content of acetate esters in comparison with the green-fruited species. We show that the difference in volatile ester content between the red- and green-fruited species is associated with insertion of a retrotransposon adjacent to the most enzymatically active member of a family of esterases. This insertion causes higher expression of the esterase, resulting in the reduced levels of multiple esters that are negatively correlated with human preferences for tomato. The insertion was evolutionarily fixed in the red-fruited species, suggesting that high expression of the esterase and consequent low ester content may provide an adaptive advantage in the ancestor of the red-fruited species. These results illustrate at a molecular level how closely related species exhibit major differences in volatile production by altering a volatile-associated catabolic activity.

Novel glutamate dehydrogenase genes show increased transcript and protein abundances in mature tomato fruits

NAD(P)H-glutamate dehydrogenase (GDH, EC 1.4.1.3) contributes to the control of glutamate homeostasis in all living organisms. In bacteria and animals, GDH is a homohexamer allosterically regulated, whereas in plants NADH-GDH (EC 1.4.1.2) is also found as heterohexamer of α- and β-subunits, but its regulation remains undefined. In tomato (Solanum lycopersicum), GDH activity increases during the fruit ripening along with the content of free glutamate, the most abundant amino acid of ripe fruit involved in conferring the genuine tomato flavour. In this work, novel Slgdh-NAD genes were identified in the recently deciphered tomato genome: three encoding the α-subunit (Slgdh-NAD;A1-3) and one additional gene encoding the β-subunit of GDH (Slgdh-NAD;B1) isolated from a genomic library. These genes are located in different chromosomes. Slgdh-NAD;A1-3 show conserved structures, whereas Slgdh-NAD;B1 includes a novel 5'-untranslated exon. Slgdh-NAD;A1-3 transcripts were detected in all tomato tissues examined, showing the highest levels in mature green fruits, contrasting with Slgdh-NAD;B1 transcripts which were detected mainly in roots or in mature fruits when treated with glutamate, NaCl or salicylic acid. Analyses of GDH activity and protein distribution in different tissues of the Micro-Tom cultivar showed that only the active homohexamer of GDH β-subunits was detected in roots while heterohexamers of GDH α- and β-subunits were found in fruits. These results indicate that GDH β-subunit could modulate the heteromeric isoforms of GDH in response to the environment and physiology of the tomato fruit. This information is relevant to manipulate glutamate contents in tomato fruits genetically.

Os11Gsk gene from a wild rice, Oryza rufipogon improves yield in rice

Chromosomal segments from wild rice species Oryza rufipogon, introgressed into an elite indica rice restorer line (KMR3) using molecular markers, resulted in significant increase in yield. Here we report the transcriptome analysis of flag leaves and fully emerged young panicles of one of the high yielding introgression lines IL50-7 in comparison to KMR3. A 66-fold upregulated gene Os11Gsk, which showed no transcript in KMR3 was highly expressed in O. rufipogon and IL50-7. A 5-kb genomic region including Os11Gsk and its flanking regions could be PCR amplified only from IL50-7, O. rufipogon, japonica varieties of rice-Nipponbare and Kitaake but not from the indica varieties, KMR3 and Taichung Native-1. Three sister lines of IL50-7 yielding higher than KMR3 showed presence of Os11Gsk, whereas the gene was absent in three other ILs from the same cross having lower yield than KMR3, indicating an association of the presence of Os11Gsk with high yield. Southern analysis showed additional bands in the genomic DNA of O. rufipogon and IL50-7 with Os11Gsk probe. Genomic sequence analysis of ten highly co-expressed differentially regulated genes revealed that two upregulated genes in IL50-7 were derived from O. rufipogon and most of the downregulated genes were either from KMR3 or common to KMR3, IL50-7, and O. rufipogon. Thus, we show that Os11Gsk is a wild rice-derived gene introduced in KMR3 background and increases yield either by regulating expression of functional genes sharing homology with it or by causing epigenetic modifications in the introgression line.

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.

Genetic dissection of vitamin E biosynthesis in tomato

Vegetables are critical for human health as they are a source of multiple vitamins including vitamin E (VTE). In plants, the synthesis of VTE compounds, tocopherol and tocotrienol, derives from precursors of the shikimate and methylerythritol phosphate pathways. Quantitative trait loci (QTL) for α-tocopherol content in ripe fruit have previously been determined in an Solanum pennellii tomato introgression line population. In this work, variations of tocopherol isoforms (α, β, γ, and δ) in ripe fruits of these lines were studied. In parallel all tomato genes structurally associated with VTE biosynthesis were identified and mapped. Previously identified VTE QTL on chromosomes 6 and 9 were confirmed whilst novel ones were identified on chromosomes 7 and 8. Integrated analysis at the metabolic, genetic and genomic levels allowed us to propose 16 candidate loci putatively affecting tocopherol content in tomato. A comparative analysis revealed polymorphisms at nucleotide and amino acid levels between Solanum lycopersicum and S. pennellii candidate alleles. Moreover, evolutionary analyses showed the presence of codons evolving under both neutral and positive selection, which may explain the phenotypic differences between species. These data represent an important step in understanding the genetic determinants of VTE natural variation in tomato fruit and as such in the ability to improve the content of this important nutriceutical.

The genetic architecture of branched-chain amino acid accumulation in tomato fruits

Previous studies of the genetic architecture of fruit metabolic composition have allowed us to identify four strongly conserved co-ordinate quantitative trait loci (QTL) for the branched-chain amino acids (BCAAs). This study has been extended here to encompass the other 23 enzymes described to be involved in the pathways of BCAA synthesis and degradation. On coarse mapping the chromosomal location of these enzymes, it was possible to define the map position of 24 genes. Of these genes eight co-localized, or mapped close to BCAA QTL including those encoding ketol-acid reductoisomerase (KARI), dihydroxy-acid dehydratase (DHAD), and isopropylmalate dehydratase (IPMD). Quantitative evaluation of the expression levels of these genes revealed that the S. pennellii allele of IPMD demonstrated changes in the expression level of this gene, whereas those of KARI and DHAD were invariant across the genotypes. Whilst the antisense inhibition of IPMD resulted in increased BCAA, the antisense inhibition of neither KARI nor DHAD produced a clear effect in fruit BCAA contents. The results are discussed both with respect to the roles of these specific enzymes within plant amino acid metabolism and within the context of current understanding of the regulation of plant branched-chain amino acid metabolism.

The influence of fruit load on the tomato pericarp metabolome in a Solanum chmielewskii introgression line population

It has been recently demonstrated, utilizing interspecific introgression lines of tomato, generated from the cross between Solanum lycopersicum and the wild species Solanum pennellii, that the efficiency of photosynthate partitioning exerts a considerable influence on the metabolic composition of tomato fruit pericarp. In order to further evaluate the influence of source-sink interaction, metabolite composition was determined by gas chromatography-mass spectrometry in a different population. For this purpose, we used 23 introgression lines resulting from an interspecific cross between S. lycopersicum and the wild species Solanum chmielewskii under high (unpruned trusses) and low (trusses pruned to one fruit) fruit load conditions. Following this strategy, we were able to contrast the metabolite composition of fruits from plants cultivated at both fruit loads as well as to compare the network behavior of primary metabolism in the introgression line population. The study revealed that while a greater number of metabolic quantitative trait loci were observed under high fruit load (240) than under low fruit load (128) cultivations, the levels of metabolites were more highly correlated under low fruit load cultivation. Finally, an analysis of genotype × fruit load interactions indicated a greater influence of development and cultivation than genotype on fruit composition. Comparison with previously documented transcript profiles from a subset of these lines revealed that changes in metabolite levels did not correlate with changes in the levels of genes associated with their metabolism. These findings are discussed in the context of our current understanding of the genetic and environmental influence on metabolic source-sink interactions in tomato, with particular emphasis given to fruit amino acid content.

Reciprocal insights into adaptation from agricultural and evolutionary studies in tomato

Although traditionally separated by different aims and methodologies, research on agricultural and evolutionary problems shares a common goal of understanding the mechanisms underlying functionally important traits. As such, research in both fields offers potential complementary and reciprocal insights. Here, we discuss adaptive stress responses (specifically to water stress) as an example of potentially fruitful research reciprocity, where agricultural research has clearly produced advances that could benefit evolutionary studies, while evolutionary studies offer approaches and insights underexplored in crop studies. We focus on research on Solanum species that include the domesticated tomato and its wild relatives. Integrated approaches to understanding ecological adaptation are particularly attractive in tomato and its wild relatives: many presumptively adaptive phenotypic differences characterize wild species, and the physiological and mechanistic basis of many relevant traits and environmental responses has already been examined in the context of cultivated tomato and some wild species. We highlight four specific instances where these reciprocal insights can be combined to better address questions that are fundamental both to agriculture and evolution.