Population structure and genetic differentiation associated with breeding history and selection in tomato (Solanum lycopersicum L.)

European traditional tomatoes galore: a result of farmers' selection of a few diversity-rich loci

A comprehensive collection of 1254 tomato accessions, corresponding to European traditional and modern varieties, early domesticated varieties, and wild relatives, was analyzed by genotyping by sequencing. A continuous genetic gradient between the traditional and modern varieties was observed. European traditional tomatoes displayed very low genetic diversity, with only 298 polymorphic loci (95% threshold) out of 64 943 total variants. European traditional tomatoes could be classified into several genetic groups. Two main clusters consisting of Spanish and Italian accessions showed higher genetic diversity than the remaining varieties, suggesting that these regions might be independent secondary centers of diversity with a different history. Other varieties seem to be the result of a more recent complex pattern of migrations and hybridizations among the European regions. Several polymorphic loci were associated in a genome-wide association study with fruit morphological traits in the European traditional collection. The corresponding alleles were found to contribute to the distinctive phenotypic characteristic of the genetic varietal groups. The few highly polymorphic loci associated with morphological traits in an otherwise a low-diversity population suggests a history of balancing selection, in which tomato farmers likely maintained the morphological variation by inadvertently applying a high selective pressure within different varietal types.

Atlas of phenotypic, genotypic and geographical diversity present in the European traditional tomato

The Mediterranean basin countries are considered secondary centres of tomato diversification. However, information on phenotypic and allelic variation of local tomato materials is still limited. Here we report on the evaluation of the largest traditional tomato collection, which includes 1499 accessions from Southern Europe. Analyses of 70 traits revealed a broad range of phenotypic variability with different distributions among countries, with the culinary end use within each country being the main driver of tomato diversification. Furthermore, eight main tomato types (phenoclusters) were defined by integrating phenotypic data, country of origin, and end use. Genome-wide association study (GWAS) meta-analyses identified associations in 211 loci, 159 of which were novel. The multidimensional integration of phenoclusters and the GWAS meta-analysis identified the molecular signatures for each traditional tomato type and indicated that signatures originated from differential combinations of loci, which in some cases converged in the same tomato phenotype. Our results provide a roadmap for studying and exploiting this untapped tomato diversity.

Breeding for Resistance to Fusarium Wilt of Tomato: A Review

For over a century, breeders have worked to develop tomato (Solanum lycopersicum) cultivars with resistance to Fusarium wilt (Fol) caused by the soilborne fungus Fusarium oxysporum f. sp. lycopersici. Host resistance is the most effective strategy for the management of this disease. For each of the three Fol races, resistance has been introgressed from wild tomato species, predominately in the form of R genes. The I, I-2, I-3, and I-7 R genes have each been identified, as well as the corresponding Avr effectors in the fungus with the exception of Avr7. The mechanisms by which the R gene protein products recognize these effectors, however, has not been elucidated. Extensive genetic mapping, gene cloning, and genome sequencing efforts support the development of tightly-linked molecular markers, which greatly expedite tomato breeding and the development of elite, Fol resistant cultivars. These resources also provide important tools for pyramiding resistance genes and should support the durability of host resistance.

Cultivated Tomato (Solanum lycopersicum L.) Suffered a Severe Cytoplasmic Bottleneck during Domestication: Implications from Chloroplast Genomes

In various crops, genetic bottlenecks occurring through domestication can limit crop resilience to biotic and abiotic stresses. In the present study, we investigated nucleotide diversity in tomato chloroplast genome through sequencing seven plastomes of cultivated accessions from the Campania region (Southern Italy) and two wild species among the closest (Solanum pimpinellifolium) and most distantly related (S. neorickii) species to cultivated tomatoes. Comparative analyses among the chloroplast genomes sequenced in this work and those available in GenBank allowed evaluating the variability of plastomes and defining phylogenetic relationships. A dramatic reduction in genetic diversity was detected in cultivated tomatoes, nonetheless, a few de novo mutations, which still differentiated the cultivated tomatoes from the closest wild relative S. pimpinellifolium, were detected and are potentially utilizable as diagnostic markers. Phylogenetic analyses confirmed that S. pimpinellifolium is the closest ancestor of all cultivated tomatoes. Local accessions all clustered together and were strictly related with other cultivated tomatoes (S. lycopersicum group). Noteworthy, S. lycopersicum var. cerasiforme resulted in a mixture of both cultivated and wild tomato genotypes since one of the two analyzed accessions clustered with cultivated tomato, whereas the other with S. pimpinellifolium. Overall, our results revealed a very reduced cytoplasmic variability in cultivated tomatoes and suggest the occurrence of a cytoplasmic bottleneck during their domestication.

Differential Regulation of Phytoene Synthase PSY1 During Fruit Carotenogenesis in Cultivated and Wild Tomato Species (Solanum section Lycopersicon)

In plants, carotenoids define fruit pigmentation and are involved in the processes of photo-oxidative stress defense and phytohormone production; a key enzyme responsible for carotene synthesis in fruit is phytoene synthase 1 (PSY1). Tomatoes (Solanum section Lycopersicon) comprise cultivated (Solanum lycopersicum) as well as wild species with different fruit color and are a good model to study carotenogenesis in fleshy fruit. In this study, we identified homologous PSY1 genes in five Solanum section Lycopersicon species, including domesticated red-fruited S. lycopersicum and wild yellow-fruited S. cheesmaniae and green-fruited S. chilense, S. habrochaites and S. pennellii. PSY1 homologs had a highly conserved structure, including key motifs in the active and catalytic sites, suggesting that PSY1 enzymatic function is similar in green-fruited wild tomato species and preserved in red-fruited S. lycopersicum. PSY1 mRNA expression directly correlated with carotenoid content in ripe fruit of the analyzed tomato species, indicating differential transcriptional regulation. Analysis of the PSY1 promoter and 5′-UTR sequence revealed over 30 regulatory elements involved in response to light, abiotic stresses, plant hormones, and parasites, suggesting that the regulation of PSY1 expression may affect the processes of fruit senescence, seed maturation and dormancy, and pathogen resistance. The revealed differences between green-fruited and red-fruited Solanum species in the structure of the PSY1 promoter/5′-UTR, such as the acquisition of ethylene-responsive element by S. lycopersicum, could reflect the effects of domestication on the transcriptional mechanisms regulating PSY1 expression, including induction of carotenogenesis during fruit ripening, which would contribute to red coloration in mature fruit.

ddRAD sequencing-based genotyping for population structure analysis in cultivated tomato provides new insights into the genomic diversity of Mediterranean 'da serbo' type long shelf-life germplasm

Double digest restriction-site associated sequencing (ddRAD-seq) is a flexible and cost-effective strategy for providing in-depth insights into the genetic architecture of germplasm collections. Using this methodology, we investigated the genomic diversity of a panel of 288 diverse tomato (Solanum lycopersicum L.) accessions enriched in 'da serbo' (called 'de penjar' in Spain) long shelf life (LSL) materials (152 accessions) mostly originating from Italy and Spain. The rest of the materials originate from different countries and include landraces for fresh consumption, elite cultivars, heirlooms, and breeding lines. Apart from their LSL trait, 'da serbo' landraces are of remarkable interest for their resilience. We identified 32,799 high-quality SNPs, which were used for model ancestry population structure and non-parametric hierarchical clustering. Six genetic subgroups were revealed, clearly separating most 'da serbo' landraces, but also the Spanish germplasm, suggesting a subdivision of the population based on type and geographical provenance. Linkage disequilibrium (LD) in the collection decayed very rapidly within <5 kb. We then investigated SNPs showing contrasted minor frequency allele (MAF) in 'da serbo' materials, resulting in the identification of high frequencies in this germplasm of several mutations in genes related to stress tolerance and fruit maturation such as CTR1 and JAR1. Finally, a mini-core collection of 58 accessions encompassing most of the diversity was selected for further exploitation of key traits. Our findings suggest the presence of a genetic footprint of the 'da serbo' germplasm selected in the Mediterranean basin. Moreover, we provide novel insights on LSL 'da serbo' germplasm as a promising source of alleles for tolerance to stresses.

Mediterranean Long Shelf-Life Landraces: An Untapped Genetic Resource for Tomato Improvement

The Mediterranean long shelf-life (LSL) tomatoes are a group of landraces with a fruit remaining sound up to 6-12 months after harvest. Most have been selected under semi-arid Mediterranean summer conditions with poor irrigation or rain-fed and thus, are drought tolerant. Besides the convergence in the latter traits, local selection criteria have been very variable, leading to a wide variation in fruit morphology and quality traits. The different soil characteristics and agricultural management techniques across the Mediterranean denote also a wide range of plant adaptive traits to different conditions. Despite the notorious traits for fruit quality and environment adaptation, the LSL landraces have been poorly exploited in tomato breeding programs, which rely basically on wild tomato species. In this review, we describe most of the information currently available for Mediterranean LSL landraces in order to highlight the importance of this genetic resource. We focus on the origin and diversity, the main selective traits, and the determinants of the extended fruit shelf-life and the drought tolerance. Altogether, the Mediterranean LSL landraces are a very valuable heritage to be revalued, since constitutes an alternative source to improve fruit quality and shelf-life in tomato, and to breed for more resilient cultivars under the predicted climate change conditions.

Trichome Independent Resistance against Western Flower Thrips in Tomato

Western flower thrips (WFT) are a major pest on many crops, including tomato. Thrips cause yield losses, not only through feeding damage, but also by the transmission of viruses of which the Tomato Spotted Wilt Virus is the most important one. In cultivated tomato, genetic diversity is extremely low, and all commercial lines are susceptible to WFT. Several wild relatives are WFT resistant and these resistances are based on glandular trichome-derived traits. Introgression of these traits in cultivated lines did not lead to WFT resistant commercial varieties so far. In this study, we investigated WFT resistance in cultivated tomato using a F2 population derived from a cross between a WFT susceptible and a WFT resistant cultivated tomato line. We discovered that this WFT resistance is independent of glandular trichome density or trichome-derived volatile profiles and is associated with three QTLs on chromosomes 4, 5 and 10. Foliar metabolic profiles of F3 families with low and high WFT feeding damage were clearly different. We identified α-tomatine and a phenolic compound as potential defensive compounds. Their causality and interaction need further investigation. Because this study is based on cultivated tomato lines, our findings can directly be used in nowadays breeding programs.

Genetic variants associated with the root system architecture of oilseed rape (Brassica napus L.) under contrasting phosphate supply

Breeding crops with ideal root system architecture for efficient absorption of phosphorus is an important strategy to reduce the use of phosphate fertilizers. To investigate genetic variants leading to changes in root system architecture, 405 oilseed rape cultivars were genotyped with a 60K Brassica Infinium SNP array in low and high P environments. A total of 285 single-nucleotide polymorphisms were associated with root system architecture traits at varying phosphorus levels. Nine single-nucleotide polymorphisms corroborate a previous linkage analysis of root system architecture quantitative trait loci in the BnaTNDH population. One peak single-nucleotide polymorphism region on A3 was associated with all root system architecture traits and co-localized with a quantitative trait locus for primary root length at low phosphorus. Two more single-nucleotide polymorphism peaks on A5 for root dry weight at low phosphorus were detected in both growth systems and co-localized with a quantitative trait locus for the same trait. The candidate genes identified on A3 form a haplotype ‘BnA3Hap’, that will be important for understanding the phosphorus/root system interaction and for the incorporation into Brassica napus breeding programs.

Bayesian QTL mapping using genome-wide SSR markers and segregating population derived from a cross of two commercial F1 hybrids of tomato

Using newly developed euchromatin-derived genomic SSR markers and a flexible Bayesian mapping method, 13 significant agricultural QTLs were identified in a segregating population derived from a four-way cross of tomato. So far, many QTL mapping studies in tomato have been performed for progeny obtained from crosses between two genetically distant parents, e.g., domesticated tomatoes and wild relatives. However, QTL information of quantitative traits related to yield (e.g., flower or fruit number, and total or average weight of fruits) in such intercross populations would be of limited use for breeding commercial tomato cultivars because individuals in the populations have specific genetic backgrounds underlying extremely different phenotypes between the parents such as large fruit in domesticated tomatoes and small fruit in wild relatives, which may not be reflective of the genetic variation in tomato breeding populations. In this study, we constructed F₂ population derived from a cross between two commercial F₁ cultivars in tomato to extract QTL information practical for tomato breeding. This cross corresponded to a four-way cross, because the four parental lines of the two F₁ cultivars were considered to be the founders. We developed 2510 new expressed sequence tag (EST)-based (euchromatin-derived) genomic SSR markers and selected 262 markers from these new SSR markers and publicly available SSR markers to construct a linkage map. QTL analysis for ten agricultural traits of tomato was performed based on the phenotypes and marker genotypes of F₂ plants using a flexible Bayesian method. As results, 13 QTL regions were detected for six traits by the Bayesian method developed in this study.

Solar Radiation-Associated Adaptive SNP Genetic Differentiation in Wild Emmer Wheat, Triticum dicoccoides

Whole-genome scans with large number of genetic markers provide the opportunity to investigate local adaptation in natural populations and identify candidate genes under positive selection. In the present study, adaptation genetic differentiation associated with solar radiation was investigated using 695 polymorphic SNP markers in wild emmer wheat originated in a micro-site at Yehudiyya, Israel. The test involved two solar radiation niches: (1) sun, in-between trees; and (2) shade, under tree canopy, separated apart by a distance of 2-4 m. Analysis of molecular variance showed a small (0.53%) but significant portion of overall variation between the sun and shade micro-niches, indicating a non-ignorable genetic differentiation between sun and shade habitats. Fifty SNP markers showed a medium (0.05 ≤ FST ≤ 0.15) or high genetic differentiation (FST > 0.15). A total of 21 outlier loci under positive selection were identified by using four different FST -outlier testing algorithms. The markers and genome locations under positive selection are consistent with the known patterns of selection. These results suggested that genetic differentiation between sun and shade habitats is substantial, radiation-associated, and therefore ecologically determined. Hence, the results of this study reflected effects of natural selection through solar radiation on EST-related SNP genetic diversity, resulting presumably in different adaptive complexes at a micro-scale divergence. The present work highlights the evolutionary theory and application significance of solar radiation-driven natural selection in wheat improvement.

Genetic diversity and association mapping in the Colombian Central Collection of Solanum tuberosum L. Andigenum group using SNPs markers

The potato (Solanum tuberosum L.) is the fourth most important crop food in the world and Colombia has one of the most important collections of potato germplasm in the world (the Colombian Central Collection-CCC). Little is known about its potential as a source of genetic diversity for molecular breeding programs. In this study, we analyzed 809 Andigenum group accessions from the CCC using 5968 SNPs to determine: 1) the genetic diversity and population structure of the Andigenum germplasm and 2) the usefulness of this collection to map qualitative traits across the potato genome. The genetic structure analysis based on principal components, cluster analyses, and Bayesian inference revealed that the CCC can be subdivided into two main groups associated with their ploidy level: Phureja (diploid) and Andigena (tetraploid). The Andigena population was more genetically diverse but less genetically substructured than the Phureja population (three vs. five subpopulations, respectively). The association mapping analysis of qualitative morphological data using 4666 SNPs showed 23 markers significantly associated with nine morphological traits. The present study showed that the CCC is a highly diverse germplasm collection genetically and phenotypically, useful to implement association mapping in order to identify genes related to traits of interest and to assist future potato genetic breeding programs.

Diverse responses of wild and cultivated tomato to BABA, oligandrin and Oidium neolycopersici infection

Background and Aims

Current strategies for increased crop protection of susceptible tomato plants against pathogen infections include treatment with synthetic chemicals, application of natural pathogen-derived compounds or transfer of resistance genes from wild tomato species within breeding programmes. In this study, a series of 45 genes potentially involved in defence mechanisms was retrieved from the genome sequence of inbred reference tomato cultivar Solanum lycopersicum 'Heinz 1706'. The aim of the study was to analyse expression of these selected genes in wild and cultivated tomato plants contrasting in resistance to the biotrophic pathogen Oidium neolycopersici , the causative agent of powdery mildew. Plants were treated either solely with potential resistance inducers or by inducers together with the pathogen.

Methods

The resistance against O. neolycopersici infection as well as RT-PCR-based analysis of gene expression in response to the oomycete elicitor oligandrin and chemical agent β-aminobutyric acid (BABA) were investigated in the highly susceptible domesticated inbred genotype Solanum lycopersicum 'Amateur' and resistant wild genotype Solanum habrochaites .

Key Results

Differences in basal expression levels of defensins, germins, β-1,3-glucanases, heveins, chitinases, osmotins and PR1 proteins in non-infected and non-elicited plants were observed between the highly resistant and susceptible genotypes. Moreover, these defence genes showed an extensive up-regulation following O. neolycopersici infection in both genotypes. Application of BABA and elicitin induced expression of multiple defence-related transcripts and, through different mechanisms, enhanced resistance against powdery mildew in the susceptible tomato genotype.

Conclusions

The results indicate that non-specific resistance in the resistant genotype S. habrochaites resulted from high basal levels of transcripts with proven roles in defence processes. In the susceptible genotype S. lycopersicum 'Amateur', oligandrin- and BABA-induced resistance involved different signalling pathways, with BABA-treated leaves displaying direct activation of the ethylene-dependent signalling pathway, in contrast to previously reported jasmonic acid-mediated signalling for elicitins.

Exploring New Alleles Involved in Tomato Fruit Quality in an Introgression Line Library of Solanum pimpinellifolium

We have studied a genomic library of introgression lines from the Solanum pimpinellifolium accession TO-937 into the genetic background of the "Moneymaker" cultivar in order to evaluate the accession's breeding potential. Overall, no deleterious phenotypes were observed, and the plants and fruits were phenotypically very similar to those of "Moneymaker," which confirms the feasibility of translating the current results into elite breeding programs. We identified chromosomal regions associated with traits that were both vegetative (plant vigor, trichome density) and fruit-related (morphology, organoleptic quality, color). A trichome-density locus was mapped on chromosome 10 that had not previously been associated with insect resistance, which indicates that the increment of trichomes by itself does not confer resistance. A large number of quantitative trait loci (QTLs) have been identified for fruit weight. Interestingly, fruit weight QTLs on chromosomes 1 and 10 showed a magnitude effect similar to that of QTLs previously defined as important in domestication and diversification. Low variability was observed for fruit-shape-related traits. We were, however, able to identify a QTL for shoulder height, although the effects were quite low, thus demonstrating the suitability of the current population for QTL detection. Regarding organoleptic traits, consistent QTLs were detected for soluble solid content (SSC). Interestingly, QTLs on chromosomes 2 and 9 increased SSC but did not affect fruit weight, making them quite promising for introduction in modern cultivars. Three ILs with introgressions on chromosomes 1, 2, and 10 increased the internal fruit color, making them candidates for increasing the color of modern cultivars. Comparing the QTL detection between this IL population and a recombinant inbred line population from the same cross, we found that QTL stability across generations depended on the trait, as it was very high for fruit weight but low for organoleptic traits. This difference in QTL stability may be due to a predominant additive gene action for QTLs involved in fruit weight, whereas epistatic and genetic background interactions are most likely important for the other traits.

Novel resistance mechanisms of a wild tomato against the glasshouse whitefly

The glasshouse whitefly, Trialeurodes vaporariorum, is an important pest of many crop plants including tomato, Solanum lycopersicum. Many wild tomato species exhibit a higher resistance to whiteflies. Therefore, locating the source of this enhanced resistance and breeding it into commercial tomato species is an important strategy to reduce the impact of pests on crops. Here, we assessed the pest resistance of Lycopersicon pimpinellifolium by comparing oviposition and feeding data from T. vaporariorum on this wild tomato species with data collected from a susceptible commercial tomato, S. lycopersicum var. 'Elegance'. The location of resistance factors was examined by use of electrical penetration graph (EPG) studies on these tomato species. Results show that whiteflies preferentially settled on the commercial tomato more often in 80 % of the replicates when given free choice between the two tomato species and laid significantly fewer eggs on L. pimpinellifolium. Whiteflies exhibited a shorter duration of the second feeding bout, reduced pathway phase probing, longer salivation in the phloem and more non-probing activities in the early stages of the EPG on the wild tomato species compared to the commercial tomato. These findings evidence that a dual mode of resistance is present in this wild tomato against T. vaporariorum: a post-penetration, pre-phloem resistance mechanism and a phloem-located factor, which to the best of our knowledge is the first time that evidence for this has been presented. These findings can be used to inform future breeding strategies to increase the resistance of commercial tomato varieties against this important pest.

Identification of Genetic Differentiation between Waxy and Common Maize by SNP Genotyping

Waxy maize (Zea mays L. var. ceratina) is an important vegetable and economic crop that is thought to have originated from cultivated flint maize and most recently underwent divergence from common maize. In this study, a total of 110 waxy and 110 common maize inbred lines were genotyped with 3072 SNPs to evaluate the genetic diversity, population structure, and linkage disequilibrium decay as well as identify putative loci that are under positive selection. The results revealed abundant genetic diversity in the studied panel and that genetic diversity was much higher in common than in waxy maize germplasms. Principal coordinate analysis and neighbor-joining cluster analysis consistently classified the 220 accessions into two major groups and a mixed group with mixed ancestry. Subpopulation structure in both waxy and common maize sets were associated with the germplasm origin and corresponding heterotic groups. The LD decay distance (1500–2000 kb) in waxy maize was lower than that in common maize. Fourteen candidate loci were identified as under positive selection between waxy and common maize at the 99% confidence level. The information from this study can assist waxy maize breeders by enhancing parental line selection and breeding program design.

Population Structure, Genetic Variation, and Linkage Disequilibrium in Perennial Ryegrass Populations Divergently Selected for Freezing Tolerance

Low temperature is one of the abiotic stresses seriously affecting the growth of perennial ryegrass (Lolium perenne L.), and freezing tolerance is a complex trait of major agronomical importance in northern and central Europe. Understanding the genetic control of freezing tolerance would aid in the development of cultivars of perennial ryegrass with improved adaptation to frost. The plant material investigated in this study was an experimental synthetic population derived from pair-crosses among five European perennial ryegrass genotypes, representing adaptations to a range of climatic conditions across Europe. A total number of 80 individuals (24 of High frost [HF]; 29 of Low frost [LF], and 27 of Unselected [US]) from the second generation of the two divergently selected populations and an unselected (US) control population were genotyped using 278 genome-wide SNPs derived from perennial ryegrass transcriptome sequences. Our studies investigated the genetic diversity among the three experimental populations by analysis of molecular variance and population structure, and determined that the HF and LF populations are very divergent after selection for freezing tolerance, whereas the HF and US populations are more similar. Linkage disequilibrium (LD) decay varied across the seven chromosomes and the conspicuous pattern of LD between the HF and LF population confirmed their divergence in freezing tolerance. Furthermore, two F st outlier methods; finite island model (fdist) by LOSITAN and hierarchical structure model using ARLEQUIN, both detected six loci under directional selection. These outlier loci are most probably linked to genes involved in freezing tolerance, cold adaptation, and abiotic stress. These six candidate loci under directional selection for freezing tolerance might be potential marker resources for breeding perennial ryegrass cultivars with improved freezing tolerance.

