Am. J. Bot. Cross-Journal Searching
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 (16)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Allphin, L.
Right arrow Articles by Harper, K. T.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Allphin, L.
Right arrow Articles by Harper, K. T.
Agricola
Right arrow Articles by Allphin, L.
Right arrow Articles by Harper, K. T.
(American Journal of Botany. 1998;85:1251-1261.)
© 1998 Botanical Society of America, Inc.


Genetic diversity and gene flow in the endangered dwarf bear poppy, Arctomecon humilis (Papaveraceae)1

Loreen Allphin5,a, Michael D. Windhamb and Kimball T. Harperc

a Department of Biology, University of Utah, Salt Lake City, Utah 84112; b Garrett Herbarium, Utah Museum of Natural History, University of Utah, Salt Lake City, Utah 84112; and c Department of Botany and Range Science, Brigham Young University, Provo, Utah 84602


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Arctomecon humilis is a critically endangered species endemic to the Moenkopi shale of Washington County, Utah. Recovery plans for the species would be improved by an understanding of genetic diversity and gene flow among its remaining populations. Ten variable isozyme loci were used to calculate genetic diversity statistics for study populations. Westerly populations possessed higher levels of genetic variability than other populations at the same isozyme loci. Three of the populations exhibited significant deviations from Hardy-Weinberg expectations. No correlation existed between genetic distance and geographic distance. Most of the genetic diversity was distributed among populations with little gene flow between populations, suggesting that observed genetic differences may arise from genetic drift. For the westerly populations, similar genotypes were observed in the seedling and old age classes, while intermediate age classes typically possessed an alternate set of genotypes at Pgi-2. Mean heterozygosity increased with age class across populations. Westerly populations of A. humilis shared more alleles with the nearest geographic population of A. californica than other populations. Since the westerly populations contained more genetic variability and more alleles in common with a near relative, they may be relictual. Other populations may contain less genetic diversity due to founder effects and/or genetic drift.

Key Words: Arctomecon • bear poppy • heterozygosity • gene flow • genetic diversity • Papaveraceae


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
An array of ecological factors and life history traits may affect the amount and distribution of genetic diversity in and among plant populations (Hamrick, Linhart, and Mitton, 1979; Loveless and Hamrick, 1984; Nevo, Beiles, and Ben-Shlomo, 1984; Hamrick and Godt, 1990). Genetic variability is known to increase fitness in populations of many plant and animal species (Hamrick, Linhart, and Mitton, 1979; Wills, 1981; Danzmann et al., 1986; Ledig, 1986). Increased levels of heterozygosity have been shown to increase longevity in long-lived perennial plant species (Schaal and Levin, 1976; Hamrick, Linhart, and Mitton, 1979). In addition, greater levels of genetic variation may buffer genotypes against environmental challenges. Conversely, the loss of genetic variability could render populations more vulnerable to extinction in cases of habitat perturbation, reproductive bottlenecks, etc. (Wright, 1933; O'Brien et al., 1985; Lacy, 1987; Simberloff, 1988; Barrett and Kohn, 1991).

According to Wright's (1943, 1946) isolation by distance models, mating is dependent on the distance between individuals and their ability to disperse propagules. Based on these models, inbreeding may occur if isolation increases and prevents gene flow between small patch populations. Habitat fragmentation and loss of critical habitat are becoming increasingly more common for many species, especially those in and around urban centers. Rare and/or endangered plant species often occur in small, disjunct populations with reduced genetic diversity due to increased habitat fragmentation and inbreeding (Stebbins, 1980; Simberloff, 1988; Barrett and Kohn, 1991; Godt, Johnson, and Hamrick, 1996; Sun, 1996).

Arctomecon humilis Coville (Papaveraceae), the dwarf bear poppy, is a short-lived perennial herb (average life span ~5 yr) restricted to gypsiferous shale substrates in Washington County, Utah (Atwood, 1977). The species occupies isolated sites within a 10-km radius of the city of St. George, Utah. During the last decade, urban development in the St. George area has eliminated or fragmented much of the critical habitat for the species. From 1970 to 1995, Washington County, Utah, grew by an average annual rate of 6.6%, which was by far the fastest growth rate of any county in the state of Utah (Christiansen and Robson, 1997). Washington County increased its population by 45% from 1990 to 1994 making it the fifth fastest growing county in the United States (Wright, 1996). By the year 2000, the city of St. George, Utah, has been estimated to grow to a population level of over 82 000 people, in addition to the estimated 2 million tourists that visit annually (Governor's Office of Planning and Budget, 1997).

The distribution of A. humilis is typically restricted to the white-gray (Shinob Kibe) or sometimes red gypsiferous members of the Moenkopi Formation, substrates that support only sparse vegetation (Nelson and Harper, 1991). These Moenkopi shale substrates are attractive to off-road vehicle users that not only destroy living plants, but also accelerate soil compaction and erosion. Much of this gypsum substrate is covered by mining claims. The 1872 General Mining Law required annual labor ($100 per claim) for the maintenance of these gypsum mining claims (MacDonnell, 1993). This work was often very disruptive and destructive. Many claim holders, for example, simply used a bulldozer to move soil around at a claim site, thus fulfilling the annual labor requirement. Today these activities are much reduced by mandatory fees (MacDonnell, 1993). Nevertheless, these labor requirements, in the past, have marred significant bear poppy habitat.

In 1979, the U.S. Fish and Wildlife Service listed A. humilis as an endangered species (U.S. Fish and Wildlife Service, 1979). In response to increasing threats, several hundred acres of poppy habitat were closed to off-road vehicle use by the Bureau of Land Management and the state of Utah (U.S. Bureau of Land Management, 1987). These closures include three populations (Red Bluff, Beehive Dome, and Warner Ridge), which are the largest extant populations of the species. Although these areas are posted, the closures have not completely curtailed use by off-road vehicles. In addition, housing development continues in the St. George area. These forces threaten the existence of the 11 extant populations of the species, especially the eight unprotected populations. Because of the narrow endemism and decline due to many human related factors, bear poppy has become one of the most critically endangered species in Utah.

Land managers faced with the difficult task of developing recovery plans for A. humilis have had to rely on information from habitat studies and demographic surveys (Nelson, 1989; Nelson and Harper, 1991). Recent analysis of DNA markers in these known populations using randomly amplified DNA procedures showed little variation among populations but offered no insight into genetic variation within populations, since the analysis was based on bulked samples from the individual populations (Van Buren and Harper, 1996). Management decisions regarding the future of these populations would be much improved by more detailed information concerning the genetics of the species in its various subpopulations (Falk and Holsinger, 1991).

A recent reproductive study of A. humilis reveals that the species has a mixed mating system (V. J. Tepedino, Utah State University, personal communication). The species requires an insect vector for successful fertilization of its ovules. The species is partially self-compatible. However, only approximately half of geitonogamously crossed flowers set fruit, while all of xenogamously crossed flowers set fruit (V. J. Tepedino, personal communication). The species is primarily bee pollinated (Tepedino and Griswold, 1989). The flowers are visited by an array of native, ground-nesting bees (often, a different suite of bees for each population) and Apis mellifera, the domestic honey bee. The native bees are often more efficient than the honey bees at pollination in A. humilis (visiting when stigmas are more receptive) and several of these native bee species are quite rare (V. J. Tepedino, personal communication; Tepedino and Griswold, 1989). Seed dispersal in this species is short-distance, primarily by ants attracted to the seeds by arils (Nelson, 1989).

Conservation and management suggestions for several rare plant and animal species have been made based on recent genetic analyses (Odasz and Savolainen, 1996; Sheely and Meagher, 1996; Cummings et al., 1997; Godt, Walker, and Hamrick, 1997; Ramirez and Froehlig, 1997; Sipes and Wolf, 1997). It is important that land mangers understand how genetic diversity is distributed among the remaining natural populations of A. humilis, if they are to preserve its remaining genetic diversity. This is especially critical considering the high levels of habitat loss and fragmentation due to housing development, gypsum mining claims, and off-road vehicle use, and considering the short-distance dispersal of pollen and seeds of this species. Short-distance dispersal of propagules, in association with increased isolation of populations due to habitat fragmentation, may result in inbreeding (Wright, 1943, 1946). In order to facilitate preservation of genetic diversity within the species, genetic variability and gene flow within and among the populations of A. humilis were evaluated in this study.

The research addressed in this manuscript explored the following questions. Does gene flow exist among these populations? How much genetic diversity exists within populations? How is genetic diversity distributed among populations and age classes within a given population? How does genetic variability in A. humilis compare with that of a near relative A. californica Torr.&Frem.?


    MATERIALS AND METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Study sites and sampling regime
For the allozyme analyses, leaf samples were taken from six populations of A. humilis from near St. George, Utah, and two populations of a near relative A. californica (Van Buren and Harper, 1996) found near Lake Mead, Arizona and Nevada (Fig. 1). The sampled populations of A. humilis encompass the total range of the species and are relatively isolated from one another. Most of the additional (unsampled) populations of A. humilis are located in heavily developed areas and are currently represented by small numbers of highly scattered individuals. The two populations of A. californica selected for this study are located on opposite sides of Lake Mead (Fig. 1).



View larger version (25K):
[in this window]
[in a new window]
 
Fig. 1. Localities of the six study populations of A. humilis in Washington County, Utah, and two study populations of A. californica from southeastern Nevada and northeastern Arizona. Circles represent western populations of A. humilis.

 
Leaves were collected from 24 to 30 individuals in each of the A. humilis populations (SE Boomer Hills, Warner Ridge, Shinob Kibe, Red Bluff, Beehive Dome, and White Dome; Table 1). Sampled individuals from each population were chosen from three different age classes (seedling, intermediate, and old) across the entire population. Age classes were assigned based upon phenotypic variability among tagged individuals of known age that have been followed demographically in the Red Bluff population for 10 yr: (1) plants that have not yet reached reproductive maturity represent the seedling age class, (2) individuals reproductively mature but showing no signs of senescence represent the intermediate age class, and (3) large individuals that show senescence represent the old age class. Prereproductive individuals are distinguished from those that skip a year of flowering by size and the absence of inflorescence stalk remnants from the previous year. Numbers of individuals sampled for each population and age class are given in Table 1. Voucher specimens were collected at each population and deposited in the Garret Herbarium, University of Utah.


View this table:
[in this window]
[in a new window]
 
Table 1. Estimated population size and numbers of individuals sampled by populations and age class for A. humilis and A. californica in this study. Population names are accompanied by abbreviations to be used throughout the remainder of the manuscript.

 
Electrophoresis
An electrophoretic survey was used to estimate genetic variability within and among selected populations of A. humilis. Leaf material of all sampled individuals was ground in a standard PVP-phosphate grinding buffer (Soltis et al., 1983). Ground material was absorbed into Whatman 3MM filter paper wicks and stored in an ultra-cold freezer (-70°C) until electrophoresis. Twenty enzymes were surveyed for variable loci in A. humilis and A. californica using a variety of gel and electrode conditions (Soltis et al., 1983). Allozyme markers for 10 enzyme loci were variable and provided consistent, interpretable results (Table 2).


View this table:
[in this window]
[in a new window]
 
Table 2. Enzymes used in isozyme analysis of A. humilis and A. californica.

 
Statistical analyses
Genetic data collected from the ten variable loci were used to calculate diversity statistics following Hamrick, Linhart, and Mitton (1979), Hedrick (1985), and Hamrick and Godt (1990) using Lynsprog, a statistical genetics program written by M. D. Loveless (College of Wooster, Wooster, Ohio). These statistics included: mean observed heterozygosity (Ho; direct estimate), polymorphic index (PI; mean proportion of polymorphic loci), Hardy-Weinberg expected heterozygosity (He; Nei and Li, 1979), mean number alleles per locus (A), mean number of effective alleles per locus (Ae; number of equally frequent alleles in an ideal population that would produce the same homozygosity as the actual population; Hartl, 1981), Nei's genetic identity (I), and Nei's genetic distance (D; Nei, 1972). A Mantel matrix randomization test (Mantel, 1967) was performed to determine whether there was a significant positive relationship between geographic distance and genetic distance (Rohlf, 1992).

Distribution of genetic variation among populations was estimated using Nei's (1973) genetic diversity statistics employing Lynsprog, a statistical genetics program written by M. D. Loveless (College of Wooster, Wooster, Ohio). These included: total genetic diversity (Ht), intrapopulation genetic diversity (Hs), diversity between populations (Dst), and the proportion of genetic variation distributed among populations (Gst). Gene flow (Nm, number of migrants per generation) among populations was calculated based upon Gst, using methods of Wright (1951).

The fixation index (F) was computed for each locus and each population using methods of Wright (1965). Chi-square analysis of Wright's fixation indices for study populations was used to test significance of deviations from Hardy-Weinberg expectations. Fixation indices and deviations from Hardy-Weinberg expectations were computed following Hamrick, Linhart, and Mitton (1979), Hedrick (1985), and Hamrick and Godt (1990) using Lynsprog, a statistical genetics program written by M. D. Loveless (College of Wooster, Wooster, Ohio).

Average linkage clustering using the unweighted pair-group method (UPGMA) was performed for study populations using Nei's genetic identities (Romesburg, 1984) in order to graphically portray genetic similarity among populations. Percentage similarity was computed among A. humilis populations for abiotic (soil) characteristics (Nelson and Harper, 1991) using Ruzicka's (1958) similarity index. A UPGMA phenogram was generated from these similarity indices (Romesburg, 1984). The UPGMA phenogram generated from abiotic characteristics of A. humilis populations (Nelson and Harper, 1991) was compared with the UPGMA phenogram constructed based on Nei's genetic identities, using only the populations considered in the abiotic UPGMA phenogram.

Percent similarity was computed among age-classes using genotype frequencies at the Pgi-2 locus and Ruzicka's (1958) similarity index. Mean Ho was calculated over all loci for each age class among all populations (Nei, 1973). Populations were bulked in this analysis because sample sizes by age class were small within populations (only 8–10 individuals). A one-way ANOVA and Tukey-Kramer multiple means comparisons were performed to test for significant differences among age-classes (SAS, 1994).


    RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Of the A. humilis populations, two of the westerly populations (Fig. 1), RB and SEBH exhibited the highest levels of genetic variability for all allozyme loci (Table 3). SEBH had the highest mean observed heterozygosity (Ho = 0.185) and the highest polymorphic index (PI = 0.172). This population also exhibited the highest mean number of alleles (A = 1.70) and the highest mean number of effective alleles per locus (Ae = 1.63). The population at WD exhibited the lowest genetic variability (Ho = 0.038; PI = 0.070). Red Bluff was the population with the greatest number of unique genotypes (3). All other populations had only one or no unique genotypes. The mean Ho for the species is 0.100.


View this table:
[in this window]
[in a new window]
 
Table 3. Levels of intrapopulations isozyme variability across populations of A. humilis and A. californica. Values are given as means across ten enzyme loci.

 
Arctomecon californica populations exhibited higher mean polymorphic indices (0.269), higher mean observed heterozygosities (0. 175), higher mean number of alleles per locus (A, 2.05), and a greater number of unique genotypes in a population (6–7) than A. humilis populations (Table 3). Allele frequencies used to generate these diversity data in A. humilis and A. californica are located in the Appendix of this manuscript.

Most of the total gene diversity (Ht = 0.339) in A. humilis was distributed among populations (Dst = 0.209). Gst (the proportion of genetic variation distributed among populations) was estimated as 0.617, thus 61.7% of the genetic variability was distributed among populations. Total gene flow across all populations was low (Nm = 0.155).

Fixation index (F) was determined for all loci in each population of A. humilis. Chi-square analysis demonstrated significant deviation from Hardy-Weinberg expectations for some of the populations and loci (Table 4). SE Boomer Hills, the most genetically variable of the study populations, showed significant deviation from Hardy-Weinberg at two enzyme loci. Warner Ridge and WD also exhibited significant deviation (P <= 0.05) from Hardy-Weinberg expectations at Tpi-3 (Table 4). Since six populations were sampled and ten loci, there were a total of 60 determinations of F deviations from Hardy-Weinberg equilibrium. For A. humilis, four of the 60 determinations (~7%) deviated significantly (P <= 0.05) from Hardy-Weinberg. This is only slightly more than you would randomly expect to be significant by chance alone (at the 5 % level). However, only two populations of A. californica were sampled for this study. Of the 20 possible determinations of F deviations from Hardy-Weinberg expectations, five showed significant deviation (Pgm-1 and Tpi-1 for both populations and Nadh-dp for the OVER population; Table 4). This is much greater than might be expected by chance alone. In addition, the He values for these populations are much higher than the Ho values (Table 3), thus suggesting that the two A. californica populations sampled are significantly deviating from Hardy-Weinberg expectations.


View this table:
[in this window]
[in a new window]
 
Table 4. Chi-square values generated using Lynsprog, following Hedrick (1985), for selected enzyme loci. Populations exhibiting significant deviations of F (fixation index) from zero are followed by a symbol.

 
Nei's genetic distances (D) among populations of A. humilis ranged from 0.008 to 0.212 and geographic distances among all A. humilis populations ranged from 5.1 to 21.3 km (Table 5). A Mantel randomization test showed no significant positive relationship between geographic distance and genetic distance (P = 0.4057; r = -0.0686). Genetic identity values among study populations were obtained by subtracting Nei's genetic distances from 1 (Nei, 1972). Population clustering based upon Nei's genetic identities among populations is shown in Fig. 2. The Beehive Dome and WD populations were the most genetically similar (I = 0.992). Red Bluff and SEBH were also similar with respect to genetic identity derived from allozyme markers (I = 0.981; Fig. 2). The Shinob Kibe population was the most genetically unique, exhibited the lowest genetic identity values with other sampled populations (I = 0.865), of the six A. humilis populations used in this study. This was especially evident at the Pgi-2 locus for which SK exhibited a high allele frequency for an allele (Pgi-2B) that was rarely found in other populations (Fig. 3; Appendix).


View this table:
[in this window]
[in a new window]
 
Table 5. Geographic distances and genetic distances between sampled populations of A. humilis and A. californica. Geographic distances (km) are given above the diagonal and Nei's genetic distances (D ) are given below the diagonal.

 


View larger version (12K):
[in this window]
[in a new window]
 
Fig. 2. Phenogram showing relationships of sampled populations of A. humilis and A. californica based upon Nei's (1972) genetic identities among populations calculated from allele frequencies.

 


View larger version (42K):
[in this window]
[in a new window]
 
Fig. 3. Gel photograph illustrating allozyme variability at the Pgi-2 locus for all sampled populations of Arctomecon.

 
The genetic groupings of the sampled populations were similar to those generated from abiotic (soil) characteristics for the same population (Nelson and Harper, 1991; Fig. 4A). The two phenograms differed, however, in their placement of WR and SK populations. Red Bluff (RB) was genetically more similar to SEBH than WR. White Dome (WD) was more similar genetically to RB, SEBH, and WR than the SK, although they share abiotic similarities (Fig. 4B).



View larger version (14K):
[in this window]
[in a new window]
 
Fig. 4. (A) Phenogram showing percentage age similarity among populations of A. humilis based upon abiotic (site) characteristics, modified and redrawn from Nelson (1989). (B) Phenogram of population clusterings of same populations based upon genetic identities (Nei, 1972).

 
When the data were considered collectively for all populations of A. humilis and A. californica, the two sampled populations of A. californica had a low genetic identity (I = 0.867) in comparison to relatively higher genetic identities exhibited among the A. humilis populations (I = 0.927). Arctomecon humilis is genetically very different from A. californica. The mean genetic identity between A. humilis and A. californica populations was low (I = 0.311). The western populations of A. humilis (RB and SEBH) had a higher average genetic identity with A. californica (I = 0.357) than did the eastern populations (WD, BD, WR, and SK; I = 0.288). The western populations shared additional alleles in common with A. californica (Pgi-2D and Pgm-1B) that were not present in the eastern populations (Fig. 3; Appendix).

The genetically diverse western populations of A. humilis (RB and SEBH) also exhibited patterns of temporal (age class) variability evident only at the Pgi-2 locus. Similar genotypes for this locus were observed in the seedling and old age classes, while intermediate age classes typically possessed an alternate set of genotypes at Pgi-2 (Table 6). However, while the temporal variability was evident at both populations, RB and SEBH had different distributions of genotype frequencies among the age classes. Percentage age similarity (Ruzicka, 1958) was calculated among age classes for both the RB and SEBH populations based on genotype frequencies at Pgi-2. Seedlings and old individuals are 56.2% similar at Pgi-2 for RB and SEBH. Seedlings and intermediate-aged individuals are only 23.2% similar and intermediates and old individuals are 24.6% similar. Statistical analyses were not performed due to the small sample sizes by age-class for each population and the inability to bulk the populations due to the different distribution of genotype frequencies within the two populations (Table 6).


View this table:
[in this window]
[in a new window]
 
Table 6. Genotypic frequencies at the Pgi-2 locus across three age classes from the RB and SEBH populations of A. humilis.

 
Across all selected populations of A. humilis, mean observed heterozygosity increased with age (Fig. 5). The old age class had the highest mean observed heterozygosity of 0.238, while the seedlings exhibited a significantly lower mean observed heterozygosity of 0.184 (F ratio = 37.22; P = 0.001). All populations were bulked in this analysis due to small sample size for each of the age classes per population.



View larger version (121K):
[in this window]
[in a new window]
 
Fig. 5. Mean observed heterozygosity among three age classes for all sampled populations of A. humilis at ten loci (all populations are pooled due to the small sample size by age in each population). Seedlings and old age classes are significantly different at P = 0.001.

 

    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
The results presented in this manuscript have many implications for managers of the dwarf bear poppy. The two westernmost populations (RB and SEBH) exhibit higher levels of genetic variability than the other sampled populations. Other populations show low levels of genetic variability for the same allozyme loci, possibly due to founder effects or genetic drift. The population exhibiting the highest mean observed heterozygosity and polymorhpic index, SEBH, is not one of the three populations closed to off-road vehicle use. One-half of this population is located on land managed by the Bureau of Land Management (BLM), while the other is on private land. Cattle currently graze both the private and BLM lands containing this population.

The most genetically unique of the sampled populations is Shinob Kibe. This population lies on privately owned land, which is being heavily developed. The Nature Conservancy is now attempting to purchase ~28 ha of land adjacent to the Shinob Kibe Butte. If purchased, the land would be fenced to exclude off-road vehicles and casual foot traffic. We previously mentioned a recent study performed on the genus Arctomecon using RAPD (randomly amplified polymorphic DNA) markers (Van Buren and Harper, 1996). The RAPD study suggested that the Shinob Kibe (SK), the most northern population, had a genetic make-up somewhat different from other populations because it had the fewest number of RAPD markers shared with other populations. These data are consistent with our findings based on genetic identity values among populations. Shinob Kibe had the lowest mean genetic identity value with other sampled populations (Fig. 2).

Most of the total genetic diversity (Ht) of the species, using data from this study, is distributed between populations and gene flow between populations is low. Populations of A. humilis are becoming progressively more fragmented and isolated from one another due to housing developments and disturbance from off-road vehicles. However, Van Buren and Harper (1996) stated that A. humilis showed little RAPD variation among populations. In this manuscript, we find that 61.7% of all the variability in A. humilis is distributed among populations. Van Buren and Harper (1996) bulked leaf material from all sampled individuals in a population into a single sample for each population. This type of bulking prevents any interpretation on variability within populations. This type of sampling also prevents an accurate assessment of the distribution of genetic variability among populations. Markers are typically noted for a particular population without any consideration of their frequency within the population.

Distribution of genetic diversity may be affected by spatial and historical factors (Kimura and Maruyama, 1971; Slatkin, 1987; Sheely and Meagher, 1996) and/or ecological and life history traits (Hamrick et al., 1979; Karron, 1987; Karron et al., 1988; Hamrick and Godt, 1990). The historical range and past population sizes in A. humilis are relatively unknown. However, one of the authors of this manuscript (K. T. Harper) has observed a decline in the number of individuals in several of the populations over the past 20 yr.

Arctomecon humilis is confined to the Shinob Kibe member of the Moenkopi formation. This geologic formation is often discontinuous across the landscape surface. Some of the sampled populations, therefore, may have been isolated for a long period of time, because the substrate was laid down during the Holocene (Hintze, 1973). The SK population is located on the edge of the Shinob Kibe Butte that is more geologically separated from the other populations, which might explain its genetic uniqueness.

Some of the populations (RB, SEBH, and WD), on the other hand, are growing on different parts of a, at one time, continuous member (Shinob Kibe) of the Moenkopi formation. Little development has occurred between the RB and SEBH populations. It is, therefore, likely that gene flow events may occasionally occur between these two populations. However, much urban development has occurred between the RB and WD populations, including the construction of a major interstate (Interstate 15) and a suburban development (Bloomington). This urban development might explain the genetic differentiation between the WD and the RB and SEBH populations. Off-road vehicle users have also heavily hit the WD population. The authors have seen a significant reduction in the number of individuals in this population over the past 10 yr. The WD population was also the least genetically variable of the A. humilis populations.

Since genetic distance could not be correlated with geographic distance in this study, historic and geologic distribution of populations may be responsible for the observed distribution of genetic diversity among populations, rather than a disruption of gene flow due to habitat fragmentation. However, additional study is necessary. We feel that this manuscript provides baseline genetic information for future studies. In the future, allozyme variation should be reassessed in these populations to look at the long-term effects of human induced habitat fragmentation and habitat degradation on genetic diversity and genetic structure in this species.

Some populations of A. humilis, especially SEBH, appear to be inbreeding (F values varied significantly from zero). This population was the smallest of the study populations (<100 individuals, Table 1). Inbreeding may result in reduced fitness for some of these populations (Schemske and Lande, 1985; Simberloff, 1988; Levin, 1991). However, more compelling are the data on F deviations from zero in A. californica. Because only two populations of this more common species were sampled in this study, we will refrain from discussing possible reasons for the H-W deviations. We note that a graduate student (L. Hitchcock, Utah State University) is assessing allozyme diversity in A. californica. It will be interesting to compare the results of the two studies upon completion.

Mean observed heterozygosity is significantly lower in A. humilis (Ho = 0.100) than other outcrossing plant species (Ho = 0.167) and dicots (Ho = 0.136), but only slightly lower than other short-lived herbaceous perennials (Hamrick and Godt, 1990; Table 7). Genetic diversity in A. humilis is consistent with the level of variation observed in other narrow endemic species (Ho = 0.096; Hamrick and Godt, 1990). Arctomecon humilis populations exhibit less genetic variability on average than the sampled populations of A. californica. Both species are endemic to gypsic soils of Mojave Desert regions. However, A. californica is geographically more widespread and occupies an ecologically broader range of habitats than A. humilis (Nelson and Welsh, 1993). These results are consistent with examples in the literature of geographically restricted species having less genetic variation than more widespread congeners (Karron, 1987; Karron et al., 1988; Hamrick et al., 1991; Hamrick and Godt, 1996; Godt, Walker, and Hamrick, 1997). Arctomecon californica exhibited significantly higher number of effective alleles per locus than other comparative species groups (Hamrick and Godt, 1990; Table 7). The species' mean observed heterozygosity (0.175) was higher than other dicots, short-lived, herbaceous perennials, and outcrossing (requiring an animal vector) and endemic species (Hamrick and Godt, 1990; Table 7).


View this table:
[in this window]
[in a new window]
 
Table 7. Level of allozyme variation at species level for species in different categories. The data were extracted from Hamrick and Godt (1990). These data are compared with allozyme variation at the species level for both A. humilis and A. californica. Standard error values are given in parentheses below the means. Statistical differences among means for species or groups were determined using a multiple means comparison (Snedecor and Cochran, 1967). means followed by the same letter in a column are not significantly different at P < 0.05.

 
The genetic identities among the six study populations of A. humilis (I = 0.937) were within the range expected of conspecific populations (Crawford, 1983). However, genetic identity values between the two taxa were extremely low (I = 0.311), lower than expected for closely related congeners (Crawford, 1983). More alleles were shared in common between the two westerly populations of A. humilis and the nearest geographic population of A. californica in Nevada than other populations of A. humilis. However, it is uncertain whether this similarity represents a historic gene flow event or alleles shared through common ancestry. Since the westerly populations of A. humilis contain more genetic variability and more alleles in common with a near relative, they may be relictual. Other populations may contain less genetic diversity due to founder effects and/or genetic drift. Because the RAPD study (Van Buren and Harper, 1996) included both of the taxa that were used in this study, comparisons can be made between the two studies. For example, RAPD data support A. humilis and A. californica as distinct species. This is consistent with our findings based on genetic identity. However, further study is suggested for all of the Arctomecon species using sensitive molecular techniques.

Allozyme markers have been considered as "neutral" genes reflecting evolutionary processes affecting the entire genome (Avise, 1994). Mitton (1994), however, postulates that environment may occasionally select for allozyme markers. The phenogram resulting from cluster analysis of the populations based upon Nei's genetic identities is similar to that based on soil analyses for the same populations (Nelson and Harper, 1991). These results suggest that abiotic conditions may have selected for the allozyme markers observed in A. humilis. Before this hypothesis can be accepted, however, genetic variability within and among populations must be assessed using a molecular technique with known neutrality of markers.

The patterns of genetic variability among age classes in the Red Bluff and SE Boomer Hill populations were only observed at the Pgi-2 locus, which is a rapidly evolving and highly variable enzyme locus (Gillespie and Kojima, 1968). The observed genetic differences among the three age classes indicate that age structure in existing populations is important and should be considered in preservation of genetic diversity within A. humilis populations.

Seeds of the Papaveraceae are dispersed before they mature and require a prolonged period of after-ripening followed by a specific weather regime (cool, wet springs) to germinate (Nelson, 1989). This may require years of residence in the seed bank before germination. The time necessary for completion of development and germination is unknown. Therefore, the seed bank undoubtedly has a major impact on a population's gene pool. Seed banks are believed to buffer populations against changes in genetic composition and increase the effective population size (Epling, Harlan, and Ball, 1959). They may also serve as a source of novel genetic variation because of mutations in seeds as they age (Levin, 1990). Soil erosion and compaction, however, could eliminate a viable seed bank.

Age classes of A. humilis plants differ with respect to mean heterozygosity. Heterozygosities and polymorphic indices have been demonstrated to increase with age in some trees and long-lived perennials (Schaal and Levin, 1976; Hamrick, 1979). It appears, however, that even in a short-lived perennial (average longevity <5 yr) such as A. humilis, heterozygosity may increase survivorship and thus longevity. Such a conclusion is supported by recent work with two herbaceous perennials (Cabin, 1996; Hossaert-McKey et al., 1996).

Recovery plans for endangered species, especially A. humilis, should not only consider the genetic variability among populations, but variability within populations. This research provides exceptions to common generalizations. First, it is risky to generalize that the largest populations of any species are the most variable genetically. Several of the large populations of A. humilis possessed very little genetic variability. Second, protection for an endangered species should include not only all ages within a given population, but the seed bank as well. In order to preserve genetic variability in species where populations are declining as with A. humilis, managers should attempt to preserve individuals representative of all age classes. In situations where areas are reseeded, seed should be collected from all age classes in a given population, as well as from the seed bank.

We recognize, however, the limitation of our small sample sizes per population (24–30 individuals). Interpretation of within-population and age-class variation is difficult when only 30 individuals out of several hundred are sampled. Therefore, future genetic studies in A. humilis should also include larger sample sizes for a better understanding of the distribution of variation among age classes within populations in A. humilis.

Because the historical range and past population sizes in A. humilis are relatively unknown, interpretation of the data presented in this manuscript is difficult. However, it is evident that populations are becoming fragmented and lost. Off-road vehicles and development still pose a serious threat to the existence of the dwarf bear poppy. If these current human activities continue, it is likely that gene flow among populations will continue to be interrupted and genetic differentiation and drift among populations will therefore increase. We also feel that the data in this manuscript illustrate the importance of maintaining the age structure within populations. Since the seed bank is vital to the species, it is critical that land managers enforce restrictions to off-road vehicle use. Additional dwarf bear poppy populations will likely go extinct. However, land managers and groups, like the Nature Conservancy, can use these genetic data to ensure that the genetic diversity within A. humilis is preserved.


    FOOTNOTES
 
1 The authors thank Dr. George Yatskievych, Missouri Botanical Garden, for his assistance with the statistical genetic data analyses; Wade Woolstenhulme, Monica Brady, and Molly Windham for their assistance with field tissue collection; and Heidi Christy and Elizabeth Bodmer for their assistance with the enzyme electrophoresis. This research was partly supported by grants from the U.S. Fish and Wildlife Service and the Bureau of Land Management.


View this table:
[in this window]
[in a new window]
 
Appendix. Allele frequencies at ten polymorphic loci for sampled A. humilis and A. californica populations.

 
Back 5 Author for correspondence, current address: Department of Botany and Range Science, Brigham Young University, Provo, UT 84602 (email:LAWoolst{at}bioag.byu.edu). Back


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Atwood, N. D.1977.The dwarf bear-poppy. Mentzelia 3: 6–7.

Avise, J. C.1994.Molecular markers, natural history and evolution. Chapman and Hall, New York, NY.

Barrett, S. C., and J. R. Kohn.1991.Genetic and evolutionary consequences of small population size in plants: implications for conservation. In D. A. Falk and K. E. Holsinger [eds.], Genetics and conservation of rare plants, 3–30. Oxford University Press, New York, NY.

Cabin, R. J.1996.Genetic comparisons of seed bank and seedling populations of a perennial mustard, Lesquerella fendleri. Evolution 50:1830–1841.

Christiansen, M. E., and J. Robson.1997.Regional/national comparisons, population urbanization trends in Utah. In M. L. Wood [ed.], 1997 Economic report to the governor. Governor's Office of Planning and Budget, Salt Lake City, UT.

Crawford, D. J.1983.Phylogenetic and systematic inferences from electrophoretic studies. In S. D. Tanksley and J. J. Orton [eds.], Isozymes in plant breeding, part A, 257–287. Elsevier, Amsterdam.

Cummings, S. A., E. L. Brannon, K. J. Adams, and G. H. Thorgaard.1997.Genetic analyses to establish captive breeding priorities for endangered Snake River sockeye salmon. Conservation Biology 11: 662–669. [CrossRef][ISI]

Danzmann, R. G., M. M. Ferguson, F. W. Allendorf, and K. L. Knudsen.1986.Heterozygosity and developmental rate in a strain of rainbow trout. Evolution 40: 86–93. [CrossRef][ISI]

Epling, C., L. Harlan, and F. M. Ball.1959.The breeding group and seed storage: a study in population dynamics. Evolution 14: 238–255. [CrossRef][ISI]

Falk, D. A., and K. Holsinger.1991.Genetics and conservation of rare plants. Oxford University Press, New York, NY.

Gillespie, J. H., and K. Kojima.1968.The degree of polymorphism in enzymes involved in energy production compared to that in nonspecific enzymes in two Drosophila ananassae populations. Proceedings of the National Academy of Sciences, USA 61: 582–585. [Free Full Text]

Godt, M. J. W., B. R. Johnson, and J. L. Hamrick.1996.Genetic diversity and population size in four rare southern Appalachian plant species. Conservation Biology 10: 796–805. [CrossRef][ISI]

———, J. Walker, and J. L. Hamrick.1997.Genetic diversity in the endangered lily Harperocallis flava and a close relative, Tofieldia racemosa. Conservation Biology 11: 361–366.

Governor's Office of Planning and Budget.1997.Utah's long-term projections. In M. L. Wood [ed.], 1997 Economic report to the governor. Utah Office of Planning and Budget, Salt Lake City, UT.

Hamrick, J. L.1979.Genetic variation and longevity. In O. Solbrig, S. Jain, G. Johnson, and P. Raven [eds.], Topics in plant reproductive biology, 84–113. Columbia University Press, New York, NY.

———, and M. J. W. Godt.1990.Allozyme diversity in plant species. In A. H. D. Brown, M. T. Clegg, A. L. Kahler, and B. S. Weir [eds.], Plant population genetics, breeding, and genetic resources, 43–63. Sinauer, Sunderland, MA.

———, ———, D. A. Murawski, and M. D. Loveless.1991. Correlations between species traits and allozyme diversity: implications for conservation biology. In D. A. Falk and K. E. Holsinger [eds.], Genetics and conservation of rare plants, 75–86. Oxford University Press, New York, NY.

———, and ———.1996.Conservation genetics of endemic plant species. In J. C. Avise and J. L. Hamrick [eds.], Conservation genetics: case histories from nature, 281–304. Chapman and Hall, New York, NY.

———, Y. B. Linhart, and J. B. Mitton.1979.Relationships between life history characteristics and electrophoretically detectable genetic variation in plants. Annual Review of Ecology and Systematics 10: 173–200.

Hartl, D. L.1981.A primer of population genetics. Sinauer, Sunderland, MA.

Hedrick, P. W.1985.Genetics of populations. Jones and Bartlett, Boston, MA.

Hintze, L.1973.Geologic history of Utah. Brigham Young University Geologic Studies 20(3): 164–181.

Hossaert-McKey, M. M., M. Valero, D. Magda, M. Jarry, J. Cuguen, and P. Vernet.1996.The evolving genetic history of a population of Lathyrus sylvestrus: evidence from temporal and spatial genetic structure. Evolution 50: 1808–1812. [CrossRef][ISI]

Karron, J. D.1987.A comparison of levels of genetic polymorphism and self-compatibility in geographically restricted and widespread plant congeners. Evolutionary Ecology 1: 47–58.

———, Y. B. Linhart, C. A. Chaulk, and C. A. Robertson.1988.Genetic structure of populations of geographically restricted and widespread species of Astragalus (Fabaceae). American Journal of Botany 75: 1114–1119. [CrossRef][ISI]

Kimura, M., and T. Muruyama.1971.Pattern of neutral polymorphism in a geographically structured population. Genetical Research, Cambridge 18: 125–130.

Lacy, R. C.1987.Loss of genetic diversity from managed populations: interacting effects of drift, mutation, immigration, selection, and population subdivision. Conservation Biology 1: 143–158.

Ledig, F. T.1986.Heterozygosity, heterosis, and fitness in outbreeding plants. In M. E. Soulé [ed.], Conservation biology: the science of scarcity and diversity, 77–104. Sinauer, Sunderland, MA.

Levin, D. A.1990.The seed bank as a source of genetic novelty in plants. American Naturalist 135: 563–572. [CrossRef][ISI]

———.1991.The effect of inbreeding on survivorship in Phlox. Evolution 45: 1047–1049. [CrossRef]

Loveless, M. D., and J. L. Hamrick.1984.Ecological determinants of genetic structure in plant populations. Annual Review of Ecology and Systematics 15: 65–95.

MacDonnell, L. J.1993.Mineral law in the United States: a study in legal change. In L. J. MacDonnell and S. F Bates [eds.], Natural resources policy and law: trends and directions, 66–93. Island Press, Washington, DC.

Mantel, N.1967.The detection of disease clustering and a generalized regression approach. Cancer Research 27: 209–220. [ISI][Medline]

Mitton, J. G.1994.Molecular approach to population biology. Annual Review of Ecology and Systematics 25: 45–69.

Nei, M.1972.Genetic distance between populations. American Naturalist 106: 283–292. [CrossRef][ISI]

———.1973.Analysis of gene diversity in subdivided populations. Proceedings of the National Academy of Sciences, USA 70: 3321–3323. [Abstract/Free Full Text]

———, and W. H. Li.1979.Mathematical models for studying genetic variation in terms of restriction endonucleases. Proceedings of the National Academy of Sciences, USA 76: 191–193.

Nelson, D. R.1989.Site characteristics and habitat requirements of the endangered dwarf bear-claw poppy (Arctomecon humilis Coville, Papaveraceae). Master's thesis, Brigham Young University, Provo, UT.

———, and K. T. Harper.1991.Site characteristics and habitat requirements of the endangered dwarf bear-claw poppy (Arctomecon humilis Coville, Papaveraceae). Great Basin Naturalist 51: 167–175.

———, and S. L. Welsh.1993.Taxonomic revision of Arctomecon Torr.&Frem. Rhodora 95: 197–215. [ISI]

Nevo, E., A. Beiles, and R. Ben-Shlomo.1984.The evolutionary significance of genetic diversity: ecological, demographic, and life-history correlates. In G. S. Mani [ed.], Evolutionary dynamics of genetic diversity, 13–213. Springer-Verlag, Berlin.

O'Brien, S. J., M. E. Roelke, M. A. Newman, C. A. Winkler, D. Meltzer, L. Colly, J. F. Eveman, M. Bush, and D. E. Wildt.1985.Genetic basis for vulnerability in the cheetah. Science 227: 1428–1434. [Abstract/Free Full Text]

Odasz, A. M., and O. Savolainen.1996.Genetic variation in populations of the arctic perennial Pedicularis dasyantha (Scrophulariaceae), on Svalbard, Norway. American Journal of Botany 83: 1379–1385. [CrossRef][ISI]

Ramirez, M. G., and J. L. Froehlig.1997.Minimal genetic variation in a coastal dune arthropod: the trapdoor spider Aptostichus simus (Cyrtaucheniidae). Conservation Biology 11: 256–259.[CrossRef][ISI]

Rohlf, F. J.1992.NTSYS-pc: Numerical taxonomy and multivariate analysis system, version 1.70. Applied Biostatistics, Setauket, NY.

Romesburg, H. C.1984.Cluster analysis for researchers. Lifetime Learning Publications, Belmont, CA.

Ruzicka, M.1958.Awndung matematisch-statishticher methoden in der geobotanik (syshetische bearbutung von autnahme). Biologia, Bratisl 13: 647–661.

SAS.1994.JMP Statistics for the Apple Macintosh, version 3. SAS Institute, Cary, NC.

Schaal, B. A., and D. A. Levin.1976.The demographic genetics of Liatris cylindracea MICHX. (Compositae). American Naturalist 110: 191–205. [CrossRef][ISI]

Schemske, D. W., and R. Lande.1985.The evolution of self-fertilization and inbreeding depression in plants. II. Empirical observations. Evolution 39: 41–52. [CrossRef][ISI]

Sheely, D. L., and T. R. Meagher.1996.Genetic diversity in Micronesian island populations of the tropical tree Campnosperma brevipetiolata (Anacardiaceae). American Journal of Botany 83: 1571–1579. [CrossRef][ISI]

Simberloff, D. W.1988.The contribution of population and community biology to conservation science. Annual Review of Ecology and Systematics 19: 473–511.

Sipes, S. D., and P. G. Wolf.1997.Clonal structure and patterns of allozyme diversity in the rare endemic Cycladenia humilis var. jonesii (Apocynaceae). American Journal of Botany 84: 401–409. [Abstract]

Slatkin, M.1987.Gene flow and the geographic structure of natural populations. Science 236: 787–792. [Abstract/Free Full Text]

Snedecor, G. W., and W. G. Cochran.1967.Statistical methods, 6th ed. Iowa State University Press, Ames, IA.

Soltis, D. E., D. H. Haufler, D. C. Darrow, and J. 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]

Stebbins, G. L.1980.Rarity of plant species: a synthetic viewpoint. Rhodora 82: 77–86. [ISI]

Sun, M.1996.Effects of population size, mating system, and evolutionary origin on genetic diversity in Spiranthes sinenesis and S. hongkongensis. Conservation Biology 10: 785–795.

Tepedino, V. J., and T. L. Griswold.1989.Pollination biology of endangered plants. In Annual report, grasshopper integrated pest management project, 74–81. USDA, Animal and Plant Health Inspection Service, Plant Protection and Quarantine, Boise, ID.

U.S. Bureau of Land Management.1987.Off-road vehicle closure, Utah. Federal Register 52(49): 7941.

U.S. Fish and Wildlife Service.1979.Determination that Arctomecon humilis is an endangered species. Federal Register 44(26): 64250–64252.

Van Buren, R., and K. T. Harper.1996.Genetic variation among populations of Arctomecon. Southwestern rare and endangered plants: proceedings of the second conference. USDA Forest Service General Technical Report RM-GTR-283: 77–85.

Wills, C.1981.Genetic variability. Clarendon Press, Oxford.

Wright, L.1996.St. George: blooming flower of spreading weed. Salt Lake Tribune, April 7, page B-1.

Wright, S.1933.Inbreeding and homozygosis. Proceedings of National Academy of Sciences, USA 19: 411–433. [Free Full Text]

———.1943.Isolation-by-distance. Genetics 28: 114–138. [Free Full Text]

———.1946.Isolation by distance under diverse systems of mating. Genetics 31: 39–59. [Free Full Text]

———.1951.The genetic structure of populations. Annals of Eugenics 15: 323–354. [ISI]

———.1965.The interpretation of population structure by F-statistics with special regard to systems of mating. Evolution 19: 395–420. [CrossRef][ISI]




This article has been cited by other articles:


Home page
Am. J. Bot.Home page
M. C. Neel and N. C. Ellstrand
Patterns of allozyme diversity in the threatened plant Erigeron parishii (Asteraceae)
Am. J. Botany, May 1, 2001; 88(5): 810 - 818.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Bot.Home page
M. A. Gitzendanner and P. S. Soltis
Patterns of genetic variation in rare and widespread plant congeners
Am. J. Botany, June 1, 2000; 87(6): 783 - 792.
[Abstract] [Full Text]


Home page
Am. J. Bot.Home page
K. Schmidt and K. Jensen
Genetic structure and AFLP variation of remnant populations in the rare plant Pedicularis palustris (Scrophulariaceae) and its relation to population size and reproductive components
Am. J. Botany, May 1, 2000; 87(5): 678 - 689.
[Abstract] [Full Text]


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 (16)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Allphin, L.
Right arrow Articles by Harper, K. T.
Right arrow Search for Related Content
PubMed