HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS		QUICK SEARCH:   [advanced] Author: Keyword(s): Year:  Vol:  Page:

This Article Abstract Full Text (PDF) Supplementary Data Submit a response Alert me when this article is cited Alert me when eLetters are posted Alert me if a correction is posted Services Email this article to a friend Similar articles in this journal Similar articles in ISI Web of Science Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via ISI Web of Science (50) Citing Articles via Google Scholar Google Scholar Articles by Peralta, I. E. Articles by Spooner, D. M. Search for Related Content PubMed Articles by Peralta, I. E. Articles by Spooner, D. M. Agricola Articles by Peralta, I. E. Articles by Spooner, D. M.

Granule-bound starch synthase (GBSSI) gene phylogeny of wild tomatoes (Solanum L. section Lycopersicon [Mill.] Wettst. subsection Lycopersicon)¹

Vegetable Crops Research Unit, USDA, Agricultural Research Service, Department of Horticulture, University of Wisconsin, 1575 Linden Drive, Madison, Wisconsin 53706-1590 USA

Received for publication August 11, 2000. Accepted for publication March 8, 2001.

	ABSTRACT

TOP ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS DISCUSSION LITERATURE CITED

 
Eight wild tomato species are native to western South America and one to the Galapagos Islands. Different classifications of tomatoes have been based on morphological or biological criteria. Our primary goal was to examine the phylogenetic relationships of all nine wild tomato species and closely related outgroups, with a concentration on the most widespread and variable tomato species Solanum peruvianum, using DNA sequences of the structural gene granule-bound starch synthase (GBSSI, or waxy). Results show some concordance with previous morphology-based classifications and new relationships. The ingroup comprised a basal polytomy composed of the self-incompatible green-fruited species S. chilense and the central to southern Peruvian populations of S. peruvianum, S. habrochaites, and S. pennellii. A derived clade contains the northern Peruvian populations of S. peruvianum (also self-incompatible, green-fruited), S. chmielewskii, and S. neorickii (self-compatible, green-fruited), and the self-compatible and red- to orange- to yellow-fruited species S. cheesmaniae, S. lycopersicum, and S. pimpinellifolium. Outgroup relationships are largely concordant with prior chloroplast DNA restriction site phylogenies, support S. juglandifolium and S. ochranthum as the closest outgroup to tomatoes with S. lycopersicoides and S. sitiens as basal to these, and support allogamy, self-incompatibility, and green fruits as primitive in the tomato clade.

Key Words: GBSSI • Lycopersicon • phylogeny • Solanaceae • Solanum sect. Lycopersicon • tomato • waxy gene

	INTRODUCTION

TOP ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS DISCUSSION LITERATURE CITED

 
Wild tomatoes (Solanum L. subsect. Lycopersicon, an autonym of Solanum section Lycopersicon (Mill.) Wettst.) are entirely American in distribution, growing in the western South American Andes from central Ecuador through Peru to northern Chile and in the Galapagos Islands, where the endemic species S. cheesmaniae grows (Fig. 1; see Table 1 for taxonomic authorities). In addition, S. lycopersicum, the wild ancestor of cultivated tomatoes, is more widespread and perhaps more recently distributed into Mexico, Colombia, Bolivia, and other South American countries (Rick and Holle, 1990 ). Wild tomatoes grow in a variety of habitats, from near sea level to >3300 m in elevation (Rick, 1973 ; Taylor, 1986 ). These habitats include the arid Pacific coast to the mesic uplands in the high Andes. All species are diploid (2n = 24; Rick, 1979 ). Breeding systems vary from allogamous self-incompatible to facultative allogamous and self-compatible to autogamous and self-compatible (Rick, 1963, 1979, 1984, 1986 ; Table 1). The self-incompatibility system in tomatoes is gametophytic and controlled by a single, multiallelic S locus (Tanksley and Loaiza-Figueroa, 1985 ).

View larger version (48K): [in this window] [in a new window]   Fig. 1. Distribution of the wild tomato accessions used in this study. Map numbers correspond to generalized map localities in Table 1 . A line drawn at 10° S indicates the northern (N) and southern (S) edges of distribution of "north" and "south" Solanum peruvianum populations that reflect the phylogenetic relationships showed in Fig. 5

The two most recent and complete taxonomic treatments of tomatoes are those of Müller (1940) and Luckwill (1943 ; Fig. 2). Müller (1940) considered Lycopersicon as a distinct genus and divided it into two subgenera, Eulycopersicon C. H. Müll., with two species possessing glabrous and red- to orange-colored fruits, flat seeds, bractless inflorescences, and leaves without pseudostipules, and Eriopersicon C. H. Müll., with four species bearing pubescent green fruits, thick seeds, bracteate inflorescences, and leaves usually with pseudostipules. Three years later, Luckwill (1943) adopted the same supraspecific categories but recognized different infraspecific taxa and five species in the subgenus Eriopersicon (Fig. 2). These treatments have become outdated as the number of species and races collected from South America has increased (Rick, 1971, 1991 ; Holle, Rick, and Hunt, 1978, 1979 ; Taylor, 1986 ).

View larger version (28K): [in this window] [in a new window]   Fig. 2. Comparison of classifications of Solanum L. section Lycopersicon [Mill.] Wettst. subsection Lycopersicon. Species in boldface have red to orange to yellow fruits and those without boldface have green fruits. The numbers in parentheses represent numbers of infraspecific ranks (subspecies, varieties, and forms). The lines connect synonymous taxa

In a recent taxonomic treatment, Child (1990 ; Fig. 2) placed the tomatoes in Solanum subgenus Potatoe (G. Don) D'Arcy, section Lycopersicon, subsection Lycopersicon, and segregated them into three series: Eriopersicon (C. H. Müll.) Child, Lycopersicon, and Neolycopersicon (Correll) Child. He also hypothesized that Solanum subsection Lycopersicoides Child and section Juglandifolium (Rydb.) Child are the closest relatives of section Lycopersicon. Different criteria used in classification, morphology, and crossability, therefore, have led to different numbers of species, subspecies, and varieties and to conflicting hypotheses of interspecific relationships (Fig. 2).

Molecular data provide numerous and independent characters useful for plant phylogenetic inference (Sytsma, 1990 ). Single or low-copy nuclear genes are of great interest in molecular systematics because they are potentially most informative for species-level comparisons. Nuclear gene introns, such as in the genes Phy, Pgi, adh-1, and GapA, evolve faster and have been considered useful for angiosperm phylogenetic studies (Soltis and Soltis, 1998 ). In the perennial relatives of soybean (Glycine) sequences from the histone H3-D intron provided useful characters for assessing relationships at the species level (Doyle, Kanazin, and Shoemaker, 1996 ).

The phylogenetic utility of the structural nuclear gene for granule-bound starch synthase (GBSSI or waxy gene) was first explored by Mason-Gamer and Kellogg (1996) in grasses. The intron–exon regions revealed an interesting structure with potential utility at different taxonomic levels. The GBSSI introns provided more phylogenetically informative characters, and higher levels of variation, than the ITS region in grasses (Mason-Gamer, Weil, and Kellogg, 1998 ). A similar sequence variation was found among Ipomoea species comparing ITS with three exons (9, 10, 11; Fig. 3) and two introns, comprising in total 651 nucleotide sites, of the 3' end of GBSSI (Miller, Raush, and Manos, 1999 ). Nevertheless, ITS sequences showed more variation than GBSSI sequences within the genus Physalis, although only 700 base pairs (bp) of GBSSI were analyzed (Whitson and Manos, 1999 ).

View larger version (13K): [in this window] [in a new window]   Fig. 3. Schematic of the structural potato GBBSI (waxy) gene analyzed in this study. The promoter sequence (TACAAAT), translational start (ATG), stop codon (TAA), and polyadenylation sites (polyA) are indicated. The black box represents an untranslated exon, shaded boxes the 13 translated exons and interconnecting lines introns. The fragment obtained by the polymerase chain reaction comprises 1525 base pairs (from the first to the eighth exon) and is indicated with a line below the schematic. Arrow heads indicate position and direction of primers 5' and 3' that were used to amplify the fragment and forward primers A and B and reverse primers CR and DR that were used for cycle sequencing. The scale represents the total length of complete genomic nucleotide sequence of potato GBBSI (modified from van der Leij et al., 1991 )

The main goal of our research was to examine the phylogenetic relationships of all nine wild tomato species and closely related outgroups, with a concentration on the most widespread and variable tomato species, Solanum peruvianum, using DNA sequences of the structural gene granule-bound starch synthase. A secondary goal was to compare these results to previous systematic treatments in wild tomatoes that use morphology (Müller, 1940 ; Luckwill, 1943 ; Child, 1990 ), crossing relationships (Rick, 1963, 1979, 1986 ), isozyme (Rick and Fobes, 1975a, b ; Rick and Tanksley, 1981 ; Rick, 1983, 1986 ; Bretó, Asins, and Carbonell, 1993 ), or molecular data (Palmer and Zamir, 1982 ; McClean and Hanson, 1986 ; Miller and Tanksley, 1990 ).

	MATERIALS AND METHODS

TOP ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS DISCUSSION LITERATURE CITED

 
Plants
In the ingroup (Solanum section Lycopersicon subsection Lycopersicon) three accessions per species were analyzed for S. cheesmaniae, S. chmielewskii, S. chilense, S. habrochaites, S. lycopersicum var. cerasiforme, S. neorickii, S. pennellii; five accessions for S. pimpinellifolium; and 39 for S. peruvianum. Many more accessions of S. peruvianum were analyzed because this species is the most widespread and has the greatest degree of morphological, molecular, and crossability variation (Rick, 1963, 1979, 1986 ; Miller and Tanksley, 1990 ). One accession of a cultivated tomato, S. lycopersicum CV ‘Mountain Spring,’ was also included. The outgroups were selected according to the cpDNA results of Spooner, Anderson, and Jansen (1993) . Like the ingroup, they are all diploid at 2n = 24 (Table 2). Two accessions per species were included in the outgroup taxa S. juglandifolium, S. lycopersicoides, S. ochranthum, S. sitiens and one accession per species in the outgroup taxa S. bulbocastanum, S. etuberosum, S. jamesii, S. muricatum and S. palustre. In total, 79 accessions were used for the phylogenetic analysis (Table 1). All the tomato (Solanum sect. Lycopersicon subsection Lycopersicon) accessions and outgroups in sect. Lycopersicon subsection Lycopersicoides and sect. Juglandifolium were obtained from the Tomato Genetics Resource Center, Department of Vegetable Crops, University of California, Davis. All potato accessions (Solanum sect. Petota) and members of sect. Etuberosum were obtained from the United States Potato Introduction Station, National Research Support Program-6 (NRSP-6), Sturgeon Bay, Wisconsin, USA (Bamberg et al., 1996 ); one accession of sect. Basarthrum was obtained from Dr. Gregory Anderson, University of Connecticut, Storrs, Connecticut, USA. Dr. Charles Rick, curator of the Tomato Genetics Resource Center, provided advice for the choice of tomato accessions based on geographic distribution, morphology, genetic diversity, and breeding behavior. Vouchers for these accessions were deposited at DAV, MERL, PTIS, and WIS.

Generation of GBSSI sequences
The oligonucleotide primers were developed based on the GBSSI sequences of S. tuberosum obtained from GenBank. Two oligonucleotides were designed as external primers (Fig. 3) for the polymerase chain reaction (PCR), in the first exon region at position 136 (5'-GATGGGCTCCAATCAAGAACTAAT-3') and in the eighth exon at position 1639 (5'-GCCATTCACAATCCCAGTTATGC-3') of the published sequence of GBSSI (van der Leij et al., 1991 ), using the software program Oligo (Rychlik, 1998 ). Within this region (1311 bp) there are seven introns (from second to eighth introns) of lengths 84, 84, 101, 83, 80, 90, and 83 bp, respectively (total of 602 bp; van der Leij et al., 1991 ). This region was amplified using PCR for all accessions. In all cases, a single PCR product was obtained (Fig. 3).

Amplification reactions were conducted in a total volume of 50 µL with a Perkin Elmer programmable thermal controller. Each reaction comprised Promega (Madison, Wisconsin, USA) reaction buffer (1x final concentration), MgCl (2.5 mmol/L), dNTPs (200 µmol/L each), the two external primers (0.4 mmol/L each), bovine serum albumin (0.4 mmol/L), template DNA (50–100 ng/µL), and Promega Taq DNA polymerase (2.5 units/µL).

The amplification reaction was conducted for 4 min at 94°C, then 5 cycles of 30 sec at 94°C, 1 min at 50°C, and 1.5 min at 72°C, then 30 cycles of 30 sec at 94°C, 1 min at 55°C, and 1.5 min at 72°C, with a 10-min extension at 72°C, ending with 4°C hold. To determine whether amplification was successful, samples were run on a 1.0–1.2% agarose minigel stained with ethidium bromide and viewed with ultraviolet light. As a control, a lambda clone of the GBSSI gene from potato, kindly provided by Dr. Richard Visser (Wageningen, The Netherlands), was used. A DNA ladder (Gibco BRL low mass DNA ladder; Life Technologies, Rockville, Maryland, USA) was included to estimate the length and concentration of the amplified products.

Successful reactions were purified for cycle sequencing using cleaning columns (QIAgen Valencia, California, USA). Concentration of amplified cleaned product was checked with a fluorimeter and adjusted to 10–30 ng per 20 µL reaction. Cycle sequencing reactions were carried out using the ABI Prism Dye Terminator Cycle Sequencing Kit Reaction with AmpliTaq DNA Polymerase FS, and with Big Dye (Perkin-Elmer, Applied Biosystems, Foster City, California, USA) according to manufacturer's instructions.

Four internal primers (Fig. 3) were designed based on potato Genbank sequence to match conserved exon regions, beginning at position 174 (GBSSI A: 5'-CAAGATGGCATCCAGAACTGAGA-3'), 622 (GBSSI B: 5'-CACTGCTATAAACGTGGGGTTGA-3'), 1555 (GBSSI CR: 5'-GGCATAGTATGGGCTCACTGTAA-3'), and 1126 (GBSSI DR 5'-GGAATGAGAGCTGTGTGCCAATC-3') of the published sequence of GBSSI (van der Leij et al., 1991 ). Sequences were obtained from both the forward (A, B) and reverse (C_R and D_R) primers.

Cycle sequencing was conducted for 25 cycles, each cycle consisting of 96°C for 30 sec, 46°C for 15 sec, and 60°C for 4 min. Sequencing products were precipitated in 50 µL of 95% ethanol and 2 µL of 3 mol/L sodium acetate, and then vacuum-dried. The samples were resuspended in formamide loading buffer (84% formamide, 4 mmol/L EDTA, 8 mg/mL Blue Dextran) and electrophoresed on a Perkin-Elmer Applied Biosystems 377XL automated DNA sequencer.

Some of the GBSSI sequences showed polymorphisms, particularly in the third intron, which was difficult to align. A similar protocol was used to amplify GBSSI in the problematic accessions, but in this case using Pfu polymerase (Promega, Madison, Wisconsin, USA), which produces less error during the PCR amplification than Taq polymerase. Cleaned PCR products were then cloned into Promega's pGEM-T Easy vector. The vector contain an A 3' terminal thymidine overhang, and T7 and SP6 RNA polymerase promoters flanking multiple cloning regions including the -peptide coding region of the enzyme ß-galactosidase and multiple restriction sites. Because Pfu polymerase produces blunt ends, an A-tailing reaction was performed using Taq DNA polymerase before the ligation. Ligation, transformation, and plating were done according to the protocols provided by the manufacturer. Insertional inactivation of the -peptide allows recombinant clones to be directly identified by color screening on indicator plates. Plasmid preparations were done with alkaline lysis (Wizard plus plasmid purification systems; Promega, Madison, Wisconsin, USA). Purified plasmid DNA was digested with EcoRI, and cut and uncut products were checked on 0.9% agarose gel for presence of inserts. A lambda control of GBSSI was included to check the size of the band. Two molecular ladders were also included (Gibco BRL low and high mass DNA ladders; Life Technologies, Rockville, Maryland, USA), and concentrations of the inserts were estimated by visual comparisons of band intensity. Plasmid minipreps containing the GBSSI fragment were cycle sequenced as described above. Six to 12 clones were sequenced for each individual plant per accession (Table 2).

Sequence alignment
DNA sequences initially were edited using Sequencher software (Gene Codes Corporation, Ann Arbor, Michigan, USA). These initially aligned sequences were then imported into PAUP* version 4.0d64 (Swofford, 1998 ) and were further aligned visually.

Phylogenetic analysis
Parsimony analyses of the DNA sequences were conducted with PAUP* 4.0d64 (Swofford, 1998 ). Solanum muricatum was used to root the trees based on the results of Spooner, Anderson, and Jansen (1993) and Olmstead and Palmer (1997) . The first analyses were performed using Fitch parsimony, which allows free reversibility (change of characters in any direction among the nucleotide states). All the characters were included and were equally weighted, and gaps were considered as missing data. An initial tree was generated by jackknife (36.8% deletion, 10 000 replicates, and fast stepwise addition) (Farris et al., 1996 ). This tree was used as a starting tree for further heuristic searches using tree-bisection-reconnection (TBR) swapping algorithm, with retention of all equally parsimonious trees (MULPARS) with stepwise addition and steepest descent. The amount of homoplasy was evaluated with the consistency index (CI; Kluge and Farris, 1969 ) and the retention index (RI; Farris, 1989 ). Consensus trees were generated from all most parsimonious trees. Bootstrap analysis (Felsenstein, 1985 ), using 1000 replicates, and jackknife, also with 1000 replicates (36.8% deletions; Farris et al., 1996 ), were conducted to estimate the internal relative support for each branch. Also, a decay analysis (Bremer, 1988 ; Donoghue et al., 1992 ) was calculated using an heuristic TBR search. Bootstrap and decay values are placed on the strict consensus tree (Fig. 5).

View larger version (37K): [in this window] [in a new window]   Fig. 5. Strict consensus tree of 15 000 most parsimonious trees (indels treated as characters) based on GBSSI sequences. Percentage of 1000 bootstrap replicates is given above branches. Decay values between 1 and 3 are in parentheses; the other branches in the outgroup are at least 13 steps longer. Distinct GBSSI sequence types from different clones obtained from single individuals of seven accessions are indicated in italic boldface and a number (see Table 2 ). Map localities correspond to Fig. 1 and Table 1 ; G represents accessions from the Galapagos Islands, and S represents unmapped continental accessions. "North per" refers to accessions from the northern range of S. peruvianum, while "south per" refers to accessions from central and southern populations of S. peruvianum, as shown in Fig. 1

	RESULTS

TOP ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS DISCUSSION LITERATURE CITED

 
GBSSI sequence variation within individuals
About 15% of the 79 accessions showed some sites with two, and sometimes three, overlapping nucleotide peaks in the same position in the GBSSI sequence. These polymorphic sites were always found within the intron regions, mainly at the beginning of the third intron. To ascertain the nature of these variants, the GBSSI DNA fragments were cloned and sequenced from seven of these individuals that displayed infraindividual polymorphism (S. cheesmaniae LA1450, S. chilense LA1930, S. peruvianum LA110, LA111, LA1336, LA1556, LA2152). Table 2 lists the number of GBSSI clones that were sequenced from these seven individuals, the number of distinct sequence types found, and their nucleotide variation and position. A putative microsatellite (AT)_n was found in the third intron, being n = 1 in potato and n = 2 up to n = 9 in the other accessions.

For some accessions (e.g., S. chilense LA1930 and S. cheesmaniae LA1450) all clones were identical, yet these differed from the sequences of the PCR product generated using Taq polymerase. We assumed that the sequences obtained from the clones were the most accurate, because Pfu polymerase produces less error. These sequences were included in the cladistic analysis. Sequences from five accessions of S. peruvianum differed by substitutions in the exon and intron regions and deletions/insertions only in the intron regions (Table 2). We initially included all sequence types in a cladistic analysis, but most were so similar (differing only by 1–3 base pairs) that they made no difference in the cladistic results. Others differed by up to 3 base pairs and were entered separately in the final analysis (indicated in italic boldface, Figs. 4, 5). In all cases they remained in the same part of the cladogram. In addition, cloned sequences of S. pennellii LA1376 has complex polymorphisms suggesting to us it is a possible hybrid; we do not include these results here but are investigating this further.

View larger version (20K): [in this window] [in a new window]   Fig. 4. One representative phylogram of 15 000 equally parsimonious trees (indels treated as characters), length = 347, consistency index = 0.88, retention index = 0.95. Different sequence types are indicated in italic boldface. Numbers indicate branch lengths, and unnumbered branches have a length of 1

Bootstrap analysis using 1000 replicates gave an excellent support for the outgroups (86–100%). Almost identical support values were found with the jackknife method. The ingroup had much less resolution and consisted of a basal polytomy, from which some branches were resolved, and a terminal clade. The basal polytomy was composed entirely of the self-incompatible green-fruited species S. chilense (map localities 29, 36, 38; Fig. 1), the central to southern Peruvian populations of S. peruvianum (24 populations, plus three additional sequence types, map localities 18, 19, 21, 22, 23, 24, 28, 29, 30, 32, 33, 34, 35, 36, 37), S. habrochaites (map localities 3, 11, 24), and S. pennellii (map localities 19, 28, 31). The latter two species are resolved within this basal polytomy, with bootstrap values of 86 and 96%, respectively. Three other branches of S. peruvianum were supported within this basal polytomy with bootstrap values of 62–99%. Two of these branches contain populations of S. peruvianum in the same generalized geographic area (map location 21) or contiguous generalized areas (21, 22, and 23).

The derived clade (70% bootstrap value) contains 14 northern Peruvian populations of S. peruvianum (also self-incompatible, green-fruited, map localities 6, 9, 10, 11, 12, 13, 15), and one central Peruvian population of S. peruvianum (LA1609, map locality 22), S. chmielewskii and S. neorickii (self-compatible, green-fruited, map localities 25, 26, 27 and 7, 17, 20, respectively), and the self-compatible and red- to orange-fruited species S. cheesmaniae, S. lycopersicum (map localities 5, 21, 32), and S. pimpinellifolium (map localities 2, 7, 8, 23). The latter three species cluster together as a polytomy within this advanced clade, but include one northern accession of S. peruvianum (LA2163, green-fruited, self-incompatible). The derived clade includes two groups of S. peruvianum with bootstrap values higher than 50%, and a branch corresponding to the red- to orange- to yellow-fruited self-compatible species, with a bootstrap value of 24%, containing S. cheesmaniae, S. lycopersicum, S. pimpinellifolium, and one accession of S. peruvianum. For branches with <50% bootstrap support the decay index = 1, for branches between 50 and 70% the decay index = 2, and for branches between 70 and 95% the decay index = 3 (Fig. 5).

	DISCUSSION

TOP ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS DISCUSSION LITERATURE CITED

 
The GBSSI sequences provided good resolution at the outgroup level and helped clarify relationships among five of the nine sections recognized by Child (1990) within Solanum subgenus Potatoe. Those results are entirely concordant with the cpDNA phylogeny of Solanum (Spooner, Anderson, and Jansen, 1993 ) regarding interspecific relationships of Solanum section Lycopersicon subsection Lycopersicoides, sect. Juglandifolium, sect. Basarthrum, sect. Petota, and sect. Etuberosum. The only difference is the cladistic separation of the three tuber-bearing members of sect. Petota, S. tuberosum on one branch and S. bulbocastanum + S. jamesii on a different branch, but this likely is a taxa sampling-density problem.

The GBSSI sequences provided less phylogenetic information at the ingroup than at the outgroup level, but still provided insights into taxonomic relationships among tomatoes. These results support S. juglandifolium and S. ochranthum as the closest outgroup to tomatoes, with S. lycopersicoides and S. sitiens as basal to these taxa, and support allogamy, self-incompatibility, and green fruits as primitive.

Mating system have played an important role in the evolution of wild tomatoes. Self-incompatibility has been considered as an ancestral condition to self-compatibility in tomatoes and probably never reversed to self-incompatibility (Rick, 1982 ). In this evolutionary context, changes from self-incompatibility to self-compatibility are expected to arise frequently and independently in different lineages (Rick, 1982 ). Interestingly, S. habrochaites and S. pennellii have both self-incompatible and self-compatible populations. The self-incompatible populations have higher genetic variation and occupy the center of their species geographic distribution. Self-compatible populations occur toward the northern and southern edges of distribution and have less genetic variation and smaller flower parts (Rick, Fobes, and Tanksley, 1979 ; Rick and Tanksley, 1981 ; Rick, 1984 ). The change from self-incompatibility to self-compatibility has been reported in only one population of S. peruvianum LA2157 from Chotano, Peru (Rick, 1986 ). Self-compatibility and the red to orange to yellow fruit color in S. cheesmaniae, S. lycopersicum, and S. pimpinellifolium have been considered as derived synapomorphies for this monophyletic group (Palmer and Zamir, 1982 ). The two closely related species, S. chmielewskii and S. neorickii, are also self-compatible. The latter species is exclusively autogamous with low population genetic variation and small flowers with stigmas included in the anther tube. In contrast, the facultative allogamous S. chmielewskii has larger flower parts and higher levels of heterozygosity. It has been postulated that S. neorickii evolved from S. chmielewskii (Rick et al., 1976 ). All the populations of S. chilense, and the most closely related species outgroups, S. lycopersicoides, S. sitiens, S. ochranthum, and S. juglandifolium, are exclusively self-incompatible. In Solanum sect. Petota, sect. Etuberosum and sect. Basarthrum self-compatibility most probably arose independently.

Within the ingroup, the phylogenetic analyses of GBSSI sequences identified S. pennellii and S. habrochaites as well-supported branches within the basal polytomy. Solanum peruvianum is supported as paraphyletic, with the taxa separating into two geographical groups, one in northern Peru ("northern group") and another into the central and southern Peru ("southern group"). All the red- to orange- to yellow-fruited self-compatible autogamous species (S. cheesmaniae, S. lycopersicum, and S. pimpinellifolium) and the green-fruited facultative allogamous species (S. chmielewskii and S. neorickii) are derived and on the same clade as these northern populations of S. peruvianum.

The GBSSI sequences presented very low variation in the autogamous species, compared with the considerable diversity found in all allogamous species. These results are in concordance with isozyme data (Rick and Fobes, 1975a, b ; Rick and Tanksley, 1981 ; Rick, 1983, 1986 ; Bretó, Asins, and Carbonell, 1993 ) and nRFLP data that showed self-incompatible species to harbor, on average, ten times more genetic variation within accessions than the self-compatible ones (Miller and Tanksley, 1990 ). In fact, more variation was found within a single accession of the self-incompatible species S. peruvianum than among all accessions tested of any one of the self-compatible species (Miller and Tanksley, 1990 ). The estimated amount of DNA polymorphism reflects breeding systems, consequently less within-species variation was found in the selfing species than in the outcrossing ones (Wolfgang and Langley, 1998 ). Similarly, low levels of DNA polymorphism were found in molecular linkage maps based on a cross between two closely related autogamous species, S. lycopersicum and S. pimpinellifolium, compared with the higher level of polymorphisms found in maps based on intraspecific crosses between S. lycopersicum and other allogamous species (Chen and Foolad, 1998 ). Results from different molecular markers demonstrated the difficulty of finding polymorphisms among self-compatible species. More variable markers, such as microsatellites used to identify tomato cultivars (Broun and Tanksley, 1996 ; Smulders et al., 1997 ), may be better suited to compare relationships among autogamous tomato species.

The low level of genetic variation found within self-compatible species may be explained by the role of autogamy that drives the loss of variation and fixation of alleles (Rick, 1984 ). The low variation among the self-compatible species showed by cpDNA (Palmer and Zamir, 1982 ), allozyme (Bretó, Asins, and Carbonell, 1993 ), nRFLP (Miller and Tanksley, 1990 ), molecular linkage maps (Chen and Foolad, 1998 ), and GBSSI sequences may reflect recent lineage divergence (Moniz de Sá and Drouin, 1996 ).

The maternal phylogenetic results from cpDNA restriction sites of Palmer and Zamir (1982) showed that the red- to orange- to yellow-fruited self-compatible autogamous species (S. cheesmaniae, S. lycopersicum, and S. pimpinellifolium) formed a clearly derived monophyletic group. The nRFLP showed the same results (Miller and Tanksley, 1990 ). The same group has been classified as series Lycopersicon in a morphological treatment by Child (1990 ; Fig. 2). The GBSSI data also were concordant in supporting the derived nature of the red-fruited autogamous species, and the monophyly of this group was confirmed (but with a low bootstrap value of 24%), with the exception of inclusion the single accession of S. peruvianum LA2163 from Yamalúc (Cajamarca Province, Peru). This northern accession came from the same area and is closely related genetically to the only self-compatible population of this species, S. peruvianum LA2157 from Tunel Chotano (Cajamarca Province, Peru). These northern Peruvian populations have been considered ancestral stocks (Rick, 1986 ), and S. peruvianum LA2163 may represent a progenitor population of this derived clade.

Critical results revealed by the GBSSI data are the paraphyly of S. peruvianum and the close relationship of the northern populations with the self-compatible taxa, also supported by nRLFP (Miller and Tanksley, 1990 ). The morphological and crossability differences between the northern and southern populations of S. peruvianum have long been recognized. The northern populations were separated as Lycopersicon [Solanum] peruvianum var. humifusum C. H. Müll. based on plants with decumbent stems, very short trichomes, and leaves with few leaflets (Müller, 1940 ). Rick (1963, 1986) demonstrated that these northern races were partially crossable among themselves but with reduced crossability to the southern populations that grow south of 8° S latitude (Rick, 1986 ). In addition, Miller and Tanksley (1990) showed that these two groups could be distinguished by nuclear restriction fragment length polymorphism data, although only one population of the northern group was analyzed, S. peruvianum var. humifusum LA2150 (not included in our GBSSI study). Within the northern S. peruvianum clade, only one central Peruvian population, S. peruvianum LA1609, was also found. Interestingly, this accession is geographically and morphologically related to other two coastal races of S. peruvianum, LA107 and LA1373 (not included in our GBSSI study), that have been used as crossing bridge populations between the northern and the southern populations of S. peruvianum (Rick, 1963, 1986 ). The inclusion of S. peruvianum LA1609 within the northern population clade suggests a possibility of naturally occurring gene exchange between northern and southern groups.

The GBSSI gene sequence results also showed the close relationship of two self-compatible species, S. chmielewskii and S. neorickii, with the northern S. peruvianum group. Similarly, the nRFLP supported a northern population, S. peruvianum var. humifusum LA2150, as the most closely related to S. chmielewskii and S. neorickii. In the maternal phylogenetic results from cpDNA restriction sites, S. chmielewskii also appeared closely related to a northern population, S. peruvianum LA1032 (not included in our GBSSI study), and Palmer and Zamir (1982) proposed the inclusion of S. chmielewskii as subspecies of S. peruvianum. In the phylogenetic GBSSI gene sequence results, the allogamous self-incompatible S. chilense was included within the basal polytomy with the southern races of S. peruvianum. A similar relationship was found with cpDNA restriction sites, and Palmer and Zamir (1982) also proposed the inclusion of S. chilense as subspecies of S. peruvianum.

Mitochondrial DNA (mtDNA) restriction fragment length studies were used to address relationships of nine Lycopersicon species and two closely related Solanum species (McClean and Hanson, 1986 ). The divergence estimates were based on the ratio between the number of shared fragments and all fragments evaluated. The mtDNA divergence is higher than that in cpDNA, indicating that the DNA of the two organelles is evolving at different rates (McClean and Hanson, 1986 ). Tomato mitochondrial genomes exhibit a small amount of divergence, and the variability was more limited within the self-compatible S. lycopersicum (0.04%) than between this species and the self-incompatible S. sitiens (2.7%). The inferred phylogeny from the GBSSI data confirms the position of this latter taxa as an outgroup species. Nevertheless, cpDNA, nRFLPs, and GBSSI results were not in agreement with the species relationships showed by mitochondrial results. These differences can be explained by the possibility of mitochondrial genome rearrangements that could result in estimates of divergences greater than the true values (McClean and Hanson, 1986 ). The incorrect assessment of homology of the mtDNA bands is another possible explanation that may account for the species relationship incongruence among data sets.

In agreement with crossing data, GBSSI gene sequence results distinguished the northern and southern S. peruvianum groups and also the close relationship of this species with S. chilense (Rick, 1963, 1979, 1986 ). Nevertheless, the molecular results of this study do not support the placement of S. pennellii and S. habrochaites within the Esculentum group (Rick, 1979 ), and S. chmielewskii and S. neorickii were more closely related to the northern S. peruvianum group than to the red- to orange- to yellow-fruited self-compatible species. The GBSSI results support a close relationship among green-fruited species, in agreement with morphological treatments of wild tomatoes (Müller, 1940 ; Luckwill, 1943 ).

Multiple GBSSI clones obtained from a single individual had identical sequences for S. chilense LA1930 and S. cheesmaniae LA1450 (Table 2). Nevertheless, multiple clones obtained from a single PCR performed on five different individuals, each one representing five S. peruvianum accessions, differed by nucleotide substitution and deletion/insertion events (Table 2). In four cases (S. peruvianum LA110, LA111, LA1336, and LA2152) more than two GBSSI sequence variants were identified. Theoretically, in diploid wild tomato species, a single nuclear gene as GBSSI has up to two different alleles in the heterozygote. True allelic variation at a single locus could explain the presence of two distinct sequences of GBSSI. However, the presence of more than two variants could be caused by artifactual errors that may occur as a consequence of nucleotide misincorporation during the PCR amplification of the GBSSI gene (Shimada and Tada, 1991 ). Another possibility to explain the observed intraindividual variation is the presence of more than one GBSSI locus (Salehuzzaman, Jacobsen, and Visser, 1993 ; Evans et al., 2000 ). In spite of the variation found among clones, all the variants used in with the cladistic analysis fell within the same clade.

The GBSSI DNA data, like other molecular markers used in tomatoes to date, were not able to completely differentiate very closely related, and most probably recent divergent species, especially the self-compatible ones. The GBSSI provide useful information about intraspecific variation in S. peruvianum and about character evolution regarding breeding systems and fruit color. A future paper will use these same accessions to investigate the morphological support for species boundaries and relationships.

	FOOTNOTES

² Current address: Department of Biological Sciences, National University of Cuyo, Almirante Brown 500, C. C. 7, 5505 Chacras de Coria, Lujan, Mendoza, Argentina (e-mail: iperalta{at}lab.cricyt.edu.ar ).

³ Author for reprint requests (e-mail: dspooner{at}facstaff.wisc.edu ).

	LITERATURE CITED

TOP ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS DISCUSSION LITERATURE CITED

 
Bamberg J. B. M. W. Martin J. J. Schartner D. M. Spooner 1996 Inventory of tuber-bearing Solanum species. Catalog of potato germplasm. Potato Introduction Station, National Research Support Program-6 (NRSP-6), Sturgeon Bay, Wisconsin, USA

Baum D. A. K. J. Sytsma P. C. Hoch 1994 A phylogenetic analysis of Epilobium (Onagraceae) based on nuclear ribosomal DNA sequences. Systematic Botany 19: 363-388[CrossRef][ISI]

Bohs L. R. G. Olmstead 1997 Phylogenetic relationships in Solanum (Solanaceae) base on ndhF sequences. Systematic Botany 22: 5-17

Bohs L. R. G. Olmstead 1999 Solanum phylogeny inferred from chloroplast DNA sequence data. In M. Nee, D. E. Symon, R. N. Lester, and J. P. Jessop [eds.], Solanaceae IV, advances in biology and utilization, 97–110. Royal Botanic Gardens, Kew, UK

Börner C. 1912 Solanum L. und Solanopsis gen. nov. Abhandlungen herausgegeben vom Naturwissenschaftlichen Verein zu Bremen 21: 282

Bremer K. 1988 The limits of amino acid sequence data in angiosperm phylogenetic reconstruction. Evolution 42: 795-803[CrossRef][ISI]

Bretó M. P. M. J. Asins E. A. Carbonell 1993 Genetic variability in Lycopersicon species and their genetic relationships. Theoretical and Applied Genetics 86: 113-120[ISI]

Broun P. S. D. Tanksley 1996 Characterization and mapping of simple repeat sequences in the tomato genome. Molecular and General Genetics 250: 39-49

Chen F. Q. M. R. Foolad 1998 A molecular linkage map of tomato based on a cross between Lycopersicon esculentum and L. pimpinellifolium and its comparison with other molecular maps of tomato. Genome 42: 94-103[CrossRef][ISI]

Child A. 1990 A synopsis of Solanum subgenus Potatoe (G. Don) D'Arcy [Tuberarium (Dun.) Bitter (s.l.)]. Feddes Repertorium 101: 209-235

Clark J. R. M. Robertson C. C. Ainsworth 1991 Nucleotide sequence of a wheat (Triticum aestivum L.) cDNA clone encoding the waxy protein. Plant Molecular Biology 16: 1099-1101[CrossRef][ISI][Medline]

D'Arcy W. G. 1987 The circumscription of Lycopersicon. Solanaceae Newsletter 2: 60-61

Donoghue M. J. R. G. Olmstead J. F. Smith J. D. Palmer 1992 Phylogenetic relationships of Dipsacales based on rbcL sequences. Annals of the Missouri Botanical Garden 79: 333-345

Doyle J. J. V. Kanazin R. C. Shoemaker 1996 Phylogenetic utility of histone H3 intron sequences in the perennial relatives of soybean (Glycine; Leguminosae) . Molecular Phylogenetics and Evolution 6: 438-447[CrossRef][ISI][Medline]

Dry I. A. Smith A. Edwards M. Bhattacharyya P. Dunn C. Martin 1992 Characterization of cDNAs encoding two isoforms of granule-bound starch synthase which show differential expression in developing storage organs of pea and potato. Plant Journal 2: 193-202[ISI][Medline]

Evans R. C. L. A. Alice C. S. Campbell E. A. Kellogg T. A. Dickinson 2000 The granule-bound starch synthase (GBSSI) gene in the Rosaceae: multiple loci and phylogenetic utility. Molecular Phylogenetics and Evolution 17: 388-400[CrossRef][ISI][Medline]

Farris J. S. 1989 The retention index and homoplasy excess. Systematic Zoology 38: 406-407[CrossRef]

Farris J. S. V. A. Albert M. Källersjö D. Lipscomb A. G. Kluge 1996 Parsimony jackknifing outperforms neighbor-joining. Cladistics 12: 99-124[CrossRef][ISI]

Felsenstein J. 1985 Confidence limits on phylogenies: an approach using the bootstrap. Evolution 39: 783-791[CrossRef][ISI]

Fosberg F. R. 1987 New nomenclatural combination for Galapagos plant species. Phytologia 62: 181-183

Heine H. 1976 Flora de la Nouvelle Caledonie, vol. 7. Musée Nacional D'Histoire Naturelle, Paris, France

Holle M. C. M. Rick D. G. Hunt 1978 Catalog of collections of green-fruited Lycopersicon species and Solanum pennellii found in watersheds of Peru. Tomato Genetics Cooperative Report 28: 49-78

Holle M. C. M. Rick D. G. Hunt 1979 Catalog of collections of green-fruited Lycopersicon species and Solanum pennellii found in watersheds of Peru. Tomato Genetics Cooperative Report 29: 63-91

Klösgen R. B. A. Gierl Z. S. Schwarz-Sommer H. Saedler 1986 Molecular analysis of the waxy locus of Zea mays. Molecular and General Genetics 203: 237-244[CrossRef]

Kluge A. G. J. S. Farris 1969 Qualitative phyletics and the evolution of anurans. Systematic Zoology 18: 1-32

Knapp S. D. M. Spooner 1999 A new name for a common Ecuatorian and Peruvian wild tomato species. Novon 9: 375-376[CrossRef][ISI]

Lester R. N. 1991 Evolutionary relationships of tomato, potato and pepino, and wild species of Lycopersicon and Solanum. In J. G. Hawkes, R. N. Lester, M. Nee, and N. Estrada [eds.], Solanaceae III: taxonomy, chemistry, evolution, 283–301. Royal Botanic Gardens, Kew, UK

Linnaeus C. 1753 Species Plantarum, 1st ed. Holmiae, Stockholm, Sweden

Luckwill L. C. 1943 The genus Lycopersicon: an historical, biological, and taxonomic survey of the wild and cultivated tomatoes. Aberdeen University Studies 120: 1-44

MacBride J. F. 1962 Flora of Peru: Solanaceae. Field Museum of Natural History, Botanical Series 13: 3-267

Mason-Gamer R. J. E. A. Kellogg 1996 Potential utility of the nuclear gene waxy for plant phylogenetic analysis. American Journal of Botany (supplement) 83: 178

Mason-Gamer R. J. C. F. Weil E. A. Kellogg 1998 Granule-bound starch synthase: structure, function, and phylogenetic utility. Molecular Biology and Evolution 15: 1658-1673[Abstract]

McClean P. E. M. R. Hanson 1986 Mitochondrial DNA sequence divergence among Lycopersicon and related Solanum species. Genetics 112: 649-667[Abstract/Free Full Text]

Mérida A. J. M. RodrÍguez-Galán C. Vincent J. M. Romero 1999 Expression of the granule-bound starch synthase I (Waxy) gene from snapdragon is developmentally and circadian clock regulated. Plant Physiology 120: 401-409[Abstract/Free Full Text]

Miller J. C. S. D. Tanksley 1990 RFLP analysis of phylogenetic relationships and genetic variation in the genus Lycopersicon. Theoretical and Applied Genetics 80: 437-448[ISI]

Miller P. 1754 The gardeners dictionary, abridged 4th ed. London, UK

Miller R. E. M. D. Rausher P. S. Manos 1999 Phylogenetic systematics of Ipomoea (Convolvulaceae) based on ITS and Waxy. Systematic Botany 24: 209-227[CrossRef][ISI]

Moniz de sá M. G. Drouin 1996 Phylogeny and substitution rates of angiosperm actin genes. Molecular Biology and Evolution 13: 1198-1212[Abstract]

Müller C. H. 1940 A revision of the genus Lycopersicon. United States Department of Agriculture, Miscellaneous Publication 382: 1-28

Olmstead R. G. J. D. Palmer 1997 Implications for phylogeny, classification, and biogeography of Solanum from cpDNA restriction site variation. Systematic Botany 22: 19-29

Olmstead R. G. J. A. Sweere R. E. Spangler L. Bohs J. D. Palmer 1999 Phylogeny and provisional classification of the Solanaceae based on chloroplast DNA. In M. Nee, D. E. Symon, R. N. Lester, and J. P. Jessop [eds.], Solanaceae IV, advances in biology and utilization, 111–137. Royal Botanic Gardens, Kew, UK

Palmer J. D. D. Zamir 1982 Chloroplast DNA evolution and phylogenetic relationships in Lycopersicon. Proceedings of the National Academy of Sciences, USA 79: 5006-5010[Abstract/Free Full Text]

Peralta I. E. 2000 Phylogeny of wild tomatoes (Solanum L. section Lycopersicum [Mill.] Wettst. subsection Lycopersicum) based on morphology and waxy gene sequences. Ph.D. dissertation, Plant Breeding and Plant Genetics Program, University of Wisconsin, Madison, Wisconsin, USA

Peralta I. E. D. M. Spooner 2000 Classification of wild tomatoes: a review. Kurtziana 28: 45-54

Rick C. M. 1960 Hybridization between Lycopersicon esculentum and Solanum pennellii: phylogenetic and cytogenetic significance. Proceedings of the National Academy of Sciences, USA 46: 78-82[Free Full Text]

Rick C. M. 1963 Barriers to interbreeding in Lycopersicon peruvianum. Evolution 17: 216-232[CrossRef][ISI]

Rick C. M. 1971 Lycopersicon. In I. L. Wiggins and D. M. Porter [eds.], Flora of the Galapagos Islands, 468–471. Stanford University Press, Palo Alto, California, USA

Rick C. M. 1973 Potential genetic resources in tomato species: clues from observations in native habitats. In A. M. Srb [ed.], Genes, enzymes, and populations, 255–269. Plenum, New York, New York, USA

Rick C. M. 1979 Biosystematic studies in Lycopersicon and closely related species of Solanum. In J. G. Hawkes, R. N. Lester, and A. D. Skelding [eds.], The biology and taxonomy of Solanaceae. Linnean Society Symposium Series 7, 667–677. Academic Press, New York, New York, USA

Rick C. M. 1982 Genetic relationships between self-incompatibility and floral traits in the tomato species. Biologisches Zentralblatt 101: 185-198[ISI]

Rick C. M. 1983 Tomato (Lycopersicon). In S. D. Tanksley and T. J. Orton [eds.], Isozymes in plant genetics and breeding part B, 147–165. Elsevier Science, Amsterdam, The Netherlands

Rick C. M. 1984 Evolution of mating systems: evidence from allozyme variation. In V. L. Chopra, B. C. Joshi, R. P. Sharma, and H. C. Bansal [eds.], Genetics: new frontiers. Proceedings of the XV International Congress of Genetics, vol. 4, Applied genetics, 215–221. Bowker, Epping, Essex, UK

Rick C. M. 1986 Reproductive isolation in the Lycopersicon peruvianum complex. In W. G. D'Arcy [ed.], Solanaceae, biology and systematics, 477–495. Columbia University Press, New York, New York, USA

Rick C. M. 1991 Tomato resources of South America reveal many genetic treasures. Diversity 7: 54-56

Rick C. M. J. F. Fobes 1975a Allozymes of Galapagos tomatoes: polymorphism, geographic distribution and affinities. Evolution 29: 443-457[CrossRef][ISI]

Rick C. M. J. F. Fobes 1975b Allozyme variation in the cultivated tomato and closely related species. Bulletin of the Torrey Botanical Club 102: 376-384[CrossRef][ISI]

Rick C. M. J. F. Fobes S. D. Tanksley 1979 Evolution of mating systems in Lycopersicon hirsutum as deduced from genetic variation in electrophoretic and morphological characters. Plant Systematics and Evolution 132: 279-298[CrossRef][ISI]

Rick C. M. M. Holle 1990 Andean Lycopersicon esculentum var. cerasiforme: genetic variation and its evolutionary significance. Economic Botany 44: 69-78[ISI]

Rick C. M. E. Kesicki J. F. Fobes M. Holle 1976 Genetic and biosystematic studies on two new sibling species of Lycopersicon from Interandean Peru. Theoretical and Applied Genetics 47: 55-68