|
|
||||||||
Population Biology |
2Department of Plant Biology, University of Georgia, Athens, Georgia 30602 USA; 3Invasive Plant Research Laboratory, USDA/ARS, Ft. Lauderdale, Florida, 33314 USA; 4Institute of Ecology, University of Georgia, Athens, Georgia 30602 USA; 5Departments of Plant Biology and Genetics, University of Georgia, Athens, Georgia 30602 USA
Received for publication December 1, 2003. Accepted for publication May 6, 2004.
| ABSTRACT |
|---|
|
|
|---|
Key Words: allozymes clonal structure coastal plants gene flow population genetics sea oats southeastern United States Uniola paniculata
| INTRODUCTION |
|---|
|
|
|---|
Large-scale population structuring may also occur if genetic patterns differ between the center of a species' range and the margins of its distribution. There is currently some debate as to whether populations near the center should display more or less genetic variability than those at the edge of the range. The classical view, often referred to as the CentralMarginal Model, is that central populations are more genetically variable than marginal populations because of fewer ecological niches at the margins (da Cunha and Dobzhansky, 1954
) or a greater importance of heterosis in central populations (Carson, 1959
; Wallace, 1984
). In contrast, Brussard (1984)
proposed that selection at environmentally variable margins may promote genetic diversity. In reviewing the literature on Drosophila species, Brussard (1984)
found declines in inversion polymorphisms, which are presumably under selection, from the center to the margins of the species' ranges, but no change in supposedly neutral allozyme variation. Patterns of genetic variation between central and marginal populations may thus depend on the degree to which the markers are under selection, rates of gene flow into marginal populations, and the amount of habitat variability and rates of population growth in central and marginal areas.
Coastal plants are particularly appropriate for addressing questions regarding patterns of spatial genetic structuring and comparisons of genetic variation between central and marginal areas because of their long, narrow, and naturally fragmented habitat. How and whether the unique geographic distribution of coastal plants affects patterns of gene flow remains largely unknown because of the paucity of population genetic studies in these species. In addition, many coastal dune plants reproduce both asexually and sexually (Ranwell, 1972
). While vegetative propagation (Pemadasa and Lovell, 1976
) and flowering (Pemadasa and Lovell, 1974
) have been examined separately for some species, the relative importance of asexual vs. sexual reproduction in coastal dune plants is generally unknown. The objective of this study was to investigate the population genetics of the coastal dune grass Uniola paniculata L. at three spatial scales: among regions, within and among populations, and among clones within patches.
Uniola paniculata (sea oats) is the dominant dune-stabilizing grass in southeastern USA. This native plant occurs on the primary dunes of beaches and barrier islands from Virginia, south along the Atlantic Ocean, into the Bahamas, and along the Gulf of Mexico coast to Veracruz, Mexico (Wagner, 1964
). This range spans more than 5000 km of coastline. However, the width of the habitat within a dune is quite narrowgenerally less than 200 m. Thus, this species has an essentially linear distribution that is naturally fragmented: patches of beach occur on barrier islands or on mainland beaches that are isolated by waterways. Additional fragmentation occurs both by natural disturbances such as hurricanes and anthropogenic disturbances such as coastal housing developments.
Uniola is outcrossing and wind-pollinated, and its seeds (caryopses) are dispersed by wind and water (Harper and Seneca, 1974
; Hester and Mendelssohn, 1987
). While these factors should promote gene flow and decrease genetic differences among populations, its linear, fragmented distribution should limit gene flow and increase genetic structure. In addition, Seneca (1972)
found that Virginia and North Carolina populations of Uniola required cold stratification for germination, while Florida populations did not. This suggests that selection may differ among climatically distinct regions of the range of Uniola, further increasing among-population differences. Based on this information, we hypothesized that Uniola would have more genetic structure than other outcrossing species, that genetic distance would increase with geographic distance, that genetically similar populations would cluster within regions, and that genetic diversity would be lowest at the margin of the species' range.
In addition to sexual reproduction, Uniola reproduces clonally by the production of rhizomes (Wagner, 1964
). Clonal reproduction is thought to be extensive for this species (Harper and Seneca, 1974
), and rates of germination and seedling emergence and establishment in the field tend to be low (Wagner, 1964
; Franks, 2003
). We hypothesized that extensive ramet production and low seed germination and establishment would cause populations to be dominated by a few large clones.
| MATERIALS AND METHODS |
|---|
|
|
|---|
|
|
Electrophoresis
Leaf tissue samples were kept cool, taken back to the laboratory, and crushed for protein extraction using liquid nitrogen, sea sand, and Camellia buffer (Wendel and Parks, 1982
). The extracted material was absorbed onto Whatman 3MM filter paper wicks. These wicks were stored at 70°C until needed for electrophoresis. Extracted proteins were run on 9.5% starch gels. Seventeen enzymes were examinedPGI, MNR, PGM and ADH on buffer system 34; MDH, IDH, AK, SKDH, MPI and F16 on buffer system 11; DIA, AAT, TPI, FE, GDH, and ME on buffer system 8; and 6-PGD on buffer system 4 (Soltis et al., 1983
for all buffers and stains except AAT, ADH, AK, DIA, MNR, MPI [Cheliak and Pitel, 1984
]). Gels were then stained with enzyme-specific stains according to the protocols above and scored for allozyme banding patterns. Interpretations of the banding patterns follow the recommendations of Weeden and Wendel (1989)
and Kephart (1990)
for polyploid species.
Analyses
Population genetic statistics were calculated using the program LYNSPROG developed by M. D. Loveless (College of Wooster, Wooster, Ohio, USA) and A. F. Schnabel (University of Indiana, South Bend, Indiana, USA). At the species (subscript s) and population (subscript p) levels, we estimated percentage polymorphic loci (p), mean number of alleles per polymorphic locus (AP), observed heterozygosity (Ho), and expected heterozygosity (He = 1
pi2), which estimates genetic diversity. Standard errors for within-population parameters were obtained by averaging across all populations.
Nei's (1972)
genetic distance and identity statistics were calculated for each pair of populations. Genetic distance was used to perform a cluster analysis using the unweighted pairwise groups method using arithmetic averages (UPGMA) with the program Phylip (version 3.2; Felsenstein, 1989
). Geographic distance was calculated between all pairs of populations. For adjacent populations, we calculated the distance between the sites using the latitude and longitude coordinates of each site. For nonadjacent populations, we calculated the distance along the coast between them by summing the distances connecting adjacent populations between them. This method assumes that the relevant geographic distance among populations is the distance along the coast, because noncoastal habitat is unsuitable for this species. To examine the relationship between geographic distance and Rousset's (1997)
genetic distance (FST/[1 FST]), we performed Mantel tests and reduced major axis (RMA) regression on natural log-transformed data using the program IBD (version 2.03) by A. J. Bohonak (Bohonak, 2002
). We used Rousset's genetic distance and RMA regression because these statistics are appropriate for genetic/geographic distance comparisons (Bohonak, 2002
).
To partition genetic diversity into geographical components, we conducted a hierarichal analysis of genetic structure. Genetic structure, GST, was calculated following Nei (1987)
. Populations were grouped by two coasts (Atlantic and Gulf) or by five smaller regions, and total GST was partitioned into components among regions and among populations within regions.
We examined fine-scale genetic structure by investigating both clonal structure and patterns of relatedness among individuals within a patch. To assess clonal structure, we first determined the multilocus genotypes for each individual within each patch. We assessed genetic diversity first as the number of genets per number of ramets sampled (G/N). Our second measure of genotypic diversity is Simpson's index: D = 1
([ni{ni 1}]/[N{N 1}]); ni = number of ramets of the ith genet and N = total number of ramets (Pielou, 1969
). We also mapped multilocus genotypes to examine spatial patterns of clonal structure. To examine fine-scale patterns of genetic relatedness, we estimated the coefficient of coancestry (Cockerham, 1969
) between all pairs of genetic individuals within patches using the program fijAnal3D developed by J. Nason (Department of Botany, Iowa State University, Ames, Iowa, USA). Coancestry (
ij) was estimated following Parker et al. (2001)
. We generated 95% confidence intervals around the hypothesis of no spatial genetic structure by bootstrapping using the program BS_fij3D developed by J. Nason.
| RESULTS |
|---|
|
|
|---|
|
|
A similar lack of strong regional genetic pattern was found by organizing the populations into a UPGMA tree (Fig. 2). This clustering dendrogram based on Nei's genetic distances revealed that populations spatially near each other tended to be genetically similar. However, there were notable exceptions. For example, Virginia Beach, Virginia, and South Padre Island, Texas, were the two most widely separated populations geographically, yet are fairly similar genetically.
|
Clonal study
Clonal structure and diversity varied widely among the eight patches sampled (Table 4). One patch (Cape Hatteras A) contained only two multilocus genotypes for 60 ramets sampled. Another patch (Emerald Coast A) had 13 genotypes. The average number of ramets/genet was 11.8 (G/N = 0.13), and mean clonal diversity (measured as Simpson's diversity index) was 0.60 (Table 4). The scale of spatial autocorrelation tended to be fairly small (<3 m) based on spatial analysis of coancestry coefficients. At distances greater than 3 m, ramets were not any more likely to be genetically similar than expected by chance. The probability that two individuals shared a multilocus genotype by chance alone (Harada et al., 1997
) was fairly small (approximately 0.097 for both sites).
|
| DISCUSSION |
|---|
|
|
|---|
In one of the few genetic studies of coastal plants, Clausing et al. (2000)
examined genetic variation in two coastal species (Cakile maritima and Eryngium maritimum) in Europe. For these species, genetic distance and geographic distance along the coast were well correlated. They also found regional differences between the Atlantic and Mediterranean populations of both species, which contrasts with the lack of regional differences found for Uniola.
Regional genetic patterns
Based on the CentralMarginal Model (da Cunha and Dobzhansky, 1954
), we predicted that genetic diversity should increase southward from the northern edge of the species range. Our results did not support this prediction. Instead, some populations near the northern margin of the species range (e.g., Cape Hatteras) were among the most diverse. Also, diversity was low in all populations from the Gulf of Mexico west of the Florida panhandle, possibly due to natural or anthropogenic disturbance.
The results of this study contrast with predictions of the CentralMarginal Model as well as with some empirical work. For example, Robinson et al. (2002)
showed a decline in genetic diversity toward the edge of the range of a mite species. A possible explanation for this discrepancy is that the Central Marginal Model, developed for mobile insect species, does not hold for plants that are sessile except for dispersal of pollen and seeds. It is also possible that only loci under selection increase in diversity away from the margin. This hypothesis is supported by Brussard (1984)
, who found that in Drosophila, inversion polymorphisms, which are presumably under selection, supported the model but that allozyme data did not. Other theoretical work (Mayr, 1970
; Kirkpatrick and Barton, 1997
) suggests that gene flow into marginal populations may offset drift and directional selection, minimizing genetic differences between central and marginal populations. Our results are consistent with the hypothesis of relatively high gene flow into marginal populations.
Previous studies have examined genetic patterns for several vertebrate and invertebrate species specifically in the region of the Florida peninsula. Many of these studies have found genetic disjunctions between populations that are located on the Gulf of Mexico and the Atlantic coasts and between south Florida and more northern populations, though other studies failed to find such a disjunction (reviewed in Avise, 2000
, pages 224242). Avise (2000)
argues that historical biogeographic factors were responsible for establishing barriers to gene flow and creating the patterns of genetic disjunctions seen in this region. In our study, we could explain interpopulation genetic variation by differences among regions, but we did not find a regional disjunction between Atlantic and Gulf coasts or between northern and southern areas, as has been found for other species. There are several possible reasons for this discrepancy. First, studies in this region cited in Avise (2000)
were all of animals, while our study was of a coastal plant. Although many of the animals studied were highly mobile, the outcrossing and wind-dispersed plant we examined may be better able to overcome potential geographic barriers to gene flow than animal species. Second, our study utilized biparentally inherited markers (allozymes) rather than uniparentally inherited markers, such as mtDNA, used in the animal studies. Uniparentally inherited markers allow additional power to trace genetic lineages, so it is possible that a study using mtDNA or cpDNA would find genetic disjunctions in Uniola that were not detected using allozymes. Third, while the overall genetic patterns did not show a strong GulfAtlantic or northsouth disjunction for Uniola, at least one locus, DIA, did show a geographic cline in allele frequencies.
Fine-scale genetic structure
The results of the clonal study indicate that patterns of clonal diversity and structure can vary widely within a relatively small area. Both sites (Cape Hatteras, North Carolina, and Grayton Beach, Florida) had patches with both low and high clonal diversity. Overall, clonal diversity tended to be higher than expected, averaging 0.600 (Simpson's diversity index, Table 4). In most areas, there were many different genotypes, and ramets spaced more than 4 m apart were unlikely to belong to the same clone. This result indicates that both clonal and sexual reproduction are likely important for this species.
The finding that many of the patches had a high diversity of clones was somewhat surprising in light of the fact that clonal reproduction in this species is known to be extensive (Harper and Seneca, 1974
) and that, while germination rates can be high under laboratory conditions (Seneca, 1972
; Hester and Mendelssohn, 1987
), seedling emergence in the field is often very low (Wagner, 1964
; Franks, 2003
). There are, however, occasional periods in which a large number of seedlings establish. Such a localized flush of seedlings has been observed on Keywadin Island (S. J. Franks, personal observation). This pattern of sporadic and patchy establishment could also help to explain geographic patterns of genetic structure, with rare long-distance dispersal events reducing genetic structure but clonal reproduction and recruitment of seedlings into localized patches maintaining differences among populations.
Widen et al. (1994)
surveyed clonal structure in 45 plant species and found an average of 26% multilocus genotypes, which is lower than was found in this study for Uniola. Other studies have shown that clonal structure varies widely both among species and among populations within species (Kudoh et al., 1999
; Reusch et al., 2000
; Stehlik and Holderegger, 2000
; Suyama et al., 2000
). Variation in clonal diversity can be related to such biotic factors as breeding system and dispersal patterns and to such abiotic factors as nutrient heterogeneity and disturbance (Kudoh et al., 1999
; Stehlik and Holderegger, 2000
; Suyama et al., 2000
; Parker et al., 2001
). Kudoh et al. (1999)
found that for the temperate woodland herb Uvularia perfoliata, clonal structure was strongly dependent on disturbance patterns, with patches in canopy gap habitats containing a few large clones but many different genotypes and patches in closed canopy habitats comprised of a single multilocus genotype. In contrast, Parker et al. (2001)
found that although fire history altered population age structure and genetic structure in seedlings of Pinus clausa, disturbance did not detectably affect adult genetic structure. Disturbance may play a role in the clonal structure of Uniola, especially because such factors as burial, erosion, and storm overwash are important in its distribution and population dynamics (Wagner, 1964
; Franks and Peterson, 2003
). However, causes of variation in clonal structure among Uniola patches could not be determined from this study and remain uncertain.
In summary, several conclusions can be drawn from the patterns of genetic variation and structure in Uniola paniculata. First, the degree of genetic differentiation among populations and the relationship between geographic and genetic distance indicate that this species experiences moderate barriers to gene flow, which is consistent with expectation for a species with a linear, fragmented distribution. Second, there were no discreet disjunctions among regions, and patterns of genetic similarity did not always follow geographic location, potentially due to occasional long-distance dispersal, disturbance events, or to anthropogenic factors. Third, genetic diversity was not lowest at the margins of the species' range, which does not support the predictions of the CentralMarginal Model. Finally, fine-scale genetic patterns indicate that clonal diversity is highly variable and that both asexual reproduction via ramets and sexual reproduction via seeds are likely important for regeneration in this species. These results aid our understanding of the population genetics of a species with a unique distribution and can be used in conservation and restoration efforts involving this important coastal plant. Experiments designed to elucidate the degree of local adaptation and the optimum genotypes for use in restoration would complement the findings of this study and aid coastal management efforts.
| FOOTNOTES |
|---|
| LITERATURE CITED |
|---|
|
|
|---|
Berg E. E. J. L. Hamrick 1997 Quantification of genetic diversity at allozyme loci. Canadian Journal of Forest Research 27: 415-424
Bohonak A. J. 2002 IBD (Isolation by Distance): a program for analyses of isolation by distance. Journal of Heredity 93: 153-154
Brussard P. F. 1984 Geographic patterns and environmental gradients: the CentralMarginal Model in Drosophila revisited. Annual Review of Ecology and Systematics 15: 25-64[CrossRef][ISI]
Carson H. L. 1959 Genetic conditions which promote or retard the formation of species. Cold Spring Harbor Symposium on Quantitative Biology 24: 87-103[ISI][Medline]
Cheliak W. M. J. A. Pitel 1984 Techniques for starch gel electrophoresis of enzymes from forest tree species. Petawawa National Forestry Institute, Canadian Forestry Service, Chalk River, Ontario, Canada, Information Report P1-X-42
Clausing G. K. Vickers J. W. Kadereit 2000 Historical biogeography in a linear system: genetic variation of sea rocket (Cakile maritima) and sea holly (Eryngium maritimum) along European coasts. Molecular Ecology 9: 1823-1833[CrossRef][Medline]
Cockerham C. C. 1969 Variance of gene frequencies. Evolution 23: 72-84[CrossRef][ISI]
Da Cunha A. B. T. Dobzhansky 1954 A further study of chromosomal polymorphism in Drosophila willistoni in relation to environment. Evolution 8: 119-134[CrossRef][ISI]
Felsenstein J. 1989 PHYLIPPhylogeny Inference Package. Cladistics 5: 164-166
Franks S. J. 2003 Facilitation in multiple life-history stages: evidence for nucleated succession in coastal dunes. Plant Ecology 168: 1-11[CrossRef][ISI]
Franks S. J. C. J. Peterson 2003 Burial disturbance leads to facilitation among coastal dune plants. Plant Ecology 168: 13-21[CrossRef][ISI]
Futuyma D. 1986 Evolutionary biology. Sinauer, Sunderland, Massachusetts, USA
Godt M. J. W. J. L. Hamrick 1998 Allozyme diversity in the grasses. In G. P. Cheplick [ed.], Population biology of the grasses, 1129. Cambridge University Press, Cambridge, UK
Hamrick J. L. M. J. W. Godt 1996 Effects of life history traits on genetic diversity in plant species. Philosophical Transactions of the Royal Society London, series B 351: 1291-1298[CrossRef]
Harada Y. S. Kawano Y. Iwasa 1997 Probability of clonal identity: inferring the relative success of sexual versus clonal reproduction from spatial genetic patterns. Journal of Ecology 85: 591-600[CrossRef]
Harper J. L. E. D. Seneca 1974 A preliminary study of flowering in Uniola paniculata along the North Carolina coast. Bulletin of the Torrey Botanical Club 101: 7-13[CrossRef][ISI]
Hedrick P. W. 2000 Genetics of populations. Jones & Bartlett, Sudbury, Massachusetts, USA
Hester M. W. I. A. Mendelssohn 1987 Seed production and germination response of four Louisiana populations of Uniola paniculata (Gramineae). American Journal of Botany 74: 1093-1101[CrossRef][ISI]
Kephart S. R. 1990 Starch gel electrophoresis of plant isozymes: a comparative analysis of techniques. American Journal of Botany 77: 693-712[CrossRef][ISI]
Kirkpatrick M. N. H. Barton 1997 Evolution of a species' range. American Naturalist 150: 1-23[CrossRef][ISI]
Kudoh H. H. Shibaike H. Takasu D. F. Whigham S. Kawano 1999 Genetic structure and determinants of clonal structure in a temperate deciduous woodland herb, Uvularia perfoliata. Journal of Ecology 87: 244-257[CrossRef]
Mayr E. 1970 Populations, species, and evolution. Belknap Press, Cambridge, Massachusetts, USA
Nassar J. M. J. L. Hamrick T. H. Fleming 2001 Genetic variation and population structure of the mixed-mating cactus, Melocactus curvispinus (Cactaceae). Heredity 87: 69-79[CrossRef][ISI][Medline]
Nei M. 1972 Genetic distance between populations. American Naturalist 106: 283-292[CrossRef][ISI]
Nei M. 1987 Molecular evolutionary genetics. Columbia University Press, New York, New York, USA
Parker K. C. J. L. Hamrick A. J. Parker J. D. Nason 2001 Fine-scale genetic structure in Pinus clausa (Pinaceae) populations: effects of disturbance history. Heredity 87: 99-113[CrossRef][ISI][Medline]
Pemadasa M. A. P. H. Lovell 1974 Factors controlling flowering time of some dune annuals. Journal of Ecology 62: 869-880[CrossRef]
Pemadasa M. A. P. H. Lovell 1976 Effects of timing of life-cycle on vegetative growth of some dune annuals. Journal of Ecology 64: 213-222[CrossRef]
Pielou E. C. 1969 An introduction to mathematical ecology. Wiley Interscience, New York, New York, USA
Ranwell D. S. 1972 Ecology of salt marshes and sand dunes. Chapman and Hall, New York, New York, USA
Reusch T. B. H. W. T. Stam J. L. Olsen 2000 A microsatellite-based estimation of clonal diversity and population subdivision in Zostera marina, a marine flowering plant. Molecular Ecology 9: 127-140[CrossRef][Medline]
Robinson M. T. A. R. Weeks A. A. Hoffman 2002 Geographic patterns of clonal diversity in the earth mite species Penthaleus major with particular emphasis on species margins. Evolution 56: 1160-1167[ISI][Medline]
Rousset F. 1997 Genetic differentiation and estimation of gene flow from F-statistics under isolation by distance. Genetics 145: 1219-1228[Abstract]
Seneca E. D. 1972 Germination and seedling response of Atlantic and Gulf Coast populations of Uniola paniculata. American Journal of Botany 59: 290-296[CrossRef][ISI]
Soltis D. E. C. H. Haufler D. C. Darrow G. J. Gastony 1983 Starch gel electrophoresis of ferns: a compilation of grinding buffer, gel and electrode buffers, and staining schedules. American Fern Journal 73: 9-27[CrossRef][ISI]
Stehlik I. R. Holderegger 2000 Spatial genetic structure and clonal diversity of Anemone nemorosa in late successional deciduous woodlands of Central Europe. Journal of Ecology 88: 424-435[CrossRef]
Suyama Y. K. Obayashi I. Hayashi 2000 Clonal structure in a dwarf bamboo (Sasa senanensis) population inferred from amplified fragment length polymorphism (AFLP) fingerprints. Molecular Ecology 9: 901-906[CrossRef][Medline]
Wagner R. H. 1964 The ecology of Uniola paniculata L. in the dune-strand habitat of North Carolina. Ecological Monographs 34: 79-96[CrossRef]
Wallace B. 1984 A possible explanation for observed differences in the geographical distributions of chromosomal arrangements of plants and Drosophila. Egyptian Journal of Genetics and Cytology 13: 121-136
Weeden N. F. J. F. Wendel 1989 Genetics of plant isozymes. In D. E. Soltis and P. S. Soltis [eds.], Isozymes in plant biology, 4672. Dioscorides Press, Portland, Oregon, USA
Wendel J. F. C. R. Parks 1982 Genetic control of isozyme variation in Camellia japonica. Journal of Heredity 73: 197-204
Widen B. N. Cronberg M. Widen 1994 Genotypic diversity, molecular markers, and spatial distribution of genets in clonal plants: a literature survey. Folia Geobotanica et Phytotaxonomica 29: 245-263
This article has been cited by other articles:
![]() |
E. PEREZ-COLLAZOS and P. CATALAN Palaeopolyploidy, Spatial Structure and Conservation Genetics of the Narrow Steppe Plant Vella pseudocytisus subsp. paui (Vellinae, Cruciferae) Ann. Bot., April 1, 2006; 97(4): 635 - 647. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. MANDAK, K. BIMOVA, I. PLACKOVA, V. MAHELKA, and J. CHRTEK Loss of Genetic Variation in Geographically Marginal Populations of Atriplex tatarica (Chenopodiaceae) Ann. Bot., October 1, 2005; 96(5): 901 - 912. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |