|
|
||||||||
a Ecology Program, Department of Biological Sciences, Texas Tech University, Lubbock, Texas 794093131
| ABSTRACT |
|---|
|
|
|---|
Key Words: allozyme clone Fagaceae genetic variation isozyme population structure Quercus
| INTRODUCTION |
|---|
|
|
|---|
|
| MATERIALS AND METHODS |
|---|
|
|
|---|
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 |
|---|
|
|
|---|
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.
|
|
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 14 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).
|
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.
|
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).
|
|
|
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 |
|---|
|
|
|---|
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 |
|---|
|
|
| REFERENCES |
|---|
|
|
|---|
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: 714.
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: 307313. [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: 5161.
Ellstrand, N. C., and M. L. Roose.1987Patterns of genotypic diversity in clonal plant species. American Journal of Botany 74: 123131. [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: 339351. [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, 4363. Sinauer, Sunderland, MA.
, , and S. L. Sherman-Broyles.1992Factors influencing levels of genetic diversity in woody plant species. New Forests 5: 95124.
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: 26032607.
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: 10041012. [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: 743755. [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: 15531564. [Abstract]
Muller, C. H.1951The significance of vegetative reproduction in Quercus. Madroño 11: 129137.
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, 44124. Sinauer, Sunderland, MA.
Nei, M.1978Estimation of average heterozygosity and genetic heterozygosity and genetic distance from a small number of individuals. Genetics 89: 583590.
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: 10181025.
Pettit, R. D.1979Effects of picloram and tebuthiuron pellets on sand shinnery oak communities. Journal of Range Management 32: 196200. [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: 240251. [CrossRef][ISI]
Scifres, C. J.1972Herbicide interactions in control of sand shinnery oak. Journal of Range Management 25: 386389. [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: 4990.
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, 107152, 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: 401409. [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: 927. [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): 128136.
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: 281283.
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: 30223031.
Wahlund, S.1928Zusammensetzung von Population und Korrelationsersceinung vom Standpunkt der Vererbungslehre aus Betrachtet. Hereditas 11: 65105. [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, 545. Dioscorides Press, Portland, OR.
Werth, C. R.1985Implementing an isozyme laboratory at a field station. Virginia Journal of Science 36: 5376.
.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:446461.
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: 13501355. [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: 407423. [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.
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |