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.

Imputation of ancient canid genomes reveals inbreeding history over the past 10,000 years

The multi-millenia long history between dogs and humans has placed them at the forefront of archeological and genomic research. Despite ongoing efforts including the analysis of ancient dog and wolf genomes, many questions remain regarding their geographic and temporal origins, and the microevolutionary processes that led to the diversity of breeds today. Although ancient genomes provide valuable information, their use is hindered by low depth of coverage and post-mortem damage, which inhibits confident genotype calling. In the present study, we assess how genotype imputation of ancient dog and wolf genomes, utilising a large reference panel, can improve the resolution provided by ancient datasets. Imputation accuracy was evaluated by down-sampling high coverage dog and wolf genomes to 0.05-2x coverage and comparing concordance between imputed and high coverage genotypes. We measured the impact of imputation on principal component analyses and runs of homozygosity. Our findings show high (R2>0.9) imputation accuracy for dogs with coverage as low as 0.5x and for wolves as low as 1.0x. We then imputed a dataset of 90 ancient dog and wolf genomes, to assess changes in inbreeding during the last 10,000 years of dog evolution. Ancient dog and wolf populations generally exhibited lower inbreeding levels than present-day individuals. Interestingly, regions with low ROH density maintained across ancient and present-day samples were significantly associated with genes related to olfaction and immune response. Our study indicates that imputing ancient canine genomes is a viable strategy that allows for the use of analytical methods previously limited to high-quality genetic data.

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.

IntroUNET: identifying introgressed alleles via semantic segmentation

A growing body of evidence suggests that gene flow between closely related species is a widespread phenomenon. Alleles that introgress from one species into a close relative are typically neutral or deleterious, but sometimes confer a significant fitness advantage. Given the potential relevance to speciation and adaptation, numerous methods have therefore been devised to identify regions of the genome that have experienced introgression. Recently, supervised machine learning approaches have been shown to be highly effective for detecting introgression. One especially promising approach is to treat population genetic inference as an image classification problem, and feed an image representation of a population genetic alignment as input to a deep neural network that distinguishes among evolutionary models (i.e. introgression or no introgression). However, if we wish to investigate the full extent and fitness effects of introgression, merely identifying genomic regions in a population genetic alignment that harbor introgressed loci is insufficient-ideally we would be able to infer precisely which individuals have introgressed material and at which positions in the genome. Here we adapt a deep learning algorithm for semantic segmentation, the task of correctly identifying the type of object to which each individual pixel in an image belongs, to the task of identifying introgressed alleles. Our trained neural network is thus able to infer, for each individual in a two-population alignment, which of those individual's alleles were introgressed from the other population. We use simulated data to show that this approach is highly accurate, and that it can be readily extended to identify alleles that are introgressed from an unsampled "ghost" population, performing comparably to a supervised learning method tailored specifically to that task. Finally, we apply this method to data from Drosophila, showing that it is able to accurately recover introgressed haplotypes from real data. This analysis reveals that introgressed alleles are typically confined to lower frequencies within genic regions, suggestive of purifying selection, but are found at much higher frequencies in a region previously shown to be affected by adaptive introgression. Our method's success in recovering introgressed haplotypes in challenging real-world scenarios underscores the utility of deep learning approaches for making richer evolutionary inferences from genomic data.

Genetic and cultural adaptations underlie the establishment of dairy pastoralism in the Tibetan Plateau

BACKGROUND

Domestication and introduction of dairy animals facilitated the permanent human occupation of the Tibetan Plateau. Yet the history of dairy pastoralism in the Tibetan Plateau remains poorly understood. Little is known how Tibetans adapted to milk and dairy products.

RESULTS

We integrated archeological evidence and genetic analysis to show the picture that the dairy ruminants, together with dogs, were introduced from West Eurasia into the Tibetan Plateau since ~ 3600 years ago. The genetic admixture between the exotic and indigenous dogs enriched the candidate lactase persistence (LP) allele 10974A > G of West Eurasian origin in Tibetan dogs. In vitro experiments demonstrate that - 13838G > A functions as a LP allele in Tibetans. Unlike multiple LP alleles presenting selective signatures in West Eurasians and South Asians, the de novo origin of Tibetan-specific LP allele - 13838G > A with low frequency (~ 6-7%) and absence of selection corresponds - 13910C > T in pastoralists across eastern Eurasia steppe.

CONCLUSIONS

Results depict a novel scenario of genetic and cultural adaptations to diet and expand current understanding of the establishment of dairy pastoralism in the Tibetan Plateau.

Genomic analysis of Tibetan ground tits identifies molecular adaptations associated with cooperative breeding

Cooperative breeding is a sophisticated altruistic social behavior that helps social animals to adapt to harsh environments. The Tibetan ground tit, Pseudopodoces humilis, is a high-altitude bird endemic to the Tibetan plateau. Recently, it has become an exciting system for studying the evolution of facultative cooperative breeding. To test for molecular adaptations associated with cooperative breeding, we resequenced the whole genome of ground tits from 6 wild populations that display remarkable variation in the frequency of cooperative breeding. Population structure analyses showed that the 6 populations were divided into 4 lineages, which is congruent with the major geographical distribution of the sampling sites. Using genome-wide selective sweep analysis, we identified putative positively selected genes (PSGs) in groups of tits that displayed high and low cooperative breeding rates. The total number of PSGs varied from 146 to 722 in high cooperative breeding rate populations, and from 272 to 752 in low cooperative breeding rate populations. Functional enrichment analysis of these PSGs identified several significantly enriched ontologies related to oxytocin signaling, estrogen signaling, and insulin secretion. PSGs involved in these functional ontologies suggest that molecular adaptations in hormonal regulation may have played important roles in shaping the evolution of cooperative breeding in the ground tit. Taken together, our study provides candidate genes and functional ontologies involved in molecular adaptations associated with cooperative breeding in Tibetan ground tits, and calls for a better understanding of the genetic roles in the evolution of cooperative breeding.

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.

The Australasian dingo archetype: de novo chromosome-length genome assembly, DNA methylome, and cranial morphology

BACKGROUND

One difficulty in testing the hypothesis that the Australasian dingo is a functional intermediate between wild wolves and domesticated breed dogs is that there is no reference specimen. Here we link a high-quality de novo long-read chromosomal assembly with epigenetic footprints and morphology to describe the Alpine dingo female named Cooinda. It was critical to establish an Alpine dingo reference because this ecotype occurs throughout coastal eastern Australia where the first drawings and descriptions were completed.

FINDINGS

We generated a high-quality chromosome-level reference genome assembly (Canfam_ADS) using a combination of Pacific Bioscience, Oxford Nanopore, 10X Genomics, Bionano, and Hi-C technologies. Compared to the previously published Desert dingo assembly, there are large structural rearrangements on chromosomes 11, 16, 25, and 26. Phylogenetic analyses of chromosomal data from Cooinda the Alpine dingo and 9 previously published de novo canine assemblies show dingoes are monophyletic and basal to domestic dogs. Network analyses show that the mitochondrial DNA genome clusters within the southeastern lineage, as expected for an Alpine dingo. Comparison of regulatory regions identified 2 differentially methylated regions within glucagon receptor GCGR and histone deacetylase HDAC4 genes that are unmethylated in the Alpine dingo genome but hypermethylated in the Desert dingo. Morphologic data, comprising geometric morphometric assessment of cranial morphology, place dingo Cooinda within population-level variation for Alpine dingoes. Magnetic resonance imaging of brain tissue shows she had a larger cranial capacity than a similar-sized domestic dog.

CONCLUSIONS

These combined data support the hypothesis that the dingo Cooinda fits the spectrum of genetic and morphologic characteristics typical of the Alpine ecotype. We propose that she be considered the archetype specimen for future research investigating the evolutionary history, morphology, physiology, and ecology of dingoes. The female has been taxidermically prepared and is now at the Australian Museum, Sydney.

The Australasian dingo archetype: De novo chromosome-length genome assembly, DNA methylome, and cranial morphology

Background

One difficulty in testing the hypothesis that the Australasian dingo is a functional intermediate between wild wolves and domesticated breed dogs is that there is no reference specimen. Here we link a high-quality de novo long read chromosomal assembly with epigenetic footprints and morphology to describe the Alpine dingo female named Cooinda. It was critical to establish an Alpine dingo reference because this ecotype occurs throughout coastal eastern Australia where the first drawings and descriptions were completed.

Findings

