Population genetics of Ice Age brown bears
Abstract
The Pleistocene was a dynamic period for Holarctic mammal species, complicated by episodes of glaciation, local extinctions, and intercontinental migration. The genetic consequences of these events are difficult to resolve from the study of present-day populations. To provide a direct view of population genetics in the late Pleistocene, we measured mitochondrial DNA sequence variation in seven permafrost-preserved brown bear (Ursus arctos) specimens, dated from 14,000 to 42,000 years ago. Approximately 36,000 years ago, the Beringian brown bear population had a higher genetic diversity than any extant North American population, but by 15,000 years ago genetic diversity appears similar to the modern day. The older, genetically diverse, Beringian population contained sequences from three clades now restricted to local regions within North America, indicating that current phylogeographic patterns may provide misleading data for evolutionary studies and conservation management. The late Pleistocene phylogeographic data also indicate possible colonization routes to areas south of the Cordilleran ice sheet.
The major climatic changes that occurred toward the end of the late Pleistocene had an important influence on the evolution and distribution of extant taxa. However, the genetic consequences of these events have been difficult to determine (1–6). The brown bear is a large, Holarctic carnivore whose distribution was dramatically altered by late Pleistocene events. In North America, the brown bear has had a limited history, appearing in eastern Beringia only 50–70,000 years ago and spreading into the contiguous United States about 13,000 years ago (7, 8). Previous research on mitochondrial control region sequences from 317 extant brown bears found that they defined one European (I) and three distinct North American (II, IIIa/IIIb, IV) clades (9) (Figs. (Figs.11 and and22a). Bears with sequences from the basal and divergent clade II were restricted to the Admiralty, Baranof, and Chichagof (ABC) islands. Sequences from populations in northern Canada and Alaska fell in a clade comprising two closely related groups (IIIa and IIIb), whereas populations in southern Canada and the contiguous United States belonged to a southern clade (IV). Of 22 localities surveyed, none had sequences from more than one clade, with one exception, where IIIa and IIIb cooccurred (ref. 9, Fig. Fig.22a). Sequences from the western Alaskan clade (IIIa) also were found in Asia and northern and eastern Europe, suggesting a recent connection across the Bering Land Bridge (9, 10).
Maximum likelihood trees of brown bear control region and cytochrome b sequences. For simplicity, we use the clade I–IV nomenclature that was established previously (9), although in our smaller dataset several clades appear paraphyletic, particularly in the less well resolved cytochrome b tree. Extant sequences are in color [W (9), TS (10), TB (11)], and permafrost sequences (P1–P7) are in black. Reliability values (above nodes) and bootstrap percentages (below nodes) with values greater than 50% are indicated. Corrected radiometric dates (12) and reference numbers are as follows: P1, UCR3742/CAMS-51806 (15,370 ± 60 bp); P2, UCR3741/CAMS-51805 (14,980 ± 60 bp); P3, UCR3743/CAMS-54128 (13,760 ± 50 bp); P4, UCR3745/CAMS-54129 (15,680 ± 50 bp); P5, UCR3746/CAMS-54130 (42,850 ± 850); P6, UCR3744/CAMS-51808 (35,970 ± 660); and P7, Beta16162 (36,500 ± 1,150).
Current (A) and past geographic distribution of brown bear control region sequence clades 15,000 (B) and 35,000–45,000 (C) years ago. Locations of modern samples (9) are indicated by triangles. The asterisk indicates a brown bear fossil dated to greater than 40,000 years ago on Prince of Wales Island (16), and no fossils are known in the contiguous United States before 13,000 years. The approximate extent of glacial ice sheets 15,000 ago is superimposed on current continental boundaries in hatching (17).
The recent appearance of brown bears in North America and their division into four genetically distinct populations with highly divergent mitochondrial control region sequences implied separate invasions of North America from long-isolated Old World populations (9, 10). Clade II was postulated to represent the earliest invasion, followed by clades IV and IIIa/IIIb. Furthermore, the precise correspondence between geographic and phylogenetic divisions suggests three evolutionary significant units for conservation (9, 13) (clades II, III, and IV). In contrast, nuclear microsatellite data in Alaskan brown bears do not support long-term genetic isolation of the three Alaskan clades, suggesting the possibility of sex-biased dispersal (14).
Mammalian remains have been recovered from vast permafrost deposits in central Alaska and northwestern Canada as a consequence of placer gold-mine operations (8). More than 100,000 skeletal elements are preserved in collections at the American Museum of Natural History (AMNH), New York, and the Canadian Museum of Nature (CMN), Ottawa. Mitochondrial and nuclear DNA sequences have been amplified from several individual permafrost specimens, and cold conditions are thought to favor DNA survival (15). Therefore, museum collections of permafrost-preserved bones potentially represent an extensive genetic record of mammal populations over the past 50,000 years. To directly record changes in genetic diversity through the late Pleistocene, we examined nine unassociated, permafrost-preserved brown bear bones from localities about 400 km apart, near Fairbanks, Alaska, and Sixtymile, Yukon Territory (Fig. (Fig.22a).
Acknowledgments
We thank the CMN and AMNH for permission to sample their collections. We are grateful to Charles Marshall and the Oxford University Museum of Natural History for laboratory facilities. C. R. Harington, R. Tedford, B. Van Valkenburgh, L. Waits, and R. Ward provided helpful comments and assistance. This research was supported by the National Science Foundation, Natural Environment Research Council, and Royal Society.
Abbreviations
| CMN | Canadian Museum of Nature |
| AMNH | American Museum of Natural History |
| LGM | last glacial maximum |
Note Added in Proof
Recent studies of permafrost remains have demonstrated that it is possible to amplify single copy nuclear sequences, including nuclear copies of mitochondrial genes (35). The latter is unlikely to be a factor in the current study as several different primer pairs and sequence regions produced consistent results.
Footnotes
This paper was submitted directly (Track II) to the PNAS office.
Data deposition: The sequences reported in this paper have been deposited in the GenBank database (accession nos. {"type":"entrez-nucleotide","attrs":{"text":"AF225566","term_id":"7107377","term_text":"AF225566"}}AF225566–{"type":"entrez-nucleotide","attrs":{"text":"AF225572","term_id":"7107385","term_text":"AF225572"}}AF225572).
Article published online before print: Proc. Natl. Acad. Sci. USA, 10.1073/pnas.040453097.
Article and publication date are at www.pnas.org/cgi/doi/10.1073/pnas.040453097
References
- 1. Bermingham E, Rohwer S, Freeman S, Wood C. Proc Natl Acad Sci USA. 1992;89:6624–6628.
- 2. Hewitt G M. Biol J Linn Soc. 1996;58:247–276.[PubMed]
- 3. Klicka J, Zink R M. Science. 1997;277:1666–1669.[PubMed]
- 4. Avise J C, Walker D. Proc R Soc London Ser B. 1998;265:457–463.
- 5. Klicka J, Zink R M. Proc R Soc London Ser B. 1999;266:695–700.[PubMed]
- 6. Wooding S, Ward R. Mol Biol Evol. 1997;14:1096–1105.[PubMed]
- 7. Kurtén B, Anderson E Pleistocene Mammals of North America. New York: Columbia Univ. Press; 1980. [PubMed][Google Scholar]
- 8. Guthrie R D Frozen Fauna of the Mammoth Steppe. Chicago: Univ. of Chicago Press; 1990. [PubMed][Google Scholar]
- 9. Waits L P, Talbot S L, Ward R H, Shields G F. Cons Biol. 1998;12:408–417.[PubMed]
- 10. Talbot S L, Shields G F. Mol Phylogenet Evol. 1996;5:477–494.[PubMed]
- 11. Taberlet P, Bouvet J. Proc Roy Soc London Ser B. 1994;255:195–200.[PubMed]
- 12. Stuiver M, Polach H A. Radiocarbon. 1977;19:355–363.[PubMed]
- 13. Moritz C. Trends Ecol Evol. 1994;9:373–375.[PubMed]
- 14. Paetkau D, Shields G F, Strobeck C. Mol Ecol. 1998;7:1283–1292.[PubMed]
- 15. Lindahl T. Nature (London) 1993;362:709–715.[PubMed]
- 16. Heaton T H, Talbot S L, Shields G F. Quat Res. 1996;46:186–192.[PubMed]
- 17. Dawson A G Ice Age Earth. London: Routledge; 1992. [PubMed][Google Scholar]
- 18. Hanni C, Laudet V, Stehelin D, Taberlet P. Proc Natl Acad Sci USA. 1994;91:12336–12340.
- 19. Krings M, Stone A, Schmitz R W, Krainitski H, Stoneking M, Pääbo S. Cell. 1997;90:19–30.[PubMed]
- 20. Taylor R E Radiocarbon Dating: An Archaeological Perspective. New York: Academic; 1987. [PubMed][Google Scholar]
- 21. Kirner D L, Taylor R E, Southon J R. Radiocarbon. 1995;37:697–704.[PubMed]
- 22. Swofford, D. L. (1998) paup4d65: Phylogenetic Analysis Using Parsimony and Other Methods (Sinauer, Sunderland, MA), Test Version.
- 23. Strimmer K, Von Haeseler A. Mol Biol Evol. 1996;13:964–969.[PubMed]
- 24. Tajima F. Genetics. 1993;135:599–607.
- 25. Pielou E C After the Ice Age: The Return of Life to Glaciated North America. Chicago: Univ. of Chicago Press; 1991. [PubMed][Google Scholar]
- 26. Anderson B G, Borns H W., Jr . The Ice Age World. Oslo: Scandinavian Univ. Press; 1997. [PubMed]
- 27. Avise J C, Neigel J E, Arnold J. J Mol Evol. 1984;20:99–105.[PubMed]
- 28. Hudson R R. Oxford Surv Evol Biol. 1990;7:1–44.[PubMed]
- 29. Templeton A R, Crandall K A, Sing C F. Genetics. 1995;132:619–633.
- 30. Masuda R, Murata K, Aiurzaniin A, Yoshida M C. Hereditas. 1998;128:277–280.[PubMed]
- 31. Matsuhashi T, Masuda R, Mano T, Yoshida M C. Mol Biol Evol. 1999;16:676–684.[PubMed]
- 32. Byun S A, Koop B F, Reimchen T E. Evolution. 1997;51:1647–1653.[PubMed]
- 33. Mann D H, Hamilton T D. Quat Sci Rev. 1995;14:449–471.[PubMed]
- 34. Harington C R. Can J Earth Sci. 1975;12:903–919.[PubMed]
- 35. Greenwood A D, Capelli C, Possnert G, Pääbo S. Mol Biol Evol. 1999;16:1466–1473.[PubMed]

