Grey wolf genomic history reveals a dual ancestry of dogs

Rapid evolution of prehistoric dogs from wolves by natural and sexual selection emerges from an agent-based model

Wolves are among the earliest animals to be domesticated. However, the mechanism by which ancient wolves were domesticated into modern dogs is unknown. Most of the prevailing domestication hypotheses posit that humans selectively bred wolves that were more docile. However, a competing hypothesis states that wolves that were less hostile towards humans would essentially domesticate themselves by naturally selecting for tamer wolves since that would allow for easier access to food from human settlements. A major critique of the latter hypothesis is whether evolution by this natural selective pathway could have occurred in a sufficiently short time span. Simulating the process would help demonstrate if such an objection is sufficient to dismiss this hypothesis. Thus, we constructed an agent-based model of the evolution of a single trait, a measure of human tolerance, in canines to test the merit of the time constraint objection. We tested scenarios both with and without mate preference to provide a potential sexual selective force. We used fecundity and mortality rates from the literature for validation. Hartigan's dip test of unimodality was used to measure if and when divergence of populations occurred. Our results indicate that the proto-domestication hypothesis cannot be rejected on the basis of time constraints.

Ancient genomics and the origin, dispersal, and development of domestic sheep

The origins and prehistory of domestic sheep (Ovis aries) are incompletely understood; to address this, we generated data from 118 ancient genomes spanning 12,000 years sampled from across Eurasia. Genomes from Central Türkiye ~8000 BCE are genetically proximal to the domestic origins of sheep but do not fully explain the ancestry of later populations, suggesting a mosaic of wild ancestries. Genomic signatures indicate selection by ancient herders for pigmentation patterns, hornedness, and growth rate. Although the first European sheep flocks derive from Türkiye, in a notable parallel with ancient human genome discoveries, we detected a major influx of Western steppe-related ancestry in the Bronze Age.

[Thousands of years of human-dog relationship]

During recent years, much progress has been made in understanding the origin and evolution of the dog. Thanks to the collaboration between zooarchaeology, genomics and paleogenetics, researchers were able to hypothesize scenarios regarding the origins of the canine lineages present in Europe at the end of the Pleistocene and the beginning of the Holocene. Research has also shown a correlation between human and canine migration across time and space, highlighting a strong relationship between man and his best friend. This proximity between the two species is also illustrated by the adaptation of this species to anthropogenic selective pressures, particularly in parallel with cultural transitions. Although the history of this species still requires much exploration to be fully understood, these results provide new theoretical bases for understanding the interplay between humans and dogs.

Phenotypic and Complete Reference Whole Genome Sequence Analyses of Two Paenibacillus spp. Isolates from a Gray Wolf (Canis lupus) Gastrointestinal Tract

Inflammatory bowel disease (IBD) is increasing among mammals around the world, and domestic dogs are no exception. There is no approved cure for canine IBD with limited treatment options. Novel probiotic bacteria discovery from free-ranging animals for the treatment of IBD in domestic pets can likely yield promising probiotic candidates. Consequently, the overall aim was to isolate bacteria from free-ranging animals that could potentially be utilized as novel probiotics. Two bacteria identified as unique Paenibacillus spp. strains by small ribosomal RNA (16S) gene sequencing were isolated from the gastrointestinal tract of a North American Gray Wolf (Canis lupus). The bacteria were typed as Gram-variable, and both were catalase/oxidase positive as well as sensitive to commonly used antibiotics. The bacteria digested complex carbohydrates and lipids by standard assays. The isolated bacteria also inhibited the growth of Staphylococcus aureus and Micrococcus luteus. The whole genome sequence (WGS) length of bacterial isolate ClWae17B was 6,939,193 bp, while ClWae19 was 7,032,512 bp, both similar in size to other Paenibacillus spp. The genomes of both bacteria encoded enzymes involved with the metabolism of complex starches and lipids, such as lyases and pectinases, along with encoding antimicrobials such as lanthipeptides, lasso peptides, and cyclic-lactone-autoinducers. No pernicious virulence genes were identified in the WGS of either bacterial isolate. Phylogenetically, the most closely related bacteria based on 16S gene sequences and WGS were P. taichungensis for ClWae17B and P. amylolyticus for ClWae19. WGS analyses and phenotypic assays supported the hypothesis that the isolates described constitute two novel candidate probiotic bacteria for potential use in dogs.

A sedimentary ancient DNA perspective on human and carnivore persistence through the Late Pleistocene in El Mirón Cave, Spain

Caves are primary sites for studying human and animal subsistence patterns and genetic ancestry throughout the Palaeolithic. Iberia served as a critical human and animal refugium in Europe during the Last Glacial Maximum (LGM), 26.5 to 19 thousand years before the present (cal kya). Therefore, it is a key location for understanding human and animal population dynamics during this event. We recover and analyse sedimentary ancient DNA (sedaDNA) data from the lower archaeological stratigraphic sequence of El Mirón Cave (Cantabria, Spain), encompassing the (1) Late Mousterian period, associated with Neanderthals, and (2) the Gravettian (c. 31.5 cal kya), Solutrean (c. 24.5-22 cal kya), and Initial Magdalenian (d. 21-20.5 cal kya) periods, associated with anatomically modern humans. We identify 28 animal taxa including humans. Fifteen of these taxa had not been identified from the archaeozoological (i.e., faunal) record, including the presence of hyenas in the Magdalenian. Additionally, we provide phylogenetic analyses on 70 sedaDNA mtDNA genomes of fauna including the densest Iberian Pleistocene sampling of C. lupus. Finally, we recover three human mtDNA sequences from the Solutrean levels. These sequences, along with published data, suggest mtDNA haplogroup continuity in Iberia throughout the Solutrean/Last Glacial Maximum period.

GenomeDelta: detecting recent transposable element invasions without repeat library

We present GenomeDelta, a novel tool for identifying sample-specific sequences, such as recent transposable element (TE) invasions, without requiring a repeat library. GenomeDelta compares high-quality assemblies with short-read data to detect sequences absent from the short reads. It is applicable to both model and non-model organisms and can identify recent TE invasions, spatially heterogeneous sequences, viral insertions, and hotizontal gene transfers. GenomeDelta was validated with simulated and real data and used to discover three recent TE invasions in Drosophila melanogaster and a novel TE with geographic variation in Zymoseptoria tritici.

Late Pleistocene onset of mutualistic human/canid (Canis spp.) relationships in subarctic Alaska

Large canids (wolves, dogs, and coyote) and people form a close relationship in northern (subarctic and arctic) socioecological systems. Here, we document the antiquity of this bond and the multiple ways it manifested in interior Alaska, a region key to understanding the peopling of the Americas and early northern lifeways. We compile original and existing genomic, isotopic, and osteological canid data from archaeological, paleontological, and modern sites. Results show that in contrast to canids recovered in non-anthropic contexts, canids recovered in association with human occupations are markedly diverse. They include multiple species and intraspecific lineages, morphological variation, and diets ranging from terrestrial to marine. This variation is expressed along both geographic and temporal gradients, starting in the terminal Pleistocene with canids showing high marine dietary estimates. This paper provides evidence of the multiple ecological relationships between canids and people in the north-from predation, probable commensalism, and taming, to domestication-and of their early onset.

A map of canine sequence variation relative to a Greenland wolf outgroup

For over 15 years, canine genetics research relied on a reference assembly from a Boxer breed dog named Tasha (i.e., canFam3.1). Recent advances in long-read sequencing and genome assembly have led to the development of numerous high-quality assemblies from diverse canines. These assemblies represent notable improvements in completeness, contiguity, and the representation of gene promoters and gene models. Although genome graph and pan-genome approaches have promise, most genetic analyses in canines rely upon the mapping of Illumina sequencing reads to a single reference. The Dog10K consortium, and others, have generated deep catalogs of genetic variation through an alignment of Illumina sequencing reads to a reference genome obtained from a German Shepherd Dog named Mischka (i.e., canFam4, UU_Cfam_GSD_1.0). However, alignment to a breed-derived genome may introduce bias in genotype calling across samples. Since the use of an outgroup reference genome may remove this effect, we have reprocessed 1929 samples analyzed by the Dog10K consortium using a Greenland wolf (mCanLor1.2) as the reference. We efficiently performed remapping and variant calling using a GPU-implementation of common analysis tools. The resulting call set removes the variability in genetic differences seen across samples and breed relationships revealed by principal component analysis are not affected by the choice of reference genome. Using this sequence data, we inferred the history of population sizes and found that village dog populations experienced a 9-13 fold reduction in historic effective population size relative to wolves.

Epithelial barrier dysfunction and associated diseases in companion animals: Differences and similarities between humans and animals and research needs

Since the 1960s, more than 350,000 new chemicals have been introduced into the lives of humans and domestic animals. Many of them have become part of modern life and some are affecting nature as pollutants. Yet, our comprehension of their potential health risks for both humans and animals remains partial. The "epithelial barrier theory" suggests that genetic predisposition and exposure to diverse factors damaging the epithelial barriers contribute to the emergence of allergic and autoimmune conditions. Impaired epithelial barriers, microbial dysbiosis, and tissue inflammation have been observed in a high number of mucosal inflammatory, autoimmune and neuropsychiatric diseases, many of which showed increased prevalence in the last decades. Pets, especially cats and dogs, share living spaces with humans and are exposed to household cleaners, personal care products, air pollutants, and microplastics. The utilisation of cosmetic products and food additives for pets is on the rise, unfortunately, accompanied by less rigorous safety regulations than those governing human products. In this review, we explore the implications of disruptions in epithelial barriers on the well-being of companion animals, drawing comparisons with humans, and endeavour to elucidate the spectrum of diseases that afflict them. In addition, future research areas with the interconnectedness of human, animal, and environmental well-being are highlighted in line with the "One Health" concept.

Filtering out the noise: metagenomic classifiers optimize ancient DNA mapping

Contamination with exogenous DNA presents a significant challenge in ancient DNA (aDNA) studies of single organisms. Failure to address contamination from microbes, reagents, and present-day sources can impact the interpretation of results. Although field and laboratory protocols exist to limit contamination, there is still a need to accurately distinguish between endogenous and exogenous data computationally. Here, we propose a workflow to reduce exogenous contamination based on a metagenomic classifier. Unlike previous methods that relied exclusively on DNA sequencing reads mapping specificity to a single reference genome to remove contaminating reads, our approach uses Kraken2-based filtering before mapping to the reference genome. Using both simulated and empirical shotgun aDNA data, we show that this workflow presents a simple and efficient method that can be used in a wide range of computational environments-including personal machines. We propose strategies to build specific databases used to profile sequencing data that take into consideration available computational resources and prior knowledge about the target taxa and likely contaminants. Our workflow significantly reduces the overall computational resources required during the mapping process and reduces the total runtime by up to ~94%. The most significant impacts are observed in low endogenous samples. Importantly, contaminants that would map to the reference are filtered out using our strategy, reducing false positive alignments. We also show that our method results in a negligible loss of endogenous data with no measurable impact on downstream population genetics analyses.

Survey on the Occurrence of Zoonotic Bacterial Pathogens in the Feces of Wolves (Canis lupus italicus) Collected in a Protected Area in Central Italy

Previous investigations have explored the involvement of wolves in parasitic and viral diseases, but data on the zoonotic bacteria are limited. The aim of this study was to assess the occurrence of bacterial zoonotic agents in 16 wolf (Canis lupus italicus) fecal samples collected in a protected area in Central Italy. Campylobacter spp., Salmonella spp., Yersinia spp., Listeria monocytogenes, and Shiga Toxin-Producing Escherichia coli (STEC) were investigated by culture, while polymerase chain reaction (PCR) was employed to detect Coxiella burnetii, Mycobacterium spp., Brucella spp., and Francisella tularensis. The presence of Extended Spectrum β-Lactamase (ESBL)- and carbapenemase-producing Enterobacteriaceae was also evaluated, using selective isolation media and detection of antimicrobial resistance genes. All samples were negative for Campylobacter spp., Salmonella spp., C. burnetii, Mycobacterium spp., Brucella spp., F. tularensis, and carbapenemase-producing Enterobacteriaceae. One sample tested positive for Yersinia aldovae and three for Yersinia enterocolitica BT1A. One L. monocytogenes (serogroup IIa) and one STEC, carrying the stx1 gene, were isolated. Two ESBL isolates were detected: one Serratia fonticola, carrying blaFONA-3/6 gene, and one Escherichia coli, carrying blaCTX-M-1 gene. Both ESBL isolates were resistant to different antimicrobials and therefore classified as multi-drug-resistant. Our data suggest that wolves are potential carriers of zoonotic bacteria and may contribute to the environmental contamination through their feces.

Breed function and behaviour correlate with endocranial volume in domestic dogs

Domestic dogs display a remarkable diversity of functions, morphologies and cognitive abilities. Using data from 1682 dogs representing 172 breeds, we tested for variation in relative endocranial volume (REV), a proxy for brain size and a basic measure of cognitive ability, in relation to breed function, phylogenetic classification, cranial shape, cooperative behaviour and temperament. Function, body size, phylogenetic clade and cranial shape correlate with REV. Toy dogs, functioning mainly as companions, have the largest endocranial volumes relative to their body size. Working dogs, bred to perform complex human-assistance skills and reportedly possessing higher cognitive abilities, have the smallest. Our results thus show that complex skills and cooperative behaviour-a hallmark of social cognition-do not predict larger REV in dogs. However, REV increases with fear and aggression, attention-seeking and separation anxiety and decreases with trainability. Significant correlations between REV and behavioural traits underscore the evolutionary plasticity of mammalian brain size under domestication and artificial selection and provide support for hypotheses linking the modulation of fear and aggression to brain size change under domestication.

Coyotes can do 'puppy dog eyes' too: comparing interspecific variation in Canis facial expression muscles

Facial expressions are critical for non-verbal communication. The Canis genus epitomizes the interplay between behaviour and morphology in the evolution of non-verbal communication. Recent work suggests that the levator anguli oculi medialis (LAOM) muscle is unique to dogs (Canis familiaris) within the Canis genus and evolved due to domestication. The LAOM raises the inner eyebrows, resulting in the 'puppy dog eyes' expression. Here, we test whether the LAOM is a derived trait in dogs by (i) examining the facial expression muscles of a closely related and ancestral wild Canis species, the coyote (C. latrans) and (ii) comparing our results with other Canis and canid taxa. We discover that coyotes have a well-developed LAOM like dogs, which differs from the modified/absent LAOM in grey wolves. Our findings challenge the hypothesis that the LAOM developed due to domestication. We suggest that the LAOM is a basal trait that was lost in grey wolves. Additionally, we find inter- and intraspecific variations in the size of the muscles of the outer ear, forehead, lips and rostrum, indicating potential adaptations related to sensory perception, communication and individual-level functional variations within canids. Together, this research expands our knowledge of facial expressions, their evolution and their role in communication.

Estimation of spatial demographic maps from polymorphism data using a neural network

A fundamental goal in population genetics is to understand how variation is arrayed over natural landscapes. From first principles we know that common features such as heterogeneous population densities and barriers to dispersal should shape genetic variation over space, however there are few tools currently available that can deal with these ubiquitous complexities. Geographically referenced single nucleotide polymorphism (SNP) data are increasingly accessible, presenting an opportunity to study genetic variation across geographic space in myriad species. We present a new inference method that uses geo-referenced SNPs and a deep neural network to estimate spatially heterogeneous maps of population density and dispersal rate. Our neural network trains on simulated input and output pairings, where the input consists of genotypes and sampling locations generated from a continuous space population genetic simulator, and the output is a map of the true demographic parameters. We benchmark our tool against existing methods and discuss qualitative differences between the different approaches; in particular, our program is unique because it infers the magnitude of both dispersal and density as well as their variation over the landscape, and it does so using SNP data. Similar methods are constrained to estimating relative migration rates, or require identity-by-descent blocks as input. We applied our tool to empirical data from North American grey wolves, for which it estimated mostly reasonable demographic parameters, but was affected by incomplete spatial sampling. Genetic based methods like ours complement other, direct methods for estimating past and present demography, and we believe will serve as valuable tools for applications in conservation, ecology and evolutionary biology. An open source software package implementing our method is available from https://github.com/kr-colab/mapNN.

Testing times: disentangling admixture histories in recent and complex demographies using ancient DNA

Our knowledge of human evolutionary history has been greatly advanced by paleogenomics. Since the 2020s, the study of ancient DNA has increasingly focused on reconstructing the recent past. However, the accuracy of paleogenomic methods in resolving questions of historical and archaeological importance amidst the increased demographic complexity and decreased genetic differentiation remains an open question. We evaluated the performance and behavior of two commonly used methods, qpAdm and the f3-statistic, on admixture inference under a diversity of demographic models and data conditions. We performed two complementary simulation approaches-firstly exploring a wide demographic parameter space under four simple demographic models of varying complexities and configurations using branch-length data from two chromosomes-and secondly, we analyzed a model of Eurasian history composed of 59 populations using whole-genome data modified with ancient DNA conditions such as SNP ascertainment, data missingness, and pseudohaploidization. We observe that population differentiation is the primary factor driving qpAdm performance. Notably, while complex gene flow histories influence which models are classified as plausible, they do not reduce overall performance. Under conditions reflective of the historical period, qpAdm most frequently identifies the true model as plausible among a small candidate set of closely related populations. To increase the utility for resolving fine-scaled hypotheses, we provide a heuristic for further distinguishing between candidate models that incorporates qpAdm model P-values and f3-statistics. Finally, we demonstrate a significant performance increase for qpAdm using whole-genome branch-length f2-statistics, highlighting the potential for improved demographic inference that could be achieved with future advancements in f-statistic estimations.

Performance of qpAdm -based screens for genetic admixture on admixture-graph-shaped histories and stepping-stone landscapes

qpAdm is a statistical tool that is often used in exploratory archaeogenetic studies for finding optimal admixture models of population history. Despite its popularity, qpAdm remains untested on histories in the form of admixture graphs of random topology or stepping-stone landscapes. We analyzed data from such simulations and found that while for admixture-graph-shaped histories there exist simple solutions (temporal stratification) for minimizing false findings of gene flow, in the case of stepping-stone landscapes the method generates results that do not appear suspect but are misleading: feasible qpAdm models are either accurate but simplistic in the context of landscapes, or highly inaccurate in the case of multi-component models. This is largely is due to two reasons: 1) because of complex migration networks that violate the assumptions of the method, there is poor correlation between qpAdm p -values and model optimality in many sections of the parameter space; 2) admixture fraction estimates between 0 and 1 are largely restricted to symmetric source configurations around targets, hence popular [0, 1] model plausibility criteria confound analyses of landscape-type demographies, unless their interpretations are explicitly spatial. For many species/regions/periods archaeogenetic sampling is very sparse and may be random with respect to population density of ancient individuals. In this situation only a specific combination of landscape properties and feasibility criteria allows to efficiently reject highly asymmetric non-optimal models most abundant in random deme sets. This problem may obscure useful signal (rare optimal models) and might be responsible for some claims about rapid long-distance migrations in the literature.

Comparative Population Genomics of Arctic Sled Dogs Reveals a Deep and Complex History

Recent evidence demonstrates genomic and morphological continuity in the Arctic ancestral lineage of dogs. Here, we use the Siberian Husky to investigate the genomic legacy of the northeast Eurasian Arctic lineage and model the deep population history using genome-wide single nucleotide polymorphisms. Utilizing ancient dog-calibrated molecular clocks, we found that at least two distinct lineages of Arctic dogs existed in ancient Eurasia at the end of the Pleistocene. This pushes back the origin of sled dogs in the northeast Siberian Arctic with humans likely intentionally selecting dogs to perform different functions and keeping breeding populations that overlap in time and space relatively reproductively isolated. In modern Siberian Huskies, we found significant population structure based on how they are used by humans, recent European breed introgression in about half of the dogs that participate in races, moderate levels of inbreeding, and fewer potentially harmful variants in populations under strong selection for form and function (show, sled show, and racing populations of Siberian Huskies). As the struggle to preserve unique evolutionary lineages while maintaining genetic health intensifies across pedigreed dogs, understanding the genomic history to guide policies and best practices for breed management is crucial to sustain these ancient lineages and their unique evolutionary identity.

Neolithic Yersinia pestis infections in humans and a dog

Yersinia pestis has been infecting humans since the Late Neolithic (LN). Whether those early infections were isolated zoonoses or initiators of a pandemic remains unclear. We report Y. pestis infections in two individuals (of 133) from the LN necropolis at Warburg (Germany, 5300-4900 cal BP). Our analyses show that the two genomes belong to distinct strains and reflect independent infection events. All LN genomes known today (n = 4) are basal in the phylogeny and represent separate lineages that probably originated in different animal hosts. In the LN, an opening of the landscape resulted in the introduction of new rodent species, which may have acted as Y. pestis reservoirs. Coincidentally, the number of dogs increased, possibly leading to Y. pestis infections in canines. Indeed, we detect Y. pestis in an LN dog. Collectively, our data suggest that Y. pestis frequently entered human settlements at the time without causing significant outbreaks.

Multiple ancestries and shared gene flow among modern livestock guarding dogs

Livestock guarding dogs (LGDs) have been used to protect livestock for millennia. While previous works suggested a single origin of modern LGDs, the degree and source of shared ancestry have not been tested. To address this, we generated genome-wide SNP data from 304 LGDs and combined it with public genomic data from 2,183 modern and 22 ancient dogs. Our findings reveal shared ancestry and extensive gene flow among modern LGD breeds which we attribute to historical livestock migrations. Additionally, admixture between LGDs and free-ranging dogs argues against reproductive isolation as a core mechanism for maintaining the specialized skills of LGDs. Finally, we identify two lineages within modern LGDs and uncover multiple ancestries tracing back to distinct Eurasian ancient dogs, concordant with the absence of a single ancestor. Overall, our work explores the complex evolutionary history of LGDs, offering valuable insights into how human and livestock co-migrations shaped this functional group.

Large-scale genomic analysis of the domestic dog informs biological discovery

Recent advances in genomics, coupled with a unique population structure and remarkable levels of variation, have propelled the domestic dog to new levels as a system for understanding fundamental principles in mammalian biology. Central to this advance are more than 350 recognized breeds, each a closed population that has undergone selection for unique features. Genetic variation in the domestic dog is particularly well characterized compared with other domestic mammals, with almost 3000 high-coverage genomes publicly available. Importantly, as the number of sequenced genomes increases, new avenues for analysis are becoming available. Herein, we discuss recent discoveries in canine genomics regarding behavior, morphology, and disease susceptibility. We explore the limitations of current data sets for variant interpretation, tradeoffs between sequencing strategies, and the burgeoning role of long-read genomes for capturing structural variants. In addition, we consider how large-scale collections of whole-genome sequence data drive rare variant discovery and assess the geographic distribution of canine diversity, which identifies Asia as a major source of missing variation. Finally, we review recent comparative genomic analyses that will facilitate annotation of the noncoding genome in dogs.

Ancient genomes reveal over two thousand years of dingo population structure

Dingoes are culturally and ecologically important free-living canids whose ancestors arrived in Australia over 3,000 B.P., likely transported by seafaring people. However, the early history of dingoes in Australia-including the number of founding populations and their routes of introduction-remains uncertain. This uncertainty arises partly from the complex and poorly understood relationship between modern dingoes and New Guinea singing dogs, and suspicions that post-Colonial hybridization has introduced recent domestic dog ancestry into the genomes of many wild dingo populations. In this study, we analyzed genome-wide data from nine ancient dingo specimens ranging in age from 400 to 2,746 y old, predating the introduction of domestic dogs to Australia by European colonists. We uncovered evidence that the continent-wide population structure observed in modern dingo populations had already emerged several thousand years ago. We also detected excess allele sharing between New Guinea singing dogs and ancient dingoes from coastal New South Wales (NSW) compared to ancient dingoes from southern Australia, irrespective of any post-Colonial hybrid ancestry in the genomes of modern individuals. Our results are consistent with several demographic scenarios, including a scenario where the ancestry of dingoes from the east coast of Australia results from at least two waves of migration from source populations with varying affinities to New Guinea singing dogs. We also contribute to the growing body of evidence that modern dingoes derive little genomic ancestry from post-Colonial hybridization with other domestic dog lineages, instead descending primarily from ancient canids introduced to Sahul thousands of years ago.

Estimation of spatial demographic maps from polymorphism data using a neural network

A fundamental goal in population genetics is to understand how variation is arrayed over natural landscapes. From first principles we know that common features such as heterogeneous population densities and barriers to dispersal should shape genetic variation over space, however there are few tools currently available that can deal with these ubiquitous complexities. Geographically referenced single nucleotide polymorphism (SNP) data are increasingly accessible, presenting an opportunity to study genetic variation across geographic space in myriad species. We present a new inference method that uses geo-referenced SNPs and a deep neural network to estimate spatially heterogeneous maps of population density and dispersal rate. Our neural network trains on simulated input and output pairings, where the input consists of genotypes and sampling locations generated from a continuous space population genetic simulator, and the output is a map of the true demographic parameters. We benchmark our tool against existing methods and discuss qualitative differences between the different approaches; in particular, our program is unique because it infers the magnitude of both dispersal and density as well as their variation over the landscape, and it does so using SNP data. Similar methods are constrained to estimating relative migration rates, or require identity by descent blocks as input. We applied our tool to empirical data from North American grey wolves, for which it estimated mostly reasonable demographic parameters, but was affected by incomplete spatial sampling. Genetic based methods like ours complement other, direct methods for estimating past and present demography, and we believe will serve as valuable tools for applications in conservation, ecology, and evolutionary biology. An open source software package implementing our method is available from https://github.com/kr-colab/mapNN .

Genomic analysis of wolves from Pakistan clarifies boundaries among three divergent wolf lineages

Among the three main divergent lineages of gray wolf (Canis lupus), the Holarctic lineage is the most widespread and best studied, particularly in North America and Europe. Less is known about Tibetan (also called Himalayan) and Indian wolf lineages in southern Asia, especially in areas surrounding Pakistan where all three lineages are thought to meet. Given the endangered status of the Indian wolf in neighboring India and unclear southwestern boundary of the Tibetan wolf range, we conducted mitochondrial and genome-wide sequencing of wolves from Pakistan and Kyrgyzstan. Sequences of the mitochondrial D-loop region of 81 wolves from Pakistan indicated contact zones between Holarctic and Indian lineages across the northern and western mountains of Pakistan. Reduced-representation genome sequencing of eight wolves indicated an east-to-west cline of Indian to Holarctic ancestry, consistent with a contact zone between these two lineages in Pakistan. The western boundary of the Tibetan lineage corresponded to the Ladakh region of India's Himalayas with a narrow zone of admixture spanning this boundary from the Karakoram Mountains of northern Pakistan into Ladakh, India. Our results highlight the conservation significance of Pakistan's wolf populations, especially the remaining populations in Sindh and Southern Punjab that represent the highly endangered Indian lineage.

Improving data archiving practices in ancient genomics

Ancient DNA is producing a rich record of past genetic diversity in humans and other species. However, unless the primary data is appropriately archived, its long-term value will not be fully realised. I surveyed publicly archived data from 42 recent ancient genomics studies. Half of the studies archived incomplete datasets, preventing accurate replication and representing a loss of data of potential future use. No studies met all criteria that could be considered best practice. Based on these results, I make six recommendations for data producers: (1) archive all sequencing reads, not just those that aligned to a reference genome, (2) archive read alignments too, but as secondary analysis files, (3) provide correct experiment metadata on samples, libraries and sequencing runs, (4) provide informative sample metadata, (5) archive data from low-coverage and negative experiments, and (6) document archiving choices in papers, and peer review these. Given the reliance on destructive sampling of finite material, ancient genomics studies have a particularly strong responsibility to ensure the longevity and reusability of generated data.

A comparative analysis reveals the genomic diversity among 8 Muscovy duck populations

The Muscovy duck (Cairina moschata) is a waterfowl indigenous to the neotropical regions of Central and South America. It has low demand for concentrated feed and strong adaptability to different rearing conditions. After introduced to China through Eurasian commercial trade, Muscovy ducks have a domestication history of around 300 years in the Fujian Province of China. In the 1990s, the commodity Muscovy duck breed "Crimo," cultivated in Europe, entered the Chinese market for consumption and breeding purposes. Due to the different selective breeding processes, Muscovy ducks have various populational traits and lack transparency of their genetic background. To remove this burden in the Muscovy duck breeding process, we analyzed genomic data from 8 populations totaling 83 individuals. We identify 11.24 million single nucleotide polymorphisms (SNPs) and categorized these individuals into the Fujian-bred and the Crimo populations according to phylogenetic analyses. We then delved deeper into their evolutionary relationships through assessing population structure, calculating fixation index (FST) values, and measuring genetic distances. Our exploration of runs of homozygosity (ROHs) and homozygous-by-descent (HBD) uncovered genomic regions enriched for genes implicated in fatty acid metabolism, development, and immunity pathways. Selective sweep analyses further indicated strong selective pressures exerted on genes including TECR, STAT2, and TRAF5. These findings provide insights into genetic variations of Muscovy ducks, thus offering valuable information regarding genetic diversity, population conservation, and genome associated with the breeding of Muscovy ducks.

Dogs and their genes: what ever will they think of next?

The Edward Novitski Prize recognizes creativity and intellectual ingenuity in the solution of problems in genetics research. The prize honors scientific experimental work-either a single experimental accomplishment or a body of work. Ostrander is recognized for work developing the domestic dog as an experimental system for solving fundamental biological problems and identifying genetic sequences of relevance to human health and disease. Including work on disease and behavioral health, Ostrander has shown a dedication to creative methods for understanding canine genetics and the value of translating research organisms to human genetics.

Historic dog Furs Unravel the Origin and Artificial Selection of Modern Nordic Lapphund and Elkhound dog Breeds

The origins and extreme morphological evolution of the modern dog breeds are poorly studied because the founder populations are extinct. Here, we analyse eight 100 to 200 years old dog fur samples obtained from traditional North Swedish clothing, to explore the origin and artificial selection of the modern Nordic Lapphund and Elkhound dog breeds. Population genomic analysis confirmed the Lapphund and Elkhound breeds to originate from the local dog population, and showed a distinct decrease in genetic diversity in agreement with intense breeding. We identified eleven genes under positive selection during the breed development. In particular, the MSRB3 gene, associated with breed-related ear morphology, was selected in all Lapphund and Elkhound breeds, and functional assays showed that a SNP mutation in the 3'UTR region suppresses its expression through miRNA regulation. Our findings demonstrate analysis of near-modern dog artifacts as an effective tool for interpreting the origin and artificial selection of the modern dog breeds.

Protein folding, cellular stress and cancer

Proteins are acknowledged as the phenotypical manifestation of the genotype, because protein-coding genes carry the information for the strings of amino acids that constitute the proteins. It is widely accepted that protein function depends on the corresponding "native" structure or folding achieved within the cell, and that native protein folding corresponds to the lowest free energy minimum for a given protein. However, protein folding within the cell is a non-deterministic dissipative process that from the same input may produce different outcomes, thus conformational heterogeneity of folded proteins is the rule and not the exception. Local changes in the intracellular environment promote variation in protein folding. Hence protein folding requires "supervision" by a host of chaperones and co-chaperones that help their client proteins to achieve the folding that is most stable according to the local environment. Such environmental influence on protein folding is continuously transduced with the help of the cellular stress responses (CSRs) and this may lead to changes in the rules of engagement between proteins, so that the corresponding protein interactome could be modified by the environment leading to an alternative cellular phenotype. This allows for a phenotypic plasticity useful for adapting to sudden and/or transient environmental changes at the cellular level. Starting from this perspective, hereunder we develop the argument that the presence of sustained cellular stress coupled to efficient CSRs may lead to the selection of an aberrant phenotype as the resulting adaptation of the cellular proteome (and the corresponding interactome) to such stressful conditions, and this can be a common epigenetic pathway to cancer.

Archaeology meets environmental genomics: implementing sedaDNA in the study of the human past

Sedimentary ancient DNA (sedaDNA) has become one of the standard applications in the field of paleogenomics in recent years. It has been used for paleoenvironmental reconstructions, detecting the presence of prehistoric species in the absence of macro remains and even investigating the evolutionary history of a few species. However, its application in archaeology has been limited and primarily focused on humans. This article argues that sedaDNA holds significant potential in addressing key archaeological questions concerning the origins, lifestyles, and environments of past human populations. Our aim is to facilitate the integration of sedaDNA into the standard workflows in archaeology as a transformative tool, thereby unleashing its full potential for studying the human past. Ultimately, we not only underscore the challenges inherent in the sedaDNA field but also provide a research agenda for essential enhancements needed for implementing sedaDNA into the archaeological workflow.

Cervidae antlers exploited to manufacture prehistoric tools and hunting implements as a reliable source of ancient DNA

Antler is one of the primary animal raw materials exploited for technical purposes by the hunter-gatherer groups of the Eurasian Upper Palaeolithic (UP) all over the ecological range of deers, and beyond. It was exhaustively employed to produce one of the most critical tools for the survival of the UP societies: hunting weapons. However, antler implements can be made from diverse deer taxa, with different ecological requirements and ethological behaviours. Identifying the antler's origin at a taxonomic level is thus essential in improving our knowledge of humans' functional, practical and symbolic choices, as well as the human-animal interface during Prehistoric times. Nevertheless, palaeogenetics analyses have focused mainly on bone and teeth, with genetic studies of antler generally focused on modern deer conservation. Here we present the results of the first whole mitochondrial genome ancient DNA (aDNA) analysis by means of in-solution hybridisation capture of antlers from pre-Holocene archaeological contexts. We analysed a set of 50 Palaeolithic and Neolithic (c. 34-8ka) antler and osseous objects from South-Western Europe, Central Europe, South-Western Asia and the Caucasus. We successfully obtained aDNA, allowing us to identify the exploited taxa and demonstrate the archaeological relevance of those finds. Moreover, as most of the antlers were sampled using a minimally-invasive method, further analyses (morphometric, technical, genetic, radiometric and more) remain possible on these objects.

Ancient Mitogenomes Reveal the Maternal Genetic History of East Asian Dogs

Recent studies have suggested that dogs were domesticated during the Last Glacial Maximum (LGM) in Siberia, which contrasts with previous proposed domestication centers (e.g. Europe, the Middle East, and East Asia). Ancient DNA provides a powerful resource for the study of mammalian evolution and has been widely used to understand the genetic history of domestic animals. To understand the maternal genetic history of East Asian dogs, we have made a complete mitogenome dataset of 120 East Asian canids from 38 archaeological sites, including 102 newly sequenced from 12.9 to 1 ka BP (1,000 years before present). The majority (112/119, 94.12%) belonged to haplogroup A, and half of these (55/112, 49.11%) belonged to sub-haplogroup A1b. Most existing mitochondrial haplogroups were present in ancient East Asian dogs. However, mitochondrial lineages in ancient northern dogs (northeastern Eurasia and northern East Asia) were deeper and older than those in southern East Asian dogs. Results suggests that East Asian dogs originated from northeastern Eurasian populations after the LGM, dispersing in two possible directions after domestication. Western Eurasian (Europe and the Middle East) dog maternal ancestries genetically influenced East Asian dogs from approximately 4 ka BP, dramatically increasing after 3 ka BP, and afterwards largely replaced most primary maternal lineages in northern East Asia. Additionally, at least three major mitogenome sub-haplogroups of haplogroup A (A1a, A1b, and A3) reveal at least two major dispersal waves onto the Qinghai-Tibet Plateau in ancient times, indicating eastern (A1b and A3) and western (A1a) Eurasian origins.

From wolves to humans: oral microbiome resistance to transfer across mammalian hosts

The mammalian mouth is colonized by complex microbial communities, adapted to specific niches, and in homeostasis with the host. Individual microbes interact metabolically and rely primarily on nutrients provided by the host, with which they have potentially co-evolved along the mammalian lineages. The oral environment is similar across mammals, but the diversity, specificity, and evolution of community structure in related or interacting mammals are little understood. Here, we compared the oral microbiomes of dogs with those of wild wolves and humans. In dogs, we found an increased microbial diversity relative to wolves, possibly related to the transition to omnivorous nutrition following domestication. This includes a larger diversity of Patescibacteria than previously reported in any other oral microbiota. The oral microbes are most distinct at bacterial species or strain levels, with few if any shared between humans and canids, while the close evolutionary relationship between wolves and dogs is reflected by numerous shared taxa. More taxa are shared at higher taxonomic levels including with humans, supporting their more ancestral common mammalian colonization followed by diversification. Phylogenies of selected oral bacterial lineages do not support stable human-dog microbial transfers but suggest diversification along mammalian lineages (apes and canids). Therefore, despite millennia of cohabitation and close interaction, the host and its native community controls and limits the assimilation of new microbes, even if closely related. Higher resolution metagenomic and microbial physiological studies, covering a larger mammalian diversity, should help understand how oral communities assemble, adapt, and interact with their hosts.IMPORTANCENumerous types of microbes colonize the mouth after birth and play important roles in maintaining oral health. When the microbiota-host homeostasis is perturbed, proliferation of some bacteria leads to diseases such as caries and periodontitis. Unlike the gut microbiome, the diversity of oral microbes across the mammalian evolutionary space is not understood. Our study compared the oral microbiomes of wild wolves, dogs, and apes (humans, chimpanzees, and bonobos), with the aim of identifying if microbes have been potentially exchanged between humans and dogs as a result of domestication and cohabitation. We found little if any evidence for such exchanges. The significance of our research is in finding that the oral microbiota and/or the host limit the acquisition of exogenous microbes, which is important in the context of natural exclusion of potential novel pathogens. We provide a framework for expanded higher-resolution studies across domestic and wild animals to understand resistance/resilience.

Japanese wolves are most closely related to dogs and share DNA with East Eurasian dogs

Although the domestic dog's origin is still unclear, this lineage is believed to have been domesticated from an extinct population of gray wolves, which is expected to be more closely related to dogs than to other populations of gray wolves. Here, we sequence the whole genomes of nine Japanese wolves (7.5-100x: Edo to Meiji periods) and 11 modern Japanese dogs and analyze them together with those from other populations of dogs and wolves. A phylogenomic tree shows that, among the gray wolves, Japanese wolves are closest to the dog, suggesting that the ancestor of dogs is closely related to the ancestor of the Japanese wolf. Based on phylogenetic and geographic relationships, the dog lineage has most likely originated in East Asia, where it diverged from a common ancestor with the Japanese wolf. Since East Eurasian dogs possess Japanese wolf ancestry, we estimate an introgression event from the ancestor of the Japanese wolf to the ancestor of the East Eurasian dog that occurred before the dog's arrival in the Japanese archipelago.

Behavioural Synchronisation between Dogs and Humans: Unveiling Interspecific Motor Resonance?

Dogs' behavioural synchronisation with humans is of growing scientific interest. However, studies lack a comprehensive exploration of the neurocognitive foundations of this social cognitive ability. Drawing parallels from the mechanisms underlying behavioural synchronisation in humans, specifically motor resonance and the recruitment of mirror neurons, we hypothesise that dogs' behavioural synchronisation with humans is underpinned by a similar mechanism, namely interspecific motor resonance. Based on a literature review, we argue that dogs possess the prerequisites for motor resonance, and we suggest that interspecific behavioural synchronisation relies on the activation of both human and canine mirror neurons. Furthermore, interspecific behavioural studies highlight certain characteristics of motor resonance, including motor contagion and its social modulators. While these findings strongly suggest the potential existence of interspecific motor resonance, direct proof remains to be established. Our analysis thus paves the way for future research to confirm the existence of interspecific motor resonance as the neurocognitive foundation for interspecific behavioural synchronisation. Unravelling the neurocognitive mechanisms underlying this behavioural adjustment holds profound implications for understanding the evolutionary dynamics of dogs alongside humans and improving the day-to-day management of dog-human interactions.

Population structure and adaptability analysis of Schizothorax o'connori based on whole-genome resequencing

BACKGROUND

Schizothorax o'connori is an endemic fish distributed in the upper and lower reaches of the Yarlung Zangbo River in China. It has experienced a fourth round of whole gene replication events and is a good model for exploring the genetic differentiation and environmental adaptability of fish in the Qinghai-Tibet Plateau. The uplift of the Qinghai-Tibet Plateau has led to changes in the river system, thereby affecting gene exchange and population differentiation between fish populations. With the release of fish whole genome data, whole genome resequencing has been widely used in genetic evolutionary analysis and screening of selected genes in fish, which can better elucidate the genetic basis and molecular environmental adaptation mechanisms of fish. Therefore, our purpose of this study was to understand the population structure and adaptive characteristics of S. o'connori using the whole-genome resequencing method.

RESULTS

The results showed that 23,602,746 SNPs were identified from seven populations, mostly distributed on chromosomes 2 and 23. There was no significant genetic differentiation between the populations, and the genetic diversity was relatively low. However, the Zangga population could be separated from the Bomi, Linzhi, and Milin populations in the cluster analysis. Based on historical dynamics analysis of the population, the size of the ancestral population of S. o'connori was affected by the late accelerated uplift of the Qinghai Tibet Plateau and the Fourth Glacial Age. The selected sites were mostly enriched in pathways related to DNA repair and energy metabolism.

CONCLUSION

Overall, the whole-genome resequencing analysis provides valuable insights into the population structure and adaptive characteristics of S. o'connori. There was no obvious genetic differentiation at the genome level between the S. o'connori populations upstream and downstream of the Yarlung Zangbo River. The current distribution pattern and genetic diversity are influenced by the late accelerated uplift of the Qinghai Tibet Plateau and the Fourth Ice Age. The selected sites of S. o'connori are enriched in the energy metabolism and DNA repair pathways to adapt to the low temperature and strong ultraviolet radiation environment at high altitude.

The history of Coast Salish "woolly dogs" revealed by ancient genomics and Indigenous Knowledge

Ancestral Coast Salish societies in the Pacific Northwest kept long-haired "woolly dogs" that were bred and cared for over millennia. However, the dog wool-weaving tradition declined during the 19th century, and the population was lost. In this study, we analyzed genomic and isotopic data from a preserved woolly dog pelt from "Mutton," collected in 1859. Mutton is the only known example of an Indigenous North American dog with dominant precolonial ancestry postdating the onset of settler colonialism. We identified candidate genetic variants potentially linked with their distinct woolly phenotype. We integrated these data with interviews from Coast Salish Elders, Knowledge Keepers, and weavers about shared traditional knowledge and memories surrounding woolly dogs, their importance within Coast Salish societies, and how colonial policies led directly to their disappearance.

Anthropogenic Change and the Process of Speciation

Anthropogenic impacts on the environment alter speciation processes by affecting both geographical contexts and selection patterns on a worldwide scale. Here we review evidence of these effects. We find that human activities often generate spatial isolation between populations and thereby promote genetic divergence but also frequently cause sudden secondary contact and hybridization between diverging lineages. Human-caused environmental changes produce new ecological niches, altering selection in diverse ways that can drive diversification; but changes also often remove niches and cause extirpations. Human impacts that alter selection regimes are widespread and strong in magnitude, ranging from local changes in biotic and abiotic conditions to direct harvesting to global climate change. Altered selection, and evolutionary responses to it, impacts early-stage divergence of lineages, but does not necessarily lead toward speciation and persistence of separate species. Altogether, humans both promote and hinder speciation, although new species would form very slowly relative to anthropogenic hybridization, which can be nearly instantaneous. Speculating about the future of speciation, we highlight two key conclusions: (1) Humans will have a large influence on extinction and "despeciation" dynamics in the short term and on early-stage lineage divergence, and thus potentially speciation in the longer term, and (2) long-term monitoring combined with easily dated anthropogenic changes will improve our understanding of the processes of speciation. We can use this knowledge to preserve and restore ecosystems in ways that promote (re-)diversification, increasing future opportunities of speciation and enhancing biodiversity.

Testing Times: Challenges in Disentangling Admixture Histories in Recent and Complex Demographies

Paleogenomics has expanded our knowledge of human evolutionary history. Since the 2020s, the study of ancient DNA has increased its focus on reconstructing the recent past. However, the accuracy of paleogenomic methods in answering questions of historical and archaeological importance amidst the increased demographic complexity and decreased genetic differentiation within the historical period remains an open question. We used two simulation approaches to evaluate the limitations and behavior of commonly used methods, qpAdm and the f3-statistic, on admixture inference. The first is based on branch-length data simulated from four simple demographic models of varying complexities and configurations. The second, an analysis of Eurasian history composed of 59 populations using whole-genome data modified with ancient DNA conditions such as SNP ascertainment, data missingness, and pseudo-haploidization. We show that under conditions resembling historical populations, qpAdm can identify a small candidate set of true sources and populations closely related to them. However, in typical ancient DNA conditions, qpAdm is unable to further distinguish between them, limiting its utility for resolving fine-scaled hypotheses. Notably, we find that complex gene-flow histories generally lead to improvements in the performance of qpAdm and observe no bias in the estimation of admixture weights. We offer a heuristic for admixture inference that incorporates admixture weight estimate and P-values of qpAdm models, and f3-statistics to enhance the power to distinguish between multiple plausible candidates. Finally, we highlight the future potential of qpAdm through whole-genome branch-length f2-statistics, demonstrating the improved demographic inference that could be achieved with advancements in f-statistic estimations.

Limited historical admixture between European wildcats and domestic cats

Domestic cats were derived from the Near Eastern wildcat (Felis lybica), after which they dispersed with people into Europe. As they did so, it is possible that they interbred with the indigenous population of European wildcats (Felis silvestris). Gene flow between incoming domestic animals and closely related indigenous wild species has been previously demonstrated in other taxa, including pigs, sheep, goats, bees, chickens, and cattle. In the case of cats, a lack of nuclear, genome-wide data, particularly from Near Eastern wildcats, has made it difficult to either detect or quantify this possibility. To address these issues, we generated 75 ancient mitochondrial genomes, 14 ancient nuclear genomes, and 31 modern nuclear genomes from European and Near Eastern wildcats. Our results demonstrate that despite cohabitating for at least 2,000 years on the European mainland and in Britain, most modern domestic cats possessed less than 10% of their ancestry from European wildcats, and ancient European wildcats possessed little to no ancestry from domestic cats. The antiquity and strength of this reproductive isolation between introduced domestic cats and local wildcats was likely the result of behavioral and ecological differences. Intriguingly, this long-lasting reproductive isolation is currently being eroded in parts of the species' distribution as a result of anthropogenic activities.

Modeling of African population history using f-statistics is biased when applying all previously proposed SNP ascertainment schemes

f-statistics have emerged as a first line of analysis for making inferences about demographic history from genome-wide data. Not only are they guaranteed to allow robust tests of the fits of proposed models of population history to data when analyzing full genome sequencing data-that is, all single nucleotide polymorphisms (SNPs) in the individuals being analyzed-but they are also guaranteed to allow robust tests of models for SNPs ascertained as polymorphic in a population that is an outgroup in a phylogenetic sense to all groups being analyzed. True "outgroup ascertainment" is in practice impossible in humans because our species has arisen from a substructured ancestral population that does not descend from a homogeneous ancestral population going back many hundreds of thousands of years into the past. However, initial studies suggested that non-outgroup-ascertainment schemes might produce robust enough results using f-statistics, and that motivated widespread fitting of models to data using non-outgroup-ascertained SNP panels such as the "Affymetrix Human Origins array" which has been genotyped on thousands of modern individuals from hundreds of populations, or the "1240k" in-solution enrichment reagent which has been the source of about 70% of published genome-wide data for ancient humans. In this study, we show that while analyses of population history using such panels work well for studies of relationships among non-African populations and one African outgroup, when co-modeling more than one sub-Saharan African and/or archaic human groups (Neanderthals and Denisovans), fitting of f-statistics to such SNP sets is expected to frequently lead to false rejection of true demographic histories, and failure to reject incorrect models. Analyzing panels of SNPs polymorphic in archaic humans, which has been suggested as a solution for the ascertainment problem, has limited statistical power and retains important biases. However, by carrying out simulations of diverse demographic histories, we show that bias in inferences based on f-statistics can be minimized by ascertaining on variants common in a union of diverse African groups; such ascertainment retains high statistical power while allowing co-analysis of archaic and modern groups.

The extinct Sicilian wolf shows a complex history of isolation and admixture with ancient dogs

The Sicilian wolf remained isolated in Sicily from the end of the Pleistocene until its extermination in the 1930s-1960s. Given its long-term isolation on the island and distinctive morphology, the genetic origin of the Sicilian wolf remains debated. We sequenced four nuclear genomes and five mitogenomes from the seven existing museum specimens to investigate the Sicilian wolf ancestry, relationships with extant and extinct wolves and dogs, and diversity. Our results show that the Sicilian wolf is most closely related to the Italian wolf but carries ancestry from a lineage related to European Eneolithic and Bronze Age dogs. The average nucleotide diversity of the Sicilian wolf was half of the Italian wolf, with 37-50% of its genome contained in runs of homozygosity. Overall, we show that, by the time it went extinct, the Sicilian wolf had high inbreeding and low-genetic diversity, consistent with a population in an insular environment.

Genome sequencing of 2000 canids by the Dog10K consortium advances the understanding of demography, genome function and architecture

BACKGROUND

The international Dog10K project aims to sequence and analyze several thousand canine genomes. Incorporating 20 × data from 1987 individuals, including 1611 dogs (321 breeds), 309 village dogs, 63 wolves, and four coyotes, we identify genomic variation across the canid family, setting the stage for detailed studies of domestication, behavior, morphology, disease susceptibility, and genome architecture and function.

RESULTS

We report the analysis of > 48 M single-nucleotide, indel, and structural variants spanning the autosomes, X chromosome, and mitochondria. We discover more than 75% of variation for 239 sampled breeds. Allele sharing analysis indicates that 94.9% of breeds form monophyletic clusters and 25 major clades. German Shepherd Dogs and related breeds show the highest allele sharing with independent breeds from multiple clades. On average, each breed dog differs from the UU_Cfam_GSD_1.0 reference at 26,960 deletions and 14,034 insertions greater than 50 bp, with wolves having 14% more variants. Discovered variants include retrogene insertions from 926 parent genes. To aid functional prioritization, single-nucleotide variants were annotated with SnpEff and Zoonomia phyloP constraint scores. Constrained positions were negatively correlated with allele frequency. Finally, the utility of the Dog10K data as an imputation reference panel is assessed, generating high-confidence calls across varied genotyping platform densities including for breeds not included in the Dog10K collection.

CONCLUSIONS

We have developed a dense dataset of 1987 sequenced canids that reveals patterns of allele sharing, identifies likely functional variants, informs breed structure, and enables accurate imputation. Dog10K data are publicly available.

Seriously cilia: A tiny organelle illuminates evolution, disease, and intercellular communication

The borders between cell and developmental biology, which have always been permeable, have largely dissolved. One manifestation is the blossoming of cilia biology, with cell and developmental approaches (increasingly complemented by human genetics, structural insights, and computational analysis) fruitfully advancing understanding of this fascinating, multifunctional organelle. The last eukaryotic common ancestor probably possessed a motile cilium, providing evolution with ample opportunity to adapt cilia to many jobs. Over the last decades, we have learned how non-motile, primary cilia play important roles in intercellular communication. Reflecting their diverse motility and signaling functions, compromised cilia cause a diverse range of diseases collectively called "ciliopathies." In this review, we highlight how cilia signal, focusing on how second messengers generated in cilia convey distinct information; how cilia are a potential source of signals to other cells; how evolution may have shaped ciliary function; and how cilia research may address thorny outstanding questions.

Multiple Origins and Genomic Basis of Complex Traits in Sighthounds

Sighthounds, a distinctive group of hounds comprising numerous breeds, have their origins rooted in ancient artificial selection of dogs. In this study, we performed genome sequencing for 123 sighthounds, including one breed from Africa, six breeds from Europe, two breeds from Russia, and four breeds and 12 village dogs from the Middle East. We gathered public genome data of five sighthounds and 98 other dogs as well as 31 gray wolves to pinpoint the origin and genes influencing the morphology of the sighthound genome. Population genomic analysis suggested that sighthounds originated from native dogs independently and were comprehensively admixed among breeds, supporting the multiple origins hypothesis of sighthounds. An additional 67 published ancient wolf genomes were added for gene flow detection. Results showed dramatic admixture of ancient wolves in African sighthounds, even more than with modern wolves. Whole-genome scan analysis identified 17 positively selected genes (PSGs) in the African population, 27 PSGs in the European population, and 54 PSGs in the Middle Eastern population. None of the PSGs overlapped in the three populations. Pooled PSGs of the three populations were significantly enriched in "regulation of release of sequestered calcium ion into cytosol" (gene ontology: 0051279), which is related to blood circulation and heart contraction. In addition, ESR1, JAK2, ADRB1, PRKCE, and CAMK2D were under positive selection in all three selected groups. This suggests that different PSGs in the same pathway contributed to the similar phenotype of sighthounds. We identified an ESR1 mutation (chr1: g.42,177,149 T > C) in the transcription factor (TF) binding site of Stat5a and a JAK2 mutation (chr1: g.93,277,007 T > A) in the TF binding site of Sox5. Functional experiments confirmed that the ESR1 and JAK2 mutation reduced their expression. Our results provide new insights into the domestication history and genomic basis of sighthounds.

Dog breeds and conformations in the UK in 2019: VetCompass canine demography and some consequent welfare implications

INTRODUCTION

Growing concerns over health and welfare impacts from extreme phenotypes in dogs have created an urgent need for reliable demographic information on the national breed structures of dogs.

METHODS

This study included all dogs under primary veterinary care in the UK during 2019 at practices participating in VetCompass. Demographic data on these dogs were analysed to report on the frequency of common breeds and also to report on conformation, bodyweight, sex and neuter associations with these breeds.

RESULTS

The study included 2,237,105 dogs under UK veterinary care in 2019. Overall, 69.4% (n = 1,551,462) were classified as purebred, 6.7% (149,308) as designer-crossbred and 24.0% (536,335) as nondesigner-crossbred. Across 800 unique breed names, the most frequent breeds at any age were nondesigner-crossbred (n = 536,335, 24.0%), Labrador Retriever (154,222, 6.9%) and Jack Russell Terrier (101,294, 4.5%). Among 229,624 (10.3%) dogs aged under one year, the most frequent breeds were nondesigner-crossbred (n = 45,995, 20.0%), French Bulldog (16,036, 7.0%) and Cockapoo (14,321, 6.2%). Overall, based on breed characteristics, 17.6% (395,739) were classified as brachycephalic, 43.1% (969,403) as mesaticephalic and 8.3% (186,320) as dolichocephalic. Of 1,551,336 dogs that were classifiable based on breed, 52.6% (815,673) were chondrodystrophic. Of 1,462,925 dogs that were classifiable, there were 54.6% (n = 798,426) short haired, 32.6% (476,883) medium haired and 12.8% (186,934) long haired. Of 1,547,653 dogs that were classifiable for ear carriage, 24.5% (n = 379,581) were erect, 28.1% (434,273) were semi-erect, 19.7% (305,475) were v-shaped drop and 27.7% (428,324) were pendulous. Overall, there was a 1.09:1.00 ratio of male (n = 1,163,512; 52.2%) to female dogs (n = 1,067,552; 47.8%).

CONCLUSIONS

Health and welfare issues linked to popular breeds with extreme phenotypes suggest that there is much work to do to help owners to make more welfare-friendly decisions when choosing which type of dog to own.

Origin and population structure of native dog breeds in the Korean peninsula and East Asia

To study the ancestry and phylogenetic relationships of native Korean dog breeds to other Asian dog populations, we analyzed nucleotide variations in whole-genome sequences of 205 canid individuals. Sapsaree, Northern Chinese indigenous dog, and Tibetan Mastiff were largely related to West Eurasian ancestry. Jindo, Donggyeongi, Shiba, Southern Chinese indigenous (SCHI), Vietnamese indigenous dogs (VIET), and Indonesian indigenous dogs were related to Southeast and East Asian ancestry. Among East Asian dog breeds, Sapsaree presented the highest haplotype sharing with German Shepherds, indicating ancient admixture of European ancestry to modern East Asian dog breeds. SCHI showed greater haplotype sharing with New Guinea singing dogs, VIET, and Jindo than with other Asian breeds. The predicted divergence time of East Asian populations from their common ancestor was approximately 2,000 to 11,000 years ago. Our results expand understanding of the genetic history of dogs in the Korean peninsula to the Asian continent and Oceanic region.

Possible origins and implications of atypical morphologies and domestication-like traits in wild golden jackals (Canis aureus)

Deciphering the origins of phenotypic variations in natural animal populations is a challenging topic for evolutionary and conservation biologists. Atypical morphologies in mammals are usually attributed to interspecific hybridisation or de-novo mutations. Here we report the case of four golden jackals (Canis aureus), that were observed during a camera-trapping wildlife survey in Northern Israel, displaying anomalous morphological traits, such as white patches, an upturned tail, and long thick fur which resemble features of domesticated mammals. Another individual was culled under permit and was genetically and morphologically examined. Paternal and nuclear genetic profiles, as well as geometric morphometric data, identified this individual as a golden jackal rather than a recent dog/wolf-jackal hybrid. Its maternal haplotype suggested past introgression of African wolf (Canis lupaster) mitochondrial DNA, as previously documented in other jackals from Israel. When viewed in the context of the jackal as an overabundant species in Israel, the rural nature of the surveyed area, the abundance of anthropogenic waste, and molecular and morphological findings, the possibility of an individual presenting incipient stages of domestication should also be considered.

Being a Dog: A Review of the Domestication Process

The process of canine domestication represents certainly one of the most interesting questions that evolutionary biology aims to address. A "multiphase" view of this process is now accepted, with a first phase during which different groups of wolves were attracted by the anthropogenic niche and a second phase characterized by the gradual establishment of mutual relationships between wolves and humans. Here, we provide a review of dog (Canis familiaris) domestication, highlighting the ecological differences between dogs and wolves, analyzing the molecular mechanisms which seem to have influenced the affiliative behaviors first observed in Belyaev's foxes, and describing the genetics of ancient European dogs. Then, we focus on three Mediterranean peninsulas (Balkan, Iberian and Italian), which together represent the main geographic area for studying canine domestication dynamics, as it has shaped the current genetic variability of dog populations, and where a well-defined European genetic structure was pinpointed through the analysis of uniparental genetic markers and their phylogeny.

Tracing Eastern Wolf Origins From Whole-Genome Data in Context of Extensive Hybridization

Southeastern Canada is inhabited by an amalgam of hybridizing wolf-like canids, raising fundamental questions regarding their taxonomy, origins, and timing of hybridization events. Eastern wolves (Canis lycaon), specifically, have been the subject of significant controversy, being viewed as either a distinct taxonomic entity of conservation concern or a recent hybrid of coyotes (C. latrans) and grey wolves (C. lupus). Mitochondrial DNA analyses show some evidence of eastern wolves being North American evolved canids. In contrast, nuclear genome studies indicate eastern wolves are best described as a hybrid entity, but with unclear timing of hybridization events. To test hypotheses related to these competing findings we sequenced whole genomes of 25 individuals, representative of extant Canadian wolf-like canid types of known origin and levels of contemporary hybridization. Here we present data describing eastern wolves as a distinct taxonomic entity that evolved separately from grey wolves for the past ∼67,000 years with an admixture event with coyotes ∼37,000 years ago. We show that Great Lakes wolves originated as a product of admixture between grey wolves and eastern wolves after the last glaciation (∼8,000 years ago) while eastern coyotes originated as a product of admixture between "western" coyotes and eastern wolves during the last century. Eastern wolf nuclear genomes appear shaped by historical and contemporary gene flow with grey wolves and coyotes, yet evolutionary uniqueness remains among eastern wolves currently inhabiting a restricted range in southeastern Canada.

Companion animal organoid technology to advance veterinary regenerative medicine

First year medical and veterinary students are made very aware that drugs can have very different effects in various species or even in breeds of one specific species. On the other hand, the “One Medicine” concept implies that therapeutic and technical approaches are exchangeable between man and animals. These opposing views on the (dis)similarities between human and veterinary medicine are magnified in regenerative medicine. Regenerative medicine promises to stimulate the body's own regenerative capacity via activation of stem cells and/or the application of instructive biomaterials. Although the potential is enormous, so are the hurdles that need to be overcome before large scale clinical implementation is realistic. It is in the advancement of regenerative medicine that veterinary regenerative medicine can play an instrumental and crucial role. This review describes the discovery of (adult) stem cells in domesticated animals, mainly cats and dogs. The promise of cell-mediated regenerative veterinary medicine is compared to the actual achievements, and this will lead to a set of unanswered questions (controversies, research gaps, potential developments in relation to fundamental, pre-clinical, and clinical research). For veterinary regenerative medicine to have impact, either for human medicine and/or for domesticated animals, answering these questions is pivotal.