|
|
||||||||
1 Department of Biology, University of South Dakota, Vermillion, South Dakota 57069, USA
2 Department of Biological Sciences, St. Cloud State University, St. Cloud, Minnesota 56301, USA
3 Corresponding author (email: hbritten{at}usd.edu)
| ABSTRACT |
|---|
|
|
|---|
Key words: Black-tailed prairie dog, gene diversity, heterozygosity, plague, population genetic structure, RAPD-PCR, Yersinia pestis.
| INTRODUCTION |
|---|
|
|
|---|
Sylvatic plague, caused by the bacterium Yersinia pestis, was introduced to North America from Asia approximately 100 years ago. Plague is a generalist bacterium spread by fleas or direct contact and is now found in over 200 species of mammals (Pollitzer and Meyer, 1961; Poland and Barnes, 1979; Biggens and Kosoy, 2001; Parkhill et al., 2001). In North America, the significant and widespread declines of black-tailed prairie dogs (Cynomys ludovicianus), which prompted candidate listing under the Endangered Species Act, are attributed, in part, to the devastating effects of plague (Gober, 2000). Sylvatic plague is the only disease known to cause widespread fatalities of black-tailed prairie dogs, and no evidence has been found confirming resistance (Barnes, 1993; Cully and Williams, 2001). High-density communal living and clumped distributions of prairie dog colonies favor the spread of plague. Plague epizootics lead to substantial reductions in colony size (often 99%), local extirpations, and increased isolation among surviving colonies (Biggens and Kosoy, 2001; Cully and Williams, 2001; Lomolino and Smith, 2001).
Population bottlenecks caused by sylvatic plague are expected to decrease heterozygosity and increase the loss of alleles in affected prairie dog colonies (Hartl and Clark, 1997; Frankham et al., 2002). Moreover, recovering colonies, particularly those that are small and isolated, are subject to genetic drift and inbreeding, processes that contribute further to the loss of variability. However, gene flow via migration among recovering colonies can potentially restore genetic variability lost during bottlenecks. Sylvatic plague epizootics provide a unique opportunity to test the relative importance of genetic drift, inbreeding, and migration in prairie dogs.
Plague greatly reduced the range of black-tailed prairie dogs in Phillips County, Montana, USA, during the early 1990s (Fauna West Wildlife Consultants, 1999; Gober, 2000; Potera, 2000). Occupied range of prairie dogs decreased by 45%, from approximately 10,600 to 5,800 ha. The smallest occupied range was recorded in 1996, but since then the black-tailed prairie dog population has slowly increased. Our goal was to identify the effects of sylvatic plague and geographic isolation on genetic variability in prairie dog colonies. We hypothesized that colonies recovering from plague lost genetic variability by both the initial reduction in the number of individuals and subsequent drift and inbreeding. Furthermore, we expected to detect higher levels of genetic variability in clumped colonies as they should have received more gene flow, which would have ameliorated the effects of drift regardless of plague history. The combination of plague and isolation was expected to contribute to substantial decreases in genetic variability within colonies and genetic differentiation among colonies.
| MATERIALS AND METHODS |
|---|
|
|
|---|
We extracted DNA from 0.1 g of liver tissue using the Gen EluteTM Mammalian Genomic DNA Kit (Sigma, St. Louis, Missouri, USA). Randomly chosen oligonucleotide 10-mer primers (Operon Technologies, Alemeda, California, USA) were used in randomly amplified polymorphic DNA (RAPD) analysis (Palumbi, 1996) following initial screening confirming polymorphism. Ready-To-GoTM RAPD Analysis Beads (Amersham Pharmacia Biotech, Piscataway, New Jersey, USA) were used in PCR reactions following manufacturers protocols. To ensure consistent results, we used negative controls in all RAPD-PCR reactions and gel runs, and we reanalyzed samples for correspondence (Silva and Russo, 2000). We electrophoresed amplified samples on 1% agarose gels in 1 x TBE buffer for approximately 2 hr at 100 V using a Power-Pac 300 power supply (Bio-Rad, Hercules, California, USA). AmpliSizeTM Molecular Rulers were included as size standards (Bio-Rad). Gels were stained with ethidium bromide, visualized with ultraviolet light, and recorded digitally. Only primers exhibiting polymorphism upon screening were chosen for amplification, and all loci consistently reproduced were scored. We scored bands as 1 or 0 indicating band presence and absence, respectively. Loci were assumed to be biallelic and representative of only one chromosomal location.
Because RAPD markers are dominant, the presence of a band on a gel can represent either a homozygous or heterozygous genotype at any given locus. Therefore, it is necessary to estimate the frequencies of recessive alleles at any locus using the Hardy-Weinberg equilibrium equations and their attendant assumptions. A previous study using allozymes (codominant markers for which all heterozygotes can be directly identified) at four of the current study colonies (UL Bend, B01, B49, and B98) showed overall conformance to Hardy-Weinberg expectations (Britten, unpublished data), which suggested that this approach was appropriate. However, comparisons across studies using RAPDs and other markers that provide multilocus genotypes (e.g., isozymes and microsatellites) should be interpreted with caution because the same set of multilocus markers is rarely used in more than one study and because of the need to estimate frequencies of recessive alleles when using dominant markers (Sunnucks, 2000).
Number of alleles per locus, number of effective alleles per locus, allelic diversity, percent polymorphism, and Neis (1987) gene diversity (h, an analog of Hardy-Weinberg expected heterozygosity) were estimated for all six study colonies using POPGENE (version 1.31; Yeh et al., 1999). Mean per locus and per colony gene diversities were determined for all six colonies. We pooled per locus mean gene diversities for colonies that had experienced plague versus colonies without plague and clumped versus peripheral colonies. We used Wilcoxon paired-sample tests to determine differences in average per locus effective number of alleles and gene diversity between both colonies with and without plague and clumped and peripheral groups.
We used POPGENE to estimate GST, a measure of genetic differentiation between colonies. Arlequin (version 2.000; Schneider et al., 2000) was used to perform analyses of molecular variance (AMOVA) to examine population genetic structuring in sampled prairie dog colonies. Analysis of molecular variance is a hierarchical analysis of population substructure that partitions allele frequency variance components among all colonies and within colonies (Excoffier et al., 1992; Schneider et al., 2000). We then generated genetic differentiation statistics (F-statistics). Analysis of molecular variance and the F-statistic variants (GST and
-statistics) all estimate levels of variance in allele frequencies within vs. among colonies. Finally, we used POPGENE to calculate the number of genetically effective migrants among colonies per generation (Nm) from
ST and GST.
| RESULTS |
|---|
|
|
|---|
The percentage of polymorphic loci per colony ranged from 58.3% to 95.8% (Table 1
). Colonies that experienced plague (h=0.27±0.18) had significantly lower average per locus pooled gene diversity than colonies without plague (h=0.36±0.14, T=40, P<0.001, n=24; Table 1
). Similarly, colonies that experienced plague had a lower average effective number of alleles per locus (ne=1.44±0.35) than colonies without plague (ne=1.62±0.31; T=56, P<0.01, n=24; Table 1
). The difference in pooled gene diversity between clumped (h=0.29±0.17) and peripheral colonies (h=0.32±0.17) was not statistically significant (T=113, P>0.20, n=24; Table 1
). Likewise, peripheral colonies exhibited a slight but nonsignificantly higher average effective number of alleles per locus (ne=1.55±0.35) than clumped colonies (ne=1.50±0.35, T=117, P>0.05, n=24; Table 1
).
|
ST=0.232, P<0.0001; GST=0.194, P<0.0001). Numbers of genetically effective migrants per generation (Nm) were 0.873 and 1.039 estimated from
ST and GST, respectively. | DISCUSSION |
|---|
|
|
|---|
We predicted that degree of geographic isolation would result in a negative correlation with genetic variability. However, peripheral prairie dog colonies exhibited higher, yet statistically nonsignificant, average gene diversity than clumped colonies. There are two possible explanations for this result. First, designation of peripheral colonies in this study may not have captured the scale of isolation necessary to distinguish clumped from peripheral colonies because black-tailed prairie dogs migrate up to 10 km (Knowles, 1985; Garrett and Franklin, 1988; Table 1
). Thus, gene flow may account for the maintenance of variability in peripheral colonies. Second, Fennell (2002) found that colonies farther than 7 km from their nearest neighbor with plague persisted during an epizootic on the Northern Cheyenne Indian Reservation in central Montana, whereas colonies with closer neighbors with plague were unlikely to persist. The peripheral colony that did not experience plague had the highest genetic variability of colonies we sampled and was located approximately 50 km from the nearest colony with plague. Thus, geographic isolation appears to protect colonies from contracting plague, allowing them to remain large and to retain genetic variability.
A large proportion of total genetic differentiation (22.3%) was attributed to variation among colonies. Estimates of fixation indices support significant genetic structuring among the six sampled colonies, and our indices (
ST=0.223; GST=0.194) are higher than reported in the literature for other black-tailed prairie dog colonies (Chesser, 1983, FST=0.103; Daley, 1992, FST=0.115; Foltz and Hoogland, 1983, FST=0.028; Roach et al., 2001,
=0.118). Relatively high fixation indices observed for prairie dog colonies in Montana may represent lower levels of gene flow among colonies than observed in other studies, although comparisons across studies should be viewed with caution because of the variety of molecular markers (isozymes and microsatellites) used in the previous studies.
Number of genetically effective prairie dog migrants per generation ranged from 0.87 to 1.04 and are low enough that population genetic differentiation may be expected (Mills and Allendorf, 1996). Estimates of genetically effective numbers of migrants per generation are derived from F-statistics and their analogs as a means of comparing gene flow levels among populations, not as a means of comparing the actual number of dispersing individuals. Mills and Allendorf (1996) and Frankham et al. (2002) discuss the assumptions used to estimate gene flow. Intermediate levels of gene flow will not prevent divergence in allele frequencies among colonies, but may maintain levels of genetic variability within each disjunctive colony (Stenseth and Lidicker, 1992; Frankham et al., 2002). Although we observed a significant reduction in genetic diversity within colonies that experience plague, gene flow may ameliorate the effects of the sylvatic plague by reintroducing variability to genetically depauperate postplague colonies. Given enough time, gene flow will erase the effects of plague on genetic variability, assuming that colony size is stable or increases following recovery.
Finally, our results suggest that management of clumped black-tailed prairie dog colonies in north-central Montana should focus on preserving gene flow among colonies to maintain genetic variability within the system following plague-induced bottlenecks. This form of metapopulation structure would also facilitate demographic rescue of colonies that are affected by plague. Implementation of this strategy, however, must be balanced by the threat of increased mobility of plague-infected fleas during intercolony movement of prairie dogs. Peripheral colonies benefit from demographic and genetic input from dispersers. However, peripheral colonies may also act as relatively persistent reservoirs of genetic variability during plague epizootics because of their isolation. As such, care in managing peripheral colonies may be critical to the persistence of black-tailed prairie dogs in a given area during a plague epizootic.
| ACKNOWLEDGMENTS |
|---|
| LITERATURE CITED |
|---|
|
|
|---|
BIGGINS, D. E., AND M. Y. KOSOY. 2001. Influences of introduced plague on North American mammals: Implications from ecology of plague in Asia. Journal of Mammalogy 82: 906916.
CHESSER, R. K. 1983. Genetic variability within and among populations of the black-tailed prairie dog. Evolution 37: 320331.
CULLY, J. F., AND E. S. WILLIAMS. 2001. Interspecific comparisons of sylvatic plague in prairie dogs. Journal of Mammalogy 82: 894905.
DALEY, J. G. 1992. Population reductions and genetic variability in black-tailed prairie dogs. Journal of Wildlife Management 56: 212220.
EXCOFFIER, L., P. SMOUSE, AND J. M. QUATTRO. 1992. Analysis of molecular variance inferred from metric distances among DNA haplotypes: Application to human mitochondrial DNA restriction data. Genetics 131: 479491.[Abstract]
FAUNAWEST WILDLIFE CONSULTANTS. 1999. Status of the black-and white-tailed prairie dogs in Montana. Montana Fish, Wildlife, and Parks, Helena, Montana, pp. 128.
FENNELL, J. D. 2002. Recovery of black-tailed prairie dog colonies following a sylvatic plague epizootic. M.S. Thesis, Department of Biology, Montana State University, Bozeman, Montana, 29 pp.
FOLTZ, D. W., AND J. L. HOOGLAND. 1983. Genetic evidence of outbreeding in the black-tailed prairie dog (Cynomys ludovicianus). Evolution 37: 273281.
FRANKEL, O. H., AND M. E. SOULé . 1981. Conservation and evolution. Cambridge University Press, Cambridge, UK, 327 pp.
FRANKHAM, R., J. D. BALLOU, AND D. A. BRISCOE. 2002. Introduction to conservation genetics. Cambridge University Press, New York, New York, 613 pp.
GARRETT, M. G., AND W. L. FRANKLIN. 1988. Behavioral ecology of dispersal in the black-tailed prairie dog. Journal of Mammalogy 69: 236250.
GOBER, P. 2000.12-Month administrative finding, black-tailed prairie dog. Federal Register 65: 54765488.
HARTL, D. L., AND A. G. CLARK. 1997. Principles of population genetics, 3rd Edition. Sinauer Associates, Inc., Sunderland, Massachusetts, 542 pp.
KALINOWSKI, S. T., AND R. S. WAPLES. 2002. Relationship of effective to census size in fluctuating populations. Conservation Biology 16: 129136.
KNOWLES, C. J. 1985. Observations on prairie dog dispersal in Montana. Prairie Naturalist 17: 3339.
LOMOLINO, M. V., AND G. A. SMITH. 2001. Dynamic biogeography of prairie dog (Cynomys ludovicianus) towns near the edge of their range. Journal of Mammalogy 82: 937945.
MILLS, S. L., AND F. W. ALLENDORF. 1996. The one-migrant-per-generation rule in conservation and management. Conservation Biology 10: 15091518.
NEI, M. 1987. Molecular evolutionary genetics. Columbia University Press, New York, New York, 512 pp.
OBRIEN, S. J., AND J. F. EVERMANN. 1988. Interactive influence of infectious disease and genetic diversity in natural populations. Trends in Ecology and Evolution 3: 254259.
, J. S. MARTENSON, S. MITHTHAPALA, D. JANCZEWSKI, J. PECON-SLATERY, W. JOHNSON, D. A. GILBERT, M. ROELKE, C. PACKER, M. BUSH, AND D. E. WILDT. 1996. Conservation genetics of the felidae. In Conservation genetics: Case histories from nature, J. C. Avise and J. L. Hamrick (eds.). Chapman and Hall, New York, New York, pp. 5074.
PALUMBI, S. R. 1996. Nucleic acids II: The polymerase chain reaction. In Molecular Systematics, 2nd Edition, D. M. Hillis, C. Moritz, and B. K.Mable (eds.). Sinauer Associates, Sunderland, Massachusetts, pp. 204247.
PARKHILL, J., B. W. WREN, N. R. THOMSON, R. W.TITBALL, M. T. G. HOLDEN, M. B. PRENTICE, M. SEBAIHIA, K. D. JAMES, C. CHURCHER, K. L. MUNGALL, S. BAKER, D. BASHAM, S. D. BENTLEY, K. BROOKS, A. M. CERDENO-TáRRAGA, T. CHILLINGWORTH, A. CRONIN, R. M. DAVIES, P. DAVIS, G. DOUGAN, T. FELTWELL, N. HAMLIN, S. HOLROYD, K. JAGELS, A. V. KARLYSHEV, S. LEATHER, S. MOULE, P. C. F. OYSTON, M. QUAIL, K. RUTHERFORD, M. SIMMONDS, J. SKELTON, K. STEVENS, S. WHITEHEAD, AND B. G. BARRELL. 2001. Genome sequence of Yersinia pestis, the causative agent of plague. Nature 413: 523527.[Medline]
POLAND, J. D., AND A. M. BARNES. 1979. Plague. In CRC handbook series in zoonoses, H. Stoenner, W. Kaplan, and M. Torten (eds.). CRC Press, Inc., Boca Raton, Florida, pp. 515559.
POLLITZER, R., AND K. F. MEYER. 1961. The ecology of plague. In Studies in disease ecology, J. M. May (ed.). Hafner Publishing Company, Inc., New York, New York, pp. 433590.
POTERA, C. 2000. Prairie dogs plagued by YERSINIA pestis. ASM News 66: 718719.
ROACH, J. L., P. STAPP, B. VAN HORNE, AND M. F. ANTOLIN. 2001. Genetic structure of a meta-population of black-tailed prairie dogs. Journal of Mammalogy 82: 946959.
SCHNEIDER, S., D. ROESSLI, AND L. EXCOFFIER. 2000. Arlequin: A software for population genetics data analysis, Version 2.000. Genetics and Biometry Laboratory, Dept. of Anthropology, University of Geneva, http:lgb.unige.ch/arlequin/. Accessed 9 October 2003.
SEDDON, J. M., AND P. R. BAVERSTOCK. 1999. Variation on islands: Major histocompatibility complex (Mhc) polymorphism in population of the Australian bush rat. Molecular Ecology 8: 20712079.[Medline]
SILVA, E. P., AND C. A. M. RUSSO. 2000. Techniques and statistical data analysis in molecular population genetics. Hydrobiologia 420: 119135.
STENSETH, N. C., AND W. Z. LIDICKER. 1992. Animal dispersal: Small mammals as a model. Chapman and Hall, New York, New York, 365 pp.
SUNNUCKS, P. 2000. Efficient genetic markers for population biology. Trends in Ecology and Evolution. 15: 199203.
TRAVIS, S. E., C. N. SLOBODCHIKOFF, AND P. KEIM. 1997. DNA fingerprinting reveals low genetic diversity in Gunnisons prairie dog (Cynomys gunnisoni). Journal of Mammalogy 78: 725732.
VUCETICH, J. A., T. A. WAITE, AND L. NUNNEY. 1997. Fluctuating population size and the ratio of effective to census population size. Evolution 51: 20172021.
WESTEMEIER, R. L., J. D. BRAWN, S. A. SIMPSON, T.L. ESKER, R. W. JANSEN, J. W. WALK, E. L. KIR-SHNER, J. L. BOUZAT, AND K. N. PAIGE. 1998. Tracking the long-term decline and recovery of an isolated population. Science 282: 16951698.
YEH, F. C., R. C. YANG, T. J. B. BOYLE, Z. H. YE, AND J. X. MAO. 1999. POPGENE, the user-friendly shareware for population genetic analysis. Molecular Biology and Biotechnology Centre, University of Alberta, Canada, http://www.ualberta.ca/~fyeh/. Accessed 9 October 2003.
Received for publication 5 April 2003.
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |