Dog10K: an international sequencing effort to advance studies of canine domestication, phenotypes and health

Potential Adaptive Introgression From Dogs in Iberian Grey Wolves (Canis lupus)

Invading species along with increased anthropogenization may lead to hybridization events between wild species and closely related domesticates. As a consequence, wild species may carry introgressed alleles from domestic species, which is generally assumed to yield adverse effects in wild populations. The opposite evolutionary consequence, adaptive introgression, where introgressed genes are positively selected in the wild species, is possible but has rarely been documented. Grey wolves (Canis lupus) are widely distributed across the Holarctic and frequently coexist with their close relative, the domestic dog (C. familiaris). Despite ample opportunity, hybridization rarely occurs in most populations. Here we studied the geographically isolated grey wolves of the Iberian Peninsula, who have coexisted with a large population of loosely controlled dogs for thousands of years in a human-modified landscape. We assessed the extent and impact of dog introgression on the current Iberian grey wolf population by analysing 150 whole genomes of Iberian and other Eurasian grey wolves as well as dogs originating from across Europe and western Siberia. We identified almost no recent introgression and a small (< 5%) overall ancient dog ancestry. Using a combination of single scan statistics and ancestry enrichment estimates, we identified positive selection on six genes (DAPP1, NSMCE4A, MPPED2, PCDH9, MBTPS1, and CDH13) for which wild Iberian wolves carry alleles introgressed from dogs. The genes with introgressed and positively selected alleles include functions in immune response and brain functions, which may explain some of the unique behavioural phenotypes in Iberian wolves such as their reduced dispersal compared to other wolf populations.

Going back to "basics": Harlow's learning set task with wolves and dogs

To survive and reproduce, animals need to behave adaptively by adjusting their behavior to their environment, with learning facilitating some of these processes. Dogs have become a go-to model species in comparative cognition studies, making our understanding of their learning skills paramount at multiple levels, not only with regards to basic research on their cognitive skills and the effects of domestication, but also with applied purposes such as training. In order to tackle these issues, we tested similarly raised wolves and dogs in a serial learning task inspired by Harlow's "learning set." In Phase 1, different pairs of objects were presented to the animals, one of which was baited while the other was not. Both species' performance gradually improved with each new set of objects, showing that they "learnt to learn," but no differences were found between the species in their learning speed. In Phase 2, once subjects had learned the association between one of the objects and the food reward, the contingencies were reversed and the previously unrewarded object of the same pair was now rewarded. Dogs' performance in this task seemed to be better than wolves', albeit only when considering just the first session of each reversal, suggesting that the dogs might be more flexible than wolves. Further research (possibly with the aid of refined methods such as computer-based tasks) would help ascertain whether these differences between wolves and dogs are persistent across different learning tasks.

Socio-economic evaluation of human-dog coexistence: A 40,000 years history

The human-animal relationship and its benefits have been a popular issue in recent decades. This study focuses on the dog as a paradigmatic example of this complex bond. The relationship between human and dog evolved over thousands of years through the domestication process. During this period, climatic, cultural, geographical and social changes have had a direct impact on this bond. Dogs have a remarkable capacity to play a variety of roles that mirror the evolution of human societies. Subsequently, questions arise regarding the important and symbolic roles that canines played in ancient cultures, religions and mythological traditions. Specifically, how have these roles changed through history? This research aimed to explore the socio-economic aspects of the human-dog relationship, spanning from ancient civilizations to modern societies. Through a systematic search methodology, this investigation seeks to provide insights into the dynamics of the complex human-dog relationship, including historical backgrounds, socio-cultural dynamics, economic implications and associated challenges. Additionally, the present study addresses the environmental concerns that emerge alongside rising stray dog populations. This research emphasizes the importance of strategies to balance the benefits of this coexistence that promote human and animal welfare.

Dog10K: an integrated Dog10K database summarizing canine multi-omics

The diversity observed in canine breed phenotypes, together with their risk for heritabily disorders of relevance to dogs and humans, makes the species an ideal subject for studies aimed at understanding the genetic basis of complex traits and human biomedical models. Dog10K is an ongoing international collaboration that aims to uncover the genetic basis of phenotypic diversity, disease, behavior, and domestication history of dogs. To best present and make the extensive data accessible and user friendly, we have established the Dog10K (http://dog10k.kiz.ac.cn/) database, a comprehensive-omics resource summarizing multiple types of data. This database integrates single nucleotide variants (SNVs) from 1987 canine genomes, de-novo mutations (DNMs) from 43 dog breeds with >40× sequence, RNA-seq data of 105057 single nuclei from hippocampus, 74067 single cells from leukocytes and 30 blood samples from published canid studies. We provide clear visualization, statistics, browse, searching, and downloading functions for all data. We have integrated three analysis tools, Selscan, LiftOver and AgeConversion, to aid researchers in custom exploration of the comprehensive-omics data. The Dog10K database will serve as a foundational platform for analyzing, presenting and utilizing canine multi-omics data.

Neuronal ceroid lipofuscinosis in a Schapendoes dog is caused by a missense variant in CLN6

Neuronal ceroid lipofuscinosis (NCL) is a group of neurodegenerative disorders that occur in humans, dogs, and several other species. NCL is characterised clinically by progressive deterioration of cognitive and motor function, epileptic seizures, and visual impairment. Most forms present early in life and eventually lead to premature death. Typical pathological changes include neuronal accumulation of autofluorescent, periodic acid-Schiff- and Sudan black B-positive lipopigments, as well as marked loss of neurons in the central nervous system. Here, we describe a 19-month-old Schapendoes dog, where clinical signs were indicative of lysosomal storage disease, which was corroborated by pathological findings consistent with NCL. Whole genome sequencing of the affected dog and both parents, followed by variant calling and visual inspection of known NCL genes, identified a missense variant in CLN6 (c.386T>C). The variant is located in a highly conserved region of the gene and predicted to be harmful, which supports a causal relationship. The identification of this novel CLN6 variant enables pre-breeding DNA-testing to prevent future cases of NCL6 in the Schapendoes breed, and presents a potential natural model for NCL6 in humans.

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.

The value of hybrid genomes: Building two highly contiguous reference genome assemblies to advance Canis genomic studies

Previous studies of canid population and evolutionary genetics have relied on high-quality domestic dog reference genomes that have been produced primarily for biomedical and trait mapping studies in dog breeds. However, the absence of highly contiguous genomes from other Canis species like the gray wolf and coyote, that represent additional distinct demographic histories, may bias inferences regarding interspecific genetic diversity and phylogenetic relationships. Here, we present single haplotype de novo genome assemblies for the gray wolf and coyote, generated by applying the trio-binning approach to long sequence reads generated from the genome of a female first-generation hybrid produced from a gray wolf and coyote mating. The assemblies were highly contiguous, with contig N50 sizes of 44.6 and 42.0 Mb for the wolf and coyote, respectively. Genome scaffolding and alignments between the two Canis assemblies and published dog reference genomes showed near complete collinearity, with one exception: a coyote-specific chromosome fission of chromosome 13 and fusion of the proximal portion of that chromosome with chromosome 8, retaining the Canis-typical haploid chromosome number of 2n = 78. We evaluated mapping quality for previous RADseq data from 334 canids and found nearly identical mapping quality and patterns among canid species and regional populations regardless of the genome used for alignment (dog, coyote, or gray wolf). These novel wolf and coyote genome reference assemblies will be important resources for proper and accurate inference of Canis demography, taxonomic evaluation, and conservation genetics.

Large variance in inbreeding within the Iberian wolf population

The gray wolf (Canis lupus) population on the Iberian Peninsula was the largest in western and central Europe during most of the 20th century, with its size apparently never under a few hundred individuals. After partial legal protection in the 1970s in Spain, the northwest Iberian population increased to about 300 to 350 packs and then stabilized. In contrast to many current European wolf populations, which have been connected through gene flow, the Iberian wolf population has been isolated for decades. Here, we measured changes in genomic diversity and inbreeding through the last decades in a geographic context. We find that the level of genomic diversity in Iberian wolves is low compared with other Eurasian wolf populations. Despite population expansion in the last 50 years, some modern wolves had very high inbreeding, especially in the recently recolonized and historical edge areas. These individuals contrast with others with low inbreeding within the same population. The high variance in inbreeding despite population expansion seems associated with small-scale fragmentation of the range that is revealed by the genetic similarity between modern and historical samples from close localities despite being separated by decades, remaining differentiated from other individuals that are just over 100 km away, a small distance for a species with great dispersal capacity inhabiting a continuous range. This illustrates that, despite its demographically stable condition, the population would probably benefit from favoring connectivity within the population as well as genetic exchange with other European wolf populations to avoid excessive fragmentation and local inbreeding depression.

The Addis Ababa Lions: Whole-Genome Sequencing of a Rare and Precious Population

Lions are widely known as charismatic predators that once roamed across the globe, but their populations have been greatly affected by environmental factors and human activities over the last 150 yr. Of particular interest is the Addis Ababa lion population, which has been maintained in captivity at around 20 individuals for over 75 yr, while many wild African lion populations have become extinct. In order to understand the molecular features of this unique population, we conducted a whole-genome sequencing study on 15 Addis Ababa lions and detected 4.5 million distinct genomic variants compared with the reference African lion genome. Using functional annotation, we identified several genes with mutations that potentially impact various traits such as mane color, body size, reproduction, gastrointestinal functions, cardiovascular processes, and sensory perception. These findings offer valuable insights into the genetics of this threatened lion population.

A cautionary tale of low-pass sequencing and imputation with respect to haplotype accuracy

BACKGROUND

Low-pass whole-genome sequencing and imputation offer significant cost savings, enabling substantial increases in sample size and statistical power. This approach is particularly promising in livestock breeding, providing an affordable means of screening individuals for deleterious alleles or calculating genomic breeding values. Consequently, it may also be of value in companion animal genomics to support pedigree breeding. We sought to evaluate in dogs the impact of low coverage sequencing and reference-guided imputation on genotype concordance and association analyses.

RESULTS

DNA isolated from saliva of 30 Labrador retrievers was sequenced at low (0.9X and 3.8X) and high (43.5X) coverage, and down-sampled from 43.5X to 9.6X and 17.4X. Genotype imputation was performed using a diverse reference panel (1021 dogs), and two subsets of the former panel (256 dogs each) where one had an excess of Labrador retrievers relative to other breeds. We observed little difference in imputed genotype concordance between reference panels. Association analyses for a locus acting as a disease proxy were performed using single-marker (GEMMA) and haplotype-based (XP-EHH) tests. GEMMA results were highly correlated (r ≥ 0.97) between 43.5X and ≥ 3.8X depths of coverage, while for 0.9X the correlation was lower (r ≤ 0.8). XP-EHH results were less well correlated, with r ranging from 0.58 (0.9X) to 0.88 (17.4X). Across a random sample of 10,000 genomic regions averaging 17 kb in size, we observed a median of three haplotypes per dog across the sequencing depths, with 5% of the regions returning more than eight haplotypes. Inspection of one such region revealed genotype and phasing inconsistencies across sequencing depths.

CONCLUSIONS

We demonstrate that saliva-derived canine DNA is suitable for whole-genome sequencing, highlighting the feasibility of client-based sampling. Low-pass sequencing and imputation require caution as incorrect allele assignments result when the subject possesses alleles that are absent in the reference panel. Larger panels have the capacity for greater allelic diversity, which should reduce the potential for imputation error. Although low-pass sequencing can accurately impute allele dosage, we highlight issues with phasing accuracy that impact haplotype-based analyses. Consequently, if accurately phased genotypes are required for analyses, we advocate sequencing at high depth (> 20X).

The P10K database: a data portal for the protist 10 000 genomes project

Protists, a highly diverse group of microscopic eukaryotic organisms distinct from fungi, animals and plants, exert crucial roles within the earth's biosphere. However, the genomes of only a small fraction of known protist species have been published and made publicly accessible. To address this constraint, the Protist 10 000 Genomes Project (P10K) was initiated, implementing a specialized pipeline for single-cell genome/transcriptome assembly, decontamination and annotation of protists. The resultant P10K database (https://ngdc.cncb.ac.cn/p10k/) serves as a comprehensive platform, collating and disseminating genome sequences and annotations from diverse protist groups. Currently, the P10K database has incorporated 2959 genomes and transcriptomes, including 1101 newly sequenced datasets by P10K and 1858 publicly available datasets. Notably, it covers 45% of the protist orders, with a significant representation (53% coverage) of ciliates, featuring nearly a thousand genomes/transcriptomes. Intriguingly, analysis of the unique codon table usage among ciliates has revealed differences compared to the NCBI taxonomy system, suggesting a need to revise the codon tables used for these species. Collectively, the P10K database serves as a valuable repository of genetic resources for protist research and aims to expand its collection by incorporating more sequenced data and advanced analysis tools to benefit protist studies worldwide.

Recent Advances in the Nutrition and Metabolism of Dogs and Cats

Domestic dogs (facultative carnivores) and cats (obligate carnivores) have been human companions for at least 12,000 and 9000 years, respectively. These animal species have a relatively short digestive tract but a large stomach volume and share many common features of physiological processes, intestinal microbes, and nutrient metabolism. The taste buds of the canine and feline tongues can distinguish sour, umami, bitter, and salty substances. Dogs, but not cats, possess sweet receptors. α-Amylase activity is either absent or very low in canine and feline saliva, and is present at low or substantial levels in the pancreatic secretions of cats or dogs, respectively. Thus, unlike cats, dogs have adapted to high-starch rations while also consuming animal-sourced foods. At metabolic levels, both dogs and cats synthesize de novo vitamin C and many amino acids (AAs, such as Ala, Asn, Asp, Glu, Gln, Gly, Pro, and Ser) but have a very limited ability to form vitamin D3. Compared with dogs, cats have higher requirements for AAs, some B-complex vitamins, and choline; greater rates of gluconeogenesis; a higher capacity to tolerate AA imbalances and antagonism; a more limited ability to synthesize arginine and taurine from glutamine/proline and cysteine, respectively; and a very limited ability to generate polyunsaturated fatty acids (PUFAs) from respective substrates. Unlike dogs, cats cannot convert either β-carotene into vitamin A or tryptophan into niacin. Dogs can thrive on one large meal daily and select high-fat over low-fat diets, whereas cats eat more frequently during light and dark periods and select high-protein over low-protein diets. There are increasing concerns over the health of skin, hair, bone, and joints (specialized connective tissues containing large amounts of collagen and/or keratin); sarcopenia (age-related losses of skeletal-muscle mass and function); and cognitive function in dogs and cats. Sufficient intakes of proteinogenic AAs and taurine along with vitamins, minerals, and PUFAs are crucial for the normal structures of the skin, hair, bone, and joints, while mitigating sarcopenia and cognitive dysfunction. Although pet owners may have different perceptions about the feeding and management practice of their dogs and cats, the health and well-being of the companion animals critically depend on safe, balanced, and nutritive foods. The new knowledge covered in this volume of Adv Exp Med Biol is essential to guide the formulation of pet foods to improve the growth, development, brain function, reproduction, lactation, and health of the companion animals.

Leading the pack: Best practices in comparative canine cancer genomics to inform human oncology

Pet dogs develop spontaneous cancers at a rate estimated to be five times higher than that of humans, providing a unique opportunity to study disease biology and evaluate novel therapeutic strategies in a model system that possesses an intact immune system and mirrors key aspects of human cancer biology. Despite decades of interest, effective utilization of pet dog cancers has been hindered by a limited repertoire of necessary cellular and molecular reagents for both in vitro and in vivo studies, as well as a dearth of information regarding the genomic landscape of these cancers. Recently, many of these critical gaps have been addressed through the generation of a highly annotated canine reference genome, the creation of several tools necessary for multi-omic analysis of canine tumours, and the development of a centralized repository for key genomic and associated clinical information from canine cancer patients, the Integrated Canine Data Commons. Together, these advances have catalysed multidisciplinary efforts designed to integrate the study of pet dog cancers more effectively into the translational continuum, with the ultimate goal of improving human outcomes. The current review summarizes this recent progress and provides a guide to resources and tools available for comparative study of pet dog cancers.

Domestication and breeding objective did not shape the interpretation of physical and social cues in goats (Capra hircus)

Artificial selection by humans, either through domestication or subsequent selection for specific breeding objectives, drives changes in animal cognition and behaviour. However, most previous cognitive research comparing domestic and wild animals has focused on companion animals such as canids, limiting any general claims about the effects of artificial selection by humans. Using a cognitive test battery, we investigated the ability of wild goats (non-domestic, seven subjects), dwarf goats (domestic, not selected for milk production, 15 subjects) and dairy goats (domestic, selected for high milk yield, 18 subjects) to utilise physical and social cues in an object choice task. To increase the heterogeneity of our test samples, data for domestic goats were collected by two experimenters at two research stations (Agroscope; Research Institute for Farm Animal Biology). We did not find performance differences between the three groups in the cognitive test battery for either physical or social cues. This indicates that for a domestic non-companion animal species, domestication and selection for certain breeding objectives did not measurably shape the physical and cognitive skills of goats.

A frameshift-deletion mutation in Reelin causes cerebellar hypoplasia in White Swiss Shepherd dogs

Cerebellar hypoplasia is a heterogeneous neurological condition in which the cerebellum is smaller than usual or not completely developed. The condition can have genetic origins, with Mendelian-effect mutations described in several mammalian species. Here, we describe a genetic investigation of cerebellar hypoplasia in White Swiss Shepherd dogs, where two affected puppies were identified from a litter with a recent common ancestor on both sides of their pedigree. Whole genome sequencing was conducted for 10 dogs in this family, and filtering of these data based on a recessive transmission hypothesis highlighted five protein-altering candidate variants - including a frameshift-deletion of the Reelin (RELN) gene (p.Val947*). Given the status of RELN as a gene responsible for cerebellar hypoplasia in humans, sheep and mice, these data strongly suggest the loss-of-function variant as underlying these effects. This variant has not been found in other dog breeds nor in a cohort of European White Swiss Shepherds, suggesting a recent mutation event. This finding will support the genotyping of a more diverse sample of dogs, and should aid future management of the harmful allele through optimised mating schemes.

GWAS using low-pass whole genome sequence reveals a novel locus in canine congenital idiopathic megaesophagus

Congenital idiopathic megaesophagus (CIM) is a gastrointestinal disorder of dogs wherein the esophagus is dilated and swallowing activity is reduced, causing regurgitation of ingesta. Affected individuals experience weight loss and malnourishment and are at risk for aspiration pneumonia, intussusception, and euthanasia. Great Danes have among the highest incidences of CIM across dog breeds, suggesting a genetic predisposition. We generated low-pass sequencing data for 83 Great Danes and used variant calls to impute missing whole genome single-nucleotide variants (SNVs) for each individual based on haplotypes phased from 624 high-coverage dog genomes, including 21 Great Danes. We validated the utility of our imputed data set for genome-wide association studies (GWASs) by mapping loci known to underlie coat phenotypes with simple and complex inheritance patterns. We conducted a GWAS for CIM with 2,010,300 SNVs, identifying a novel locus on canine chromosome 1 (P-val = 2.76 × 10-10). Associated SNVs are intergenic or intronic and are found in two clusters across a 1.7-Mb region. Inspection of coding regions in high-coverage genomes from affected Great Danes did not reveal candidate causal variants, suggesting that regulatory variants underlie CIM. Further studies are necessary to assess the role of these non-coding variants.

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.

Whole Animal Genome Sequencing: user-friendly, rapid, containerized pipelines for processing, variant discovery, and annotation of short-read whole genome sequencing data

Advancements in massively parallel short-read sequencing technologies and the associated decreasing costs have led to large and diverse variant discovery efforts across species. However, processing high-throughput short-read sequencing data can be challenging with potential pitfalls and bioinformatics bottlenecks in generating reproducible results. Although a number of pipelines exist that address these challenges, these are often geared toward human or traditional model organism species and can be difficult to configure across institutions. Whole Animal Genome Sequencing (WAGS) is an open-source set of user-friendly, containerized pipelines designed to simplify the process of identifying germline short (SNP and indel) and structural variants (SVs) geared toward the veterinary community but adaptable to any species with a suitable reference genome. We present a description of the pipelines [adapted from the best practices of the Genome Analysis Toolkit (GATK)], along with benchmarking data from both the preprocessing and joint genotyping steps, consistent with a typical user workflow.

Comparative neurogenetics of dog behavior complements efforts towards human neuropsychiatric genetics

Domestic dogs display a wide array of heritable behaviors that have intermediate genetic complexity thanks to a long history of human-influenced selection. Comparative genetics in dogs could address the scarcity of non-human neurogenetic systems relevant to human neuropsychiatric disorders, which are characterized by mental, emotional, and behavioral symptoms and involve vastly complex genetic and non-genetic risk factors. Our review describes the diverse behavioral "phenome" of domestic dogs, past and ongoing sources of behavioral selection, and the state of canine behavioral genetics. We highlight two naturally disordered behavioral domains that illustrate how dogs may prove useful as a comparative, forward neurogenetic system: canine age-related cognitive dysfunction, which can be examined more rapidly given the attenuated lifespan of dogs, and compulsive disorders, which may have genetic roots in purpose-bred behaviors. Growing community science initiatives aimed at the companion dog population will be well suited to investigating such complex behavioral phenotypes and offer a comparative resource that parallels human genomic initiatives in scale and dimensionality.

A loss-of-function variant in canine GLRA1 associates with a neurological disorder resembling human hyperekplexia

Hereditary hyperekplexia is a rare neuronal disorder characterized by an exaggerated startle response to sudden tactile or acoustic stimuli. In this study, we present a Miniature Australian Shepherd family showing clinical signs, which have genetic and phenotypic similarities with human hereditary hyperekplexia: episodes of muscle stiffness that could occasionally be triggered by acoustic stimuli. Whole genome sequence data analysis of two affected dogs revealed a 36-bp deletion spanning the exon-intron boundary in the glycine receptor alpha 1 (GLRA1) gene. Further validation in pedigree samples and an additional cohort of 127 Miniature Australian Shepherds, 45 Miniature American Shepherds and 74 Australian Shepherds demonstrated complete segregation of the variant with the disease, according to an autosomal recessive inheritance pattern. The protein encoded by GLRA1 is a subunit of the glycine receptor, which mediates postsynaptic inhibition in the brain stem and spinal cord. The canine GLRA1 deletion is located in the signal peptide and is predicted to cause exon skipping and subsequent premature stop codon resulting in a significant defect in glycine signaling. Variants in GLRA1 are known to cause hereditary hyperekplexia in humans; however, this is the first study to associate a variant in canine GLRA1 with the disorder, establishing a spontaneous large animal disease model for the human condition.

Common Orthopedic Traits and Screening for Breeding Programs

Orthopedic diseases are complex traits, meaning genetics and environmental factors affect risk, making identification of genetic associations difficult. In the United States, hip and elbow scores, patellar luxation scores, Legg-Calvé-Perthes disease, and shoulder osteochondrosis affectedness are available in the Orthopedic Foundation for Animals registry. Distraction indices and extended, ventrodorsal hip conformation scores are recorded by PennHIP. Application of estimated breeding values for hip and elbow dysplasia in breeder selection reduces the severity and prevalence of these traits. Genomic prediction and whole-genome sequence technologies and methods should improve knowledge of genetics underlying orthopedic diseases, leading to improved canine orthopedic genetic quality.

2022 updates to the Rat Genome Database: a Findable, Accessible, Interoperable, and Reusable (FAIR) resource

The Rat Genome Database (RGD, https://rgd.mcw.edu) has evolved from simply a resource for rat genetic markers, maps, and genes, by adding multiple genomic data types and extensive disease and phenotype annotations and developing tools to effectively mine, analyze, and visualize the available data, to empower investigators in their hypothesis-driven research. Leveraging its robust and flexible infrastructure, RGD has added data for human and eight other model organisms (mouse, 13-lined ground squirrel, chinchilla, naked mole-rat, dog, pig, African green monkey/vervet, and bonobo) besides rat to enhance its translational aspect. This article presents an overview of the database with the most recent additions to RGD's genome, variant, and quantitative phenotype data. We also briefly introduce Virtual Comparative Map (VCMap), an updated tool that explores synteny between species as an improvement to RGD's suite of tools, followed by a discussion regarding the refinements to the existing PhenoMiner tool that assists researchers in finding and comparing quantitative data across rat strains. Collectively, RGD focuses on providing a continuously improving, consistent, and high-quality data resource for researchers while advancing data reproducibility and fulfilling Findable, Accessible, Interoperable, and Reusable (FAIR) data principles.

Genomic Diversity and Runs of Homozygosity in Bernese Mountain Dogs

Bernese mountain dogs are a large dog breed formed in the early 1900s in Switzerland. While originally farm dogs that were used for pulling carts, guarding, and driving cattle, today they are considered multi-purpose companion and family dogs. The breed is predisposed to several complex diseases, such as histiocytic sarcoma, degenerative myelopathy, or hip dysplasia. Using whole-genome sequencing (WGS) data, we assessed the genomic architecture of 33 unrelated dogs from four countries: France, Sweden, Switzerland, and the United States. Analysis of runs of homozygosity (ROH) identified 12,643 ROH with an average length of 2.29 Mb and an average inbreeding coefficient of 0.395. Multidimensional scaling analysis of the genetic relatedness revealed limited clustering of European versus USA dogs, suggesting exchanges of breeding stock between continents. Furthermore, only two mtDNA haplotypes were detected in the 33 studied dogs, both of which are widespread throughout multiple dog breeds. WGS-based ROH analyses revealed several fixed or nearly fixed regions harboring discreet morphological trait-associated as well as disease-associated genetic variants. Several genes involved in the regulation of immune cells were found in the ROH shared by all dogs, which is notable in the context of the breed’s strong predisposition to hematopoietic cancers. High levels of inbreeding and relatedness, strongly exaggerated in the last 30 years, have likely led to the high prevalence of specific genetic disorders in this breed.

Amino acid nutrition and metabolism in domestic cats and dogs

Domestic cats and dogs are carnivores that have evolved differentially in the nutrition and metabolism of amino acids. This article highlights both proteinogenic and nonproteinogenic amino acids. Dogs inadequately synthesize citrulline (the precursor of arginine) from glutamine, glutamate, and proline in the small intestine. Although most breeds of dogs have potential for adequately converting cysteine into taurine in the liver, a small proportion (1.3%-2.5%) of the Newfoundland dogs fed commercially available balanced diets exhibit a deficiency of taurine possibly due to gene mutations. Certain breeds of dogs (e.g., golden retrievers) are more prone to taurine deficiency possibly due to lower hepatic activities of cysteine dioxygenase and cysteine sulfinate decarboxylase. De novo synthesis of arginine and taurine is very limited in cats. Thus, concentrations of both taurine and arginine in feline milk are the greatest among domestic mammals. Compared with dogs, cats have greater endogenous nitrogen losses and higher dietary requirements for many amino acids (e.g., arginine, taurine, cysteine, and tyrosine), and are less sensitive to amino acid imbalances and antagonisms. Throughout adulthood, cats and dogs may lose 34% and 21% of their lean body mass, respectively. Adequate intakes of high-quality protein (i.e., 32% and 40% animal protein in diets of aging dogs and cats, respectively; dry matter basis) are recommended to alleviate aging-associated reductions in the mass and function of skeletal muscles and bones. Pet-food grade animal-sourced foodstuffs are excellent sources of both proteinogenic amino acids and taurine for cats and dogs, and can help to optimize their growth, development, and health.

Bayesian mixed model analysis uncovered 21 risk loci for chronic kidney disease in boxer dogs

Chronic kidney disease (CKD) affects 10% of the human population, with only a small fraction genetically defined. CKD is also common in dogs and has been diagnosed in nearly all breeds, but its genetic basis remains unclear. Here, we performed a Bayesian mixed model genome-wide association analysis for canine CKD in a boxer population of 117 canine cases and 137 controls, and identified 21 genetic regions associated with the disease. At the top markers from each CKD region, the cases carried an average of 20.2 risk alleles, significantly higher than controls (15.6 risk alleles). An ANOVA test showed that the 21 CKD regions together explained 57% of CKD phenotypic variation in the population. Based on whole genome sequencing data of 20 boxers, we identified 5,206 variants in LD with the top 50 BayesR markers. Following comparative analysis with human regulatory data, 17 putative regulatory variants were identified and tested with electrophoretic mobility shift assays. In total four variants, three intronic variants from the MAGI2 and GALNT18 genes, and one variant in an intergenic region on chr28, showed alternative binding ability for the risk and protective alleles in kidney cell lines. Many genes from the 21 CKD regions, RELN, MAGI2, FGFR2 and others, have been implicated in human kidney development or disease. The results from this study provide new information that may enlighten the etiology of CKD in both dogs and humans.

Domestic dog lineages reveal genetic drivers of behavioral diversification

Selective breeding of domestic dogs has generated diverse breeds often optimized for performing specialized tasks. Despite the heritability of breed-typical behavioral traits, identification of causal loci has proven challenging due to the complexity of canine population structure. We overcome longstanding difficulties in identifying genetic drivers of canine behavior by developing a framework for understanding relationships between breeds and the behaviors that define them, utilizing genetic data for over 4,000 domestic, semi-feral, and wild canids and behavioral survey data for over 46,000 dogs. We identify ten major canine genetic lineages and their behavioral correlates and show that breed diversification is predominantly driven by non-coding regulatory variation. We determine that lineage-associated genes converge in neurodevelopmental co-expression networks, identifying a sheepdog-associated enrichment for interrelated axon guidance functions. This work presents a scaffold for canine diversification that positions the domestic dog as an unparalleled system for revealing the genetic origins of behavioral diversity.

Across-breed genetic investigation of canine hip dysplasia, elbow dysplasia, and anterior cruciate ligament rupture using whole-genome sequencing

Here, we report the use of genome-wide association study (GWAS) for the analysis of canine whole-genome sequencing (WGS) repository data using breed phenotypes. Single-nucleotide polymorphisms (SNPs) were called from WGS data from 648 dogs that included 119 breeds from the Dog10K Genomes Project. Next, we assigned breed phenotypes for hip dysplasia (Orthopedic Foundation for Animals (OFA) HD, n = 230 dogs from 27 breeds; hospital HD, n = 279 dogs from 38 breeds), elbow dysplasia (ED, n = 230 dogs from 27 breeds), and anterior cruciate ligament rupture (ACL rupture, n = 279 dogs from 38 breeds), the three most important canine spontaneous complex orthopedic diseases. Substantial morbidity is common with these diseases. Previous within- and between-breed GWAS for HD, ED, and ACL rupture using array SNPs have identified disease-associated loci. Individual disease phenotypes are lacking in repository data. There is a critical knowledge gap regarding the optimal approach to undertake categorical GWAS without individual phenotypes. We considered four GWAS approaches: a classical linear mixed model, a haplotype-based model, a binary case-control model, and a weighted least squares model using SNP average allelic frequency. We found that categorical GWAS was able to validate HD candidate loci. Additionally, we discovered novel candidate loci and genes for all three diseases, including FBX025, IL1A, IL1B, COL27A1, SPRED2 (HD), UGDH, FAF1 (ED), TGIF2 (ED & ACL rupture), and IL22, IL26, CSMD1, LDHA, and TNS1 (ACL rupture). Therefore, categorical GWAS of ancestral dog populations may contribute to the understanding of any disease for which breed epidemiological risk data are available, including diseases for which GWAS has not been performed and candidate loci remain elusive.

Balancing at the Borderline of a Breed: A Case Study of the Hungarian Short-Haired Vizsla Dog Breed, Definition of the Breed Profile Using Simple SNP-Based Methods

The aim of this study was to determine the breed boundary of the Hungarian Short-haired Vizsla (HSV) dog breed. Seventy registered purebred HSV dogs were genotyped on approximately 145,000 SNPs. Principal Component Analysis (PCA) and Admixture analysis certified that they belong to the same population. The outer point of the breed demarcation was a single Hungarian Wire-haired Vizsla (HWV) individual, which was the closest animal genetically to the HSV population in the PCA analysis. Three programs were used for the breed assignment calculations, including the widely used GeneClass2.0 software and two additional approaches developed here: the 'PCA-distance' and 'IBS-central' methods. Both new methods calculate a single number that represents how closely a dog fits into the actual reference population. The former approach calculates this number based on the PCA distances from the median of HSV animals. The latter calculates it from identity by state (IBS) data, measuring the distance from a central animal that is the best representative of the breed. Having no mixed-breed dogs with known HSV genome proportion, admixture animals were simulated by using data of HSV and HWV individuals to calibrate the inclusion/exclusion probabilities for the assignment. The numbers generated from these relatively simple calculations can be used by breeders and clubs to keep their populations under genetic supervision.

Intrinsic base substitution patterns in diverse species reveal links to cancer and metabolism

Analyses of large-scale cancer sequencing data have revealed that mutagenic processes can create distinctive patterns of base substitutions, called mutational signatures. Interestingly, mutational patterns resembling some of these signatures can also be observed in normal cells. To determine whether similar patterns exist more generally, we analyzed large data sets of genetic variation, including mutations from 7 model species and single nucleotide polymorphisms in 42 species, totaling >1.9 billion variants. We found that base substitution patterns for most species closely match single base substitution (SBS) mutational signature 5 in the Catalog of Somatic Mutations in Cancer (COSMIC) database. SBS5 is ubiquitous in cancers and also present in normal human cells, suggesting that similar patterns of genetic variation across so many species are likely due to conserved biochemistry. We investigated the mechanistic origins of the SBS5-like mutational pattern in Saccharomyces cerevisiae, and show that translesion DNA synthesis and sugar metabolism are directly linked to this form of mutagenesis. We propose that conserved metabolic processes in cells are coupled to continuous generation of genetic variants, which can be acted upon by selection to drive the evolution of biological entities.

A single-nucleus transcriptomic atlas of the dog hippocampus reveals the potential relationship between specific cell types and domestication

The process of domestication has led to dramatic differences in behavioral traits between domestic dogs and gray wolves. Whole-genome research found that a class of putative positively selected genes were related to various aspects of learning and memory, such as long-term potentiation and long-term depression. In this study, we constructed a single-nucleus transcriptomic atlas of the dog hippocampus to illustrate its cell types, cell lineage and molecular features. Using the transcriptomes of 105 057 nuclei from the hippocampus of a Beagle dog, we identified 26 cell clusters and a putative trajectory of oligodendrocyte development. Comparative analysis revealed a significant convergence between dog differentially expressed genes (DEGs) and putative positively selected genes (PSGs). Forty putative PSGs were DEGs in glutamatergic neurons, especially in Cluster 14, which is related to the regulation of nervous system development. In summary, this study provides a blueprint to understand the cellular mechanism of dog domestication.

Canine Saliva as a Possible Source of Antimicrobial Resistance Genes

While the One Health issues of intensive animal farming are commonly discussed, keeping companion animals is less associated with the interspecies headway of antimicrobial resistance. With the constant advance in veterinary standards, antibiotics are regularly applied in companion animal medicine. Due to the close coexistence of dogs and humans, dog bites and other casual encounters with dog saliva (e.g., licking the owner) are common. According to our metagenome study, based on 26 new generation sequencing canine saliva datasets from 2020 and 2021 reposited in NCBI SRA by The 10,000 Dog Genome Consortium and the Broad Institute within Darwin’s Ark project, canine saliva is rich in bacteria with predictably transferable antimicrobial resistance genes (ARGs). In the genome of potentially pathogenic Bacteroides, Capnocytophaga, Corynebacterium, Fusobacterium, Pasteurella, Porphyromonas, Staphylococcus and Streptococcus species, which are some of the most relevant bacteria in dog bite infections, ARGs against aminoglycosides, carbapenems, cephalosporins, glycylcyclines, lincosamides, macrolides, oxazolidinone, penams, phenicols, pleuromutilins, streptogramins, sulfonamides and tetracyclines could be identified. Several ARGs, including ones against amoxicillin–clavulanate, the most commonly applied antimicrobial agent for dog bites, were predicted to be potentially transferable based on their association with mobile genetic elements (e.g., plasmids, prophages and integrated mobile genetic elements). According to our findings, canine saliva may be a source of transfer for ARG-rich bacteria that can either colonize the human body or transport ARGs to the host bacteriota, and thus can be considered as a risk in the spread of antimicrobial resistance.

Recent, full-length gene retrocopies are common in canids

Gene retrocopies arise from the reverse transcription and insertion into the genome of processed mRNA transcripts. Although many retrocopies have acquired mutations that render them functionally inactive, most mammals retain active LINE-1 sequences capable of producing new retrocopies. New retrocopies, referred to as retro copy number variants (retroCNVs), may not be identified by standard variant calling techniques in high-throughput sequencing data. Although multiple functional FGF4 retroCNVs have been associated with skeletal dysplasias in dogs, the full landscape of canid retroCNVs has not been characterized. Here, retroCNV discovery was performed on a whole-genome sequencing data set of 293 canids from 76 breeds. We identified retroCNV parent genes via the presence of mRNA-specific 30-mers, and then identified retroCNV insertion sites through discordant read analysis. In total, we resolved insertion sites for 1911 retroCNVs from 1179 parent genes, 1236 of which appeared identical to their parent genes. Dogs had on average 54.1 total retroCNVs and 1.4 private retroCNVs. We found evidence of expression in testes for 12% (14/113) of the retroCNVs identified in six Golden Retrievers, including four chimeric transcripts, and 97 retroCNVs also had significantly elevated F ST across dog breeds, possibly indicating selection. We applied our approach to a subset of human genomes and detected an average of 4.2 retroCNVs per sample, highlighting a 13-fold relative increase of retroCNV frequency in dogs. Particularly in canids, retroCNVs are a largely unexplored source of genetic variation that can contribute to genome plasticity and that should be considered when investigating traits and diseases.

Prime editor-mediated correction of a pathogenic mutation in purebred dogs

Canine hip dysplasia (HD) is a multifactorial disease caused by interactions between genetic and environmental factors. HD, which mainly occurs in medium- to large-sized dogs, is a disease that causes severe pain and requires surgical intervention. However, the procedure is not straight-forward, and the only way to ameliorate the situation is to exclude individual dogs with HD from breeding programs. Recently, prime editing (PE), a novel genome editing tool based on the CRISPR-Cas9 system, has been developed and validated in plants and mice. In this study, we successfully corrected a mutation related to HD in Labrador retriever dogs for the first time. We collected cells from a dog diagnosed with HD, corrected the mutation using PE, and generated mutation-corrected dogs by somatic cell nuclear transfer. The results indicate that PE technology can potentially be used as a platform to correct genetic defects in dogs.

Generation of genome-edited dogs by somatic cell nuclear transfer

Background

Canine cloning technology based on somatic cell nuclear transfer (SCNT) combined with genome-editing tools such as CRISPR-Cas9 can be used to correct pathogenic mutations in purebred dogs or to generate animal models of disease.

Results

We constructed a CRISPR-Cas9 vector targeting canine DJ-1. Genome-edited canine fibroblasts were established using vector transfection and antibiotic selection. We performed canine SCNT using genome-edited fibroblasts and successfully generated two genome-edited dogs. Both genome-edited dogs had insertion-deletion mutations at the target locus, and DJ-1 expression was either downregulated or completely repressed.

Conclusion

SCNT successfully produced genome-edited dogs by using the CRISPR-Cas9 system for the first time.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12896-022-00749-3.

[Dogs, wolves and genes]

The most significant polymorphism associated with dog size occurs in the region of the IGF1 gene and concerns a single base change in a neighbouring lncRNA. The "small" (C) allele of this SNP is mostly found in small modern breeds and canids (foxes, coyotes, jackals) while the "large" (T) dominates in wolves and large dogs. However, the small allele is also present at low level in ancient wolves and is shown to represent the ancestral allele in canids, which has been recently selected in small dog breeds obtained by human selection.

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.

Dogs lacking Apolipoprotein E show advanced atherosclerosis leading to apparent clinical complications

Atherosclerotic cardiovascular disease resulting from dysregulated lipid metabolism is the leading cause of morbidity and mortality worldwide. Apolipoprotein E (ApoE) plays a critical role in cholesterol metabolism. Knockouts in lipid-metabolizing proteins including ApoE in multiple model organisms such as mice and rats exhibiting elevated levels of cholesterol have been widely used for dissecting the pathology of atherosclerosis, but few of these animal models exhibit advanced atherosclerotic plaques leading to ischemia-induced clinical symptoms, limiting their use for translational studies. Here we report hypercholesterolemia and severe atherosclerosis characterized by stenosis and occlusion of arteries, together with clinical manifestations of stroke and gangrene, in ApoE knockout dogs generated by CRISPR/Cas9 and cloned by somatic cell nuclear transfer technologies. Importantly, the hypercholesterolemia and atherosclerotic complications in F0 mutants are recapitulated in their offspring. As the ApoE-associated atherosclerosis and clinical manifestations in mutant dogs are more similar to that in human patients compared with those in other animal models, these mutant dogs will be invaluable in developing and evaluating new therapies, including endovascular procedures, against atherosclerosis and related disorders.

Imprinting of canine IGF2 and H19

Genomic imprinting occurs in therian mammals and is a phenomenon whereby the two alleles of a gene are differentially expressed, based on the sex of the parent from whom the alleles were inherited. The allelic differences in expression are the consequence of different epigenetic modifications that are established in the sperm or oocyte during gametogenesis and transmitted at fertilization to offspring. A small minority of genes is regulated in this way but they have important biological functions, and aberrant regulation of imprinted genes contributes to disease aetiology in humans and other animals. The factors driving the evolution of imprinted genes are also of considerable interest, as these genes appear to forego the benefits of diploidy. To broaden the phylogenetic analysis of genomic imprinting, we began a study of imprinted genes in the domestic dog, Canis familiaris. In this report, we show that canine IGF2 and H19 are imprinted, with parent-of origin-dependent monoallelic expression patterns in neonatal umbilical cord. We identify a putative imprint control region associated with the genes, and provide evidence for differential methylation of this region in a somatic tissue (umbilical cord) and for its hypermethylation in the male germline. Canis familiaris is fast becoming a highly informative system for elucidating disease processes and evolution, and the study of imprinted genes in this species may help in understanding how these genes contribute to the generation of morphological and behavioral diversity.

Best practices for analyzing imputed genotypes from low-pass sequencing in dogs

Although DNA array-based approaches for genome-wide association studies (GWAS) permit the collection of thousands of low-cost genotypes, it is often at the expense of resolution and completeness, as SNP chip technologies are ultimately limited by SNPs chosen during array development. An alternative low-cost approach is low-pass whole genome sequencing (WGS) followed by imputation. Rather than relying on high levels of genotype confidence at a set of select loci, low-pass WGS and imputation rely on the combined information from millions of randomly sampled low-confidence genotypes. To investigate low-pass WGS and imputation in the dog, we assessed accuracy and performance by downsampling 97 high-coverage (> 15×) WGS datasets from 51 different breeds to approximately 1× coverage, simulating low-pass WGS. Using a reference panel of 676 dogs from 91 breeds, genotypes were imputed from the downsampled data and compared to a truth set of genotypes generated from high-coverage WGS. Using our truth set, we optimized a variant quality filtering strategy that retained approximately 80% of 14 M imputed sites and lowered the imputation error rate from 3.0% to 1.5%. Seven million sites remained with a MAF > 5% and an average imputation quality score of 0.95. Finally, we simulated the impact of imputation errors on outcomes for case-control GWAS, where small effect sizes were most impacted and medium-to-large effect sizes were minorly impacted. These analyses provide best practice guidelines for study design and data post-processing of low-pass WGS-imputed genotypes in dogs.

Wolves, dogs and humans in regular contact can mutually impact each other's skin microbiota

In contrast to humans and dogs, the skin microbiota of wolves is yet to be described. Here, we investigated the skin microbiota of dogs and wolves kept in outdoor packs at the Wolf Science Center (WSC) via 16S rRNA gene amplicon sequencing. Skin swab samples were also collected from human care takers and their pet dogs. When comparing the three canine groups, representing different degrees of human contact to the care takers and each other, the pet dogs showed the highest level of diversity. Additionally, while human skin was dominated by a few abundant phylotypes, the skin microbiota of the care takers who had particularly close contact with the WSC animals was more similar to the microbiota of dogs and wolves compared to the humans who had less contact with these animals. Our results suggest that domestication may have an impact on the diversity of the skin microbiota, and that the canine skin microbiota can be shared with humans, depending on the level of interaction.

Genome-Wide Association Studies Reveal Neurological Genes for Dog Herding, Predation, Temperament, and Trainability Traits

Genome-wide association study (GWAS) using dog breed standard values as phenotypic measurements is an efficient way to identify genes associated with morphological and behavioral traits. As a result of strong human purposeful selections, several specialized behavioral traits such as herding and hunting have been formed in different modern dog breeds. However, genetic analyses on this topic are rather limited due to the accurate phenotyping difficulty for these complex behavioral traits. Here, 268 dog whole-genome sequences from 130 modern breeds were used to investigate candidate genes underlying dog herding, predation, temperament, and trainability by GWAS. Behavioral phenotypes were obtained from the American Kennel Club based on dog breed standard descriptions or groups (conventional categorization of dog historical roles). The GWAS results of herding behavior (without body size as a covariate) revealed 44 significantly associated sites within five chromosomes. Significantly associated sites on CFA7, 9, 10, and 20 were located either in or near neuropathological or neuronal genes including THOC1, ASIC2, MSRB3, LLPH, RFX8, and CHL1. MSRB3 and CHL1 genes were reported to be associated with dog fear. Since herding is a restricted hunting behavior by removing killing instinct, 36 hounds and 55 herding dogs were used to analyze predation behavior. Three neuronal-related genes (JAK2, MEIS1, and LRRTM4) were revealed as candidates for predation behavior. The significantly associated variant of temperament GWAS was located within ACSS3 gene. The highest associated variant in trainability GWAS is located on CFA22, with no variants detected above the Bonferroni threshold. Since dog behaviors are correlated with body size, we next incorporate body mass as covariates into GWAS; and significant signals around THOC1, MSRB3, LLPH, RFX8, CHL1, LRRTM4, and ACSS3 genes were still detected for dog herding, predation, and temperament behaviors. In humans, these candidate genes are either involved in nervous system development or associated with mental disorders. In conclusion, our results imply that these neuronal or psychiatric genes might be involved in biological processes underlying dog herding, predation, and temperament behavioral traits.

Improving the resolution of canine genome-wide association studies using genotype imputation: A study of two breeds

Genotype imputation using a reference panel that combines high-density array data and publicly available whole genome sequence consortium variant data is potentially a cost-effective method to increase the density of extant lower-density array datasets. In this study, three datasets (two Border Collie; one Italian Spinone) generated using a legacy array (Illumina CanineHD, 173 662 SNPs) were utilised to assess the feasibility and accuracy of this approach and to gather additional evidence for the efficacy of canine genotype imputation. The cosmopolitan reference panels used to impute genotypes comprised dogs of 158 breeds, mixed breed dogs, wolves and Chinese indigenous dogs, as well as breed-specific individuals genotyped using the Axiom Canine HD array. The two Border Collie reference panels comprised 808 individuals including 79 Border Collies and 426 326 or 426 332 SNPs; and the Italian Spinone reference panel comprised 807 individuals including 38 Italian Spinoni and 476 313 SNPs. A high accuracy for imputation was observed, with the lowest accuracy observed for one of the Border Collie datasets (mean R² = 0.94) and the highest for the Italian Spinone dataset (mean R² = 0.97). This study’s findings demonstrate that imputation of a legacy array study set using a reference panel comprising both breed-specific array data and multi-breed variant data derived from whole genomes is effective and accurate. The process of canine genotype imputation, using the valuable growing resource of publicly available canine genome variant datasets alongside breed-specific data, is described in detail to facilitate and encourage use of this technique in canine genetics.

A three-dimensional musculoskeletal model of the dog

The domestic dog is interesting to investigate because of the wide range of body size, body mass, and physique in the many breeds. In the last several years, the number of clinical and biomechanical studies on dog locomotion has increased. However, the relationship between body structure and joint load during locomotion, as well as between joint load and degenerative diseases of the locomotor system (e.g. dysplasia), are not sufficiently understood. Collecting this data through in vivo measurements/records of joint forces and loads on deep/small muscles is complex, invasive, and sometimes unethical. The use of detailed musculoskeletal models may help fill the knowledge gap. We describe here the methods we used to create a detailed musculoskeletal model with 84 degrees of freedom and 134 muscles. Our model has three key-features: three-dimensionality, scalability, and modularity. We tested the validity of the model by identifying forelimb muscle synergies of a walking Beagle. We used inverse dynamics and static optimization to estimate muscle activations based on experimental data. We identified three muscle synergy groups by using hierarchical clustering. The activation patterns predicted from the model exhibit good agreement with experimental data for most of the forelimb muscles. We expect that our model will speed up the analysis of how body size, physique, agility, and disease influence neuronal control and joint loading in dog locomotion.

Dogs as Pets and Pests: Global Patterns of Canine Abundance, Activity, and Health

Dogs (Canis familiaris) were the first domesticated species and, at an estimated population of 1 billion individuals, are globally ubiquitous today. Describing the tremendous morphometric diversity and evolutionary origins of dogs is a scientific endeavor that predates Darwin, yet our interdisciplinary understanding of the species is just beginning. Here, I present global trends in dog abundance, activity, and health. While the human-dog relationship has for millennia been close, it is also complicated. As pets, companion dogs are often treated as family members and constitute the largest sector of the ever-growing >$200 billion USD global pet care industry. As pests, free-roaming dogs are an emerging threat to native species via both predation and nonconsumptive effects (e.g., disturbance, competition for resources, and hybridization). Furthermore, I briefly discuss mounting evidence of dogs as not only infectious disease reservoirs but also as bridges for the transmission of pathogens between wild animals and humans in zoonotic spillover events, triggering intensive dog population management strategies such as culling. Dog mobility across the urban-wildland interface is an important driver for this and other adverse effects of canines on wildlife populations and is an active topic of disease ecologists and conservation biologists. Other canine scientists, including veterinary clinicians and physiologists, study more mechanistic aspects of dog mobility: the comparative kinetics, kinematics, and energetics of dog locomotor health. I outline the prevalent methodological approaches and breed-specific findings within dog activity and health research, then conclude by recognizing promising technologies that are bridging disciplinary gaps in canine science.

Copy number variation underlies complex phenotypes in domestic dog breeds and other canids

Extreme phenotypic diversity, a history of artificial selection, and socioeconomic value make domestic dog breeds a compelling subject for genomic research. Copy number variation (CNV) is known to account for a significant part of inter-individual genomic diversity in other systems. However, a comprehensive genome-wide study of structural variation as it relates to breed-specific phenotypes is lacking. We have generated whole genome CNV maps for more than 300 canids. Our data set extends the canine structural variation landscape to more than 100 dog breeds, including novel variants that cannot be assessed using microarray technologies. We have taken advantage of this data set to perform the first CNV-based genome-wide association study (GWAS) in canids. We identify 96 loci that display copy number differences across breeds, which are statistically associated with a previously compiled set of breed-specific morphometrics and disease susceptibilities. Among these, we highlight the discovery of a long-range interaction involving a CNV near MED13L and TBX3, which could influence breed standard height. Integration of the CNVs with chromatin interactions, long noncoding RNA expression, and single nucleotide variation highlights a subset of specific loci and genes with potential functional relevance and the prospect to explain trait variation between dog breeds.

Identification of common predisposing loci to hematopoietic cancers in four dog breeds

Histiocytic sarcoma (HS) is a rare but aggressive cancer in both humans and dogs. The spontaneous canine model, which has clinical, epidemiological, and histological similarities with human HS and specific breed predispositions, provides a unique opportunity to unravel the genetic basis of this cancer. In this study, we aimed to identify germline risk factors associated with the development of HS in canine-predisposed breeds. We used a methodology that combined several genome-wide association studies in a multi-breed and multi-cancer approach as well as targeted next-generation sequencing, and imputation We combined several dog breeds (Bernese mountain dogs, Rottweilers, flat-coated retrievers, and golden retrievers), and three hematopoietic cancers (HS, lymphoma, and mast cell tumor). Results showed that we not only refined the previously identified HS risk CDKN2A locus, but also identified new loci on canine chromosomes 2, 5, 14, and 20. Capture and targeted sequencing of specific loci suggested the existence of regulatory variants in non-coding regions and methylation mechanisms linked to risk haplotypes, which lead to strong cancer predisposition in specific dog breeds. We also showed that these canine cancer predisposing loci appeared to be due to the additive effect of several risk haplotypes involved in other hematopoietic cancers such as lymphoma or mast cell tumors as well. This illustrates the pleiotropic nature of these canine cancer loci as observed in human oncology, thereby reinforcing the interest of predisposed dog breeds to study cancer initiation and progression.

Because of specific breed structures and artificial selection, pet dogs suffer from numerous genetic diseases, including cancers and represent a unique spontaneous model of human cancers. Here, we focused on histiocytic sarcoma (HS), a rare and highly aggressive cancer in humans. In this study, we have used spontaneous affected and unaffected dogs from three predisposed dog breeds to identify loci involved in HS predisposition. Through genetic analyses, we showed that these canine cancer predispositions are due to the additive effect of several risk haplotypes also involved in the predisposition of other hematopoietic cancers. The corresponding chromosomal regions in humans are involved in the predisposition of several cancers and are also associated with immune traits. This study demonstrates the pleiotropic nature of these canine cancer loci as observed in human oncology, thereby reinforcing the interest of predisposed dog breeds to study mechanisms involved in cancer initiation.

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.