Association Analysis for Bacterial Spot Resistance in a Directionally Selected Complex Breeding Population of Tomato

Bacterial spot of tomato is caused by at least four species of Xanthomonas with multiple physiological races. We developed a complex breeding population for simultaneous discovery of marker-trait linkage, validation of existing quantitative trait loci (QTL), and pyramiding of resistance. Six advanced accessions with resistance from distinct sources were crossed in all combinations and their F1 hybrids were intercrossed. Over 1,100 segregating progeny were evaluated in the field following inoculation with X. euvesicatoria race T1 strains. We selected 5% of the most resistant and 5% of the most susceptible progeny for evaluation as plots in two subsequent replicated field trials inoculated with T1 and T3 (X. perforans) strains. The estimated heritability of T1 resistance was 0.32. In order to detect previously reported resistance genes, as well as novel QTL, we explored methods to correct for population structure and analysis based on single markers or haplotypes. Both single-point and haplotype analyses identified strong associations in the genomic regions known to carry Rx-3 (chromosome 5) and Rx-4/Xv3 (chromosome 11). Accounting for kinship and structure generally improved the fit of statistical models. Detection of known loci was improved by adding kinship or a combination of kinship and structure using a Q matrix from model-based clustering. Additional QTL were detected on chromosomes 1, 4, 6, and 7 for T1 resistance and chromosomes 2, 4, and 6 for T3 resistance (P < 0.01). Haplotype analysis improved our ability to trace the origin of positive alleles. These results demonstrate that both known and novel associations can be identified using complex breeding populations that have experienced directional selection.

Exploring a Tomato Landraces Collection for Fruit-Related Traits by the Aid of a High-Throughput Genomic Platform

During its evolution and domestication Solanum lycopersicum has undergone various genetic ‘bottlenecks’ and extreme inbreeding of limited genotypes. In Europe the tomato found a secondary centre for diversification, which resulted in a wide array of fruit shape variation given rise to a range of landraces that have been cultivated for centuries. Landraces represent a reservoir of genetic diversity especially for traits such as abiotic stress resistance and high fruit quality. Information about the variation present among tomato landrace populations is still limited. A collection of 123 genotypes from different geographical areas was established with the aim of capturing a wide diversity. Eighteen morphological traits were evaluated, mainly related to the fruit. About 45% of morphological variation was attributed to fruit shape, as estimated by the principal component analysis, and the dendrogram of relatedness divided the population in subgroups mainly on the basis of fruit weight and locule number. Genotyping was carried out using the tomato array platform SolCAP able to interrogate 7,720 SNPs. In the whole collection 87.1% markers were polymorphic but they decreased to 44–54% when considering groups of genotypes with different origin. The neighbour-joining tree analysis clustered the 123 genotypes into two main branches. The STRUCTURE analysis with K = 3 also divided the population on the basis of fruit size. A genomic-wide association strategy revealed 36 novel markers associated to the variation of 15 traits. The markers were mapped on the tomato chromosomes together with 98 candidate genes for the traits analyzed. Six regions were evidenced in which candidate genes co-localized with 19 associated SNPs. In addition, 17 associated SNPs were localized in genomic regions lacking candidate genes. The identification of these markers demonstrated that novel variability was captured in our germoplasm collection. They might also provide a viable indirect selection tool in future practical breeding programs.

Exploring a Tomato Landraces Collection for Fruit-Related Traits by the Aid of a High-Throughput Genomic Platform

During its evolution and domestication Solanum lycopersicum has undergone various genetic 'bottlenecks' and extreme inbreeding of limited genotypes. In Europe the tomato found a secondary centre for diversification, which resulted in a wide array of fruit shape variation given rise to a range of landraces that have been cultivated for centuries. Landraces represent a reservoir of genetic diversity especially for traits such as abiotic stress resistance and high fruit quality. Information about the variation present among tomato landrace populations is still limited. A collection of 123 genotypes from different geographical areas was established with the aim of capturing a wide diversity. Eighteen morphological traits were evaluated, mainly related to the fruit. About 45% of morphological variation was attributed to fruit shape, as estimated by the principal component analysis, and the dendrogram of relatedness divided the population in subgroups mainly on the basis of fruit weight and locule number. Genotyping was carried out using the tomato array platform SolCAP able to interrogate 7,720 SNPs. In the whole collection 87.1% markers were polymorphic but they decreased to 44-54% when considering groups of genotypes with different origin. The neighbour-joining tree analysis clustered the 123 genotypes into two main branches. The STRUCTURE analysis with K = 3 also divided the population on the basis of fruit size. A genomic-wide association strategy revealed 36 novel markers associated to the variation of 15 traits. The markers were mapped on the tomato chromosomes together with 98 candidate genes for the traits analyzed. Six regions were evidenced in which candidate genes co-localized with 19 associated SNPs. In addition, 17 associated SNPs were localized in genomic regions lacking candidate genes. The identification of these markers demonstrated that novel variability was captured in our germoplasm collection. They might also provide a viable indirect selection tool in future practical breeding programs.

Genomic variation in tomato, from wild ancestors to contemporary breeding accessions

BACKGROUND

Domestication modifies the genomic variation of species. Quantifying this variation provides insights into the domestication process, facilitates the management of resources used by breeders and germplasm centers, and enables the design of experiments to associate traits with genes. We described and analyzed the genetic diversity of 1,008 tomato accessions including Solanum lycopersicum var. lycopersicum (SLL), S. lycopersicum var. cerasiforme (SLC), and S. pimpinellifolium (SP) that were genotyped using 7,720 SNPs. Additionally, we explored the allelic frequency of six loci affecting fruit weight and shape to infer patterns of selection.

RESULTS

Our results revealed a pattern of variation that strongly supported a two-step domestication process, occasional hybridization in the wild, and differentiation through human selection. These interpretations were consistent with the observed allele frequencies for the six loci affecting fruit weight and shape. Fruit weight was strongly selected in SLC in the Andean region of Ecuador and Northern Peru prior to the domestication of tomato in Mesoamerica. Alleles affecting fruit shape were differentially selected among SLL genetic subgroups. Our results also clarified the biological status of SLC. True SLC was phylogenetically positioned between SP and SLL and its fruit morphology was diverse. SLC and "cherry tomato" are not synonymous terms. The morphologically-based term "cherry tomato" included some SLC, contemporary varieties, as well as many admixtures between SP and SLL. Contemporary SLL showed a moderate increase in nucleotide diversity, when compared with vintage groups.

CONCLUSIONS

This study presents a broad and detailed representation of the genomic variation in tomato. Tomato domestication seems to have followed a two step-process; a first domestication in South America and a second step in Mesoamerica. The distribution of fruit weight and shape alleles supports that domestication of SLC occurred in the Andean region. Our results also clarify the biological status of SLC as true phylogenetic group within tomato. We detect Ecuadorian and Peruvian accessions that may represent a pool of unexplored variation that could be of interest for crop improvement.

Analysis of wild-species introgressions in tomato inbreds uncovers ancestral origins

BACKGROUND

Decades of intensive tomato breeding using wild-species germplasm have resulted in the genomes of domesticated germplasm (Solanum lycopersicum) being intertwined with introgressions from their wild relatives. Comparative analysis of genomes among cultivated tomatoes and wild species that have contributed genetic variation can help identify desirable genes, such as those conferring disease resistance. The ability to identify introgression position, borders, and contents can reveal ancestral origins and facilitate harnessing of wild variation in crop breeding.

RESULTS

Here we present the whole-genome sequences of two tomato inbreds, Gh13 and BTI-87, both carrying the begomovirus resistance locus Ty-3 introgressed from wild tomato species. Introgressions of different sizes on chromosome 6 of Gh13 and BTI-87, both corresponding to the Ty-3 region, were identified as from a source close to the wild species S. chilense. Other introgressions were identified throughout the genomes of the inbreds and showed major differences in the breeding pedigrees of the two lines. Interestingly, additional large introgressions from the close tomato relative S. pimpinellifolium were identified in both lines. Some of the polymorphic regions were attributed to introgressions in the reference Heinz 1706 genome, indicating wild genome sequences in the reference tomato genome.

CONCLUSIONS

The methods developed in this work can be used to delineate genome introgressions, and subsequently contribute to development of molecular markers to aid phenotypic selection, fine mapping and discovery of candidate genes for important phenotypes, and for identification of novel variation for tomato improvement. These universal methods can easily be applied to other crop plants.

Genes involved in floral meristem in tomato exhibit drastically reduced genetic diversity and signature of selection

BACKGROUND

Domestication and selection of crops have notably reshaped fruit morphology. With its large phenotypic diversity, tomato (Solanum lycopersicum) illustrates this evolutive trend. Genes involved in flower meristem development are known to regulate also fruit morphology. To decipher the genetic variation underlying tomato fruit morphology, we assessed the nucleotide diversity and selection footprints of candidate genes involved in flower and fruit development and performed genome-wide association studies.

RESULTS

Thirty candidate genes were selected according to their similarity with genes involved in meristem development or their known causal function in Arabidopsis thaliana. In tomato, these genes and flanking regions were sequenced in a core collection of 96 accessions (including cultivated, cherry-type and wild relative accessions) maximizing the molecular diversity, using the Roche 454 technology. A total amount of 17 Mb was sequenced allowing the discovery of 6,106 single nucleotide polymorphisms (SNPs). The annotation of the 30 gene regions identified 231 exons carrying 517 SNPs. Subsequently, the nucleotide diversity (π) and the neutral evolution of each region were compared against genome-wide values within the collection, using a SNP array carrying 7,667 SNPs mainly distributed in coding sequences.About half of the genes revealed footprints of selection and polymorphisms putatively involved in fruit size variation by showing negative Tajima's D and nucleotide diversity reduction in cultivated tomato compared to its wild relative. Among the candidates, FW2.2 and BAM1 sequences revealed selection footprints within their promoter regions suggesting their potential involvement in their regulation. Two associations co-localized with previously identified loci: LC (locule number) and Ovate (fruit shape).

CONCLUSION

Compared to whole genome genotypic data, a drastic reduction of nucleotide diversity was shown for several candidate genes. Strong selection patterns were identified in 15 candidates highlighting the critical role of meristem maintenance genes as well as the impact of domestication on candidates. The study highlighted a set of polymorphisms putatively important in the evolution of these genes.

Exploring genetic variation in the tomato (Solanum section Lycopersicon) clade by whole-genome sequencing

