Am. J. Bot. Join the BSA
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Submit a response
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via ISI Web of Science (5)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Silva-Montellano, A.
Right arrow Articles by Eguiarte, L. E.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Silva-Montellano, A.
Right arrow Articles by Eguiarte, L. E.
Agricola
Right arrow Articles by Silva-Montellano, A.
Right arrow Articles by Eguiarte, L. E.
(American Journal of Botany. 2003;90:700-706.)
© 2003 Botanical Society of America, Inc.


Population Biology

Geographic patterns in the reproductive ecology of Agave lechuguilla (Agavaceae) in the Chihuahuan desert. II. Genetic variation, differentiation, and inbreeding estimates1

Arturo Silva-Montellano2 and Luis E. Eguiarte

Departamento de Ecología Evolutiva, Instituto de Ecología, Apartado postal 70-275, C.U., Universidad Nacional Autónoma de Mexico, CP 04510, D.F., Mexico

Received for publication July 5, 2002. Accepted for publication December 5, 2002.


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Plants with natural variation in their floral traits and reproductive ecology are ideal subjects for analyzing the effects of natural selection and other evolutionary forces on genetic structure of natural populations. Agave lechuguilla shows latitudinal changes in floral morphology, color, and nectar production along its distribution through north-central Mexico. Both the type and abundance of its pollinators also change with latitude. Using starch electrophoresis, we examined the levels and patterns of variation of 13 polymorphic allozyme loci in 11 populations of A. lechuguilla. The overall level of genetic variability was high (He = 0.394), but the levels of genetic variation had no geographic pattern. However, the southern populations exhibited an excess of heterozygotes in relation to expectations for Hardy-Weinberg equilibrium, whereas the northern populations had an excess of homozygotes. Total differentiation among populations was low ({theta} = 0.083), although gene flow estimates (Nm) varied among groups of populations: southern populations had the lowest levels of genetic differentiation, suggesting high levels of gene flow; northern populations had greater levels of genetic differentiation ({theta} = 0.115), suggesting low gene flow among them. The patterns and inferences of the genetic structure of the population at the molecular level is consistent with variation in floral traits and pollinator visitation rates across the range of the species.

Key Words: Agavaceae • Agave lechuguilla • Chihuahuan desert • gene flow • population differentiation • population genetics


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
The reproductive ecology and genetic structure of natural populations are deeply interconnected. The reproductive ecology of plants determines the amount of self- and cross-pollination as well as the distance that pollen and seeds move within populations, factors that ultimately influence the levels of inbreeding within populations (Schoen, 1982a , b ; Handel, 1983 ; Schoen and Clegg, 1985 ; Hamrick and Godt, 1989 ; Hamrick and Loveless, 1989 ; Eguiarte et al., 1993 ). Additionally, reproductive ecology can influence the amount of gene flow among populations, their effective population sizes, and the potential for natural selection to occur (Hamrick and Godt, 1989 ; Hamrick and Loveless, 1989 ; Eguiarte et al., 1993 ; Parra et al., 1993 ; Slatkin, 1994 ). On the other hand, genetic structure affects the evolution of reproductive traits, which in turn depends upon the availability of useful genetic variation and on effective population sizes (Hedrick, 1983 ; Parra et al., 1993 ). Although the literature on plant reproductive ecology and genetic structures is vast, few studies have explicitly studied the joint variation of reproductive ecology and genetic structure geographically (i.e., Schoen, 1982a , b ; Ohara et al., 1996 ; Wong and Sun, 1999 ).

Studies that examine differences in the reproductive ecology and genetic structure of several populations with contrasting environmental conditions serve to corroborate predictions about the effect of variation in ecological characteristics or selfing-outcrossing rates on molecular genetic diversity. For instance, Schoen (1982a , b ) studied Gilia achilleifolia in California and found correlations between pollinator identity, reproductive biology traits, outcrossing rates, and levels of genetic variation among several populations. Similar patterns were found in Trillium kamtschaticum by Ohara et al. (1996) , in which variation in the breeding system and genetic diversity of populations is related to the abundance of different pollinators in distinct areas on the island of Hokkaido, Japan.

We studied the relationship between reproductive ecology and genetic structure in several populations of Agave lechuguilla along a latitudinal gradient in the Mexican portion of the Chihuahuan desert. The reproductive ecology of the species was reported in a separate paper (Silva-Montellano and Eguiarte, 2003 ). This work demonstrated that along a latitudinal gradient, floral traits of the species clearly vary: southern populations have long, tubular, pale flowers, which produce large amounts of diluted nectar and are intensively visited by nocturnal and diurnal pollinators. In contrast, northern populations have shorter, more open reddish flowers, which produce smaller quantities of concentrated nectar and are less intensely visited (Cadaval, 1999 ; Silva-Montellano, 2001 ). The efficiency of fruit production is also higher in southern than in northern populations (Silva-Montellano and Eguiarte, 2003 ). Here we analyze the genetic structure of A. lechuguilla along its distributional range and investigate possible correlations between interpopulation differences in floral biology and differences in levels of genetic variation, inbreeding, and differentiation among populations. Because flowers are less frequently visited in the north, we expected to see higher levels of inbreeding and greater genetic differentiation in northern populations than in southern ones.


    MATERIALS AND METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Study species
Agave lechuguilla Torr. (Agavaceae, subgenus Littaea, Marginateae group; Gentry, 1982 ) is one of the most abundant Agave species. This small century plant is abundant throughout the Chihuahuan desert. Silva-Montellano and Eguiarte (2003) and Cadaval (1999) describe the general characteristics and floral biology of A. lechuguilla.

Studied populations
We studied 11 populations along a latitudinal gradient that extends through most of the Mexican portion of the Chihuahuan desert (Table 1). These populations did not statistically differ in density of rosettes and inflorescences as described in detail in Silva-Montellano and Eguiarte (2003) . Samples for genetic analyses were collected from June to September in 1996. Forty different individuals were collected in each population from an area of approximately 2 ha. To avoid repeatedly sampling members of clones, we collected tissues from individuals that were more than 2 m apart (see Trame et al., 1995 ). We sampled the younger expanded leaf of each rosette and collected a small portion from the base of the rosette. Samples were transported in liquid N2 and stored at –80°C until used.


View this table:
[in this window]
[in a new window]
 
Table 1. Location and genetic diversity for the studied populations of Agave lechuguilla and number of loci in each population that differs significantly from Hardy-Weinberg equilibria. N = mean number of plants assayed per locus (SE); A = mean number of alleles per locus; Ae = effective number of alleles per locus; P = proportion of polymorphic loci (most common alelle); He = expected heterozygosity (SE; unbiased estimator; Nei, 1978); Ho = observed heterozygosity (SE). H-W: + = positively different; – = negatively different

 
Allozyme procedures
We followed standard procedures for starch-gel electrophoresis (Soltis et al., 1983 ; Martínez-Palacios et al., 1999 ). Samples were mechanically disrupted with an electric drill, using 15–20 drops of the extraction buffer described in Aguirre-Planter et al. (2000) . Extracts were absorbed in 12 x 1.5 mm paper wicks and stored at –80°C. Allozymes were run in 12% starch gels for 7 h at 60 mA, using buffer system 8 of Soltis et al. (1983) with gel buffer at pH 7.6. We analyzed 10 enzymes, giving a total of 13 loci: phosphoglucose isomerase (PGI, E.C. [Enzyme Commission number] 5.3.1.9, two loci), glutamate oxaloacetate transaminase (GOT, E.C. 2.6.1.1, one locus), leucine aminopeptidase (LAP, E.C. 3.4.11.1, one locus), esterase (EST, E.C. 3.1.1.1, one locus), diaphorase (DIA, E.C. 1.6.4.3, two loci), malic enzyme (ME, E.C.1.1.1.40, one locus), phosphoglucomutase (PGM, E.C. 2.7.5.1, two loci), acid phosphatase (ACPH, E.C.3.1.3.2, one locus), anodic peroxidase (APX, E.C. 1.11.1.7, one locus), and cathodic peroxidase (CPX, E.C. 1.11.1.7, one locus). Staining procedures correspond to those of Soltis et al. (1983) , with minor modifications available from the authors.

The fastest loci and alleles were scored as 1, the second 2, and so on, following genetic interpretations of the closely related species, Agave victoriae-reginae (Martínez-Palacios et al., 1999 ) for the nine loci common to both species. One population of A. victoriae-reginae (population 7 from Martínez-Palacios et al., 1999 ) was used in the analysis both as an outgroup (see below) and for comparison of the loci and alleles, including the 13 loci of the present study.

Statistical analyses
Most of the statistical genetic analyses were done using TFPGA (Miller, 1997 ), unless indicated otherwise. We obtained the allelic frequencies for each locus, and from these frequencies, we calculated the proportion of polymorphic alleles (P), the observed (Ho) and expected (He) heterozygosities, the average number of alleles per locus (A) and the effective number of alleles (Ae) (Hedrick, 1983 ). To analyze the deviations from Hardy-Weinberg equilibrium, we estimated the fixation index (f) for each locus from each population with BIOSYS (Swofford and Selander, 1989 ), using a {chi}2 test to evaluate if deviations were different from zero (Li and Horvitz, 1953 ). Clinal variation of expected and observed heterozygosities per locus and fixation index per locus along the latitudinal gradient were tested using regression for repeated measures of Y (Sokal and Rohlf, 1995 ), which includes an ANOVA for differences among populations.

Wright's F statistics were obtained following Weir and Cockerham (1984 ; F, f, and {theta}, equivalent to FIT, FIS, and FST, respectively) procedures. Single locus values for each index were tested if different from zero using {chi}2 tests (Li and Horvitz, 1953 ; Workman and Niswander, 1970 ) and averaged by means of a jackknife procedure (Weir, 1990 ). Confidence intervals of 95% for each statistic were obtained from 1000 bootstrap samples. Nm was aproximated from {theta} using Crow and Aoki's (1984) formula.

