Burying Dogs in Ancient Cis-Baikal, Siberia: Temporal Trends and Relationships with Human Diet and Subsistence Practices

Abstract

The first objective of this study is to examine temporal patterns in ancient dog burials in the Lake Baikal region of Eastern Siberia. The second objective is to determine if the practice of dog burial here can be correlated with patterns in human subsistence practices, in particular a reliance on terrestrial mammals. Direct radiocarbon dating of a suite of the region’s dog remains indicates that these animals were given burial only during periods in which human burials were common. Dog burials of any kind were most common during the Early Neolithic (∼7–8000 B.P.), and rare during all other time periods. Further, only foraging groups seem to have buried canids in this region, as pastoralist habitation sites and cemeteries generally lack dog interments, with the exception of sacrificed animals. Stable carbon and nitrogen isotope data indicate that dogs were only buried where and when human diets were relatively rich in aquatic foods, which here most likely included river and lake fish and Baikal seal (Phoca sibirica). Generally, human and dog diets appear to have been similar across the study subregions, and this is important for interpreting their radiocarbon dates, and comparing them to those obtained on the region’s human remains, both of which likely carry a freshwater old carbon bias. Slight offsets were observed in the isotope values of dogs and humans in our samples, particularly where both have diets rich in aquatic fauna. This may result from dietary differences between people and their dogs, perhaps due to consuming fish of different sizes, or even different tissues from the same aquatic fauna. This paper also provides a first glimpse of the DNA of ancient canids in Northeast Asia.

About the Authors:

Robert J. Losey

* E-mail: [email protected]

Affiliation: Department of Anthropology, University of Alberta, Edmonton, Alberta, Canada

Sandra Garvie-Lok

Affiliation: Department of Anthropology, University of Alberta, Edmonton, Alberta, Canada

Jennifer A. Leonard

Affiliation: Conservation and Evolutionary Genetics Group, Estación Biológica de Doñana (Consejo Superior de Investigaciones Cientificas), Seville, Spain

M. Anne Katzenberg

Affiliation: Department of Archaeology, University of Calgary, Calgary, Alberta, Canada

Mietje Germonpré

Affiliation: Department of Palaeontology, Royal Belgian Institute of Natural Sciences, Brussels, Belgium

Tatiana Nomokonova

Affiliation: Department of Anthropology, University of British Columbia, Vancouver, BC, Canada

Mikhail V. Sablin

Affiliation: Zoological Institute, Russian Academy of Science, Saint-Petersburg, Russia

Olga I. Goriunova

Affiliation: Department of Archaeology and Ethnography, Irkutsk State University, Irkutsk, Russia

Natalia E. Berdnikova

Affiliation: Department of Archaeology and Ethnography, Irkutsk State University, Irkutsk, Russia

Nikolai A. Savel’ev

Affiliation: Department of Archaeology and Ethnography, Irkutsk State University, Irkutsk, Russia

Introduction

During the Late Pleistocene, people in Eurasia began ritually burying domesticated dogs, indicating that humans’ close personal interactions with these animals have a very deep, and likely a very complex, history [1]. Skeletal remains of ancient dogs from around the globe are increasingly being studied to better understand these animals’ evolution, life histories, and relationships with humans. Cis-Baikal, the region encompassing the southern and western shores of Eastern Siberia’s Lake Baikal and the lands surrounding the upper portions of the Angara and Lena rivers west of this lake (Figure 1), has rich collections of archaeological canid remains, including dog and wolf skeletons from graves and disarticulated remains from habitation sites [2]. This region also has one of North Asia’s most thoroughly studied Holocene archaeological records, which provides a firm culture history framework within which human-dog interaction can be explored [3].

[Figure omitted. See PDF.]

Figure 1. Map of the study area.

Figure 1a. Location of study area within the Russian Federation. Figure 1b. Site locations mentioned in the text. Numbers 1–8 are in the Angara River/South Baikal region, and numbers 9–17 are in Priol’khon’e: 1. Shamanka II, 2. Lokomotiv, 3. Kitoi, 4. Ust’-Khaita, 5. Ust’-Belaia, 6. Pad Lenkovka, 7. Pad Kalashnikova, 8. Ust’-Ida, 9. Bugul’deika II, 10. Khotoruk, 11. Ulan-Khada, 12. Uliarba II, 13. Sarminskii Mys, 14. Khuzhir-Nuge XIV, 15. Kurma XI, 16. Shamanskii Mys, 17. Todakta I.

https://doi.org/10.1371/journal.pone.0063740.g001

In a previous examination of canid remains from the region, we proposed that Cis-Baikal foragers were burying select dogs and wolves because these animals had attained near-human status, and because the culturally appropriate way to care for such individuals at death that emerged during the Middle Holocene was interment in cemeteries [2]. At the time only two of the region’s canids had been directly dated, both to the Middle Holocene, and it was unknown if dogs only were buried by the foragers living here during this period, or if they were interred throughout the Holocene, including by pastoralists, who first arrived ∼3400 B.P (all dates are in calibrated years before present). Further, no detailed studies had been made on the diets of ancient dogs in Cis-Baikal and how these might have varied temporally, geographically, and in relation to human dietary and subsistence patterns. As a result, it was unknown if the practice of burying dogs could be related to the particular ways in which dogs were utilized by people in the past. For example, dog burials potentially could be most common among foraging groups that were heavily dependent on terrestrial game like deer, which can be effectively hunted with dogs. Correspondingly, dog interment then might be uncommon among foraging groups relying mostly on fishing, which likely would not have directly involved the use of dogs. Alternatively, dog burials might be widespread among pastoralists who became emotionally tied to certain dogs through their use in herding, livestock protection, and hunting. Long-term practices of dog burial in relation to such factors rarely have been directly examined with archaeological data.

This paper provides the first in-depth analysis of the archaeological context, age, and diets of an array of ancient Cis-Baikal dogs in order to asses how the practice of dog burial relates to the chronology of human mortuary practices, dietary patterns, and subsistence strategies over the course of the last ∼12,000 years. To do this we provide new radiocarbon dates, stable isotope values, and ancient mitochondrial DNA (mtDNA) data for the canid specimens, and integrate previously published data on other regional canids. The contextual information and new technical data point to a series of previously unrecognized geographical and temporal patterns in dog burials in this region, but also a first glimpse at the DNA of early dogs in Northeast Asia. To best understand these datasets, the general culture history of the region is first discussed, including brief notes on the previous identification of dogs from each period, and trends in human dietary and subsistence practices are outlined.

Background

The Early Holocene in Cis-Baikal, referred to as the Mesolithic, is characterized by the use of microblade technologies and the absence of pottery and cemeteries (Table 1). One wolf burial and one unburied dog have been reported from this period [4]. Forager cemeteries first appear here ∼8000 B.P. at the start of the Early Neolithic period, which also is said to be marked by the use of pottery and ground stone technologies (but not agriculture; [3], [5]). Several dogs have been identified from the Early Neolithic, including burials in cemeteries [2]. Around 7000 to 6800 B.P. the Early Neolithic mortuary traditions cease–nearly no humans are interred here for ∼1000 years. At 6000 to 5800 B.P., with the advent of the Late Neolithic period, human burials again become common, but are of different mortuary traditions. Dog remains have not been reported for the Late Neolithic.

[Figure omitted. See PDF.]

Table 1. Simplified Holocene culture history model for Cis-Baikal.

https://doi.org/10.1371/journal.pone.0063740.t001

The time between the Early and Late Neolithic mortuary traditions has been referred to as the ‘hiatus’ or Middle Neolithic, and genetic research has demonstrated that the latter human populations were genetically discontinuous with the Early Neolithic mortuary populations [6], [7]. No dog remains have been identified from the hiatus period. A third major period of forager cemetery use spans from 5200/5000 to 4000 B.P., or the Early Bronze age, and these groups may be culturally and genetically derived from local Late Neolithic populations. Dog remains previously assigned to this period include burials and isolated elements [2]. After ∼3400 B.P., pastoralists primarily herding sheep, goats, cattle, and horses arrive and establish a suite of new burial traditions [8], [9]. The chronology of Cis-Baikal’s Late Holocene culture history is subdivided here into the Late Bronze Age (∼3400 to 2250 B.P), Early Iron Age (∼2250 to 1350 B.P.), Late Iron Age (1350 to 850 B.P.), and Early Mongolian (850 to 550 B.P.) periods. No dog remains have been described from any of these periods.

Human dietary behaviors during the Holocene in Cis-Baikal have been assessed using stable carbon and nitrogen isotope analyses of human skeletal remains, and through the study of faunal remains from habitation sites. For the Mesolithic, fauna at habitation sites suggest a focus on hunting large ungulates, including deer [10]. Faunal remains and stable carbon and nitrogen isotope values for the region’s Middle Holocene human foragers indicate that these groups relied primarily on terrestrial game such as red deer (Cervus elaphus) and roe deer (Capreolus pygargus) and increased but variable amounts of the region’s freshwater fauna, including riverine and lake fish and Lake Baikal seal (Phoca sibirica) [11]–[16]. Overall, protein diets at this time were least based on aquatic foods among foragers living on the Upper Lena River, and more aquatically focused among those living on the shores of Lake Baikal and on the Angara River. No stable isotope data is available for Cis-Baikal’s Late Holocene human remains, but zooarchaeological data indicates subsistence economies based on herding, hunting (including for deer), sealing, and some fishing, at least along the lake shore [9].

Importantly, the Holocene culture history chronology for Cis-Baikal mostly has been built through radiocarbon dating human skeletal remains, and it is clear that its precision needs to be reevaluated because of the recent documentation of an old carbon offset in Lake Baikal’s aquatic fauna [17]. Humans and dogs that consumed aquatic foods from the lake and its outlet, the Angara River, also likely carry this offset to some extent, and potentially have produced radiocarbon age assessments centuries older than their true ages. While this project was not designed to provide a model for correcting radiocarbon dates on Cis-Baikal dog or human remains, understanding dogs’ dependence on aquatic foods relative to local humans should help to assess the relative extent of bias in the dates made on dog remains, allowing us to more firmly assign them to culture history period.

The canid specimens analyzed here come from two subregions of Cis-Baikal, the Angara River/South Baikal region, and Priol’khon’e (Figure 1). The third Cis-Baikal subregion where extensive archaeological research has occurred, the Upper Lena, has yet to produce clear evidence of dog remains from any period; only a single unidentified and undated canid has been reported [18], but was unavailable to us for study. The Angara River/South Baikal region and Priol’khon’e have quite different ecologies, and these differences in ecology likely shaped local subsistence patterns and stable isotope ecologies. The Angara River is Lake Baikal’s only outlet and is ecologically influenced by the lake, especially in terms of the fishery in its upper sections. Black grayling (Thymallus arcticus baicalensis) and lenok (Brachymystax lenok), for example, move from the lake into the river as far as Irkutsk for spawning [19]. As the Angara progresses further downstream to the northwest, its fish species take on a more typical Siberia boreal river composition. The surrounding landscape is characterized by rolling hills, a mix of forest types, and in some areas, patches of forest-steppe. The adjacent South Baikal region includes the lake’s south shoreline and the Irkut River valley, which has tributaries that drain several high-altitude mountain ranges. It generally is somewhat wetter than the Angara region, and some areas are densely forested. The second study subregion, Priol’khon’e, which is the area on the west coast of Lake Baikal around the Little Sea, is markedly more arid than first subregion, with steppe and forest-steppe vegetation being dominant. The southern end of the Little Sea, the body of water between Ol’khon Island and the west coast, is one of the lake’s larger stretches of relatively shallow (<5 m) water, and is a rich littoral fishery [20]. Hills, mountains, and cliffs characterize the coastline on Ol’khon’s east coast and that immediately north and south of the Little Sea. These shoreline regions are adjacent to some of the lake’s deepest sections and water temperatures typically are cold, even in summer. Here the density of fishes in the upper portion of the water column generally is low; Baikal seals congregate in these areas of the lake, particularly when it is ice-covered.

Materials and Methods

Archaeological and Reference Material

The canid remains analyzed in this study derive from excavations carried out as early as the 1950s. For some individuals, only small element fragments or crania now survive, limiting the extent of analyses possible. Because all specimens analyzed in this study were excavated during other projects, no permits were needed to conduct the analyses reported here. Contextual information for the canid remains was derived from original site reports, publications, and photographs, and is summarized in Tables 2 and 3.

[Figure omitted. See PDF.]

Table 2. Context information and radiocarbon dating data for Angara/South Baikal canid remains.

https://doi.org/10.1371/journal.pone.0063740.t002

[Figure omitted. See PDF.]

Table 3. Context information and radiocarbon dating data for Priol’khon’e canid remains.

https://doi.org/10.1371/journal.pone.0063740.t003

All newly reported radiocarbon dates are AMS dates obtained through the Oxford Radiocarbon Accelerator Unit [21], [22]. Radiocarbon results are summarized in Tables 2 and 3, and all were calibrated using Oxcal (v4.1) and the INTCAL09 dataset [23]. Site locations are shown on Figure 1.

Where possible, the prehistoric canids from Cis-Baikal were compared with reference sets of recent and prehistoric dogs and wolves. Standard measurements on the Cis-Baikal canid crania, mandibles, dentition, and long bones are presented in Table 4.

[Figure omitted. See PDF.]

Table 4. Metric data for Cis-Baikal canids. VDD measurements follow [73], MY measurements are from [74]. Long bone measurements are greatest length, following [73].

https://doi.org/10.1371/journal.pone.0063740.t004

The first comparative reference group we utilize consists of European Paleolithic dogs [24], [25] (Tables 5, 6, 7). The second reference group is made up of Chukotka dogs from the 19th and 20th century, hereafter termed Siberian Huskies. The third reference group consists of recent North American Arctic dogs, all from the 20th century. Reference groups of wolves measured by us include Pleistocene (Weichselian) wolves from Belgium, the Czech Republic, the Ukraine and Russia [24], and recent Siberian and Tibetan wolves. From the literature we added Weichselian wolves from southern France [26], recent North American wolves [27], and European Mesolithic dogs from England [28], Germany [29], Portugal [30], and the Baltic region [31]. Estimates of the canid’s shoulder heights, based on their long bone lengths (following [32]), are given in Table 7, also for comparative purposes.

[Figure omitted. See PDF.]

Table 5. Comparison of cranial lengths of Cis-Baikal canids with modern and ancient canids.

https://doi.org/10.1371/journal.pone.0063740.t005

[Figure omitted. See PDF.]

Table 6. Comparison of teeth size of the Ust’-Khaita canid with select canid populations. All measurements in mm.

https://doi.org/10.1371/journal.pone.0063740.t006

[Figure omitted. See PDF.]

Table 7. Comparison of Cis-Baikal canid long bone lengths with those of modern and ancient wolves, and large northern dogs. Shoulder heights are estimates based on regression equations on long bone lengths [32].

https://doi.org/10.1371/journal.pone.0063740.t007

Genetic Analyses

Bone preparation, DNA isolation and polymerase chain reaction (PCR) set-up were all performed in a dedicated, spatially isolated ancient DNA laboratory using all standard ancient DNA precautions. Negative controls were included with each extraction and PCR. DNA extraction was either with phenol-chloroform as described in [2] or with the silica column based DNeasy kit (Qiagen). Fragments of the mitochondrial DNA were amplified in 25 µl reactions including 3 µl of extract. Reactions contained 1× buffer, 2.5 mM MgCl2, 0.2 mM dNTPs, 1 µM of each primer and 2.5 U AmpliTaq Gold (Applied Biosystems). Primers which target the ca. 425 base pairs at the 5′ end of the mitochondrial control region in three overlapping fragments were used from [33] (size uncertainty due to presence of indels). The PCR conditions were an initial 5 min denaturation step at 95°C followed by 20 cycles of 95°C for 30 s, touchdown 0.5°C/step from 60 to 50 for 1 min, and extension at 72°C for 1 min followed by 40 cycles of 95°C for 30s, 48°C for 1 min and 72°C for 1 min with a final elongation step of 7 min at 72°C, following [2]. Amplification was attempted from each extract at least twice for each fragment. All PCRs included negative PCR controls and the extraction negative controls. All reactions were checked on agarose gels stained with CyberSafe (Invitrogen), and all bands were directly (Sanger) sequenced in both directions with the same primers as were used in the PCR.

A BLAST search (http://www.ncbi.nlm.nih.gov/genbank) was made with each clean sequence to determine taxa of origin, to ensure it was canid and not a common reagent contaminant such as human or cow [34]. Sequences of canid origin were checked against replicates of the same fragment and overlapping fragments, and concatenated in Geneious [35]. Degredation via deamination of the DNA in ancient material can yield apparent mutations. These should occur randomly, so they can be identified and resolved through replication. In the two cases where a mismatch was identified between replicates, the base pair was resolved through generation of sequences from at least one additional independent amplification. Sequences were aligned with a representative subset of previously published dog and wolf sequences [2], [33], [36], [37], [38] in MUSCLE [39] and a phylogenetic analysis was performed using a maximum likelihood algorithm and an HKY85 model of sequence evolution (as in [2]), four rate categories, and an estimated gamma parameter and transition/transversion ratio in the program PHYML [40]. Coyote (Canis latrans) haplotypes were used to root the phylogeny [37]. The minimum spanning network of the dog clade I haplotypes, as identified in the phylogeny, was constructed using TCS v. 1.21 [41].

Stable Isotope Analyses

The stable carbon (δ13C) and nitrogen (δ15N) isotope ecology of Cis-Baikal was previously established through the analyses of 251 modern and archaeological faunal specimens [11]–[15]. These analyses were conducted at the University of Calgary, using the methods outlined in [14]. Briefly, samples were cleaned ultrasonically, demineralized in 1% HCl, and soaked in 0.125 NaOH before being freeze-dried, ground and analyzed. The new canid stable isotope values presented below were obtained through the Oxford Radiocarbon Accelerator Unit using the methods outlined in [21]. These were the same samples utilized for radiocarbon dating, and the C/N ratios reported in Tables 2 and 3 were used to help evaluate sample preservation for both the radiocarbon dates and the stable isotope results. For analyses, these samples were ground, demineralized in 0.5 M HCl and treated with 0.1 M NaOH before being gelatinized and ultrafiltered to collect the >30-kD gelatin fraction, which was freeze-dried and analyzed. To test for inter-comparability of human and canid stable isotope results, three canid samples were analyzed by both laboratories, with the average difference in δ13C being 0.4‰, and δ15N 0.13‰.

Results

Osteometric Identification

Seven complete adult canid skulls are available from Holocene sites in Cis-Baikal. The Shamanka II dog skull and the Lokomotiv wolf skull were studied previously [2]. The five remaining complete skulls, three from Ust’-Belaia and two from Pad’ Kalashnikova, have lengths ranging from 165.8–213 mm (Tables 4, 5). They are all shorter than the skulls of the recent and Pleistocene wolves in our reference groups. The Cis-Baikal skulls fit in the size range of the recent Siberian husky skulls, except for the Ust’-Belaia 4 skull, which is shorter than these modern specimens (Table 5). Except for the Ust’-Belaia 1 skull, all are shorter than those of the recent North American Arctic dogs, and all five are shorter than those of the European Palaeolithic dogs. The upper carnassials of the five Cis-Baikal skulls (Table 4) are comparable in length to the carnassials of the recent dogs in our data set and are clearly smaller than the carnassials of the fossil and recent wolves (Table 6). In other words, all this size data points to the specimens being dogs, which is unsurprising, given that all date to Middle Holocene, long after the initial domestication of dogs had began.

While the cranium of the Ust’-Khaita canid was mostly complete, it was a juvenile at the time of death, likely 5–8 months of age, making it difficult to compare the shape and size of its cranium to adult specimens of known taxonomic status. All its permanent teeth are fully erupted, and we thus focus on the sizes of its teeth in our attempts to identify this specimen (Table 6). The length of the P3 falls inside the ranges of the P3 lengths of Palaeolithic dogs and recent North American Arctic dogs. This premolar is slightly longer than the P3 from the Cis-Baikal Middle Holocene dogs and European prehistoric dogs, but it is clearly smaller than those of the Pleistocene wolves. The length of the upper carnassial of the Ust’-Khaita canid falls in the range of the tooth lengths of the recent North American Arctic dogs, outside the range of the other dogs, but is smaller than the upper carnassials of all wolves, fossil and recent, in our data set. The upper first molar length falls inside the range of the lengths of this tooth in the European Palaeolithic dogs, the Cis-Baikal Middle Holocene dogs and the recent North American Arctic dogs; it is somewhat longer than this molar of recent Siberian Huskies. It falls just inside the range of recent Siberian wolves, but is smaller than this tooth in the Pleistocene wolves. Finally, the length of the lower carnassial of the Ust’-Khaita canid is similar to that of Palaeolithic dogs and European prehistoric dogs, but is longer than that of the other prehistoric and recent dogs. It falls at the lower limit of the range of the Tibetan wolves but is clearly smaller than this tooth of Pleistocene wolves and recent Siberian wolves. Its post-canine teeth also exhibit crowding, a trait seen elsewhere in early dogs that is relatively rare in wolves, but this could be due to its cranium and mandible not being fully developed [42], [43]. Perhaps also relevant to its identification is the fact that its cranium was punctured near the suture between the right parietal and temporal bones (Figure 2); this partially healed prior to death, perhaps indicating it had received care from humans. Also, the blade of the right ilium was fractured and it too had partially healed. It is possible both injuries occurred in the same event. Given this life history and size data, we only identify this canid as a probable dog. Note that the original analysis of this specimen also identified it as a dog [4]. Clearly, however, additional analyses will be needed in the future to more specifically identify this specimen.

[Figure omitted. See PDF.]

Figure 2. The cranium of the Ust’-Khaita canid with the location of the partially healed puncture to the cranium indicated.

https://doi.org/10.1371/journal.pone.0063740.g002

Four fragmented canid crania are available from Cis-Baikal, including specimens from Ust’-Belaia, Uliarba II, Ulan-Khada and Todakta I. However, the lengths of their carnassials are clearly shorter than those of all the wolves from our reference groups, and are similar to the lengths of the dog carnassials (Tables 4, 6). Furthermore, the size of the long bones of several of these canids fall in the range of North American Arctic dogs and are clearly smaller than the same elements in wolves (Table 7). Based on the size of their carnassials and long bones, these canid remains also are identified as dogs.

The Bugul’deika canid is represented by only a single fragmented mandible, but its carnassial is shorter than that of all of the other Cis-Baikal canids, except that of the Late Iron Age Todakta I dog, and is clearly far shorter than any wolf carnassials in our dataset (Tables 4, 6). We thus identify this specimen as a dog. Small postcranial fragments are all that remain from the Shamanskii Mys canid skeletons, but the original site reports identify all of these specimens as dogs that were similar in size and shape to modern Siberian Huskies [44]; identification methods are not provided. For this paper we also consider these specimens as dogs. Finally, the Khotoruk canid was represented only by fragments of postcranial elements and was not previously analyzed; it could not be more specifically identified through osteometric analyses.

Genetic Analyses

Of the 14 canid remains analyzed, complete sequences were obtained from five individuals: two from Ust’-Belaia (2010–021 and 2010-019), and one each from Ulan-Khada (2010-001), Pad’ Kalashnikova (2010-023) and Khotoruk (1997.282), and those sequences have been submitted to GenBank (KC776175-KC776179). The first four of these came from remains that were osteometrically identified as dogs and were either the same as a dog in Genbank, or one base pair different. The animals from Ulan-Khada and Pad’ Kalashnikova had the same haplotype, which is widespread and common in living dogs, and has been documented in America (GenBank accession HQ126706; [45]) and Africa (GQ375187; [46]). The two dogs from Ust’-Belaia had sequences that differed from this common sequence by one mutation, and which were not present in the data base (Figure 3). All three of the dog haplotypes identified here belong to dog clade I, the most common and diverse clade (Figure 4). The last individual, from Khotoruk, was represented by fragments only and could not be osteometically identified. It had a sequence that matched a Russian wolf sequence (GenBank accession GQ376507), and its isotopic signature is consistent with this identification (see below), so we identify it as a wolf.

[Figure omitted. See PDF.]

Figure 3. Minimum spanning network of haplotypes in dog clade I.

Previously published sequences are in white and the ancient Baikal dogs presented here are in the shaded ovals. The American clade has been collapsed and is represented by the oval labeled “American clade”. Each link represents a single mutation, and bars across them additional mutations. The black circle represents a hypothetical haplotype. Not all alternative links are shown.

https://doi.org/10.1371/journal.pone.0063740.g003

[Figure omitted. See PDF.]

Figure 4. Maximum likelihood phylogeny of recent and ancient dogs and wolves (prefix lu), rooted with coyotes (prefix la).

The dogs form the previously identified four clades (labeled dog clades I–IV) within the diversity of wolves. Specimen 2010-19 is Ust’-Belaia dog 2 in the text, 2010-021 is Ust’-Belaia dog 3, 2010-023 Pad’ Kalashnikova dog 1, 1997.282 Khotoruk, and 2010-001 Ulan-Khada.

https://doi.org/10.1371/journal.pone.0063740.g004

Temporal Variability in Dog Burials

Several patterns emerge when examining the radiocarbon dates and contextual information for the Angara River/South Baikal canids (Tables 1, 2, 3). First, the earliest possible dog in this subregion is at Ust’-Khaita, a habitation site, and this animal was not buried upon its death. Three previous radiocarbon dates on unidentified animal bones suggested the layer containing the canid skeleton dated to the Early Holocene [4], [10]. However, a radiocarbon date (Ox23873; 10375±45) on a fragment of the canid cranium indicates it died around 12,380 to 12,135 B.P., at the end of the Pleistocene. Notably, all other faunal remains from the layer, which include specimens of Bovidae, Cervidae, and Equidae, consist of scattered and fragmented skeletal elements. The canid is the only animal in the layer represented by a nearly whole skeleton. If these canid remains truly are from a dog, it would be among the earliest yet documented in East Asia, with specimens identified as dogs and of similar age previously being reported in China and Kamchatka [47]; a far older probable dog recently has been documented in the Altai region west of Cis-Baikal [48], [49].

Second, canid burials first appear in Angara River/South Baikal subregion near the Mesolithic/Early Neolithic transition, beginning with the wolf interred in the Lokomotiv cemetery (Tables 1, 2). Dog burials from the Early Neolithic are well evidenced and include the previously described whole skeleton in the Shamanka II cemetery [2], and the two primary burials at Pad’ Kalashnikova, a site which contains both habitation site materials and human burials. Multiple artifacts were buried with the body of the pit #2 dog at Pad’ Kalashnikova, and a pebble was placed in its mouth (Figure 5). The second dog burial at this site was interred in the sitting or crouched position but was not associated with artifacts (Figure 6). Perhaps the most interesting dog burial was found in the Ust’-Belaia habitation and cemetery site. This dog (Ust’-Belaia dog 1) was interred wearing a necklace of eight red deer canine tooth pendants and also associated with its skeleton were a Bovidae scapula and horn core, two whole roe deer antlers, and other unidentified bones (Figure 7); this dog appears to have been buried during the Early to Middle Neolithic transition. Ust’-Belaia also contained one dog dating to the Early Neolithic which appears to have been discarded in a trash pit, and one dog of unknown provenience within the site also dating to this period. In addition, one dog dating to the Middle Neolithic was present in Ust’-Belaia, this one also apparently discarded in a pit and found partially disarticulated. Third and finally, our studies provide no evidence for dog burials during the Middle Neolithic, and also absent are dog remains from any contexts dating to the Late Neolithic or later periods in this subregion.

[Figure omitted. See PDF.]

Figure 5. Pad’ Kalashnikova dog from pit #2 under excavation.

Stone and bone implements are present near and under the cranium, and a round pebble is visible within the mouth.

https://doi.org/10.1371/journal.pone.0063740.g005

[Figure omitted. See PDF.]

Figure 6. Pad’ Kalasnikova dog from pit #1.

This dog was buried in a crouched or sitting position.

https://doi.org/10.1371/journal.pone.0063740.g006

[Figure omitted. See PDF.]

Figure 7. Illustration of dog burial #1 at Ust’-Belaia.

The dog was interred wearing a necklace of four red deer canine teeth pendants and placed beside its body were various faunal remains. Redrawn from [67].

https://doi.org/10.1371/journal.pone.0063740.g007

The temporal variability in dog burial practices in Priol’khon’e follows a similar pattern to that seen in the South Baikal/Angara River region, but also provides some indication of dog postmortem treatment during the latter portions of the Holocene (Tables 1, 3). First, the earliest canid burials in Priol’khon’e are found at the Shamanskii Mys habitation and cemetery site, and these date to the Early Neolithic period. These include two primary burials (only one of which is directly dated) of dogs within an elaborate grave containing the remains of a human adult male (which also is directly dated to the Early Neolithic; Figure 8), and a single primary interment of a probable dog in its own pit near another human grave. This latter dog burial also contained multiple artifacts and fragmented faunal remains.

[Figure omitted. See PDF.]

Figure 8. Illustration of grave #3 (1972) at Shamanksii Mys.

Figure shows the skeleton of one dog placed above and to the right of the human burial, and a second dog skeleton present to the left of the human body. Note that the sewn birch bark sheet separating the human remains and the dog skeletons is not visible in this illustration, but several of the birch poles that were below this sheet are depicted on the chest of the underlying human burial. Redrawn from [44].

https://doi.org/10.1371/journal.pone.0063740.g008

Second, dog remains from any context are absent from the Middle and Late Neolithic periods in Priol’khon’e, but are evidenced in later periods. A secondary dog burial was found within the upper portion of a human grave at the Uliarba II cemetery (Figure 9), and like the human remains at this site, dates to the Early Bronze Age. Scattered and unburied remains of another Early Bronze Age dog were found at the Ulan-Khada habitation site, while an isolated dog mandible dating to the Late Bronze age was present at the Bugul’deika II habitation site. The scattered post-cranial remains of the wolf found within the upper portions of an undated grave at the Khotoruk cemetery proved to date to the Early Iron Age, far post-dating the radiocarbon dated human remains at this site, all of which are from the Early Neolithic. Third and finally, the dog from the Todakta I cemetery dates to the Late Iron Age, and this primary burial was found within its own pit near an undated human grave typologically assigned to the Early Mongolian period (Figure 10). Directly under the skeleton of this juvenile dog were the cranium, mandibles, front lower limbs, and caudal vertebrae of a juvenile cow, and smaller fragments of sheep or goat bone. This group of remains may be from a sacrifice, where a young dog was killed and buried with remnants of several eaten animals, including perhaps the hide of the cow with the head and feet still attached.

[Figure omitted. See PDF.]

Figure 9. Remains of a secondary dog burial within the upper portion of a human grave at Uliarba II.

https://doi.org/10.1371/journal.pone.0063740.g009

[Figure omitted. See PDF.]

Figure 10. The Todakta I dog interment under excavation [72].

Just below the dog were the cranium, mandible and lower leg bones of a calf. These calf remains may have been left attached to the hide and used to wrap the dog, which appears to have been a sacrificed animal.

https://doi.org/10.1371/journal.pone.0063740.g010

Stable Isotope Analyses

The stable isotope ecology of Cis-Baikal is well documented but quite complex (Figure 11). In general, fish and Baikal seal here are more enriched in the heavier nitrogen isotope (δ15N mean for arch. seal = 13.79±1.19‰; modern fish = 11.08±2.02‰) than terrestrial mammals such as deer (δ15N mean for archaeological specimens = 5.14±1.05‰). Variability in the nitrogen isotope values of humans and canids thus indicates the relative contribution of aquatic foods to their diets. Fish in Lake Baikal range widely on the δ13C scale (−28.6‰ to −9.6‰), reflecting variability in algae at the base of the aquatic food web. That variability results from different sources of carbon (dissolved organic and inorganic carbon as well as atmospheric CO2) and variation in photosynthesis related to water temperature and pH [11], [50]. Fish from shallower waters tend to be enriched in the heavier isotope of carbon relative to fish from open waters. While there is some evidence for a trophic level effect in Baikal of about 1‰ for δ13C, the majority of variation in δ13C appears to relate to where the fishes were foraging within the lake, not their trophic level. Riverine fish, by contrast, show little variation in δ13C since multiple carbon sources are constantly mixed. Human and canid variability on the δ13C scale should then relate primarily to the relative reliance upon fish from these varying environments.

[Figure omitted. See PDF.]

Figure 11. Plot of δ13C and δ15N values for select Cis-Baikal mammals and fish.

Green squares represent archaeological ungulates; green hexagons are modern Eurasian ground squirrel and hare; white circles are modern Angara fish; blue triangles are modern fish from the open Baikal coast; blue squares modern fish from the Little Sea of Lake Baikal; white triangles modern Lena River fish; pink triangles archaeological Baikal seals. Modern samples are adjusted by 1.0‰ on the carbon scale to compensate for the industrial isotope effect.

https://doi.org/10.1371/journal.pone.0063740.g011

Stable isotope values also are available for 282 Middle Holocene human skeletons from South Baikal (n = 62), the Angara River (n = 119) and Priol’khon’e (n = 101; [14], [15]; Figures 12, 13, 14). Previous studies have reported stable isotope values for six Cis-Baikal canids (including two from Shamanskii Mys and one from Khotoruk, also analyzed here), all of which at that time were undated, and three (from Khotoruk, Sagan-Zaba II, and Obkhoi) never identified in the original site reports as specifically dog or wolf [14], [15]. Two canids in these earlier studies (from Sagan-Zaba II and Obkhoi) are excluded here because they remain undated and unidentified. We include the Khotoruk canid isotope values because we were able to genetically identify this specimen as a wolf (see above), and also discussed here are the isotope values from the three Shamanskii Mys canids, all previously identified as dogs in the literature [44]. All canid stable isotope values used here are listed in Tables 8 and 9.

[Figure omitted. See PDF.]

Figure 12. Plot of δ13C and δ15N values for humans and a dog buried at the Shamanka II cemetery.

https://doi.org/10.1371/journal.pone.0063740.g012

[Figure omitted. See PDF.]

Figure 13. Plot of δ13C and δ15N values for humans and canids from the Angara River area.

Linear regression for human values (y = 0.929x +28.692, R2 = 0.846; N = 119) and canid values (y = 0.399x +17.958, R2 = 0.864; N = 8) indicated.

https://doi.org/10.1371/journal.pone.0063740.g013

[Figure omitted. See PDF.]

Figure 14. Plot of δ13C and δ15N values for humans and dogs from the Priol’khon’e area.

https://doi.org/10.1371/journal.pone.0063740.g014

[Figure omitted. See PDF.]

Table 8. Stable carbon and nitrogen isotope values for Angara River/South Baikal area canids analyzed in this study.

https://doi.org/10.1371/journal.pone.0063740.t008

[Figure omitted. See PDF.]

Table 9. Stable carbon and nitrogen isotope values for Priol’kon’e canids analyzed in this study.

https://doi.org/10.1371/journal.pone.0063740.t009

Moving to results by subregion, the human stable isotope values at the Shamanka II cemetery in South Baikal, all from Early Neolithic individuals, have relatively high δ13C and δ15N values (Figure 12; δ13C mean = −16.27±0.7‰; δ15N mean = 14.55±1.0‰), indicating diets with substantial aquatic components from the lake [15]. The Early Neolithic dog buried at this site has a δ13C value (−16.1‰) very near the site mean, and a δ15N value (13.0‰) well below the site mean (Figure 12). This pattern, where dogs are similarly enriched in δ13C and lower in δ15N than local humans eating aquatic food-rich diets, generally occurs across all of our study samples.

On the Angara River, human stable isotope values are available for Early and Late Neolithic and Early Bronze Age humans, and these fall into a clear linear pattern when plotted (y = 0.929x +28.692, R2 = 0.846; n = 119; Figure 13). The Early Neolithic humans from Lokomotiv (n = 72) are more enriched in δ13C and δ15N (δ13C mean = −15.75±0.8‰; δ15N mean = 14.15±0.7‰; Figure 13) while the Late Neolithic and Early Bronze Age individuals from Ust’-Ida (n = 41) about 115 km downstream are less enriched in both isotopes (δ13C mean = −18.14±0.9‰; δ15N mean = 11.60±0.8‰). Ust’-Belaia is located between these two cemeteries on the Angara River, and the stable isotope values for four humans buried there (three Early Neolithic, one Early Bronze) mostly plot between those of the two larger cemetery samples. Values for one Early Bronze Age Lokomotiv human and one Early Neolithic individual from the Kitoi cemetery also are available. Overall, there is a clear geographical and temporal trend in human diets on the Angara during the Middle Holocene: Early Neolithic individuals’ diets had greater aquatic content than did the diets of Late Neolithic or Early Bronze Age groups, and diets generally were less aquatic as one moved downstream from Lake Baikal.

The canid stable isotope values from the Angara region include those from dogs dating to the Upper Paleolithic/Mesolithic transition, and the Early and Middle Neolithic periods, as well as the Lokomotiv wolf, which dates to the very late Mesolithic (Figure 13). The early Ust’-Khaita probable dog, the Lokomotiv wolf, and the Ust’-Belaia dog buried wearing a necklace all have relatively low isotope values, indicating relatively low aquatic content in their diets. The remaining Angara dogs, including buried and non-buried individuals, have higher δ13C and δ15N values; the higher δ15N values suggest their diets included relatively greater quantities of aquatic foods. The stable isotope values for all the canids fall in a linear pattern when plotted (y = 0.399x+17.958, R2 = 0.864; n = 8), and this pattern differs from that seen in the Angara human isotopes values. Specifically, as canid and human isotope values increase in δ13C, the canids tend to show lower δ15N when compared to humans with the same δ13C values.

For Priol’khon’e, stable isotope values are available for four humans dating to the Early Neolithic, with the remainder being Late Neolithic (n = 8) and Early Bronze Age (n = 89; [14], [15]; Figure 14). These data indicate substantial variability in reliance upon aquatic foods from the lake, with δ15N values ranging from 10.3‰ to 17.4‰. The δ13C values range from −20.1‰ to −16.9‰, with one individual with a δ13C value of 15.0‰ identified previously as an outlier [15]. While clearly less linear than the Angara River human dataset, the Priol’khon’e human δ13C and δ15N values are still positively correlated (with outlier removed, R2 = 0.542; y = 1.485*x +41.506). The dataset indicates no significant changes in human diets across the three time periods represented.

Stable isotope values are available for eight Priol’khon’e canids, including three Early Neolithic dogs, two Early Bronze Age dogs, one Late Bronze Age dog, and one probable wolf and one dog from the Iron Age (Figure 14). The Neolithic and Early Bronze Age dogs, all presumably living with foragers, have only a slightly higher average δ15N value (mean = 13.2±0.4‰) than that of the two later dating dogs (mean = 12.4±0.5‰), suggesting these later dating pastoralists’ dogs relied marginally less on aquatic foods than the older dating dogs. The range in δ13C values for all dogs (−18.6‰ to –17.0‰) overlaps well with upper portion of the range of human δ13C values; these humans appear to have had the greatest reliance on aquatic foods within the Priol’khon’e population. The dogs also tend to have lower δ15N values than humans with similar δ13C values, similar to the pattern seen among some of the Angara River dogs; these dogs appear to have had diets rich in protein from aquatic foods. The Khotoruk specimen’s isotope values mark it as a clear outlier, and indicate that it made little use of aquatic foods, consistent with its genetic identification as a wolf.

Discussion

The earliest dog potentially evidenced in Cis-Baikal’s archaeological record is represented by scattered skeletal remains at Ust’-Khaita, and dates to the Late Pleistocene, a period in which human subsistence practices appear to have been focused on large ungulates, including deer. The first canid burial, which dates to just over 8000 years ago, is the wolf interred with human remains at the Lokomotiv cemetery. Dog burials, some with human skeletons and some in their own burial pits, appear several centuries later in both study subareas (but not the Upper Lena) within the Early Neolithic period, at sites Shamanskii Mys, Shamanka II, and Pad’ Kalashnikova. This same period is characterized by the widespread appearance of human graves and cemeteries across Cis-Baikal, and other ritualized treatments of animal remains (Table 1).

Human diets in the Angara/South Baikal region during the Early Neolithic were variable, but generally included some aquatic foods, and faunal assemblages from this subregion and time period are dominated by deer and fish [10], [51]. The buried dogs found there, all dating to the Early Neolithic period, had variable diets, with some relying heavily on aquatic foods, others much less so. Human diets in Priol’khon’e during the Middle Holocene were variable, primarily in regard to individuals’ reliance on local aquatic foods [14], [15]. Forager habitation sites on the Little Sea coast of Priol’khon’e are dominated by littoral fishes, with terrestrial mammals and seals being present in trace amounts only [16], [52]. Site faunal assemblages on the open coast of Priol’khon’e from any part of the Holocene are dominated by seal remains, but domesticated ungulates are common in deposits post-dating ∼3000 B.P [9], [53]. Small quantities of deer remains are present in these sites, particularly on the open coast. All Middle Holocene dogs analyzed from this subregion, including buried and non-buried individuals, had diets with substantial aquatic content, with a mean δ13C value similar to that of the Priol’khon’e foragers with aquatic-rich diets, and a mean δ15N value slightly below this group’s average–they were eating much like the local humans who were most reliant on aquatic foods.

Dog burials appear to have been absent in the third Cis-Baikal subregion, the Upper Lena, and Middle Holocene foragers buried there had the lowest mean δ15N values across all of Cis-Baikal (mean = 11.0‰; [15]), suggesting a greater reliance on terrestrial mammals than in the other two study subregions. Further, there also is no evidence of dog burials during the Early Holocene, when subsistence adaptations seem to be focused on ungulates. The advent and continued practice of foragers burying dogs in Cis-Baikal then is clearly temporally correlated with the practice of regularly burying humans (and sometimes other animals) in this region, but not with a single dietary pattern, and certainly not with a human dietary focus on terrestrial mammals. The opposite is true–dogs were most commonly buried in areas where foragers tended to have diets rich in aquatic foods. In Priol’khon’e and South Baikal, these aquatic foods likely included littoral fishes and Baikal seals, the latter of which dogs could have been used to help procure. On the Angara, however, aquatic fauna are largely limited to riverine fishes, and dogs likely were of little use in obtaining these. Clearly, no single human subsistence practice or diet is correlated with the practice of burying dogs in Cis-Baikal. If anything, dogs were more commonly buried where diets were broader as a result of use of both terrestrial and aquatic fauna.

Offsets between the δ13C and δ15N values of dogs and humans have been noted in previous studies, with the most common pattern being one where dogs show lower δ15N values, or both lower δ13C and lower δ15N (see Table S1). In the Angara River dataset, dog δ15N values are offset from those of humans only when both species are eating diets rich in aquatic foods (Figure 13). In Priol’khon’e, all dogs had diets rich in protein from aquatic sources, and on average have lower δ15N values than humans with similar δ13C values (Figure 14). The most commonly suggested cause of lower δ15N values in dogs than in humans is that dogs were consuming human feces, which was thought to be depleted in 15N compared to the human diet (see references in Table S1). However, recent studies of δ13C and δ15N values in the feces of humans fed controlled diets found no significant difference in the stable isotope ratios of the feces and the diets [54]. This suggests that dogs eating feces is unlikely to account for the δ15N offsets seen in our data. Another proposed explanation for isotopic offsets between canids and humans has been bone collagen ingestion [55]. In mammals, the δ15N value of bone collagen is essentially the same as the muscle, while the δ13C value of bone collagen is enriched by ∼4‰ relative to the muscle [56]. This is a poor fit to our data because dogs consuming more collagen from mammal bones than humans would show higher δ13C, not lower δ15N, and the offset should also have affected dogs eating diets with little aquatic protein content; such a pattern clearly is not evident in our samples. Less is known about the isotopic values of fish bone relative to muscle. At least one study indicates that fish bone has a lower δ15N value than fish muscle while another suggests it is higher [57], [58]. An alternative explanation for the offsets in our samples may be that when fish were eaten, humans preferentially kept and consumed larger fish, which would have higher δ15N values, and gave the smaller fish to their dogs. Although humans and canids likely shared many foods, their diets would not have been identical in almost any circumstance. Humans may have retained some foods for themselves and canid diets presumably included a number of items such as mammal bones, fish bones and viscera that humans consumed less frequently, potentially leading to isotopic differences.

While some Lake Baikal fauna consumed by humans clearly carry an old carbon offset when radiocarbon dated [9], what factors contribute to this offset, its potential variability within the lake fauna, and its magnitude in fish in the lake’s tributaries remain unknown. Regardless, the degree to which humans and dogs relied on local aquatic foods likely will correlate to some degree with the magnitude of offset in their radiocarbon ages. The general similarity in Cis-Baikal dog and human diets suggests both would carry an old carbon bias when radiocarbon dated, with dietary differences between the two species potentially producing different magnitudes of age offset. The Lokomotiv wolf grave provides a useful example of maximum magnitude of offset to expect between humans and canids. Humans from Lokomotiv had diets rich aquatic-foods (δ13C mean = −15.75±0.8‰; δ15N mean = 14.15±0.7‰), while the wolf was eating very little of these foods (δ13C = −20.2‰; δ15N = 9.7‰). Correspondingly, the two radiocarbon dates on human remains buried with the wolf were ∼500 radiocarbon years older than the dates on the canid’s skeleton [2]; this magnitude of offset might characterize most of the radiocarbon dates on human remains in this cemetery. The Angara River dogs in this study, all but one of which date to the Early or Middle Neolithic (the Ust’-Khaita canid, excluded here), also have lower mean δ15N and δ13C values (δ13C mean = −16.8±1.5‰; δ15N = 11.3±0.6‰) than the Angara Early Neolithic human group (δ13C mean = −15.8±0.8‰; δ15N average = 14.2±0.7‰), but these dietary differences are of far lower magnitude than that between the Lokomotiv wolf and humans (Figure 13). This suggests that the Angara River dogs’ age offsets are of lower magnitude than that inferred for the human remains for this region. In Priol’khon’e, dog and human diets also are relatively similar. However, the dogs here appear to have eaten at slightly higher trophic levels than at least some locally-buried humans, and at lower levels than many others, suggesting there is a more complex pattern in the errors to be expected from the dates on human and canid remains in this region. Regardless, we believe the data suggest the chronology for the Cis-Baikal dogs is not biased to the extent that the overall temporal patterns in dog burial practices cannot still be inferred with reasonable confidence. Precise corrections to Cis-Baikal’s radiocarbon dates is probably only possible through future paired dating of directly associated remains of terrestrial herbivores and dog or human skeletons.

Proceeding with the assumption that the canid radiocarbon dates are providing a reasonably reliable indication of their place in the culture history model, additional temporal patterns can be examined. No dog remains examined in this study date to the Late Neolithic period, and even if all dog radiocarbon dates were decreased in age by 3–4 centuries, only a single specimen, Ust’-Belaia dog #3, would fall within this period as it is currently defined, and then only marginally so. Unadjusted, this specimen dates firmly within the Middle Neolithic period, when human burials are very rare; this specimen also is not from a formal burial. It is possible that dogs were rare during this period in general, or perhaps the pattern is a result of curation and sampling issues. Within Cis-Baikal, the Angara River area appears to have contained the greatest number of Late Neolithic sites, but most of these were excavated in the 1970s or earlier and the collections are now lost [59]. Human burials have been dated to this period, but are notably far fewer in number than those from the Early Neolithic or Early Bronze Age [3]. Perhaps the lack of Late Neolithic dogs relates to the overall dearth of Late Neolithic archaeological materials now available.

We also have no primary dog burials dating to the Early Bronze Age, although a possible secondary dog burial was present within a human grave at Uliarba II, and other dog remains in the region clearly date to this period. Further, an Early Bronze Age primary dog burial at the Pad’ Lenkovka cemetery on the Angara River previously has been reported [60] (this specimen could not be relocated for study), suggesting that at least some dogs were buried during this period. Overall, dog burials were most common during the Early Neolithic, absent in the Middle Neolithic, and rare in the Late Neolithic and Early Bronze Age. This trend occurs on the Angara, where human diets had apparently shifted to less reliance on fish by the start of the Late Neolithic, and is observed in Priol’khon’e, where there is no evidence for a change in human dietary patterns across the entire Middle Holocene.

The Late Holocene in Cis-Baikal, which is marked in part by the presence of pastoral groups and various human mortuary traditions, also lacks clearly identifiable dog burials in human cemeteries. The possible exception is the Todakta I dog, which we have interpreted as part of the sacrifice of several domesticated animals associated with human burials at this site. Pastoralist groups by definition lived in close association with many domesticated species, here most often sheep, goats, horses, and cattle. These groups regularly rode, sold, traded, sacrificed, and consumed these animals, and considered them property. Perhaps because dogs in these societies no longer had the unique position of being humans’ only cohabitant animal, and because people’s relationships to animals more broadly had changed with the emergence of pastoralism, dogs no longer were considered to have spiritual equivalency with humans and were no longer considered eligible for burial in human cemeteries.

The mtDNA recovered from five of the Cis-Baikal canids provides a first picture of the genetics of ancient dogs from East Asia, and also another wolf. The dog specimens yielding mtDNA all belonged to dog clade I, which also includes Middle and Late Holocene archaeological dogs from Europe, West Asia, Alaska, and Central and South America [33], [47], [61] and modern dogs on several continents [45], [46]. Our samples show this clade was present in East Asia since at least the Middle Holocene. Notably, most Pre-Columbian American dogs also belong to clade I [33], [62], [63]. Further, the mtDNA haplogroup identified in an Early Neolithic dog from Pad’ Kalashnikova was also found in an Early Bronze Age dog from Ulan-Khada. This may suggest at least some continuity in dog populations here during the Middle Holocene, a period in which discontinuity in human haplogroups has been documented [6], [7], although this haplotype is currently widespread and so could have been common and/or widespread in the past as well.

The wolf haplotypes identified in our samples do not appear to have contributed genetically to the local dog population, but additional DNA analyses of ancient wolves in the area are needed to further verify this. The haplotype in the Khotoruk wolf matches a haplotype found in living Russian wolves, indicating some level of population continuity for the last 1400 years or so.

Conclusion

Dogs appear to have lived among humans in Cis-Baikal throughout much of the Holocene, and perhaps since the very Late Pleistocene. Middle Holocene foraging groups here seem to have been the only societies that provided dogs formal burials within cemeteries otherwise used for human interment. In other words, these groups treated select dogs much like humans upon their deaths, interring some in their own graves with various accoutrements, and others in graves also bearing humans. This occurred where human and dog diets were relatively broad and included significant aquatic components, and not where and when they were primarily focused on ungulates. Dogs would seemingly have been useful in hunting these terrestrial mammals, but there is no evidence for them being buried in any period on the Upper Lena, nor are there any documented dog burials during the Early or Late Holocene, the former characterized by a human focus on hunting ungulates, the latter by pastoralism, ungulate hunting, and fishing. Even where and when dogs were buried by foragers, simple relationships between human dietary patterns and dog burials is not evidenced, as dogs were buried where groups were particularly reliant on riverine fish and ungulates, and also among groups who were hunting seals, fishing in Lake Baikal, and taking some deer. Human emotional attachment to dogs alone cannot account for these patterns, as both humans and dogs almost certainly would have had the capacity for such close personal relationships across all of Cis-Baikal and throughout the last 12,000 years. The evidence presented here indicated that the foragers in Cis-Baikal burying dogs in cemeteries did this only when they were regularly burying their human dead (the Middle Holocene) in such settings. Dogs were the only domesticated animals living with humans at this time, and it appears that dogs and animals such as bears [64] were considered by foraging groups here to be spiritually similar to humans, as were many animals among historic northern indigenous groups [65]. When these broader beliefs were combined with intimate personal relationships with dogs, which here involved sharing many of the same foods, and the broader practice of burying one’s group members in cemeteries, some dogs were given ‘human’ mortuary rites.

How the Cis-Baikal dogs were being used by their human counterparts remains unclear, but it is apparent from the way that at least some of these animals were treated at death that they considered more than mere bio-technologies–they were parts of people’s emotional and social lives. Dogs in such settings likely held a number of roles, even within their individual lifetimes, ranging from companionship or guardian to fellow hunter and burden carrier. Further details of these roles can only emerge when the life histories of ancient dogs are explored through detailed studies of their skeletons’ chemistry and morphology, and with the integration of this data with archaeological contextual information. Stable isotope studies of dog diets are one step in this process, and help move the study of domestication solely from a search for deep origins to an exploration of human-animal interaction.

Supporting Information

[Figure omitted. See PDF.]

Table S1.

1a. Comparison of human and dog collagen stable carbon and nitrogen isotope values from other studies. Criterion for inclusion in the table is that the data represent multiple dogs and humans clearly from the same site and time period. Table S1b. Studies excluded in table S1a but which were reviewed because they have multiple isotope values for both dogs and humans.

https://doi.org/10.1371/journal.pone.0063740.s001

(DOCX)

Acknowledgments

The authors would like to acknowledge the helpful comments on an earlier draft of this document provided by Torben Rick and Greger Larson and two anonymous reviewers. Thanks are also offered to Andzej Weber, German Medvedev, Vladimir Bazaliiskii, and Sergei Cheremisin, who provided advice, access to materials, and logistical support. Irina Uvarova of the Archive of the Institute of Archaeology, Russian Academy of Science graciously provided access to the Pad’ Kalashnikova field report and photographs. Finally, we would like to acknowledge Natal’ia Kasprishina, who illustrated the line drawings used in the paper. Logistical support was provided by Laboratorio de Ecología Molecular, Estación Biológica de Doñana, CSIC (LEM-EBD). All necessary permits were obtained for the described study, which complied with all relevant regulations.

Author Contributions

Conceived and designed the experiments: RJL SGL JAL TN. Performed the experiments: RJL JAL MAK MG. Analyzed the data: RJL SGL JAL MAK MG TN. Contributed reagents/materials/analysis tools: MVS OIG NEB NAS. Wrote the paper: RJL SGL JAL MAK MG TN.

Citation: Losey RJ, Garvie-Lok S, Leonard JA, Katzenberg MA, Germonpré M, Nomokonova T, et al. (2013) Burying Dogs in Ancient Cis-Baikal, Siberia: Temporal Trends and Relationships with Human Diet and Subsistence Practices. PLoS ONE 8(5): e63740. https://doi.org/10.1371/journal.pone.0063740

References

1. Morey D (2010) Dogs: Domestication and the Development of a Social Bond Cambridge: Cambridge University Press. p. 356.

2. Losey RJ, Bazaliiskii VI, Garvie-Lok S, Germonpre M, Leonard J, et al. (2011) Canids as persons: Early Neolithic dog and wolf burials, Cis-Baikal, Siberia. J Anthropol Archaeol 30: 174–189.

3. Weber AW, Katzenberg MA, Schurr TG (2010) Prehistoric Hunter-gatherers of the Baikal Region, Siberia. Philadelphia: University of Pennsylvania Museum of Archaeology and Anthropology. p. 319.

4. Klement’ev AM, Igumnova EA, Savel’ev NA (2005) Khishchniki (Carnivora, Mammalia) Ust’-Khaitinskogo arkheologicheskogo mestonakhozhdeniia. Istoki, Formirovanie i Razvitie Evraziiskoi Polikul’turnosti. Kul’tury i Obshchestva Severnoi Azii v Istoricheskom Proshlom i Sovremennosti. Irkutsk: ‘‘Radian’’. 26–29.

5. Weber AW, Bettinger RL (2010) Middle Holocene hunter–gatherers of Cis-Baikal, Siberia: an overview for the new century. J Anthropol Archaeol 29: 491–506.

6. Mooder KP, Schurr TG, Bamforth FJ, Bazaliiski VI, Savel’ev NA (2006) Population Affinities of Neolithic Siberians: A Snapshot From Prehistoric Lake Baikal. Am J Phys Anthropol 129: 349–361.

7. Mooder KP, Weber AW, Schurr TG, Bamforth FJ, Bazaliiski VI, et al. (2005) Matrilineal affinities and prehistoric Siberian mortuary practices: a case study from Neolithic Lake Baikal. J Archaeolog Sci 4: 619–634.

8. Kharinskii AV (2001) Predbaikal’e v kon. I tys. do n.e. – ser. II tys n.e.: Genesis Kultur i Ikh Periodizatsiia. Irkutsk: Izdatel’stvo IrGTU. p.199.

9. Nomokonova T, Losey RJ, Weber AW, Goriunova OI (2010) Late Holocene subsistence practices among Cis-Baikal pastoralists, Siberia: zooarchaeological insights from Sagan-Zaba II. Asian Perspec 49: 157–179.

10. Savel’ev NA, Teten’kin AV, Igumnova ES, Abdulov TA, Ineshin EM, et al.. (2001) Mnogosloinyi geoarkheologicheskii ob’’ekt Ust’-Khaita – predvaritel’nye dannye. In Sovremennye Problemy Evraziiskogo Paleolitovedeniia. Izdatel’stvo Instituta Arkheologii i Etnografii (SO RAN, Novosibirsk). 338–352.

11. Katzenberg MA, Weber AW (1999) Stable isotope ecology and paleodiet in the Lake Baikal region of Siberia. J Archaeolog Sci 26: 651–659.

12. Katzenberg MA, Goriunova OI, Weber A (2009) Paleodiet reconstruction of Bronze Age Siberians from the mortuary site of Khuzhir-Nuge XIV, Lake Baikal. J Archaeolog Sci 36: 663–674.

13. Katzenberg MA, Bazaliiskii VI, Goriunova OI, Savel’ev NA, Weber AW (2010) Diet reconstruction of prehistoric hunter-gatherers in the Lake Baikal region. In Weber AW, Katzenberg MA, Schurr TG, editors. Prehistoric Hunter-gatherers of the Baikal Region, Siberia. Philadelphia: University of Pennsylvania Museum of Archaeology and Anthropology. 175–191.

14. Katzenberg A, McKenzie HG, Losey RJ, Goriunova OI, Weber AW (2011) Prehistoric dietary adaptations among hunter-fisher-gatherers from the Little Sea of Lake Baikal, Siberia, Russian Federation. J Archaeolo Sci 39: 2612–2626.

15. Weber AW, White D, Bazaliiskii VI, Goriunova OI, Savel’ev NA, et al. (2011) Hunter–gatherer foraging ranges, migrations, and travel in the middle Holocene Baikal region of Siberia: Insights from carbon and nitrogen stable isotope signatures. J Anthropol Archaeol 30: 523–548.

16. Losey RJ, Nomokonova T, Goriunova OI (2008) Fishing Ancient Lake Baikal, Siberia: Inferences from the Reconstruction of Harvested Perch (Perca fluviatilis) Size. J Archaeol Sci 38: 577–590.

17. Nomokonova T, Losey RJ, Goriunova OI, Weber AW (2013) A freshwater old carbon offset in Lake Baikal, Siberia and problems with the radiocarbon dating of archaeological sediments: Evidence from the Sagan-Zaba II site. Quatern Int 290–1: 110–125.

18. Weber AW, Link DW, Katzenberg MA (2002) Hunter–gatherer culture change and continuity in the Middle Holocene Cis-Baikal, Siberia. J Anthropol Archaeol 21: 230–299.

19. Khozov MM (1950) Presnye vody Vostochnoi Sibiri (bassein Baikala, Angary, Vitima, verkhnego technical Leny I Nizhnei Tunguski). Irkutsk: Irkutskoe Oblastnoe Gosudarstvennoe Izdatel’stvo. p.368.

20. Kozhova OM, Izmest’eva LR (1998) Lake Baikal: Evolution and Diversity Leiden: Backhuys Publishers. p.418.

21. Bronk Ramsey C, Higham T, Bowles A, Hedges R (2004) Improvements to the pretreatment of bone at Oxford. Radiocarbon 46: 155–63.

22. Brock F, Higham TFG, Ditchfield P, Bronk Ramsey C (2010) Current pretreatment methods for AMS radiocarbon dating at the Oxford Radiocarbon Accelerator Unit (ORAU). Radiocarbon 52: 103–112.

23. Reimer PJ, Baillie MGL, Bard E, Bayliss A, Beck JW, et al. (2009) IntCal09 and Marine09 radiocarbon age calibration curves, 0–50,000 years cal BP. Radiocarbon 51: 1111–1150.

24. Germonpré M, Lázničková-Galetová M, Sablin M (2012) Palaeolithic dog skulls at the Gravettian Předmostí site, the Czech Republic. J Archaeol Sci 39: 184–202.

25. Germonpré M, Sablin MV, Despres V, Hofreiter M, Laznickova-Galetova M, et al. (2013) Palaeolithic dogs and the early domestication of the wolf: a reply to the comments of Crockford and Kuzmin (2012). J Archaeol Sci 40: 786–792.

26. Boudadi-Magali M (2010) Les Canis pléistocènes du sud de la France : approche biosystématique, evolutive et biochronologique. Thèse pour obtenir le grade de Docteur, Université Bordeaux 1, Spécialité: Préhistoire et Géologie du Quaternaire.

27. Meachen JA, Samuels JX (2012) Evolution in coyotes (Canis latrans) in response to the megafaunal extinctions. P Natl Acad Sci USA 109: 4191–4196.

28. Degerbøl M (1961) On a find of a Preboreal domestic dog (Canis familiaris L.) from Star Carr, Yorkshire, with remarks on other Mesolithic dogs. P Prehist Soc 27: 35–55.

29. Street M (1989) Ein frühmesolithischer Hund und Hundverbiss an Knochen vom Fundplatz Bedburg-Königshoven, Niederrhein. Archäologische Informationen 12: 203–215.

30. Detry C, Cardoso JL (2010) On some remains of dog (Canis familiaris) from the Mesolithic shell-middens of Muge, Portugal. J Archaeol Sci 37: 2762–2774.

31. Benecke N (1994) Archäozoologische Studien zur Entwicklung der Haustierhaltung. Berlin: Akademie Verlag. p.451.

32. Harcourt RA (1974) The dog in prehistoric and early historic Britain. J Archaeol Sci 1: 151–175.

33. Leonard JA, Wayne RK, Wheeler J, Valadez R, Guillen S, et al. (2002) Ancient DNA evidence for Old World origin of New World dogs. Science 298: 1613–1616.

34. Leonard JA, Shanks O, Hofreiter M, Kreuz E, Hodges L, et al. (2007) Animal DNA in PCR reagents plagues ancient DNA research. J Archaeolog Sci 34: 1361–1366.

35. Drummond AJ, Ashton B, Buxton S, Cheung M, Cooper A, et al.. (2011) Geneious v5.4.

36. Vilà C, Savolainen P, Maldonado JE, Amorim IR, Rice JE, et al. (1997) Multiple and ancient origins of the domestic dog. Science 276: 1687–1689.

37. Vilà C, Amorim IR, Leonard JA, Posada D, Castroviejo J, et al. (1999) Mitochondrial DNA phylogeography and population history of the grey wolf Canis lupus. Mol Ecol 8: 2089–2103.

38. Sharma DK, Maldonado JE, Jhala YV, Fleischer RC (2004) Ancient wolf lineages in India. Proc R Soc Lond (Suppl 3): S1–S4.

39. Edgar RC (2004) MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res 32: 1792–1797.

40. Guindon S, Gascuel O (2003) A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol 52: 696–704.

41. Clement M, Posada D, Crandall KA (2000) TCS: a computer program to estimate gene genealogies. Mol Ecol 9: 1657–1659.

42. Lawrence B (1967) Early domestic dogs. Z Säugetierkunde 32: 44–59.

43. Clutton-Brock J (1995) Origins of the dog: domestication and early history. In Serpell J, editor. Domestic Dog: its Evolution, Behavior and Interactions with People. Cambridge: Cambridge University Press. 2–19.

44. Konopatskii AK (1982) Drevnie Kul’tury Baikala. Novosibirsk: Nauka. p. 175.

45. Castroviejo-Fisher S, Skoglund P, Valadez R, Vilà C, Leonard JA (2011) Vanishing native American dog lineages. BMC Evol Biol 11: 73.

46. Boyko AR, Boyko RH, Boyko CM, Parker HG, Castelhano M, et al. (2009) Complex population structure in African village dogs and its implications for inferring dog domestication history. Proc Natl Acad Sci USA 106: 13903–13908.

47. Larson G, Karlsson EK, Perri A, Webster MT, Ho SY, et al. (2012) Rethinking dog domestication by integrating genetics, archeology, and biogeography. Proc Natl Acad Sci USA 109: 8878–8883.

48. Ovodov ND, Crockford SJ, Kuzmin YV, Higham TF, Hodgins GW, et al. (2011) A 33,000-year-old incipient dog from the Altai Mountains of Siberia: evidence of the earliest domestication disrupted by the Last Glacial Maximum. Plos One 6: e22821.

49. Druzkhova AS, Thalmann O, Trifonov VA, Leonard JA, Vorobieva NV, et al. (2013) Ancient DNA analysis affirms the canid from Altai as a primitive dog. Plos One 8: e57754.

50. Hecky RE, Hesslein RH (1995) Contributions of benthic algae to lake food webs as revealed by stable isotope analysis. J N Am Benthol Soc 14: 631–653.

51. Medvedev GI (1969) Results of the investigation of the Mesolithic in the stratified settlement of Ust’-Belaia 1957–1964. Arctic Anthropol 6: 61–73.

52. Nomokonova T, Goriunova OI, Losey RJ, Savel’ev NA (2011) Use of Ulan-Khada bay at Lake Baikal in the Holocene (based on fauna remains). Archaeolog Ethnolog Anthropolog Eurasia 39: 30–36.

53. Goriunova OI, Ovodov ND, Novikov AG (2007) Analiz faunisticheskikh materiialov s mnogosloinogo poseleniia Tyshkine III (oz. Baikal). Severnaia Evraziia v Antropogene: Chelovek, Paleotekhnologii, Geoekologiia, Etnologiia i Antropologiia, Vol.1. Irkutsk: “Ottisk”. 168–174.

54. Kuhnle GGC, Joosen AMCP, Kneale CJ, O’Connell TC (2012) Carbon and nitrogen isotopic ratios of urine and faeces as novel nutritional biomarkers of mean and fish intake. Eur J Nutr. DOI: 10.1007/s00394-012-0328-2.

55. Rick TC, Culleton BJ, Smith CB, Johnson JR, Kennet DJ (2011) Stable isotope analysis of dog, fox and human diets at a Late Holocene Chumash village (CA-SRI-2) on Santa Rosa Island, California. J Archaeolog Sci 38: 1385–1393.

56. Ambrose SH, Norr L (1993) Experimental evidence for the relationship of the carbon isotope ratios of whole diet and dietary protein to those of bone collagen and carbonate. In Lambert JB, Grupe G, editors. Prehistoric Human Bone: Archaeology at the Molecular Level. New York: Springer-Verlag. 1–3.

57. Sholto-Douglas AD, Field JG, James AG, van der Merwe NJ (1991) 13C/12C and 15N/14N isotope ratios in the Southern Benguela Ecosystem: indicators of food web relationships among different size-classes of plankton and pelagic fish; differences between fish muscle and bone collagen tissues. Mar Ecol Prog Ser 78: 23–31.

58. Richards MP, Hedges REM (1999) Stable isotope evidence for similarities in the types of marine foods used by Late Mesolithic humans at sites along the Atlantic coast of Europe. J Archaeolog Sci 26: 717–722.

59. Weber AW (1995) The neolithic and early bronze age of the lake Baikal Region, Siberia: a review of recent research. J World Prehist 9: 99–165.

60. Okladnikov AP (1955) Neolit i bronzovyi vek Prebaikal’ia (chast’ III). Materialy i issledovaniia po arkheologii SSSR, vol. 43. Moscow: Akademiia nauk SSSR. p.208.

61. Verginelli F, Capelli C, Coia V, Musiani M, Falchetti M, et al. (2005) Mitochondrial DNA from prehistoric canids highlights relationship between dogs and South-east European wolves. Mol Biol Evol 22: 2541–2551.

62. Tito RY, Belknap III SL, Sobolik KD, Ingraham RC, Cleeland LM, et al. (2011) Brief Communication: DNA from early Holocene dog. Am J Phys Anthropol 145: 653–657.

63. Brown SK, Darwent CM, Sacks BN (2012) Ancient DNA evidence for genetic continuity in arctic dogs. J Archaeolog Sci. DOI: 10.1016/j.jas.2012.09.010.

64. Losey RJ, Bazaliiskii VI, Lieverse A, Waters-Rist A, Faccia K, Weber A (2013) The bear-able likeness of being: ursine remains at the Shamanka II cemetery, Lake Baikal, Siberia. In Watts C, editor. Relational Archaeologies. London: Routledge Press.

65. Ingold T (2000) The Perception of the Environment. London: Routledge. p.465.

66. Bazaliiskiy VI, Savelyev NA (2003) The wolf of Baikal: the ‘‘Lokomotiv’’ early Neolithic cemetery in Siberia (Russia). Antiquity 77: 20–30.

67. Krizhevskaia LIA (1978) Neolit poseleniia v ust’e r. Belaia (po materialam raskopok 1957 i 1959 gg.). Drevnie Kul’tury Priangar’ia, ed Vasilevskii PS. Novosibirsk: Nauka. 69–95.

68. Berdnikova NE (2001) Geoarkheologicheskii ob’’ekt Ust’-Belaia. Kul’turnye kompleksy. Kamennyi Vek Iuzhnogo Priangar’ia, Bel’skii Geoarkheologicheskii Raion, Vol. 2. Irkutsk: IGU. 113–146.

69. Larichev VE (1958) Otchet ob Arkheologicheskikh Raskopkakh Poseleniia v Padi Kalashnikovoi (g. Svirsk) v 1958 g. Arkhiv IA RAN, Saint-Petersburg. Holding-1, Section-1, Collection 1, #1788.

70. Goriunova OI, Novikov AG, Ziablin LP, Smotrova VI (2004). Drevnie Pogrebeniia Mogil’nika Uliarba na Baikale (Neolit-Paleometall) Novosibirsk: IAiE SO RAN. p.86.

71. Nomokonova T, Goriunova OI, Losey RJ, Savel’ev NA (2011) Use of Ulan-Khada bay at Lake Baikal in the Holocene (based on fauna remains). Archaeol Ethnol Anthropol Eurasia 39: 30–36.

72. Pavlutskaia VV (1990) Baikal’skoe Priol’khon’e v Epokhu Pozdnikh Kochevnikov (Po Materialam Mogil’nikov XI–XIX vv.). Ph.D. Summary. Leningrad: LGU. p.16.

73. Von den Driesch A (1976) A Guide to the Measurement of Animal Bones from Archaeological Sites. Cambridge MA: Peabody Museum Press. p.149.

74. Morey DF (1992) Size, shape and development in the evolution of the domestic dog. J Archaeol Sci 19: 181–204.