Skip to main content

Vanishing native American dog lineages

Abstract

Background

Dogs were an important element in many native American cultures at the time Europeans arrived. Although previous ancient DNA studies revealed the existence of unique native American mitochondrial sequences, these have not been found in modern dogs, mainly purebred, studied so far.

Results

We identified many previously undescribed mitochondrial control region sequences in 400 dogs from rural and isolated areas as well as street dogs from across the Americas. However, sequences of native American origin proved to be exceedingly rare, and we estimate that the native population contributed only a minor fraction of the gene pool that constitutes the modern population.

Conclusions

The high number of previously unidentified haplotypes in our sample suggests that a lot of unsampled genetic variation exists in non-breed dogs. Our results also suggest that the arrival of European colonists to the Americas may have led to an extensive replacement of the native American dog population by the dogs of the invaders.

Background

Dogs colonized the Americas with early human groups from Asia [1] and were widespread by the time Europeans arrived late in the 15th century [2]. Most of the dogs around the world today have mitochondrial DNA (mtDNA) control region sequences that form a well-defined phylogenetic clade (Clade I) [3]. Genetic characterization of ancient American dogs (samples obtained from human settlements that had not been in contact with Europeans, hereafter referred to as native American dogs) revealed a unique set of mtDNA sequences that clustered within this clade, but have not been observed in extant dogs [e.g. [36]]. Most notably, a subclade (Ia) has so far only been identified in ancient American dogs and had a frequency of 62% in ancient Latin American animals [1]. However, most genetic studies are based on purebred dogs, and since most internationally recognized breeds today are primarily European or Asian in origin, it is possible that American dog lineages have been excluded.

To determine contemporary distribution and frequency of the native American dog lineages in the Americas (both North America and South America), we analyzed the fragment of the mtDNA control region comparable to available genetic data for ancient native American dogs in 400 village and non-breed dogs from Alaska to Patagonia. These included dogs living in small isolated settlements in Canada, Mexico, Central America, the Caribbean, the Orinoco Llanos, the Andes, the Amazon basin and Patagonia (Figure 1). This novel data set allowed us to compare the genetic composition of past and present populations, and to use statistical population genetic models to estimate the maximum possible contribution of pre-Columbian dogs to the extant population.

Figure 1
figure 1

Map illustrating geographic distribution of modern American dog samples included in the study. Numbers indicate number sample size per country.

Results and discussion

In the total sample of 400 modern American dogs (Genbank numbers in Table S1), we identified 40 unique mtDNA haplotypes, of which 23 (57.5%) had not been identified in previous studies that included samples from around the world. Haplotypes were widespread across sampled localities and we did not detect any geographic pattern. This shows that significant undescribed diversity is likely to be present in non-purebred dogs [5]. In addition, the level of mtDNA nucleotide diversity in the present Latin American population (θπ = 0.017, = 40,000) indicates an effective population size similar to that in ancient America (θπ = 0.015, = 35,000).

More than half of the haplotypes (23/40) and individuals (259/400) belonged to the common Clade I. However, none of the sequences was identical to nor clustered within the ancient native American dog Clade Ia (Figure 2). This striking difference in haplogroup frequencies between ancient and modern dogs suggests that very little of the mtDNAs present in the extant population trace back to the native American population This could be due to either genetic drift or population discontinuity between the two time points.

Figure 2
figure 2

Maximum likelihood phylogeny of Clade I haplotypes. Ancient American dog haplotypes are on red branches. Clade Ia is highlighted with a label and the dashed rectangle highlights the only extant haplotype inferred to be derived from an ancient dog sequence. Tip labels correspond to non-redundant haplotypes from Table S1 or [1, 35, 7].

We investigated the probability of this result under different demographic models and found that the observation of a private haplogroup in the ancient sample with a frequency as high as Ia allows rejection of complete continuity between ancient and present populations for an extensive range of assumptions. For the constant size population model (see Methods), we found that all models with an assumed N e > 3,000 could be rejected at the 5% significance level, and the population expansion scenario was similarly rejected for initial N e > 2,100 (P < 0.05).

