Ancient DNA analysis reveals divergence of the cave bear, Ursus spelaeus, and brown bear, Ursus arctos, lineages

The cave bear, Ursus spelaeus, represents one of the most frequently found paleontological remains from the Pleistocene in Europe. The species has always been confined to Europe and was contemporary with the brown bear, Ursus arctos. Relationships between the cave bear and the two lineages of brown bears defined in Europe, as well as the origins of the two species, remain controversial, mainly due to the wide morphological diversity of the fossil remains, which makes interpretation difficult [1, 2]. Sequence analysis of ancient DNA is a useful tool for resolving such problems because it provides an independent source of data [3]. We previously amplified a short DNA fragment of the mitochondrial DNA control region (mt control region) of a 40,000-year-old Ursus spelaeus sample [4]. In this paper, we describe the DNA analysis of two mtDNA regions, the control region and the cytochrome b gene. Control region sequences were obtained from ten samples of cave bears ranging from 130,000 to 20,000 years BP, and one particularly well-conserved sample gave a complete cyt b sequence. Our data demonstrate that cave bears split largely before the lineages of brown bears around 1.2 million years ago. Given its abundance, its wide distribution in space and time, and its large morphological diversity, the cave bear is a promising model for direct observation of the evolution of sequences throughout time, extinction periods, and the differentiation of populations shaped by climatic fluctuations during the Pleistocene.     

A three-dimensional analysis of tooth-root morphology in living bears and implications for feeding behaviour in the extinct cave bear

ABSTRACT

The morphology of both crowns and tooth-roots reflects dietary specialisation in mammalian carnivores. In this article, we analyse the tooth-root morphology of maxillary teeth from CT scans of living bears (Ursus arctos, Ursus americanus, Ursus maritimus, Ursus thibetanus, Melursus ursinus, Helarctos malayanus, Tremarctos ornatus and Ailuropoda melanoleuca) in order to make inferences about the diet and feeding behaviour of the extinct cave bear (Ursus spelaeus sensu lato). Specifically, we investigate two major mitochondrial clades of extinct cave bears recognized by previous authors: Ursus ingressus and Ursus spelaeus (U. spelaeus spelaeus, U. spelaeus ladinicus, U. spelaeus eremus). Our results indicate a close association between tooth-root surface area and feeding behaviour in all living bear species. Tooth-root surface area values of cave bears suggest that they relied more on vegetative matter than living brown bears (Ursus arctos) but subtle differences between these species/subspecies could also indicate different feeding strategies among the members of cave bear complex.     

Ancestral Polymorphisms and Sex-Biased Migration Shaped the Demographic History of Brown Bears and Polar Bears

Recent studies have reported discordant gene trees in the evolution of brown bears and polar bears. Genealogical histories are different among independent nuclear loci and between biparentally inherited autosomal DNA (aDNA) and matrilineal mitochondrial DNA (mtDNA). Based on multi-locus genomic sequences from aDNA and mtDNA, we inferred the population demography of brown and polar bears and found that brown bears have 6 times (aDNA) or more than 14 times (mtDNA) larger population sizes than polar bears and that polar bear lineage is derived from within brown bear diversity. In brown bears, the effective population size ratio of mtDNA to aDNA was at least 0.62, which deviated from the expected value of 0.25, suggesting matriarchal population due to female philopatry and male-biased migration. These results emphasize that ancestral polymorphisms and sex-biased migration may have contributed to conflicting branching patterns in brown and polar bears across aDNA genes and mtDNA.     

Sequencing mtDNA of the cave bearUrsus spelaeus from the Bavarian Alps is feasible by nested and touchdown PCR

Ancient DNA from bones of the extinctUrsus spelaeus Rosamueller et Hainroth, 1794 found in the Bavarian Alps has been amplified by PCR. Two out of five samples yielded a distinct band of 135 bp originating from the mtDNA control region. A combination of nested and touchdown PCR supported the amplification. Analysis of the nucleotide sequences revealed four transitions compared to the French cave bear sequence, the only cave bear data known so far (Genbank database X80259, AF121779). The consensus distant matrix tree clustered the two cave bears next to the brown bearUrsus arctos Linnaeus, 1758.     

Ancient DNA evidence for the loss of a highly divergent brown bear clade during historical times

The genetic diversity of present-day brown bears (Ursus arctos) has been extensively studied over the years and appears to be geographically structured into five main clades. The question of the past diversity of the species has been recently addressed by ancient DNA studies that concluded to a relative genetic stability over the last 35,000 years. However, the post-last glacial maximum genetic diversity of the species still remains poorly documented, notably in the Old World. Here, we analyse Atlas brown bears, which became extinct during the Holocene period. A divergent brown bear mitochondrial DNA lineage not present in any of the previously studied modern or ancient bear samples was uncovered, suggesting that the diversity of U. arctos was larger in the past than it is now. Specifically, a significant portion (with respect to sequence divergence) of the intraspecific diversity of the brown bear was lost with the extinction of the Atlas brown bear after the Pleistocene/Holocene transition.     

Genetics and conservation of European brown bears Ursus arctos

• 1 We review the genetics research that has been conducted on the European brown bear Ursus arctos, one of the genetically best‐studied mammalian species.
• 2 The first genetics studies on European brown bears were on phylogeography, as a basis for proposed population augmentations. Two major mitochondrial DNA lineages, western and eastern, and two clades within the western lineage were found. This led to a hypothesis that brown bears had contracted to southern refugia during the last glacial maximum. More recent results suggest that gene flow among brown bears blurred this structure and they survived north of these putative refugia. Thus, today's structure might be a result of population fragmentation caused by humans.
• 3 The nuclear diversity of European brown bears is similar in range to that in North American bears: low levels occur in the small populations and high levels in the large populations.
• 4 Many non‐invasive genetic methods, developed during research on brown bears, have been used for individual identification, censusing populations, monitoring migration and gene flow, and testing methods that are easier to use in endangered populations and over large areas.
• 5 Genetics has been used to study many behavioural and population ecological questions that have relevance for the conservation and management of brown bears.
• 6 The European brown bear has served, and will continue to serve, as a model for the development of methods, analyses and hypotheses in conservation genetics.

     

Hominid exploitation of the environment and cave bear populations. The case of Ursus spelaeus Rosenmüller-Heinroth in Amutxate cave (Aralar, Navarra-Spain)

Cave bears (Ursus deningeri and U. spelaeus) and hominids (Homo heidelbergensis, H. neanderthalensis, and H. sapiens) were potential competitors for environmental resources (subterranean and open air). Here, we examined the age at death of cave bear (Ursus spelaeus Rosenmüller-Heinroth) specimens from Amutxate cave in order to shed light on the effect of resource sharing between cave bears and hominids. After studying dental wear of the deciduous and permanent dentitions, the ontogenetic development of mandibles, and incremental layers of cement (annuli), we defined five age groups differentiated by marked development and size gaps. Our findings indicate that after hibernating, bears abandoned the den, thereby leaving the subterranean environment (caves) free for temporary hominid occupation-this would explain the subtle traces of hominid presence in many dens. However, a simple calculation based on age at death of subadult and adult cave bear specimens in Amutxate cave, extrapolated to the whole cave area, showed that the area surrounding this cave hosted bears for at least 9,000 years. This length of habitation, quite similar to the time-span derived from amino acid racemization and electron spin resonance, indicates that bear populations in the Amutxate cave constituted a serious constraint for hominid exploitation of the environment.     

Genetic diversity of endangered brown bear (Ursus arctos) populations at the crossroads of Europe, Asia and Africa

Aim  Middle East brown bears ( Ursus arctos syriacus Hemprich and Ehrenberg, 1828) are presently on the edge of extinction. However, little is known of their genetic diversity. This study investigates that question as well as that of Middle East brown bear relationships to surrounding populations of the species.
Location  Middle East region of south-western Asia.
Methods  We performed DNA analyses on 27 brown bear individuals. Twenty ancient bone samples (Late Pleistocene to 20th century) from natural populations and seven present-day samples obtained from captive individuals were analysed.
Results  Phylogenetic analyses of the mitochondrial sequences obtained from seven ancient specimens identify three distinct maternal clades, all unrelated to one recently described from North Africa. Brown bears from Iran exhibit striking diversity (three individuals, three haplotypes) and form a unique clade that cannot be linked to any extant one. Individuals from Syria belong to the Holarctic clade now observed in Eastern Europe, Turkey, Japan and North America. Specimens from Lebanon surprisingly appear as tightly linked to the clade of brown bears now in Western Europe. Moreover, we show that U. a. syriacus in captivity still harbour haplotypes closely linked to those found in ancient individuals.
Main conclusion  This study brings important new information on the genetic diversity of brown bear populations at the crossroads of Europe, Asia and Africa. It reveals a high level of diversity in Middle East brown bears and extends the historical distribution of the Western European clade to the East. Our analyses also suggest the value of a specific breeding programme for captive populations.     

First DNA sequences from Asian cave bear fossils reveal deep divergences and complex phylogeographic patterns

Until recently, cave bears were believed to have only inhabited Europe. However, recent morphological evidence suggests that cave bears' geographic range extended as far east as Transbaikalia, Eastern Siberia. These Asian cave bears were morphologically distinct from European cave bears. However, how they related to European lineages remains unclear, stressing the need to assess the phylogenetic and phylogeographic relationship between Asian cave bears and their European relatives. In this work, we address this issue using a 227 base-pair fragment of the mitochondrial control region obtained from nine fossil bone samples from eight sites from the Urals, Caucasus, Altai Mountains, Ukraine and Yana River region in Eastern Siberia. Results of the phylogenetic analyses indicate that (i) the cave bear from the Yana River is most closely related to cave bears from the Caucasus region; (ii) the Caucasus/Yana group of bears is genetically very distinct from both European cave bears and brown bears, suggesting that these bears could represent an independent species; and (iii) the Western European cave bear lineage reached at least temporarily to the Altai Mountains, 7000 km east of their known centre of distribution. These results suggest that the diversity of cave bears was greater than previously believed, and that they could survive in a much wider range of ecological conditions than previously assumed. They also agree with recent studies on other extinct and extant species, such as wolves, hyenas and steppe bison, which have also revealed higher genetic and ecological diversity in Pleistocene populations than previously known.     

Recent history of the brown bear in the Maghreb

Bones of the brown bear (Ursus arctos, mammalia, Carnivora) found in a cave of the Akouker massif (Djurdjura, Algeria) have been dated according to the 14C method as belonging to the historical times (420-600 A.D.). The bone and teeth measurements correspond to a small-sized animal, the smallest ever found in the Maghreb. A review of fossil bears in the Quaternary faunas of North Africa clearly shows that the area of distribution, which was initially wide, shrank at the end of the Upper Pleistocene. The bears had temporarily taken to mountainous areas difficult to access. The bone remains discovered up to now prove that the brown bear was represented by individuals or populations showing a large diversity of size.     

Nuclear DNA microsatellite analysis of genetic diversity and gene flow in the Scandinavian brown bear (Ursus arctos)

In the 1930s, the Scandinavian brown bear was close to extinction due to vigorous extermination programmes in Norway and Sweden. Increased protection of the brown bear in Scandinavia has resulted in the recovery of four subpopulations, which currently contain close to 1000 individuals. Effective conservation and management of the Scandinavian brown bear requires knowledge of the current levels of genetic diversity and gene flow among the four subpopulations. Earlier studies of mitochondrial DNA (mtDNA) diversity revealed extremely low levels of genetic variation, and population structure that grouped the three northern subpopulations in one genetic clade and the southernmost subpopulation in a second highly divergent clade. In this study, we extended the analysis of genetic diversity and gene flow in the Scandinavian brown bear using data from 19 nuclear DNA microsatellite loci. Results from the nuclear loci were strikingly different than the mtDNA results. Genetic diversity levels in the four subpopulations were equivalent to diversity levels in nonbottlenecked populations from North America, and significantly higher than levels in other bottlenecked and isolated brown bear populations. Gene flow levels between subpopulations ranged from low to moderate and were correlated with geographical distance. The substantial difference in results obtained using mtDNA and nuclear DNA markers stresses the importance of collecting data from both types of genetic markers before interpreting data and making recommendations for the conservation and management of natural populations. Based on the results from the mtDNA and nuclear DNA data sets, we propose one evolutionarily significant unit and four management units for the brown bear in Scandinavia.     

Diet reconstruction in cave bears from craniodental morphology: past evidences,new results and future directions

ABSTRACT

The diet of the cave bear (Ursus spelaeus) is a controversial topic, as different paleobiological approaches (e.g. dental wear, isotopic biochemistry, skull morphometrics) result in different dietary inferences for the cave bear, ranging from carnivory to pure herbivory. Here, we review the main results obtained from these approaches, with special emphasis on those obtained from the morphometric analyses of the cave bear craniodental skeleton. Then, we compute a between-group Principal Components Analysis from a set of 3D-landmarks digitized on 103 mandibles of living bears and extinct cave bears and using a phylomorphospace approach. Moreover, we also reconstructed the evolutionary trajectory of the cave bear mandible from the hypothetical shape of its inferred ancestor. Our results indicate that the mandible of the cave bear possess specific traits indicative of a highly-herbivorous diet or, at least, more herbivorous than their closest living relative, the brown bear (Ursus arctos). Moreover, we also propose new directions for future research to obtain more detailed inferences on the potential food resources consumed by the cave bear being crucial to understand the ‘life and death’ of this vanished animal.     

Mitogenetic structure of brown bears (Ursus arctos L.) in northeastern Europe and a new time frame for the formation of European brown bear lineages

We estimated the phylogenetic relationships of brown bear maternal haplotypes from countries of northeastern Europe (Estonia, Finland and European Russia), using sequences of mitochondrial DNA (mtDNA) control region of 231 bears. Twenty-five mtDNA haplotypes were identified. The brown bear population in northeastern Europe can be divided into three haplogroups: one with bears from all three countries, one with bears from Finland and Russia, and the third composed almost exclusively of bears from European Russia. Four haplotypes from Finland and European Russia matched exactly with haplotypes from Slovakia, suggesting the significance of the current territory of Slovakia in ancient demographic processes of brown bears. Based on the results of this study and those from the recent literature, we hypothesize that the West Carpathian Mountains have served either as one of the northernmost refuge areas or as an important movement corridor for brown bears of the Eastern lineage towards northern Europe during or after the last ice age. Bayesian analyses were performed to investigate the temporal framework of brown bear lineages in Europe. The molecular clock was calibrated using Beringian brown bear sequences derived from radiocarbon-dated ancient samples, and the estimated mutation rate was 29.8% (13.3%-47.6%) per million years. The whole European population and Western and Eastern lineages formed about 175,000, 70,000 and 25,000 years before present, respectively. Our approach to estimating the time frame of brown bear evolution demonstrates the importance of using an appropriate mutation rate, and this has implications for other studies of Pleistocene populations.     

Conservation genetics of the European brown bear - a study using excremental PCR of nuclear and mitochondrial sequences

In the Brenta area of northern Italy, a brown bear Ursus arctos population is rapidly going extinct. Restocking of the population is planned. In order to study the genetics of this highly vulnerable population with a minimum of stress to the animals we have developed a PCR-based method that allows the study of mitochondrial and nuclear gene sequences from droppings collected in the field. This method is generally applicable to animals in the wild. Using excremental as well as hair samples, we show that the Brenta population is monomorphic for one mitochondrial lineage and that female as well as male bears exist in the area. In addition, 70 samples from other parts of Europe were studied. As others have previously reported, the mitochondrial gene pool of European bears is divided into two major clades, one with a western and the other with an eastern distribution. Whereas populations generally belong to either one or the other mitochondrial clade, the Romanian population contains both clades. The bears in the Brenta belong to the western clade. The implications for the management of brown bears in the Brenta and elsewhere in Europe are discussed.     

Genomic Evidence for Island Population Conversion Resolves Conflicting Theories of Polar Bear Evolution

Despite extensive genetic analysis, the evolutionary relationship between polar bears (Ursus maritimus) and brown bears (U. arctos) remains unclear. The two most recent comprehensive reports indicate a recent divergence with little subsequent admixture or a much more ancient divergence followed by extensive admixture. At the center of this controversy are the Alaskan ABC Islands brown bears that show evidence of shared ancestry with polar bears. We present an analysis of genome-wide sequence data for seven polar bears, one ABC Islands brown bear, one mainland Alaskan brown bear, and a black bear (U. americanus), plus recently published datasets from other bears. Surprisingly, we find clear evidence for gene flow from polar bears into ABC Islands brown bears but no evidence of gene flow from brown bears into polar bears. Importantly, while polar bears contributed <1% of the autosomal genome of the ABC Islands brown bear, they contributed 6.5% of the X chromosome. The magnitude of sex-biased polar bear ancestry and the clear direction of gene flow suggest a model wherein the enigmatic ABC Island brown bears are the descendants of a polar bear population that was gradually converted into brown bears via male-dominated brown bear admixture. We present a model that reconciles heretofore conflicting genetic observations. We posit that the enigmatic ABC Islands brown bears derive from a population of polar bears likely stranded by the receding ice at the end of the last glacial period. Since then, male brown bear migration onto the island has gradually converted these bears into an admixed population whose phenotype and genotype are principally brown bear, except at mtDNA and X-linked loci. This process of genome erosion and conversion may be a common outcome when climate change or other forces cause a population to become isolated and then overrun by species with which it can hybridize.     

Ancient DNA reveals differences in behaviour and sociality between brown bears and extinct cave bears

Ancient DNA studies have revolutionized the study of extinct species and populations, providing insights on phylogeny, phylogeography, admixture and demographic history. However, inferences on behaviour and sociality have been far less frequent. Here, we investigate the complete mitochondrial genomes of extinct Late Pleistocene cave bears and middle Holocene brown bears that each inhabited multiple geographically proximate caves in northern Spain. In cave bears, we find that, although most caves were occupied simultaneously, each cave almost exclusively contains a unique lineage of closely related haplotypes. This remarkable pattern suggests extreme fidelity to their birth site in cave bears, best described as homing behaviour, and that cave bears formed stable maternal social groups at least for hibernation. In contrast, brown bears do not show any strong association of mitochondrial lineage and cave, suggesting that these two closely related species differed in aspects of their behaviour and sociality. This difference is likely to have contributed to cave bear extinction, which occurred at a time in which competition for caves between bears and humans was likely intense and the ability to rapidly colonize new hibernation sites would have been crucial for the survival of a species so dependent on caves for hibernation as cave bears. Our study demonstrates the potential of ancient DNA to uncover patterns of behaviour and sociality in ancient species and populations, even those that went extinct many tens of thousands of years ago.     

Sudden expansion of a single brown bear maternal lineage across northern continental Eurasia after the last ice age: a general demographic model for mammals?

The brown bear has proved a useful model for studying Late Quaternary mammalian phylogeography. However, information is lacking from northern continental Eurasia, which constitutes a large part of the species' current distribution. We analysed mitochondrial DNA sequences (totalling 1943 bp) from 205 bears from northeast Europe and Russia in order to characterize the maternal phylogeography of bears in this region. We also estimated the formation times of the sampled brown bear lineages and those of its extinct relative, the cave bear.
Four closely related haplogroups belonging to a single mitochondrial subclade were identified in northern continental Eurasia. Several haplotypes were found throughout the whole study area, while one haplogroup was restricted to Kamchatka. The haplotype network, estimated divergence times and various statistical tests indicated that bears in northern continental Eurasia recently underwent a sudden expansion, preceded by a severe bottleneck. This brown bear population was therefore most likely founded by a small number of bears that were restricted to a single refuge area during the last glacial maximum. This pattern has been described previously for other mammal species and as such may represent one general model for the phylogeography of Eurasian mammals. Bayesian divergence time estimates are presented for different brown and cave bear clades. Moreover, our results demonstrate the extent of substitution rate variation occurring throughout the phylogenetic tree, highlighting the need for appropriate calibration when estimating divergence times.     

Phylogenetic Relationships Among Fragmented Asian Black Bear (Ursus Thibetanus) Populations in Western Japan

The number of Asian black bears (Ursus thibetanus) in Japan has been reduced and their habitats fragmented and isolated because of human activities. Our previous study examining microsatellite DNA loci revealed significant genetic differentiation among four local populations in the western part of Honshu. Here, an approximate 700-bp nucleotide sequence of mitochondrial DNA (mtDNA) control region was analysed in 119 bears to infer the evolutionary history of these populations. Thirteen variable sites and variation in the number of Ts at a T-repeat site were observed among the analysed sequences, which defined 20 mtDNA haplotypes with the average sequence divergence of 0.0051 (SD = 0.00001). The observed haplotype frequencies differed significantly among the four populations. Phylogeographic analysis of the haplotypes suggested that black bears in this region have gone through two different colonisation histories, since the observed haplotypes belonged to two major monophyletic lineages and the lineages were distributed with an apparent border. The spatial genetic structure revealed by using mtDNA was different from that observed using microsatellite DNA markers, probably due to female philopatry and male-biased dispersal. Since nuclear genetic diversity will be lost in the three western populations because of the small population size and genetic isolation, their habitats need to be preserved, and these four populations should be linked to each other by corridors to promote gene flow from the easternmost population with higher nuclear genetic diversity.     

Using a reference population yardstick to calibrate and compare genetic diversity reported in different studies: an example from the brown bear

In species with large geographic ranges, genetic diversity of different populations may be well studied, but differences in loci and sample sizes can make the results of different studies difficult to compare. Yet, such comparisons are important for assessing the status of populations of conservation concern. We propose a simple approach of using a single well-studied reference population as a ‘yardstick'' to calibrate results of different studies to the same scale, enabling comparisons. We use a well-studied large carnivore, the brown bear (Ursus arctos), as a case study to demonstrate the approach. As a reference population, we genotyped 513 brown bears from Slovenia using 20 polymorphic microsatellite loci. We used this data set to calibrate and compare heterozygosity and allelic richness for 30 brown bear populations from 10 different studies across the global distribution of the species. The simplicity of the reference population approach makes it useful for other species, enabling comparisons of genetic diversity estimates between previously incompatible studies and improving our understanding of how genetic diversity is distributed throughout a species range.