Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2016 Feb;26(2):163-73.
doi: 10.1101/gr.197517.115. Epub 2015 Dec 17.

Worldwide patterns of genomic variation and admixture in gray wolves

Affiliations

Worldwide patterns of genomic variation and admixture in gray wolves

Zhenxin Fan et al. Genome Res. 2016 Feb.

Abstract

The gray wolf (Canis lupus) is a widely distributed top predator and ancestor of the domestic dog. To address questions about wolf relationships to each other and dogs, we assembled and analyzed a data set of 34 canine genomes. The divergence between New and Old World wolves is the earliest branching event and is followed by the divergence of Old World wolves and dogs, confirming that the dog was domesticated in the Old World. However, no single wolf population is more closely related to dogs, supporting the hypothesis that dogs were derived from an extinct wolf population. All extant wolves have a surprisingly recent common ancestry and experienced a dramatic population decline beginning at least ∼30 thousand years ago (kya). We suggest this crisis was related to the colonization of Eurasia by modern human hunter-gatherers, who competed with wolves for limited prey but also domesticated them, leading to a compensatory population expansion of dogs. We found extensive admixture between dogs and wolves, with up to 25% of Eurasian wolf genomes showing signs of dog ancestry. Dogs have influenced the recent history of wolves through admixture and vice versa, potentially enhancing adaptation. Simple scenarios of dog domestication are confounded by admixture, and studies that do not take admixture into account with specific demographic models are problematic.

PubMed Disclaimer

Figures

Figure 1.
Figure 1.
Sample distribution. Solid circles are samples sequenced in this study. Open circles indicate sequences from Zhang et al. (2014). Triangles and boxes indicate sequences from Wang et al. (2013) and Freedman et al. (2014), respectively. Species memberships are indicated by color: gray wolf (red), domestic dog (blue), coyote (green), and golden jackal (yellow). The reference dog genome is from a boxer.
Figure 2.
Figure 2.
Total length of runs of homozygosity (ROHs) and heterozygosity. The black line is the total length of ROHs (Mb) in each genome, and the blue and red bars are the genome-wide heterozygosity with and without ROHs, respectively.
Figure 3.
Figure 3.
The maximum likelihood tree of 30 sequences. Numbers represent node support inferred from 100 bootstrap repetitions. The reference genome boxer was not included. The Israeli golden jackal is the outgroup.
Figure 4.
Figure 4.
Principal component analyses. (A) PC1 and PC2 of dogs and 20 wolves; (B) PC1 and PC2 of dogs and 18 wolves, excluding the Tibetan wolf 1 and Qinghai wolf 1; (C) PC3 and PC4 of dogs and 20 wolves; (D) PC3 and PC4 of dogs and 18 wolves, excluding the Tibetan wolf 1 and Qinghai wolf 1. (□) Highland Asian wolves; (▵) lowland Asian wolves; (○) Middle Eastern wolves; (■) European wolves; (▲) dogs; (●) North American wolves.
Figure 5.
Figure 5.
Demographic history inferred using PSMC. Following Freedman et al. (2014) and Zhang et al. (2014), we used a generation time = 3 and a mutation rate = 1.0 × 10−8 per generation. The Tibetan wolf 1 and Inner Mongolia wolf 4 are shown in all the plots for comparison purposes. (A) All the Asian wolves; (B) all the European wolves, Middle Eastern wolves, and Indian wolf; (C) dogs; (D) Mexican wolf and Yellowstone wolves.
Figure 6.
Figure 6.
Demographic model inferred using G-PhoCS. Estimates of divergence times and effective population sizes (Ne) inferred by applying a Bayesian demography inference method (G-PhoCS) to sequence data from 13,647 putative neutral loci in a subset of 22 canid genomes (because of limitations in computational power). Estimates were obtained in four separate analyses (Methods; Supplemental Table 6). Ranges of Ne are shown and correspond to 95% Bayesian credible intervals. Estimates are calibrated by assuming a per-generation mutation rate of μ = 10−8. Mean estimates (vertical lines) and ranges corresponding to 95% Bayesian credible intervals are provided at select nodes. Scales are given in units of years by assuming an average generation time of 3 yr and two different mutation rates: μ = 10−8 (dark blue) and μ = 4 × 10−9 (brown). The model also considered gene flow between different population groups (see Table 1).

Similar articles

Cited by

References

    1. Aggarwal RK, Ramadevi J, Singh L. 2003. Ancient origin and evolution of the Indian wolf: evidence from mitochondrial DNA typing of wolves from Trans-Himalayan region and Pennisular India. Genome Biol 4: P6.
    1. Aldenderfer M. 2011. Peopling the Tibetan plateau: insights from archaeology. High Alt Med Biol 12: 141–147. - PubMed
    1. Anderson TM, vonHoldt BM, Candille SI, Musiani M, Greco C, Stahler DR, Smith DW, Padhukasahasram B, Randi E, Leonard JA, et al. 2009. Molecular and evolutionary history of melanism in North American gray wolves. Science 323: 1339–1343. - PMC - PubMed
    1. Boyko AR, Quignon P, Li L, Schoenebeck J, Degenhardt JD, Lohmueller KE, Zhao K, Brisbin A, Parker HG, vonHoldt BM, et al. 2010. A simple genetic architecture underlies morphological variation in dogs. PLoS Biol 8: e1000451. - PMC - PubMed
    1. Brantingham PJ, Rhode D, Madsen DB. 2010. Archaeology augments Tibet's genetic history. Science 329: 1467. - PubMed

Publication types

Substances

LinkOut - more resources