We explored genetic variation by sequencing a selection of 84 tomato accessions and related wild species representative of the Lycopersicon, Arcanum, Eriopersicon and Neolycopersicon groups, which has yielded a huge amount of precious data on sequence diversity in the tomato clade. Three new reference genomes were reconstructed to support our comparative genome analyses. Comparative sequence alignment revealed group-, species- and accession-specific polymorphisms, explaining characteristic fruit traits and growth habits in the various cultivars. Using gene models from the annotated Heinz 1706 reference genome, we observed differences in the ratio between non-synonymous and synonymous SNPs (dN/dS) in fruit diversification and plant growth genes compared to a random set of genes, indicating positive selection and differences in selection pressure between crop accessions and wild species. In wild species, the number of single-nucleotide polymorphisms (SNPs) exceeds 10 million, i.e. 20-fold higher than found in most of the crop accessions, indicating dramatic genetic erosion of crop and heirloom tomatoes. In addition, the highest levels of heterozygosity were found for allogamous self-incompatible wild species, while facultative and autogamous self-compatible species display a lower heterozygosity level. Using whole-genome SNP information for maximum-likelihood analysis, we achieved complete tree resolution, whereas maximum-likelihood trees based on SNPs from ten fruit and growth genes show incomplete resolution for the crop accessions, partly due to the effect of heterozygous SNPs. Finally, results suggest that phylogenetic relationships are correlated with habitat, indicating the occurrence of geographical races within these groups, which is of practical importance for Solanum genome evolution studies.

The potential of pale flax as a source of useful genetic variation for cultivated flax revealed through molecular diversity and association analyses

Pale flax (Linum bienne Mill.) is the wild progenitor of cultivated flax (Linum usitatissimum L.) and represents the primary gene pool to broaden its genetic base. Here, a collection of 125 pale flax accessions and the Canadian flax core collection of 407 accessions were genotyped using 112 genome-wide simple sequence repeat markers and phenotyped for nine traits with the aim of conducting population structure, molecular diversity and association mapping analyses. The combined population structure analysis identified two well-supported major groups corresponding to pale and cultivated flax. The L. usitatissimum convar. crepitans accessions most closely resembled its wild progenitor, both having dehiscent capsules. The unbiased Nei's genetic distance (0.65) confirmed the strong genetic differentiation between cultivated and pale flax. Similar levels of genetic diversity were observed in both species, albeit 430 (48 %) of pale flax alleles were unique, in agreement with their high genetic differentiation. Significant associations were identified for seven and four traits in pale and cultivated flax, respectively. Favorable alleles with potentially positive effect to improve yield through yield components were identified in pale flax. The allelic frequencies of markers associated with domestication-related traits such as capsular dehiscence indicated directional selection with the most common alleles in pale flax being absent or rare in cultivated flax and vice versa. Our results demonstrated that pale flax is a potential source of novel variation to improve multiple traits in cultivated flax and that association mapping is a suitable approach to screening pale flax germplasm to identify favorable quantitative trait locus alleles.

In silico identification and experimental validation of insertion-deletion polymorphisms in tomato genome

Comparative analysis of the genome sequences of Solanum lycopersicum variety Heinz 1706 and S. pimpinellifolium accession LA 1589 using MUGSY software identified 145 695 insertion-deletion (InDel) polymorphisms. A selected set of 3029 candidate InDels (≥2 bp) across the entire tomato genome were subjected to PCR validation, and 82.4% could be verified. Of 2272 polymorphic InDels between LA 1589 and Heinz 1706, 61.6, 45.2, and 31.6% were polymorphic in 8 accessions of S. pimpinellifolium, 4 accessions of S. lycopersicum var. cerasiforme, and 10 varieties of S. lycopersicum, respectively. Genetic distance was 0.216 in S. pimpinellifolium, 0.202 in S. lycopersicum var. cerasiforme, and 0.108 in S. lycopersicum. The data suggested a reduction of genetic variation from S. pimpinellifolium to S. lycopersicum var. cerasiforme and S. lycopersicum. Cluster analysis showed that the 8 accessions of S. pimpinellifolium were in one group, whereas 4 accessions of S. lycopersicum var. cerasiforme and 10 varieties of S. lycopersicum were in the same group.

Patchwork sequencing of tomato San Marzano and Vesuviano varieties highlights genome-wide variations

Background

Investigation of tomato genetic resources is a crucial issue for better straight evolution and genetic studies as well as tomato breeding strategies. Traditional Vesuviano and San Marzano varieties grown in Campania region (Southern Italy) are famous for their remarkable fruit quality. Owing to their economic and social importance is crucial to understand the genetic basis of their unique traits.

Results

Here, we present the draft genome sequences of tomato Vesuviano and San Marzano genome. A 40x genome coverage was obtained from a hybrid Illumina paired-end reads assembling that combines de novo assembly with iterative mapping to the reference S. lycopersicum genome (SL2.40). Insertions, deletions and SNP variants were carefully measured. When assessed on the basis of the reference annotation, 30% of protein-coding genes are predicted to have variants in both varieties. Copy genes number and gene location were assessed by mRNA transcripts mapping, showing a closer relationship of San Marzano with reference genome. Distinctive variations in key genes and transcription/regulation factors related to fruit quality have been revealed for both cultivars.

Conclusions

The effort performed highlighted varieties relationships and important variants in fruit key processes useful to dissect the path from sequence variant to phenotype.

SNP genotyping reveals genetic diversity between cultivated landraces and contemporary varieties of tomato

Background

The tomato (Solanum lycopersium L.) is the most widely grown vegetable in the world. It was domesticated in Latin America and Italy and Spain are considered secondary centers of diversification. This food crop has experienced severe genetic bottlenecks and modern breeding activities have been characterized by trait introgression from wild species and divergence in different market classes.

Results

With the aim to examine patterns of polymorphism, characterize population structure and identify putative loci under positive selection, we genotyped 214 tomato accessions (which include cultivated landraces, commercial varieties and wild relatives) using a custom-made Illumina SNP-panel. Most of the 175 successfully scored SNP loci were found to be polymorphic. Population structure analysis and estimates of genetic differentiation indicated that landraces constitute distinct sub-populations. Furthermore, contemporary varieties could be separated in groups (processing, fresh and cherry) that are consistent with the recent breeding aimed at market-class specialization. In addition, at the 95% confidence level, we identified 30, 34 and 37 loci under positive selection between landraces and each of the groups of commercial variety (cherry, processing and fresh market, respectively). Their number and genomic locations imply the presence of some extended regions with high genetic variation between landraces and contemporary varieties.

Conclusions

Our work provides knowledge concerning the level and distribution of genetic variation within cultivated tomato landraces and increases our understanding of the genetic subdivision of contemporary varieties. The data indicate that adaptation and selection have led to a genomic signature in cultivated landraces and that the subpopulation structure of contemporary varieties is shaped by directed breeding and largely of recent origin. The genomic characterization presented here is an essential step towards a future exploitation of the available tomato genetic resources in research and breeding programs.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-14-835) contains supplementary material, which is available to authorized users.

Comparative analyses identify the contributions of exotic donors to disease resistance in a barley experimental population

Introgression of novel genetic variation into breeding populations is frequently required to facilitate response to new abiotic or biotic pressure. This is particularly true for the introduction of host pathogen resistance in plant breeding. However, the number and genomic location of loci contributed by donor parents are often unknown, complicating efforts to recover desired agronomic phenotypes. We examined allele frequency differentiation in an experimental barley breeding population subject to introgression and subsequent selection for Fusarium head blight resistance. Allele frequency differentiation between the experimental population and the base population identified three primary genomic regions putatively subject to selection for resistance. All three genomic regions have been previously identified by quantitative trait locus (QTL) and association mapping. Based on the degree of identity-by-state relative to donor parents, putative donors of resistance alleles were also identified. The successful application of comparative population genetic approaches in this barley breeding experiment suggests that the approach could be applied to other breeding populations that have undergone defined breeding and selection histories, with the potential to provide valuable information for genetic improvement.

SNP-revealed genetic diversity in wild emmer wheat correlates with ecological factors

Background

Patterns of genetic diversity between and within natural plant populations and their driving forces are of great interest in evolutionary biology. However, few studies have been performed on the genetic structure and population divergence in wild emmer wheat using a large number of EST-related single nucleotide polymorphism (SNP) markers.

Results

In the present study, twenty-five natural wild emmer wheat populations representing a wide range of ecological conditions in Israel and Turkey were used. Genetic diversity and genetic structure were investigated using over 1,000 SNP markers. A moderate level of genetic diversity was detected due to the biallelic property of SNP markers. Clustering based on Bayesian model showed that grouping pattern is related to the geographical distribution of the wild emmer wheat. However, genetic differentiation between populations was not necessarily dependent on the geographical distances. A total of 33 outlier loci under positive selection were identified using a F_ST-outlier method. Significant correlations between loci and ecogeographical factors were observed.

Conclusions

Natural selection appears to play a major role in generating adaptive structures in wild emmer wheat. SNP markers are appropriate for detecting selectively-channeled adaptive genetic diversity in natural populations of wild emmer wheat. This adaptive genetic diversity is significantly associated with ecological factors.

Genetic diversity revealed by single nucleotide polymorphism markers in a worldwide germplasm collection of durum wheat

Evaluation of genetic diversity and genetic structure in crops has important implications for plant breeding programs and the conservation of genetic resources. Newly developed single nucleotide polymorphism (SNP) markers are effective in detecting genetic diversity. In the present study, a worldwide durum wheat collection consisting of 150 accessions was used. Genetic diversity and genetic structure were investigated using 946 polymorphic SNP markers covering the whole genome of tetraploid wheat. Genetic structure was greatly impacted by multiple factors, such as environmental conditions, breeding methods reflected by release periods of varieties, and gene flows via human activities. A loss of genetic diversity was observed from landraces and old cultivars to the modern cultivars released during periods of the Early Green Revolution, but an increase in cultivars released during the Post Green Revolution. Furthermore, a comparative analysis of genetic diversity among the 10 mega ecogeographical regions indicated that South America, North America, and Europe possessed the richest genetic variability, while the Middle East showed moderate levels of genetic diversity.

High-density SNP genotyping of tomato (Solanum lycopersicum L.) reveals patterns of genetic variation due to breeding

The effects of selection on genome variation were investigated and visualized in tomato using a high-density single nucleotide polymorphism (SNP) array. 7,720 SNPs were genotyped on a collection of 426 tomato accessions (410 inbreds and 16 hybrids) and over 97% of the markers were polymorphic in the entire collection. Principal component analysis (PCA) and pairwise estimates of F(st) supported that the inbred accessions represented seven sub-populations including processing, large-fruited fresh market, large-fruited vintage, cultivated cherry, landrace, wild cherry, and S. pimpinellifolium. Further divisions were found within both the contemporary processing and fresh market sub-populations. These sub-populations showed higher levels of genetic diversity relative to the vintage sub-population. The array provided a large number of polymorphic SNP markers across each sub-population, ranging from 3,159 in the vintage accessions to 6,234 in the cultivated cherry accessions. Visualization of minor allele frequency revealed regions of the genome that distinguished three representative sub-populations of cultivated tomato (processing, fresh market, and vintage), particularly on chromosomes 2, 4, 5, 6, and 11. The PCA loadings and F(st) outlier analysis between these three sub-populations identified a large number of candidate loci under positive selection on chromosomes 4, 5, and 11. The extent of linkage disequilibrium (LD) was examined within each chromosome for these sub-populations. LD decay varied between chromosomes and sub-populations, with large differences reflective of breeding history. For example, on chromosome 11, decay occurred over 0.8 cM for processing accessions and over 19.7 cM for fresh market accessions. The observed SNP variation and LD decay suggest that different patterns of genetic variation in cultivated tomato are due to introgression from wild species and selection for market specialization.