Open Access

Barking up the wrong tree: Modern northern European dogs fail to explain their origin

  • Helena Malmström1Email author,
  • Carles Vilà1,
  • M Thomas P Gilbert3,
  • Jan Storå4,
  • Eske Willerslev3,
  • Gunilla Holmlund2 and
  • Anders Götherström1
BMC Evolutionary Biology20088:71

DOI: 10.1186/1471-2148-8-71

Received: 05 July 2007

Accepted: 28 February 2008

Published: 28 February 2008

Abstract

Background

Geographic distribution of the genetic diversity in domestic animals, particularly mitochondrial DNA, has often been used to infer centers of domestication. The underlying presumption is that phylogeographic patterns among domesticates were established during, or shortly after the domestication. Human activities are assumed not to have altered the haplogroup frequencies to any great extent. We studied this hypothesis by analyzing 24 mtDNA sequences in ancient Scandinavian dogs. Breeds originating in northern Europe are characterized by having a high frequency of mtDNA sequences belonging to a haplogroup rare in other populations (HgD). This has been suggested to indicate a possible origin of the haplogroup (perhaps even a separate domestication) in central or northern Europe.

Results

The sequences observed in the ancient samples do not include the haplogroup indicative for northern European breeds (HgD). Instead, several of them correspond to haplogroups that are uncommon in the region today and that are supposed to have Asian origin.

Conclusion

We find no evidence for local domestication. We conclude that interpretation of the processes responsible for current domestic haplogroup frequencies should be carried out with caution if based only on contemporary data. They do not only tell their own story, but also that of humans.

Background

Mitochondrial phylogeography is a useful tool for the study of wild populations [1]. But applying phylogeography to domestic species is more complicated. For example, the arrival of dogs into the New World was not a simple expansion of Asian dog populations, but a consequence of human migrations [2]. Humans have a history of population movements and social changes over the last ten thousand years [36], and this, together with the human-mediated breeding history has affected the genetic composition of livestock animals [7]. Therefore, human influence cannot be neglected when studying patterns of genetic diversity in domesticates.

A classic example of the use of mitochondrial DNA (mtDNA) diversity to infer the history of domestication refers to dogs (Canis familiaris). Four to six mitochondrial haplogroups (Hg) have been described in genetic studies of modern dogs [8], indicating recurrent domestication or backcrosses between domestic dogs and wild wolves (Canis lupus). Three of the major Hgs are distributed throughout the world, whereas one (D) is restricted to Europe, especially in breeds originating in Scandinavia [811]. Similar patterns of fragmented genetic diversity have been used to argue for local domestication in other species [12]. Such scenario could apply to dogs as they appear as early as 9000 years ago in Scandinavia [13, 14], and as dogs and wolf remains have been found on the same sites [15, 16].

We have genetically typed mtDNA from ancient Scandinavian dog remains to investigate whether the high frequency of a rare dog Hg in northern Europe is evidence of a local domestication event. We predict that an early domestication event in Scandinavia would be visible in a notable frequency of Hg D in ancient Scandinavian dog remains.

Results and Discussion

Reproducible sequence data for the full 219 bp sequence was recovered from 18 samples. The remaining six specimens only yielded reproducible sequence data for a shorter 107 bp product (DQ860843-3 DQ860864, AY673648–AY673672, Table 1). Phylogenetic reconstruction including the 18 prehistoric dogs that yielded complete sequences and 543 modern dog sequences [11] revealed that all prehistoric dogs grouped within modern Hgs (Figure 1). Seven of the 18 sequences belonged to Hg A, the group that encompasses more than 70% of all tested modern dogs around the world [11], while the remaining eleven belonged to the now rare Hg C. Previously published data on 273 Swedish dogs indicated frequencies of the four major Hgs: A 69.2%, B 18.3%, C 7.0%, and D 5.5% [9]. When only native Scandinavian breeds were considered (n = 54), the haplotype frequency for Hg D increased to 33.3% [9]. A comparison of the 219 bp sequences obtained from the ancient dog specimens (n = 18) and the contemporary native breeds indicates that the frequency of Hg D is significantly reduced in the ancient data set (0% in the ancient dataset, 33% in the modern dataset, p < 0.001). The remaining six ancient Scandinavian dogs that yielded only 107 bp of reproducible sequence data, also indicated an absence of Hg D (Table 2). In total, eleven sequences had substitution patterns indicative for Hg A, while the remaining 13 had a substitution indicative for Hg C (Table 1).
Table 1

Description of material and Hg belonging.

Sample

Locality

Element

Age

Haplogroup

21,2

Korsnäs#

Bone

Neolithic

C

31

Korsnäs#

Bone

Neolithic

A*

41,2

Korsnäs#

Bone

Neolithic

A

121,3

Bergsgraven

Bone

Neolithic

C

131,2

Ajvide

Bone

Neolithic

C

151,2

Ajvide#

Bone

Neolithic

C

161,2

Ajvide#

Bone

Neolithic

A

171,2

Ajvide

Bone

Neolithic

C

181,2

Ajvide

Bone

Neolithic

C

191,2

Ajvide#

Bone

Neolithic

C

211,2

Ajvide

Teeth

Neolithic

C

221,2

Ajvide#

Teeth

Neolithic

C

231,2

Ajvide

Teeth

Neolithic

C

251,3

Ajvide

Tooth

Neolithic

C

261,3

Skara A

Bone

Medieval

A*

271,2

Skara B

Bone

Medieval

A

281

Stockholm

Bone

Medieval

C*

302,3

Visby

Teeth

Neolithic

A

323

Eketorp

Bone

Medieval

C*

333

Eketorp

Bone

Medieval

A*

343

Eketorp

Bone

Medieval

A*

12433

Sunnerby

Tooth

Medieval

A

2293

Sunnerby

Bone

Medieval

A

2313

Sunnerby

Bone

Medieval

A

Scandinavian dog remains yielding DNA. Sample name, locality, element (tooth or bone), age, and haplogroup are specified. Those independently replicated are marked with #, those only yielding the 107 bp fragment are marked with * (otherwise, the total length of the studied sequence is 219 bp). Sequences are from Malmström et al 20051 (107 bp), Malmström et al 20072 (112 bp), and from this study3 (either or both fragments).

Table 2

Description of genetic data.

 

15595

15611

15625

15626

15638

15639

15650

15652

15665

Hg

Ref

C

T

T

A

T

T

T

G

T

-

2

·

·

·

·

G

·

C

·

·

C

3

·

·

·

·

A

·

·

·

C

A

4

·

·

·

G

A

·

·

·

C

A

12

·

·

·

·

G

·

C

·

·

C

13

·

C

·

·

G

·

C

·

·

C

15

·

C

·

·

G

·

C

·

·

C

16

·

·

·

·

A

·

·

·

C

A

17

·

C

·

·

G

·

C

·

·

C

18

·

C

·

·

G

·

C

·

·

C

19

·

C

·

·

G

·

C

·

·

C

21

·

C

·

·

G

·

C

·

·

C

22

·

C

·

·

G

·

C

·

·

C

23

·

C

·

·

G

·

C

·

·

C

25

·

C

·

·

G

·

C

·

·

C

26

·

·

·

G

A

·

·

·

·

A

27

·

·

·

G

A

·

·

·

·

A

28

·

C

·

·

G

·

C

·

·

C

30

·

·

·

·

A

·

·

·

C

A

32

·

C

·

·

G

·

C

·

·

C

33

·

·

·

G

A

·

·

·

·

A

34

·

·

·

·

A

·

·

·

·

A

1243

·

·

·

·

·

·

·

·

·

A

229

·

·

·

·

·

·

·

·

·

A

231

·

·

·

·

·

·

·

·

·

A

Haplogroup-specific substitutions in a 107 bp fragment of the D-loop (15561–15668) and their occurrence in 24 ancient Scandinavian dogs. The reference sequence is from Kim et al. 1998.

https://static-content.springer.com/image/art%3A10.1186%2F1471-2148-8-71/MediaObjects/12862_2007_Article_640_Fig1_HTML.jpg
Figure 1

Reduced median network based on 561 sequences (216 bp in length) including 18 ancient sequences (full sequences, not only those in Table 2) produced in this study and a coyote ( Canis latrans ) outgroup. After removing tri-status characters, 36 haplotypes remain; those in red contain at least one ancient Scandinavian individual. A NJ tree based on Kimura-2 and bootstrapped with 1000 pseudoreplicates yielded a similar topology where Hg A was not supported, Hg B received 81% support, Hg C 84%, and Hg D 65%. Four medieval and three Neolithic dogs are distributed among the Hg A sequences, and 11 Neolithic sequences are distributed within the C clade.

When the complete ancient dataset was compared to the contemporary dataset, it proved equally improbable (p < 0.001) that Hg D was present in the oldest specimens in a frequency similar to that observed in modern dogs. Using a binomial distribution, we estimate that the frequency of Hg D must have been lower than 0.118 in the original population for us not to have noticed it in 24 samples (p < 0.05), and lower than 0. 182 during the Neolithic for us not to have noticed it in 15 samples (p < 0.05). Lastly, the χ2 test between the Hg frequencies of the complete ancient dataset with that of the complete Scandinavian (not only considering native breeds) demonstrated a significant deviation between the modern and ancient frequencies (p < 0.0388, χ2 = 4.27, df = 1).

Our results indicate that Hg frequencies have been altered in Scandinavian dogs since their first arrival. Interestingly, while Hg C is overrepresented in our ancient material, there is a complete lack of the Scandinavian group D in our ancient dataset. Hg D is the one that could support a Scandinavian origin whereas Hg C is suggested to be of Asian origin [11]. Thus, we find no obvious evidence for prehistoric canid domestication in Scandinavia. An external origin of Scandinavian dogs is supported by morphologic data, as even the oldest remains of dogs in Scandinavia were of smaller size than those of prehistoric and extant wolves [15, 17]. While canid domestication may have occurred in other parts of Europe [18], Scandinavian dogs were likely imported and had experienced a long period of morphological change under human control before they reached the Scandinavian peninsula.

Conclusion

In a wider context, our data calls for caution when using modern sequences of domestic animals to interpret the history of domestication. There is a profound difference on selective forces on wild and domestic animals. Domesticates are subject to constant manipulation and follow the same historic vicissitudes as the associated human population. For example, genetic diversity may change rapidly due to intense selective breeding, backcrossing with the wild ancestor [1922], migration of human groups with genetically different domesticates [2], and decimation of local populations by invaders [23]. Hg frequency could change over a short period of few generations. While diversity and branching pattern of mitochondrial DNA may be used to study the number and possible antiquity of domestication events, phylogeographical patterns in domestic species may indeed not be suitable to base discussion on their geographical origin. Close collaboration between zooarchaeologists and molecular biologists, involving the study of ancient specimens, may be the only way to verify the continuity through time of the species that have set the basis of our culture.

Methods

We analysed sequences from 24 archaeological domestic dogs. The skeletal material represents eight Scandinavian sites: four Neolithic (≈5300–4500 BP) and four Medieval (≈1000–500 BP) (Table 1, Figure 2). DNA was extracted, PCR-amplified and sequenced according to previously published protocols [24, 25] (Table 1). A subset of the sample sequences (n = 7) was replicated in an independent laboratory (Center for Ancient Genetics, University of Copenhagen) confirming the original results in all cases. The new data was combined with previously published data [24, 25] for analysis (see Table 1). Most of the Neolithic samples, those from Ajvide, Visby, and Korsnäs (Figure 2), were recovered as isolated findings in cultural layers. One sample was from a complete skeleton found in a burial in Central Sweden (Bergsgraven in Linköping). The Medieval samples originated from isolated skeletal remains recovered in urban (Skara and Stockholm) or rural contexts (Eketorp, Sunnerby).
https://static-content.springer.com/image/art%3A10.1186%2F1471-2148-8-71/MediaObjects/12862_2007_Article_640_Fig2_HTML.jpg
Figure 2

Samples were collected from eight different archaeological localities in southern Scandinavia, four Neolithic (empty circles) and four medieval (black dots).

We used polymorphisms found within a 219 bp fragment of the dog D-loop region and comparison with previously published sequences to determine Hgs. To portray the relationship between sequences we constructed a reduced median network, including all 219 bp, using Network 4.1 [26]; threshold was set to 1. The default parameters for reduction were used, where parallel mutations are assumed to be more frequent between existing sequences, than between an existing sequence and a median vector. As the reduction was intense, a NJ tree was constructed to support the network using the Kimura 2-parameter model of sequence evolution. Support for the nodes was assessed with 1000 bootstrap pseudoreplicates.

The Hg frequencies within the ancient datasets were compared to those in modern Scandinavian dogs [9] using χ2 tests. To ensure a minimum sample of 5 for all cells in χ2 tests, we pooled Hgs A and D, and C and B (following the best natural grouping according to the network, Figure 1). Finally, we used the binomial distribution to estimate the probability of observing as many sequences from Hg D as found in the ancient samples assuming that the relative frequency was the same as in modern native Scandinavian dogs.

Declarations

Acknowledgements

We thank L. Drenzel at the Museum of National Antiquities, M. Douglas at Östergötland County Museum, M. Vretemark at Västergötlands Museum and A. Nyqvist Thorsson at Väner Museum for supplying bone samples, L. Tomasson for constructing Figure 2 and J. Leonard for valuable comments. The Swedish Research Council provided financial support for the study.

Authors’ Affiliations

(1)
Department of Evolutionary Biology, Uppsala University
(2)
Department of Forensic Genetics and Forensic Toxicology, National Board of Forensic Medicine
(3)
Centre for Ancient Genetics, Biological Institute, University of Copenhagen
(4)
Osteoarchaeological Research Laboratory, Stockholm University

References

  1. Hewitt GM: Genetic consequences of climatic oscillations in the Quaternary. Philos Trans R Soc Lond B Biol Sci. 2004, 359 (1442): 183-95; discussion 195. 10.1098/rstb.2003.1388.PubMed CentralView ArticlePubMedGoogle Scholar
  2. Leonard JA, Wayne RK, Wheeler J, Valadez R, Guillen S, Vila C: Ancient DNA evidence for Old World origin of New World dogs. Science. 2002, 298 (5598): 1613-1616. 10.1126/science.1076980.View ArticlePubMedGoogle Scholar
  3. Karlsson AO, Wallerstrom T, Gotherstrom A, Holmlund G: Y-chromosome diversity in Sweden - a long-time perspective. Eur J Hum Genet. 2006, 14 (8): 963-970. 10.1038/sj.ejhg.5201651.View ArticlePubMedGoogle Scholar
  4. Pakendorf B, Stoneking M: Mitochondrial DNA and human evolution. Annu Rev Genomics Hum Genet. 2005, 6: 165-183. 10.1146/annurev.genom.6.080604.162249.View ArticlePubMedGoogle Scholar
  5. Semino O, Passarino G, Oefner PJ, Lin AA, Arbuzova S, Beckman LE, De Benedictis G, Francalacci P, Kouvatsi A, Limborska S, Marcikiae M, Mika A, Mika B, Primorac D, Santachiara-Benerecetti AS, Cavalli-Sforza LL, Underhill PA: The genetic legacy of Paleolithic Homo sapiens sapiens in extant Europeans: a Y chromosome perspective. Science. 2000, 290 (5494): 1155-1159. 10.1126/science.290.5494.1155.View ArticlePubMedGoogle Scholar
  6. Thomas MG, Parfitt T, Weiss DA, Skorecki K, Wilson JF, le Roux M, Bradman N, Goldstein DB: Y chromosomes traveling south: the cohen modal haplotype and the origins of the Lemba--the "Black Jews of Southern Africa". Am J Hum Genet. 2000, 66 (2): 674-686. 10.1086/302749.PubMed CentralView ArticlePubMedGoogle Scholar
  7. Notter DR: The importance of genetic diversity in livestock populations of the future. J Anim Sci. 1999, 77 (1): 61-69.PubMedGoogle Scholar
  8. Vila C, Savolainen P, Maldonado JE, Amorim IR, Rice JE, Honeycutt RL, Crandall KA, Lundeberg J, Wayne RK: Multiple and ancient origins of the domestic dog. Science. 1997, 276 (5319): 1687-1689. 10.1126/science.276.5319.1687.View ArticlePubMedGoogle Scholar
  9. Angleby H, Savolainen P: Forensic informativity of domestic dog mtDNA control region sequences. Forensic Sci Int. 2005, 154 (2-3): 99-110. 10.1016/j.forsciint.2004.09.132.View ArticlePubMedGoogle Scholar
  10. Pires AE, Ouragh L, Kalboussi M, Matos J, Petrucci-Fonseca F, Bruford MW: Mitochondrial DNA sequence variation in Portuguese native dog breeds: diversity and phylogenetic affinities. J Hered. 2006, 97 (4): 318-330. 10.1093/jhered/esl006.View ArticlePubMedGoogle Scholar
  11. Savolainen P, Zhang YP, Luo J, Lundeberg J, Leitner T: Genetic evidence for an East Asian origin of domestic dogs. Science. 2002, 298 (5598): 1610-1613. 10.1126/science.1073906.View ArticlePubMedGoogle Scholar
  12. Jansen T, Forster P, Levine MA, Oelke H, Hurles M, Renfrew C, Weber J, Olek K: Mitochondrial DNA and the origins of the domestic horse. Proc Natl Acad Sci U S A. 2002, 99 (16): 10905-10910. 10.1073/pnas.152330099.PubMed CentralView ArticlePubMedGoogle Scholar
  13. Arnesson Westerdahl A: Mesolitiskt hundfynd vid Hornborgasjön. Värdefulla fauna- och domestikationshistoriska fynd. Flora och Fauna. 1983, 78: 211-216.Google Scholar
  14. Magnell O: Tracking Wild Boar and Hunters. Osteology of Wild Boar in Mesolithic South Scandinavia. Acrhaeology. 2005, Lund , LundGoogle Scholar
  15. Benecke N: Studies on Early Dog Remains from Northern Europe. Journal of Archaeological Science. 1987, 14 (31-49):
  16. Pluskowski A: Where are the Wolves? Investigating the Scarcity of Eurasian Grey Wolf (Canis lupus lupus) Remains in Medieval Archaeological Contexts and its Implications. International Journal of Osteoarchaeology. 2006, 16: 279-295. 10.1002/oa.824.View ArticleGoogle Scholar
  17. Jonsson L: The Vertebrate Faunal Remains from the Late Atlantic Settlement Skateholm in Scania, South Sweden. The Skateholm Project I, Man and Environment. Edited by: Larsson L. 1988, Stockholm , Almqvist & Wiksell International., 56-89.Google Scholar
  18. Verginelli F, Capelli C, Coia V, Musiani M, Falchetti M, Ottini L, Palmirotta R, Tagliacozzo A, De Grossi Mazzorin I, Mariani-Costantini R: Mitochondrial DNA from prehistoric canids highlights relationships between dogs and South-East European wolves. Mol Biol Evol. 2005, 22 (12): 2541-2551. 10.1093/molbev/msi248.View ArticlePubMedGoogle Scholar
  19. Gotherstrom A, Anderung C, Hellborg L, Elburg R, Smith C, Bradley DG, Ellegren H: Cattle domestication in the Near East was followed by hybridization with aurochs bulls in Europe. Proc Biol Sci. 2005, 272 (1579): 2345-2350. 10.1098/rspb.2005.3243.PubMed CentralView ArticlePubMedGoogle Scholar
  20. Larson G, Dobney K, Albarella U, Fang M, Matisoo-Smith E, Robins J, Lowden S, Finlayson H, Brand T, Willerslev E, Rowley-Conwy P, Andersson L, Cooper A: Worldwide phylogeography of wild boar reveals multiple centers of pig domestication. Science. 2005, 307 (5715): 1618-1621. 10.1126/science.1106927.View ArticlePubMedGoogle Scholar
  21. Svensson EM, Anderung C, Baubliene J, Persson P, Malmstrom H, Smith C, Vretemark M, Daugnora L, Gotherstrom A: Tracing genetic change over time using nuclear SNPs in ancient and modern cattle. Anim Genet. 2007,Google Scholar
  22. Vila C, Seddon J, Ellegren H: Genes of domestic mammals augmented by backcrossing with wild ancestors. Trends Genet. 2005, 21 (4): 214-218. 10.1016/j.tig.2005.02.004.View ArticlePubMedGoogle Scholar
  23. Valadez R, Mestre G: Historia del Xoloitzcuintel en México. 1999, México City , Universidad nacional Autónoma de MéxicoGoogle Scholar
  24. Malmstrom H, Stora J, Dalen L, Holmlund G, Gotherstrom A: Extensive human DNA contamination in extracts from ancient dog bones and teeth. Mol Biol Evol. 2005, 22 (10): 2040-2047. 10.1093/molbev/msi195.View ArticlePubMedGoogle Scholar
  25. Malmstrom H, Svensson EM, Gilbert MT, Willerslev E, Gotherstrom A, Holmlund G: More on contamination: the use of asymmetric molecular behavior to identify authentic ancient human DNA. Mol Biol Evol. 2007, 24 (4): 998-1004. 10.1093/molbev/msm015.View ArticlePubMedGoogle Scholar
  26. Bandelt HJ, Macaulay V, Richards M: Median networks: speedy construction and greedy reduction, one simulation, and two case studies from human mtDNA. Mol Phylogenet Evol. 2000, 16 (1): 8-28. 10.1006/mpev.2000.0792.View ArticlePubMedGoogle Scholar

Copyright

© Malmström et al; licensee BioMed Central Ltd. 2008

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Advertisement