We generated a high-quality chromosome-level reference genome assembly (Canfam_ADS) using a combination of Pacific Bioscience, Oxford Nanopore, 10X Genomics, Bionano, and Hi-C technologies. Compared to the previously published Desert dingo assembly, there are large structural rearrangements on Chromosomes 11, 16, 25 and 26. Phylogenetic analyses of chromosomal data from Cooinda the Alpine dingo and nine previously published de novo canine assemblies show dingoes are monophyletic and basal to domestic dogs. Network analyses show that the mtDNA genome clusters within the southeastern lineage, as expected for an Alpine dingo. Comparison of regulatory regions identified two differentially methylated regions within glucagon receptor GCGR and histone deacetylase HDAC4 genes that are unmethylated in the Alpine dingo genome but hypermethylated in the Desert dingo. Morphological data, comprising geometric morphometric assessment of cranial morphology place dingo Cooinda within population-level variation for Alpine dingoes. Magnetic resonance imaging of brain tissue show she had a larger cranial capacity than a similar-sized domestic dog.

Conclusions

These combined data support the hypothesis that the dingo Cooinda fits the spectrum of genetic and morphological characteristics typical of the Alpine ecotype. We propose that she be considered the archetype specimen for future research investigating the evolutionary history, morphology, physiology, and ecology of dingoes. The female has been taxidermically prepared and is now at the Australian Museum, Sydney.

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) lived^1–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.

DNA from ancient wolves spanning 100,000 years sheds light on wolves’ evolutionary history and the genomic origin of dogs.

Genome-wide investigations reveal the population structure and selection signatures of Nigerian cattle adaptation in the sub-Saharan tropics

Background

Cattle are considered to be the most desirable livestock by small scale farmers. In Africa, although comprehensive genomic studies have been carried out on cattle, the genetic variations in indigenous cattle from Nigeria have not been fully explored. In this study, genome-wide analysis based on genotyping-by-sequencing (GBS) of 193 Nigerian cattle was used to reveal new insights on the history of West African cattle and their adaptation to the tropical African environment, particularly in sub-Saharan region. 

Results

The GBS data were evaluated against whole-genome sequencing (WGS) data and high rate of variant concordance between the two platforms was evident with high correlated genetic distance matrices genotyped by both methods suggestive of the reliability of GBS applicability in population genetics. The genetic structure of Nigerian cattle was observed to be homogenous and unique from other African cattle populations. Selection analysis for the genomic regions harboring imprints of adaptation revealed genes associated with immune responses, growth and reproduction, efficiency of feeds utilization, and heat tolerance. Our findings depict potential convergent adaptation between African cattle, dogs and humans with adaptive genes SPRY2 and ITGB1BP1 possibly involved in common physiological activities.

Conclusion

The study presents unique genetic patterns of Nigerian cattle which provide new insights on the history of cattle in West Africa based on their population structure and the possibility of parallel adaptation between African cattle, dogs and humans in Africa which require further investigations.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08512-w.

Natural and human-driven selection of a single non-coding body size variant in ancient and modern canids

Domestic dogs (Canis lupus familiaris) are the most variable-sized mammalian species on Earth, displaying a 40-fold size difference between breeds.¹ Although dogs of variable size are found in the archeological record,^2-4 the most dramatic shifts in body size are the result of selection over the last two centuries, as dog breeders selected and propagated phenotypic extremes within closed breeding populations.⁵ Analyses of over 200 domestic breeds have identified approximately 20 body size genes regulating insulin processing, fatty acid metabolism, TGFβ signaling, and skeletal formation.^6-10 Of these, insulin-like growth factor 1 (IGF1) predominates, controlling approximately 15% of body size variation between breeds.⁸ The identification of a functional mutation associated with IGF1 has thus far proven elusive.^6,10,11 Here, to identify and elucidate the role of an ancestral IGF1 allele in the propagation of modern canids, we analyzed 1,431 genome sequences from 13 species, including both ancient and modern canids, thus allowing us to define the evolutionary history of both ancestral and derived alleles at this locus. We identified a single variant in an antisense long non-coding RNA (IGF1-AS) that interacts with the IGF1 gene, creating a duplex. While the derived mutation predominates in both modern gray wolves and large domestic breeds, the ancestral allele, which predisposes to small size, was common in small-sized breeds and smaller wild canids. Our analyses demonstrate that this major regulator of canid body size nearly vanished in Pleistocene wolves, before its recent resurgence resulting from human-imposed selection for small-sized breed dogs.

Whole-Genome Sequencing Reveals Lactase Persistence Adaptation in European Dogs

Coexistence and cooperation between dogs and humans over thousands of years have supported convergent evolutionary processes in the two species. Previous studies found that Eurasian dogs evolved into a distinct geographic cluster. In this study, we used the genomes of 242 European dogs, 38 Southeast Asian indigenous (SEAI) dogs, and 41 gray wolves to identify adaptation of European dogs . We report 86 unique positively selected genes in European dogs, among which is LCT (lactase). LCT encodes lactase, which is fundamental for the digestion of lactose. We found that an A-to-G mutation (chr19:38,609,592) is almost fixed in Middle Eastern and European dogs. The results of two-dimensional site frequency spectrum (2D SFS) support that the mutation is under soft sweep . We inferred that the onset of positive selection of the mutation is shorter than 6,535 years and behind the well-developed dairy economy in central Europe. It increases the expression of LCT by reducing its binding with ZEB1, which would enhance dog's ability to digest milk-based diets. Our study uncovers the genetic basis of convergent evolution between humans and dogs with respect to diet, emphasizing the import of the dog as a biomedical model for studying mechanisms of the digestive system.

Evolutionary history of the extinct Sardinian dhole

The Sardinian dhole (Cynotherium sardous)¹ was an iconic and unique canid species that was endemic to Sardinia and Corsica until it became extinct at the end of the Late Pleistocene.^2-5 Given its peculiar dental morphology, small body size, and high level of endemism, several extant canids have been proposed as possible relatives of the Sardinian dhole, including the Asian dhole and African hunting dog ancestor.^3,6-9 Morphometric analyses^3,6,8-12 have failed to clarify the evolutionary relationship with other canids.We sequenced the genome of a ca-21,100-year-old Sardinian dhole in order to understand its genomic history and clarify its phylogenetic position. We found that it represents a separate taxon from all other living canids from Eurasia, Africa, and North America, and that the Sardinian dhole lineage diverged from the Asian dhole ca 885 ka. We additionally detected historical gene flow between the Sardinian and Asian dhole lineages, which ended approximately 500-300 ka, when the land bridge between Sardinia and mainland Italy was already broken, severing their population connectivity. Our sample showed low genome-wide diversity compared to other extant canids-probably a result of the long-term isolation-that could have contributed to the subsequent extinction of the Sardinian dhole.

Genome-scale sequencing and analysis of human, wolf, and bison DNA from 25,000-year-old sediment

Cave sediments have been shown to preserve ancient DNA but so far have not yielded the genome-scale information of skeletal remains. We retrieved and analyzed human and mammalian nuclear and mitochondrial environmental "shotgun" genomes from a single 25,000-year-old Upper Paleolithic sediment sample from Satsurblia cave, western Georgia:first, a human environmental genome with substantial basal Eurasian ancestry, which was an ancestral component of the majority of post-Ice Age people in the Near East, North Africa, and parts of Europe; second, a wolf environmental genome that is basal to extant Eurasian wolves and dogs and represents a previously unknown, likely extinct, Caucasian lineage; and third, a European bison environmental genome that is basal to present-day populations, suggesting that population structure has been substantially reshaped since the Last Glacial Maximum. Our results provide new insights into the Late Pleistocene genetic histories of these three species and demonstrate that direct shotgun sequencing of sediment DNA, without target enrichment methods, can yield genome-wide data informative of ancestry and phylogenetic relationships.

A shared pathogen: Babesia rossi in domestic dogs, black-backed jackals (Canis mesomelas) and African wild dogs (Lycaon pictus) in South Africa

In sub-Saharan Africa, babesiosis in domestic dogs is caused primarily by Babesia rossi. Black-backed jackals (Canis mesomelas), which are subclinical carriers of B. rossi, were a likely reservoir host from which infection passed to domestic dogs. The role of other indigenous canids, e.g. African wild dogs (Lycaon pictus), as reservoirs of B. rossi has not been elucidated. The question also arises whether genetic differences have arisen between B. rossi infecting domestic dogs and "ancestral" B. rossi in jackals. In a previous study we found that nearly one-third (27 of 91) of jackals were infected with B. rossi; this was confirmed by 18S rDNA sequence analysis. In this study, the near full-length B. rossi 18S rRNA gene was successfully amplified from 6 domestic dogs and 3 black-backed jackals. The obtained recombinant sequences were identical (100 %) to previously described B. rossi sequences of black-backed jackals in South Africa, and 99 % similar to B. rossi from dogs in South Africa and the Sudan. Although blood specimens from 5 (10 %) of 52 free-ranging African wild dogs (from Kruger National Park, South Africa, reacted with the B. rossi probe on RLB hybridisation, the presence of B. rossi could not be confirmed by amplification and sequencing, nor by multiplex, real-time PCR. Although African wild dogs they can be infected with B. rossi without showing clinical signs, our findings suggest that they are apparently not important reservoir hosts of B. rossi.

Whole genome analyses reveal significant convergence in obsessive-compulsive disorder between humans and dogs

Obsessive-compulsive disorder (OCD) represents a heterogeneous collection of diseases with diverse levels of phenotypic, genetic, and etiologic variability, making it difficult to identify the underlying genetic and biological mechanisms in humans. Domestic dogs exhibit several OCD-like behaviors. Using continuous circling as a representative phenotype for OCD, we screened two independent dog breeds, the Belgian Malinois and Kunming Dog and subsequently sequenced ten circling dogs and ten unaffected dogs for each breed. Using population differentiation analyses, we identified 11 candidate genes in the extreme tail of the differentiated regions between cases and controls. These genes overlap significantly with genes identified in a genome wide association study (GWAS) of human OCD, indicating strong convergence between humans and dogs. Through gene expressional analysis and functional exploration, we found that two candidate OCD risk genes, PPP2R2B and ADAMTSL3, affected the density and morphology of dendritic spines. Therefore, changes in dendritic spine may underlie some common biological and physiological pathways shared between humans and dogs. Our study revealed an unprecedented level of convergence in OCD shared between humans and dogs, and highlighted the importance of using domestic dogs as a model species for many human diseases including OCD.

Pleistocene climate fluctuations drove demographic history of African golden wolves (Canis lupaster)

Pleistocene climate change impacted entire ecosystems throughout the world. In the northern hemisphere, the distribution of Arctic species expanded during glacial periods, while more temperate and mesic species contracted into climatic refugia, where isolation drove genetic divergence. Cycles of local cooling and warming in the Sahara region of northern Africa caused repeated contractions and expansions of savannah-like environments which connected mesic species isolated in refugia during interglacial times, possibly driving population expansions and contractions; divergence and geneflow in the associated fauna. Here, we use whole genome sequences of African golden wolves (Canis lupaster), a generalist mesopredator with a wide distribution in northern Africa to estimate their demographic history and past episodes of geneflow. We detect a correlation between divergence times and cycles of increased aridity-associated Pleistocene glacial cycles. A complex demographic history with responses to local climate change in different lineages was found, including a relict lineage north of the High Atlas Mountains of Morocco that has been isolated for more than 18,000 years, possibly a distinct ecotype.

Origins and genetic legacy of prehistoric dogs

Dogs were the first domestic animal, but little is known about their population history and to what extent it was linked to humans. We sequenced 27 ancient dog genomes and found that all dogs share a common ancestry distinct from present-day wolves, with limited gene flow from wolves since domestication but substantial dog-to-wolf gene flow. By 11,000 years ago, at least five major ancestry lineages had diversified, demonstrating a deep genetic history of dogs during the Paleolithic. Coanalysis with human genomes reveals aspects of dog population history that mirror humans, including Levant-related ancestry in Africa and early agricultural Europe. Other aspects differ, including the impacts of steppe pastoralist expansions in West and East Eurasia and a near-complete turnover of Neolithic European dog ancestry.

Genetic signature of hybridization between Chinese spot-billed ducks and domesticated ducks

In this study, we analyzed 93 whole genomes from Chinese spot-billed ducks (CSB), meat-type ducks (MET), and egg and dual purpose-type ducks (EDT) to characterize the genetic material flowing between the CSB and modern ducks. Using a frequency of shared identical-by-descent method, approximately 10.9 Mb introgression segments containing 140 genes were identified showing the signatures of introgression between CSB and EDT. Meanwhile, nearly 10.6 M introgression regions containing 149 genes were identified between CSB and MET. Based on the haplotypes tree of each segment, we found that the introgression between CSB and domesticated ducks was asymmetric with a high level of gene flow from domestic to CSB and a low level of migration in the opposite direction. Moreover, we identified several genes that were introgressions from CSB and showed the signature of positive selection, which may contribute to the breeding of modern ducks. Our results provide new insight into the evolution and breeding history of domestic ducks and may be useful for the future management of wild and domestic duck populations.

Deciphering the puzzles of dog domestication

The domestic dog, as a highly successful domestication model, is well known as a favored human companion. Exploring its domestication history should provide great insight into our understanding of the prehistoric development of human culture and productivity. Furthermore, investigation on the mechanisms underpinning the morphological and behavioral traits associated with canid domestication syndrome is of significance not only for scientific study but also for human medical research. Current development of a multidisciplinary canine genome database, which includes enormous omics data, has substantially improved our understanding of the genetic makeup of dogs. Here, we reviewed recent advances associated with the original history and genetic basis underlying environmental adaptations and phenotypic diversities in domestic dogs, which should provide perspectives on improving the communicative relationship between dogs and humans.

Whole genome sequencing for the investigation of canine mammary tumor inheritance - an initial assessment of high-risk breast cancer genes reveal BRCA2 and STK11 variants potentially associated with risk in purebred dogs

Background

Although, in general, cancer is considered a multifactorial disease, clustering of particular cancers in pedigrees suggests a genetic predisposition and could explain why some dog breeds appear to have an increased risk of certain cancers. To our knowledge, there have been no published reports of whole genome sequencing to investigate inherited canine mammary tumor (CMT) risk, and with little known about CMT genetic susceptibility, we carried out whole genome sequencing on 14 purebred dogs diagnosed with mammary tumors from four breed-specific pedigrees. Following sequencing, each dog’s data was processed through a bioinformatics pipeline. This initial report highlights variants in orthologs of human breast cancer susceptibility genes.

Results

The overall whole genome and exome coverage averages were 26.0X and 25.6X, respectively, with 96.1% of the genome and 96.7% of the exome covered at least 10X. Of the average 7.9 million variants per dog, initial analyses involved surveying variants in orthologs of human breast cancer susceptibility genes, BRCA1, BRCA2, CDH1, PTEN, STK11, and TP53, and identified 19 unique coding variants that were validated through PCR and Sanger sequencing. Statistical analyses identified variants in BRCA2 and STK11 that appear to be associated with CMT, and breed-specific analyses revealed the breeds at the highest risk. Several additional BRCA2 variants showed trends toward significance, but have conflicting interpretations of pathogenicity, and correspond to variants of unknown significance in humans, which require further investigation. Variants in other genes were noted but did not appear to be associated with disease.

Conclusions

Whole genome sequencing proves to be an effective method to elucidate risk of CMT. Risk variants in orthologs of human breast cancer susceptibility genes have been identified. Ultimately, these whole genome sequencing efforts have provided a plethora of data that can also be assessed for novel discovery and have the potential to lead to breakthroughs in canine and human research through comparative analyses.

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.

Long-term follow-up of owned, free-roaming dogs in South Africa naturally exposed to Babesia rossi

Babesia rossi is an important, tick-borne intraerythrocytic protozoan parasite; however, its natural history and epidemiology is poorly understood. Babesia rossi is the most virulent Babesia sp. in domestic dogs and is generally considered to cause severe babesiosis, which is fatal if left untreated. However, subclinical infections and mild disease from B. rossi have been reported, although the clinical progression of these cases was not reported. Therefore, to better understand B. rossi under field conditions, we evaluated its clinical progression and seroprevalence in an owned, free-roaming dog population in Zenzele, South Africa, where the parasite is endemic and prevention is not routine. The entire dog population in Zenzele was monitored intensively at the individual level from March 2008 until April 2014, primarily for a longitudinal study on rabies control. Subsequent evaluation of B. rossi comprised analyses of clinical and laboratory data collected from the Zenzele dog population during the 6 year study period. A substantial proportion (31% (n = 34)) of 109 dogs (randomly selected from every available dog in February/March 2010 older than ~6-8 weeks (n = 246)) tested by Indirect Fluorescent Antibody Test had seroconverted strongly to B. rossi. All 34 dogs were generally consistently healthy adults, determined from regular clinical examinations between March 2008 and April 2014. Blood smear examinations at multiple time points between July 2009 and February 2011 were also undertaken for almost all of these (34) seropositive dogs and all those tested were consistently negative for Babesia spp. Subclinical infections and mild disease were also the main findings for a separate group of 18 dogs positive for Babesia spp. on blood smear examination and confirmed to be infected with B. rossi by Polymerase Chain Reaction - Reverse Line Blot. Almost all of these dogs were positive at only one time point from repeat blood smear examinations between July 2009 and February 2011. We suggest that these observations are consistent with immunity acquired from repeated, low-level exposure to the parasite, generating transient subclinical infections or mild disease. Should this be the case, the use of tick control, particularly in adult dogs in free-roaming populations in B. rossi endemic regions, should be carefully considered.

Demographic expansion of an African opportunistic carnivore during the Neolithic revolution

The diffusion of Neolithic technology together with the Holocene Climatic Optimum fostered the spread of human settlements and pastoral activities in North Africa, resulting in profound and enduring consequences for the dynamics of species, communities and landscapes. Here, we investigate the demographic history of the African wolf (Canis lupaster), a recently recognized canid species, to understand if demographic trends of this generalist and opportunistic carnivore reflect the increase in food availability that emerged after the arrival of the Neolithic economy in North Africa. We screened nuclear and mitochondrial DNA in samples collected throughout Algeria and Tunisia, and implemented coalescent approaches to estimate the variation of effective population sizes from present to ancestral time. We have found consistent evidence supporting the hypothesis that the African wolf population experienced a meaningful expansion concurring with a period of rapid population expansion of domesticates linked to the advent of agricultural practices.

Genomic Approaches Reveal an Endemic Subpopulation of Gray Wolves in Southern China

Although gray wolves (Canis lupus) are one of the most widely distributed terrestrial mammals, their origins in China are not well understood. We sequenced six specimens from wolf skins, showing that gray wolves from Southern China (SC) derive from a single lineage, distinct from gray wolves from the Tibetan Plateau and Northern China, suggesting that SC gray wolves may form a distinct subpopulation. Of SC gray wolves, one wolf from Zhejiang carries a genetic component from a canid and had gene flow from a population related to or further diverged from wolves than the dhole. This may indicate that interspecific gene flow likely played an important role in shaping the speciation patterns and population structure in the genus Canis. Our study is the first to survey museum gray wolves' genomes from Southern China, highlighting how sequencing the paleogenome from museum specimens can help us to study extinct species.

Structural variation during dog domestication: insights from gray wolf and dhole genomes

Several processes like phenotypic evolution, disease susceptibility and environmental adaptations, which fashion the domestication of animals, are largely attributable to structural variations (SVs) in the genome. Here, we present high-quality draft genomes of the gray wolf (Canis lupus) and dhole (Cuon alpinus) with scaffold N50 of 6.04 Mb and 3.96 Mb, respectively. Sequence alignment comprising genomes of three canid species reveals SVs specific to the dog, particularly 16 315 insertions, 2565 deletions, 443 repeats, 16 inversions and 15 translocations. Functional annotation of the dog SVs associated with genes indicates their enrichments in energy metabolisms, neurological processes and immune systems. Interestingly, we identify and verify at population level an insertion fully covering a copy of the AKR1B1 (Aldo-Keto Reductase Family 1 Member B) transcript. Transcriptome analysis reveals a high level of expression of the new AKR1B1 copy in the small intestine and liver, implying an increase in de novo fatty acid synthesis and antioxidant ability in dog compared to gray wolf, likely in response to dietary shifts during the agricultural revolution. For the first time, we report a comprehensive analysis of the evolutionary dynamics of SVs during the domestication step of dogs. Our findings demonstrate that retroposition can birth new genes to facilitate domestication, and affirm the importance of large-scale genomic variants in domestication studies.

Interspecific Gene Flow Shaped the Evolution of the Genus Canis

The evolutionary history of the wolf-like canids of the genus Canis has been heavily debated, especially regarding the number of distinct species and their relationships at the population and species level [1-6]. We assembled a dataset of 48 resequenced genomes spanning all members of the genus Canis except the black-backed and side-striped jackals, encompassing the global diversity of seven extant canid lineages. This includes eight new genomes, including the first resequenced Ethiopian wolf (Canis simensis), one dhole (Cuon alpinus), two East African hunting dogs (Lycaon pictus), two Eurasian golden jackals (Canis aureus), and two Middle Eastern gray wolves (Canis lupus). The relationships between the Ethiopian wolf, African golden wolf, and golden jackal were resolved. We highlight the role of interspecific hybridization in the evolution of this charismatic group. Specifically, we find gene flow between the ancestors of the dhole and African hunting dog and admixture between the gray wolf, coyote (Canis latrans), golden jackal, and African golden wolf. Additionally, we report gene flow from gray and Ethiopian wolves to the African golden wolf, suggesting that the African golden wolf originated through hybridization between these species. Finally, we hypothesize that coyotes and gray wolves carry genetic material derived from a "ghost" basal canid lineage.