Open access peer-reviewed chapter

The Endangered Species Dioon edule in the Sierra Madre Oriental in San Luis Potosí: Demography and Genetic Diversity

Written By

Gabriel Rubio-Méndez, Alberto Prado, Jacqueline C. Bede, José Arturo De-Nova, Joel Flores, Juan Antonio Reyes-Agüero and Laura Yáñez-Espinosa

Submitted: September 17th, 2020 Reviewed: February 3rd, 2021 Published: March 9th, 2021

DOI: 10.5772/intechopen.96372

Chapter metrics overview

266 Chapter Downloads

View Full Metrics


The distribution of the endangered species Dioon edule is in populations scattered throughout the Sierra Madre Oriental in San Luis Potosí, Mexico. Its habitat is tropical dry forests at lower elevations and oak forests at higher elevations, mainly disturbed by anthropic activities. We determined and analyzed nine populations’ demographic structure and explored the genetic diversity of five using SSR markers. The population density averaged 2050 individuals ha-1 and have an aggregated distribution pattern. Differences in the plants’ size among localities are due to site quality, based on their adaptation capacity and response to climate and soil traits. Most populations have the highest mortality in the early stages of life, with a low mortality rate for those who survive this bottleneck. Two populations have a relatively constant mortality rate, attributed to disturbance of the habitat. The populations show low genetic diversity and an excess of homozygotes. Their similarity is probably related to the formation of natural corridors favoring connectivity between populations. The deterioration and fragmentation of the habitat have severe effects on the populations’ viability, like reducing gene flow, which has led to inbreeding and genetic drift.


  • Zamiaceae
  • cycads
  • gymnosperms
  • ecotone
  • dioecious
  • habitat disturbance
  • inbreeding
  • Mexico
  • population
  • gene flow

1. Introduction

The relevance of the species Dioon edule (Zamiaceae) is evolutive because it belongs to a group of plants whose origin can be traced back to 250 million years ago, allowing us to understand the adaptations developed to survive the present day. It is also ecologically relevant as it is scattered in a wide range of ecosystems, from coastal dunes in the Gulf of Mexico’s coastal plain to pine forest in the Sierra Madre Oriental, and its interactions with symbiotic micro-organisms, pollinating insects, and predator-dispersing wild and domestic fauna. The region’s cultural relevance is defined by the Xi’iuy ethnic group in San Luis Potosí, collecting the seeds for food and the leaves for ceremonies [1].

This species is distributed in the Sierra Madre Oriental (SMO), within the geographic-cultural region known as La Huasteca, including the south of Tamaulipas, east of San Luis Potosi, Hidalgo, Queretaro, and Veracruz. The elevation goes from sea level to 1525 meters above sea level (m.a.s.l.) [2, 3]. It has stems up to 6 m high in some occasions semi-prostrated, presenting reproductive events with an interval of 10 to 52 years in female individuals and 2.8 to 8.8 in male individuals [4]. Its strobilus is pollinated only by specialist insects (beetles) associated with them [5].

The species is currently classified as near threatened (NT) by International Union for Conservation of Nature and endangered by the NOM-059-SEMARNAT-2010 in Mexico. This category is mainly because of land-use change to agriculture and livestock, causing habitat fragmentation that adversely affects population dynamics and a possible reduction in gene flow among populations [4, 6].

It is necessary to integrate information from different areas, mainly demography and population genetics, to develop and implement effective conservation strategies [6] to build general conclusions from the observations derived from each of these areas.

Demographic studies provide the basis for knowledge of populations and allow us to assess populations’ viability, identify factors that reduce the population, and predict the persistence of small, isolated populations and evaluate alternative conservation scenarios [7]. Population genetics allows us to describe the genetic composition of natural populations and predict their changes in response to various evolutionary forces operating on them, allowing us to identify some risks of a genetic nature that affect the persistence of species such as habitat fragmentation and loss of adaptive potential associated with the decline of genetic diversity and inbreeding [8, 9]. Diverse studies related to population genetics and genomics of cycads have been published based on molecular markers such as Aloenzymes, RADseq, chloroplast markers, ISSR [4, 10, 11, 12, 13, 14, 15].

The objective of this research was to study the demographic and genetic status of Dioon edule Lindl. (Zamiaceae) in representative natural populations of the Sierra Madre Oriental in the state of San Luis Potosi, analyzing their population structure concerning their density, age structure, sex ratio, as well as genetic diversity.


2. Materials and methods

2.1 Study area

The Sierra Madre Oriental is fundamentally a group of minor mountain ranges formed by marine sedimentary rocks, mainly limestone, sandstone, and shales. The maximum elevations are up to 3180 m.a.s.l., with small intermontane valleys limited by high mountain ranges with steep slopes (Figure 1). In this region, there are soils of alluvial origin formed in the great plains with sediments mostly coming from limestone and lutites and residual and colluvial origin in the high and lower parts of the mountain ranges constituted by the same type of rock. In general, the soil is medium texture, shallow, and rich in organic matter and nutrients depending on the climate and vegetation [16]. The area’s climate is semi-warm, characterized by average annual temperatures above 18 °C and 600 mm of rainfall, mostly during the summer [16].

Figure 1.

Study area in the Sierra Madre oriental of San Luis Potosí, Mexico.

In the state of San Luis Potosi, three important vegetation regions for cycads stand out [16]: Tamasopo with tropical rain forest and oak forest (Quercus spp.); Aquismón with tropical rain forest; and, Rayón and Santa Catarina with oak forest, and piedmont scrub.

2.2 Demography

The study of demography was conducted in nine populations of Dioon edule distributed in the sub-province of Carso Huasteco (Figure 1) that covers 15.52% of the total surface of San Luis Potosi. It is formed almost exclusively by mountain ranges with soils of residual origin, shallow (less than 10 cm), with abundant rocky outcrops and Lithic phase. We do not provide detailed information on populations’ locations to avoid possible looting, so each one has a key with two letters.

The sampling consists of modifying the line interception method [17] and used by [18] in transects of 30 m long and 10 m wide. Rectangular units provide the advantage of evaluating variables by walking in a straight line without moving sideways, starting at one end of the line, and considering only those plants within the plot [19]. Due to the harsh characteristics of the area with rocky slopes of up to 70°, we traced two linear transects 100 m long by 2 m wide, parallel to the slope, locating them where the most significant number of individuals were concentrated, considering that D. edule populations present an aggregate distribution [18]. Each transect was divided into four 25 m long subunits to carry out the plants’ counts and measurements.

The different stages of the cycads’ life cycle were defined according to the number of leaves, the number of leaflets and stem height [7] because it has not yet been possible to define the age of the cycads due to their growth characteristics [18]. Therefore, individuals were classified as follows: seedlings, individuals that present one or more leaves without these forming a crown; juveniles, individuals that present one or more crowns, but with a non-visible stem; adults, individuals that present one or more crowns and have a visible stem. Adults were classified as male, female, or non-reproductive according to the presence of reproductive structures.

Once the cycads’ life cycle stages were defined, each stage’s survival probabilities calculation as the proportion of individuals that survived to enter the next category, relative to the initial number of seeds (lx = nx/n0) [20]. A survival curve describes the survival pattern of individuals in a population over time, plotting the lx on the y-axis, in logarithmic scale, with life cycle stages on the x-axis There are three basic survival curves, the survival curve type I represents a population in which most of the organisms die in the older age categories. The type II curve represents a population in which the mortality rate is almost constant, and type III represents a population in which most individuals die during the first stages of the life cycle, and very few reach the intermediate and late age categories [20].

The stem’s diameter and the plant’s total height were measured, the total number of leaves was counted, and the length of those on the last crown was measured.

Considering that the species D. edule is dioecious and produce only one cone per reproductive event, we evaluated the number of plants with male and female cones present as a measure of the potential reproductive capacity of the population per reproductive event [21].

The data on population density and potential reproductive capacity was analyzed using a one-way Analysis of Variance (ANOVA) to determine the statistical difference among the nine populations studied (localities) with two transects in each (n = 18).

Data on plant size and number of leaves of individuals in reproductive and non-reproductive states were analyzed through a nested ANOVA (random-effects model II) to determine the statistical difference among localities and categories of plants’ reproductive states in all localities and each location. When statistically significant differences between means occurred, a Tukey multiple comparison test was applied (p < 0.05). All analyzes were performed with XLSTAT software (Addinsoft Inc., Paris, France).

2.3 Genetic diversity

The genetic diversity was studied in five populations (Figure 1), covering an altitudinal gradient from 388 to 1050 m.a.s.l. In each population, leaf tissue was collected from 21 individuals (seven from each age category). The DNA extraction was done with the DNeasy Plant minikit® Kit (Quiagen, Valencia, California, USA), following the kit’s protocol and using approximately 100 mg of leaf tissue from which the cuticle was previously detached. A total of seven microsatellite markers designed by [22] (Table 1) labeled M-13 (−21pb) to obtain fluorescent products [23, 24], and PCR products were separated and genotyped with the LI-COR 4300 DNA analyzer (LI-COR Biosciences, Lincoln, Nebraska, USA) in a 6.5% polyacrylamide gel (LI-COR Biosciences, Lincoln, Nebraska, USA), at a wavelength of 700 nm.

LocusPrimer sequences (5′–3′)RepeatNSize

Table 1.

SSR primers used. N, number of alleles; size, allele size range (bp).

It was calculated the average number of alleles per locus, percentage of polymorphic locus, the expected and observed heterozygosity, the fixation index and it was estimated the deviation of the Hardy–Weinberg Equilibrium using the GenAlex software [17], for all the samples and grouping the individuals according to their life stage. An Analysis of Molecular Variance (AMOVA) was performed on genetic distances. Wright’s F statistics (Fis, Fst, Fit) was calculated for the whole set of populations and each age category with the GenAlex software [25], and also the degree of genetic differentiation between pairs of populations (Fst) was estimated for all the samples and grouped by life stage.


3. Results

3.1 Demography

The density of the nine locations averaged 2050 individuals ha−1 ± 293.36. Density ranges from 3775 individuals ha−1 in SA to 775 individuals ha−1 in MO. The populations with the highest density of individuals are the localities of SA and GA, and the lowest density is MO, although there are no differences between localities in the number of individuals ha−1 of the populations (ANOVA, df = 8; F = 2.287; P > 0.05).

The stages of the life cycle (Table 2) showed that seedlings were the most abundant, with an average of 49.27% ± 5.00, varying from the lowest in RN to the highest in SA. The proportion of young plants averaged 23.05% ± 2.62, with CH standing out with the lowest proportion and RN with the highest. The proportion of non-reproductive adults was 11.39% ± 3.12, being absent in SA and RN localities, unlike CH, where they represent a third of the population. Reproductive adults represented 16.29% ± 4.48, being scarce in SA and RN, but in the last, they are almost half of the population.

Non- reproductive

Table 2.

Frequency of plants per stage of their life cycle in the locations.

The survival curve based on the states of the life cycle of the populations is an inverted “J” (type III), well defined in most locations (Figure 2F,G,H,I), and scarcely defined in other localities (Figure 2D,E). The localities of SA (Figure 2A), GA (Figure 2B) and AN (Figure 2C) showed a curve like type II in the populations.

Figure 2.

Survival curves according to the life cycle stages of the nine populations of Dioon edule. (A) SA, (B) GA, (C) AN, (D) CH, (E) PO, (F) AB, (G) EJ, (H) RN, (I) MO.

The population’s potential reproductive capacity was 250 ha−1 ± 45.26 individuals in the localities with reproductive structures, of which 158 ha−1 ± 35.36 were male and 92 ha−1 ± 17.18 female. The highest number of male individuals with the reproductive structure were GA and RN, and with the lowest number was SA (Figure 3A).

Figure 3.

Ratio of (A) the density of male (solid) and female (hollow) Dioon edule plants (± SE) and (B) the sex ratio in each location. Different letters mean significant differences (Tukey <0.05).

Likewise, the localities with the highest number of female individuals were RN and CH, and those with the lowest number were MO and AN (Figure 3A). The average sex ratio for the locations was 2: 1 (male: female), with the highest ratio of 4: 1 in GA and the lowest of 0.5: 1 in CH (Figure 3B). The ANOVA showed that there are no differences between mean localities in the number of female individuals ha−1 (n = 18; df = 8; F = 1.234; P > 0.05) and male individuals ha-1 (n = 18; df = 8; F = 1.006; P > 0.05) of the populations.

An average of 35.70% ±2.79 of individuals was found in the transect up to 25 m, mainly seedlings and juveniles in the localities. From 25 to 50 m was found on average 34.39% ± 5.60, in the 50 to 75 m 19.89% ± 3.41 and in the 75 to 100 m, 10.02% ± 2.44. The 70% of the individuals were concentrated in the first 50 m from the transect’s starting point and particularly in RN 100%. At sites with a steeper slope, there were plants widely dispersed, while at sites with flatter slope (less than 20%, e.g. SA), there were large groups of individuals aggregated around adult plants or in the shade of oaks (Quercus laeta) and palms (Brahea dulcis).

Considering the plants’ size and the number of leaves as parameters associated with the cycads’ age, the longest and largest plants are in the communities of PO, CH, GA, AN, and SA. If it was impossible to measure the adult plants’ height and diameter because they had a subterranean stem, only the number of leaves was considered, obtaining similar results.

The average height for non-reproductive individuals was 23.4 cm ± 3.4, with a maximum average height of 52.2 cm and the lowest of 14.5 cm. For reproductive plants, the maximum recorded height was 56.3 cm, and the lowest was 9.5 cm, with an average of 28.8 cm ± 8.3. However, the ANOVA showed that there are differences among means in plant height between localities (df = 8; F = 5.428; P < 0.001) and between categories of reproductive status in localities (df = 8; F = 4.992; P < 0.05), but not between categories of reproductive status in each locality (df = 1; F = 2.633; P > 0.05) (Figure 4).

Figure 4.

Size of Dioon edule plants (± SE) in each location. (A) Non reproductive and, (B) reproductive. Different letters mean significant differences (Tukey <0.05).

The average stem diameter was 12.2 cm ± 2.5 in non-breeding individuals, with a width of 23.6 cm to 15.0 cm (Figure 4A), and in reproductive individuals, the average was 16.5 cm ± 4.4 with a maximum of 31.2 cm and a minimum of 5.5 cm. The ANOVA showed that there are no differences in the mean diameter between localities (df = 8; F = 1627; P > 0.05) or between categories of reproductive status in each locality (df = 1; F = 2848; P > 0.05), but there are differences between the categories of reproductive status of individuals in the localities (df = 8; F = 4034; P < 0.05) (Figure 4).

The average number of leaves was 17.69 ± 2.18 in all individuals. The average number of leaves was 17.69 ± 2.18 in all individuals. In nonreproductive adults it was 14.3 with a maximum of 21 and a minimum of 3.5 (Figure 4A), for reproductive adults it was 15 with a maximum of 22.5 leaves. The ANOVA showed that there are differences in the mean number of leaves between localities (df = 8; F = 3.166; P < 0.05) and between categories of reproductive status in each locality (df = 1; F = 7.312; P < 0.05) but not between the categories of reproductive status of individuals in the localities (df = 8; F = 1.381; P > 0.05) (Figure 4).

3.2 Genetic diversity

Only three of the first seven proposed were amplified, ED9, TOM5, and CAP5. All the loci were polymorphic, but the population AN have the lowest proportion of polymorphic loci (33.33%) in the seedlings (Table 3). The average number of alleles per locus was 3.66; in the seedling stage, the lowest value corresponds to Los AN (1.33) (Table 3). There were no unique alleles within the adult stage. However, in the juvenile and seedling, four and three unique alleles were found, respectively, and ED9 was the locus that presented the highest number of unique alleles.

LocationLife cycle stageNNaPolymorphic Loci (%)HoHeF

Table 3.

Genetic diversity index for five Dioon edule populations in San Luis Potosi, Mexico.

The average heterozygosity (Ho) was 0.207. The highest value was in the seedlings in PO with 0.333, and the lowest in the seedlings in PO (Ho = 0.0), and the adult stage in CH (0.048) (Table 3). For the expected Heterozygosity, the average was 0.504 and the highest value was 0.634 in the juvenile category for the PO locality.

In the populations, only two loci resulted in Hardy–Weinberg equilibrium (H-W), ED9 in the populations of RN, PO, and AN, and TOM5 in the population of RN (Table 4). In the seedlings, the ED9 locus resulted in H-W balance in the populations of CH, PO, and SA and TOM5 resulted in balance in RN. The loci CAP5 in PO and SA populations and ED9 in SA, were monomorphic, and CAP5 in PO and SA populations and ED9 in SA. In the juveniles, ED9 resulted in H-W equilibrium in all populations. However, in adults, this same locus resulted in equilibrium in most populations except CH (Table 4).


Table 4.

Hardy–Weinberg equilibrium test for three microsatellite loci for five Dioon edule populations in San Luis Potosi, Mexico.

The AMOVA showed that the categories’ differences are explained by the genetic variation between individuals for all populations and each age category (Seedlings = 50%; juveniles: 59%; adults: 62%: Total: 58%) (Table 4). The inbreeding coefficient (FIS) generally shows a deficiency of heterozygotes in all populations and age categories (Fis ≈ 0.666–0.708) (Table 5).

Life cycle stageSourcedfSSMSVAR%

Table 5.

Analysis of molecular variance (AMOVA) based on three microsatellite loci for five Dioon edule populations in San Luis Potosi, Mexico. BP = between populations; BI = between individuals; WI = within individuals.

The fixation index (Fst) indicates a moderate differentiation in all populations, as well as in the age categories (Fst ≈ 0.11–0.115). The populations that share more alleles are CH and PO located in the south of the study area (Fst = 0.057; p = 0.006); those that share fewer alleles are PO and SA also in the south (Fst = 0.265; p = 0.001) (Table 5). The seedlings have a high genetic differentiation between pairs of populations (Fst > 0.26; p < 0.05), while in the juveniles, the differentiation is moderate (Fst > 0.110; p < 0.05), and in adults from moderate to high (Fst > 1.15; p < 0.05).


4. Discussion

In the study area the average density was 2103 plants ha−1, lower than the 3200–4633 plants ha−1 of Dioon edule in the center of Veracruz [4, 18], but higher than the 37 plants ha−1 of D. purpusii [26], 116 to 167 plants ha−1 of D. merolae [27] and 589 plants ha−1 of D. holmgrenii [28].

The absence of statistically significant differences between locations is mainly due to the parameters’ wide variation observed. However, the populations with lower density coincide with other studies [27, 28], where the habitat’s quality has decreased. The low quality of the habitat is due to the recurrent disturbance caused by agriculture, mainly the establishment of sugar cane fields in RN, pastures for livestock in EJ and MO, and the extraction of complete plants or their parts, like seeds and leaf harvest for tamale production [1, 29, 30]. In general, plant density and size are population attributes that indicate habitat quality [7].

Most populations presented the survival inverted “J” curve (type III) typical of shade-tolerant species [20] and most populations of D. edule [18], D. purpusii [26], and D. angustifolium [31]. These populations have a higher density of seedlings and juveniles with high mortality and few very long-lived adults with reduced mortality. Fewer populations have a relatively constant mortality rate throughout their life cycle, different from most Dioon populations. This modification could be attributed to the reduction of populations caused by habitat disturbing [32] since the natural forest was replaced with cattle pastures.

The plant’s reproductive capacity is of great interest in conservation since it significantly impacts management and recovery strategies [33]. It is remarkable the high number of female plants present compared with other studies. The imbalance of females is probably related to climatic factors, soil quality, or the occurrence of fires, which in the study area is very common due to agricultural activities [18, 26, 34, 35, 36, 37, 38]. The plants’ sex ratio average was 2: 1 (male: female) in all locations. However, the ideal ratio for Dioon species is 3:1 to ensure an abundant amount of pollen during female strobili receptivity [13]. There are five populations of concern that could present pollen adequacy problems during fertilization, mainly the CH locality that presented a ratio of 1:2.

The Dioon edule populations’ distribution pattern is aggregated, which has long been attributed to the strobilus remaining attached to the plant’s stem during dispersal, causing many seeds to remain close and around. Seed dispersal agents are currently being studied to identify them in the study area, including reptiles, birds, and mammals. It is relevant the evidence of asexual multiplication, which means that sprouts could be an effective strategy when sexual reproduction is difficult or in response to some stress. Like the cycad Cycas armstrongii, which produces sprouts after fires around the charred stems [36], facilitating rapid propagation and ensuring populations’ maintenance [39].

The differences in height between reproductive and non-reproductive individuals have been attributed to the lowest growth rate of reproductive individuals due to the compression of stem tissues by the increase in weight of stored water and the production of heavy strobili [18, 40]. However, in this study, no significant differences were observed in any location, indicating that this difference does not always apply to all cycad populations.

Considering that there were significant plants’ height differences among populations, it would be an indicator of the site quality. The sites of CH, GA and AB that presented the tallest plants would have the best quality since they present the lowest proportion of organic matter and smooth slopes. Diference in height and diameter of Ceratozamia matudae [41], Dioon purpusii [26] and D. holmgrenii [28] is attributed to disturbance conditions.

The plants’ diameter did not present significant differences between localities. It has already been observed in Cycas micronesica [40] between reproductive and non-reproductive individuals, which can be attributed to the manoxylic wood type with a lot of succulent tissue [42]. In the same way, the difference in the number of leaves present in reproductive and non-reproductive individuals at each location and the difference between locations can be attributed that during the adult stage, the incidence of solar radiation correlates with leaf production, but the production of strobili in female individuals stops the production of new leaves, due to the greater allocation of resources for the development of the strabo [43]. The number of leaves in reproductive adults was lower than in Dioon holmgrenii, but much higher in non-reproductive ones [28].

As highlighted in this study, another determining factor in the development of Dioon edule populations is the deterioration of forests caused by agricultural activities [44]. These activities cause a lower rate of reproduction, capacity to adapt to changes, and less ability to establish themselves during the passage from one phenological stage to another, altering the age structure of the population [7].

The genetic diversity analysis tried to clear the picture of the medium and long-term repercussions of adverse factors identified when analyzing the species’ population demography. However, we consider that the number of populations and molecular markers should be increased since Dioon edule is ditributed in various habitats with different microenvironmental conditions, causing phenotypical differences that should be analyzed [3, 31].

The populations analyzed show low genetic diversity and an excess of homozygotes in the seedlings, juveniles, and adults. Contrary to that expected for adult plants, which were believed to have become established before the disturbances. Similar behavior was present in a study conducted in Veracruz, mentioning that populations outside the Pleistocene flora and fauna refuges located in Veracruz and Oaxaca probably suffered a bottleneck due to exposure to adverse environmental conditions generating a significant decrease in population size reflected in the low levels of observed heterozygosity detected [10].

Previous studies suggest that life history issues, such as geographic distribution area, reproduction systems, and pollen and seed dispersal systems, have essential consequences on genetic structure and genetic variability levels in natural populations [4, 6, 10, 45, 46]. These factors are especially relevant in cycads because of their distribution in particular areas. They also present vegetative reproduction forming inside the populations [47].

In addition to the patches of clonal sprouts, the seedling banks that form around adult plants of Dioon edule and trees of Quercus laeta and Bursera simaruba, act as an extensive reservoir of genetic variation. However, this eventually is lost due to the high mortality rates that this species presents in the establishment’s early stages. The mortality occurs naturally in populations because seedlings and juveniles slowly develop features that allow them to tolerate adverse conditions, such as prolonged drought or insolation [48]. These conditions can be aggravated by habitat fragmentation [6, 7, 44], causing gene drift in the population and making natural selection less effective by expressing deleterious or poorly tolerant genotypes to current environmental conditions [49].

Deviations in H-W were at most of the loci assessed, all tending towards a deficit of heterozygotes. The deviations may be due to the lag in the reproductive season, which can vary from four to 52 years in female plants [4]. This produces a time barrier that prevents individuals’ panmictic crossing, causing groups to form that can function as independent populations.

Five unique alleles were identified with low frequencies in the ED9 and TOM5 locus in seedlings and juveniles. These alleles’ appearance could be attributed to recurrent mutations or gene flow by pollen dispersal at great distances [50]. These alleles may benefit from genetic drift promoted by population size reduction.

AMOVA indicates a high differentiation in seedlings (FST = 0.26), probably due to the reduction in gene flow, which has led to inbreeding within the populations. Simultaneously, in adults and juveniles, the values are very similar and may be maintained by self-incompatibility that creates a dependence on pollinators to facilitate reproduction [46].

Some populations share some alleles, suggesting that they were connected by natural corridors some time ago, mainly in the vicinity of the Sierra Gorda in the northern of Querétaro and south of San Luis Potosí. However, in the current generations of seedlings, the Fst values indicate that the populations analyzed are becoming isolated. The isolation is probably caused by physical barriers associated with anthropic activities increasing in the area. This situation diminishes the habitat conditions and leads to higher mortality, affecting pollinators and favoring inbreeding and gene drift [6].


5. Conclusions

Some Dioon edule populations have declined substantially due to land-use change, primarily for sugar cane fields and cattle pastures, but the overall decline is relatively low. Populations genetic diversity suggests that at some time, they were connected by natural corridors. However, when combining the own high mortality of the species in the seedling and juvenile stages, the reduction and fragmentation of its habitat, and the decrease of pollinators, it favors endogamy and genetic drift, whose effects are more pronounced in small populations.



GRM thanks the support of CONACYT (Scholarship for Postgraduate Studies 290685) and the PMPCA-UASLP. The Cycad Society financed the research through the project “The chamal (Dioon edule Lindl.) in the state of San Luis Potosí” and CONACYT (104739), granted to LYE. IICA financed the project “The Chamal, a toxic but also edible cycad” granted to AP. We thank Raymundo Mora for his support in the field and Marc Hersh for his laboratory support.


  1. 1. Yáñez-Espinosa L. Chamal (Dioon edule Lindl.) in the State of San Luis Potosí, México. The Cycad Newsletter. 2009; 32:19-21
  2. 2. de Luca P, Sabato S, Vazquez-Torres M. Distribution and variation of Dion edule (Zamiaceae). Brittonia. 1982; 34 355-362.
  3. 3. Whitelock L. Variation in the mexican cycad Dioon edule (Zamiaceae). The Botanical Review. 2004;70 (2):240-49. doi:10.1663/0006-8101(2004)070[0240:VITMCD]2.0.CO;2
  4. 4. Octavio-Aguilar P. Demographics and population genetics for three locations in Dioon edule Lindl. (Zamiaceae) in the state of Veracruz [thesis]. Xalapa, Veracruz, Mexico: Instituto de Ecologia, A.C.;2009. [Spanish]
  5. 5. Maldonado-Ruiz MF, Flores-Vázquez JC. Entomofauna asociada a Dioon sp. nov. y actividad de los insectos polinizadores en San Jerónimo Taviche, Oaxaca, México. In: Proceedings of Cycad 2008 The 8th International Conference on Cycad Biology; 13-15 January 2008; Panama City, Panama: Memoirs of The New York Botanical Garden, 106 ; 2012. p. 295-300 [Spanish]
  6. 6. Cabrera-Toledo D, González-Astorga J, Vovides AP, Casas A, Vargas-Ponce O, Carrillo-Reyes P, Nolasco-Soto J, Vega E. Surviving background extinction: Inferences from historic and current dynamics in the contrasting population structures of two endemic Mexican cycads. Population Ecology. 2019:61:62-73.
  7. 7. Lopez-Gallego C. Demographic variation in cycad populations inhabiting contrasting forest fragment. Biodiversity Conservation. 2008;17:1213-1225.
  8. 8. Godoy JA. Genetics, molecular markers, and species conservation.Ecosistemas.2009;18:23-33. [Spanish]
  9. 9. Griffiths A, Miller J, Suzuki D, Lewontin D, Gelbart W. Genetics. 7th ed. Madrid: McGraw-Hill-Interamericana; 2002. 849p. [Spanish]
  10. 10. González-Astorga J, Vovides AP, Ferrer M, Iglesias C. Population Genetics of Dioon edule Lindl. (Zamiaceae, Cycadales): Biogeographical and Evolutionary Implications. Botanical Journal of the Linnean Society. 2003;80(3):457-467. https://doi.org10.1046/j.1095-8312.2003.00257.x
  11. 11. Xie J, Jian S, Liu N. Genetic variation in the endemic plant Cycas debaoensis on the basis of ISSR analysis. Australian Journal of Botany. 2005;53:141-145.
  12. 12. Yang, Y, Li, Y, Li, L-F, Ge, X-J, Gong, X. Isolation and characterization of microsatellite markers for Cycas debaoensis YC Zhong et CJ Chen (Cycadaceae). Molecular Ecology Resources. 2008;8:913-915.
  13. 13. Cabrera-Toledo D, González-Astorga J, Vovides AP. Heterozygote excess in ancient populations of the critically endangered Dioon caputoi (Zamiaceae, Cycadales) from central Mexico. Botanical Journal of the Linnean Society. 2008;158:436-447.
  14. 14. Zhan QQ , Wang JF, Gong X, Peng H. Patterns of chloroplast DNA variation in Cycas debaoensis (Cycadaceae): conservation implications. Conservation Genetics. 2011;12:959-970.
  15. 15. Clugston JAR, Kenicer GJ, Milne R, Overcast I, Wilson TC, Nagalingum NS. RADseq as a valuable tool for plants with large genomes—A case study in cycads. Molecular Ecology Resources. 2019;19;1610-1622. DOI: 10.1111/1755-0998.13085
  16. 16. Instituto Nacional de Estadística, Geografía e Informatica [INEGI]. Geographical synthesis of the state of San Luis Potosi. México, D.F.:1985. [Spanish]
  17. 17. Franco-López J, de la Cruz AG, Navarrete SN, Flores MG, Kato ME, Sánchez CS, Abarca ALG, Bedia, CM. Handbook of Ecology. 2nd ed. México: Trillas;2005. 266 p. [Spanish]
  18. 18. Vovides AP. Spatial distribution, survival, and fecundity of Dioon edule (Zamiaceae) in a tropical deciduous forest in Veracruz, Mexico, with notes on its habitat. American Journal of Botany. 1990;77:1532-1543.
  19. 19. Matteucci S, Colma A. Methodology for the study of vegetation. Serie de biología, monografía no. 22. Washington, D. C.: Secretaría General Organización de los Estados Americanos. Programa Regional de Desarrollo Científico y Tecnológico. 1982. [Spanish]
  20. 20. Carabias J, Meave JA, Valverde T, Cano-Santana Z. Ecology and environment in the 21st century. 1st ed. Mexico: Pearson Educación de México, S. A., 2009. 267 p. [Spanish]
  21. 21. Tremblay R L. The effect of population structure, plant size, herbivory and reproductive potential on effective population size in the temperate epiphytic orchid, Sarcochilus australis. Cunninghamia. 2006;9:529-535
  22. 22. Moynihan J, Meerow A, Francisco-Ortega J. Isolation, characterization and cross-species amplification of microsatellite loci in the cycad genus Dioon (Zamiaceae). Potential utilization in population genetics studies of Dioon edule. Molecular Ecology Notes. 2007;7:72-74. doi:10.1111/j.1471-8286.2006.01531.x
  23. 23. Shuelke M. An economic method for the fluorescent labeling of PCR fragments. Nature Biotechnology. 2000;18:233-34.
  24. 24. Blacket MJ, Robin C, Good RT, Lee SF, Miller D. Universal primers for fluorescent labelling of PCR fragments—an efficient and cost-effective approach to genotyping by fluorescence. Molecular Ecology Resources. 2012:12: 456-463
  25. 25. Peakall R, Smouse P. Genalex 6: Genetic Analysis in Excel. Population Genetic Software for Teaching and Research. Molecular Ecology Notes. 2006;6 (1):288-95. doi:10.1111/j.1471-8286.2005.01155.x
  26. 26. Yáñez-Espinosa L, Sosa-Sosa F. Population structure of Dioon purpusii Rose in Oaxaca, México. Neotropical Biology and Conservation. 2007;2:46-54
  27. 27. Lázaro-Zermeño JM, González-Espinosa M, Mendoza A, Martínez-Ramos M, Quintana-Ascencio PF. Individual growth, reproduction and population dynamics of Dioon merolae (Zamiaceae) under different leaf harvest histories in central Chiapas, México. Forest Ecology and Management. 2011:261:427-439.
  28. 28. Velasco-García MV, Valdez-Hernández JI, Ramírez-Herrera C, Hernández-Hernández ML, López-Upton J, López-Mata L, López-Sánchez, H. Structure, stage heterogeneity and spatial dispersion pattern of Dioon holmgrenii (Zamiaceae). Botanical Sciences. 2016;94:75-87. [Spanish]
  29. 29. Rubio-Méndez G, Chávez-Acuña IJ, Yáñez-Espinosa L, Fortanelli-Martínez J. Effect of seed collection on the structure and dynamics of wild populations of Dioon edule (Zamiaceae). The Journal of the Torrey Botanical Society. (in press)
  30. 30. Tristán-Martínez E, Fortanelli-Martínez J, Bonta M. Toxic Harvest: Chamal Cycad (Dioon edule) Food Culture in Xi'Iuy Indigenous Communities of San Luis Potosi, Mexico. Journal of Ethnobiology. 2020:4:519-534.
  31. 31. Rubio-Méndez G, Yáñez-Espinosa L, Salinas-Rodríguez MM, Hernández-Hernández KA, Balderas-González D. Demographics and population structure of Dioon angustifolium (Zamiaceae) in northeastern Mexico. Botanical Sciences. 2019;97(4):685-690. [Spanish]
  32. 32. Keppel G. Notes on the natural history of Cycas seemannii (Cycadaceae). South Pacific Journal of Natural Science. 2001;19:35-41.
  33. 33. Gaol ML, Fox ED. Reproductive potential of Acacia species in the central wheatbelt: variation between years. Conservation Science Western Australia Journal. 2002;4:147-152
  34. 34. Alejandre-Rosas J, Sanchez-Tinoco MY, Vazquez-Torres M. Population structure of Ceratozamia mexicana Brong. (Zamiaceae) in a forest in central Veracruz. La Ciencia y el Hombre. 1990;1:93-112. [Spanish]
  35. 35. Ornduff R. Size classes, reproductive behavior, and insect associates of Cycas media (Cycadaceae) in Australia. Botanical Gazette. 1991;152:203-107.
  36. 36. Watkinson AR, Powell JC. The life history and population structure of Cycas armstrongii in monsoonal northern Australia. Oecologia. 1997;111:341-349.
  37. 37. Nicolalde-Morejón, F. Ecology and Taxonomy of Zamia in Ecuador. Quito, Ecuador: Fundación Ecuatoriana para la Investigación y Desarrollo de la Botánica, Boletín 9; 2001. [Spanish]
  38. 38. Pérez-Farrera MA, Vovides AP. Spatial distribution, population structure, and fecundity of Ceratozamia matudae (Zamiaceae) in El Triunfo Biosphere Reserve, Chiapas, México. The Botanical Review. 2004;70:299-311.[0299:sdpsaf];2
  39. 39. Brenes-Cambronero L, Di Stefano JF. Possible influence of ramets on the population structure and distribution of the giant tree Warszewicsia uxpanapensis, cordillera Tilarán, Costa Rica. Revista Biología Tropical. 2006;54:1179-1188. [Spanish]
  40. 40. Niklas K, Marler TE. Sex population differences in the allometry of an endangered cycad species Cycas micronesica (Cycadales). International Journal of Plant Sciences. 2008;169: 659-665.
  41. 41. Pérez-Farrera MA, Quintana-Ascencio PF, Salvatierra-Izaba B, Vovides AP. Population dynamic of Ceratozamia matudai in el Triunfo Biosphere Reserve, Chiapas, México. Journal of the Torrey Botanical Society. 2000;127:291-299.
  42. 42. Stevenson DW. Radial growth in the Cycadales. American Journal of Botany. 1980;67:465-475.
  43. 43. Clark D, Clark DB. Leaf production and the cost of reproduction in the Neotropical rain forest cycad Zamia skinneri. Journal of Ecology. 1987;76:1153-1163.
  44. 44. Donaldson, JS, editor. 2003. Cycads status, survey and conservation action plan. Cambridge, UK: IUCN/SSC Cycad Specialist Group. IUCN. 2003. ix+86 pp
  45. 45. Bohonak A. Dispersal, gene flow and population structure. The Quarterly Review of Biology. 1999;74(1):21-45.
  46. 46. Mathiasen P, Rovere A, Premoli A. Genetic Structure and Early Effects of Inbreeding in Fragmented Temperate Forest of a Self-Incompatible Tree, Embothrium coccineum. Conservation Biology. 2006;21 (1):232-40. doi:10.1111/j.1523-1739.2006.00565.x
  47. 47. Cabrera-Toledo D, González-Astorga J, Flores-Vázquez J. Fine-scale spatial genetic structure in two mexican cycad species Dioon caputoi and Dioon merolae (Zamiaceae, Cycadales): Implications for conservation. Biochemical Systematics and Ecology. 2012;40:6.
  48. 48. Augspurger CK. Seedling survival of tropical tree species: Interactions of dispersal distance, light gaps, and pathogens. Ecology. 1984;65:1705-1712.
  49. 49. Allendorf FW, Luikart GH. Conservation and the Genetics of Populations. 2nd ed. UK:Blackwell Publishing. 2007 ISBN: 978-0-470-67146-7
  50. 50. Acuña W, Yalta C, Veli E. Cross species transferability of microsatellite markers from Anas platyrhynchos to Peruvian Muscovy Duck Cairina moschata domestica. Revista Peruana de Biología. 2020:27:255-260.

Written By

Gabriel Rubio-Méndez, Alberto Prado, Jacqueline C. Bede, José Arturo De-Nova, Joel Flores, Juan Antonio Reyes-Agüero and Laura Yáñez-Espinosa

Submitted: September 17th, 2020 Reviewed: February 3rd, 2021 Published: March 9th, 2021