Grey wolf genomic history reveals a dual ancestry of dogs.

The grey wolf (Canis lupus) was the first species to give rise to a domestic population, and they remained widespread throughout the last Ice Age when many other large mammal species went extinct. Little is known, however, about the history and possible extinction of past wolf populations or when and where the wolf progenitors of the present-day dog lineage (Canis familiaris) lived1-8. Here we analysed 72 ancient wolf genomes spanning the last 100,000 years from Europe, Siberia and North America. We found that wolf populations were highly connected throughout the Late Pleistocene, with levels of differentiation an order of magnitude lower than they are today. This population connectivity allowed us to detect natural selection across the time series, including rapid fixation of mutations in the gene IFT88 40,000-30,000 years ago. We show that dogs are overall more closely related to ancient wolves from eastern Eurasia than to those from western Eurasia, suggesting a domestication process in the east. However, we also found that dogs in the Near East and Africa derive up to half of their ancestry from a distinct population related to modern southwest Eurasian wolves, reflecting either an independent domestication process or admixture from local wolves. None of the analysed ancient wolf genomes is a direct match for either of these dog ancestries, meaning that the exact progenitor populations remain to be located.

Identifying unknown Indian wolves by their distinctive howls: its potential as a non-invasive survey method.

Previous studies have posited the use of acoustics-based surveys to monitor population size and estimate their density. However, decreasing the bias in population estimations, such as by using Capture-Mark-Recapture, requires the identification of individuals using supervised classification methods, especially for sparsely populated species like the wolf which may otherwise be counted repeatedly. The cryptic behaviour of Indian wolf (Canis lupus pallipes) poses serious challenges to survey efforts, and thus, there is no reliable estimate of their population despite a prominent role in the ecosystem. Like other wolves, Indian wolves produce howls that can be detected over distances of more than 6 km, making them ideal candidates for acoustic surveys. Here, we explore the use of a supervised classifier to identify unknown individuals. We trained a supervised Agglomerative Nesting hierarchical clustering (AGNES) model using 49 howls from five Indian wolves and achieved 98% individual identification accuracy. We tested our model's predictive power using 20 novel howls from a further four individuals (test dataset) and resulted in 75% accuracy in classifying howls to individuals. The model can reduce bias in population estimations using Capture-Mark-Recapture and track individual wolves non-invasively by their howls. This has potential for studies of wolves' territory use, pack composition, and reproductive behaviour. Our method can potentially be adapted for other species with individually distinctive vocalisations, representing an advanced tool for individual-level monitoring.

Analysis of the Mitochondrial Genomes of Japanese Wolf Specimens in the Siebold Collection, Leiden.

The taxonomic status of extinct Japanese or Honshu wolves (Canis lupus hodophilax) has been disputed since the name hodophilax was first proposed by Temminck in 1839 on the basis of specimens stored in Leiden, the Netherlands. Points of controversy include whether the type specimen of hodophilax (Jentink c: RMNH.MAM.39181) and the other two specimens from Leiden (Jentink a: RMNH.MAM.39182 and Jentink b: RMNH.MAM.39183) represent different varieties or subspecies of Japanese wolves or not. Two Japanese names, ookami and jamainu, used to describe wild Canis species, further complicate the issue. In this study, the taxonomic status of Japanese wolves was clarified using mitochondrial DNA of the three specimens stored at the Naturalis Biodiversity Center in Leiden, in addition to three Japanese wolf specimens stored at the Museum für Naturkunde in Berlin and five new samples from Japan. The mitochondrial genomes of the type specimen of hodophilax (Jentink c) and another sample from Leiden (Jentink b) as well as Berlin specimens were included in the cluster of Japanese wolves distinct from other grey wolves. However, the other sample from Leiden (Jentink a) was identified as a domestic dog. A mitochondrial genome analysis suggested that Japanese wolves could be categorized into two distinct clusters. Studies of nuclear genomes are needed to further clarify the taxonomic status, divergence time, and population genetic structure of Japanese wolves.

A novel canis lupus familiaris reference genome improves variant resolution for use in breed-specific GWAS.

Reference genome fidelity is critically important for genome wide association studies, yet most vary widely from the study population. A typical whole genome sequencing approach implies short-read technologies resulting in fragmented assemblies with regions of ambiguity. Further information is lost by economic necessity when genotyping populations, as lower resolution technologies such as genotyping arrays are commonly used. Here, we present a phased reference genome for Canis lupus familiaris using high molecular weight DNA-sequencing technologies. We tested wet laboratory and bioinformatic approaches to demonstrate a minimum workflow to generate the 2.4 gigabase genome for a Labrador Retriever. The de novo assembly required eight Oxford Nanopore R9.4 flowcells (∼23X depth) and running a 10X Genomics library on the equivalent of one lane of an Illumina NovaSeq S1 flowcell (∼88X depth), bringing the cost of generating a nearly complete reference genome to less than $10K (USD). Mapping of short-read data from 10 Labrador Retrievers against this reference resulted in 1% more aligned reads versus the current reference (CanFam3.1, P < 0.001), and a 15% reduction of variant calls, increasing the chance of identifying true, low-effect size variants in a genome-wide association studies. We believe that by incorporating the cost to produce a full genome assembly into any large-scale genotyping project, an investigator can improve study power, decrease costs, and optimize the overall scientific value of their study.

Dire wolves were the last of an ancient New World canid lineage.

Dire wolves are considered to be one of the most common and widespread large carnivores in Pleistocene America1, yet relatively little is known about their evolution or extinction. Here, to reconstruct the evolutionary history of dire wolves, we sequenced five genomes from sub-fossil remains dating from 13,000 to more than 50,000 years ago. Our results indicate that although they were similar morphologically to the extant grey wolf, dire wolves were a highly divergent lineage that split from living canids around 5.7 million years ago. In contrast to numerous examples of hybridization across Canidae2,3, there is no evidence for gene flow between dire wolves and either North American grey wolves or coyotes. This suggests that dire wolves evolved in isolation from the Pleistocene ancestors of these species. Our results also support an early New World origin of dire wolves, while the ancestors of grey wolves, coyotes and dholes evolved in Eurasia and colonized North America only relatively recently.

Revisiting the Woolly wolf (Canis lupus chanco) phylogeny in Himalaya: Addressing taxonomy, spatial extent and distribution of an ancient lineage in Asia.

Of the sub-species of Holarctic wolf, the Woolly wolf (Canis lupus chanco) is uniquely adapted to atmospheric hypoxia and widely distributed across the Himalaya, Qinghai Tibetan Plateau (QTP) and Mongolia. Taxonomic ambiguity still exists for this sub-species because of complex evolutionary history anduse of limited wild samples across its range in Himalaya. We document for the first time population genetic structure and taxonomic affinity of the wolves across western and eastern Himalayan regions from samples collected from the wild (n = 19) using mitochondrial control region (225bp). We found two haplotypes in our data, one widely distributed in the Himalaya that was shared with QTP and the other confined to Himachal Pradesh and Uttarakhand in the western Himalaya, India. After combining our data withpublished sequences (n = 83), we observed 15 haplotypes. Some of these were shared among different locations from India to QTP and a few were private to geographic locations. A phylogenetic tree indicated that Woolly wolves from India, Nepal, QTP and Mongolia are basal to other wolves with shallow divergence (K2P; 0.000-0.044) and high bootstrap values. Demographic analyses based on mismatch distribution and Bayesian skyline plots (BSP) suggested a stable population over a long time (~million years) with signs of recent declines. Regional dominance of private haplotypes across its distribution range may indicate allopatric divergence. This may be due to differences in habitat characteristics, availability of different wild prey species and differential deglaciation within the range of the Woolly wolf during historic time. Presence of basal and shallow divergence within-clade along with unique ecological requirements and adaptation to hypoxia, the Woolly wolf of Himalaya, QTP, and Mongolian regions may be considered as a distinct an Evolutionary Significant Unit (ESU). Identifying management units (MUs) is needed within its distribution range using harmonized multiple genetic data for effective conservation planning.

Genomic regions under selection in the feralization of the dingoes.

Dingoes are wild canids living in Australia, originating from domestic dogs. They have lived isolated from both the wild and the domestic ancestor, making them a unique model for studying feralization. Here, we sequence the genomes of 10 dingoes and 2 New Guinea Singing Dogs. Phylogenetic and demographic analyses show that dingoes originate from dogs in southern East Asia, which migrated via Island Southeast Asia to reach Australia around 8300 years ago, and subsequently diverged into a genetically distinct population. Selection analysis identifies 50 positively selected genes enriched in digestion and metabolism, indicating a diet change during feralization of dingoes. Thirteen of these genes have shifted allele frequencies compared to dogs but not compared to wolves. Functional assays show that an A-to-G mutation in ARHGEF7 decreases the endogenous expression, suggesting behavioral adaptations related to the transitions in environment. Our results indicate that the feralization of the dingo induced positive selection on genomic regions correlated to neurodevelopment, metabolism and reproduction, in adaptation to a wild environment.

Global Phylogeographic and Admixture Patterns in Grey Wolves and Genetic Legacy of An Ancient Siberian Lineage.

The evolutionary relationships between extinct and extant lineages provide important insight into species' response to environmental change. The grey wolf is among the few Holarctic large carnivores that survived the Late Pleistocene megafaunal extinctions, responding to that period's profound environmental changes with loss of distinct lineages and phylogeographic shifts, and undergoing domestication. We reconstructed global genome-wide phylogeographic patterns in modern wolves, including previously underrepresented Siberian wolves, and assessed their evolutionary relationships with a previously genotyped wolf from Taimyr, Siberia, dated at 35 Kya. The inferred phylogeographic structure was affected by admixture with dogs, coyotes and golden jackals, stressing the importance of accounting for this process in phylogeographic studies. The Taimyr lineage was distinct from modern Siberian wolves and constituted a sister lineage of modern Eurasian wolves and domestic dogs, with an ambiguous position relative to North American wolves. We detected gene flow from the Taimyr lineage to Arctic dog breeds, but population clustering methods indicated closer similarity of the Taimyr wolf to modern wolves than dogs, implying complex post-divergence relationships among these lineages. Our study shows that introgression from ecologically diverse con-specific and con-generic populations was common in wolves' evolutionary history, and could have facilitated their adaptation to environmental change.

Identifying suitable habitat and corridors for Indian Grey Wolf (Canis lupus pallipes) in Chotta Nagpur Plateau and Lower Gangetic Planes: A species with differential management needs.

Different Biogeographic provinces and environmental factors are known to influence the dispersibility of long-ranging carnivores over the landscape. However, lack of empirical data on long-ranging carnivores may lead to erroneous decisions in formulating management plans. The Indian Grey wolf (Canis lupus pallipes) is known to be distributed in the vast areas of the Indian subcontinent. However, the actual population estimates are available only for Gujarat, Karnataka, Rajasthan and Bihar. Whereas, its distribution, population and habitat ecology is poorly known from the eastern region. Hence, this article aimed to evaluate the habitat suitability along with landscape connectivity for the species over the two major biogeographic provinces of India, i.e., Lower Gangetic Plains (7b) and Chhota Nagpur Plateau (6b). The present model with significantly higher Area under the curve (AUC) value of 0.981, indicates its accuracy in predicting the suitable habitats and identifying biological corridors by using environmental, topological and anthropogenic variables. Precipitation of the driest quarter and the precipitation of seasonality were the two best performing variables in our model, capable of explaining about 26% and 22.4% variation in the data respectively. Out of the total area i.e. 4,16,665 Km2, about 18,237 Km2 (4.37%) was found to be highly suitable area and about 3,16,803 Km2 (76.03%) areas as least suitable. The corridor analysis indicated that the habitat connectivity was highest in the border line area of the two biotic provinces located in the south-eastern zone via districts of Purba Singhbhum and Paschim Singhbhum of Jharkhand state and Bankura and West Midnapore districts of West Bengal state. Among the Protected Areas (PAs), natural corridors exist connecting the Simlipal National Park (NP)-Satkosia Wildlife Santuray (WLS), Dalma ranges of Chotta Nagpur plateau along with Badrama WLS, Khulasuni WLS and Debrigarh WLS. Differential management through landscape level planning may be helpful in securing the future of the species in the landscape.

Natural re-colonization and admixture of wolves (Canis lupus) in the US Pacific Northwest: challenges for the protection and management of rare and endangered taxa.

Admixture resulting from natural dispersal processes can potentially generate novel phenotypic variation that may facilitate persistence in changing environments or result in the loss of population-specific adaptations. Yet, under the US Endangered Species Act, policy is limited for management of individuals whose ancestry includes a protected taxon; therefore, they are generally not protected under the Act. This issue is exemplified by the recently re-established grey wolves of the Pacific Northwest states of Washington and Oregon, USA. This population was likely founded by two phenotypically and genetically distinct wolf ecotypes: Northern Rocky Mountain (NRM) forest and coastal rainforest. The latter is considered potentially threatened in southeast Alaska and thus the source of migrants may affect plans for their protection. To assess the genetic source of the re-established population, we sequenced a ~ 300 bp portion of the mitochondrial control region and ~ 5 Mbp of the nuclear genome. Genetic analysis revealed that the Washington wolves share ancestry with both wolf ecotypes, whereas the Oregon population shares ancestry with NRM forest wolves only. Using ecological niche modelling, we found that the Pacific Northwest states contain environments suitable for each ecotype, with wolf packs established in both environmental types. Continued migration from coastal rainforest and NRM forest source populations may increase the genetic diversity of the Pacific Northwest population. However, this admixed population challenges traditional management regimes given that admixture occurs between an adaptively distinct ecotype and a more abundant reintroduced interior form. Our results emphasize the need for a more precise US policy to address the general problem of admixture in the management of endangered species, subspecies, and distinct population segments.

The evolutionary history of dogs in the Americas.

Dogs were present in the Americas before the arrival of European colonists, but the origin and fate of these precontact dogs are largely unknown. We sequenced 71 mitochondrial and 7 nuclear genomes from ancient North American and Siberian dogs from time frames spanning ~9000 years. Our analysis indicates that American dogs were not derived from North American wolves. Instead, American dogs form a monophyletic lineage that likely originated in Siberia and dispersed into the Americas alongside people. After the arrival of Europeans, native American dogs almost completely disappeared, leaving a minimal genetic legacy in modern dog populations. The closest detectable extant lineage to precontact American dogs is the canine transmissible venereal tumor, a contagious cancer clone derived from an individual dog that lived up to 8000 years ago.

A bioinformatic pipeline for identifying informative SNP panels for parentage assignment from RADseq data.

The development of high-throughput sequencing technologies is dramatically increasing the use of single nucleotide polymorphisms (SNPs) across the field of genetics, but most parentage studies of wild populations still rely on microsatellites. We developed a bioinformatic pipeline for identifying SNP panels that are informative for parentage analysis from restriction site-associated DNA sequencing (RADseq) data. This pipeline includes options for analysis with or without a reference genome, and provides methods to maximize genotyping accuracy and select sets of unlinked loci that have high statistical power. We test this pipeline on small populations of Mexican gray wolf and bighorn sheep, for which parentage analyses are expected to be challenging due to low genetic diversity and the presence of many closely related individuals. We compare the results of parentage analysis across SNP panels generated with or without the use of a reference genome, and between SNPs and microsatellites. For Mexican gray wolf, we conducted parentage analyses for 30 pups from a single cohort where samples were available from 64% of possible mothers and 53% of possible fathers, and the accuracy of parentage assignments could be estimated because true identities of parents were known a priori based on field data. For bighorn sheep, we conducted maternity analyses for 39 lambs from five cohorts where 77% of possible mothers were sampled, but true identities of parents were unknown. Analyses with and without a reference genome produced SNP panels with ≥95% parentage assignment accuracy for Mexican gray wolf, outperforming microsatellites at 78% accuracy. Maternity assignments were completely consistent across all SNP panels for the bighorn sheep, and were 74.4% consistent with assignments from microsatellites. Accuracy and consistency of parentage analysis were not reduced when using as few as 284 SNPs for Mexican gray wolf and 142 SNPs for bighorn sheep, indicating our pipeline can be used to develop SNP genotyping assays for parentage analysis with relatively small numbers of loci.

Genome-Wide Analysis of SNPs Is Consistent with No Domestic Dog Ancestry in the Endangered Mexican Wolf (Canis lupus baileyi).

The Mexican gray wolf (Canis lupus baileyi) was historically distributed throughout the southwestern United States and northern Mexico. Extensive predator removal campaigns during the early 20th century, however, resulted in its eventual extirpation by the mid 1980s. At this time, the Mexican wolf existed only in 3 separate captive lineages (McBride, Ghost Ranch, and Aragón) descended from 3, 2, and 2 founders, respectively. These lineages were merged in 1995 to increase the available genetic variation, and Mexican wolves were reintroduced into Arizona and New Mexico in 1998. Despite the ongoing management of the Mexican wolf population, it has been suggested that a proportion of the Mexican wolf ancestry may be recently derived from hybridization with domestic dogs. In this study, we genotyped 87 Mexican wolves, including individuals from all 3 captive lineages and cross-lineage wolves, for more than 172000 single nucleotide polymorphisms. We identified levels of genetic variation consistent with the pedigree record and effects of genetic rescue. To identify the potential to detect hybridization with domestic dogs, we compared our Mexican wolf genotypes with those from studies of domestic dogs and other gray wolves. The proportion of Mexican wolf ancestry assigned to domestic dogs was only between 0.06% (SD 0.23%) and 7.8% (SD 1.0%) for global and local ancestry estimates, respectively; and was consistent with simulated levels of incomplete lineage sorting. Overall, our results suggested that Mexican wolves lack biologically significant ancestry with dogs and have useful implications for the conservation and management of this endangered wolf subspecies.

Analysis of the canid Y-chromosome phylogeny using short-read sequencing data reveals the presence of distinct haplogroups among Neolithic European dogs.

BACKGROUND: Most genetic analyses of ancient and modern dogs have focused on variation in the autosomes or on the mitochondria. Mitochondrial DNA is more easily obtained from ancient samples than nuclear DNA and mitochondrial analyses have revealed important insights into the evolutionary history of canids. Utilizing a recently published dog Y-chromosome reference, we analyzed Y-chromosome sequence across a diverse collection of canids and determined the Y haplogroup of three ancient European dogs. RESULTS: We identified 1121 biallelic Y-chromosome SNVs using whole-genome sequences from 118 canids and defined variants diagnostic to distinct dog Y haplogroups. Similar to that of the mitochondria and previous more limited studies of Y diversity, we observe several deep splits in the Y-chromosome tree which may be the result of retained Y-chromosome diversity which predates dog domestication or post-domestication admixture with wolves. We find that Y-chromosomes from three ancient European dogs (4700-7000 years old) belong to distinct clades. CONCLUSIONS: We estimate that the time to the most recent comment ancestor of dog Y haplogroups is 68-151 thousand years ago. Analysis of three Y-chromosomes from the Neolithic confirms long stranding population structure among European dogs.

Postcranial diversity and recent ecomorphic impoverishment of North American gray wolves.

Recent advances in genomics and palaeontology have begun to unravel the complex evolutionary history of the gray wolf, Canis lupus Still, much of their phenotypic variation across time and space remains to be documented. We examined the limb morphology of the fossil and modern North American gray wolves from the late Quaternary (

Assessing SNP genotyping of noninvasively collected wildlife samples using microfluidic arrays.

Noninvasively collected samples are a common source of DNA in wildlife genetic studies. Currently, single nucleotide polymorphism (SNP) genotyping using microfluidic arrays is emerging as an easy-to-use and cost-effective methodology. Here we assessed the performance of microfluidic SNP arrays in genotyping noninvasive samples from grey wolves, European wildcats and brown bears, and we compared results with traditional microsatellite genotyping. We successfully SNP-genotyped 87%, 80% and 97% of the wolf, cat and bear samples, respectively. Genotype recovery was higher based on SNPs, while both marker types identified the same individuals and provided almost identical estimates of pairwise differentiation. We found that samples for which all SNP loci were scored had no disagreements across the three replicates (except one locus in a wolf sample). Thus, we argue that call rate (amplification success) can be used as a proxy for genotype quality, allowing the reduction of replication effort when call rate is high. Furthermore, we used cycle threshold values of real-time PCR to guide the choice of protocols for SNP amplification. Finally, we provide general guidelines for successful SNP genotyping of degraded DNA using microfluidic technology.

Population expansion and individual age affect endoparasite richness and diversity in a recolonising large carnivore population.

The recent recolonisation of the Central European lowland (CEL) by the grey wolf (Canis lupus) provides an excellent opportunity to study the effect of founder events on endoparasite diversity. Which role do prey and predator populations play in the re-establishment of endoparasite life cycles? Which intrinsic and extrinsic factors control individual endoparasite diversity in an expanding host population? In 53 individually known CEL wolves sampled in Germany, we revealed a community of four cestode, eight nematode, one trematode and 12 potential Sarcocystis species through molecular genetic techniques. Infections with zoonotic Echinococcus multilocularis, Trichinella britovi and T. spiralis occurred as single cases. Per capita endoparasite species richness and diversity significantly increased with population size and changed with age, whereas sex, microsatellite heterozygosity, and geographic origin had no effect. Tapeworm abundance (Taenia spp.) was significantly higher in immigrants than natives. Metacestode prevalence was slightly higher in ungulates from wolf territories than from control areas elsewhere. Even though alternative canid definitive hosts might also play a role within the investigated parasite life cycles, our findings indicate that (1) immigrated wolves increase parasite diversity in German packs, and (2) prevalence of wolf-associated parasites had declined during wolf absence and has now risen during recolonisation.

Ungulate predation and ecological roles of wolves and coyotes in eastern North America.

Understanding the ecological roles of species that influence ecosystem processes is a central goal of ecology and conservation biology. Eastern coyotes (Canis latrans) have ascended to the role of apex predator across much of eastern North America since the extirpation of wolves (Canis spp.) and there has been considerable confusion regarding their ability to prey on ungulates and their ecological niche relative to wolves. Eastern wolves (C. lycaon) are thought to have been the historical top predator in eastern deciduous forests and have previously been characterized as deer specialists that are inefficient predators of moose because of their smaller size relative to gray wolves (C. lupus). We investigated intrinsic and extrinsic influences on per capita kill rates of white-tailed deer (Odocoileus virginianus) and moose (Alces alces) during winter by sympatric packs of eastern coyotes, eastern wolves, and admixed canids in Ontario, Canada to clarify the predatory ability and ecological roles of the different canid top predators of eastern North America. Eastern coyote ancestry within packs negatively influenced per capita total ungulate (deer and moose combined) and moose kill rates. Furthermore, canids in packs dominated by eastern coyote ancestry consumed significantly less ungulate biomass and more anthropogenic food than packs dominated by wolf ancestry. Similar to gray wolves in previous studies, eastern wolves preyed on deer where they were available. However, in areas were deer were scarce, eastern wolves killed moose at rates similar to those previously documented for gray wolves at comparable moose densities across North America. Eastern coyotes are effective deer predators, but their dietary flexibility and low kill rates on moose suggest they have not replaced the ecological role of wolves in eastern North America.

Computed tomography examination and mitochondrial DNA analysis of Japanese wolf skull covered with skin.

A Canis skull, right half of the mandible and part of the left half of the mandible were subjected to three-dimensional (3D) computed tomography (CT) observation and mitochondrial DNA (mtDNA) analysis in order to determine whether the specimens belonged to the extinct Japanese wolf, Canis lupus hodophilax (Temminck, 1839). Osteometric analysis of the skull and right half of the mandible revealed that the material (JW275) was indeed typical of the Japanese wolf. Sequence analysis of a 600-bp mtDNA region revealed that the JW275 belonged to haplotype Group B, which is characterized by an 8-bp deletion in the mtDNA control region. The findings of this study suggest that 3D CT analysis is well suited to examining fragile and valuable biological samples, as it removes the need for destructive sampling.