Phenotypic diversity within cultivated tomato (Solanum lycopersicum) is particularly evident for fruit shape and size. Four genes that control tomato fruit shape have been cloned. SUN and OVATE control elongated shape whereas FASCIATED (FAS) and LOCULE NUMBER (LC) control fruit locule number and flat shape. We investigated the distribution of the fruit shape alleles in the tomato germplasm and evaluated their contribution to morphology in a diverse collection of 368 predominantly tomato and tomato var. cerasiforme accessions. Fruits were visually classified into eight shape categories that were supported by objective measurements obtained from image analysis using the Tomato Analyzer software. The allele distribution of SUN, OVATE, LC, and FAS in all accessions was strongly associated with fruit shape classification. We also genotyped 116 representative accessions with additional 25 markers distributed evenly across the genome. Through a model-based clustering we demonstrated that shape categories, germplasm classes, and the shape genes were nonrandomly distributed among five genetic clusters (P < 0.001), implying that selection for fruit shape genes was critical to subpopulation differentiation within cultivated tomato. Our data suggested that the LC, FAS, and SUN mutations arose in the same ancestral population while the OVATE mutation arose in a separate lineage. Furthermore, LC, OVATE, and FAS mutations may have arisen prior to domestication or early during the selection of cultivated tomato whereas the SUN mutation appeared to be a postdomestication event arising in Europe.

The concurrent development of high-throughput genotyping platforms and next generation sequencing (NGS) has increased the number and density of genetic markers, the efficiency of constructing detailed linkage maps, and our ability to overlay recombination and physical maps of the genome. We developed an array for tomato with 8,784 Single Nucleotide Polymorphisms (SNPs) mainly discovered based on NGS-derived transcriptome sequences. Of the SNPs, 7,720 (88%) passed manufacturing quality control and could be scored in tomato germplasm. The array was used to generate high-density linkage maps for three interspecific F₂ populations: EXPEN 2000 (Solanum lycopersicum LA0925 x S. pennellii LA0716, 79 individuals), EXPEN 2012 (S. lycopersicum Moneymaker x S. pennellii LA0716, 160 individuals), and EXPIM 2012 (S. lycopersicum Moneymaker x S. pimpinellifolium LA0121, 183 individuals). The EXPEN 2000-SNP and EXPEN 2012 maps consisted of 3,503 and 3,687 markers representing 1,076 and 1,229 unique map positions (genetic bins), respectively. The EXPEN 2000-SNP map had an average marker bin interval of 1.6 cM, while the EXPEN 2012 map had an average bin interval of 0.9 cM. The EXPIM 2012 map was constructed with 4,491 markers (1,358 bins) and an average bin interval of 0.8 cM. All three linkage maps revealed an uneven distribution of markers across the genome. The dense EXPEN 2012 and EXPIM 2012 maps showed high levels of colinearity across all 12 chromosomes, and also revealed evidence of small inversions between LA0716 and LA0121. Physical positions of 7,666 SNPs were identified relative to the tomato genome sequence. The genetic and physical positions were mostly consistent. Exceptions were observed for chromosomes 3, 10 and 12. Comparing genetic positions relative to physical positions revealed that genomic regions with high recombination rates were consistent with the known distribution of euchromatin across the 12 chromosomes, while very low recombination rates were observed in the heterochromatic regions.

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

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.

We used molecular markers to identify quantitative trait loci (QTL) that contribute to resistance to bacterial canker of tomato caused by Clavibacter michiganensis subsp. michiganensis. Resistance was first identified as a marker-trait association in an inbred backcross (IBC) population derived from crossing Lycopersicon hirsutum accession (LA407) with L. esculentum. Single-marker QTL analysis suggested that at least two loci originating from L. hirsutum LA407, Rcm 2.0 on chromosome 2 and Rcm 5.1 on chromosome 5, contribute to resistance in replicated trials. Two segregating F₂populations were developed by crossing resistant inbred backcross lines (IBLs) to elite L. esculentum lines and used to confirm QTL associations detected in the IBC population. In these populations, realized heritability estimates were higher for selection based on maximal disease than for selection based on disease progression. Realized heritability in the population carrying Rcm 2.0 was 0.63 and 0.14, respectively, for each selection criteria. Realized heritability estimates were 0.85 for selection based on maximal disease and 0.37 for selection based on disease progression in a population carrying Rcm 5.1. The disease response of F₃families selected for resistance suggested that both Rcm 2.0 and Rcm 5.1 confer resistance to bacterial strains in the repetitive sequence-based polymerase chain reaction DNA fingerprint classes A and C. Markers linked to Rcm 2.0 explained up to 56% of the total phenotypic variation for resistance in one population, and markers linked to Rcm 5.1 explained up to 73% of the total phenotypic variation for resistance in a separate population.