Print version ISSN 0187-7380
Rev. fitotec. mex vol.34 n.3 Chapingo Jul./Sep. 2011
Genetic diversity of Agave cupreata Trel. & Berger. Considerations for its conservation
Diversidad genética de Agave cupreata Trel. & Berger. Consideraciones para su conservación
Alejandro MartínezPalacios1*, Juan M. GómezSierra1, Cuauhtémoc SáenzRomero1, Nidia PérezNasser2 y Nahúm SánchezVargas1
1 Instituto de Investigaciones Agropecuarias y Forestales, Universidad Michoacana de San Nicolás de Hidalgo. Km 9.5 Carretera MoreliaZinapécuaro. 58880, Tarímbaro, Michoacán, México. Phone +(52)(443) 3340475 ext. 119, fax ext. 200. * Corresponding author (email@example.com)
2 Centro de Investigaciones en Ecosistemas, Universidad Nacional Autónoma de México. Antigua Carretera a Pátzcuaro No. 8701, Col ExHacienda de San Miguel de la Huerta. 58190, Morelia, Michoacán.
Recibido: 07 de Julio del 2010.
Aceptado: 15 de Junio del 2011.
Agave cupreata Trel. & Berger is an endemic plant naturally distributed in the Balsas Depression, a semiarid region in the states of Guerrero and Michoacán in Southwestern México. Their populations are heavily decimated because mature individuals just before their single life flowering period are harvested to produce mescal, an alcoholic beverage. The genetic variation among and within 12 natural populations was examined for nine isozyme loci. Results indicate high average proportion of polymorphic loci (93 %) and expected heterozygosity (H = 0.467), with an excess of observed heterozygotes in relation to HardyWeinberg expectations (Ho = 0.521, F = 0.1179). These results represent the largest heterozygosity reported for Agave species endemic to México. There is also a statistically significant genetic differentiation among populations (FST = 0.042). An UPGMA dendrogram reveals the absence of a geographic pattern, as confirmed by a Mantel test (r = 0.110, P = 0.769), which did not show significant isolation by distance. Estimated minimum viable effective population size was very large (Ne =16,165), larger than in any other known natural population. To protect the natural genetic variation, it is suggested to design and manage A. cupreata natural populations as forest genetic resource conservation units (FGRCUs) using realistic and modest Ne sizes, perhaps between 500 and 5000 plants, ideally with intermediate plantations that could serve as pollinator corridors. Commercial plantations and ex situ FGRCUs need to be established to gradually develop a sustainable management, perhaps at higher altitudes than current locations, as a management measure for adaptation to the climatic change.
Key words: Agavaceae, population genetics, minimum viable effective population size.
Agave cupreata es una planta endémica que se distribuye en la Depresión del Balsas, una región semiárida de los Estados de Guerrero y Michoacán, al suroccidente de México. Sus poblaciones están muy diezmadas debido a que los individuos maduros se cosechan justo antes de su floración, que ocurre sólo una vez en su vida, para producir mezcal, una bebida alcohólica. La variación genética entre y dentro de 12 poblaciones naturales fue examinada en nueve loci isoenzimáticos. Los resultados indican un promedio muy elevado de proporción de loci polimórficos (93%) y de heterocigosidad esperada (He = 0.467), con un exceso de heterocigotos observados en relación con lo esperado bajo el equilibrio de HardyWeinberg (Ho = 0.521, FIS = 0.1179). Estos resultados de heterocigosidad son la más altos reportados en especies de Agave endémicas de México. También hay una diferenciación genética entre poblaciones estadísticamente significativa (FST = 0.042). Un dendrograma UPGMA revela la ausencia de un patrón geográfico, confirmada por una prueba de Mantel (r = 0.110, P = 0.769), que demostró que no hay un significativo aislamiento por distancia. El tamaño mínimo efectivo viable de la población fue muy grande (Ne = 16,165), más grandes que en cualquier otra población natural conocida. Para proteger la variación genética natural de A. cupreata se sugiere diseñar y gestionar la conservación de poblaciones naturales como unidades de conservación de recursos genéticos forestales (UCRGFs) con tamaños de Ne realistas y modestos, tal vez entre 500 y 5000 plantas, idealmente con plantaciones intermedias que podrían servir como corredores de polinizadores. Se necesita establecer plantaciones comerciales y UCRGF ex situ para desarrollar gradualmente un manejo sustentable, tal vez a altitudes mayores que a la que se encuentran las localidades actuales, como una medida de manejo para adaptarse al cambio climático.
Palabras clave: Agavaceae, genética de poblaciones, tamaño efectivo mínimo de población viable.
Agave cupreata Trel. & Berger is an endemic plant of the Agavaceae family (Crenatea group, subgenus Agave), now classified as a subfamily of the expanded Asparagalean families (Chase et al., 2004). It is naturally distributed in the interior uphill slopes (1220 to 1850 m of altitude) of the Balsas Depression, a semiarid region in the states of Guerrero and Michoacán in Southwestern México (Gentry, 1982). Many of the natural populations are heavily decimated because mature whole individuals are harvested just before flowering, when the sugar concentration is at its maximum, to produce mescal in local artisan distilleries, a traditional liquor similar to the well known Mexican tequila (ColungaGarcíaMarín and ZizumboVillarreal, 2007; ZizumboVillarreal and ColungaGarcíaMarín, 2007). Harvesting prevents the species reproduction because A. cupreata is a semelparous plant (it can reproduce sexually only once in their life time) and it does not reproduce vegetatively. The effect of harvesting in the population demography is aggravated because plants need between 7 and 15 years to reach their sexual maturity (Illsley et al., 2007). Apparently, the overexploitation of the natural populations is more severe in the state of Michoacán than that in Guerrero. There are a few examples of well organized Indian communities to conduct some management of the A. cupreata natural populations in the region of Chilapa, CentralEast part of the state of Guerrero, with specific rules that contribute to a more sustainable management, such as harvesting no more than 80 % of the individuals at flowering stage, reintroducing seedlings propagated in nurseries and setting aside populations designated as reserves (IllsleyGranich C. 2004; Personal Comm.)1.
Local producers of mescal, mainly from the state of Michoacán, recently recognized that their use of natural populations of A. cupreata is not sustainable, and gradually began to organize themselves in mescal producer associations that have established approximately 800 ha of plantations (SEDRU, 2009) from seedlings produced in rustic greenhouses without control of seed source. Although such actions are an encouraging start, a more sophisticated management program is required which should combine a sustainable harvesting with a biological conservation program to protect the genetic variation of A. cupreata (Eguiarte and Souza, 2007). For reaching this goal, the patterning of genetic variation of natural populations needs to be determined.
Genetic diversity among and within populations of several species of Agave endemic to México has been studied using molecular markers. In general, high levels of expected heterozygosity (He ) were observed, and high genetic differentiation among populations (FST or GST) in one case, and medium to low in most cases, has been documented: He = 0.335 and FST = 0.236 for Agave victoriaereginae T. Moore (MartínezPalacios et al., 1999); He = 0.394 and θ= 0.083 for Agave lechuguilla (SilvaMontellano and Eguiarte, 2003b); He = 0.12 to He = 0.29 and GST = 0.11±0.02 for the complex Agave deserti Engelm., A. cerulata Trel. and A. subsimplex Trel. (NavarroQuezada et al., 2003); and He= 0.261 for Agave angustifolia Haw. (BarrazaMorales et al., 2006). However, to our present knowledge there are no previously published studies of the genetic variability of A. cupreata.
The objectives of the present study are: a) estimate levels of variation and genetic structure among and within natural populations of Agave cupreata using isozyme markers, and b) suggest measures of conservation of the genetic diversity of A. cupreata.
MATERIALS AND METHODS
Studied populations and material collection
Leaves of 37 randomly selected mature healthy individuals of each of 12 natural populations (8 from the state of Michoacán and 4 from Guerrero) were collected along the entire distribution range of A. cupreata (Table 1). The groups of individual plants represented by the samples are termed populations while the location where a population was collected is called provenance. Leaf samples were transported to the lab in coolers (6 °C), and 3 cm x 3 cm samples were stored at 70 °C until processed.
Starch gel electrophoresis and staining protocols were conducted following standard methods (Conkle et al., 1982; Cheliak and Pitel, 1984), with modifications as described by MartínezPalacios et al. (1999). Sections of fresh, young leaves were crushed into 0.5 mL of the grinding buffer [3:1 (v:v) mixture of buffer YO from Yeh and O'Malley (1980) and VeglI from Pitel and Cheliak (1984)]. The extracts were adsorbed in Whatman No. 17 paper wicks that were applied later to starch gels 11 %. Two systems of electrode and gel buffers, and seven enzymes with good resolution were used. A total of 9 loci were analyzed. System LiOH 8 (Soltis et al., 1983; MartínezPalacios et al., 1999) resolved glutamate oxaloacetate transaminase (Got1), malic enzymes (Me2), menadione reductase (Mnr1), phosphoglucoisomerase (Pgi1 and Pgi2). System C (Stuber et al., 1988) resolved acid phosphatase (Acph1), glutamate dehydrogenase (Gdh1), phosphoglucomutase (Pgm1 and Pgm2). Gels were electrophoresed at 60 mA for 7 to 8 h.
Enzymatic systems with more than one locus were numbered according to their mobility relative to the alleles of standard individuals present in all gels and systems. Allelic variants were numbered sequentially from anode to cathode. The numbers of loci and alleles genetically controlling the enzyme activity were inferred from the observed banding patterns and from data on quaternary structure (Wendel and Weeden, 1989). Initially we essayed other loci (Dial, Dia2, Estl and Est2) that worked with Agave victoriaereginae (MartínezPalacios et al., 1999), but we discarded them because they showed poor resolution with A. cupreata extracts. However, since we included Gdhl, Pgml and Pgm2, our final number of examined loci (9) was not much different than that used by A. victoriaereginae (10 loci).
We estimated the number of polymorphic loci (P, criterion of 95 %), average number of alleles per locus (A), observed (Ho) and expected (He) heterozygosity under HardyWeinberg equilibrium Chi2 test (Snedecor and Cochran, 1967), exact test of Haldane (1954), Wright's F statistics (Wright, 1965) and their 95 % confidence intervals (using 10 000 iterations), genetic distances among populations (Nei, 1978), and a dendrogram was constructed through the unweighted pair group method with arithmetic mean (UPGMA) using the tools of the population genetics program (TFPGA) (Miller, 1997). A Mantel test (Manly, 1987) was performed between genetic and geographic distances. Minimum viable effective population size (Ne) was estimated based on the average He using the linear model fitted by regressing Ne values against He values, as used by Millar and Libby (1991) in several conifer species. Ne values estimated by Millar and Libby (1991) assumed a mutation rate of 105 and are designed to maintain given He values. The regression model proposed by ViverosViveros et al. (2010) was used:
Ne = 984.58 + (36723 He)
RESULTS AND DISCUSSION
Variation and genetic structure
High polymorphism was found, with an average among populations of 93 % polymorphic loci and 2.9 alleles per locus. The heterozygosity was also very large, with an average expected value of He = 0.467, and an even larger value of observed heterozygosity of Ho = 0.521 (Table 1). These values are the largest ever reported for an Agave species endemic to México (MartínezPalacios et al., 1999; NavarroQuezada et al., 2003; SilvaMontellano and Eguiarte, 2003a; Demey et al., 2004). This result should be taken with caution because we used only 9 loci; there is a possibility that by examining a larger number of loci the estimated heterozygosity could be lower, if by chance our sample could be biassed. Haldane test showed that 5 out of 9 loci had significant excess of heterozygotes, in comparison to HardyWeinberg expectations (P < 0.03).
Although moderate, there is a significant genetic differentiation among populations with FST = 0.042 (95 % confidence interval from 0.0613 to 0.0281) (Table 2). Such differentiation is lower than the values reported for other Agave species: FST = 0.08 for A. subsimplex (NavarroQuezada et al., 2003); FST = 0.083 for A. lechuguilla (SilvaMontellano and Eguiarte, 2003a); FST = 0.135 for A. deserti; and FST = 0.236 for A. victoriaereginae (MartínezPalacios et al., 1999).
It is likely that the moderate genetic differentiation is produced by the counteracting effect of a relevant gene flow, estimated as Nm=5.68 (Table 2), which might be due to the action of efficient pollinators such as bees, birds and bats (Gentry, 1982; Eguiarte et al., 2000).
Average genetic distance between pairs of populations was relatively low (D = 0.04). The largest genetic distance (D = 0.09) was found between the Escobas and Llanitos populations, both from Michoacán, while the lowest genetic distance (D = 0.008) was between Tixtla, Guerrero and Escalera, Michoacán (Table 3). The UPGMA dendrogram does not follow a geographic pattern since geographically distant populations clustered together (for example: Limón, Michoacán grouped with Mesones, Guerrero, both separated by 322 km), while some geographically close populations did not clustered together (Figure 1). This lack of association between genetic distances and geographic locations was confirmed by the Mantel test (r = 0.110, P = 0.769), which showed no significant isolation by distance.
Despite the heavy exploitation pressure to which A. cupreata populations are subjected, there is no evidence of inbreeding within populations: FIS = 0.1179 (significantly different from zero, 95 % confidence interval from 0.0072 to 0.2054), a value that indicates a significant excess of heterozygotes. This can be related to the mating system of A. cupreata that only reproduces sexually, their flowers display herkogamy (spatial separation of anthers and stamens) and dichogamy (temporal separation of male and female reproductive organs), and specifically protandry (earlier maturation of anthers than that of pistils); these characteristics are favoring outcrossing and apparently making self pollination nearly impossible (Illsley et al., 2007). Thus, this outcrossing system apparently plays a relevant role to prevent the accumulation of inbreeding, even when the population sizes are reduced. However, interspecific crossing is seemingly possible, since seedlings that appear to be hybrids of A. cupreata with A. attenuata and A. inaequidens have been observed in established plantations of A. cupreata (MartínezPalacios A. 2009; Personal comm.)2.
Estimated minimum viable effective population size was Ne =16 165. This value is very large, as a result of the very large average heterozygosity value. In general, for conservation of genetic variability the suggested Ne sizes range between 500 (Franklin, 1980; Frankham et al., 2002) and 5000 individuals (Lande 1995; Frankham et al., 2002). In our case it would be more realistic to use much more modest Ne sizes than the estimated value (16 165), considering that natural populations of A. cupreata have substantially lower census sizes, perhaps a maximum of Ne ≈ 500 (reproductive individuals flowering the same year). Thus, in order to designate A. cupreata natural populations as forest genetic resource conservation units (sensu SaenzRomero et al., 2003, based on Ledig , Millar and Libby ), the largest populations would need to be located and protected.
Considerations for conservation
A sustainable management program of A. cupreata that allows continuous harvesting of mature individuals for mescal distillation needs to gradually satisfy the demand of this beverage by relying more on commercial plantations, rather than in natural populations. In parallel, it would be needed to protect natural populations and to establish additional ones to be managed as forest genetic resource conservation units (FGRCUs). Provided that the estimated minimum viable effective population size (Ne) for such FGRCUs is larger than the actual largest natural populations, it would be required to establish plantations between large natural populations, to serve as corridors for pollinators, in order to maintain the genetic flow among populations and also for increasing Ne, counterbalancing in that way the loss of individuals that could have occurred due to past harvesting. On this regard, although the recent establishment of commercial plantations of A. cupreata aided by subsidies from the federal, state and municipality governments seems to be a measure in the correct direction, it might be needed to establish additional plantations with the main goal of conservation, and not only for commercial harvesting.
The predicted decoupling between genotypes and climate due to factual climatic change could be diminished if new commercial plantations are established at altitudes higher than those of present populations, in order to match genotypes fitness to predicted future climates. Although it is generally suggested to plant at an elevation 300 m higher than the seed source, to match the climate predicted for 2030 (SaenzRomero et al., 2010), it is recommended to make decisions based on estimations for contemporary and future climate for specific locations, using available resources, including web based tools (e.g. Crookstone, 2011).
Agave cupreata populations have a very large expected heterozygosity (He = 0.467), with an excess of observed heterozygotes (Ho = 0.521, FIS = 0.1179). The genetic differentiation among populations is relatively low (FST = 0.042), yet statistically significant. There is not significant isolation by distance. Estimated minimum viable effective population size was very large (Ne =16 165), due to the large heterozygosity values.
Financial support was provided to Alejandro MartínezPalacios by the joint research fund for basic research between the Mexican National Council of Science and Technology (CONACYT) and the Mexican Ministery of Education (SEP, project 2004P47777Z), the Scientific Research Council (CIC) of the Universidad Michoacana de San Nicolás de Hidalgo (UMSNH, project 5.6) and the State of Michoacán Council of Science and Technology (COECYTMichoacán, Project CB07021228). Three anonymous reviewers contributed to significantly improve the manuscript.
BarrazaMorales A,FL SánchezTeyer, M Robert, M Esqueda, A Gardea (2006) Variabilidad genética en Agave angustifolia Haw. de la sierra Sonorense, México, determinada con marcadores AFLP. Rev. Fitotec. Mex. 29:18. [ Links ]
Chase M W, J L Reveal, M F Fay (2009) A subfamilial classification for the expanded asparagalean families Amaryllidaceae, Asparagaceae and Xanthorrhoeaceae. Bot. J. Linnean Soc. 161:132136. [ Links ]
Cheliak W M, J A Pitel (1984) Electrophoretic identification of clones in trembling aspen. Can. J. For. Res. 14:740743. [ Links ]
ColungaGarcíaMarín P, D ZizumboVillarreal (2007) El tequila y otros mezcales del centrooccidente de México: domesticación, diversidad y conservación de germoplasma. In: En lo Ancestral hay Futuro: del Tequila, los Mezcales y otros Agaves. P ColungaGarcíaMarín, L Eguiarte, A Larqué, D ZizumboVillarreal (eds). CICYCONACYTCONABIOSEMARNATINE. México D.F. pp:113131. [ Links ]
Conkle M T, P D Hodgskiss, L B Nunnally, S C Hunter (1982) Starch gel electrophoresis of conifer seeds: a laboratory manual. General Technical Report PSW64. Pacific Southwest Forest and Range xperiment Station. United States Forest Service. Berkeley, California. 18 p. [ Links ]
Crookstone N L (2011) Custom Data Requests. In: Research on Forest Climate Change: Potential Effects of Global Warming on Forests and Plant Climate Relationships in Western North America and Mexico. Moscow Forestry Sciences Laboratory, Rocky Mountain Research Station, USDA Forest Service. Available at: http://forest.moscowfsl.wsu.edu/climate/customData/ (May 2011). [ Links ]
Demey J R, E Gamez, S Molina, D Infante (2004) Comparative study of the discriminating capacity of AFLP and ISTR markers for genetic analysis of Agave fourcroydes. Plant Mol. Biol. Rep. 22:2935. [ Links ]
Eguiarte L E, V Souza, A SilvaMontellano (2000) Evolución de la familia Agavaceae: filogenia, biología reproductiva y genética de poblaciones. Bol. Soc. Bot. Méx. 66:131150. [ Links ]
Eguiarte L E, V Souza (2007) Historia natural del Agave y sus parientes: Evolución y ecología. In: En lo Ancestral hay Futuro: del Tequila, los Mezcales y otros Agaves. P ColungaGarcíaMarín, L Eguiarte, A Larqué, D ZizumboVillarreal (eds). CICYCONACYTCONABIOSEMARNATINE. México D.F. pp:321. [ Links ]
Franklin I R (1980) Evolutionary change in small populations. In: Conservation Biology: an EvolutionaryEcological Perspective. M E Soulé, B A Wilcox (eds). Sinauer. Sunderland, MA, U.S.A. pp:135150. [ Links ]
Frankham R, J D Ballou, DA Briscoe (2002) Introduction to Conservation Genetics. Cambridge University Press. Cambridge, UK. 617 p. [ Links ]
Gentry S H (1982) Agaves of Continental North America. Univ. Arizona Press. Tucson, AZ. 669 p. [ Links ]
Haldane J B S (1954) An exact test for randomness of mating. J. Genetics 52:631635. [ Links ]
Illsley C, E Vega, I Pisanty, A Tlacotempa, P García, P Morales, G Rivera, J García, V Jiménez, F Castro, M Calzada (2007) Maguey papalote: hacia el manejo campesino de un recurso colectivo en el trópico seco de Guerrero, México. In: En lo Ancestral hay Futuro: del Tequila, los Mezcales y otros Agaves. P ColungaGarcíaMarín, L Eguiarte, A Larqué, D ZizumboVillarreal (eds). CICYCONACYTCONABIOSEMARNATINE. México D.F. pp:319338. [ Links ]
Lande R (1995) Mutation and conservation. Conserv. Biol. 9:782791. [ Links ]
Ledig F T (1988) The conservation of diversity in forest trees. Why and how should genes be conserved? Bioscience 38:471479. [ Links ]
Manly B F J (1987) The Statistics of Natural Selection on Animal Populations. Chapman and Hall. New York, N.Y., USA. 484 p. [ Links ]
MartínezPalacios A, L E Eguiarte, G R Furnier (1999) Genetic diversity of endangered endemic Agave victoriaereginae (Agavaceae) in Chihuahuan desert. Am. J. Bot. 86:10931098. [ Links ]
Millar C I, W J Libby (1991) Strategies for conserving clinal, ecotypic, and disjunct population diversity in widespread species. In: Genetics and Conservation of Rare Plants. D A Falk, K E Holsinger (eds). Oxford University Press. New York. pp:149170. [ Links ]
Miller M P (1997) Tools for Population Genetics Analyses (TFPGA) 1.3: A windows program for the analysis of allozyme and molecular population genetic data. Computer software distributed by author. Available at: http://www.marksgeneticsoftware.net/_vti_bin/shtml.exe/tfpga.htm (May 2010). [ Links ]
NavarroQuezada A, R GonzálezChauvet, F MolinaFreaner, L E Eguiarte (2003) Genetic differentiation in the Agave deserti (Agavaceae) complex of the Sonoran desert. Heredity 90:220227. [ Links ]
Nei M (1978) The theory of genetic distance and evolution of human races. J. Human Genet. 23:341369. [ Links ]
Pitel JA,WM Cheliak (1984) Effect of extraction buffers on characterization of isoenzymes from vegetative tissues of five conifers species: a user's manual. Information Report PIX34, Petawawa National Forestry Institute. Canadian Forestry Service. Chalk River, Ontario, Canada. 64 p. [ Links ]
SáenzRomero C, A E Snively, R LindigCisneros (2003) Conservation and restoration of pine forest genetic resources in México. Silvae Genet. 52:233237. [ Links ]
SáenzRomero C, G E Rehfeldt, N L Crookston, P Duval, R StAmant, J Beaulieu, B A Richardson (2010) Spline models of contemporary, 2030, 2060 and 2090 climates for México and their use in understanding climatechange impacts on the vegetation. Climatic Change 102:595623. [ Links ]
SEDRU, Secretaría de Desarrollo Rural (2009) Exigen productores de mezcal en Michoacán la denominación de origen de su producto. Secretaría de Desarrollo Rural. Gobierno del Estado de Michoacán. Morelia, Mich. Available at: http://www.michoacan.gob.mx/SEDRU/Exigen_productores_de_mezcal_en_Michoacan_la_ denominacion_de_origen_de_su_producto (July 2009). [ Links ]
SilvaMontellano A, L E Eguiarte (2003a) Geographic patterns in the reproductive ecology of Agave lechuguilla (Agavaceae) in the Chihuahuan desert. I. Floral characteristics, visitors and fecundity. Am. J. Bot. 90:377387. [ Links ]
SilvaMontellano A, L E Eguiarte (2003b) Geographic patterns in the reproductive ecology of Agave lechuguilla (Agavaceae) in the Chihuahuan desert. II. Genetic variation, differentiation and inbreeding estimates. Am. J. Bot. 90:700706. [ Links ]
Snedecor G W, W G Cochran (1967) Statistical Methods, 6th ed. The Iowa State University Press. Ames, IA. 593 p. [ Links ]
Soltis D E, C H Haufler, D C Darrow, G J Gastrony (1983) Starch gel electrophoresis of ferns: a compilation of grinding buffers, and staining schedules. Amer. Fern J. 73:927. [ Links ]
Stuber C W, J F Wendel, M M Goodman, J S C Smith (1988) Techniques and scoring procedures for starch gel electrophoresis of enzymes from maize (Zea mays L.). Technical Bulletin 286. North Carolina State University. Raleigh, N.C., EUA. 87 p. [ Links ]
ViverosViveros H, B L TapiaOlivares, C SáenzRomero, J J VargasHernández, J LópezUpton, A SantacruzVarela, G RamírezValverde (2010) Isoenzymatic variation of Pinus hartwegii Lindl. along an altitudinal gradient in Michoacán, México. Agrociencia 44:723733. [ Links ]
Wendel J F, N F Weeden (1989) Visualization and interpretation of plant isozymes. In: Isozymes in Plant Biology. D E Soltis, P S Soltis (eds). Chapman and Hall. London, UK. pp:545. [ Links ]
Wright S (1965) The interpretation of population structure by Fstatistics with special regard to systems of mating. Evolution 19:395420. [ Links ]
Yeh F C, D M O'Malley (1980) Enzyme variation in natural populations of Douglasfir, (Pseudotsuga menziessii (Mirb.) Franco) from British Columbia. I. Genetic variation patterns in coastal populations. Silvae Genet. 29:8392. [ Links ]
ZizumboVillarreal D, P ColungaGarcíaMarín (2007) La introducción de la destilación y el origen de los mezcales en el occidente de México. In: En lo Ancestral hay Futuro: del Tequila, los Mezcales y otros Agaves. P ColungaGarcíaMarín, L Eguiarte, A Larqué, D ZizumboVillarreal (eds). CICYCONACYTCONABIOSEMARNATINE. México D.F., pp:85112. [ Links ]
1 IIlsleyGranich C (2004) Manejo campesino sustentable del maguey papalote de Chilapa. Informe final de proyecto. Grupo de Estudios Ambientales AC y Comisión Nacional para el Conocimiento y Uso de la Biodiversidad, México D.F. 43 p.
2 MartínezPalacios A, S Chávez Mendoza, M Gómez Sierra, Ma Sierra Yxta, R Cárdenas Navarro (2009) Management and Conservation of Agave cupreata (Agavaceae), (Resumenes en extenso en CD). First International Conference on Sustainable Cities. 26 al 29 de Octubre de 2009, Morelia Michoacán, México. pp:111.