|
|
||||||||
a Dep. of Botany, Univ. of Wisconsin, 430 Lincoln Drive, Madison, WI 53706-1381
b Vegetable Crops Research Unit, USDA-ARS, Dep. of Horticulture, Univ. of Wisconsin, 1575 Linden Drive, Madison, WI 53706-1590
* Corresponding author (dspooner{at}facstaff.wisc.edu)
| ABSTRACT |
|---|
|
|
|---|
Abbreviations: bp, base pairs cpDNA, chloroplast DNA IAM, infinite allele model PAUP, Phylogenetic Analysis Using Parsimony RAPD, randomly amplified polymorphic DNA RFLP, restriction fragment length polymorphism SMM, stepwise mutation model Tm, melting temperature
| INTRODUCTION |
|---|
|
|
|---|
The Andean and Chilean tetraploid cultivated potatoes have been treated as two separate species, S. andigenum Juz. and Bukasov and S. tuberosum L.; as subspecies of S. tuberosum; and as two cultivar groups within S. tuberosum (Huamán and Spooner, 2002). Cultivar groups are taxonomic categories used by the International Code of Nomenclature of Cultivated Plants (ICNCP) to associate cultivated plants with traits that are of use to agriculturists, but imply no phylogenetic differences between the groups. Species and subspecies, in contrast, are treated by the International Code of Botanical Nomenclature (ICBN) and generally assume phylogenetic differences (Spooner et al., 2002). Subspecies tuberosum is indigenous to Chiloé Island, the Chonos Archipelago to the south and adjacent areas in south-central mainland Chile. Subspecies andigenum is indigenous to Andean South America and ranges from Venezuela to northern Chile and Argentina (Hawkes, 1990).
Swaminathan and Magoon (1961) hypothesized that subsp. andigenum arose through autoploidy of the cultivated diploid species S. stenotomum Juz. and Bukasov. Cribb and Hawkes (1986) advocated that amphiploid subsp. andigenum originated from a cross between diploid cultivated S. stenotomum and the wild S. sparsipilum (Bitter) Juz. and Bakasov Matsubayashi (1991) suggested that subsp. andigenum originated from two diploid cultivated species S. phureja Juz. and Bukasov and S. stenotomum.
The origins of Chilean populations of subsp. tuberosum are speculative. Juzepczuk and Bukasov (1929) proposed that subsp. tuberosum originated from indigenous tetraploid cultivated species S. fonckii Phil. ex Reich. (nomen nudeum, no locallity other than Chile designated), S. leptostigma Juz. ex Bukasov (type from tubers collected in Chiloé Island), and S. molinae Juz. (type from tubers collected on Chiloé Island). Hawkes (1990) treated these three taxa as subsp. tuberosum. Juzepczuk and Bukasov (1929) suggested that S. palustre Schltdl. (then treated as S. brevidens Phil.) may be a progenitor. Hawkes (1990) proposed that subsp. andigenum evolved into subsp. tuberosum after transport to Chile, parallel with the evolution of subsp. tuberosum from subsp. andigenum in Europe. On the basis of starch grains, Ugent et al. (1987) proposed the wild species S. maglia Schltdl. as a progenitor of subsp. tuberosum. Grun (1990) hypothesized that subsp. tuberosum evolved from a cross between subsp. andigenum and an unidentified wild species.
The genetic relationships and extent of genetic differentiation between S. tuberosum subsp. tuberosum and subsp. andigenum is unresolved. On the basis of cytoplasmic sterility factors, geographical isolation, and ecological differences, Grun (1990) suggested that subsp. tuberosum was distinct from subsp. andigenum. Hawkes (1990) distinguished the two subspecies by subsp. tuberosum having fewer stems with foliage aligned at a broad angle to the stem and having less-dissected leaves with wider leaflets and thicker pedicels. Chloroplast DNA (cpDNA) restriction site data documented five chloroplast genotypes (A, C, S, T, and W types) in S. tuberosum, which included subsp. tuberosum and andigenum. Subspecies andigenum has all five types and native Chilean subsp. tuberosum has three types: A, T, and W (Hosaka and Hanneman, 1988). The most frequently observed type in Chilean subsp. tuberosum is T, which is characterized by a 241-base-pair deletion (Kawagoe and Kikuta, 1991).
Microsatellite or simple sequence repeats present a relatively new and promising technique to examine the genetic difference of subsp. tuberosum and subsp. andigenum. Microsatellites are tandemly repeated short sequences of 1 to 6 base pairs (bp) in length, spread throughout the genome, highly polymorphic, and have been used to investigate relationships among closely related taxa (Goldstein and Pollock, 1997). This study is designed to test differences between the two subspecies of S. tuberosum using the mapped potato microsatellite loci of Milbourne et al. (1998).
| MATERIALS AND METHODS |
|---|
|
|
|---|
|
|
|
Seeds from the U.S. Potato Genebank were planted in a greenhouse in Madison, WI, in May 1997. Tubers collected in Bolivia were planted in quarantine facilities in Glenn Dale, MD, in May 1999. Plants from both seeds and tubers were grown in pots in potting soil under ambient light. Young leaves were harvested after three months of growth. Leaves from plants grown in Maryland were sent on ice to Madison, WI, for DNA extraction. Deoxyribonucleic acid was extracted using the minipreparation methods outlined in Ballard et al. (1998). DNA concentrations were estimated by gel electrophoresis on 1% agarose gels with a low mass DNA ladder (Gibco BRL Low DNA Mass Ladder, Life Technologies, Rockford, MD)1.
Amplification and Fragment Analysis
Primers of 23 mapped microsatellite loci from Milbourne et al. (1998) (stm0007, stm0010, stm0013, stm0019, stm0025, stm0037, stm0038, stm0051, stm0052, stm1003, stm1006, stm1008, stm1017, stm1020, stm1029, stm1031, stm1049, stm1055, stm1069, stm1100, stm1104, stm2020, and stm3016) were labeled with fluorescent dyes (FAM, HEX, TETRA from Applied Biosystems, Foster City, CA) and used to amplify DNA from one seedling per accession (Table 1). Reaction conditions were optimized as described in Provan et al. (1996). Conditions for a 25-µL reaction were as follows: 1XPCR Buffer II (Perkin Elmer, Norwalk, CT), 1.5 mM MgCl2, 0.2 mM deoxyribonucleoside triphosphates, 0.4 µM of each primer pair (labeled forward and unlabeled reverse), 1 unit of AmpliTaq Gold (Perkin Elmer), and 10 to 20 ng of DNA. All reactions were amplified in a Perkin Elmer 9600 thermocycler set at the following times and temperatures: 1 cycle of 10 min at 94°C, 2 min at melting temperature (Tm) (Milbourne et al., 1998), 5 min at 72°C, followed by 29 cycles of 1 min at 94°C, 45 s at Tm, 5 min at 72°C, and ending with 45 min hold at 72°C. The Tm for each primer followed Milbourne et al. (1998).
Usually, three microsatellite products from one individual were pooled by adding 3 µL of the FAM, 6 µL of TETRA, and 12 to 15 µL of HEX, depending on the strength of the product. An aliquot of 1.5 µL of each pooled sample was mixed with 1.8 µL of sample loading buffer [80% formamide, 10 mg mL-1 blue dextran, 5 mM Disodium(ethylene dinitrilo)tetraacetic acid dihydrate pH 8.0] and 0.45 µL molecular weight standard (TAMRA 5000 standard from PE Biosystems, Wellesley, MA). Samples were heated for 3 min at 95°C, then chilled on ice. Approximately 1 µL of the chilled sample was loaded on a 5% LongRanger (FMC Bioproducts, Rockland, ME) polyacrylamide/6 M urea gel in a PE Biosystems 377XL DNA sequencing apparatus. Samples were run at 3000 V at 2400 scans hr-1 in 36-cm well-to-read plates. Data were collected using the DNA Sequencer Data Collection v. 2.0 (PE Biosystems) and analyzed with GeneScan v. 2.1 (PE Biosystems).
Fluorescent peaks were labeled as fragment sizes by using Genotyper v. 2.1 (PE Biosystems). All peaks were manually edited using Genotyper's manual click option. Peaks were scored by allele size and as presence or absence. One replicate of identical DNA (subsp. tuberosum, PI 595455) was included as a standard, as well as two accessions with DNA extracted from separate seedlings of the same accession.
Data Analysis
The allelic bp scores were entered into Microsat version 1.5d (Minch et al., 1997) to generate a distance matrix using
µ2, the stepwise mutation model (SMM) of Goldstein et al. (1995). This model assumes that new mutations occur by stepping up or down in size within a short region of DNA, and thus could lead to a previously existing allele. We also generated distance matrices by an infinite allele model [IAM, Nei (1972) in Phylogenetic Analysis Using Parsimony (PAUP)] that assumes each mutation forms a new allele of infinite size.
Both distance matrices were entered into PAUP version 4.0d65 (Swofford, 1998), where branching trees were built using the tree building methods of unweighted pair group method with arithmetic means (UPGMA; Sokal and Sneath, 1963), neighbor-joining (Saitou and Nei, 1987), and unweighted least-squares methods (Fitch and Margoliash, 1967). All trees were built using a random input order. The unweighted least-squares tree was generated using a heuristic stepwise-addition search for 10 replicates. The distance matrices were also entered into KITSCH, a program within PHYLIP version 3.573c (Felsentein, 1995), to obtain a Fitch-Margoliash method with contemporaneous taxa (power = 2.0; 10 random replicates). Significant differences in allele distributions were determined in JMP statistical software (SAS Institute, 1995) by the Tukey-Kramer honestly significant difference (HSD) test.
| RESULTS |
|---|
|
|
|---|
Of the 18 loci scored, a total of 208 alleles were detected for all 94 individuals. Both subspecies of S. tuberosum accounted for 46.5% of the total missing data and remaining accessions for 53.5%. Each primer produced one to four peaks per individual scored as numbers of bp for the tetraploid individuals. The diploid individuals had three identifiable peaks, with the source of the third peak not known. The number of heterozygous individuals (more than two peaks scored for one primer) across all loci ranged from 20 individuals in stm1029 to 77 individuals in stm2020, with a mean of 46.8 per locus. The size ranged from 11 bp in stm1008 to 122 bp in stm0019. The number of alleles per locus ranged from 6 in stm1069 to 19 in stm0019, with a mean of 11.6 alleles per locus.
Allelic distributions varied by ploidy, taxonomy, and taxonomic distance (ingroup vs. outgroup) (Table 2) . The two subspecies of S. tuberosum showed no difference in mean number of null alleles per accession or mean number of heterozygous alleles per accession. However, mean number of alleles per accession were significantly (P = 0.05) different between the two subspecies. All other comparisons were highly significant (P = 0.05), with the greatest difference occurring in the number of null alleles present in the ingroup and the outgroup comprised of mostly tetraploids and diploids, respectively.
|
µ2 analysis. Primers stm1104 and stm1017 showed a one-bp difference and we treated them as single nucleotide repeats in the
µ2 distance calculation. Primer stm1029, which amplifies a single nucleotide repeat region, showed alleles of two-bp differences; thus we considered it a dinucleotide repeat.
Comparison of IAM and SMM Trees
The IAM model consistently generated trees that placed the replicate sample and the two individuals from duplicate accessions near each other (Fig. 2). Conversely, all trees generated from SMM separated the replicate samples and duplicate accessions, intermixing subspecies of S. tuberosum, cultivated species, and S. brevicaule taxa.
The IAM or neighbor-joining tree (Fig. 2) separated 33 of the 35 accessions of subsp. tuberosum, intermixed with S. maglia. It placed two accessions of subsp. tuberosum with subsp. andigenum that also contained eight accessions of wild or other cultivated species. Two accessions of subsp. andigenum grouped with the outgroups. The tree placed the cultivated species S. x ajanhuiri Juz. & Bukasov, S. x curtilobum Juz. & Bukasov, S. phureja, and S. stenotomum within the ingroup, and the outgroups in the basal branch but intermixed with two accessions of S. tuberosum subsp. andigenum and seven other ingroups.
| DISCUSSION |
|---|
|
|
|---|
Although microsatellites are claimed to be useful for species-level phylogenies (Takezaki and Nei, 1996), it appears there is no consensus among researchers as to which evolutionary model is most appropriate for reconstructing phylogenies based on microsatellite data (Feldman et al., 1999; Goldstein and Pollock, 1997). Therefore, trees of S. tuberosum were constructed using both the SMM model (
µ2) (Goldstein et al., 1995) and the IAM model of Nei (1972). Both models failed to absolutely distinguish subsp. andigenum from subsp. tuberosum, or from the other cultivated species. Neither method will clearly separate subsp. andigenum from some of its ingroup relatives in the S. brevicaule complex and other cultivated species (Grun, 1990; Miller and Spooner, 1999).
It is possible that neither the SMM, which calculates distances based on the number of repeats, nor the IAM, which evaluates all possible alleles, provides the best assessment of the diversity and substructure of S. tuberosum (Di Rienzo et al., 1994). The SMM and IAM models were designed for diploids, thus neither method may best estimate the distance between the accessions examined. The microsatellite results (Fig. 2) are not totally congruent with the outgroup relationships of Spooner and Castillo (1997), suggesting that microsatellite data show homoplasy with more distant related taxa. Decreasing utility of microsatellites with increasing taxonomic distance is also clear from the much greater number of null alleles in the outgroup (Table 2), as is common in most other systems studied (Roa et al., 2000). Because the source of the primer sequences is S. tuberosum (Milbourne et al., 1998), our ingroup results are more likely to be valid.
Taxonomy and Origins of Solanum tuberosum
Various hypotheses have been proposed to explain the origins of the cultivated tetraploid potato. Infinite allele model or neighbor-joining results (Fig. 2) support the concept that most of the Chilean populations of subsp. tuberosum are distinct from subsp. andigenum, and other diploid and tetraploid wild and cultivated populations. Subspecies andigenum and tuberosum are separable but only by using a complex of morphological characters that have overlapping character states (Huamán and Spooner, 2002). They also differ in daylength, adaptation to flowering, and most populations of subsp. tuberosum have a cpDNA deletion mutation. Many studies (Grun, 1990; Huamán and Spooner, 2002) suggested that cultivated potatoes are possibly of multiple origins, and that they hybridize with wild species.
The microsatellite data place S. maglia with subsp. tuberosum, which could be interpreted to support this species as its progenitor, but in agreement with Cribb and Hawkes (1986), Hosaka and Hanneman (1988), and Hawkes (1990), we consider this unlikely. Solanum maglia is restricted to coastal Chile, 1000 km north of native S. tuberosum, with a single population in the mountains of Argentina, and most populations are sterile triploids (Hawkes, 1990). Identification of extant populations of S. maglia in southern Chile by Ugent et al. (1987) are not backed up by voucher specimens of others who have collected potatoes extensively in Chile (Contreras-M, 1987; Spooner et al., 1991). It is possible that S. maglia is a diverged population of S. tuberosum.
In traditional Andean farming systems, there is ample room for mixing wild and weedy species with cultivated potatoes. Fields of native cultivated potatoes in the Andes are a mixture of various ploides (Ochoa, 1958), and these frequently occur with weedy wild relatives. Farmers select and maintain varieties on factors such as yield, disease resistance, storage longevity, and taste (Brush et al., 1981). As a result, they have developed a large diversity of cultivated varieties that are ecologically versatile (Brush et al., 1981; Zimmerer, 1995). Because of this selection within a loose agroecosystem, the species of Andean potatoes likely cross and form a large, plastic genepool. At the genotype level, there is considerable heterogeneity and endemism of native potatoes (Zimmerer and Douches, 1991), but this variability is difficult to divide into species.
In conclusion, nuclear microsatellites developed for S. tuberosum analyzed with IAM or a neighbor-joining model (Fig. 2) support the genetic difference of most landrace populations of subsp. tuberosum and andigenum. These support their recognition at some classification level.
| ACKNOWLEDGMENTS |
|---|
| NOTES |
|---|
|
|
|---|
Received for publication May 21, 2001.
| REFERENCES |
|---|
|
|
|---|
Related articles in Crop Science:
This article has been cited by other articles:
![]() |
D. M. Spooner DNA barcoding will frequently fail in complicated groups: An example in wild potatoes Am. J. Botany, June 1, 2009; 96(6): 1177 - 1189. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. M. Spooner, J. Nunez, G. Trujillo, M. del Rosario Herrera, F. Guzman, and M. Ghislain Extensive simple sequence repeat genotyping of potato landraces supports a major reevaluation of their gene pool structure and classification PNAS, December 4, 2007; 104(49): 19398 - 19403. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. Rios, M. Ghislain, F. Rodriguez, and D. M. Spooner What Is the Origin of the European Potato? Evidence from Canary Island Landraces Crop Sci., May 31, 2007; 47(3): 1271 - 1280. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| The SCI Journals | Agronomy Journal | Vadose Zone Journal | |||
| Journal of Natural Resources and Life Sciences Education |
Soil Science Society of America Journal | ||||
| Journal of Plant Registrations | Journal of Environmental Quality |
The Plant Genome | |||