American wolves), PCDs, modern Arctic dogs,
and the CTVT founder dog. Our analyses in-
dicated that, unlike Arctic dogs, PCDs share a
number of derived alleles with coyotes and
North American wolves, indicative of admixture
(figs. S16 and S17). The CTVT founder dog also
showed some weak evidence of coyote ancestry
but did not appear to possess admixture with
North American wolves (figs. S16 and S17). Be-
cause coyotes are restricted to North America,
this finding suggests that CTVT may have orig-
inated there. As we did not ascertain the degree
of coyote ancestry in ancient PCD-related dogs in
Northern Siberia (such as the Zhokhov Island
dogs) (Fig. 1), however, this analysis does not
establish the location in which CTVT originated.
Furthermore, studies that used somatic muta-
tions to reconstruct the phylogeography of the
CTVT clone indicated deep divergence in Asia
and rece nt introduction to the Americas (21).
Altogether, these results suggest a scenario in
which CTVT originated in Asia from a dog that
was closely related to PCDs, although we cannot
exclude the possibility that the clone arose in
America and then dispersed early into Asia be-
fore being reintro duced to America.
ThelegacyofPCDsinmodernAmericancanid
populations is uncertain. It has been suggested
that some North American wolves obtained a
mutation leading to black coat color possib ly via
admixture with early American dogs (23). This
allele was not present, however, in either of the
two higher-coverage ancient PCDs in this study
(3)orinCTVT(20). Additional ancient genomes
are necessary to determine if this allele was
present in the PCD population.
In addition, previous studies have argued that
some modern American dog populations possess
a genetic signature fromindigenousAmerican
dogs (8, 9, 24). To test this hypothesis, we
analyzed nuclear data obtained from more than
5000 modern dogs (including American village
dogs) genotyped on a 180,000 SNP array (9).
We found 7 to 20% PCD ancestry in modern
American Arctic dogs (Alaskan huskies, Alaskan
malamutes, and Greenland dogs) by using f4
ratios (tables S10 and S11) (3). This result, how-
ever, may reflect ancient population substructure
in Arctic dogs rather than genuine admixture (3).
Our f4 ratio analysis did not detect a significant
admixture signal from PCDs in any modern
American dogs of European ancestry (table S10).
Our admixture analysis detected varying de-
grees (0 to 33%) of PCD/Arctic ancestry in three
individual Carolina dogs (fig. S20). This analysis,
however, could not distinguish between PCD and
Arctic ancestry, and we cannot rule out that this
signal was a result of admixture from modern
Arctic dogs and not from PCDs (3). The majority
of modern American dog populations, including
138 village dogs from South America and mul-
tiple “native” breeds (e.g., hairless dogs and
Catahoulas), possess no detectable traces of PCD
ancestry (Fig. 2A, fig. S20, and table S10), though
this analysis may suffer from ascertainment bias.
To further assess the contribution of PCDs
to modern American dog populations, we also
analyzed 590 additional modern dog mitoge-
nomes, including those from 169 village and
breed dogs that were sampled in North and
South America (21). We identified two modern
American dogs (a chihuahua and a mixed-breed
dog from Nicaragua) that carried PCD mito-
chondrial haplotypes (fig. S5), consistent with
a limited degree of PCD ancestry (<2%) in
modern American dogs. We also identified
three East Asian dogs that carried a PCD hap-
lotype, possibly as a result of ancient population
substructure or recent dog dispersal (fig. S5) (3).
Although greater degrees of PCD ancestry may
remain in American dogs that have not yet been
sampled, our results suggest that European dogs
almost completely replaced native American dog
lineages. This near disappea rance of PCDs likely
resulted from the arrival of Europeans, which led
to shifts in cultural preferences and the per-
secution of indigenous dogs (25). Introduced
European dogs may also have brought infectious
diseases to which PCDs were susceptible.
The first appearance of dogs in the North Amer-
ican archaeological record occurred ~4500 years
after the earliest evidence of human activity on
the continent (4, 11). In addition, our molecular
clock analysis indicates that the PCD lineage ap-
peared ~6500 years after North American human
lineages (Fig. 1B) (10). These discrepancies suggest
that dogs may not have arrived into the Americas
alongside
the very first human migration but
were instead potentially part of a later arrival (12)
before the flooding of the Bering land bridge
~11,000 years ago (11). This timing is compati-
ble with both the archaeological record and our
PCD divergence time estimate and suggests a
scenario in which dogs were brought to the
Americas several thousand years after the first
people arrived there.
This initial dog population entered North
America and then dispersed throughout the
Americas, where it remained isolated for at least
9000 years. Within the past 1000 years, however ,
at leas t three independent reintroductions of
dogs have occurred. The first may have consisted
of Arctic dogs that arrived with the Thule culture
~1000 years ago (6). Then, beginning in the 15th
century, Europeans brought a second wave of
dogs that appear to have almost completely re-
placed native dogs. Lastly, Siberian huskies were
introduced to the American Arctic during the
Alaskan gold rush (25). As a result of these more
recent introductions, the modern American dog
population is largely derived from Eurasian breeds,
and the closest known extant vestige of the first
American dogs now exists as a worldwide trans-
missible cancer.
REFERENCES AND NOTES
1. G. Larson et al., Proc. Natl. Acad. Sci. U.S.A. 109, 8878–8883
(2012).
2. A. Perri et al., bioRxiv 343574 [Preprint]. 11 June 2018.
3. See supplementary materials.
4. T. Goebel, M. R. Waters, D. H. O’Rourke, Science 319,
1497–1502 (2008).
5. J. A. Leonard et al., Science 298, 1613–1616 (2002).
6. S. K. Brown, C. M. Darwent, E. J. Wictum, B. N. Sacks, Heredity
115, 488–495 (2015).
7. K. E. Witt et al., J. Hum. Evol. 79, 105–118 (2015).
8. B. van Asch et al., Proc. R. Soc. London Ser. B 280, 20131142
(2013).
9. L. M. Shannon et al., Proc. Natl. Acad. Sci. U.S.A. 112,
13639–13644 (2015).
10. M. Raghavan et al., Science 349, aab3884 (2015).
11. M. Jakobsson et al., Clim. Past 13, 991–1005 (2017).
12. J. V. Moreno-Mayar et al., Nature 553, 203–207 (2018).
13. Z. Fan et al., Genome Res. 26, 163–173 (2016).
14. G.-D. Wang et al., Cell Res. 26,21–33 (2016).
15. A. H. Freedman et al.,
PLOS Genet. 10,
e1004016 (2014).
16. L. A. F. Frantz et al., Science 352, 1228–1231 (2016).
17. B. M. vonHoldt et al., Nature 464, 898–902 (2010).
18. J. K. Pickrell, J. K. Pritchard, PLOS Genet. 8, e1002967
(2012).
19. M. Raghavan et al., Science 345, 1255832 (2014).
20. E. P. Murchison et al., Science 343, 437–440 (2014).
21. A. Strakova et al., eLife 5, e14552 (2016).
22. C. A. Rebbeck, A. M. Leroi, A. Burt, Science 331, 303 (2011).
23. T. M. Anderson et al., Science 323, 1339–1343 (2009).
24. H. G. Parker et al., Cell Rep. 19, 697–708 (2017).
25. M. Derr, A Dog’s History of America: How Our Best Friend
Explored, Conquered, and Settled a Continent (Farrar, Straus
and Giroux, 2005).
26. M. Ní Leathlobhair et al., The evolutionary history of dogs in the
Americas, Dryad (2018); doi: 10.5061/dryad.s1k47j4.
AC KN OW LE D GM E NT S
We thank L. Orlando, R. K. Wayne, and D. Meltzer for their valuable
comments; B. M. Kemp, M. Masson, and J. Chupasko for support;
and J. Southon (W. M. Keck Carbon Cycle Accelerator Mass
Spectrometry Laboratory, University of California, Irvine) for the
radiocarbon date on the Port au Choix dog. We acknowledge the
University of Oxford Advanced Research Computing (ARC) facility
for providing computing time. We thank the Illinois State Museum,
the Illinois State Archaeological Survey, the Glenn A. Black
Laboratory of Archaeology at Indiana University Bloomington, the
Instituto Nacional de Antropologia e Historia, and the Ohio
Historical Society for access to material. We thank The Rooms
(Museum Division), the Board Executive, and the Government of
Newfoundland and Labrador for permission to access and sample
the Port au Choix material. We are grateful to M. Ptaszynska for
useful information and to S. Zhang for assistance with samples. We
thank the staff of the Danish National High-throughput Sequencing
Centre for assistance in data generation. Funding: L.A.F.F. was
supported by the Wellcome Trust (210119/Z/18/Z) and by Wolfson
College (University of Oxford). L.A.F.F., J.H., A.L., A.H.-B., O.L.,
K.M.D., and G.L. were supported by a European Research Council
grant (ERC-2013-StG-337574-UNDEAD) or Natural Environmental
Research Council grants (NE/K005243/1 and NE/K003259/1)
or both. M.N.L. and E.P.M. were supported by Wellcome
(102942/Z/13/A) and a Philip Leverhulme Prize awarded by the
Leverhulme Trust. A.R.P. was supported by the Max Planck
Society. E.K.I.-P. was supported by a Clarendon Fund scholarship
from the University of Oxford. M.T.P.G. was supported by a
European Research Council grant (ERC-2015-CoG-681396–
Extinction Genomics). A.M. was supported by the Muséum National
d’Histoire Naturelle. K.E.W. and R.S.M. were supported by an NSF
grant (BCS-1540336) and a Wenner-Gren grant. V.G. was
supported by a Social Sciences and Humanities Research Council
Insight grant. V.V.P., E.Y.P., and P.A.N. were supported by Russian
Science Foundation project N16-18-10265-RNF. Y.-M.K. was
supported by a Herchel Smith research fellowship. S.J.C. was
supported by Millennia Research. J.J. was supported by the Santa
Barbara Museum of Natural History. A.R.B. was supported by the
American Kennel Club and the NIH. We thank the Illinois State
Museum Society for funding. Author contributions: L.A.F.F., G.L.,
and E.P.M. conceived of the project and designed the research;
A.R.P., K.M.D., and G.L. coordinated the archaeological analyses
and sample collection efforts with input from R.S.M., C.A., A.H.-B.,
and K.E.W.; A.R.P., C.A., J.B., E.G., A.J.H., M.-H.S.S., S.J.C., M.E.,
Y.N.C., V.G., J.J., A.K.K., P.A.N., C.P.L., A.M., T.M., K.N.M., M.O.,
E.Y.P., P.S., V.V.I., C.W., and V.V.P. provided and/or collected
samples; K.E.W., A.L., J.H., O.L., S.B., A.D., E.A.D., J.E., J.-M.R.,
and M.-H.S.S. conducted the ancient laboratory work with input
from R.S.M., G.L., L.A.F.F., E.W., I.B., and M.T.P.G.; M.M., E.P.M.,
and A.S. provided and/or collected CTVT samples; M.N.L. and
Y.-M.K. conducted the CTVT analyses with input from E.P.M., K.G.,
and L.A.F.F.; M.N.L., L.A.F.F., and E.K.I.-P. conducted the analyses
of ancient data with input from S.G., A.K., A.R.B., and E.P.M.;
and L.A.F.F., G.L., E.P.M., M.N.L., and A.R.P. wrote the paper with
input from all other authors. Competing interests: A.D., J.E.,
and J.-M.R. are employees of Arbor Biosciences, which provided
target enrichment kits used in this study. J.-M.R. is also a founder
of Arbor Biosciences. A.R.B. is the founder and chief strategy
Ní Leathlobhair et al., Science 361,81–85 (2018) 6 July 2018 4of5
RESEARCH | REPORT
on December 3, 2018 http://science.sciencemag.org/Downloaded from