Tomato breeding in the genomics era: insights from a SNP array
© Víquez-Zamora et al.; licensee BioMed Central Ltd. 2013
Received: 28 February 2013
Accepted: 20 May 2013
Published: 27 May 2013
The major bottle neck in genetic and linkage studies in tomato has been the lack of a sufficient number of molecular markers. This has radically changed with the application of next generation sequencing and high throughput genotyping. A set of 6000 SNPs was identified and 5528 of them were used to evaluate tomato germplasm at the level of species, varieties and segregating populations.
From the 5528 SNPs, 1980 originated from 454-sequencing, 3495 from Illumina Solexa sequencing and 53 were additional known markers. Genotyping different tomato samples allowed the evaluation of the level of heterozygosity and introgressions among commercial varieties. Cherry tomatoes were especially different from round/beefs in chromosomes 4, 5 and 12. We were able to identify a set of 750 unique markers distinguishing S. lycopersicum ‘Moneymaker’ from all its distantly related wild relatives. Clustering and neighbour joining analysis among varieties and species showed expected grouping patterns, with S. pimpinellifolium as the most closely related to commercial tomatoesearlier results.
Our results show that a SNP search in only a few breeding lines already provides generally applicable markers in tomato and its wild relatives. It also shows that the Illumina bead array generated data are highly reproducible. Our SNPs can roughly be divided in two categories: SNPs of which both forms are present in the wild relatives and in domesticated tomatoes (originating from common ancestors) and SNPs unique for the domesticated tomato (originating from after the domestication event). The SNPs can be used for genotyping, identification of varieties, comparison of genetic and physical linkage maps and to confirm (phylogenetic) relations. In the SNPs used for the array there is hardly any overlap with the SolCAP array and it is strongly recommended to combine both SNP sets and to select a core collection of robust SNPs completely covering the entire tomato genome.
KeywordsSingle Nucleotide Polymorphisms (SNP) Custom made infinium array Tomato wild relatives
Landraces and wild relatives constitute a vast genetic resource that can be tapped to introduce novel traits into tomato breeding programmes . During the last decades, the focus has mainly been on the introduction of disease resistance genes. But, within the breeding efforts, the lack of sufficient molecular markers in tomato has been a bottle neck in genetic and linkage studies. Although all known marker systems have been applied in tomato, most of them fall short in the genomics area mostly because they are too laborious and too low throughput . These shortcomings are now being overcome by next generation sequencing projects and Single Nucleotide Polymorphisms (SNPs) identification . The importance of SNPs as bi-allelic molecular markers is now widely recognized and their use is rapidly increasing [4, 5] , since they have the advantage of being locus specific markers that can be scored co-dominantly in a flexible way. Technology has been developed for scoring single SNPs in thousands of different samples, all the way up to scoring millions of SNPs in a single sample .
Currently, the most widely used systems for high throughput SNP genotyping are the Illumina GoldenGate™, Infinium™ arrays and the KBioscience Competitive Allele‒Specific PCR genotyping system (KASPar: http://www.kbioscience.co.uk) [7–9]. The evolution of genotyping technologies has resulted in unprecedented possibilities for evaluating germplasm collections, characterizing populations, and finding markers linked to specific alleles of important genes. SNPs are also markers of choice for studying evolutionary processes . Characterization of a large set of tomato varieties with a large number of markers can show the impact of breeding on the molecular level and the extent to which these markers are useful for variety identification [11, 12].
A whole genome tomato genotyping array (custom made) using the Illumina® Infinium Beadarray technology  (http://www.illumina.com) was constructed to generate a multiplexing platform [7, 14] to analyse tomato germplasm. A set of 5528 SNPs was used to evaluate more than a thousand tomato samples. This enabled us to compare data at the level of species, varieties and segregating populations. Within the Solanaceae Coordinated Agricultural Project (SolCAP: http://solcap.msu.edu/), in 2012 Sim et al. also developed a genotyping array . However, they focused on different applications and, as we found out, with almost 100% different markers. We were interested in the question to what extent our SNP collection, which is based on a limited number of genotypes, can be applied for variety identification, phylogenetic analysis, genetic mapping, evaluation of introgressions and germplasm identification.
SNPs evaluation and distribution
Validated SNPs and their distribution over the chromosomes
454-seq on cDNA from breeding lines
Illumina Solexa on gDNA from breeding lines
Illumina Solexa on gDNA from introgression free varieties
Markers from previous analysis
Total SNPs per chromosome
However, approximately 10% of the markers could not be reliably scored mainly because of wrong automatic clustering by the GenomeStudio software. Closely linked markers in segregating populations can be used to find the correct score and the reasons for the mistakes in the automatic clustering (Additional file 1). Six percent of the SNPs resulted in NCs. These markers were removed resulting in 4072 SNPs for further analysis. Forty eight percent of the monomorphic markers within S. lycopersicum still were useful because they were polymorphic within tomato wild relatives or between wild relatives and cultivated tomato.
Constructing genetic maps
Variation among varieties
Among the cherry varieties, Gardeners Delight did not have the cherry specific Chromosome 5, and this variety has somewhat larger fruits than what we considered as cherry. In total four round tomatoes clearly had the cherry specific Chromosome 5 (including R100), but after close inspection these were catalogued as deviating from round and more plum types.
Identification of introgressions
SNPs for interspecific crosses in tomato
SNPs among tomato species
Accessions from S. cheesmaniae and S. galapagense clustered together (Figure 8). In our study only approximately 30 SNPs (from 5528) were found between the accessions of S. cheesmaniae and S. galapagense in spite of clear phenotypical differences in leaf structure and trichomes .
Quality of SNP data
The high reproducibility of the results for the 12 Heinz samples shows the robustness of the data obtained with the Infinium array. This was also evident from the comparison between cv Moneymaker and cv Moneyberg where the only differences were a few NCs. Although the data were of high quality, individual SNP calls can be wrong. Wrong calls can be recognized in dense genetic linkage maps of a species from which the sequence is known. We observed errors in 10% of the SNPs, when using the standard settings of the Illumina Genome Studio software. Such errors can be corrected manually or 10% of the SNPs can be deleted . Since the amount of data is vast, enough data remained after deleting 10% of the SNPs. Reasons for errors can be DNA quality, presence of outliers (Additional file 1) within the germplasm and, in few cases, double signalling due to duplications in the genome.
General applicability of the SNPs
Even though the SNPs were looked for in a limited number (4) of breeding lines of S. lycopersicum in combination with four introgression free varieties, they were polymorphic enough in the Solanum sect. Lycopersicon germplasm to discriminate varieties and species as well as to confirm phenetic relations. This implies that many of the SNPs originated from the time before domestication [21–23].
The S. lycopersicum specific markers must have evolved after this species separated from the others. These markers will be polymorphic in any interspecific cross (Figure 7). A relatively cheap, SNP array with a limited number (as few as 20 per chromosome) of well distributed markers will be an excellent tool for a first fast characterization of any new interspecific mapping population involving S. lycopersicum. Based on our results such an array can be easily developed.
Furthermore, it is interesting to note that so many SNPs were found among the four breeding lines. This result was quite unexpected as S. lycopersicum is considered as a species with little genetic variation .
The use of SNPs to improve the tomato genome sequence
One application of the SNP array was to compare genetic with physical positions when working with mapping populations. On the genetic linkage maps most of the markers were in the expected order; identical to the order in the assembled tomato sequence . This confirmed the accurateness of the assembled Tomato Genome.
Some unassigned markers could be mapped to specific chromosomal positions in one or more of the linkage maps that we produced (results not shown). Comparison of genetic linkage maps and the physical linkage map also pointed out a misassembly on the long arm of chromosome 12 (between 48.8 Mbp and 61.7 Mbp; Additional file 6) in version 2.4 of the tomato genome. Also the data published by Sim et al.  suggest a disruption of marker order in the same region (see their Figure 3), but the conclusion that this might be due to a misassembly was not drawn. Markers should be used to genetically validate and further improve de novo genome assemblies.
Several DNA profiling techniques have been used for variety identification . For tomato, one of the most extensive studies was done by Bredemeijer et al. in 2002  using simple sequence repeats (SSRs). They showed that 90% of the more than 500 varieties that were genotyped had a unique SSR profile using 20 markers (on average this is less than 2 markers per chromosome and one chromosome was even without markers). The SNP array covered between 150 and 900 markers per chromosome and all varieties could be distinguished, except the varieties Moneyberg and Moneymaker. That these two showed identical profiles means that they are highly related, if not identical. Both have been registered by the International Union for the Protection of New Varieties of Plants in the National Listing in Great Britain (UPOV : http://www.upov.int), so phenotypic differences must have been seen. Under the UPOV act of 1991 such varieties would likely be considered as essential derived varieties . The SNP markers developed in our study will be very useful for establishing whether varieties are essentially derived from other varieties using the protocol developed for lettuce .
The trend to exploit genes from tomato wild relatives for specific traits enlarges the variation in cultivated tomato and the differences among varieties . Such introgressions can easily be detected using the SNP array as we have shown for the Mi1.2 and TMV gene. When gene-specific (or closely linked) SNP markers are used, genotyping may substitute phenotypic assays even in variety registration as was demonstrated by Arens et al. in 2010 . The markers also allowed us to determine the level of heterozygous markers in present day varieties, which varied between zero and almost 45%. It is interesting to see that the highest numbers are found for some of the plum/cherry tomato. This is most likely because they are hybrids between round and cherry tomatoes and the 955 cherry specific SNPs will contribute to a large number of heterozygous markers (Figures 4 and 5). The high throughput SNP marker determination can be carried out at relatively low cost and is less laborious than other methods used. Therefore it is likely that SNP markers will be the markers of choice for variety identification and registration in future. However it may be anticipated that the SNP arrays will soon be replaced by complete sequencing of varieties.
Differences between round/beef and cherry tomatoes
Many of the polymorphisms located on chromosomes 4, 5 and 12 were between round/beef and cherry tomatoes. This suggests that regions on these chromosomes are essential to get the full cherry tomato phenotype and that there is selection for these regions in breeding programs for cherry tomatoes. The fact that whole chromosomes (4, 5 and 12) look to be involved is possibly due to suppression of recombination in the large pericentromeric regions [29, 30]. This is not the case on Chromosome 1 where the cherry region is a hotspot of recombination as shown in a RIL population of S. lycopersicum and S. pimpinellifolium under study (unpublished observations by the authors
Cherry type tomatoes have more SNPs in common with S. pimpinellifolium accessions than the round/beef varieties indicating that cherry tomatoes are closer to this wild relative than round and beef commercial lines. The varieties chosen for SNP selection might have been the reason that so many cherry specific markers were found. The SolCAP array also revealed different patterns of genetic variation particularly for chromosomes 2, 4, 5, 6 and 11. For chromosome 4 and 5 this is probably also due to the cherry round differences we observed. In general, relatively little is known about genomic regions distinguishing cultivated tomato gene pools .
Some regions are known to contain genes/QTLs that are related to differences between cherry and round. For instance a QTL for fruit weight and soluble solids content, is found on chromosome 2, QTLs for yield, brix, fruit weight, fruit shape, colour and epidermal reticulation have been mapped on chromosome 4 . Chromosome 5 is known to harbour QTLs for fruit colour and QTLs for viscosity traits related to total red yield and pH in chromosome 12 are known [33, 34].
SNPs in Solanum sect. Lycopersicon
Our SNP based phenetic trees were comparable to the ones made by Bretó et al. in 1993 using isozymes , Palmer & Zamir in 1982 and Spooner et al. in 1993 with chloroplast DNA , McClean & Hanson in 1986 with mitochondrial DNA , Miller & Tanksley  with genomic DNA, Marshall et al. in 2001  with internal transcribed spacer (ITS) region of nuclear ribosomal DNA sequences and, also Alvarez et al. in 2001 with microsatellite markers . Peralta et al.  performed the most extensive taxonomic study of tomato and its wild relatives and our results confirm their findings.
In our analysis we found S. pimpinellifolium as the closest wild relative to S. lycopersicum, which is similar to observations made by Grandillo et al.  and The Tomato Genome Consortium in 2012. The cherry tomato is considered either as a domesticated group or as an admixture of S. pimpinellifolium and S. lycopersicum. S. cheesmaniae and S. galapagense are also very closely related to the domesticated tomato. Introgressions in the cultivated germplasm can affect the similarity weight in the relationships between S. pimpinellifolium, S. galapagense and S. cheesmaniae on one hand and S. lycopersicum hybrids on the other hand. For phylogenetic studies it is important to define the initial germplasm and its characteristics. In the case of S. habrochaites and S. pennellii the increased number of NCs decreased the resolution.
Prospects of SNP data in tomato
For our custom made array, the SNP selection was based on commercial breeding lines. Sim et al. [15, 31] developed a large SNP genotyping array using commercial varieties. To evaluate if the same SNPs were present, the precise SNP positions from both arrays were compared (allowing a window of ± 3 base pairs). Only 98 SNPs, less than 2% of our SNPs were found in the exact same position or within the allowed window. This means that there is still a large number of SNPs to be discovered in tomato. For further comparisons among the two arrays we made the SNPs including the flanking sequences available at: http://www.plantbreeding.wur.nl/Publications/SNP/4072SNP-Sequences.xlsx.
Our results show that an SNP search in only a few breeding lines permitted the development of markers generally applicable in tomato and its wild relatives and furthermore that the Illumina bead array generated highly reproducible data. Our SNPs can be roughly divided in two categories: SNPs of which both forms are present in the wild relatives and in domesticated tomatoes and SNPs unique for the domesticated tomato. The SNPs can be used for genotyping, identification of varieties, comparison of genetic and physical linkage maps and to confirm phylogenetic relations. There is hardly any overlap with the SolCAP array and we suggest to combine both SNP sets and to select a core collection of robust SNPs completely covering the tomato genome for the development of future arrays.
Tomato germplasm was obtained from the collection of Wageningen UR Plant Breeding, The Netherlands: the Tomato Genetics Resource Center (TGRC) at University of California, Davis; the Centre for Genetic Resources (CGN), The Netherlands; and from the breeding companies Monsanto, RijkZwaan, Takii, Vilmorin & Cie (VCo), ENZA and Syngenta. The evaluated material included hybrid varieties of the project within the Centre of Biosystems Genomics (CBSG: http://www.cbsg.nl). Based on QTL model predictions, four breeding lines were chosen to obtain a large diversity in taste related characteristics . A half diallel was made with the four breeding lines resulting in six segregating populations. The parents were C74 (cherry, orange), C85 (cherry, red), R75 (round, yellow), and R104 (round, red). Further material included landraces, hybrids, commercial varieties, accessions of tomato wild relatives and mapping populations (Additional file 2: Table S1). The genotyping results of the varieties with the used SNPs can be found in Additional file 3.
DNA and RNA extraction
Genomic DNA from young leaflets was extracted following a CTAB based protocol [44, 45] adjusted for high throughput isolation. Two young leaflets were ground with a Retsch 300 mm shaker (Retsch BV, Ochten, The Netherlands) using 1 ml micronic tubes (Micronic BV, Lelystad, The Netherlands). The DNA pellets were washed in 76% EtOH with 10mM NH4Ac before re-suspending the DNA in TE buffer.
Total RNA was isolated using TRIzol reagent  according to the manufacturer’s instructions (Roche, Switzerland) and finally treated with DNaseI (Invitrogen).
SNPs identified through Roche/454-sequencing
Total RNA was isolated from the four chosen breeding lines (C74, C85, R75 and R104), and at Vertis Biotechnologie AG (Freising, Germany: http://www.vertis-biotech.com/) cDNA was made. The 454 Sequencing gave 1.3 ×106 reads of a median length of 400 base pairs. The reads were aligned to the tomato genome (v2.10) and SNPs were called using QualitySNPng after being adapted for large numbers of reads . After Tomato v2.30 was available the SNP positions were renamed based on this version.
SNPs identified through illumina/solexa-sequencing
A potential risk with the four breeding lines was that primarily interspecific SNPs would be found due to introgressed regions originating from tomato wild relatives (Additional file 2: Table S3). To include additional intraspecific (S. lycopersicum) variation four introgression free varieties were also included in the Illumina/Solexa sequencing. To reduce the complexity, genomic DNA (gDNA) of the eight different samples (C74, C85, R75, R104 and the introgression free varieties, Ailsa Craig/round, Rutgers/beef and Gardeners Delight/cherry plus the reference line Heinz/round) was digested with restriction enzyme MboI (four cutter) and the 400–600 bp fraction was cut out of a 1.5% agarose gel and purified. Theoretically, this should result in a coverage of at least 23x per fragment. After Illumina sequencing 15 × 106 fragments were blasted against the Heinz v2.10 contigs and compared. The Illumina reads of 72 basepairs were aligned with the software tool Bowtie (>95% similarity) . After alignment SNPs were called with VarScan (variant detection in massively parallel sequencing data) . All SNPs with a minimal coverage of three in a genotype were listed in Excell. A SNP was called when it was present in at least six reads in one genotype and six reads in another genotype.
Allocation of SNPs
Putative SNPs and their flanking regions were blasted against the then available contig sequences of tomato (Tomato WGS contigs v2.10) in order to choose SNPs as dispersed over the genome, when possible at least one SNP per contig. Later the availability of the tomato genome sequence (Tomato WGS chromosomes v2.3) allowed us to assign the SNPs to their physical location. A total of 6000 SNPs with two times 50 bp flanking sequences of Heinz were used for designing the oligo’s for the Illumina beadarray . After the oligo’s were synthesized, ~8% of them did not comply to the quality standards set by Illumina and were discarded leaving 5528 SNP markers per array.
Illumina® infinium bead array analysis
Solanum sp. DNA samples with a concentration of 50 ng/μl were sent to ServiceXS, Leiden, The Netherlands, where 4μl was processed according to the Infinium HD Ultra Assay protocol  and used for hybridization onto the BeadChip .
Genotyping data processing
All the SNPs were named after their position on the SL2.30 version of the tomato genome sequence published online by the International Tomato Genome Sequencing Project (http://solgenomics.net/). This version contains approximately 85% of the tomato genome sequence. The lacking sequences are mostly highly repetitive or heterochromatic regions .
The Genotyping Module 1.9.4 of the Illumina’s software GenomeStudio® V2011.1 software package was used to analyse the genotyping results under default settings. The software assigned allele calls (‘GeneCall’) according to the intensity signals obtained, resulting in a [AA], [BB], [AB] or a non-call for each SNP. Advanced assembling within each correspondent analysis was performed and manual inspection and adjustment were performed in order to optimize call rates in the case of questionable SNPs. In particular those cases, and based on the knowledge on segregation patterns within the material, clustering errors were identified and amended .
Before further analysis, markers that were more than 98% monomorphic, were removed, as well as markers with more than 25% heterozygosity in accessions or breeding lines. Finally, also markers with a large number of NCs were removed. For this two thresholds for the percentage NCs were used: more than 20% NCs among the commercial hybrids and/or more than 50% among wild relatives.
When specific populations were evaluated, synchronization of parental lines together with the corresponding offspring was performed. This means that, for each analysis alleles were sorted according to the parent lines and replaced by a specific allele designation (A or B) for each parent.
For cluster analysis the genotype calls were converted into numerical values: [AA]=1, [AB]=2, [BB]=3. Cluster analyses were done using the Jukes-Cantor similarity measure with 1000 bootstraps. Neighbour joining analysis using the Manhattan similarity measure with an out-group rooting and 1000 bootstraps was performed using the statistical package PAST version 2.12 . The BioNJ analysis was carried out using SplitsTree version 4.6 with 1000 bootstraps.
Data visualization heat maps were made in GeneMaths XT 2.12 (Applied Maths). Linkage maps were constructed using JoinMap® version 4.1 (Kyazma©) . The default calculation parameters were adjusted to cope with the large number of markers. In the similarity thresholds the option ‘show individual pairs with a similarity larger than’ was decreased from 0.95 to 0.7. Recombination frequency was used as a grouping parameter and the linkage parameters were set to take all LOD values from 0 to 100. The ‘Show strong linkages with a rec. freq. larger/smaller than’ were set to 0.5/0. The number of maximum linkages to show per locus was set to 0. As algorithm we used the ML (Maximum Likelihood) mapping option, and within the map building, the spatial sampling thresholds were set one to 0.1 the first and the rest to 0. The ‘Number of map optimization rounds per sample’ was fixed to 1. Thereafter, linkage groups were compared with chromosomal distribution in the physical maps using MapChart 2.2 .
This project was carried out within the research programme of the Centre of BioSystems Genomics (CBSG) which is part of the Netherlands Genomics Initiative / Netherlands Organization for Scientific Research.
We would like to acknowledge Fien Meijer-Dekens for collecting and maintaining all the tomato samples at WUR-Plant Breeding. The department of Bioinformatics of Wageningen UR for their support in the development of the SNP markers and in allocating them into the different versions of the Tomato Genome Sequence.
- Miller JC, Tanksley SD: RFLP analysis of phylogenetic relationships and genetic variation in the genus Lycopersicon. Theor Appl Genet. 1990, 80: 437-448.PubMed
- Agarwal M, Shrivastava N, Padh H: Advances in molecular marker techniques and their applications in plant sciences. Plant Cell Rep. 2008, 27: 617-631.View ArticlePubMed
- Tang J, Vosman B, Voorrips RE, van der Linden CG, Leunissen JAM: QualitySNP: a pipeline for detecting single nucleotide polymorphisms and insertions/deletions in EST data from diploid and polyploid species. BMC Bioinforma. 2006, 7: 438-View Article
- Labate JA, Baldo AM: Tomato SNP discovery by EST mining and resequencing. Mol Breed. 2005, 16: 343-349.View Article
- Anithakumari AM, Tang J, Van Eck HJ, Visser RGF, Leunissen JAM, Vosman B, van der Linden CG: A pipeline for high throughput detection and mapping of SNPs from EST databases. Mol Breed. 2010, 26 (1): 65-75.PubMed CentralView ArticlePubMed
- Shirasawa K, Isobe S, Hirakawa H, Asamizu E, Fukuoka H, Just D, Rothan C, Sasamoto S, Fujishiro T, Kishida Y, Kohara M, Tsuruoka H, Wada T, Nakamura Y, Sato S, Tabata S: SNP discovery and linkage map construction in cultivated tomato. DNA Res. 2010, 17 (6): 381-391.PubMed CentralView ArticlePubMed
- Gunderson KL, Steemers FJ, Ren H, Ng P, Zhou L, Tsan C, Chang W, Bullis D, Musmacker J, King C, Lebruska LL, Barker D, Oliphant A, Kuhn KM, Shen R: Whole-genome genotyping. Methods Enzymol. 2006, 410: 359-376.View ArticlePubMed
- Steemers FJ, Gunderson KL: Whole genome genotyping technologies on theBeadArray™ platform. Biotechnol J. 2007, 2: 41-49.View ArticlePubMed
- Appleby N, Edwards D, Batley J: New technologies for ultra-high throughput genotyping in plants. Methods Mol Biol. 2009, 513: 19-39.View ArticlePubMed
- Van Tienderen PH, De Haan AA, van der Linden CG, Vosman B: Biodiversity assessment using markers for ecologically important traits. Trends Ecol Evol. 2002, 17 (12): 577-582.View Article
- Bredemeijer GMM, Arens P, Wouters D, Visser D, Vosman B: The use of semi-automated fluorescent microsatellite analysis for tomato cultivar identification. Theor Appl Genet. 1998, 97: 584-590.View Article
- Bredemeijer GMM, Cooke RJ, Ganal MW, Peeters R, Isaac P, Noordijk Y, Rendell S, Jackson J, Roder MS, Wendehake K, Dijcks M, Amelaine M, Wickaert V, Bertrand L, Vosman B: Construction and testing of a microsatellite database containing more than 500 tomato varieties. Theor Appl Genet. 2002, 105: 1019-1026.View ArticlePubMed
- Illumina®: Infinium® HD Assay: Ultra Protocol Guide. 2009, California, USA: ©Illumina, Inc, 1-224. Catalog #WG-901-4007
- Gunderson KL: Whole-genome genotyping on bead arrays. Methods Mol Biol. 2009, 529: 197-213.View ArticlePubMed
- Sim S-C, Durstewitz G, Plieske J, Wieseke R, Ganal MW, Van Deynze A, Hamilton JP, Buell CR, Causse M, Wijeratne S, Francis DM: Development of a Large SNP Genotyping Array and Generation of High-Density Genetic Maps in Tomato. PLoS One. 2012, 7 (7): 1-11. e40563View Article
- The Tomato Genome Consortium: The tomato genome sequence provides insights into fleshy fruit evolution. Nature. 2012, 485: 635-641.View Article
- Labate JA, Robertson LD: Evidence of cryptic introgression in tomato (Solanum lycopersicum L.) based on wild tomato species alleles. BMC Plant Biol. 2012, 12: 133-PubMed CentralView ArticlePubMed
- Rodriguez F, Wu F, Ané C, Tanksley S, Spooner DM: Do potatoes and tomatoes have a single evolutionary history, and what proportion of the genome supports this history?. BMC Evol Biol. 2009, 9 (1): 191-PubMed CentralView ArticlePubMed
- Gascuel O: BIONJ: an improved version of the NJ algorithm based on a simple model of sequence data. Mol Biol Evol. 1997, 14: 685-695.View ArticlePubMed
- Darwin SC, Knapp S, Peralta IE: Taxonomy of tomatoes in the Galápagos Islands: native and introduced species of Solanum section Lycopersicon (Solanaceae). Syst Biodivers. 2003, 1 (1): 29-53.View Article
- Morin PA, Luikart G, Wayne RK, SNP workshop group: SNPs in ecology, evolution and conservation. Trends Ecol Evol. 2004, 19: 208-216.View Article
- Street SL, Kyes RC, Grant R, Ferguson B: Single nucleotide polymorphisms (SNPs) are highly conserved in rhesus (Macaca mulatta) and cynomolgus (Macaca fascicularis) macaques. BMC Genomics. 2007, 8: 480-PubMed CentralView ArticlePubMed
- Peralta IE, Spooner DM, Knapp S: Taxonomy of wild tomatoes and their relatives (Solanum sections Lycopersicoides, Juglandifolia, Lycopersicon; Solanaceae). Syst Bot Monogr. 2008, 84: 1-186.
- Cooke RJ: Modern methods for cultivar verification and the transgenic plant challenge. Seed Sci Technol. 1999, 27: 669-680.
- UPOV: International convention for the protection of new varieties of plants of December 2, 1961. 1961, Revised at Geneva on November 10, 1972, on October 23, 1978, and on March 19, 1991: http://www.upov.int/en/publications/conventions/1991/act1991.htm
- Vosman B, Visser D, Rouppe J, Van der Voort JR, Smulders MJM, Van Eeuwijk F: Detection of mutants and the establishment of essential derivation among rose varieties using AFLP. Theor Appl Genet. 2004, 109: 1718-1725.View ArticlePubMed
- Van Eeuwijk FA, Law JR: Statistical aspects of essential derivation, with illustrations based on lettuce and barley. Euphytica. 2004, 137: 129-137.View Article
- Arens P, Mansilla C, Deinum D, Cavellini L, Moretti A, Rolland S, van der Schoot H, Calvache D, Ponz F, Collonnier C, Mathis R, Smilde D, Caranta C, Vosman B: Development and evaluation of robust molecular markers linked to disease resistance in tomato for distinctness, uniformity and stability testing. Theor Appl Genet. 2010, 120 (3): 655-664.PubMed CentralView ArticlePubMed
- Sherman JD, Stack SM: Two-dimensional spreads of synaptonemal complexes from solanaceous plants. V. Tomato (Lycopersicon esculentum) karyotype and idiogram. Genome. 1992, 35 (2): 354-359.View Article
- Stack SM, Royer SM, Shearer LA, Chang SB, Giovannoni JJ, Westfall DH, White RA, Anderson LK: Role of fluorescence in situ hybridization in sequencing the tomato genome. Cytogenet Genome Res. 2009, 124: 339-350.View ArticlePubMed
- Sim S-C, Van Deynze A, Stoffel K, Douches DS, Zarka D, Ganal MG, Chetelat RT, Hutton SF, Scott JW, Gardner RG, Panthee DR, Mutschler M, Myers JR, Francis DM: High-Density SNP Genotyping of Tomato (Solanum lycopersicum L.) Reveals Patterns of Genetic Variation Due to Breeding. PLoS One. 2012, 7 (9): 1-18. e45520View Article
- Monforte AJ, Friedman E, Zamir D, Tanksley SD: Comparison of a set of allelic QTL-NILs for chromosome 4 of tomato: Deductions about natural variation and implications for germplasm utilization. Theor Appl Genet. 2010, 102 (4): 572-590.View Article
- Tanksley SD, Grandillo S, Fulton TM, Zamir D, Eshed Y, Petiard V, Lopez J, Beck-Bunn T: Advanced backcross QTL analysis in a cross between an elite processing line of tomato and its wild relative L. pimpinellifolium. Theor Appl Genet. 1996, 92 (2): 213-224.View ArticlePubMed
- Bernacchi D, Beck-Bunn T, Emmatty D, Eshed Y, Inai S, Lopez J, Petiard V, Sayama H, Uhlig J, Zamir D, Tanksley S: Advanced backcross QTL analysis of tomato. II. Evaluation of near-isogenic lines carrying single-donor introgressions for desirable wild QTL-alleles derived from Lycopersicon hirsutum and L. pimpinellifolium. Theor Appl Genet. 1998, 97: 170-180.View Article
- Bretó MP, Asins MJ, Carbonell EA: Genetic variability in Lycopersicon species and their genetic relationships. Theor Appl Genet. 1993, 86: 113-120.View ArticlePubMed
- Palmer JD, Zamir D: Chloroplast DNA evolution and phylogenetic relationships in Lycopersicon. Proc Natl Acad Sci USA. 1982, 79: 5006-5010.PubMed CentralView ArticlePubMed
- Spooner DM, Anderson GJ, Jansen RK: Chloroplast DNA evidence for the interrelationships of tomatoes, potatoes, and pepinos (Solanaceae). Am J Bot. 1993, 80: 676-688.View Article
- McClean PE, Hanson MR: Mitochondrial DNA sequence divergence among Lycopersicon and related Solanum species. Genetics. 1986, 112: 649-667.PubMed CentralPubMed
- Marshall JA, Knapp S, Davey MR, Power JB, Cocking EC, Bennett MD, Cox AV: Molecular systematics of Solanum section Lycopersicum (Lycopersicon) using the nuclear ITS rDNA region. Theor Appl Genet. 2001, 103: 1216-1222.View Article
- Alvarez AE, van de Wiel CCM, Smulders MJM, Vosman B: Use of microsatellites to evaluate genetic diversity and species relationships in the genus Lycopersicon. Theor Appl Genet. 2001, 103: 1283-1292.View Article
- Grandillo S, Chetelat R, Knapp S, Spooner D, Peralta I, Cammareri M, Perez O, Termolino P, Tripodi P, Chiusano ML, Ercolano MR, Frusciante L, Monti L, Pignone D: Solanum sect. Lycopersicon. Wild Crop Relatives: Genomic and Breeding Resources, Vegetables. Edited by: Kole C. 2011, Berlin Heidelberg: Springer-Verlag, 129-215.View Article
- Ranc N, Muños S, Santoni S, Causse M: A clarified position for solanum lycopersicum var. cerasiforme in the evolutionary history of tomatoes (solanaceae). BMC Plant Biol. 2008, 8: 130-PubMed CentralView ArticlePubMed
- Van Berloo R, Van Heusden S, Bovy A, Meijer-Dekens F, Lindhout P, Van Eeuwijk F: Genetic research in a public-private research consortium: prospects for indirect use of Elite breeding germplasm in academic research. Euphytica. 2008, 161: 293-300.View Article
- Steward CN, Via LE: A rapid CTAB DNA isolation technique useful for RAPD fingerprinting and other PCR applications. Biotechniques. 1993, 14: 748-750.
- Kabelka E, Franchino B, Francis DM: Two loci from Lycopersicon hirsutum LA407 confer resistance to strains of Clavibacter michiganensis subsp. michiganensis. Phytopath. 2002, 92: 504-510.View Article
- Simms D, Cizdziel PE, Chomcynski P: TRIZOL: A new reagent for optimal single-step isolation of RNA. Focus. 1993, 15: 99-102.
- Nijveen H, Van Kaauwen M, Esselink D, Hoegen B, Vosman B: QualitySNPng: a user-friendly SNP detection and visualisation tool. Nucleic Acids Res. 2013, -10.1093/nar/gkt333.
- Langmead B, Trapnell C, Pop M, Salzberg SL: Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009, 10 (3): R25-PubMed CentralView ArticlePubMed
- Koboldt DC, Chen K, Wylie T, Larson DE, McLellan MD, Mardis ER, Weinstock GM, Wilson RK, Ding L: VarScan: variant detection in massively parallel sequencing of individual and pooled samples. Bioinformatics. 2009, 25 (17): 2283-2285.PubMed CentralView ArticlePubMed
- Service XS: Project Report: Illumina iSelect Genotyping. 2010, 2333 BZ Leiden, 1-28. Report produced by company Service XS (http://www.servicexs.com)
- Teo YY, Inouye M, Small KS, Gwilliam R, Deloukas P, Kwiatkowski DP, Clark TG: A genotype calling algorithm for the Illumina BeadArrayplatform. Bioinformatics. 2007, 23: 2741-2746.PubMed CentralView ArticlePubMed
- Hammer Ø, Harper DAT, Ryan PD: PAST: Paleontological Statistics software package for education and data analysis. Palaeontol Electron. 2001, 4 (1): 1-9.
- Van Ooijen JW: Multipoint maximum likelihood mapping in a full-sib family of an outbreeding species. Genet Res. 2011, 93: 343-349.View Article
- Voorrips RE: MapChart: Software for the graphical presentation of linkage maps and QTLs. J Hered. 2002, 93 (1): 77-78.View ArticlePubMed
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.