We used a Mantel test (Manly, 1987 ) to assess the model of isolation by distance using the genetic distance for pairs of populations (Nei, 1978 ) and geographic distance among these populations. The Nei's genetic distance was also employed to obtain a UPGMA phenogram, after 1000 bootstrap samples.


    RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Levels of genetic variation and inbreeding within populations
All analyzed loci in all populations segregated in a similar patterns to those observed in A. victoriae-reginae, facilitating genotype scoring, interpretation, and comparisons. Allelic frequencies are available from the authors upon request. All loci (Table 1) were polymorphic in most (six) populations, resulting in very high levels of genetic variation. The average number of alleles per locus (A) was 2.28, the effective number of alleles per locus (Ae) was 1.72, and expected heterozygosity (He) was 0.394. Although the three southernmost populations had the highest genetic variation (Table 1; He range 0.449–0.457), populations had no statistical differences in He (ANOVA, N = 143, df = 10, F = 1.7, P = 0.087). In contrast, we found significant differences in observed heterozygosity (Ho) among populations (ANOVA, N = 143, df = 10, F = 4.56, P < 0.0001). Observed heterozygosity decreased significantly along the latitudinal gradient (Fig. 1; N = 143, R2 = 0.193, F = 33.7, P < 0.0001; slope = –0.03, t = –5.8, P < 0.0001).



View larger version (13K):
[in this window]
[in a new window]
 
Fig. 1. Observed heterozygosity (Ho, means ± 1 SE) for 13 loci along a latitudinal transect of 11 populations of Agave lechuguilla in the Chihuahuan desert in Mexico. The numbers above the symbol indicate the name of the population. Different letters below symbols indicate that the populations are statistically different (P < 0.05), contrast t test: A–B, t > 2.2, P < 0.03; A–BC, t > 3.5, P < 0.0007; A–C, t > 3.5, P < 0.0005; B–C, t > 2.1, P < 0.04. Regression formula: y = –0.03x + 1.12, R2 = 0.193, where slope has t = 6.0, P < 0.0001

 
The interpopulation differences in Ho are due to deviations from Hardy-Weinberg equilibrium, as shown by the average fixation indices per loci (Fig. 2). In the southernmost populations (1–3), several loci had significantly negative f values (Table 1; Fig. 2). In contrast, in the north several loci had significantly positive f values (Table 1, populations 5–10, including M). The average f index per loci differed among populations (ANOVA, N = 143, df = 10, F = 4.65, P < 0.0001) and increased with latitude (Fig. 2; N = 143, R2 = 0.203, F = 36.01, P < 0.0001; slope = 0.056, t = 6.0, P < 0.0001). The value of f was highest in population 8 (with eight significantly positive loci, average f = 0.542). This latitudinal pattern indicates that populations in the south (1–3) had an excess of heterozygotes compared to expected Hardy-Weinberg allelic frequencies. Populations 4 and 5 were near Hardy-Weinberg equilibrium, and in the north, several populations had an excess of homozygous individuals, particularly in populations 6, 8, 9, and 10, suggesting high levels of inbreeding.



View larger version (14K):
[in this window]
[in a new window]
 
Fig. 2. Population fixation indices (f, means ± 1 SE; for 11 populations, 13 loci) of Agave lechuguilla along the Chihuahuan desert in Mexico. The number above the symbol indicates the name of the population, and an asterisk (or a double asterisk) indicates the population is statistically different from zero ({chi}2 test = *P < 0.05, **P < 0.001). Different letters below symbols indicate that the populations are statistically different, according to contrast t test: A–B, t > 2.1, P < 0.04; A–BC, t > 3.2, P < 0.002; A–C, t > 4.9, P < 0.00001; B–C, t > 2.4, P < 0.02. Regression formula: y = 0.056x + 1.28, R2 = 0.203, where slope has t = –5.8, P < 0.0001

 
Weir and Cockerham F statistics
The mean overall f (=FIS) of 0.105 (Table 2) was not statistically different from zero. However there was variance in f among loci, and some loci exhibited significant levels of imbreeding (f > 0). {theta} (=FST) was relatively small (0.083), but it was significantly different from zero, indicating low but significant genetic differentiation among populations. As a consequence, F (=FIT) was significantly positive (0.179).


View this table:
[in this window]
[in a new window]
 
Table 2. Weir and Cockerham F statistics (1984), estimates of Nm (Crow and Aoki, 1984) per locus for all sampled populations of Agave lechuguilla, and number of populations for each locus that differs significantly from Hardy-Weinberg equilibria. Total and partial (mean) for the populations from south (1, 2, and 3), center (4, 5, 6, and M), and north (7, 8, 9, and 10) of its distribution. SD = standard deviation, jackknife; 95% CI = confidence intervals at 95%, bootstrap = 1000 iterations; all single locus estimates are significantly different from zero (P < 0.01) unless stated "ns." H-W: + = positively different; – = negatively different

 
To understand the patterns of inbreeding and differentiation, we divided the 11 populations into three geographic groups (Table 2): south (1–3), central (4–6, including M), and north (7–10). The mean value of f in the southern populations was –0.175, significantly less than zero. The mean f value of the central populations was 0.132, but it was not significantly different from zero. The northern population had a much higher mean value f at 0.355 that was greater than zero.

We detected no genetic differentiation among populations in the southern populations ({theta} = 0.009, not different from zero, according to the 95% confidence interval), but central populations had a significant amount of differentiation ({theta} = 0.059, different from zero, 95% CI) and an even greater differentiation among the northern ones ({theta} = 0.115, different from zero, 95% CI).

The single locus estimates of the mean number of migrants per generation (Nm) was always higher than one (Table 2; average Nm = 2.28), suggesting that gene flow is probably an important force in A. lechuguilla. Considering each of our geographic regions separately, Nm declined from south to north. The mean Nm was 12.23 in the south, suggesting high levels of gene flow among populations, 2.23 in the central populations, and 1.08 in the northern ones.

Allelic frequencies along the latitudinal gradient and distribution of rare alleles
The most common allele in each locus was plotted as a function of population latitude. In most cases, we were not able to detect clear patterns in the changes of allelic frequencies. Only two loci (GOT, 3 and APX, 1) had marginally significant (P {approx} 0.05) latitudinal increases or decreases in allelic frequencies (not shown).

Rare alleles (with an allelic frequency less than 0.05 within a given population) were more abundant in some populations (Fig. 3; 1, 6, M, 7, and 8), than in others (2, 3, 4, 5, 9, and 10). The frequency of rare alleles is negatively related to gene flow (Slatkin, 1985 ; Slatkin and Barton, 1989 ). Hence, the low amount of private alleles (alleles detected in only one population; one at population 8, locus DIA1-3) is congruent with the overall low {theta} and high gene flow estimated previously.



View larger version (34K):
[in this window]
[in a new window]
 
Fig. 3. Number of rare alleles (total number of alleles = 28; allelic frequency of 0.05 or less within populations) for 11 populations of Agave lechuguilla along the Chihuahuan desert in Mexico

 
Genetic distances and phenogram
We used Nei's genetic distance to explore the genetic relationships among populations (see Table 3). The average distance between all pairs of populations was low (D = 0.07). It was highest between populations 5 and 10 (D = 0.17) and lowest between populations 2 and 3 (D = 0.007). A Mantel test indicated that the relationship between the geographic (Table 3) and genetic distances was not significant (10 000 permutations; z = 2176.2, r = 0.22, P = 0.11). The reason for this lack of significance becomes clear after a UPGMA analysis using A. victoriae-reginae as an outgroup (Fig. 4). Southern populations (1–3) form a clear group (average D = 0.064), and some central populations were clustered together (4, 5, and 6, average D = 0.035). But the remaining central–northern populations have little congruence between geographic and genetic distance. The bootstrap values were sometimes low (groups 4–8 and M–9), as would be expected given the low genetic differentiation and high levels of gene flow among populations. It is noteworthy that population 10 was genetically the most differentiated (average to all the populations D = 0.13). A previous study also found that this population was divergent in general morphology and reproductive characters (Silva-Montellano and Eguiarte, 2003 ).


View this table:
[in this window]
[in a new window]
 
Table 3. Genetic (Nei, 1978) and geographic distances (in kilometers) for all pairs of sampled populations of Agave lechuguilla in the Chihuahuan desert of Mexico

 


View larger version (9K):
[in this window]
[in a new window]
 
Fig. 4. UPGMA phenogram constructed using Nei's (1978) genetic distance of 11 populations of Agave lechuguilla in the Chihuahuan desert in Mexico and the A. victoriae-reginae population as an outgroup (Martinez-Palacios et al., 1999 ). Numbers above branches indicate genetic distance, and numbers in italics are the proportion of times the bifurcation was found in 10 000 bootstrap samples

 

    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Genetic variation
Levels of genetic variation in A. lechuguilla populations are high on average (He range 0.33–0.46) relative to those of other long-lived perennials (Eguiarte, 1990 ; Hamrick et al., 1992 ; i.e., average of long-lived perennial outcrossers at the species level He = 0.211, N = 51; and at the population level He = 0.163, N = 66). Similarly high levels of He were also reported in the closely related but more geographically restricted species A. victoriae-reginae (Martínez-Palacios et al., 1999 ; He range 0.21–0.41) and in species of the A. deserti complex, using RAPDs (Navarro-Quezada, 1999 ; Gonzalez Chauvet, 2000 ; He range 0.12–0.26). Our prediction that southern populations would harbor greater genetic diversity than northern ones was met to some extent because southern populations had higher He levels, although the change along the latitudinal gradient was not significant.

The presence of similarly high levels of genetic variation in A. victoriae-reginae and in A. lechuguilla, despite the more restricted geographic range and smaller population sizes of A. victoriae-reginae, suggests that Agave species may harbor particularly high levels of genetic variation compared to other long-lived desert perennials (Martínez-Palacios et al., 1999 ; Silva-Montellano, 2001 ). Long-lived desert perennials have variable levels of genetic variation, ranging from very low, as in the palm Washingtonia filifera (McClenaghan and Beauchamp, 1986 ) to quite high, as in the cactus Echinocereus engelmannii (Neel et al., 1996 ) and other Cactaceae (Harmick et al., 2002 ), as well as in species of different families analyzed by Schuster et al. (1994 ; average He = 0.173). Nevertheless, the high levels of genetic variation in wild Agave species is not typical of all desert plants (see Comes and Abbot, 1999 ). In contrast, domesticated Agave species harbor very low levels of genetic variation (Colunga-GarcíaMarín et al., 1999 ; Gil-Vega et al., 2001 ).

Fixation indices and inbreeding within populations
Our data revealed a latitudinal trend in the deviation of genotypic frequencies from Hardy-Weinberg equilibrium. Latitude accounted for 20% of the variance in the regression of heterozygosity across latitude, decreasing from south to north. Populations in the south (1–3) had an excess of heterozygous individuals, whereas populations in the central section (4–6, including M) had genotypic frequencies near Hardy-Weinberg equilibrium and those in the north had an excess of homozygous individuals. This pattern is consistent with the higher rates of pollinator visitation and greater fruit set in the south (Silva-Montellano and Eguiarte, 2003). Lower rates of pollinator visitation may explain the lower efficiencies in fruit production and higher levels of homozygosity in northern populations. The excess of heterozygous individuals in the south could be generated by several different mechanisms (Hedrick, 1983 ; Waser, 1987 ). The most important of these is probably heterosis: more heterozygous individuals could mean greater vigor and survival (Linhart et al., 1981 ; Eguiarte et al., 1992 ). This idea is consistent with the fact that A. lechuguilla is a long-lived perennial (Eguiarte et al., 2000 ).

The low levels of genetic differentiation in the southern populations suggest very high levels of gene flow and behavior as a single panmictic population. In these very large populations, natural selection favoring the heterozygous individuals would be more efficient (Eguiarte, 1990 ; Parra et al., 1993 ). Meanwhile, in the north the situation is the reverse; populations are more isolated and show great genetic distance. Flowers in northern populations receive fewer pollinator visits and apparently generate inbreeding by geitonogamous self-pollination among flowers within a given inflorescence or by moving pollen among the different inflorescences in the same genet (Silva-Montellano, 2001 ). If we assume that all of the inbreeding in the northern populations (7–10) is due to self-pollination, the outcrossing rate can be approximated as f = 1 – t/1 + t (Haldane, 1924 ; Hedrick, 1983 ). Using this equation yields a t value of 0.467, which suggests that about half of all seeds are produced by self-pollination. The finding that outcrossing is higher in the southern populations and lower in the north was supported by a subsequent series of controlled pollination experiments and genetic analyses of outcrossing-rates between the two most contrasting populations (Silva-Montellano, 2001 ).

Genetic differentiation
We were able to detect genetic differentiation in populations of A. lechuguilla along a latitudinal transect in the Chihuahuan desert. The pattern of population differentiation along this transect was congruent with patterns of morphological and reproductive differentiation found in a previous study (Silva-Montellano, 2001 ; Silva-Montellano and Eguiarte, 2003 ). The average {theta} (0.083) found in A. lechuguilla was significantly different from zero and was lower than the value for Agave victoriae-reginae ({theta} = 0.236). The latter species is found only in 10 very localized populations in a relatively small area of the central Chihuahuan desert. These differences in {theta} may be caused by the large population sizes and wide distribution of A. lechuguilla. On the other hand, low levels of interpopulation genetic differentiation were detected in A. deserti ({theta} in this species complex ranged from 0.08 to 0.13, Navarro-Quezada, 1999 ; Gonzalez Chauvet, 2000 ). Agave deserti is relatively abundant and widespread in the Sonoran desert. We also found strong patterns of differentiation between northern and southern populations. Southern populations of A. lechuguilla were genetically similar, central populations were more differentiated (mostly due to the M population), whereas northern populations were highly differentiated. This phenomenon can be explained by differences in pollinator abundance and in taxa across the gradient. Pollinators seem to be less abundant and to have a more restricted range of movement in the north than in the south (Silva-Montellano, 2001 ; Silva-Montellano and Eguiarte, 2003 ). The movements of more abundant and more mobile pollinators would generate greater genetic connectivity among southern than among northern populations.

The UPGMA analysis suggests that the more primitive populations may be those in the north (in particular population 10). Although we are aware that some other evolutionary process may account for the observed patterns (i.e., intense genetic drift or introgression from other Agave species in some of the northern populations), our interpretation is worth considering, because it suggests a contrast to the traditional view that the Agave species originated in central Mexico, and then migrated to the north (i.e., Alvarez de Zayas, 1989 ). If our interpretation is correct, then A. lechuguilla seems to have originated in the north and more recently colonized the south. A similar pattern was detected using RAPDs in the A. deserti complex of the Sonoran desert, comprising several subspecies of A. deserti and of A. cerulata and A. subsimplex (Navarro-Quezada, 1999 ; Gonzalez Chauvet, 2000 ). If this scenario is correct, then the traits of the southern populations (larger flowers and higher outcrossing) are derived. An alternative hypothesis is that the increment in aridity to the north (Ortega, 1995 ) has led to a decline in nocturnal pollinators in the north (Silva-Montellano, 2001 ). Hence, the southern population may represent the ancestral traits before the changes occurred. In that way, northern populations have had more time to adapt to diurnal pollinators with a consequent increase in inbreeding and genetic differentiation among them. In support of this second point, while the common pollination syndrome for genus Agave is considered to be chiropterophily (Schaffer and Schaffer, 1977 ; Howell, 1979 ; Howell and Roth, 1981 ), we have evidence that the pollination syndrome in northern populations could be evolving toward diurnal pollination (Silva-Montellano and Eguiarte, 2003 ). Unfortunately, the present information of biogeography and phylogeny of the genus Agave is not enough to help solve this point (see Eguiarte et al., 2000 ). We need more detailed experiments and observations in the field to corroborate these interpretations, along with analyses of molecular evolution and population genetics of the Marginateae group of the Agave species.

In conclusion, our study of Agave lechuguilla (considering both this paper and Silva-Montellano and Eguiarte, 2003 ) clearly demonstrates that differences in pollinators not only affect the reproductive efficiency of plant populations, but also determine their genetic structure. As we already mentioned, few studies have attempted this joint approach using ecologic-genetic markers over a wide geographic range to unravel the evolutionary ecology of plant species, despite its huge potential.


    FOOTNOTES
 
1 The authors thank A. Cadaval and A. Valera for extraordinary help in the field; H. Caballero, A. Marínez-Palacios, and E. Aguirre for technical help in the laboratory; and J. Golubov, M. Rocha, F. Molina, C. Martínez del Rio, S. Good, T. Culley, and an anonymous reviewer for comments on an earlier version of the manuscript. Our research was supported by projects IN211997, Papiit-DGAPA, UNAM, and CONACyT 27983N to L.E.E. and by a DGE, UNAM scholarship to A.S.-M. The manuscript was written during a sabbatical leave to L.E.E. at the University of California at Irvine, with support from DGAPA, UNAM, and CONACyT scholarships. Back

2 Author for reprint requests (asilva{at}miranda.ecologia.unam.mx ) Back


    LITERATURE CITED
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Aguirre-Planter E. G. R. Furnier L. E. Eguiarte 2000 Low levels of genetic variation within and high levels of genetic differentiation among population of species of Abies from southern Mexico and Guatemala. American Journal of Botany 87: 362-371[Abstract/Free Full Text]

Alvarez De Zayas A. 1989 Distribución geográfica y posible origen de las Agavaceae. Revista del Jardín Botánico Nacional 10: 25-36

Cadaval A. 1999 Estudio evolutivo de los azúcares del néctar de Agave lechuguilla en el desierto de Chihuahua. Tesis de Licenciatura, Facultad de Ciencias, Universidad Nacional Autónoma de México, Ciudad de México, Distrito Federal, México

Colunga-Garcíamarín P. J. Coello-Coello L. E. Eguiarte D. Piñero 1999 Isoenzymatic variation and phylogenetic relations between henequen Agave furcroydes Lem. and its wild ancestor A. angustifolia Haw. American Journal of Botany 86: 115-123[Abstract/Free Full Text]

Comes H. R. Abbott 1999 Population genetic structure and gene flow across arid versus mesic environments: a comparative study of two parapatric Senecio species from the Near East. Evolution 53: 36-54

Crow J. K. Aoki 1984 Group selection for a phylogenetic behavioral trait: estimating the degree of population subdivision. Proceedings of the National Academy of Science, USA 81: 6073-6077[Abstract/Free Full Text]

Eguiarte L. E. 1990 Genética de poblaciones de Astrocaryum mexicanum Liebm. en Los Tuxtlas, Veracruz. Tesis Doctoral, Centro de Ecología, Universidad Nacional Autónoma de México, Ciudad de México, Distrito Federal, México

Eguiarte L. E. A. Búrquez J. Rodríguez M. Martínez-Ramos J. Sarukhán D. Piñero 1993 Direct and indirect estimates of neighborhood and effective population size in a tropical palm, Astrocaryum mexicanum. Evolution 47: 75-87[CrossRef][ISI]

Eguiarte L. E. N. Perez-Nasser D. Piñero 1992 Genetic structure, outcrossing rate and heterosis in Astrocaryum mexicanum (tropical palm): implications for evolution and conservation. Heredity 69: 217-228[ISI]

Eguiarte L. E. V. Souza A. Silva-Montellano 2000 Evolución de la familia Agavaceae: filogenia, ecología evolutiva de la reproducción y genética de poblaciones. Boletín de la Sociedad Botánica de México 66: 131-150

Gentry H. S. 1982 Agaves of continental North America. University of Arizona Press, Tucson, Arizona, USA

Gil-Vega K. M. Gonzalez-Chavira O. Martínez De La Vega J. Simpson G. Vandemark 2001 Analysis of genetic diversity in Agave tequilana var. Azul using RAPD markers. Euphytica 119: 335-341[CrossRef][ISI]

Gonzalez Chauvet R. 2000 Análisis de variación genética de Agave deserti en el desierto sonorense por medio de marcadores moleculares (RAPDs). Tesis de Licenciatura, Facultad de Ciencias, Universidad Nacional Autónoma de México, Ciudad de México, Distrito Federal, México

Haldane J. B. 1924 A mathematical theory of natural and artificial selection. II. Proceedings of the Cambridge Philosophical Society 23: 19-41[ISI]

Hamrick J. L. M. Godt 1989 Allozyme diversity in plant species. In A. Brown, M. Clegg, A. Kahler, and B. Weir [eds.], Plant population genetics, breeding and genetic resources, 43–63. Sinauer, Sunderland, Massachusetts, USA

Hamrick J. L. M. Godt S. Sherman-Broyles 1992 Factors influencing levels of genetic diversity in woody plant species. In W. T. Adams, S. H. Strauss, D. L. Copes, and A. R. Griffin [eds.], Population genetics of forest trees, 95–124. Kluwer Academic, London, UK

Hamrick J. L. M. Loveless 1989 The genetic structure of tropical tree population: association with reproductive biology. In J. Bock and Y. Linhart [eds.], The evolutionary ecology of plants, 129–146. Westview Press, Boulder, Colorado, USA

Hamrick J. L. J. D. Nason T. H. Fleming J. M. Nassar 2002 Genetic diversity in columnar cacti. In T. H. Fleming and A. Valiente-Banuet [eds.], Evolution, ecology and conservation of columnar cacti and their mutualists, 122–133. University of Arizona Press, Tucson, Arizona, USA

Handel S. 1983 Pollination ecology, plant population structure, and gene flow. In L. Real [ed.], Pollination biology, 163–211. Academic Press, Orlando, Florida, USA

Hedrick P. W. 1983 Genetics of populations. Science Books, Boston, Massachusetts, USA

Howell D. J. 1979 Flock foraging in nectar-feeding bats: advantages to the bats and the host plants. American Naturalist 144: 23-49

Howell D. J. B. Roth 1981 Sexual reproduction in agaves: the benefits of bats; the cost of semelparous advertising. Ecology 62: 1-7[CrossRef][ISI]

Li C. D. Horvitz 1953 Some methods of estimating the inbreeding coefficient. American Journal of Human Genetics 107: 511-523

Linhart Y. B. J. B. Mitton K. B. Sturgeon M. L. Davis 1981 Genetic variation in space and time in a population of Ponderosa pine. Heredity 46: 407-426[ISI]

Manly B. 1987 The statistics of natural selection on animal populations. Chapman and Hall, New York, New York, USA

Martínez-Palacios A. L. E. Eguiarte G. R. Furnier 1999 Genetic diversity of endangered endemic Agave victoriae-reginae (Agavaceae) in the Chihuahuan desert. American Journal of Botany 86: 1093-1098[Abstract/Free Full Text]

McClenaghan L. A. Beauchamp 1986 Low genetic differentiation among isolated populations of California fan palm (Washingtonia filifera). Evolution 40: 315-322[CrossRef][ISI]

Miller M. 1997 Tools for population genetic analyses (TFPGA) 1.3. Computer software. Distributed by the author. Flagstaff, Arizona, USA

Navarro-Quezada A. 1999 Estructura genética y procesos de especiación de Agave cerulata (Trel.) y A. subsimplex (Trel.) en el desierto sonorense a partir de RAPD's. Tesis de Licenciatura, Facultad de Ciencias, Universidad Nacional Autónoma de México, Ciudad de México, Distrito Federal, México

Neel M. J. Clegg N. Ellstrand 1996 Isozyme variation in Echinocereus engelmannii var. Munzii (Cactaceae). Conservation Biology 10: 622-631[CrossRef][ISI]

Nei M. 1978 Estimation of average heterozygosity and genetic distance from a small number of individuals. Genetics 89: 583-590[Abstract/Free Full Text]

Ohara M. H. Takeda Y. Ohno Y. Shimamoto 1996 Variations in the breeding system and the population genetic structure of Trillium kamtschaticum (Liliaceae). Heredity 76: 476-484[ISI]

Ortega J. R. 1995 Los paleoambientes holocénicos de la laguna de Babícora, Chihuahua, México. Geofísica Internacional 34: 107-116

Parra V. C. F. Vargas L. E. Eguiarte 1993 Reproductive biology, pollen and seed dispersal, and neighborhood size in the hummingbird-pollinated Echeveria gibbiflora (Crassulaceae). American Journal of Botany 80: 153-159[CrossRef][ISI]

Schaffer W. V. Schaffer 1977 The reproductive biology of Agavaceae I: pollen and nectar production in four Arizona agaves. Southwestern Naturalist 22: 157-168

Schoen D. 1982a The breeding system of Gilia achileifolia: variation in floral characteristics and outcrossing rate. Evolution 36: 352-360[CrossRef][ISI]

Schoen D. 1982b Genetic variation and the breeding system of Gilia achileifolia. Evolution 36: 361-370[CrossRef][ISI]

Schoen D. M. T. Clegg 1985 The influencing of flower color on outcrossing rate and male reproductive success in Ipomoea purpurea. Evolution 39: 1242-1249[CrossRef][ISI]

Schuster W. S. D. R. Sandquist S. L. Phillips J. R. Ehleringer 1994 High levels of genetic variation in populations of four dominant arid land plant species in Arizona. Journal of Arid Environments 27: 159-167

Silva-Montellano A. L. E. Aguiarte 2003 Geographic patterns in the reproductive ecology of Agave lechuguilla in the Chihuahuan desert. I. Floral characteristics, visitors, and fecundity. American Journal of Botany 90: 377-387[Abstract/Free Full Text]

Silva-Montellano J. A. 2001 Ecología reproductiva y genética de poblaciones de Agave lechuguilla (Torr.) en un gradiente latitudinal. Tesis Doctoral, Instituto de Ecología, Universidad Nacional Autónoma de México, Ciudad de México, Distrito Federal, México

Slatkin M. 1985 Rare alleles as indicators of gene flow. Evolution 39: 53-65[CrossRef][ISI]

Slatkin M. 1994 Gene flow and population structure. In L. Real [ed.], Ecological genetics, 4–17. Princeton University Press, Princeton, New Jersey, USA

Slatkin M. N. H. Barton 1989 A comparison of three indirect methods for estimating average levels of gene flow. Evolution 43: 1349-1368[CrossRef][ISI]

Sokal R. F. J. Rohlf 1995 Biometry, 3rd ed. W. H. Freeman, New York, New York, USA

Soltis D. C. Haufler D. Darrow G. Gastony 1983 Starch gel electrophoresis of ferns: a compilation of grinding buffers, gel and electrode buffers, and staining schedules. American Fern Journal 73: 9-27[CrossRef][ISI]

Swofford D. R. Selander 1989 BIOSYS-1: a computer program for the analysis of allelic variation in population genetic and biochemical systematics. Release 1.7. Natural History Survey, Champaign, Illinois, USA

Trame A. M. A. Coddinngton K. Paige 1995 Field and genetic studies testing optimal outcrossing in Agave schottii, a long-lived clonal plant. Oecologia 104: 93-100[CrossRef][ISI]

Waser N. M. 1987 Spatial genetic heterogeneity in a population of the montane perennial plant Delphinium nelsonii. Heredity 58: 249-256[ISI]

Weir B. S. 1990 Genetic data analysis. Sinauer, Sunderland, Massachusetts, USA

Weir B. S. C. Cockerham 1984 Estimating F-statistics for the analysis of population structure. Evolution 38: 1358-1370[CrossRef][ISI]

Wong K. C. M. Sun 1999 Reproductive biology and conservation genetics of Goodyera procera (Orchidaceae). American Journal of Botany 86: 1406-1413[Abstract/Free Full Text]

Workman P. J. Niswander 1970 Population studies on southwestern Indian tribes II. Local genetic differentiation in the Papago. American Journal of Human Genetics 22: 24-49[ISI][Medline]




This article has been cited by other articles:


Home page
Am. J. Bot.Home page
K. C. Parker, J. L. Hamrick, W. C. Hodgson, D. W. Trapnell, A. J. Parker, and R. K. Kuzoff
Genetic consequences of pre-Columbian cultivation for Agave murpheyi and A. delamateri (Agavaceae)
Am. J. Botany, September 1, 2007; 94(9): 1479 - 1490.
[Abstract] [Full Text] [PDF]


Home page
ANN BOT (LOND)Home page
A. E. ASSOGBADJO, T. KYNDT, B. SINSIN, G. GHEYSEN, and P. VAN DAMME
Patterns of Genetic and Morphometric Diversity in Baobab (Adansonia digitata) Populations Across Different Climatic Zones of Benin (West Africa)
Ann. Bot., May 1, 2006; 97(5): 819 - 830.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Bot.Home page
J. A. Lobo, M. Quesada, and K. E. Stoner
Effects of pollination by bats on the mating system of Ceiba pentandra (Bombacaceae) populations in two tropical life zones in Costa Rica
Am. J. Botany, February 1, 2005; 92(2): 370 - 376.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Submit a response
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via ISI Web of Science (5)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Silva-Montellano, A.
Right arrow Articles by Eguiarte, L. E.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Silva-Montellano, A.
Right arrow Articles by Eguiarte, L. E.
Agricola
Right arrow Articles by Silva-Montellano, A.
Right arrow Articles by Eguiarte, L. E.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS