Am. J. Bot.
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 ISI Web of Science (20)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Mayes, S. G.
Right arrow Articles by Werth, C. R.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Mayes, S. G.
Right arrow Articles by Werth, C. R.
Agricola
Right arrow Articles by Mayes, S. G.
Right arrow Articles by Werth, C. R.
(American Journal of Botany. 1998;85:1609-1617.)
© 1998 Botanical Society of America, Inc.

Clonal population structure and genetic variation in sand-shinnery oak,Quercus havardii (Fagaceae)1

Steven G. Mayesa, Mark A. McGinleya and Charles R. Werth2,a

a Ecology Program, Department of Biological Sciences, Texas Tech University, Lubbock, Texas 79409–3131


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
We investigated clonal population structure and genetic variation in Quercus havardii (sand-shinnery oak), a deciduous rhizomatous shrub that dominates vegetation by forming uninterrupted expanses of ground cover over sandy deposits on the plains of western Texas, western Oklahoma, and eastern New Mexico. Isozyme electrophoresis (15 loci coding 11 enzymes) was used to recognize and map clones arrayed in a 2000-m transect (50-m sample intervals) and a 200 x 190 m grid (10-m sample intervals). Ninety-four clones were discovered, 38 in the transect and 56 in the grid, resulting in an estimated density of ~15 clones per hectare. Clones varied greatly in size (~100–7000 m), shape, and degree of fragmentation. The larger clones possessed massive interiors free of intergrowth by other clones, while the smaller clones varied in degree of intergrowth. The population maintained substantial levels of genetic variation (P = 60%, A = 2.5, Hexp = 0.289) comparable to values obtained for other Quercus spp. and for other long-lived perennials. The population was outcrossing as evidenced by conformance of most loci to Hardy-Weinberg expected genotype proportions, although exceptions indicated a limited degree of population substructuring. These data indicate that despite apparent reproduction primarily through vegetative means, Q. havardii possesses conventional attributes of a sexual population.

Key Words: allozyme • clone • Fagaceae • genetic variation • isozyme • population structure • Quercus


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Among the most unusual of biotic communities in North America is the sand-shinnery oak community occurring on sandy deposits in western Texas, eastern New Mexico, and southwestern Oklahoma (Fig. 1). The overwhelming dominant in this community is the plains-endemic Quercus havardii Rydb. (sand-shinnery oak), a clonal deciduous shrub of low stature, with sparsely branched aerial shoots rarely exceeding 1 m in height. It forms large patches through horizontal growth effected by subterranean rhizomes, which may achieve a diameter of >9 cm (Muller, 1951; Wiedman, 1960; Scifres, 1972; Pettit, 1979; Dhillion et al., 1994). In the eastern portion of the species' range these patches remain more or less discrete, but in the more western sand-shinnery oak community, the patches coalesce to form a continuous ground cover with stem densities ranging from 30 to 70 stems/m (Dhillion et al., 1994; Zhang, 1996). The resulting aspect of the community is a vast homogeneous "sea" of uniformly spaced Q. havardii stems similar in height and other attributes. This vegetation type is estimated to have originally covered an area of 2.3 million hectares, but present abundance has been substantially reduced through eradication, as Q. havardii is considered a pest by ranchers (Dhillion et al., 1994).



View larger version (41K):
[in this window]
[in a new window]
 
Fig. 1. The distribution of sand-shinnery oak and the location of our field site in Yoakum County, Texas. The sand-shinnery oak community per se is predominantly in the western continuous portions of the species' range.

 
Demographic features and clonal biology of Q. havardii have been sparingly studied and remain poorly understood. Individual ramets generally live <15 yr (Wiedman, 1960; Pettit, 1979), but genets may be of considerable age and were previously estimated to expand to sizes up to ~15 m (Muller, 1951; Wiedman, 1960). Herein we address clonal population structure and genetic variation in Q. havardii. The goal of this study was to use molecular markers (allozymes) to discriminate clones and thereby facilitate perception of a population comprising discrete genets of sexual origin. From this new perspective, features of the genets such as size, shape, and intergrowth could be evaluated, and levels and patterns of genetic variation could be addressed in a manner parallel to that available for nonclonal species.


    MATERIALS AND METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Study site and sampling strategy
