*Max-Planck-Institute for Evolutionary Anthropology, Leipzig, Germany;
Institute of Legal Medicine, Catholic University of S. Cuore, Rome;
Fish and Wildlife Resources, University of Idaho;
Institut für Ur- und Frühgeschichte, Eberhard-Karls-Universität Tübingen, Tübingen, Germany;
||Institute of Palaeontology, University of Vienna, Vienna, Austria;
¶Institute of Quaternary Paleontology and Geology of the Croatian Academy of Sciences and Arts, Zagreb, Croatia
![]() |
Abstract |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
![]() |
Introduction |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Recently, ancient DNA analyses of the remains of 10 cave bears have been published (Loreille et al. 2001
), including a 282-bp piece of the control region for six specimens. Although this study convincingly showed that brown bears and cave bears are reciprocally monophyletic, it did not address questions about the genetic diversity or population history of cave bears. Here, we report the determination and analysis of 285 bp of the mitochondrial control region from 12 cave bears originating from nine different caves.
![]() |
Materials and Methods |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
|
Sequence Analysis
Three overlapping amplicons were used to reconstruct approximately 285 bp of the mtDNA hypervariable region I. Consensus sequences were constructed for each sample and aligned by eye to a brown bear sequence (GenBank accession number X75866). All 12 novel cave bear sequences determined in this study have been deposited in GenBank, accession numbers AJ300166AJ300177. For the analyses, 53 brown bear mtDNA lineages (i.e., unique mtDNA sequences) were obtained from GenBank and our own data. Six additional sequences, representing four additional cave bear mtDNA lineages, were obtained from the literature (Loreille et al. 2001
). A pyrimidine stretch of variable length was removed from all sequences, resulting in an alignment of 260 bp for both brown and cave bears. Unrooted quartet puzzling trees were reconstructed for the 13 cave bear mtDNA lineages using the PUZZLE program. The rate heterogeneity parameter was estimated from the data set, assuming gamma distributed rates with eight variable rate categories. Tree reconstructions were carried out using the TN model (Tamura and Nei 1993
) for substitution. The maximum likelihood tree found for the cave bears was compared with all possible trees where the sequences from the small-sized cave bears form a monophyletic group by a maximum likelihood ratio test (Kishino and Hasegawa 1989
) as implemented in the PUZZLE program (Strimmer and von Haeseler 1996
). A neighbor-joining tree was reconstructed with Kimura-2-parametercorrected distances for all brown and cave bear lineages with two black bear (Ursus americanus, accession numbers AF012320 and U34265) sequences as an outgroup using the program MEGA2 (Kumar et al. 2000
). Mean pairwise sequence differences (MPSDs), both within and between populations, were calculated using the program avh-2 (A. von Haeseler, personal communication). To account for a possible sampling bias, both in the cave bear (from some locations sequences could be obtained from two samples, from others only from one sample) and in the brown bear (e.g., 65 samples from Alaska gave identical sequences, whereas only few sequences are available from Asian brown bears) data, we used mtDNA lineages for estimating the MPSDs. For calculating the time to the most recent common ancestor (MRCA), a Bayesian approach (Tavare et al. 1997
) was used. We tested the constant size model using Tajima's D, Fu and Li's D* and Fu and Li's F* as implemented in the software package DNASP (Rozas and Rozas 1999
).
![]() |
Results |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
|
A total of 538 clones from 87 amplifications was used to reconstruct a 282- to 286-bp sequence of the mitochondrial control region from the 12 cave bear specimens (see supplementary information on MBE website: www.molbiolevol.org). The length variation is caused by the variable length of a pyrimidine stretch of about 28 bp, which is present in all members of the family Ursidae, except in the brown and polar bears, where it has been partly deleted. DNA sequences for six additional specimens from Belgium, France, and Spain (see fig. 1
), for which approximately 275 bp overlap with the 12 cave bear DNA sequences determined here, were obtained from the literature (Loreille et al. 2001
).
Dating
A subsample of each of the 12 samples from which the complete DNA sequence was determined was dated using TAMS. The ages range from approximately 26,500 to about 49,000 years (see supplementary information). It should be noted that carbon dates earlier than about 40,000 years are unreliable because trace amounts of contaminating modern carbon may yield later dates. Such dates should therefore be seen as minimum estimates. The ages of the samples from Loreille et al. (2001)
are in the range of 25,00045,000 years, according to stratigraphy.
Phylogenetic Analyses
Among the 18 mtDNA sequences analyzed, four pairs of sequences are identical. These four pairs come from the caves Geißenklösterle (Germany), Ramesch Cave (Austria), Sclayn (Belgium), and Cova Linares (Spain). In addition, the DNA sequences from Nixloch and Potocka Zijalka are identical after removal of length variation. The 13 cave bear mtDNA lineages (i.e., distinct mtDNA sequences) were aligned to 53 different brown and polar bear mtDNA lineages from Europe, Asia, and North America, as well as to two American black bear mtDNA lineages. A neighbor-joining tree (fig. 2
) relating these mtDNA lineages is consistent with previous findings in that polar bears fall within the variation of brown bears (Talbot and Shields 1996
; Shields et al. 2000
) and that cave bears are monophyletic with respect to brown and polar bears (Loreille et al. 2001
). Moreover, the brown bear mtDNA lineages fall into three major clades (fig. 2
). One of these contains mtDNA lineages found in western Europe (denoted west). Another clade contains mtDNA lineages from polar and brown bears from the Admirality, Baranof, and Chicagof Islands off the coast of Alaska (denoted polar/ABC), whereas the last clade contains mtDNA lineage haplotypes from brown bears in Europe, Asia, and North America (denoted east). These three mtDNA clades may reflect the outcome of population isolation in Late Pleistocene refugia (Kohn et al. 1995
, Hewitt 2000
) because the polar/ABC clade has been found in Pleistocene bear bones from the mainland of Alaska (Leonard, Wayne, and Cooper 2000
).
|
Intraspecific mtDNA Diversity
To account for the uneven geographical sampling of brown and cave bears, we used mtDNA lineages (i.e., different mtDNA sequences) to calculate the MPSD. The MPSD within all brown and polar bears is 13.2 ± 6.1, whereas it is 7.6 ± 4.0 for cave bears. To further account for the fact that we have only 13 cave bear mtDNA lineages but 53 brown and polar bear mtDNA lineages, 10,000 sets of 13 mtDNA lineages were randomly drawn from the 53 brown and polar bear mtDNA lineages, and for each set the MPSD was calculated. Only 0.4% of the distribution falls below the value of 7.6 for the cave bear (fig. 3a
). Thus, Late Pleistocene cave bears show only about half of the species-wide mtDNA diversity of extant brown bears.
|
Time to the MRCA
To estimate the time to the MRCA for cave bears, all brown bears, as well as the three brown bear mtDNA clades, we used a Bayesian approach (Tavare et al. 1997
). Prior distributions were adjusted to cover all plausible parameter values, whereas posterior distributions were estimated by 106 simulations (table 1
). Specifically, the prior generation time was set to 1017 years (Allendorf and Servheen 1986
), the divergence time between cave and brown bears to 1.21.7 Myr (Kurten 1957
; Kurten 1976
; Loreille et al. 2001
), the pairwise distance between cave bears and brown bears to 2035 differences, and the population size to 7,000425,000. The median time to the MRCA is largest for all brown bears (890,000 years), followed by the cave bears (660,000 years). For the three brown bear mtDNA clades, the times are 420,000, 350,000, and 640,000 years for polar/ABC, west, and east, respectively. Thus, cave bears had a genetic history as long as that of eastern brown bears, despite their much smaller geographical range.
|
Phylogeographic Pattern
Although brown bear mtDNA sequences fall into three large clades, bears with mtDNA sequences from different clades are sometimes found in the same area, e.g., in Romania (Kohn et al. 1995
) or on the ABC Islands (unpublished data), indicating at least some migration of females. Moreover, although some geographical structure in the frequency of mtDNA lineages is found within these clades (Waits et al. 1998
), mtDNA lineages are sometimes found over the entire geographical range of a clade. For example, one sequence from the east mtDNA clade was found in Russia and Estonia, as well as in Alaska. In the case of the cave bears, identical sequences were only found for pairs of samples derived from the same caves, and in fact in the five cases where two specimens from the same cave were analyzed, identical mtDNA sequences were found in four cases. This suggests that the gene flow between caves was restricted. On the other hand, the number of differences between mtDNA sequences correlates neither with geographical distance between the caves from which the samples were excavated nor with the age of the samples. For example, mtDNA sequences from two caves located only 3 km apart (Geißenklösterle and Hohle Fels) differ by 10 substitutions, whereas the sequences from Sclayn (Belgium) and Cova Linares (Spain) differ by only two substitutions. Similarly, the mtDNA sequences from Geißenklösterle and Hohle Fels differ by 10 substitutions but only by 3,000 years in age, whereas one of the Vindija sequences differs from the Nixloch sequence by two substitutions, despite an age difference of 6,000 years. Thus, over longer time periods, mtDNA lineages did not remain geographically restricted.
![]() |
Discussion |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Another potential problem is the presence of insertions of mitochondrial DNA in the nuclear genome. Such insertions may be mistakenly amplified when primers for mtDNA sequences are used (e.g., Zischler 1995a
, 1995b
; Greenwood and Pääbo 1999
). However, a nuclear insertion would need to fulfill several unusual requirements to contribute to the DNA sequences determined here. Firstly, nuclear insertions are expected to evolve slowly, yet the sequences determined here are highly variable. Secondly, in spite of this high polymorphism, all the specimens analyzed would come from homozygous individuals because unique sequences were determined for each individual. Thirdly, because three different primer pairs were used, they would all need to preferentially amplify the potential insertion rather than the mtDNA sequences (Krings et al. 1997
), in spite of the latter being expected to occur in much higher copy numbers than a putative insertion.
Finally, DNA damage in the template DNA can potentially cause incorrect DNA sequences to be determined from ancient remains, especially if some form of damage occurred with a high frequency at particular nucleotide positions. However, given the precautions taken, such errors are not expected to influence the results to any appreciable extent (Hofreiter et al. 2001a
).
Considering the foregoing, we conclude that the mtDNA sequences determined reflect the endogenous mitochondrial DNA of cave bears that lived during the Upper Pleistocene.
Diachronical Studies of mtDNA Diversity
With mtDNA sequences available from 18 cave bear samples, it is now possible to address questions about cave bear genetic diversity and genetic history. However, in contrast to the situation where a modern population is analyzed, ancient DNA sequences stem from different time points in the past. Thus, mutations as well as changes in allele frequencies caused by genetic drift and selection may have affected the population over the time period from which the DNA sequences are sampled.
The effect of mutations can be expected to be relatively minor in this, as well as most other sets of ancient DNA sequences. The age difference of 25,000 years between the oldest and the youngest mtDNA lineages accounts for about 4% of the time to the MRCA for all cave bears (660,000 years). Because the expected time to coalescence for a sample size of n in terms of N generations is E(TN) = 2(1 - 1/n), a tree based on n = 13 lineages is expected to coalesce to a tree with n = 9 lineages within 25,000 years. The difference in the expected length of a tree based on 13 lineages and one based on 9 lineages is about 12%. Because this is based on the maximum age difference of all possible lineage comparisons, less than 12% of the observed differences are expected to be caused by mutations.
In contrast, shifts in allele frequencies can be rapid, especially in subdivided populations, e.g., when a population goes extinct and its habitat is recolonized from another population (see Avise 2000
, pp. 232235, for an overview). Therefore, it is impossible to say whether at the time of our youngest samples, all the genetic variation observed in the entire data set was still present among cave bears. However, because mutations are not expected to influence results appreciably, it can be assumed that at the time of our oldest sample, all the variation observed was present among cave bears. Thus, the estimates of cave bear genetic diversity presented may be thought of as referring to a time point about 49,000 years ago.
Intraspecific mtDNA Diversity
It is interesting to compare the distribution of genetic diversity in brown and cave bears because, although closely related and morphologically rather similar, the two species differ markedly in distribution patterns. Whereas brown bears until recently were rather homogeneously distributed in wooded areas in the entire northern part of Eurasia and North America, the distribution of the cave bear was always patchy (Musil 1980
, pp. 510) and restricted to Europe and western Asia.
In four out of the five cases where two cave bear samples from the same cave were analyzed, their mtDNA sequences were identical and in the remaining case (Vindija) quite similar to one another. This seems to suggest that gene flow between different populations of cave bears may have been limited. Limited female-mediated gene flow has also been suggested for brown bears (Waits et al. 1998
), although it takes place over longer time periods (Leonard, Wayne, and Cooper 2000
).
On the other hand, mtDNA sequences originating from Geißenklösterle (north of the Alps), Nixloch and Potocka (in the Alps), and Vindija (south of the Alps), are all closely related and separated from the next closest sequence (Grotte Merveilleuse) by six substitutions. After removal of length variation, the sequences from Nixloch and Potocka are identical, and the sequences from Sclayn (Belgium) and Linares (Spain) differ by only two substitutions, despite the large geographical distance between the two sites. By contrast, the mtDNA sequence from Hohle Fels differs from the Geißenklösterle sequence by 10 substitutions although the two sites are located only 3 km from each other. Thus, although the sharing of two DNA sequences within, but not between, localities indicates that gene flow between cave bear populations was low over shorter time periods, migration even over long distances must have taken place.
If cave bears were subdivided into many populations with limited female genetic exchange, different mtDNA sequences could have become fixed in different populations. Because rare sequences are more likely to be lost in panmictic than in subdivided populations, this could have allowed the overall genetic diversity in cave bears to be higher than that in brown bears, which may have been more panmictic. An alternative explanation for the larger genetic diversity relative to geographic area in cave bears could be a recent geographical expansion of brown bears carrying the eastern mtDNA clade. Because no evidence for a population expansion could be detected in any of the three brown bear mtDNA clades or in all brown bear mtDNA, this possibility does not seem likely. Thus, we suggest that cave bears maintained a high level of genetic diversity relative to their geographical range because they were separated in many small populations with limited gene flow between them.
Phylogenetic Relationship of Brown and Cave Bears
The cave bears sequenced to date (Loreille et al. 2001
and this study) come from a large part of the range of this species (fig. 1
). In addition, they include a specimen from the Conturines population, which is characterized by primitive morphological features and is regarded as an early branch within cave bears (Rabeder and Nagel 2000
). Thus, the mtDNA sequences available can be expected to represent well the intraspecific variation among cave bears. In tree analyses, all cave bears fall together and outside brown bears (fig. 2
). As noted by others (Loreille et al. 2001
), this is further supported by the observation that none of the cave bear mtDNA sequences carries the partial deletion of the pyrimidine stretch found in all brown and polar bears. Thus, the theory of Mazza and Rustioni (1994)
that cave bears fall within the variation of brown bears can be rejected.
Morphological Evolution
Some cave bear populations were morphologically very different from others. The most conspicuous of these are the small-sized forms that existed at high elevations in the Alps and that have been interpreted as an adaptation to life at high altitudes (Ehrenberg 1929
; Ehrenberg 1942
; Rabeder 1999
; Rabeder, Nagel, and Pacher 2000
, pp. 6567). However, it is worth noting that the small-sized cave bears are not monphyletic with respect to their mtDNA. It is also worth noting that the ancestor of the cave bear, U. deningeri, was small in size and that extremely large-sized cave bears appeared only after 150,000 BP (Kurten 1976
, Rabeder 1999
). Thus, small size could be an ancestral feature rather than an adaptation in cave bears. This suggestion is supported by the observation that the mtDNA sequences from the Sclayn and Nixloch cave bears, which are both among the largest cave bear populations known (Kurten, 1976
, Rabeder 1999
), are separated by the second deepest divergence in the mtDNA tree (fig. 2
). It is further supported by the observation that the mtDNA sequence from the small-sized population from Ramesch falls at the base of the branch leading to the Sclayn population. Thus, we suggest that the small size of the Alpine cave bears was an ancestral trait and that large-sized bears evolved at least twice independently. Further work is necessary to clarify whether some cave bear population indeed decreased in size as an adaptation to life at high altitudes.
![]() |
Acknowledgements |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
![]() |
Footnotes |
---|
Keywords: ancient DNA
cave bear
genetic diversity
MRCA
parallel evolution
phylogeography
Address for correspondence and reprints: Michael Hofreiter, Max-Planck-Institute for Evolutionary Anthropology, Inselstr. 22, D-04103 Leipzig, Germany. hofreiter{at}eva.mpg.de
![]() |
References |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Allendorf F. W., C. Servheen, 1986 Genetics and conservation of grizzly bears Trends Ecol. Evol 1:88-89[ISI]
Avise J. C., 2000 Phylogeography. The history and formation of species Harvard University Press, Cambridge, Mass
Ehrenberg K., 1929 Die Ergebnisse der Ausgrabungen in der Schreiberwandhöhle am Dachstein Palaeont. Z. Berlin 11:261-268
. 1942 Berichte über Ausgrabungen in der Salzofenhöhle im Toten Gebirge. II. Untersuchungen über umfassendere Skelettfunde als Beitrag zur Frage der Form- und Größenverschiedenheiten zwischen Braunbär und Höhlenbär Palaeobiologica 7:531-666
Fu X. Y., W. H. Li, 1993 Statistical tests of neutrality of mutations Genetics 133:693-709
Greenwood A. D., S. Pääbo, 1999 Nuclear insertion sequences of mitochondrial DNA predominate in hair but not in blood of elephants Mol. Ecol 8:133-137[ISI][Medline]
Handt O., M. Krings, R. H. Ward, S. Pääbo, 1996 The retrieval of ancient human DNA sequences Am. J. Hum. Genet 59:368-376[ISI][Medline]
Hänni C., V. Laudet, D. Stehelin, P. Taberlet, 1994 Tracking the origins of the cave bear (Ursus spelaeus) by mitochondrial DNA sequencing Proc. Natl. Acad. Sci. USA 91:12336-12340
Hewitt G., 2000 The genetic legacy of the Quaternary ice ages Nature 405:907-913[ISI][Medline]
Höss M., P. Jaruga, T. H. Zastawny, M. Dizdaroglu, S. Pääbo, 1996 DNA damage and DNA sequence retrieval from ancient tissues Nucleic Acids Res 24:1304-1307
Höss M., S. Pääbo, 1993 DNA extraction from Pleistocene bones by a silica-based purification method Nucleic Acids Res 21:3913-3914[ISI][Medline]
Hofreiter M., V. Jaenicke, D. Serre, A. von Haeseler, S. Pääbo, 2001a. DNA sequences from multiple amplifications reveal artifacts induced by cytosine deamination in ancient DNA Nucleic Acids Res 29:4793-4799
Hofreiter M., D. Serre, H. N. Poinar, M. Kuch, S. Pääbo, 2001b. Ancient DNA Nat. Rev. Genet 2:353-359[ISI][Medline]
Kilger C., M. Krings, H. Poinar, S. Pääbo, 1997 "Colony Sequencing": direct sequencing of plasmid DNA from bacterial colonies BioTechniques 22:412-418[ISI][Medline]
Kishino H., M. Hasegawa, 1989 Evaluation of the maximum likelihood estimate of the evolutionary tree topologies from DNA sequence data, and the branching order in hominoidea J. Mol. Evol 29:170-179[ISI][Medline]
Kohn M., F. Knauer, A. Stoffella, W. Schröder, S. Pääbo, 1995 Conservation genetics of the European brown beara study using excremental PCR of nuclear and mitochondrial sequences Mol. Ecol 4:95-103[ISI][Medline]
Krings M., H. Geisert, R. W. Schmitz, H. Krainitzki, S. Pääbo, 1999 DNA sequence of the mitochondrial hypervariable region II from the Neandertal type specimen Proc. Natl. Acad. Sci. USA 96:5581-5585
Krings M., A. Stone, R. W. Schmitz, H. Krainitzki, M. Stoneking, S. Pääbo, 1997 Neanderthal DNA sequences and the origin of modern humans Cell 90:19-30[ISI][Medline]
Kumar S., K. Tamura, I. B. Jakobson, M. Nei, 2001 MEGA2: Molecular Evolutionary Genetic Analysis software Bioinformatics 17:12441245.
Kurten B., 1957 The bears and hyenas of the interglacials Quaternaria 4:69-81
. 1976 The cave bear story Columbia University Press, New York
Leonard J. A., R. K. Wayne, A. Cooper, 2000 Population genetics of Ice Age brown bears Proc. Natl. Acad. Sci. USA 97:1651-1654
Loreille O., L. Orlando, M. Patou-Mathis, M. Philippe, P. Taberlet, C. Hänni, 2001 Ancient DNA analysis reveals divergence of the cave bear, Ursus spelaeus, and brown bear, Ursus arctos, lineages Curr. Biol 11:200-203[ISI][Medline]
Mazza P., M. Rustioni, 1994 On the phylogeny of Eurasian bears Paleontogr. Abt. A 230:1-38
Musil R., 1980 Ursus spelaeusDer Höhlenbär Museum für Ur- und Frühgeschichte, Thüringens, Weimar, Germany
Poinar H. N., M. Hofreiter, W. G. Spaulding, P. S. Martin, B. A. Stankiewicz, H. Bland, R. P. Evershed, G. Possnert, S. Pääbo, 1998 Molecular coproscopy: dung and diet of the extinct ground sloth Nothotheriops shastensis Science 281:402-406
Rabeder G., 1983 Neues vom Höhlenbären: Zur Morphogenetik der Backenzähne Die Höhle 34:67-85
. 1999 Die Evolution des Höhlenbärgebisses Mitt. Komm. Quartärforsch 11:1-102
Rabeder G., D. Nagel, 2000 Alpine cave bears; phylogeny and a new taxonomy of Sixth International Cave Bear Symposium, Corona, Spain.
Rabeder G., D. Nagel, M. Pacher, 2000 Der Höhlenbär Jan Thorbecke Verlag, Stuttgart, Germany
Reisinger C., J. Hohenegger, 1998 Sexual dimorphism in limb bones of Late Pleistocene cave bear (Ursus spelaeus, Carnivora, Mammalia) from three caves in Eastern Alps (Austria and Italy) Boll. Soc. Paleontol. Ital 37:99-116
Rozas J., R. Rozas, 1999 DnaSP version 3: an integrated program for molecular population genetics and molecular evolution analysis Bioinformatics 15:174-175
Shields G. F., D. Adams, G. Garner, M. Labelle, J. Pietsch, M. Ramsay, C. Schwartz, K. Titus, S. Williamson, 2000 Phylogeography of mitochondrial DNA variation in brown bears and polar bears Mol. Phylogenet. Evol 15:319-326[ISI][Medline]
Strimmer K., A. von Haeseler, 1996 Quartet puzzling: a quartet maximum-likelihood method for reconstructing tree topologies Mol. Biol. Evol 13:964-969
Tajima F., 1989 Statistical method for testing the neutral mutation hypothesis by DNA polymorphism Genetics 123:585-595
Talbot S. L., G. F. Shields, 1996 Phylogeography of brown bears (Ursus arctos) of Alaska and paraphyly within the Ursidae Mol. Phylogenet. Evol 5:477-494[ISI][Medline]
Tamura K., M. Nei, 1993 Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees Mol. Biol. Evol 10:512-526[Abstract]
Tavare S., D. J. Balding, R. C. Griffiths, P. Donnelly, 1997 Inferring coalescence times from DNA sequence data Genetics 145:505-518
Waits L. P., S. L. Talbot, R. H. Ward, G. F. Shields, 1998 Mitochondrial DNA phylogeography of the North American brown bear and implications for conservation Conserv. Biol 12:408-417[ISI]
Zischler H., H. Geisert, A. von Haeseler, S. Pääbo, 1995a. A nuclear 'fossil' of the mitochondrial D-loop and the origin of modern humans Nature 378:489-492[ISI][Medline]
Zischler H., M. Höss, O. Handt, A. von Haeseler, A. C. van der Kuyl, J. Goudsmit, 1995b. Detecting dinosaur DNA Science 268:1192-1193[ISI][Medline]