However, complete continuity is not a realistic model since introduction of European dogs is known to have occurred. Thus, we also tested the possible contribution of the ancient population under a model where the current Latin American dog population is a mixture of Old- and New World populations (Figure 3). The simulations indicated that given the estimated mitochondrial , less than 10% of the ancestors of the modern American dog population are ancient American dogs (Figure 4). While based on inferences from a single genetic marker, this upper limit for the average genetic contribution of the ancient population is expected to apply also to the autosomal genomes of extant dogs. However, we caution that the average genome-wide native American ancestry in the extant population could be different if male and female dogs contributed an unequal number of offspring to the population (i.e. there was an unequal sex ratio), the mitochondrion has been subject to selection, or if other assumptions in the demographic model (see Methods) have been violated in the recent history of the population.

Figure 3
figure 3

Coalescent simulations used to investigate the maximum possible genetic contribution of the native American dog population to the living American dog population. A) From genealogies generated under three different models, we estimated the probability that at least 8 of 13 ancient lineages coalesce to the exclusion of all modern samples. B) Schematic illustration of the isolation-admixture model.

Figure 4
figure 4

Maximum genetic contribution of native American dogs under the isolation-admixture model. Lines represent the 95 percent limit under different assumptions about effective population size (varying across the x-axis) and generation time. The results assuming a two year generation time are shown with a solid curve and the results assuming a three year generation time with a dashed curve.

We also examined if ancient American haplotypes other than those in Clade Ia were present in the modern sample, and identified a single modern haplotype derived from a pre-Columbian haplotype (Figure 5). This previously undescribed haplotype was found in two dogs from the Maya villages of Pisté and Chan Kom, both in the State of Yucatán, Mexico (M64, M76; Additional file 1: Table S1) and is most closely related to the ancient Mexican haplotype D32 [1], from which it differed by a single substitution. Although these communities are in close contact with modern life, they still live largely according to indigenous traditions, including a lack of organized dog husbandry. The mtDNA of these two individuals may be inherited from the pre-Columbian population, and thus indicate that not all ancient lineages went completely extinct.

Figure 5
figure 5

Median joining network of Clade I haplotypes. Solid ellipse highlights Clade Ia and dashed ellipse highlights the only extant haplotype inferred to be derived from an ancient American dog haplotype and that ancient haplotype. Node size is proportional to haplotype frequency. Transversal lines indicate mutations and black dots hypothetical haplotypes. Ancient American haplotypes are red, contemporary American haplotypes are yellow and contemporary dogs from elsewhere in the world are blue. Photograph inset: dog from Chan Kom (M76 in Table S1), a likely descendant of pre-Columbian native American dogs.

Arrival by Columbus to the Americas in 1492 was quickly followed by the arrival of conquistadors, missionaries and colonists from Europe with their livestock, pets, commensals and pathogens, all of which had an important impact on native American populations and culture [7, 8]. Our results show that native American dog populations were also impacted. The extent of this impact is unexpected because of the large historical population size of dogs in the Americas and the existence of potential refugia (e.g. isolated human groups) where native lineages could have survived. Several factors might have contributed to this replacement, including direct persecution [7, 9], preference for the often larger newly arrived dogs, or susceptibility to introduced infectious diseases. Future studies including more ancient and modern dogs and more genetic markers, such as neutral autosomal markers, genes of known function and Y-chromosome markers, will contribute to a deeper understanding of the causes of and extent to which the native American dog population has changed since the arrival of Columbus.

Conclusions

Using molecular data and statistical modeling, we demonstrate that an important amount of mitochondrial haplotype diversity exists in undersampled non-breed dog populations, and that the breadth of the impact of post-Columbian colonization on the Americas has been underestimated. The extensive replacement of the native American dog population inferred from our data set illustrates that even cultural and biological elements that are not specific targets of invaders can be profoundly affected at a continental scale and in a short period of time.

Methods

We used published primers (ThrL and DLHc) and PCR conditions [1] to amplify and sequence in both directions a fragment of 425 bp from the mitochondrial control region (Table S1). Sequences were compared to available sequences in GenBank with megablast to check for nuclear insertions and new haplotypes. Additional published sequences from dogs [1, 3, 4, 10] and wolves [1014] were included to place our study in a historical and geographical framework (e.g. to root phylogenetic trees) as in [14].

We used the E-INS-i option in MAFFT [15] to align all sequences included in the study. We performed phylogenetic analysis of the complete dataset with and without removing redundant haplotypes using MacClade 4.06 [16]. We ran 1000 searches and 2000 bootstrap replicates in Garli0.96 [17] to search for a maximum likelihood topology under the model TrN+I+G selected by Modeltest 3.7 [18] using the Akaike information criterion. Sequences that clustered in Clade I [3], were further analyzed by constructing haplotype networks. We focused on Clade I because all but one of the ancient American dog sequences cluster within this clade. We built median-joining [19] haplotype networks in NETWORK 4.5 http://www.fluxus-engineering.com/ and a statistical parsimony network in TCS 1.21 [20] with gaps as an informative fifth state. Both methods yielded similar results and only the result of the former analysis is shown.

We estimated mitochondrial effective population size () for the ancient sequences described by reference [1] and a corresponding subset of our modern sequences where all samples collected north of Mexico were excluded in order to better match the geographical range of the ancient dogs from reference [1] that died before Columbus first arrived to the Americas. We used nucleotide diversity π as a direct estimate of the population-scaled mutation rate θ π and the expression θ π = 2N e μ where μ is the mutation rate per nucleotide and generation (note that N e will be affected by inclusion of multiple clades). We computed π using DNAsp v.5.00.07 [21], and used a conservative estimate of μ in the dog mitochondrial control region that assumed a divergence time between gray wolves and coyotes 2 million years ago (2.13 × 10-7 per generation) [4]. We note that assuming a 1 million year divergence time (e.g. [22]) would result in higher estimates of and thus decrease the probability of observing large differences between samples from different time points even further. We also note that mtDNA data is not optimal for accurately estimating in natural populations, but its use here is motivated by our interest for the possible magnitude of genetic drift acting on American dog mitochondria, particularly in the last ~1000 years (see below), and by the fact that no autosomal sequence data is currently available.

We investigated the probability of observing a private haplogroup in the ancient sample with a frequency as high as Clade Ia in a hypothesis-testing framework for three different population genetic models for the demographic history of American dogs using Serial Sim-Coal [23, 24], the only coalescent simulator currently available that allows modeling of both population structure and samples from different time points. By simulating 10,000 independent genealogical histories and investigating the resulting topologies with custom scripts, we approximated the probability that at least 8 of 13 lineages sampled 1000 years ago (corresponding to the ancient Latin American samples) [1] would be monophyletic with respect to 299 lineages sampled at the present (corresponding to modern Latin American samples included in the study; Figure 3A), given different assumptions about N e . The general approach of estimating the possible contribution of Native American dogs to the extant gene pool in a simulation-based hypothesis-testing framework was chosen because few analytical tools (i.e. mathematical models) have been developed to deal with probabilities under the structured coalescent with temporal samples. Simulations thus provide a flexible alternative.

First, we tested a demographic model with a single continuous population that either maintained a constant size during all of its history or started growing exponentially 200 generations ago to an N e of 15,000. To investigate the possible contribution of the native American dog population to the current population, we constructed an isolation-admixture model more in agreement with historical evidence (Figure 3B). Starting from a single population at the present, the lineages in the population were divided into two isolated populations 400 years ago (representing the post-Columbian colonization of the Americas by Europeans) with probability c and 1-c, respectively. At 1000 years ago, 13 new samples were taken from one of the subpopulations (New World, representing the data available on ancient Latin American dogs [1]). At 13,000 years further in the past the American population underwent a bottleneck, implying a tenfold reduction on its effective size, and was joined together again with the Old World population at 14,000 years ago, representing the isolation of American dogs from the Eurasian populations from which they originated [1]. The exact time American and Eurasian populations were isolated is the subject of much controversy, so we use this conservative number (14,000 years before present). However, we found that the timing of this event and the severity of the bottleneck had only marginal effects on the probability of a private ancient haplogroup compared to the admixture proportions. We investigated the effect of assuming a generation time of both 2 and 3 years.

References

  1. Leonard JA, Wayne RK, Wheeler J, Valadez R, Guillén S, Vilà C: Ancient DNA Evidence for Old World Origin of New World Dogs. Science. 2002, 298: 1613-10.1126/science.1076980.

    Article  CAS  PubMed  Google Scholar 

  2. Snyder LM, Leonard JA: Dog. Handbook of North American Indians. Enviroments, origins and populations. Edited by: Sturtevant WC. 2006, Washington: Smithsonian Institution, 3: 452-462.

    Google Scholar 

  3. Vilà 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: 1687-1689.

    Article  PubMed  Google Scholar 

  4. Savolainen P, Zhang J, Luo J, Lundeberg J, Leitner T: Genetic evidence for an East Asian origin of domestic dogs. Science. 2002, 298: 1610-1613. 10.1126/science.1073906.

    Article  CAS  PubMed  Google Scholar 

  5. Boyko AR, Boyko RH, Boyko CM, Parker HG, Castelano M, Corey L, Degenhardt J, Auton A, Hedimbl L, Kityo R, Ostrander EA, Shoenebeck J, Todhunter RJ, Jones P, Bustamante CD: Complex population structure in African village dogs and its implication for inferring dog domestication history. Proc Natl Acad Sci USA. 2009, 106: 13903-13908. 10.1073/pnas.0902129106.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Pang J-F, Kluetsch C, Zou X-J, Zhang A, Luo L-Y, Angleby H, Ardalan A, Ekström C, Sköllermo A, Lundeberg J, Matsumura S, Leitner T, Zhang Y-P, Savolainen P: mtDNA Data Indicate a Single Origin for Dogs South of Yangtze River, Less Than 16,300 Years Ago, from Numerous Wolves. Mol Biol Evol. 2009, 26: 2849-2864. 10.1093/molbev/msp195.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Crosby AW: The Columbian exchange: biological and cultural consequences of 1492. 1972, Westport: Greenwood press

    Google Scholar 

  8. Diamond J: Guns, germs, and steel: the fates of human societies. 1997, New York: Norton

    Google Scholar 

  9. Valadez R, Mendoza V: El perro como legado cultural. Nuevos Aportes, Revista de Arqueología Boliviana. 2005, 2: 14-32.

    Google Scholar 

  10. Muñoz-Fuentes V, Darimont CT, Paquet P, Leonard JA: The genetic legacy of extirpation and re-colonization in Vancouver Island wolves. Cons Gen. 2010, 11: 547-556. 10.1007/s10592-009-9974-1.

    Article  Google Scholar 

  11. Vilà C, Amorim IR, Leonard JA, Posada D, Castroviejo J, Petrucci-Fonseca F, Crandall KA, Ellegren H, Wayne RK: Mitochondrial DNA phylogeography and population history of the grey wolf Canis lupus. Mol Ecol. 1999, 8: 2089-2103.

    Article  PubMed  Google Scholar 

  12. Sharma DK, Maldonado JE, Jhala YV, Fleischer RC: Ancient wolf lineages in India. Proc R Soc Lond B. 2004, 271 (Suppl 3): S1-S4. 10.1098/rsbl.2003.0071.

    Article  CAS  Google Scholar 

  13. Musiani M, Leonard JA, Cluff HD, Gates CC, Mariani S, Paquet PC, Vilá C, Wayne RK: Differentiation of tundra/taiga and boreal coniferous forest wolves, coat color and association with migratory caribou. Mol Ecol. 2007, 16: 4149-10.1111/j.1365-294X.2007.03458.x.

    Article  CAS  PubMed  Google Scholar 

  14. Muñoz-Fuentes V, Darimont CT, Wayne RK, Paquet P, Leonard JA: Ecological factors drive differentiation in wolves from British Columbia. J Biogeogr. 2009, 36: 1516-1531. 10.1111/j.1365-2699.2008.02067.x.

    Article  Google Scholar 

  15. Katoh K, Toh H: Recent developments in the MAFFT multiple sequence alignment program. Brief Bioinform. 2008, 9: 286-10.1093/bib/bbn013.

    Article  CAS  PubMed  Google Scholar 

  16. Maddison DR, Maddison WP: MacClade: analysis of phylogeny and character evolution. 2000, Vers. 4.0 Sinauer: Sunderland

    Google Scholar 

  17. Zwickl DJ: Genetic algorithm approaches for the phylogenetic analysis of large biological sequence datasets under the maximum likelihood criterion. Ph.D. Thesis. 2006, The University of Texas at Austin

    Google Scholar 

  18. Posada D, Crandall KA: Modeltest: testing the model of DNA substitution. Bioinformatics. 1998, 14: 817-10.1093/bioinformatics/14.9.817.

    Article  CAS  PubMed  Google Scholar 

  19. Bandelt HJ, Forster P, Rohl A: Median-joining phylogenies for inferring intraspecific phylogenies. Mol Biol Evol. 1999, 16: 37-

    Article  CAS  PubMed  Google Scholar 

  20. Clement M, Posada D, Crandall K: TCS: a computer program to estimate gene genealogies. Mol Ecol. 2000, 9: 1657-10.1046/j.1365-294x.2000.01020.x.

    Article  CAS  PubMed  Google Scholar 

  21. Librado P, Rozas J: DNAsp v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009, 25: 1451-1452. 10.1093/bioinformatics/btp187.

    Article  CAS  PubMed  Google Scholar 

  22. Leonard JA, Wayne RK: Native Great Lakes wolves were not restored. Biol Lett. 2008, 4: 95-98. 10.1098/rsbl.2007.0354.

    Article  PubMed  Google Scholar 

  23. Excoffier L, Novembre J, Schneider S: SIMCOAL: a general coalescent program for simulation of molecular data in interconnected populations with arbitrary demography. J Hered. 2000, 91: 506-509. 10.1093/jhered/91.6.506.

    Article  CAS  PubMed  Google Scholar 

  24. Anderson CNK, Ramakrishnan U, Chan YL, Hadly EA: Serial SimCoal: A population genetic model for data from multiple populations and points in time. Bioinformatics. 2005, 21: 1733-1734. 10.1093/bioinformatics/bti154.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We thank all that have provided samples (Table S1), L. F. Groenelveld for help constructing the network and Mattias Jakobsson for comments on a previous draft of the manuscript. This research was supported by the Swedish Research Council, Carl Tryggers Foundation, CSIC, Universidad de los Andes, Fulbright/Ministry of Education (Spain) and Programa de Captación del Conocimiento para Andalucía (Spain).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Jennifer A Leonard.

Additional information

Authors' contributions

JAL, SC-F and CV conceived and designed the study. SC-F, RV, JAL and CV collected samples. PS and SC-F performed laboratory analyses. SC-F and PS performed sequence analyses and phylogenetic analyses. PS designed and performed population genetic analyses and simulations. SC-F, PS and JAL prepared the manuscript. All authors commented on the results, and read and approved the final manuscript.

Santiago Castroviejo-Fisher, Pontus Skoglund contributed equally to this work.

Electronic supplementary material

12862_2010_1684_MOESM1_ESM.PDF

Additional file 1: Table S1: Table describing all samples used, including GenBank number, collector, country of origin and clade to which the haplotypes belongs, as defined in[3]and[4]. (PDF 96 KB)

Authors’ original submitted files for images

Rights and permissions

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

Reprints and permissions

About this article

Cite this article

Castroviejo-Fisher, S., Skoglund, P., Valadez, R. et al. Vanishing native American dog lineages. BMC Evol Biol 11, 73 (2011). https://doi.org/10.1186/1471-2148-11-73

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2148-11-73

Keywords