Landscape genetics reveals inbreeding and genetic bottlenecks in the extremely rare short-globose cacti Mammillaria pectinifera (Cactaceae) as a result of habitat fragmentation
Maya-García Reynaa, Santiago Arizagab, Cuevas-Reyes Pabloa, Peñaloza-Ramírez Juan Manuelb, Víctor Rocha Ramírezb, Ken Oyamab     
a. Facultad de Biología, Universidad Michoacana de San Nicolas de Hidalgo, Mexico;
b. Escuela Nacional de Estudios Superiores Unidad Morelia, UNAM, Antigua Carretera a Pátzcuaro 8701, Colonia Ex-Hacienda de San José de la Huerta, 58190, Morelia, Michoacán, Mexico
Abstract: Mammillaria pectinifera is an endemic, short-globose cactus species, included in the IUCN list as a threatened species with only 18 remaining populations in the Tehuacán-Cuicatlán Valley in central Mexico. We evaluated the population genetic diversity and structure, connectivity, recent bottlenecks and population size, using nuclear microsatellites. M. pectinifera showed high genetic diversity but some evidence of heterozygote deficiency (FIS), recent bottlenecks in some populations and reductions in population size. Also, we found low population genetic differentiation and high values of connectivity for M. pectinifera, as the result of historical events of gene flow through pollen and seed dispersal. M. pectinifera occurs in sites with some degree of disturbance leading to the isolation of its populations and decreasing the levels of gene flow among them. Excessive deforestation also changes the original vegetation damaging the natural habitats. This species will become extinct if it is not properly preserved. Furthermore, this species has some ecological features that make them more vulnerable to disturbance such as a very low growth rates and long life cycles. We suggest in situ conservation to prevent the decrease of population sizes and loss of genetic diversity in the natural protected areas such as the Tehuacán-Cuicátlán Biosphere Reserve. In addition, a long-term ex situconservation program is need to construct seed banks, and optimize seed germination and plant establishment protocols that restore disturbed habitats. Furthermore, creating a supply of living plants for trade is critical to avoid further extraction of plants from nature.
Keywords: Mammillaria pectinifera    Nuclear SSR    Genetic diversity    Recent bottlenecks    Effective population size    
1. Introduction

Loss of continuous natural areas and conversion of forests to smaller areas have drastically reduced the size of plant and animal populations, increasing the isolation between populations of the same species (Young and Clarke, 2000; Lowe et al., 2005). Plant reproduction may decrease in fragmented populations because plants may receive fewer flower visitors due to a decline in the richness and abundance of pollinators, modification in species composition and limitation in movement among landscapes (Young and Clarke, 2000; Cascante et al., 2002; Goverde et al., 2002; Frankham, 2010). Disturbance also changes the abundance, richness or behavior of animal seed dispersers affecting seed dispersal (Ghazoul, 2005; Lowe et al., 2005), recruitment and establishing of seedlings (Ewers and Didham, 2006), and the spatial genetic structure of plant populations (Hamrick, 2004). Shortly following a disturbance, levels of genetic homozygosity in affected plant populations are predicted to increase, and this may lead to the expression of recessive alleles deleterious to individual fitness (Hart and Clark, 1997; Reed and Frankham, 2003). In the long term, the loss of genetic variation may induce severe population bottlenecks due to genetic drift, high population divergence and low effective population size (Zhang et al., 2004; Moreira et al., 2009).

Dry ecosystems around the world are experiencing major environmental impacts due to climate change since water is the major limiting factor for plant productivity in such environments (Fischer and Turner, 1978; Anderson et al., 1994; Hamrick et al., 2002). Other factors, such as strong drying winds, extreme temperatures, limited nutrients, high light intensity, overexploitation of natural resources, and desertification have caused habitat degradation, fragmentation and loss of biodiversity (McNeely et al., 1990; Ayyad, 2003).

The Cactaceae family includes approximately 1600 species distributed in arid and semi-arid habitats (Gibson and Nobel, 1986; Dávila et al., 2002). Mexico is the most important center of diversification with almost 850 species and 54 genera (Bravo-Hollis, 1978; Arias, 1993; Anderson et al., 1994; Hernández and Godínez, 1994). The Balsas River Basin and the Tehuacán Valley in central Mexico have a great number of cacti with approximately 81 species, of which 25% are endemic, constituting one of the principal structural and floristic elements of the dry forests (Arias et al., 1997). The species of Mammillaria predominate in this area and have biological characteristics that make them more vulnerable to the effects of disturbance in their populations such as low growth rates and long life cycles (Gibson and Nobel, 1986). Particularly, Mammillaria pectinifera is an extremely narrowly distributed species with only 18 fragmented populations remaining in the Tehuacán-Cuicatlán Valley in the states of Puebla and Oaxaca in Central Mexico (Martorell and Peters, 2005; Valverde et al., 2009). This species grows on relatively deep alkaline calcareous soils and the plants are localized on small isolated patches in xerophyllous scrub and grasslands (Arias et al., 1997). The populations of M. pectinifera fae serious risk of extinction, mainly due to massive collection of plants for commercial purposes and the degradation of their habitats (Glass, 1998; Valverde and Zavala-Hurtado, 2006). Grazing also causes the mortality of plants of M. pectinifera in its natural habitats (Meyrán, 1973; Bravo-Hollis, 1978; Bravo-Hollis and Sánchez-Mejorada, 1991). Thus, M. pectinifera is included in the category of risk of extinction in international (IUCN; CITES, 2013) and national (NOM 059 ECO1 2001; SEMARNAT, 2010) regulations.

Most of the previous population genetic studies within the Cactaceae family have been mainly conducted on columnar and short globose cacti using allozymes (Nassar et al., 2001; Hamrick et al., 2002; Moraes et al., 2005; Parra et al., 2008), RAPDs (Clark-Tapia et al., 2005) and microsatellites (Figueredo et al., 2010; Castro-Félix et al., 2013; Solórzano-Lujano et al., 2014). Cornejo-Romero et al. (2014) found that populations of species of Mammillaria were genetically structured indicating that geographic isolation and limited dispersal were the primary causes of genetic population differentiation.

Ecological studies on Mammillaria have assessed the risk of extinction of its species due to the high specificity of the habitats where they occur and the degree of disturbance of the natural environments (Contreras and Valverde, 2002; Martorell and Peters, 2005; Valverde and Zavala-Hurtado, 2006; Valverde et al., 2009; Peters et al., 2014). Many aspects about genetic structure of the populations of M. pectinifera are unknown. Therefore, it is important to conduct studies to assess the genetic diversity of their populations and provide some guidelines for conservation. The main objective of this study was to examine the levels of genetic diversity of populations of M. pectinifera and to characterize its population genetic structure, inbreeding, recent bottlenecks and fluctuations in the effective population size to identify the factors that are determining the evolution and maintenance of this endangered microendemic species, and to propose conservation strategies to preserve this species.

2. Materials and methods 2.1. Study species

M. pectinifera F.A.C.Weber in Bois (Cactaceae) is a small cactus that in habitat appears above the ground as a squat stem completely covered by flattened spines, and usually solitary. Flowers are medium-sized in a ring around the sides of the globe, white to pink in color with darker mid-strip, 20e30 mm long and stigmas are green (Arias et al., 1997). This species has a blooming season from December to March. Reproduction starts at an age of eight years. Fruits are small, red at maturity, barely emerging above the spines. Seeds, usually black, are retained among the plant tubercles and are released gradually (Arias et al., 1997).

We collected from 18 populations of M. pectinifera in the Tehuacan-Cuicatlan Biosphere Reserve in central Mexico; from 7 to 13 individuals were sampled randomly with at least 100 m of distance between each individual. In all the samples, we did cuts about 1 cm2 to get photosynthetic tissue, subsequently, the tissue was frozen at-70 ℃ until used for DNA extraction. For the DNA extraction (184 individuals), we followed the protocols of Doyle and Dickson (1987) and Doyle and Doyle (1987). Six nuclear DNA (nSSR) microsatellite loci were selected and amplified in multiplex polymerase chain reactions (PCR). Two groups of primers were arranged according to allele size and fluorescent labels. The first group was formed by the primer pairs for Mam VTC2, Mam VTC8, Mam VTC9, Mam VTC11, previously designed for the species M. crucigera by Solórzano-Lujano et al. (2009), whereas the second group included the primer pairs for P44 and P50, previously designed for Polaskia chichipe (Otero-Arnaiz, 2005). PCR was performed using the QIAGEN Multiplex PCR kit (QIAGEN) in a volume of 5 ml containing 1× Multiplex PCR Master Mix, 2 μM each primer, dH2O, and 20 ng template DNA. The thermal cycling conditions consisted of 40 cycles, each at 95 ℃ for 1 min, annealing for 1 min (i.e. first and second primer groups 58 ℃ and 60 ℃, respectively), extension for 2 min at 72 ℃ and a final extension at 72 ℃ for 10 min. Multiplex PCR products were combined with a GeneScan-500 LIZ size standard and sequenced by an ABI-PRISM 3100 Avant sequencer (Applied Biosystems). Fragments were analyzed and registered with the Peak Scanner program 1.0 (Applied Biosystems).

2.2. Genetic diversity

To test whether there was evidence of null alleles, upper allele dropout, or small genotyping errors due to stutter in the SSR data that we accumulated in 18 populations of M. pectinifera, we used the Micro-Checker v 2.2.3 program with 102 bootstrap simulations and a 95% confidence interval (Van Oosterhout et al., 2005). We calculated the values and standard error of mean number of alleles per locus (Na), mean effective number of alleles (Ne), mean observed heterozygosity (HO), mean expected heterozygosity (HE) and mean of fixation index (FIS) in 18 populations of M. pectinifera using the GENETIX 4 program (Belkhir et al., 2004).

2.3. Population structure and cluster analysis

Population genetic differentiation RST was estimated by the stepwise mutation model (SMM) performed with 104 permutations in the ARLEQUIN 3.5.1.2 software (Excoffier and Lischer, 2010). A hierarchical test of population structure was estimated using the stepwise mutation model (SMM) performed with AMOVA in ARLEQUIN 3.5. (Excoffier et al., 2005). We compared the variance distribution between groups, among populations between groups and within populations. The statistical significance was tested using 104 permutations utilizing the resulting two groups of genotypes obtained previously by STRUCTURE.

Acluster analysis using a Bayesian approach was conducted using STRUCTURE version 2.3.3 (Pritchard et al., 2000; Falush et al., 2003; Hubisz et al., 2009). In this analysis individuals are probabilistically assigned to one of the predefined K populations (gene pools) to identify the optimal number of genetic groups (Evanno et al., 2005). The optimum number of groups (K) was determined by varying the value of K from 1 to 10 and running the analysis ten times per K value, in order to determine the maximum value of posterior likelihood [LnP (D)]. Each run was performed using 105 burn-in periods and 106 Markov Chain Monte Carlo (MCMC) repetitions after burn-in. We used a model allowing for admixture with correlated allelic frequencies without any prior information. Also, we determined the most probable value of K using the maximum value of △K according to Evanno et al. (2005) implemented in the program Structure Harvester 0.6.1 (Earl and Von Holdt, 2011).

To identify probable geographic and genetic breaks among populations of M. pectinifera, we used the Monmonier's maximum difference algorithm with the software BARRIER version 2.2 (Manni et al., 2004) that creates a map of sampling locations from their geographical coordinates. Barriers are then characterized on the map by identifying the maximum values within the population pairwise genetic distance matrix. We employed an average square distances (ASD) matrix of (Goldstein et al., 1995; Slatkin, 1995) estimated for 18 populations of M. pectinifera. Resampling random subsets of individuals within populations provided 100 bootstrap replicate distances were constructed with the MSA program (Dieringer and Schlotterer, 2003) to accomplish statistical signi € ficance for the predicted barriers.

Population connectivity among M. pectinifera was scrutinized with the POPGRAPH library (http://dyerlab.bio.vcu.edu/software.html), this approaches generate a network of population linkages and defines the degree of genetic variation within populations (Dyer and Nason, 2004). In the POPGRAPH context, the set of nodes represents sampled populations, and the edges represent the multivariate measures of genetic covariance among populations. The difference in node size reflects differences in withinpopulation genetic variability, while the edge length represents the among-population component of genetic variation due to the connecting nodes (Dyer and Nason, 2004). POPGRAPH, recognize population pairs where long distance migration may have occurred by indicating population pairs with significantly greater inter-site distances (Dyer et al., 2010). It also identifies pairs of populations that are located significantly closer than predicted by inter-site separation, which suggests that a direct barrier might exist between the linear distances.

2.4. Population bottleneck and effective population size

To detect evidence of recent bottlenecks in M. pectinifera populations, we used the BOTTLENECK 1.2 software (Piry et al., 1999) to examine the two genetic groups previously obtained by STRUCTURE. Due to small sample sizes, our analysis was limited to. Also, England et al. (2006) showed that the disequilibrium method leads to severe downward bias in the estimate of Ne when sample size (S) is small relative to effective size. Therefore, the estimation by population will result in bias due to the small sample size of the population. Recent bottlenecks could be defined as a population where the rare alleles are the first to be lost, diminishing the mean number of alleles per locus. In contrast, heterozygosity is less affected, generating a transient excess in heterozygosity compared to that expected given the resulting number of alleles (Cornuet and Luikart, 1996; Luikart and Cornuet, 1998). To analyze the data set, we used 90% stepwise and 10% multistep mutations 104 iterations with the Wilcoxon signed-rank test, and the stepwise mutation (SMM), the infinite allele (IAM) and two-phase mutation (TPM) models. Moreover, we also assessed the effective population size in populations of M. pectinifera with the program LDNe (Waples and Do, 2008), which implements the bias-correction method developed by Waples (2006) to obtain Ne from each sample of S individuals. For LDNe, we used the criterion Pcrit=0.02 (alleles with frequency ˂ 0.02 are excluded), which generally provides a good balance among accuracy and bias (Waples and Do, 2008). Confidence intervals (CIs) for Ne were based in the chi-square approximation implemented by LDNe (Waples, 2006).

3. Results 3.1. Genetic diversity

In general, we did not find evidence for null alleles over all sample-loci combinations, and tests for error due to stutter and upper allele dropout were negatives in all cases. The values of the genetic diversity parameters estimated in 18 populations of M. pectinifera are reported in Table 1. We observed a wide range of values from 4.50 to 7.00 in mean number of effective alleles per locus (Na) and 0.716 to 0.950 in mean observed heterozygosity (HO) values (Table 1). Wright's inbreeding coefficients within populations (FIS) were positive in almost all cases, indicating a deficit of heterozygotes ranging from FIS=0.071 to 0.621 and some populations with negative values ranging FIS=-0.010 to 0.446 (Table 1).

Table 1 Locality name, sample size, N/W geographical coordinates (degrees), mean number of effective alleles per locus (Ne), mean observed heterozygosity (HO), mean expected heterozygosity (HE), for 18 populations of M. pectinifera in the Tehuacán-Cuicatlán Biosphere Reserve, Mexico. Standard errors are included in parenthesis.
LocalityGenetic diversity
Sample sizeCoordinatesNeHOHEFIS
1.CO1118.413/-97.4206.50 (2.07)0.893 (0.10)0.808 (0.08)0.345 (0.05)
2. CH718.619/-97.5485.83 (2.78)0.734 (0.24)0.761 (0.17)0.327 (0.07)
3.TE1018.403/-97.5487.00 (2.53)0.826(0.11)0.792 (0.14)0.237 (0.01)
4.SM1318.703/-97.6037.00 (3.03)0.845 (0.12)0.785 (0.10)-0.010 (0.05)
5. RE1118.555/-97.6356.50 (2.08)0.932 (0.08)0.808 (0.07)0.017 (0.07)
6. TC1118.865/-97.7025.83(2.31)0.842 (0.16)0.758 (0.10)0.128 (0.12)
7. SB1118.464/-97.5686.50 (2.66)0.950 (0.08)0.795 (0.08)-0.153 (0.06)
8. SI718.409/-97.4294.83 (0.98)0.849 (0.16)0.768 (0.10)-0.162 (0.09)
9.ZA1018.358/-97.4806.50 (2.07)0.866 (0.10)0.793 (0.07)0.326 (0.11)
10. FR1118.476/-97.5725.66(1.21)0.842 (0.10)0.777 (0.03)0.424 (0.01)
11. FT818.302/-97.6955.33 (2.06)0.854 (0.16)0.756 (0.12)0.071 (0.09)
12. SU1018.619/-97.4776.83 (2.40)0.796 (0.12)0.789 (0.07)0.621 (0.08)
13. EC1118.605/-97.5425.83 (1.94)0.772 (0.21)0.751 (0.10)0.225 (0.04)
14. CC1118.636/-97.4224.50 (1.37)0.771 (0.14)0.734 (0.03)0.261 (0.14)
15. LA1118.587/-97.5486.16(2.31)0.886 (0.10)0.772 (0.07)0.234 (0.06)
16. AZ1018.704/-97.4145.00 (1.26)0.948 (0.05)0.753 (0.01)-0.446 (0.19)
17. ER1018.444/-97.4635.16 (1.94)0.716(0.21)0.733 (0.12)-0.116 (0.14)
18. CQ1118.638/-97.4155.00 (1.67)0.903 (0.14)0.766 (0.07)-0.254 (0.05)
3.2. Population structure and cluster analysis

Pairwise genetic differentiation among populations indicated low to moderate genetic differentiation based on RST (Table 2). Hierarchical analysis of molecular variance (AMOVA), performed for the SMM mutation models (i.e. RST), indicated that most of the genetic variation resided within populations (ΦST=80.63%, P=0.001) followed by variation among populations within groups (ΦSC=10.32%, P=0.001), while the differentiation among groups only accounted for the remaining variation (ΦCT=9.04%, P=0.001) (Table 2).

Table 2 Pairwise population genetic differentiations among localities of M. pectinifera in the Tehuacán-Cuicatlán Biosphere Reserve using the matrix of Slatkin linearized RST (SMM), over 104 replicates in the Arlequin 3.5 software (Schneider et al., 2000).
RSTCOCHTESMREFRTCSBSIZAFTSUECCCLAAZERCQ
CO-
CH0.122-
TE0.0080.171-
SM0.0940.1570.097-
RE0.2230.1720.0500.011-
FR0.0050.1000.0020.0550.050-
TC0.0470.3050.0440.0940.2810.047-
SB0.0290.2180.0670.1180.2950.0620.018-
SI0.0580.1720.0530.0060.0860.0440.1110.186-
ZA0.1900.2560.0470.0710.0240.0240.1970.1970.104-
FT0.3370.4230.3220.2410.3950.2310.2960.1470.3810.301-
SU0.1530.3170.0060.0170.1440.0680.0470.0670.0880.0010.325-
EC0.1150.3520.0610.1430.3120.0570.0110.0340.1810.1840.1810.102-
CC0.1610.3660.0070.0250.2190.1060.0410.0870.0750.0750.2120.0680.011-
LA0.0280.2730.1300.2070.3670.0120.0670.0290.2230.3170.3900.1960.0780.246-
AZ0.2080.4110.1230.1170.3020.1580.1060.1370.1970.2270.3600.0620.0460.0300.214-
ER0.2760.4450.1990.1870.3620.1840.1490.0510.2940.2150.0930.1790.0280.0210.2780.146-
CQ0.2250.4350.1450.1910.3720.1180.1090.0780.2710.2490.1710.1000.0250.0780.1810.0100.024-

The highest posterior probability was obtained from Bayesian likelihood [LnP (D)] and the △Kapproach (Evanno et al., 2005) implemented in the program Structure Harvester (Earl and Von Holdt, 2011). This statistic determined that K=2 is the optimum value for the number of genetic clusters included in the analysis (Fig. 1). Both genetic groups showed a very widespread distribution and also some of them exhibited a great amount of exchange among groups. For instance, the first genetic group (i.e. Cluster 1, green) of populations were consistently structured and very widespread over all the distribution of M. pectinifera. In addition, some of the populations had a greater ancestry coefficient Q from 0.669 to 0.847 such as the populations CO, CH, TE, SM, RE, SI, ZA and LA (Fig. 2). Populations from the second genetic group (i.e. Cluster 2, red), such as FR, SU, EC, CC, AC, ER and CQ had a widespread distribution and were very well structured with ancestry coefficient (Q) going from 0.710 to 0.904. The rest of populations TC SB and NO showed a greater amount of exchange among genetic groups with an ancestry ranged from Q 0.410e0.590 (Fig. 2). Table 3.

Fig. 1 Mean and standard deviation values of △K for ten independent runs of STRUCTURE plotted against the number of genetic groups K. The peak indicates the most probable number of genetic groups given the data using Structure Harvester (Earl and Von Holdt, 2011).

Fig. 2 Map shows sampling localities representing 18 populations of M. pectinifera in the Tehuacán-Cuicatlán Biosphere Reserve. Green genotype and red genotype represent the genetic ancestry groups corresponding to M. pectinifera populations obtained by the program STRUCTURE. Black bars show the barriers detected with the BARRRIERS program.

Table 3 Analysis of molecular variance (AMOVA) for populations ofM. pectinifera. Analysis was performed on the nSSR data and RSTwas used for the two genetic clusters obtained by means of STRUCTURE. Asterisks indicate statistically significant values (P < 0.01). Tests were based on 104 random permutations.
Source of variationSum of squaresVariance componentsPercentage of variationFixation index
RST
Among groups10231.347.2119.04ΦCT=0.090***
Among populations within groups24313.453.88110.32ΦSC=0.010***
Within populations147309.7420.88480.63ΦST=0.080***
Total181854.5521.978

In total, four barriers were identified among the 18 populations of M. pectinifera using the Barriers analysis with 100 bootstrap replicates of ASD genetic distance matrixes (Fig. 2). The most significant barrier, with a bootstrap support of 98%, separated TC from SM in the northern part of the distribution. The second barrier, with 90% bootstrap value, indicated that there is a complex barrier that divides SM from the rest of the populations in the center and the south-western region. The third barrier, with 87% bootstrap support, separated populations in the center of the distribution, including, CH, EC, LA and SU. The fourth barrier, with a 79% bootstrap support, separated the populations from RE and NO located in the center-west of the distribution of M. pectinifera (Fig. 2).

The POPGRAPH network of populations based on nSSR genotypes had 17 edges out of the possible 22, which indicates extensive historical connectivity among most populations (Fig. 3). The most common colonization patterns occurred in a northesouth axis along the Tehuácan-Cuicatlán Biosphere Reserve, throughout southern portions of Puebla to Zapotitlán, but it was not uncommon to see eastewest gene exchanges, especially across the northeeast to the southewest part of the species range. Our data indicate occasional long-distance gene exchange across the southern portion of Tehuacán, FR, to the north toward TC and CH, where the genotypes are significantly more similar than expected based on spatial distance. Our analysis also identified a great amount of network connection between the center of the distribution towards the northwest, south, but not much connectivity with the southwest portion, which were more genetically dissimilar than predicted by spatial distance (Fig. 3).

Fig. 3 Patterns of genetic connectivity among 18 populations of M. pectinifera in the Tehuacán-Cuicatlán Biosphere Reserve, Mexico performed with the program POPGRAPH.
3.3. Population bottleneck and effective population size

The results of the analysis to test the evidence of recent bottlenecks (excess of heterozygosity) using the infinite allele model (IAM), two-phase model (TPM) and stepwise mutation models (SMM) are presented in Table 4. Results were significant (P < 0.050) for the IAM P=0.0003 and the SMM P=0.006 models observed in the loci VTCMAM11, P50 and P44. Our analysis revealed that some populations of M. pectinifera departed from neutrality, exhibiting excess heterozygosity. This suggests that the populations have recently undergone a bottleneck. The results of estimation of effective population size (Ne), performed with the program LDNe, in populations of M. pectinifera showed that both genetic groups had moderated values in the first group (Ne=455 individuals), followed by the second group (Ne=89.8); in all cases, estimates had high Jackknife support and a good confidence interval (CIs) (Table 4).

Table 4 Bottleneck analysis for M. pectinifera populations in theTehuacán-Cuicatlán Biosphere Reserve using Wilcoxon rank test under infinite allele, stepwise mutation and two phase model. Parameters for TPM: variance=10%, proportion of SMM=90%, estimation based on 104 replications. P, probability. IAM, infinite allele model; TPM, two phase model; SMM, stepwise mutation model. Also, we include the results obtained for the estimation of the population effective size for the two genetic groups, values obtained with the program LDNe.
ModelsGroup1Group2
IAM0.007***0.039**
TPM0.7810.656
SMM0.046*0.042*
LDNe455.2**89.8**
*, ** and *** Indicate significant deviation from equilibrium as value less than 0.05.
4. Discussion

Anthropogenic activities, such as massive extraction of plants and disturbance of natural habitats, may lead to a decrease in population size and changes in the genetic structure of species populations (Lowe et al., 2005). M. pectinifera has been subject to these pressures, particularly in reproductive stages, which has led to a significant reduction in their populations. Unexpectedly, the populations of M. pectinifera have been shown to have similar (López-Ortiz et al., 2013; Macías-Arrastrio, 2013) or higher (Ibarra-Suárez, 2009; Tapia-Salcido, 2011) values of genetic diversity compared to other short-globose cacti species of the same genus. However, we found significant evidence of recent bottlenecks in some populations of M. pectinifera. Lowe et al. (2005) suggested that after an episode of habitat disruption, populations that remained small for many generations suffer loss of allelic diversity and random genetic drift (Hamrick, 2004; Ellstrand and Elam, 1993). The destruction and removal of reproductive adult trees might result in fecundity variance, reduced effective population size and limited gene flow among the remnant populations (Ewers and Didham, 2006). Furthermore, M. pectinifera has life history traits that make it more vulnerable to the effects of disturbance over the course of generations such as low growth rates, long life cycles and serotiny (Gibson and Nobel, 1986; Bravo-Hollis and Sánchez-Mejorada, 1991; Valverde et al., 2009).

High inbreeding values, moderate to high heterozygosity and low to moderate effective population size are positive indicators that predict short-term population viability in M. pectinifera. While only a relatively small number of individuals may be needed in the short term to prevent inbreeding, greater genetic variation is thought to increase the likelihood of survival over much longer time scales (Frankham, 2010; Jamieson and Allendorf, 2012). Adaptation or resilience to novel or changing habitats is the most important feature (Williams and Hoffman, 2009). Thus, management strategies should be implemented that maintain, as far as possible, the full complement of genetic diversity within the Tehuacán-Cuicatlán Biosphere Reserve.

Hierarchical analysis of molecular variance revealed that most of the genetic variation is within populations, and there is low genetic differentiation among populations. This observation was confirmed by the POPGRAPH analysis that also indicated extensive historical connectivity. The connectivity values obtained for M. pectinifera were relatively high and it is likely the result of historical events of gene flow through pollen and seed dispersion by biotic agents. However, extant populations of M. pectinifera are distributed in isolated patches with few opportunities for seed dispersal across valleys and foothills (Valverde and ZavalaHurtado, 2006; Valverde et al., 2009; and Peters et al., 2014). Furthermore, populations of Mammillaria in the Tehuacán Valley present a restricted spatial distribution, high habitat specificity in restricted altitudinal ranges, soil and vegetation type (Hernandez and Godinez, 1994; Hershkovitz and Zimmer, 1997). In particular, M. pectinifera plants occurred only at altitudinal gradients from 1778 to 2100 m in mild to moderate slopes, on relatively deep alkaline calcareous soils, with high water retention capacity on stony hills (Valverde and Zavala-Hurtado, 2006; Valverde et al., 2009; Peters et al., 2014).

The massive extraction of plants, deforestation, changes in land use, and over-exploitation of areas by human activities have damaged the natural habitats and the populations of M. pectinifera (Zavala-Hurtado and Valverde, 2003; Martorell and Peters, 2005). Despite this, we observed that populations of M. pectinifera still show moderate to high levels of genetic diversity due to their strong capacity to adapt and survive in localities with high water deficits (Hernandez and Godinez, 1994; Hershkovitz and Zimmer, 1997). However, if not properly preserved, this species will become extinct in the future. For instance, according to Peters et al. (2014), 45% of all known populations of M. pectinifera are located outside of natural protected areas. For that reason, we suggest in situ conservation to prevent the decrease in population sizes and loss of genetic diversity. Natural reserves such as the Tehuacán-Cuicatlán Biosphere Reserve must play an essential role in protecting and restoring the habitat of M. pectinifera (Contreras and Valverde, 2002; Peters et al., 2014). Also, a long-term ex situ conservation program which includes collecting seed from different sources and optimizing seed germination and plant establishment protocols is needed to restore disturbed habitats. Also, a supply of living plants for trade would prevent further extraction of plants from nature.

Acknowledgements

We thank to D. Lugo-Aquino, N. Pérez-Nasser and Ana Luisa Albarrán-Lara for technical assistance and suggested improvements to the manuscript. This project was supported by SDI-UNAM (Secretaria de Desarrollo Institucional, Universidad Nacional Autónoma de México, Octavo piso Torre de Rectoría, Ciudad Universitaria).

References
Anderson, F.E., Arias, M.S., Taylor, N.P., 1994. Threatened Cacti of Mexico. Royal Botanic Gardens Kew, England, p. 135.
Arias, S. , 1993. Cactáceas: conservación y diversidad en México. Revista de la Sociedad Mexicana de Historia Natural, Volumen especial XLIV, pp. 109-115.
Arias, S. , Gama, S. , Guzmán, U. , 1997. Las cactáceas del valle de Tehuacán-Cuicatlán, Flora del Valle de Tehuacán-Cuicatlán. Fascículo 14. Instituto de Biología, UNAM, México, D. F, p. 146.
Ayyad M.A., 2003. Case studies in the conservation of biodiversity: degradation and threats. J. Arid Environ, 54: 165-182. DOI:10.1006/jare.2001.0881
Belkhir, K. , Borsa, P. , Chikhi, L. , et al. , 1996-2004. GENETIX 4. 05, logiciel sous Windows TM pour la genetique des populations. Montpellier: Laboratoire Genome, Populations, Interactions. CNRS UMR 5171. Université de Montpellier Ⅱ.
Bravo-Hollis, H. , 1978. Las cactáceas de México, vol. Ⅰ. Universidad Nacional Autonoma de México, México, D. F, p. 743.
Bravo-Hollis, H. , Sánchez-Mejorada, H. , 1991. Las cactáceas de México, vol. Ⅱ. Instituto de Biología, UNAM, Mexico, D. F, p. 404.
Cascante A., Quesada M., Lobo J.J., et al, 2002. Effects of dry tropical forest fragmentation on the reproductive success and genetic structure of the tree Samanea saman. Conserv. Biol, 16: 137-147. DOI:10.1046/j.1523-1739.2002.00317.x
Castro-Felix P., Rosas-Espinoza V.C., Díaz-Cárdenas B., et al, 2013. Genetic diversity within a declining natural population of Ferocactus histrix (DC) Lindsay. Plant Species Biol, 29: E21-E30.
Clark-Tapia R., Alfonso-Corrado C., Eguiarte L.E., et al, 2005. Clonal diversity and distribution in Stenocereus eruca (Cactaceae), a narrow endemic cactus of the Sonoran Desert. Am. J. Bot, 92: 272-278. DOI:10.3732/ajb.92.2.272
CITES, 2013. Convención sobre el Comercio Internacional de Especies Amenazadas de Fauna y Flora Silvestres. Apéndice Ⅰ, Ⅱ y Ⅲ. http://www.cites.org/eng/app/appendices.php. última consulta: 2. Ⅳ 2014.
Contreras C., Valverde T., 2002. Evaluation of the conservation status of a rare cactus (Mammillaria crucigera) through the analysis of its population dynamics. J. Arid Environ, 51: 89-102. DOI:10.1006/jare.2001.0926
Cornejo-Romero A., Medina-Sánchez J., Hernández-Hernández T., et al, 2014. Quaternary origin and genetic divergence of the endemic cactus Mammillaria pectinifera in a changing landscape in the Tehuacán Valley Mexico. Genet. Mol. Res, 13: 73-88. DOI:10.4238/2014.January.8.6
Cornuet J.M., Luikart G., 1996. Description and power analysis of two tests for detecting recent population bottlenecks from allele frequency data. Genetics, 144: 2001-2014.
Dávila P., Arizmendi M.C., Valiente-Banuet A., et al, 2002. Biological diversity in the Tehuacán-Cuicatlán Valley, Mexico. Biodivers. Conservation, 11: 421-442. DOI:10.1023/A:1014888822920
Dieringer D., Schlötterer C., 2003. Microsatellite analyser (MSA): a platform independent analysis tool for large microsatellite data set. Mol. Ecol. Notes, 3: 167-169. DOI:10.1046/j.1471-8286.2003.00351.x
Doyle J.J., Dickson E.E., 1987. Preservation of plant samples for DNA restriction endonuclease analysis. Taxon, 36: 715-722. DOI:10.2307/1221122
Doyle J.J., Doyle J.L., 1987. A rapid DNA isolation procedure for small quantities of fresh leaf tissue. Phytochem. Bull, 19: 11-15.
Dyer R.J., Nason J.D., 2004. Population graphs: the graph theoretic shape of genetic structure. Mol. Ecol, 13: 1713-1727. DOI:10.1111/mec.2004.13.issue-7
Dyer R.J., Nason J.D., Garrick R.C., 2010. Landscape modelling of gene flow: improved power using conditional genetic distance derived from the topology of population networks. Mol. Ecol, 19: 3746-3759. DOI:10.1111/j.1365-294X.2010.04748.x
Earl D.A., Von Holdt B.M., 2011. STRUCTURE HARVESTER: a website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv. Genet. Resour, 4: 359-361.
Ellstrand N.C., Ellam D.R., 1993. Population genetic consequences of small population size: implications for plant conservation. Annu. Rev. Ecol. Syst, 24: 217-242. DOI:10.1146/annurev.es.24.110193.001245
England P.R., Cornuet J.M., Berthier P., et al, 2006. Estimating effective population size from linkage disequilibrium: severe bias in small samples. Conserv. Genet, 7: 303-308. DOI:10.1007/s10592-005-9103-8
Ewers R.M., Didham R.K., 2006. Confounding factors in the detection of species responses to habitat fragmentation. Biol. Rev, 81: 117-142.
Excoffier L., Laval G., Schneider S., 2005. Arlequin (version 3. 0): an integrated software package for population genetics data analysis. Evol. Bioinforma, 1: 47-50.
Excoffier L., Lischer H.E.L., 2010. Arlequin suite ver 3. 5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol. Ecol. Resour, 10: 564-567. DOI:10.1111/men.2010.10.issue-3
Evanno G., Regnaut S., Goudet S., 2005. Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol. Ecol, 14: 2611-2620. DOI:10.1111/mec.2005.14.issue-8
Falush D., Stephens M., Pritchard J.K., 2003. Inference of population structure using multilocus genotype data: linked loci and correlated allele frequencies. Genetics, 164: 1567-1587.
Figueredo C.J., Nassar J.M., García-Rivas A.E., 2010. Population genetic diversity and structure of Pilosocereus tillianus (Cactaceae Cereeae), a columnar cactus endemic to the Venezuelan Andes. J. Arid. Environ, 74: 1392-1398. DOI:10.1016/j.jaridenv.2010.05.020
Fischer R.A., Turner N.C., 1978. Plant productivity in the arid and semiarid zones. Annu. Rev. Plant Physiol, 29: 277-317. DOI:10.1146/annurev.pp.29.060178.001425
Frankham R., 2010. Inbreeding depression Inbreeding in the wild really does matter. Heredity, 104: 124. DOI:10.1038/hdy.2009.155
Ghazoul J., 2005. Pollen and seed dispersal among dispersed plants. Biol. Rev, 80: 413-443. DOI:10.1017/S1464793105006731
Glass, C. , 1998. Guía para la identificación de Cactáceas amenazadas de México, vol. 1. Comisión Nacional para el Estudio de la Biodiversidad/Cante, México, D. F, p. 210.
Gibson B.C., Nobel P.S., 1986. The Cactus Primer. London: Harvard University Press.
Goldstein D.B., Ruiz-Linares A., Cavalli-Sforza L.L., et al, 1995. An evaluation of genetic distances for use with microsatellite loci. Genetics, 139: 463-471.
Goverde M., Schweizer K., Baur B., et al, 2002. Small-scale habitat fragmentation effects on pollinator behavior: experimental evidence from the bumblebee Bombus veteranuson calcareous grasslands. Biol. Conserv, 104: 293-299. DOI:10.1016/S0006-3207(01)00194-X
Hamrick, J. L. , Nason, J. D. , Fleming, T. H. , 2002. Genetic diversity in columnar cacti. In: Fleming, T. , Valiente-Banuet, A. (Eds. ), Evolution, Ecology and Conservation of Columnar Cacti and Their Mutualists. Arizona University Press, pp. 122-133.
Hamrick J., 2004. Response of forest trees to global environmental changes. For Ecol Manag, 197: 323-335. DOI:10.1016/j.foreco.2004.05.023
Hart, D. L. , Clark, A. G. , 1997. Principles of Population Genetics, third ed. Sinauer Associates, Sunderland, MA.
Hernandez H.M., Godinez H., 1994. Contribution to the knowledge of endangered Mexican cacti. Acta Bot. Mex, 26: 33-52.
Hershkovitz M.A., Zimmer E.Z., 1997. On the evolutionary origins of the cacti. Taxon, 46: 217-232. DOI:10.2307/1224092
Hubisz M., Falush D., Stephens M., et al, 2009. Inferring weak population structure with the assistance of sample group information. Mol. Ecol. Resour, 9: 1322-1332. DOI:10.1111/men.2009.9.issue-5
Ibarra-Suárez, A. , 2009. Estudio de la diversidad genética de cactáceas endémicas del género Mammillaria del Valle de Tehuacán-Cuicatlán, Puebla-Oaxaca, UNAM, Tlalnepantla de Baz. Thesis (Biology). Facultad de Estudios Superiores Iztacala, p. 37.
Jamieson I.G., Allendorf F.W., 2012. How does the 50/500 rule apply to MVPs?. Trends Ecol. Evol., 27: 578-584. DOI:10.1016/j.tree.2012.07.001
López-Ortiz, N. , Salas-González, P. , Dávila, P. , et al. , 2013. Diversidad y estructura genética poblacional de Mammillaria zephyranthoidesScheidwer. In: 1841 (Cactaceae) una especie endámica de México, XIX Congreso Mexicano de Botánica (2013) Tuxtla Gutiérrez, October 20-25, 2013. Chiapas, Mexico.
Lowe A.D., Boshier M.W., Bacles C., et al, 2005. Genetic resource impacts of habitat loss and degradation; reconciling empirical evidence and predicted theory for neotropical trees. Heredity, 95: 255-273. DOI:10.1038/sj.hdy.6800725
Luikart G., Cornuet J.M., 1998. Empirical evaluation of a test for identifying recently bottlenecked populations from allele frequency data. Cons. Biol, 12: 228-237. DOI:10.1046/j.1523-1739.1998.96388.x
Macías-Arrastio, F. F. , Dávila, P. , Solórzano, S. , October 20-25, 2013. Diversidad y estructura genética de Mammillaria solisioides Backeb. , especie endémica de la Mixteca de Oaxaca y Puebla. XIX Congreso Mexicano de Botánica. Tuxtla Gutiérrez, Chiapas, Mexico.
Manni F., Guerard E., Heyer E., 2004. Geographic patterns of (genetic, morphologic, linguistic) variation: how barriers can be detected by ''Monmonier's algorithm''. Hum. Biol, 76: 173-190. DOI:10.1353/hub.2004.0034
Martorell C., Peters E., 2005. The measurement of chronic disturbance and its effects on the Threatened cactus Mammillaria pectinifera. Biol. Conserv, 124: 199-207. DOI:10.1016/j.biocon.2005.01.025
McNeely, J. A. , Miller, K. R. , Reid, W. V. , et al. , 1990. Conserving the World's Biological Diversity. World Conservation Union and World Resources Institute, Gland, Switzerland, and Washington, D. C.
Meyrán, G. J. , 1973. Guía botanica de cactáceas y otras suculentas del valle de Tehuacán. Sociedad Mexicana de Cactología A. C. , México, D. F, p. 50.
Moraes E.M., Abreu A.G., Andrade S.C.S., et al, 2005. Population genetic structure of two columnar cacti with a patchy distribution in eastern Brazil. Genetica, 125: 311-323. DOI:10.1007/s10709-005-0716-0
Moreira P.A., Fernandes G.W., Collevatti R.G., 2009. Fragmentation and spatial genetic structure in Tabebuia ochracea (Bignoniaceae) a seasonally dry Neotropical tree. For. Ecol. Manag, 258: 2690-2695. DOI:10.1016/j.foreco.2009.09.037
Nassar J.M., Hamrick J., Fleming T.H., 2001. Genetic variation and population structure of the mixed-mating cactus, Melocactus curvispinus (Cactaceae). Heredity, 87: 69-79. DOI:10.1046/j.1365-2540.2001.00910.x
Otero-Arniz A., Casas A., Hamrick J.L., 2005. Direct and indirect estimates of gene flow Among Populations of wild and managed Polaskia chichipe, an endemic columnar cactus in Central Mexico. Mol. Ecol, 14: 4313-4322. DOI:10.1111/j.1365-294X.2005.02762.x
Parra F., Pérez-Nasser N., Lira R., et al, 2008. Population genetics and process of domestication of Stenocereus pruinosus (Cactaceae) in the Tehuacán Valley, México. J. Arid Enviroments, 72: 1997-2010. DOI:10.1016/j.jaridenv.2008.06.007
Peters E.M., Arizaga S., Martorell C., et al, 2014. Geographic distribution and conservation status of Mammillaria pectinifera populations. Rev. Mex. Biodivers, 85: 942-952. DOI:10.7550/rmb.36338
Piry S., Luikart G., Cornuet J.M., 1999. BOTTLENECK: a computer program for detecting recent reductions in the effective population size using allele frequency data. J. Hered, 90: 502-503. DOI:10.1093/jhered/90.4.502
Pritchard J., Stephens M., Donnelly P., 2000. Inference of population structure using multilocus genotype data. Genetics, 155: 945-959.
Reed D.H., Frankham R., 2003. The correlation between population fitness and genetic diversity. Conserv. Biol, 17: 230-237. DOI:10.1046/j.1523-1739.2003.01236.x
Schneider, S. , Roessli, D. , Excoffier, L. , 2000. Arlequin: a software for population genetics data analysis, version 2. 000. Genetics and Biometry Laboratory, Department of Anthropology, University of Geneva, Geneva (Switzerland).
Semarnat, (Secretaría de Medio Ambiente y Recursos Naturales). , 2010. Norma Oficial Mexicana NOM-059-SEMARNAT-2010, Protecci ón ambiental-Especies nativas de México de flora y fauna silvestres-Categorías de riesgo y especificaciones para su inclusión, exclusión o cambio-Lista de especies en riesgo. Diario Oficial de la Federación, 30 de Diciembre de 2010, Segunda Sección. México.
Slatkin M., 1995. A measure of population subdivision based on microsatellite allele frequencies. Genetics, 139: 457-462.
Solorzano S., Cortés-Palomec A., Ibarra A., et al, 2009. Isolation, characterization and cross-amplification of polymorphic microsatellite loci in the Threatened endemic Mammillaria crucigera (Cactaceae). Resourses Mol. Ecol, 9: 156-158. DOI:10.1111/men.2009.9.issue-1
Solórzano S., Cuevas-Alducin P.D., García-Gómez V., et al, 2014. Genetic diversity and conservation of Mammillaria huitzilopochtli and M. supertexta, two threatened species endemic of the semiarid region of central Mexico. Rev. Mex. Biodivers, 85: 565-575. DOI:10.7550/rmb.39066
Tapia-Salcido, H. J. , 2011. Análisis de la diversidad y estructura genética poblacional de dos especies del género Mammillaria, endémicas del Valle de Tehuacán-Cuicatlán. Tesis de Maestría. Facultad de Estudios Superiores Iztacala, UNAM, México.
Valverde P.L., Zavala-Hurtado J.A., Jiménez-Sierra C., et al, 2009. Assessment of extinction risk of Mammillaria pectinifera, an endemic cactus of the Tehuacán-Cuicatlán region. Rev. Mex. Biodivers, 80: 219-230.
Valverde P.L., Zavala-Hurtado J.A., 2006. Assessing the ecological status of Mammillaria pectinifera Weber (Cactaceae), a rare and threatened species endemic of the Tehuacán-Cuicatlán Region in Central Mexico. J. Arid Environ, 64: 193-208. DOI:10.1016/j.jaridenv.2005.06.001
Van-Oosterhout C., Weetman D., Hutchinson W.F., 2005. Estimation and adjustment of microsatellite null allelles in nonequilibrium. Ecol. Notes, 6: 255-256.
Waples R.S., 2006. A bias correction for estimates of effective population size based on linkage disequilibrium at unlinked gene loci. Conserv. Genet, 7: 167-184. DOI:10.1007/s10592-005-9100-y
Waples R.S., Do C., 2008. LdNe: a program for estimating effective population size from data on linkage disequilibrium. Mol. Ecol. Resour, 8: 753-756. DOI:10.1111/j.1755-0998.2007.02061.x
Williams S.E., Hoffman E.A., 2009. Minimizing genetic adaptation in captive breeding programs: a review. Biol. Conserv, 142: 2388-2400. DOI:10.1016/j.biocon.2009.05.034
Young A., Clarke G., 2000. Genetics, demography and viability of fragmented populations. Ser. Conserv. Biol, 4: 600.
Zavala-Hurtado J.A., Valverde P.L., 2003. Habitat restriction in Mammilaria pectinifera, a threatened endemic Mexican cactus. J. Veg. Sci, 14: 891-898.
Zhang X.S., Jinliang W., Hill W.G., 2004. Redistribution of gene frequency and changes of genetic variation following a bottleneck in population size. Genetics, 167: 1475-1492. DOI:10.1534/genetics.103.025874
Hamrick J.L., Godt M.J.W., 1996. Effects of life history traits on genetic diversity in plant species. Phil Trans. Biol. Sci, 351: 1291-1298. DOI:10.1098/rstb.1996.0112