This study was conducted at a permanent research site near Plains in Yoakum County, Texas (Fig. 1; Dhillion et al., 1994). The habitat consists of sand dunes dominated by Quercus havardii, which comprises >90% of the ground cover, interrupted in spots by disturbed sites and other plant species. Associated plant species included Artemisia filifolia, Yucca campestris, Schizachyrium scoparium, Andropogon hallii, Aristida purpurea, Eriogonum annuum, and Euphorbia fendleri (Dhillion et al., 1994; names of taxa based on Correll and Johnston, 1979).

As there was previously little more than anecdotal information on clone size for Quercus havardii (Muller, 1951), the most appropriate distance between sampling points was undetermined as this study commenced. To establish a first approximation of clone size, 40 samples were taken at 50-m intervals along a transect and electrophoresed (see below) to determine whether clones extended >50 m. As most samples in this transect possessed distinct genotypes and were therefore different genets (see Results), a finer sampling interval of 10 m was chosen. A grid measuring 200 x 190 m was established at the study site from which samples were obtained at 10-m intervals.

Electrophoresis of enzymes
Leaves were obtained for electrophoretic analysis from April through July 1993 and placed in sealable bags in an insulated container, and kept refrigerated until processing. Leaf tissue was homogenized within 3 d after harvest (longer delays resulted in significant loss of enzyme activity) in pH 7.5 "microbuffer" of Werth (1985) containing 5% polyvinylpyrrolidone (molecular mass 40 000 units) and 1% 2-mercaptoethanol. Homogenates were either immediately loaded into the starch gels (12% Sigma starch) or stored frozen at -85C for later electrophoresis. Electrophoresis and staining of enzymes followed standard procedures (Soltis et al., 1983; Werth, 1985; Wendel and Weeden, 1989), but employed the "zymecicle" methodology of Werth (1990). Leucine aminopeptidase (LAP) and esterase (EST) were resolved on the lithium hydroxide system (Werth, 1985); phosphoglucose isomerase (PGI), alcohol dehydrogenase (ADH), and triosephosphate isomerase (TPI) on buffer system 6 of Soltis et al. (1983); phosphoglucomutase (PGM) and 6-phosphogluconic dehydrogenase (6PGD) on the tris-citrate pH 8.0 buffer system of Selander et al. (1971); and aldolase (ALD), isocitrate dehydrogenase (IDH), malate dehydrogenase (MDH), 6PGD, and shikimate dehydrogenase (SKDH) on the morpholine-citrate buffer system (Werth, 1991).

Band interpretation and data analysis
Isozyme band patterns were interpreted using standard principles (Wendel and Weeden, 1989; Murphy et al., 1990). Alleles were designated by numbers representing the migration of the allozymes they encoded, with the lowest number assigned to the most anodal allozyme.

Genetic individuals (genets) were discriminated on the basis of their multilocus isozyme genotypes, and confidence levels for clone assignment estimated using the method of Parks and Werth (1993). Quantitative treatment of the data was carried out using BIOSYS-1 (Swofford and Selander, 1981), Statview Student for the Macintosh (1991), and NTSYS (Rohlf, 1988).


    RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
General
Fifteen loci coding 11 enzymes were consistently scorable and could be used for discrimination of genets. Of these 15 loci, 86.7% (13 loci-- Adh, Est, Idh, Lap, Mdh-1, Mdh-2, Pgi-2, Pgm-2, 6Pgd-1, 6Pgd-2, Skdh, Tpi-1, and Tpi-2) were polymorphic in the study area. The two nonvariable loci were Ald and Pgi-1. Three of the loci, Idh, Mdh-2, and 6Pgd-2, were not variable within the grid subpopulation, and one, Tpi-2, was not variable within the transect subpopulation.

In the initial transect, 38 of the 40 ramets sampled at 50-m intervals exhibited distinct genotypes and could thus be inferred to be different genets (Appendix). Two pairs of adjacent points, numbers 2/3 and 32/33 had identical genotypes, indicating the presence of single clones spanning greater than 50 m. In the two-dimensional grid, 56 distinct genotypes were observed among the 380 ramets sampled (Appendix; Fig. 2). The genotypes were numbered as encountered. In some cases the same genotype was inadvertently given two different numbers, in which case the higher number was omitted, resulting in ordinal values up to 62. Because the different sampling intervals and vastly different linear extent (2000 m vs. 200 m) of the grid and transect might influence their genetic makeup if genotypes were unevenly distributed across the population, most population genetic parameters were computed separately for these two subsamples. Allele frequencies for all loci were calculated separately for the grid and transect and then averaged to obtain population means (Table 1). For the grid subpopulation, in which closely spaced ramets resulted in repeated sampling of genets, allele frequencies were estimated both by tallying all genets (unadjusted) and by using the "round-robin" method, a subsampling routine that can avoid overestimation of frequencies for rare alleles, as described in Parks and Werth (1993). Both techniques resulted in nearly identical frequencies (Table 1), indicating that sufficient genetic variation was present to provide robust genet discrimination.



View larger version (201K):
[in this window]
[in a new window]
 
Fig. 2. Map depicting the identity of genotypes in our sampling grid (200 x 190 m). Each 10 x 10 m square was represented by a sample point in its center. Each pattern represents a different multilocus genotype, i.e., putative genet; numbers correspond to genet numbers in Appendix. Blank squares are regions devoid of Q. havardii ramets.

 

View this table:
[in this window]
[in a new window]
 
Table 1. Estimates of allele frequencies for 15 enzyme loci examined for the grid and transect subpopulations of Quercus havardii. Allele frequencies in grid were estimated using all individuals (standard estimate) and by the "round-robin" method of Parks and Werth (1993).

 
Genotype probabilities
Following Parks and Werth (1993), genotype probability Pgen, i.e., the probability that a zygote formed by combining randomly chosen gametes will have a particular genotype, was calculated for each genotype in the grid as the product across loci of Hardy-Weinberg expected frequencies, based on population mean allele frequencies. The values obtained for Pgen (Appendix), representing the chance that two consecutive ramet samples not belonging to the same genet could independently possess the same genotype, were <0.01 for all putative genets. Therefore, contiguous ramets of identical genotype could be confidently assigned to the same genet. Additionally, for each genet the probability of second encounter PSE, i.e., the chance of observing its genotype one or more additional times among the 56 genets sampled, was computed as 1 - (1 - Pgen). Values for PSE were also small, <0.01 for most genotypes, although for seven genotypes PSE > 0.05. Therefore, most noncontiguous ramets of identical genotype probably belonged to the same genet.

Clone size
Features of size and shape of Q. havardii clones in the grid were evaluated by construction of a map of genets (Fig. 2) and a histogram depicting the distribution in clone size (Fig. 3). Among the 56 clones of the grid, the majority were small (37 with only 1–4 sample points), and a few were large (five with 17 or more sample points). To rule out edge effects as a possible cause of this distribution, the clone size histogram was also plotted using only those clones that did not touch the edge of the grid, resulting in a similar distribution of clone sizes (Fig. 3). The 12 largest clones occupied >70% of the grid area, 40% being attributable to the three largest clones, clones 22, 26, and 32. The largest clone, number 22, was sampled at 70 points and thus occupied an estimated area of 7000 m (at 100 m per sample point).



View larger version (22K):
[in this window]
[in a new window]
 
Fig. 3. Clone size (number of ramets sampled) frequency distribution. The barred columns include all clones, whereas the open columns include only clones that do not come into contact with the edge of the sampling grid.

 
Clone fragmentation
Most of the 56 clones of the grid were apparently continuous and unfragmented. However, 15 of the putative genets appeared to be fragmented, i.e., including one or more ramets isolated from other ramets of the same genotype by at least one sampling point. Notable examples of putatively fragmented clones included genotype 8 (PSE = 0.030), whose two sampled ramets were ~240 m apart, genotype 33 (PSE = 0.069), with five sampled ramets in two groups ~180 m apart, genotype 47 (PSE = 0.017), with four sampled ramets in three groups spread over ~160 m, genotype 31 (PSE = 0.088), with four sampled ramets in two groups ~70 m apart, genotype 29 (PSE = 0.008), with 14 sampled ramets in two main groups ~100 m apart, and genotype 26 (PSE = 0.031), with 25 sampled ramets in two groups ~50 m apart (Fig. 2). Of these 15 putatively fragmented clones, 13 had values for PSE < 0.05 and therefore were unlikely to comprise different genets coincidentally possessing the same genotypes (but see below). Genotypes 31 and 33 had values for PSE > 0.05, such that there is a significant probability that the fragments with these genotypes could be different genets.

In evaluating the likelihood that each of these cases of separated ramet clusters represents clone fragmentation, it is necessary to consider that PSE is the probability for second encounter of an independently formed ramet with one particular genotype. The collective probability that any one or more of the detected genotypes would be reencountered is much greater than PSE because the number of "trials" (pairs of genets) is much greater. The probability of not reencountering a particular genotype is (1 - PSE). The chance that none would be reencountered is the product of (1 - PSE) for all genotypes in the grid subpopulation, which was computed as 0.329. Thus, although the chance of reencountering a particular genotype is low, the chance of reencountering one or more of the 56 genotypes is high, i.e., 0.671. Therefore, it is probable that at least one of the putatively fragmented clones may actually comprise a pair of different genets with identical genotypes, although it remains highly improbable that more than a few of the 15 cases of putative fragmentation are instead cases of distinct genet formation. It cannot be ascertained with certainty which groups of separated putative fragments with identical genotype actually comprise distinct genets. However, widely separated fragments, such as genotype 8, are especially suspect.

The number of fragments was unrelated to clone size (Fig. 4; Spearman rank correlation N = 24, rs = 0.094, p = 0.66). This result apparently deviated from the random expectation, which is a positive relationship, since larger clones should have a higher probability of becoming fragmented.



View larger version (13K):
[in this window]
[in a new window]
 
Fig. 4. Graph of the number of clone fragments as a function of clone size (number of ramets sampled) for all clones that did not come into contact with the edge of the sampling grid.

 
Clone shape
Clones exhibited irregular and variable shapes that presumably reflect their growth history (Fig. 2). Some of the larger clones, such as clones 32, 22, and 6, consisted of massive central regions, free of invading rhizomes of surrounding clones, from which radiated one or more narrower extensions winding into or between other clones. Examples of this can be seen in the northward extension of clone 22 into clone 32 and in the northeastern extension of clone 43 that divides clone 29. However, much of the perimeter of these large clones seemed to form long, fairly straight boundaries with adjacent clones, characterized by little intergrowth as perceived at the scale of sampling employed (i.e., 10-m intervals). Other clones, such as clones 13 and 26, exhibited long narrow convoluted shapes, characterized by intergrowth with adjacent clones and tending to be fragmented. The irregular boundaries of clone 13 with clones 22 and 23 illustrate this. Perhaps the most distinctive shape was exhibited by clone 39, a large C-shaped clone sampled at 18 points, whose central portion was occupied by several other clones, including clones 26 and 47. This could have resulted from the initial growth of clone 39 in competition with its neighbors, or later invasion and displacement of the interior of a circular clone.

Levels and patterns of genetic variability
Genetic variability was quantified for the grid, transect, and combined samples by computing three indices: mean expected heterozygosity (Hexp), proportion of loci polymorphic (P), and mean number of alleles per locus (A) (Table 2). Values of all three indices (Hexp = 0.289, P= 60%, and A = 2.50 for the combined sample) were higher than those reported for most other oaks that have been electrophoretically surveyed (Guttman and Weigt, 1989; Schnabel and Hamrick, 1990; Hamrick, Godt, and Sherman-Broyles, 1992; Berg and Hamrick, 1994) and were also high as compared to the mean values for flowering plants in general (Hexp = 0.113; P= 34.2%; A = 1.53) (Hamrick and Godt, 1989). These values were also high when compared to subgroups defined by life-history characters shared with Q. havardii, e.g., other dicots (Hexp = 0.096; P= 29.0%; A = 1.44), long-lived woody perennials (Hexp = 0.149; P = 50.0%; A = 1.79), temperate plants (Hexp = 0.109; P= 32.6%; A = 1.51), outcrossing wind-pollinated plants (Hexp = 0.148; P= 49.7%; A = 1.79), and other clonal plants (Hexp = 0.103; P= 29.4%; A = 1.47) (Hamrick and Godt, 1989).


View this table:
[in this window]
[in a new window]
 
Table 2. Estimates of genetic variability in Quercus havardii (± 1 SE). Values are reported separately for the grid and transect samples, and for the combined sample of all genets treated as a single population sample.

 
Conformance of genotype proportions to expected values under Hardy-Weinberg equilibrium was tested using the chi-square test (with pooling where appropriate) for each locus separately in the grid and transect subpopulations, and deviations were quantified by computing the fixation index F (GoTable 3). Most loci were in Hardy-Weinberg equilibrium in both the grid and transect subpopulations and showed low positive values for F, ranging from 0.050 to 0.192, that were not statistically different from zero. Five instances, out of 22 valid tests, of statistical departure from Hardy-Weinberg occurred: Adh (p = 0.005), Est (p = 0.046), and 6Pgd-1 (p < 0.001) in the grid subpopulation and Lap (p = 0.001) and Mdh-1 (p = 0.023) in the transect subpopulation. All five of these instances exhibited a deficiency in heterozygotes, as indicated by positive values for F (Adh, F = 0.361, Est, F = 0.160, and 6Pgd-1, F = 0.399, for the grid and Lap, F = 0.392, Mdh-1, F = 0.370, for the transect). Heterozygote deficiencies at these few loci are unlikely to reflect intensive inbreeding, which should result in heterozygote deficiencies at all loci, but instead could reflect population substructuring, i.e., the Wahlund effect (Wahlund, 1928). The possibility that the Q. havardii population might be genetically structured was evaluated in two ways, as follows.


View this table:
[in this window]
[in a new window]
 
Table 3. Values for the inbreeding coefficient F computed for each polymorphic locus in the grid and transect samples. Statistical difference of F from 0 was evaluated using chi-square analysis, and is indicated as follows: ns = not significant; * P < 0.05; ** P < 0.01; *** P < 0.001.

 
Allele-frequency heterogeneity between the grid and transect population samples was addressed by computing F statistics (Wright, 1978) for all loci, and values of FST were statistically evaluated for deviation from 0 using contingency chi-square analysis (Table 4). The amount of divergence between the two samples, as estimated by FST, was generally low with a mean for all loci of FST = 0.063, and more than half of the loci having FST < 0.02. Six loci had FST > 0.02, ranging from 0.02 to 0.359, and were significantly greater than 0 as determined by contingency chi-square analysis. Two loci, 6Pgd-1 (FST = 0.359) and Skdh (FST = 0.075), had values of FST > 0.05, a value widely considered substantial (Hartl, 1980).


View this table:
[in this window]
[in a new window]
 
Table 4. Values of F ST for all loci as an estimate of allele frequency heterogeneity between the grid and transect subsamples of the Q. havardii population. Statistical difference of F ST from 0 was evaluated using contingency chi-square analysis and is indicated as follows: ns = not significant; * P < 0.05; ** P < 0.01; *** P < 0.001.

 
Genetic structure was also evaluated by testing whether genetic relatedness of neighboring clones was greater than expected by chance. Two correlation matrices were constructed comparing each pair of clones in the transect. In one matrix the physical distances between sample points were entered, and in the other pairwise values of Rogers' similarity, computed from the isozyme genotypes, were entered as an index of genetic relatedness. The Mantel Z test (carried out using NTSYS) showed no correlation between the two matrices (r = 0.265, p 0.05). Clones from the grid were not similarly evaluated as there was no consistent method to assign distance between the genets. The results from analysis of Hardy-Weinberg proportions, FST computations, and distance/relatedness comparison combine to indicate that the population is structured (i.e., distribution of genotypes is not entirely random), but only weakly so.


    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Clonal populations often comprise an array of discrete patches that, depending on the growth strategy, may correspond to single genets (e.g., Daehler and Strong, 1994; Montalvo et al., 1997) or may each be composed of several genets (e.g., Parks and Werth, 1993; Berg and Hamrick, 1994). However, the Quercus havardii population addressed in this study comprises a single enormous patch, spread continuously over a large expanse of homogeneously suitable habitat. The use of isozymic markers has allowed for the first time a portion of this expanse to be viewed as genetically distinct individuals comprising a sexually interbreeding population. Despite rhizome growth that results in order-of-magnitude size variation among genets coupled with rarity of sexual recruitment (Muller, 1951; Wiedman, 1960; Pettit, 1979), this population appears to share features with populations of organisms with more determinate growth, i.e., large population size, an outcrossing breeding system in which mating approaches panmixia, and maintenance of a substantial level of genetic variation.

In this respect, Q. havardii parallels other species for which vegetative propagation does not substitute for sexual reproduction, but rather effects horizontal spread of sexually reproduced individuals. Such populations have been shown capable of maintaining significant genetic variation (Thomas and Dale, 1975; Ellstrand and Roose, 1987; Wolf, Haufler, and Sheffield, 1988; Wolf, Sheffield, and Haufler, 1991; Aspinwall and Christian, 1992; Berg and Hamrick, 1994; Lokker et al., 1994; Sipes and Wolf, 1997; Montalvo et al., 1997) despite possibly lowered effective population sizes in clonally reproducing species relative to species that reproduce only sexually. It is unlikely that population size constrains maintenance of genetic diversity in this species. Assuming that the density of ~15 genets per hectare in the grid is representative for populations of this species, the number of individuals across the species range (2.3 million hectares), much of which was contiguous up until the early 20th century, can be estimated as exceeding 30 million. Furthermore, the apparently low degree of genetic structure indicates that neighborhood sizes are very large, as indicated for other Quercus species (Berg and Hamrick, 1994; Montalvo et al., 1997; but see Sork, Huang, and Wiener, 1993).

Attributes of clonal architecture are highly variable among and sometimes within plant species (Ellstrand and Roose, 1987). Clone size may range from <1 m for some herbaceous species to >43 ha for a clone of Populus tremuloides Michaux (Kemperman and Barnes, 1976). Prior to the present study, the size range of Quercus havardii clones was unknown, although estimates of up to ~15 m had been made based on the sizes of discrete patches that occur in the eastern portion of the species range and of phenologically differentiated patches within the expansive communities (Muller, 1951). Here we have found that these previous estimates were conservative, as indicated by large clones such as number 22, which exceeded 150 m in maximum breadth and occupied ~7000 m.

Morphologies and growth strategies of clonal organisms vary across a broad spectrum from regularly shaped radiating circles of densely clumped ramets with minimal intergrowth (phalanx morphology) to irregularly shaped, meandering, and/or ramifying clones with more widespread ramets and a tendency for intergrowth and fragmentation (guerilla morphology) (Lovett-Doust, 1981; Silander, 1985). Quercus havardii clones, even the smaller ones, comprise numerous closely spaced ramets that occupy significant areas free of ramets from other clones (or other species) and tend not to be fragmented, thus possessing a fundamentally phalanx morphology. Phalanx morphology of Q. havardii is further suggested by anecdotal observation of its subterranean structure, usually revealed by exposure from wind erosion, rarely from purposeful excavations (J. C. Zak, Texas Tech University, personal communication). Such observations suggest a hierarchical system comprising groups of branches produced at regular intervals from a thick subterranean rhizome. However, to varying degrees individual clones also possess attributes associated with guerilla morphology. A meandering and branching form characterizes at least a portion of many of the clones, including the arms radiating from the larger ones, and some clones exhibit fragmentation. It remains uncertain, and a possible focus of future investigation, as to what kinds of spatial interactions occur between clones at a finer scale than addressed herein.

While the phenomenon of asexual proliferation through rhizome growth is readily observed in Quercus havardii, the nature of sexual reproduction and recruitment in this species is very poorly understood. Although large numbers of acorns are produced in ~3 out of every 10 yr (Muller, 1951; Pettit, 1979), Q. havardii seedlings rarely become established in the wild (Muller, 1951; Wiedman, 1960; Pettit, 1979), as virtually every available safe sight is occupied by existing genets, severely constraining the opportunity for seedling establishment. Moreover, disturbance by animals has been shown to facilitate establishment of herbaceous plant seedlings by opening sites in a deep layer of leaf-litter (Dhillion et al., 1994; Jeffery, 1998; Willig and McGinley, 1998). Perhaps disturbance may play a parallel role in recruitment of Q. havardii seedlings, but this remains to be observed directly.

Processes and histories of recruitment of long-lived perennials such as forest trees and other woody plants are difficult to evaluate because of the large time scales involved. Vegetative proliferation compounds the problem by blurring the distinction between individuals. The recognition of Q. havardii genets through isozyme fingerprints achieved in the present study narrows the universe of possibilities and clarifies the nature of questions that can be addressed in future studies. The large density of distinct genets (15 genets/ha) comprising an interbreeding sexual population with conventional attributes implicates an important and continuing, if episodic, role for sexual recruitment. The establishment of the population must have involved a complex dynamic that occurred during antiquity, the timing and progression of which can only be speculated, as the ages of the clones and the biotic community they so thoroughly dominate are unknown. The highly variable sizes and shapes of the Q. havardii clones, as they collectively occupy nearly the totality of a vast surface, must have largely been influenced by their competitive interactions as they came to abut. However, it remains unknown whether the spatial features of the genets revealed are long-enduring, having reached a stasis, are continually changing but in dynamic equilibrium, or are subject to drastic chaotic changes following episodic occurrence of environmental extremes.


    FOOTNOTES
 
1 The authors thank Kathy Stamback, Jeff Jump, and especially Thao Le for assistance with isozyme electrophoresis; Arndt Laemmerzahl, and Michele Holland for assistance with field work; John Zak and Mike Willig for input into the study design; and Mr. Johnnie Fitzgerald for permission to conduct the study on his land. This study was conducted in partial completion of a Master's degree from Texas Tech University (SGM). Financial support was provided by the Office of Research Services at Texas Tech University in the form of a summer grant-in-aid to SGM and CRW.


View this table:
[in this window]
[in a new window]
 
Appendix. Distinct multilocus genotypes detected in the grid (N = 56) and transect (N = 40) subpopulations. Individual locus genotypes are provided as a pair of single-digit numerical allele designations for the 13 polymorphic loci (see Methods; monomorphic loci Ald and Pgi-1 ). Number of loci heterozygous and probabilities for occurrence (P gen) and second encounter (P SE) were computed for the grid only. Blank entries indicate missing data.

 

View this table:
[in this window]
[in a new window]
 
Appendix. Continued.

 
Back 2 Author for correspondence. Back


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Aspinwall, N., and T. Christian.1992Clonal structure, genotypic diversity, and seed production in populations of Filipendula rubra (Rosaceae) from the northcentral United States. American Journal of Botany 79: 294–299. [CrossRef][ISI]

Berg, E. E., and J. L. Hamrick.1994Spatial and genetic structure of two sandhills oaks: Quercus laevis and Quercus margaretta (Fagaceae). American Journal of Botany 81: 7–14.

Correll, D. S., and M. C. Johnston.1979Manual of the vascular plants of Texas. University of Texas at Dallas, Richardson, TX.

Daehler, C. S., and D. R. Strong.1994Variable reproductive output among clones of Spartina alternifolia (Poaceae) invading San Francisco Bay, California: the influence of herbivory, pollination, and establishment site. American Journal of Botany 81: 307–313. [CrossRef][ISI]

Dhillion, S. S., M. A. McGinley, C. F. Friese, and J. C. Zak.1994Construction of sand-shinnery oak communities of the Llano Estacado: animal disturbances, plant community structure, and restoration. Restoration Ecology 2: 51–61.

Ellstrand, N. C., and M. L. Roose.1987Patterns of genotypic diversity in clonal plant species. American Journal of Botany 74: 123–131. [CrossRef][ISI]

Guttman, S.I., and L. A. Weigt.1989Electrophoretic evidence of relationships among Quercus (oaks) of eastern North America. Canadian Journal of Botany 67: 339–351. [CrossRef]

Hamrick, J. L., and M. J. W. Godt.1989Allozyme 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.

———, ———, and S. L. Sherman-Broyles.1992Factors influencing levels of genetic diversity in woody plant species. New Forests 5: 95–124.

Hartl, D.1980Principals of population genetics. Sinauer, Sunderland, MA.

Jeffery, J. M.1998The effect of oak-leaf litter on grass seedling germination and establishment in a sand-shinnery oak community of West Texas. Master's thesis, Texas Tech University, Lubbock, TX.

Kemperman, J. A., and B. V. Barnes.1976Clone size in American aspens. Canadian Journal of Botany 54: 2603–2607.

Lokker, C., D. Susko, L. Lovett-Doust, and J. Lovett-Doust.1994Population genetic structure of Vallisneria americana, a dioecious clonal macrophyte. American Journal of Botany 81: 1004–1012. [CrossRef][ISI]

Lovett Doust, L.1981Population dynamics and local specialization in a clonal perennial (Ranunculus repens). I. The dynamics of ramets in contrasting habitats. Journal of Ecology 69: 743–755. [CrossRef]

Montalvo, A. M., S. G. Conard, M. T. Conkle, and P. D. Hodgkiss.1997Population structure, genetic diversity, and clone formation in Quercus chrysolepis (Fagaceae). American Journal of Botany 84: 1553–1564. [Abstract]

Muller, C. H.1951The significance of vegetative reproduction in Quercus. Madroño 11: 129–137.

Murphy, R. W., J. W. Sites, Jr., D. G. Buth, and C. H. Haufler.1990Proteins I: Isozyme electrophoresis. In D. Hillis and C. Moritz [eds.], Molecular systematics, 44–124. Sinauer, Sunderland, MA.

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

Parks, J. C., and C. R. Werth.1993A study of spatial features in a population of Bracken Fern, Pteridium aquilinum (L.) Kuhn. American Journal of Botany 80: 1018–1025.

Pettit, R. D.1979Effects of picloram and tebuthiuron pellets on sand shinnery oak communities. Journal of Range Management 32: 196–200. [CrossRef][ISI]

Rohlf, F. J.1988NTSYS-pc: numerical taxonomy and multivariate analysis system. Exeter Publishing, Setauket, NY.

Schnabel, A., and J. L. Hamrick.1990Comparative analysis of population genetic structure in Quercus macrocarpa and Q. gambelii (Fagaceae). Systematic Botany 15: 240–251. [CrossRef][ISI]

Scifres, C. J.1972Herbicide interactions in control of sand shinnery oak. Journal of Range Management 25: 386–389. [CrossRef][ISI]

Selander, R. K., M. H. Smith, S. Y. Yang, W. E. Johnson, and J. B. Gentry.1971Biochemical polymorphism and systematics in the genus Peromyscus. I. Variations in the old-field mouse (P. polionotus). University of Texas Publication 103: 49–90.

Silander, J. A.1985Microevolution in clonal plants. In J. B. C. Jackson, L. W. Buss, and R. E. Cook [eds.], Population biology and evolution of clonal organisms, 107–152, Yale University Press, New Haven, CT.

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

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

Sork, V. L., S. Huang, and E. Wiener.1993Macrogeographic and fine-scale genetic structure in a North American oak species Quercus rubra L. Annales des Sciences Forestieres 50(S): 128–136.

Statview Student.1991Abacus Concepts, Inc. Berkeley, CA.

Swofford, D., and R. Selander.1981BIOSYS-1. A FORTRAN program for the comprehensive analysis of electrophoretic data in population genetics and systematics. Journal of Heredity 72: 281–283. [Abstract/Free Full Text]

Thomas, A. G., and H. M. Dale.1975The role of seed reproduction in the dynamics of established populations of Hieracium floribundum and a comparison with that of vegetative reproduction. Canadian Journal of Botany 53: 3022–3031.

Wahlund, S.1928Zusammensetzung von Population und Korrelationsersceinung vom Standpunkt der Vererbungslehre aus Betrachtet. Hereditas 11: 65–105. [ISI]

Wendel, J. F., and N. F. Weeden.1989Visualization and interpretation of plant isozymes. In D. E. Soltis and P. S. Soltis [eds.], Isozymes in plant biology, 5–45. Dioscorides Press, Portland, OR.

Werth, C. R.1985Implementing an isozyme laboratory at a field station. Virginia Journal of Science 36: 53–76.

———.1990Zymecicles: pre-prepared frozen isozyme assays. Isozyme Bulletin 23: 109.

———.1991Isozyme studies on the Dryopteris "spinulosa" complex. I. The origin of the log fern D. celsa. Systematic Botany 16:446–461.

Wiedman, V.1960Preliminary ecological study of the shinnery oak area of western Oklahoma. M. S. Thesis. Oklahoma University, Norman, OK.

Willig, M. R., and M. A. McGinley.1998Response and recovery of animal populations and communities to disturbance. In L. R. Walker [ed.], Ecosystems of disturbed ground, Elsevier Press, Amsterdam, in press.

Wolf, P., C. Haufler, and E. Sheffield.1988Electrophoretic variation and mating system of the clonal weed Pteridium aquilinum (L.) Kuhn (bracken). Evolution 42: 1350–1355. [CrossRef][ISI]

———, E. Sheffield, and C. Haufler.1991Estimates of gene flow, genetic substructure and population heterogeneity in bracken (Pteridium aquilinum). Botanical Journal of the Linnaean Society 42: 407–423. [CrossRef]

Wright, S.1978Evolution and genetics of populations, vol. 4, Variability within and among natural populations. University of Chicago Press, Chicago, IL.

Zhang, Q.1996Fungal community structure and microbial biomass in a semi-arid environment: roles in root decomposition, root growth, and soil nutrient dynamics. Ph.D. dissertation, Texas Tech University, Lubbock, TX.





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 ISI Web of Science (20)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Mayes, S. G.
Right arrow Articles by Werth, C. R.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Mayes, S. G.
Right arrow Articles by Werth, C. R.
Agricola
Right arrow Articles by Mayes, S. G.
Right arrow Articles by Werth, C. R.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS