ISOZYMATIC VARIATION AND PHYLOGENETIC

Transcription

1 American Journal of Botany 86(1): ISOZYMATIC VARIATION AND PHYLOGENETIC RELATIONSHIPS BETWEEN HENEQUÉN (AGAVE FOURCROYDES) AND ITS WILD ANCESTOR A. ANGUSTIFOLIA (AGAVACEAE) 1 PATRICIA COLUNGA-GARCÍAMARÍN, 2,3 JULIÁN COELLO-COELLO, 2 LUIS E. EGUIARTE, 4 AND DANIEL PIÑERO 4 2 Centro de Investigación Científica de Yucatán, Apartado Postal 87, Cordemex, Mérida, Yucatán, México 97310; and 4 Instituto de Ecología, UNAM. Apartado Postal Ciudad Universitaria, México, D.F. Isozymatic variation and phylogenetic relationships among extant henequén (Agave fourcroydes) germplasm and wild populations of its ancestor A. angustifolia in the Yucatan Peninsula in Mexico were analyzed. Analysis of three isozyme systems using starch gel electrophoresis indicated that while A. angustifolia populations have relatively high levels of variation, within each henequén cultivar all individuals were identical. This result corresponds to previous ethnobotanical and morphological analyses, which indicated severe loss of genetic variation of this domesticated plant as a consequence of the promotion by means of asexual propagation of only one cultivar since the middle of the last century. The three extant cultivars of henequén were distinct from each other. Two of them, Sac Ki (SK) and Yaax Ki (YK), could be matched within the progenitor, but Kitam Ki (KK) has a MDH electrophenotype not found in any of the plants growing inside the Yucatan Peninsula, but found in some A. angustifolia plants growing in the Mexican states of Oaxaca and Veracruz. A parsimony analysis of the morphological data indicated two lineages: that of SK and YK, cultivated cordage plants selected for stronger and longer fibers, whose sister group is the Tropical subdeciduous forest ecotype (SF), and that of all the other wild populations, which also included KK, the cultivated textile plants selected for finer fibers and nearly extinct in Yucatan. These results support the hypothesis of the yucatecan origin of SK and YK from the SF ecotype, as well as the hypothesis of a recent introduction of KK to the Yucatan Peninsula in a domestication trend that probably included also Chelem White (its cultivation being abandoned later). Key words: Agavaceae; Agave angustifolia; Agave fourcroydes; domestication; germplasm diversity; henequén; isozymes; phylogeny. Henequén, Agave fourcroydes Lem., is a plant cultivated for its fiber that was domesticated by the Prehispanic Maya in the Yucatan Peninsula. Its origin has been attributed to A. angustifolia Haw. (sensu Gentry, 1982), the only wild species of Agave growing in this area (Colunga-GarcíaMarín and May-Pat, 1993). Castorena-Sán- 1 Manuscript received 24 November 1997; revision accepted 22 May The authors thank Robert Bye and Abisaí García from the Jardín Botánico, UNAM, who donated plants of A. angustifolia from the states of Mexico, Veracruz, and Oaxaca to the Jardín Botánico Regional, CICY; Robert Bye for his support and advice during the course of the project; Valeria Souza, Roger Ashburner, Germán Carnevali, Ken Oyama, Marlene de la Cruz, Victor Chávez, and Paul Gepts for comments and suggestions on an early version of this paper; Karen H. Clary, David J. Boegler, and an anonymous referee for reviewing the manuscript; Filogonio May-Pat for his technical assistance in maintaining the live collection; Nidia Pérez for training the senior author in starch gel electrophoresis procedures; and Sergio Zárate for his help in the English translation.this paper is part of the Doctoral Dissertation that the senior author completed at the Instituto de Ecología of the Universidad Nacional Autónoma de México (UNAM). Financial aid was provided in part by the Comisión Nacional de Biodiversidad (CONABIO), the New York Botanic Garden through the PREBELAC program, and by UACPyP-UNAM through the PADEP program. The senior author is also grateful to the Consejo Nacional de Ciencia y Tecnología (CONACyT) for receiving support during her doctoral studies through the Programa de Cátedras Patrimoniales de Excelencia. L. E. Eguiarte received financial aid from the PAPITT project from D.G.A.P.A. UNAM, No. IN Author for correspondence ( pcolunga@cicy.cicy.mx). 115 chez, Escobedo, and Quiroz (1991) described the karyotypes of these species in Yucatan as polyploids, i.e, henequén having a chromosome number of 5x 150 and A. angustifolia with 6x 180. Within the genus, A. angustifolia has the broadest distribution, ranging along the Pacific coast, from the Mexican state of Sonora to Costa Rica, and along the Atlantic coast from Costa Rica up to Tamaulipas, at altitudes from sea level to 1500 m (Gentry, 1982). In the Yucatan Peninsula, this species exhibits a morphological variation gradient related to its geographic distribution (Orellana et al., 1985), from small plants associated with Coastal dunes, to intermediate sized plants in Tropical deciduous forests, to larger plants in Tropical subdeciduous forests. Agave angustifolia forms aggregations of individuals through vegetative propagation, either by rhizomes or bulbils, but is also endowed with a high capacity to reproduce sexually. Present morphological variation led Gentry (1982) to assume that it represented an open-pollinated complex, which has undergone considerable morphological change in response to habitat conditions, long-term climatic changes, and human activities. As far as we know, no research has been published about the reproductive biology of the species. The present paper analyzes the isozymatic variation and the phylogenetic relationships among existing henequén varieties and wild populations of A. angustifolia that grow in the Yucatan Peninsula. Results are discussed in the context of the hypotheses postulated in previous pa-

2 116 AMERICAN JOURNAL OF BOTANY [Vol. 86 TABLE 1. Number of populations and individuals per population analyzed per variant. Variant No. of populations No. of individuals per population Total no. of individuals A. angustifolia (all wild) Coastal dunes (D) 5 3, 5, 3, 4, 4 19 Tropical deciduous forest (DF) 2 2, 8 10 Tropical subdeciduous forest (SF) 3 1, 5, Chelem White (CHW) Chelem Yellow (CHY) Chelem Green (CHG) Henequén (all cultivated) 6 76 Sak Ki (SK) 5 5, 5, 5, 5, 7 27 Yaax Ki (YK) 5 7,7,2,4,5 25 Kitam Ki (KK) 3 2, 7, pers (Colunga-GarcíaMarín and May-Pat, 1993; Colunga- GarcíaMarín et al., 1993; Colunga-GarcíaMarín, Estrada- Loera and May-Pat, 1996; Colunga-GarcíaMarín and May-Pat, 1997), and in terms of alternatives for conservation of wild and cultivated germplasm. Henequén is a crop plant whose successful cultivation has been restricted to the Yucatan Peninsula, some regions of the Mexican states of Tamaulipas and Veracruz, and Cuba. In these last three areas, the germplasm was introduced from Yucatan. Propagation is strictly vegetative, by means of rhizomatous suckers. A prevailing agricultural practice in henequén cultivation is the severing of the inflorescence peduncle soon after it begins to emerge, making it unusual to encounter flowering plants. When flowers do develop, they seldom form seeds, and even then the seeds are never used for propagation. Germination percentage of henequén seeds is 9%, while that of wild plants is 73% (P. Colunga-GarcíaMarín and F. May-Pat, unpublished data). Previous studies were made of the ethnobotanical evidence on past germplasm diversity of this cultigen (Colunga-GarcíaMarín and May-Pat, 1993) and morphological variation in natural and uniform conditions of both henequén and its putative wild ancestor (Colunga- GarcíaMarín, Estrada-Loera, and May-Pat, 1996; Colunga-GarcíaMarín and May-Pat, 1997). These studies indicated the lack of evidence about the diversity of henequén during the Prehispanic Mayan era, but we assume that this diversity was at least equal or larger than that published in agronomic manuals around the beginning of the 20th century (de Echánove, 1814; Regil and Peón, 1853; Espinosa, 1860; Barba ; Bolio, 1914). Those manuals report on seven varieties of henequén and the experimental cultivation of wild plants. At present, only three of these varieties can be found: Sac Ki (SK), Yaax Ki (YK), and Kitam Ki (KK) (abbreviations also shown in Table 1). Variety SK is preferred for making cordage, because it has large, coarse, and abundant fibers. YK is morphologically very similar to SK, but it produces slightly shorter fibers and in slightly lesser amounts. KK is practically unknown to producers, and despite the fact that it is always found cultivated, it is often confused with a wild plant. It has a considerably different morphology from the other two varieties. It produces significantly less fiber, which is soft and short. Both YK and KK are found only in very small populations, particularly KK, which is only cultivated for craft-manufacturing. Ethnobotanical research (Colunga-GarcíaMarín and May-Pat, 1993) suggested that an accelerated loss of genetic diversity occurred for this crop at the beginning of this century. This loss is the result of the establishment of extensive plantations of SK exclusively dedicated to the cordage industry, accentuated by the exclusive practice of vegetative propagation. Measured coefficients of variation for the morphological characters studied (Colunga-GarcíaMarín, Estrada-Loera, and May-Pat, 1996; Colunga-GarcíaMarín and May-Pat, 1997) also suggested limited genetic diversity of the crop with respect to the wild ancestor. Morphological evidence (Colunga-GarcíaMarín, Estrada-Loera, and May-Pat, 1996; Colunga-GarcíaMarín and May-Pat, 1997) suggested that SK and YK both differ from wild plants in similar direction and magnitude, a trend that may be described as four domestication syndromes (domestication syndrome refers to a combination of characteristics with anthropocentric interest or related to the process of artificial selection). These syndromes are: gigantism, enhanced fiber content, less thorniness, and diminished reproductive capacity. There is a conspicuous correspondence among these syndromes and the anthropocentric interests that guided the process of artificial selection of henequén, at least during the former century. KK is the cultivated variety that most resembles wild plants. Ethnobotanical and morphological evidence suggests that KK was recently introduced to the Yucatan Peninsula and/or that it went through a selection process having different direction and intensity with respect to that in SK and YK (Colunga-GarcíaMarín and May-Pat, 1997). Variety SK clearly matches Gentry s (1982) diagnosis of A. fourcroydes. Given that valid and complete taxonomic descriptions for varieties YK and KK are lacking, in this paper we refer to them by their Mayan horticultural designations. With respect to the A. angustifolia populations, ethnobotanical exploration (Colunga-GarcíaMarín and May- Pat, 1993) suggested the existence of three variants of A. angustifolia, which correspond to the different habitats in which the populations grow: coastal dunes (D); Tropical deciduous forests (DF); and Tropical subdeciduous forests (SF). According to the morphological evidence (Colunga-GarcíaMarín and May-Pat, 1997), these three variants belong to two ecotypes: one including D and DF populations and the other corresponding to SF populations. This conclusion is derived from the fact that under uniform growing conditions no morphological differences were found between D and DF populations. Artisans who use wild plant fibers distinguish three variants within SF populations according to fiber quality. Of these, Chelem White (CHW) is considered more similar to cultivated variants, while both Chelem Green (CHG) and Chelem Yellow (CHY) are seen as of lower quality (in that order). This artisan classification agrees with morphological evidence (Colunga-GarcíaMarín and May-Pat, 1997). Based on ethnobotanical and morphological evidence, three hypotheses were postulated with respect to

3 January 1999] COLUNGA-GARCÍAMARÍNETAL. ISOZYMATIC VARIATION OF HENEQUÉN 117 Fig. 1. Mother-plant origin of material analyzed and location of the Jardín Botánico Regional of CICY (JBR) where the material was grown. these wild variants used by artisans (Colunga-García- Marín and May-Pat, 1993, 1997): (1) SF populations gave rise to cultivated plants; (2) CHW represents ancestral henequén populations; or (3) CHW is a clone introduced to cultivation and later abandoned during the turn of the century, a period during which attempts were made to cultivate wild populations for the textile industry. MATERIALS AND METHODS Plant material Isozyme analysis was done using 146 plants belonging to the six variants of A. angustifolia and to the three variants (SK, YK, and KK) of henequén (A. fourcroydes) found in the Yucatan Peninsula. These plants were rhizomatous suckers from plants established under uniform growth conditions in the Jardín Botánico Regional of the Centro de Investigación Científica de Yucatán (JBR) between 1985 and Morphological variation of mother plants was evaluated by Colunga-GarcíaMarín and May-Pat (1997), using eight wild A. angustifolia individuals from three Mexican States outside the Yucatan Peninsula: three from Veracruz, three from Oaxaca, and two from the State of Mexico. Each one came from a different locality. These plants were grown in the agave collection of the Jardín Botánico at UNAM, Mexico City, and then transferred and acclimated in the JBR for 1 mo previous to this study. Mother-plant origin of material analyzed is presented in Fig. 1. Table 1 shows the number of populations of each variant analyzed and the number of individuals per population included from the Yucatan Peninsula. All plants analyzed were 2 yr old at the time of sampling. To perform phylogenetic analysis we used morphological data from 136 plants grown under uniform conditions at the JBR and evaluated by Colunga-GarcíaMarín and May-Pat (1997). Isozyme analysis Tissue was collected from the leaf most recently detached from the bud in the center of the rosette and consistently sampled between 0800 and The terminal spine from this leaf was excised and 100 mg of tissue were taken from the adjacent area. Tissue samples were macerated in 100 ml of extraction buffer. The macerate was absorbed onto filter paper wicks, which were stored for 24 h at 10 C. After this, wicks were placed in starch gels (12%, containing 3% sacarose). Starch gels were run horizontally following the Wendel and Weeden (1989) protocols. In a preliminary study, 24 enzymes systems in 11 gel systems were assayed using four different extraction buffers (results available from first author). Only three of these 24 enzyme systems gave satisfactory activity, good band resolution, and consistent results and were thus selected for this study. These isozymes were acid phos-

4 118 AMERICAN JOURNAL OF BOTANY [Vol. 86 phatase (ACP) (E.C. [Enzyme Commission number] ), cathodic peroxidase (PRX) (E.C ), and malate dehidrogenase (MDH) (E.C ). Gel systems found optimal for MDH were system D (ph 6.5) designed by Stuber et al. (1988) for maize and for ACP and PRX system A (ph 5.0), also designed by Stuber et al. (1988) for maize. The best extraction buffer was that of Soltis et al. (1983). Staining protocols were for MDH that of Stuber et al. (1988) and for ACP and PRX those of Wendel and Weeden (1989), modified as follows: (1) ACP, 4 ml of 1 mol/l sodium acetate, ph 5.0, 50 ml distilled water, 100 mg magnesium chloride, 20 mg -naphtyl acid phosphate, and 75 mg of Fast Garnet (GBC) Salt, (2) PRX, 3.5 ml N-N dimethyl formamide, 35 ml distilled water, 2.5 ml 1 mol/l sodium acetate ph 5.0, 100 mg calcium chloride, 0.5 ml 3% hydrogen peroxide, and 50 mg 3-amino-9-ethyl carbazole. Due to polyploidy, activity bands could not be assigned to specific loci, given that they probably represented product combinations from different loci on homologous chromosomes. As a consequence, banding patterns were recorded as the presence-absence of individual bands: 14 bands for ACP, four bands for PRX, and nine for MDH (Fig. 2). For each isozyme system, bands were numbered consecutively, beginning with the band running closer to the origin (bottom in Fig. 2). A number was also assigned to each pattern of presence-absence of bands, or isozyme electrophenotype, and their frequency was recorded for (1) each variant analyzed, (2) pooled wild variants, and (3) pooled cultivated variants (Table 2). Phylogenetic analysis Phylogenetic Analysis Using Parsimony was done using PAUP (Swofford, 1993). An exhaustive search of the shortest tree was done based on a 66-character matrix (Colunga- GarcíaMarín, 1996). This matrix was constructed following the ideas of Thorpe (1984) and Stevens (1991) about the way of producing discrete character states from continuous variation. States were determined by the number of significant differences between taxa in a general comparision for 55 quantitative morphological characters evaluated by Colunga-GarcíaMarín and May-Pat (1997). Three Analysis of Variance (ANOVA) for unbalanced designs were performed by the minimum square method of adjustment to general linear models, using the GLM procedure of SAS, 6.04 (SAS, 1988). In the first analysis, differences were tested for all studied taxa for which data of all 21 vegetative characters included were available. In a second analysis, differences were tested among the six wild taxa and the cultivated one (SK) for all 17 flower and inflorescence characters for which data were available. In the third analysis, differences were tested among the six wild taxa for which data of the 17 fruit characters considered were available. The differences between means were analyzed with the Tukey-Kramer method (SAS, 1988). The significance level was adjusted according to the number of simultaneous comparisons made in each analysis of variance, by means of the Bonferroni inequallity (Miller, 1981). Previous to analyses, tests for normality were made on variable residuals, and these were transformed as described by Colunga-GarcíaMarín and May-Pat (1997). For coding, character states were considered to be different when means showed significant differences ( 0.05). Besides these 55 originally quantitative characters, another 11 qualitative characters were added. Of the 66 morphological characters, order was declared for 44 of these, based on the analysis of henequén domestication syndromes of Colunga-GarcíaMarín, Estrada-Loera, and May-Pat (1996) and Colunga-GarcíaMarín and May-Pat (1997), assuming that the ancestral character state was equal to 0. Characters used, their states and order are also available in Colunga-GarcíaMarín (1996). Wild variant from Dunes (D) was considered the outgroup, because its populations are perhaps the oldest in the Yucatan Peninsula. Character changes were mapped using Mac Clade 3.0 (Maddison and Maddison, 1992) to find apomorphies. Results were evaluated with a bootstrap analysis based on 1000 replicates with two random input orders and a decay analysis using the total number of possible trees, considering trees up to 10 steps longer than the most parsimonious tree. RESULTS Isozymatic variation In the 146 individuals studied from the Yucatan Peninsula, we found 23 electrophenotypes of ACP, nine of PRX, and two of MDH (Table 2). The number of electrophenotypes registered for ACP and PRX is very high for the wild species and very low for the henequén. While we registered 22 different electrophenotypes of ACP and nine of PRX in the six variants of A. angustifolia, only two electrophenotypes of each one of these isozymes were registered in the henequén varieties. The total allocation of genetic variation among wild and cultivated variants is evident; practically all the variation observed falls within the wild populations. In the case of the MDH, all plants showed the same electrophenotype, with the exception of those plants belonging to the cultivar KK. The three cultivars of henequén were distinct from each other, but no differences in enzyme activity band presence were found among individuals. This means that when genetic variation is estimated by means of the three isozyme systems analyzed all genetic variation is found among populations. The differences between SK and YK were found only on one of the PRX bands, while KK differs from them in both MDH and ACP electrophenotypes. Within A. angustifolia variants, the highest variation level was recorded in D (21 different electrophenotypes) and the lowest in CHW (five different electrophenotypes), which is the variant considered by producers to be the most similar to SK. In wild A. angustifolia populations most of the genetic variation is found within variants (Table 2). Neverthless, a significant amount of variation is present among variants, as all of them (except CHW) have electrophenotypes not present in the others. The most abundant ACP electrophenotype in A. angustifolia populations was the same as that of the SK and YK. KK had a unique electrophenotype, as did 11 A. Fig. 2. Banding patterns of the three isozyme systems analyzed. (A) Electrophoretic variation of ACP in individuals of Agave angustifolia from Dunes (D) (lanes 1 19) and Tropical deciduous forest (DF) (lanes 20 24) populations: lanes 1 3 population no. 1, lanes 4 8 population no. 2, lanes 9 11 population no. 3, lanes population no. 4, lanes population no. 5 from D; lanes population no. 1, lanes population no. 4 from DF. (B) Electrophoretic variation of PRX in individuals of A. angustifolia from D: lanes 1 3 population no. 1, lanes 4 8 population no. 2, lanes population no. 3, lanes population no. 4, lanes population no. 5, lanes are a duplicate from the individual in lane 22; lanes 9 and 17 are from an individual from DF used as a reference; lanes 10 and 18 an individual of the cultivated variant Kitam Ki (KK) also used as a reference. (C) Electrophoretic variation of MDH in individuals of the three henequén varieties: lanes 1 5 Sac Ki (SK), lanes 8 14 Yaax Ki (YK), lanes Kitam Ki (KK); lanes 6 and 15 an individual from the Chelem White variant (CHW) also used as a reference; lanes 7 and 16 an individual from the Tropical subdeciduous forest (SF) also used as a reference.

5 January 1999] COLUNGA-GARCÍAMARÍNETAL. ISOZYMATIC VARIATION OF HENEQUÉN 119

6 120 AMERICAN JOURNAL OF BOTANY [Vol. 86 TABLE 2. Number of individuals per electrophenotype per isozyme, and total number of electrophenotypes per variant per isozyme. Variant codes are as described on Table 1 for the Yucatan Peninsula. ACP acid phosphatase, PRX peroxidase anodic, and MDH malate dehydrogenase. Variant Isozyme electrophenotype ACP Total A. angustifolia inside Yucatan 22 D DF SF CHW CHG CHY Henequén 2 SK 27 1 YK 25 1 KK 24 1 A. angustifolia outside Yucatan 7 Veracruz Veracruz Veracruz Oaxaca Oaxaca Oaxaca Edo. México Edo. México angustifolia plants whose ACP electrophenotype was not repeated in other plants under study. Most of the A. angustifolia plants growing outside the Yucatan Peninsula also showed unique electrophenotypes. With respect to PRX, the most abundant electrophenotype within A. angustifolia populations was also present in SK and KK. YK had an electrophenotype with intermediate abundance in A. angustifolia. Only two wild plants from the Yucatan Peninsula had unique electrophenotypes (one of D and one of CHY). All A. angustifolia plants growing outside of the Yucatan Peninsula had similar electrophenotypes to those growing in this region. For MDH, all plants growing in the Peninsula had the same electrophenotype, with the exception of KK plants, which have an electrophenotype also found in two plants from Veracruz and two from Oaxaca. Phylogenetic relations The parsimony analysis evaluated trees in the exhaustive search. One single most parsimonious tree (total length 112 steps) was found (Fig. 2). The mean SD length of the total set of possible trees was steps, g , which indicates that our total data set is more structured than a random data set (Hillis and Huelsenbeck, 1992). Two large clades may be observed in this tree: (1) one including a monophyletic group formed by the cordage variants SK and YK, which has variant SF as its sister group (the wild variant morphologically most similar to the cultigens), and (2) a clade formed by the remaining wild variants, as well as variant KK (the cultigen most similar to the wild populations). In a bootstrap with 1000 iterations, the only strongly supported clade was that formed by the cordage variants SK and YK. Results of these iterations are consistent with the morphometric analysis previously performed (Colunga-GarcíaMarín and May- Pat, 1997). The same results were obtained with the decay analysis. In all the trees up to ten steps longer than the most parsimonious tree, the clade formed by SK and YK was found. All other clades collapsed with just one step. Apomorphies found for this clade are the same characters that comprise the domestication syndromes found with the morphological variation analysis (Colunga- GarcíaMarín and May-Pat, 1996). DISCUSSION Isozymatic variation of henequén The very small fraction of the total isozymatic variation observed in the three varieties of henequén and the complete lack of variability in each one are results to be expected in a cultigen that is always propagated vegetatively and are consistent with our previous ethnobotanical (Colunga-GarcíaMarín and May-Pat, 1993) and morphological analyses (Colunga-GarcíaMarín, Estrada-Loera, and May-Pat, 1996; Colunga-GarcíaMarín and May-Pat, 1997). These studies indicated a clear promotion of only one cultivar since the middle of the last century by means of asexual propagation, as well as a limited morphological coefficient of variation in the cultivated varieties. A similar positive correlation between genetic variation indices and morphological variation coefficients was found by van Hintum and Elings (1991). These findings reveal the need to maintain a living germplasm collection of this cultigen in

7 January 1999] COLUNGA-GARCÍAMARÍNETAL. ISOZYMATIC VARIATION OF HENEQUÉN 121 TABLE 2. Extended. Isozyme electrophenotype PRX Total MDH 1 2 Total order to preserve the little remaining diversity. Another important aspect that would contribute to the preservation of henequén diversity would be to encourage a wider pattern of utilization, beyond cordage, to promote a wider range of varieties and in this way enhance genetic variation. In traditional Maya cultivation of henequén, selection and maintenance of diverse varieties were probably closely linked to a multipurpose use of this resource, as well as with its cultivation within a wider geographic range than at present (Colunga-GarcíaMarín and May- Pat, 1993). Some potential agroindustrial uses are the extraction of sapogenins, sugar for alcohol, and cellulose for paper and chemical cellulose. The extremely small fraction of total isozymatic variation observed within cultivated varieties is our main finding. Genetic variation losses are common in cultivars (Doebley, 1989), but in none of the reported cases has such reduction been so severe relative to the wild ancestors, both within and between cultivars. Loss of genetic variation in henequén is possibly, to a certain extent, a consequence of the domestication process begun in prehispanic times. As Doebley (1992) noted, although there is little hard evidence concerning the beginnings of plant domestication, it seems a reasonable assumption that the first farmers experimented with only a small fraction of the variation present within the progenitor species of today s crops. Further, as the domestication process proceeded these farmers probably selected only the best phenotypes. For these reasons one expects a significant loss of genetic variation, but competing forces such as introgression from wild relatives and selection against loss of fitness due to inbreeding have probably counterbalanced the expected loss of genetic variation. These competing Fig. 3. Phylogenetic relations of extant henequén variants and wild A. angustifolia populations from the Yucatan Peninsula, based on 66 morphological characters. The sole most parsimonious tree derived from an exhaustive search is shown. The tree was rooted by the midpoint technique. Tree length 112 steps; Consistency Index (CI) (CI when noninformative characters were excluded); Homoplasy Index (HI) (HI when noninformative characters were excluded); Retention Index (RI) Length is shown above each branch. Numerals in a circle above each branch represent the Decay Index. Numerals in parentheses below each branch represent bootstrap values (%) based on 1000 replicates. forces are usually favored by traditional agricultural practices (Colunga-GarcíaMarín, Hernández-Xolocotzi, and Castillo-Morales, 1986; Zizumbo-Villarreal, Hernández- Xolocotzi, and Cuanalo-de la C., 1988) as has been demonstrated for species such as the common bean (Escalante et al., 1994), maize, squash, and tomato (see reviews in Doebley, 1992), but they have not been promoted in henequén cultivation today. In other vegetatively propagated cultigens, such as Opuntia, the occasional seed cultivation has been observed in traditional farming (Colunga- GarcíaMarín, Hernández-Xolocotzi, and Castillo-Morales, 1986), a practice that may result in an increase of genetic variation. Under cultivation conditions of the present century the only remaining source of genetic variation for henequén would be somatic mutations. The lack of genetic variation of other clonal cultivars has already been reported. Results of DNA oligonucleo-

8 122 AMERICAN JOURNAL OF BOTANY [Vol. 86 tide fingerprinting of Malus, Prunus, and Rubus failed to detect genetic variation within clones (Nybom, Rogstad, and Schaal, 1990). Similar results with this technique were obtained in banana, in which the only differences were detected between an induced mutant and its progenitor (Kaemmer et al., 1992). However, even if variation was not observed within clones of these cultigens, the clones belong to species having numerous varieties. This is not the case in henequén, for which the number of cultivated varieties was reduced in the present century from seven to three, two of which are nearly extinct. During the henequén cultivation boom of 1915 in the state of Yucatan, henequén was cultivated over an area of up to ha (López and García de F., 1984), and it is possible all estimated 900 million plants had a single genotype, Sac Ki. Isozymatic variation of wild A. angustifolia populations Wild populations of A. angustifolia have high levels of isozyme diversity both within and among variants, a situation that contrasts strongly with that of henequén. This is normally expected in long life-cycle outcrossing perennials having wide geographic ranges (Hamrick, Godt, and Sherman-Broyles, 1992). These high levels of genetic variation emphasize the need to preserve wild germplasm in extant populations, based on which an eventual increase in diversity of cultigens may occur through a careful breeding program. The implementation of such program would require basic research on the reproductive biology of these species. The Coastal dunes area, where we found the largest electrophenotype diversity, is threatened by intense touristic development today, and there is strong human pressure to introduce permanent agriculture to the few remaining zones of Tropical deciduous and subdeciduous forest. A conservation area within the Yucatan Peninsula is urgently needed to preserve the in situ germplasm of wild A. angustifolia. In addition, a seed bank should be established to conserve genetic variation. Phylogenetic relationships Phylogenetic relationships within wild and cultivated variants growing in the Yucatan Peninsula inferred in this study are consistent with the ethnobotanical and morphological data previously reported (Colunga-GarcíaMarín and May-Pat, 1993; Colunga-GarcíaMarín, Estrada-Loera, and May- Pat, 1996; Colunga-GarcíaMarín and May-Pat, 1997). Two different domestication trends in present-day henequén variants seem to have occurred: one of SK and YK, selected for thicker and more abundant fibers and better suited for the cordage industry (whose sister group is SF); and another of KK, which is nearly extinct in Yucatan, selected for its softer fibers and better suited for textile applications. Within this last domestication trend, CHW was probably included at the beginning of this century, with its cultivation being abandoned later on. KK is a cultivated variant that morphologically groups together with wild ones, perhaps due to the artificial selection of wild characteristics (e.g., soft fibers). Even if we only obtained a most parsimonious tree, the low bootstrap support and the decay analysis values suggest that there is a lot of homoplasy in the data. The results obtained in this study agree with the predictions of Doebley (1989) in cases of wild-cultivated derivatives when isozymes are analyzed: (1) the cultigen falls within the variation range of the putative wild progenitor (in this study SK and YK electrophenotypes are the same as those of some plants of the Yucatan wild populations of A. angustifolia), (2) the cultigen has a subset of the allelic diversity found in the wild progenitor (in our case, only one electrophenotype per cultivar per isozyme), and (3) in addition to the cultigen having less genetic variation than the wild population, this genetic variation is distributed in a different way (henequén has more variation among variants than within variants, while this is the opposite in the case of A. angustifolia). The fact that KK has a different MDH electrophenotype with respect to all other Yucatan plants studied, and similar to some A. angustifolia plants from Oaxaca and Veracruz, seems to support the hypothesis of its recent introduction (Colunga-GarcíaMarín and May-Pat, 1997). LITERATURE CITED BARBA, R El henequén en Yucatán. Boletín de la Sociedad Agrícola Mexicana Nos. 19 and 20. BOLIO, A. J. A Manual práctico del henequén, su cultivo y explotación. Empresa Editorial Católica, S.A. Mérida, Yucatán. CASTORENA-SÁNCHEZ, I., R. M. ESCOBEDO, AND A. QUIROZ New cytotaxonomical determinants recognized in six taxa of Agave in the sections Rigidae and Sisalanae. Canadian Journal of Botany 69: COLUNGA-GARCÍAMARÍN, P Origen, variación y tendencias evolutivas del henequén (Agave fourcroydes Lem.). Ph.D. dissertation, Instituto de Ecología-UACPyP/CCH, Universidad Nacional Autónoma de México, México, D.F., J. COELLO-COELLO, L. ESPEJO-PENICHE, AND L. FUENTE-MORE- NO Agave studies in Yucatan, Mexico II. Nutritional value of the inflorescence peduncle and incipient domestication. Economic Botany 47: , E. ESTRADA-LOERA, AND F. MAY-PAT Patterns of morphological variation, diversity, and domestication of wild and cultivated populations of Agave in Yucatan, Mexico. American Journal of Botany 83: , E. HERNÁNDEZ-XOLOCOTZI, AND A. CASTILLO-MORALES Variación morfológica, manejo agrícola y grados de domesticación de Opuntia spp. en el Bajío Guanajuatense. Agrociencia 65: 7 49., AND F. MAY-PAT Agave studies in Yucatan, Mexico. I. Past and present germplasm diversity and uses. Economic Botany 47: , AND Morphological variation of henequén (Agave fourcroydes, Agavaceae) germplasm and its wild ancestor (A. angustifolia) under uniform growth conditions: diversity and domestication. American Journal of Botany 84: DE ECHÁNOVE, P. A Cuadro estadístico de Yucatán en Boletín de la Sociedad Mexicana de Geografía y Estadística: DOEBLEY, J Isozymic evidence and the evolution of crop plants. In D. E. Soltis and P. S. Soltis [eds.], Isozymes in plant biology, Advances in plant sciences series, vol. 4, Dioscorides Press, Portland, OR Molecular systematics and crop evolution. In P. S. Soltis, D. E. Soltis, and J. J. Doyle [eds.], Molecular systematics of plants, Chapman and Hall, New York, NY. ESCALANTE, A. M., G. COELLO, L. E. EGUIARTE, AND D. PIÑERO Genetic structure and mating systems in wild and cultivated populations of beans Phaseolus coccineus and P. vulgaris (Fabaceae). American Journal of Botany 81: ESPINOSA, J. D Manual de mayordomos de las fincas rústicas de Yucatán. Imprenta del autor. Mérida, Yucatán. GENTRY, H. S Agaves of continental North America. University of Arizona Press, Tucson, AZ. HAMRICK, J. L., M. J. W. GODT, AND S. L. SHERMAN-BROYLES

9 January 1999] COLUNGA-GARCÍAMARÍNETAL. ISOZYMATIC VARIATION OF HENEQUÉN 123 Factors influencing levels of genetic diversity in woody plant species. New Forest 6: HILLIS, D. M., AND J. P. HUELSENBECK Signal, noise and reliability in molecular phylogenetic analyses. Journal of Heredity 192: KAEMMER, D., R. AFZA, K.WEISING, G.KAHL, AND F. J. NOVAK Oligonucleotide and amplification fingerprinting of wild species and cultivars of banana (Musa spp.). Bio/technology 10: LÓPEZ H., R., AND A. GARCÍA DE F Manual de informaciónbásica de la región henequenera de Yucatán. Centro de Investigación Científica de Yucatán, A.C. Mérida, Yucatán, México. MADDISON, W. P., AND D. R. MADDISON Mac Clade, version 3.0. Sinauer, Sunderland, MA. MILLER, R. G., JR Simultaneous statistical inference. Springer- Verlag, New York, NY. NYBOM, H., S. H. ROGSTAD, AND B. A. SCHAAL Genetic variation detected by use of the M13 DNA fingerprint probe in Malus, Prunus, and Rubus (Rosaceae). Theoretical and Applied Genetics 79: ORELLANA, R., L. VILLERS, V. FRANCO, AND L. OJEDA Algunos aspectos ecológicos de los Agaves de la Península de Yucatán. In C. Cruz, L. del Castillo, M. Robert, and R. N. Ondarza [eds.], Biología y aprovechamiento integral del Henequén y otros Agaves, Centro de Investigación Científica de Yucatán, A.C. Mérida, Yucatán, México. REGIL, J. M., AND A. M. PEÓN Estadística de Yucatán. Boletín de la Sociedad Mexicana de Geografía y Estadística 3: SAS SAS/STAT user s guide, release 6.03 edition. SAS Institute, Cary, NC. SOLTIS, D., C. HAUFLER, D. DARROW, AND G. GASTONY Starch gel electrophoresis of ferns: a compilation of grinding buffers, gel and electrode buffers, and staining schedules. American Fern Journal 73: STEVENS, P. F Character states, morphological variation, and phylogenetic analysis: a review. Systematic Botany 16: STUBER, C. W., J. F. WENDEL, M.M.GOODMAN, AND J. S. C. SMITH Techniques and scoring procedures for starch gel electrophoresis of enzymes from maize (Zea mays L.). Technical Bulletin 286, North Carolina State University, Raleigh, NC. SWOFFORD, D. L PAUP: phylogenetic analysis using parsimony, version 3.1. Computer program distributed by the Illinois Natural History Survey, Champaign, IL. THORPE, R. S Coding morphometric characters for constructing distance Wagner networks. Evolution 38: VAN HINTUM, TH. J. L., AND A. ELINGS Assessment of glutenin and phenotypic diversity of Syrian durum wheat landraces in relation to their geographic origin. Euphytica 55: WENDEL, J., AND N. WEEDEN Visualization and interpretation of plant isozymes. In D. Soltis and P. Soltis [eds.], Isozymes in plant biology, Advances in plant sciences series, vol. 4, Dioscorides Press, Portland, OR. ZIZUMBO-VILLARREAL D., E. HERNÁNDEZ-XOLOCOTZI, AND H. CUANALO- DE LA C Estrategias agrícolas tradicionales para el aprovechamiento del agua de lluvia durante el temporal: el caso de Yuriria, Guanajuato, Méx. Agrociencia 71: