Genomic Analyses Reveal the Influence of Geographic Origin, Migration, and Hybridization on Modern Dog Breed Development

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.

Archetypal Analysis for population genetics

The estimation of genetic clusters using genomic data has application from genome-wide association studies (GWAS) to demographic history to polygenic risk scores (PRS) and is expected to play an important role in the analyses of increasingly diverse, large-scale cohorts. However, existing methods are computationally-intensive, prohibitively so in the case of nationwide biobanks. Here we explore Archetypal Analysis as an efficient, unsupervised approach for identifying genetic clusters and for associating individuals with them. Such unsupervised approaches help avoid conflating socially constructed ethnic labels with genetic clusters by eliminating the need for exogenous training labels. We show that Archetypal Analysis yields similar cluster structure to existing unsupervised methods such as ADMIXTURE and provides interpretative advantages. More importantly, we show that since Archetypal Analysis can be used with lower-dimensional representations of genetic data, significant reductions in computational time and memory requirements are possible. When Archetypal Analysis is run in such a fashion, it takes several orders of magnitude less compute time than the current standard, ADMIXTURE. Finally, we demonstrate uses ranging across datasets from humans to canids.

Eye contact and sociability data suggests that Australian dingoes were never domesticated

Dogs were the first animal to become domesticated by humans, and they represent a classic model system for unraveling the processes of domestication. We compare Australian dingo eye contact and socialization with Basenji and German Shepherd dog (GSD) breeds. Australian dingoes arrived in Australia 5,000-8,000 BP, and there is debate whether they were domesticated before their arrival. The Basenji represents a primitive breed that diverged from the remaining breeds early in the domestication process, while GSDs are a breed dog selected from existing domestic dogs in the late 1800s. We conducted a 4-phase study with unfamiliar and familiar investigators either sitting passively or actively calling each canid. We found 75% of dingoes made eye contact in each phase. In contrast, 86% of Basenjis and 96% of GSDs made eye contact. Dingoes also exhibited shorter eye-gaze duration than breed dogs and did not respond to their name being called actively. Sociability, quantified as a canid coming within 1 m of the experimenter, was lowest for dingoes and highest for GSDs. For sociability duration, dingoes spent less time within 1 m of the experimenter than either breed dog. When compared with previous studies, these data show that the dingo is behaviorally intermediate between wild wolves and Basenji dogs and suggest that it was not domesticated before it arrived in Australia. However, it remains possible that the accumulation of mutations since colonization has obscured historical behaviors, and dingoes now exist in a feralized retamed cycle. Additional morphological and genetic data are required to resolve this conundrum.

Ironing out the wrinkles and folds in the epidemiology of skin fold dermatitis in dog breeds in the UK

Skin fold dermatitis (intertrigo) is an inflammatory process of closely apposing skin surfaces. Extreme conformations towards folded skin in many dog breeds are linked with higher risk. Using anonymised primary-care veterinary data from the VetCompass Programme, this study aimed to report the frequency, demographic risk factors and clinical management for skin fold dermatitis in the UK. Risk factor analysis used random effects multivariable logistic regression modelling. From a study population of 905,553 dogs, the one-year period prevalence in dogs overall was 0.37% (95% CI 0.35–0.39). Diagnosis was supported by laboratory testing in 4.21% cases. Systemic antibiosis was used in 42.30% cases. Compared with crossbreed dogs, the most highly predisposed breeds were English Bulldog (odds ratio [OR] 49.07, 95% CI 37.79–63.70), French Bulldog (OR 25.92, 95% CI 19.62–34.26,) and Pug (OR 16.27, 95% CI 12.20–21.69). The most protected breeds were Yorkshire Terrier (OR 0.14, 95% CI 0.03–0.56), Border Collie (OR 0.31, 95% CI 0.11–0.84), Jack Russell Terrier (OR 0.53, 95% CI 0.30–0.92) and Labrador Retriever (OR 0.57, 95% CI 0.35–0.93). This study adds further evidence to the welfare concerns around high popularity of dog breeds with extreme conformations. The three breeds with by far the highest odds of skin fold dermatitis represent an extreme brachycephalic conformation.

Predicting Dog Phenotypes from Genotypes

We analyze dog genotypes (i.e., positions of dog DNA sequences that often vary between different dogs) in order to predict the corresponding phenotypes (i.e., unique observed characteristics). More specifically, given chromosome data from a dog, we aim to predict the breed, height, and weight. We explore a variety of linear and non-linear classification and regression techniques to accomplish these three tasks. We also investigate the use of a neural network (both in linear and non-linear modes) for breed classification and compare the performance to traditional statistical methods. We show that linear methods generally outperform or match the performance of non-linear methods for breed classification. However, we show that the reverse is true for height and weight regression. Finally, we evaluate the results of all of these methods based on the number of input features used in the analysis. We conduct experiments using different fractions of the full genomic sequences, resulting in input sequences ranging from 20 SNPs to ∼200k SNPs. In doing so, we explore the impact of using a very limited number of SNPs for prediction. Our experiments demonstrate that these phenotypes in dogs can be predicted with as few as 0.5% of randomly selected SNPs (i.e., 992 SNPs) and that dog breeds can be classified with 50% balanced accuracy with as few as 0.02% SNPs (i.e., 40 SNPs).

Reconstructing the history of founder events using genome-wide patterns of allele sharing across individuals

Founder events play a critical role in shaping genetic diversity, fitness and disease risk in a population. Yet our understanding of the prevalence and distribution of founder events in humans and other species remains incomplete, as most existing methods require large sample sizes or phased genomes. Thus, we developed ASCEND that measures the correlation in allele sharing between pairs of individuals across the genome to infer the age and strength of founder events. We show that ASCEND can reliably estimate the parameters of founder events under a range of demographic scenarios. We then apply ASCEND to two species with contrasting evolutionary histories: ~460 worldwide human populations and ~40 modern dog breeds. In humans, we find that over half of the analyzed populations have evidence for recent founder events, associated with geographic isolation, modes of sustenance, or cultural practices such as endogamy. Notably, island populations have lower population sizes than continental groups and most hunter-gatherer, nomadic and indigenous groups have evidence of recent founder events. Many present-day groups––including Native Americans, Oceanians and South Asians––have experienced more extreme founder events than Ashkenazi Jews who have high rates of recessive diseases due their known history of founder events. Using ancient genomes, we show that the strength of founder events differs markedly across geographic regions and time––with three major founder events related to the peopling of Americas and a trend in decreasing strength of founder events in Europe following the Neolithic transition and steppe migrations. In dogs, we estimate extreme founder events in most breeds that occurred in the last 25 generations, concordant with the establishment of many dog breeds during the Victorian times. Our analysis highlights a widespread history of founder events in humans and dogs and elucidates some of the demographic and cultural practices related to these events.

A founder event occurs when small numbers of ancestral individuals give rise to a large fraction of the population. Founder events reduce genetic variation and increase the risk of recessive diseases. Despite their importance in evolutionary and disease studies, we still only have a limited comprehension of their prevalence and properties in humans and other species, as most existing methods require large sample sizes or phased genomes. Here, we present a flexible method, ASCEND, to infer the timing and the strength of founder events that is suitable for sparse datasets with few samples or limited coverage. ASCEND provides reliable estimates across a wide range of demographic scenarios. By applying it to data from two species (humans and dogs), we document a widespread history of recent founder events in both species and provide insights about the demographic processes related to these events. Our analysis helps to identify groups with strong founder events that should be prioritized for future studies as they offer a unique opportunity for biological discovery and reducing disease burden through mapping of recessive disease-causing genes and pathways, as previously shown in studies of Ashkenazi Jews and Finns.

Identification of genes associated with human-canine communication in canine evolution

The dog (Canis familiaris) was the first domesticated animal and hundreds of breeds exist today. During domestication, dogs experienced strong selection for temperament, behaviour, and cognitive ability. However, the genetic basis of these abilities is not well-understood. We focused on ancient dog breeds to investigate breed-related differences in social cognitive abilities. In a problem-solving task, ancient breeds showed a lower tendency to look back at humans than other European breeds. In a two-way object choice task, they showed no differences in correct response rate or ability to read human communicative gestures. We examined gene polymorphisms in oxytocin, oxytocin receptor, melanocortin 2 receptor, and a Williams-Beuren syndrome-related gene (WBSCR17), as candidate genes of dog domestication. The single-nucleotide polymorphisms on melanocortin 2 receptor were related to both tasks, while other polymorphisms were associated with the unsolvable task. This indicates that glucocorticoid functions are involved in the cognitive skills acquired during dog domestication.

DNA methylation clocks for dogs and humans

DNA methylation profiles have been used to develop biomarkers of aging known as epigenetic clocks, which predict chronological age with remarkable accuracy and show promise for inferring health status as an indicator of biological age. Epigenetic clocks were first built to monitor human aging, but their underlying principles appear to be evolutionarily conserved, as they have now been successfully developed for many mammalian species. Here, we describe reliable and highly accurate epigenetic clocks shown to apply to 93 domestic dog breeds. The methylation profiles were generated using the mammalian methylation array, which utilizes DNA sequences that are conserved across all mammalian species. Canine epigenetic clocks were constructed to estimate age and also average time to death. We also present two highly accurate human–dog dual species epigenetic clocks (R = 0.97), which may facilitate the ready translation from canine to human use (or vice versa) of antiaging treatments being developed for longevity and preventive medicine. Finally, epigenome-wide association studies here reveal individual methylation sites that may underlie the inverse relationship between breed weight and lifespan. Overall, we describe robust biomarkers to measure aging and, potentially, health status in canines.

Genetic Origins of the Two Canis lupus familiaris (Dog) Freight Dog Populations

Despite periodic drops in popularity, Arctic sled dogs continue to play a vital role in northern societies, providing both freight transit and recreational race activities. In this study, we selected the Mackenzie River Husky, a freight dog of complex history, and the Chinook, an American Kennel Club recognized freight dog breed whose heritage reportedly overlaps that of the MKRH, for detailed population analysis. We tested each to determine their component breeds and used admixture analysis to ascertain their population structure. We utilized haplotype analysis to identify genomic regions shared between each population and their founding breeds. Our data show that the Alaskan Malamutes and modern Greenland sled dog contributed to both populations, but there are also unexpected contributions from the German Shepherd dog and Collie. We used haplotype analysis to identify genomic regions nearing fixation in population type and identify provocative genes in each region. Finally, in response to recent reports regarding the importance of dietary lipid genes in Arctic dogs, we analyzed 8 such genes in a targeted analysis observing signatures of selection in both populations at the MLXIPL gene loci. These data highlight the genetic routes that breeds of similar function have taken toward their occupation as successful sled dogs.

Carrier Rate and Mutant Allele Frequency of GM1 Gangliosidosis in Miniature Shiba Inus (Mame Shiba): Population Screening of Breeding Dogs in Japan

GM1 gangliosidosis is a progressive, recessive, autosomal, neurodegenerative, lysosomal storage disorder that affects the brain and multiple systemic organs due to an acid β-galactosidase deficiency encoded by the GLB1 gene. This disease occurs in the Shiba Inu breed, which is one of the most popular traditional breeds in Japan, due to the GLB1:c.1649delC (p.P550Rfs*50) mutation. Previous surveys performed of the Shiba Inu population in Japan found a carrier rate of 1.02–2.94%. Currently, a miniature type of the Shiba Inu called “Mame Shiba”, bred via artificial selection to yield smaller individuals, is becoming more popular than the standard Shiba Inu and it is now one of the most popular breeds in Japan and China. The GM1 gangliosidosis mutation has yet to be surveyed in the Mame Shiba population. This study aimed to determine the frequency of the mutant allele and carrier rate of GM1 gangliosidosis in the Mame Shiba breed. Blood samples were collected from 1832 clinically healthy adult Mame Shiba Inus used for breeding across 143 Japanese kennels. The genotyping was performed using a real-time PCR assay. The survey found nine carriers among the Mame Shibas, indicating that the carrier rate and mutant allele frequency were 0.49% and 0.00246, respectively. This study demonstrated that the mutant allele has already been inherited by the Mame Shiba population. There is a risk of GM1 gangliosidosis occurrence in the Mame Shiba breed if breeders use carriers for mating. Further genotyping surveys are necessary for breeding Mame Shibas to prevent the inheritance of this disease.

The effects of age, sex, weight, and breed on canid methylomes

Unlike genomes, which are static throughout the lifespan of an organism, DNA methylomes are dynamic. To study these dynamics, we developed quantitative models that measure the effect of multiple factors on DNA methylomes including, age, sex, weight, and genetics. We conducted our study in canids, which prove to be an ideal species to assess epigenetic moderators due to their extreme variability in size and well-characterized genetic structure. We collected buccal swabs from 217 canids (207 domestic dogs and 10 grey wolves) and used targeted bisulphite sequencing to measure methylomes. We also measured genotypes at over one thousand single nucleotide polymorphisms (SNPs). As expected, we found that DNA methylomes are strongly associated with age, enabling the construction of epigenetic clocks. However, we also identify novel associations between methylomes and sex, weight, and sterilization status, leading to accurate models that predict these factors. Methylomes are also affected by genetics, and we observe multiple associations between SNP loci and methylated CpGs. Finally, we show that several factors moderate the relationship between epigenetic ages and real ages, such as body weight, which increases epigenetic ageing. In conclusion, we demonstrate that the plasticity of DNA methylomes is impacted by myriad genetics and physiological factors, and that DNA methylation biomarkers are accurate predictors of age, sex and sterilization status.

Ancestry-inclusive dog genomics challenges popular breed stereotypes

Behavioral genetics in dogs has focused on modern breeds, which are isolated subgroups with distinctive physical and, purportedly, behavioral characteristics. We interrogated breed stereotypes by surveying owners of 18,385 purebred and mixed-breed dogs and genotyping 2155 dogs. Most behavioral traits are heritable [heritability (h²) > 25%], and admixture patterns in mixed-breed dogs reveal breed propensities. Breed explains just 9% of behavioral variation in individuals. Genome-wide association analyses identify 11 loci that are significantly associated with behavior, and characteristic breed behaviors exhibit genetic complexity. Behavioral loci are not unusually differentiated in breeds, but breed propensities align, albeit weakly, with ancestral function. We propose that behaviors perceived as characteristic of modern breeds derive from thousands of years of polygenic adaptation that predates breed formation, with modern breeds distinguished primarily by aesthetic traits.

The Australian dingo is an early offshoot of modern breed dogs

Dogs are uniquely associated with human dispersal and bring transformational insight into the domestication process. Dingoes represent an intriguing case within canine evolution being geographically isolated for thousands of years. Here, we present a high-quality de novo assembly of a pure dingo (CanFam_DDS). We identified large chromosomal differences relative to the current dog reference (CanFam3.1) and confirmed no expanded pancreatic amylase gene as found in breed dogs. Phylogenetic analyses using variant pairwise matrices show that the dingo is distinct from five breed dogs with 100% bootstrap support when using Greenland wolf as the outgroup. Functionally, we observe differences in methylation patterns between the dingo and German shepherd dog genomes and differences in serum biochemistry and microbiome makeup. Our results suggest that distinct demographic and environmental conditions have shaped the dingo genome. In contrast, artificial human selection has likely shaped the genomes of domestic breed dogs after divergence from the dingo.

Patagonian sheepdog: Genomic analyses trace the footprints of extinct UK herding dogs to South America

Most modern dog breeds were developed within the last two hundred years, following strong and recent human selection based predominantly on aesthetics, with few modern breeds constructed solely to maximize their work potential. In many cases, these working breeds represent the last remnants of now lost populations. The Patagonian sheepdog (PGOD), a rare herding breed, is a remarkable example of such a population. Maintained as an isolated population for over 130 years, the PGOD offers a unique opportunity to understand the genetic relationship amongst modern herding breeds, determine key genomic structure of the founder PGOD populations, and investigate how canine genomic data can mirror human migration patterns. We thus analyzed the population structure of 159 PGOD, comparing them with 1514 dogs representing 175 established breeds. Using 150,069 SNPs from a high-density SNP genotyping array, we establish the genomic composition, ancestry, and genetic diversity of the population, complementing genomic data with the PGOD's migratory history to South America. Our phylogenetic analysis reveals that PGODs are most closely related to modern herding breeds hailing from the United Kingdom. Admixture models illustrate a greater degree of diversity and genetic heterogeneity within the very small PGOD population than in Western European herding breeds, suggesting the PGOD predates the 200-year-old construction of most pure breeds known today. We thus propose that PGODs originated from the foundational herding dogs of the UK, prior to the Victorian explosion of breeds, and that they are the closest link to a now-extinct population of herding dogs from which modern herding breeds descended.

The odds of neoplasia in dogs with and without diabetes mellitus

BACKGROUND

Increased risk of neoplasia in humans with diabetes mellitus (DM) is well documented. It is unknown if dogs with DM have increased risk of neoplasia.

OBJECTIVE

Determine if dogs with DM have an overall increased risk of neoplasia and risk for specific forms of neoplasia compared to dogs without DM.

ANIMALS

Seven hundred dogs with DM and 700 breed, age, and sex-matched dogs without DM, examined during the same years.

METHODS

Retrospective case-control study. Odds ratios (OR), corresponding 95% confidence intervals (CI), and P-values were calculated using conditional logistic regression to determine if dogs with DM had increased odds of developing neoplasia compared to dogs without DM.

RESULTS

The overall odds of developing neoplasia were not significantly different in dogs with and without DM. However, dogs with DM had significantly higher odds of developing an adrenal mass (OR, 4; 95% CI, 1.1-14.2; P = .03) compared to dogs without DM. The odds of developing a splenic mass in dogs with DM (OR, 1.2; 95% CI, 0.99-1.39) were increased compared to dogs without DM, but this difference was not significant (P = .07).

CONCLUSIONS AND CLINICAL IMPORTANCE

Dogs with DM may be at increased risk for adrenal neoplasia. Awareness of this risk can facilitate early diagnosis of this life-threatening comorbidity. Larger studies are needed to confirm these findings.

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.

Personality traits associate with behavioral problems in pet dogs

Personality traits, especially neuroticism, strongly predict psychopathology. The domestic dog (Canis lupus familiaris Linnaeus, 1758) is used as a natural model for psychiatric disorders, but the similarity between dog and human personality and the association between dog personality and unwanted behavioral traits, such as fearfulness, aggressiveness, and impulsivity/inattention, remain unknown. This study utilized structural equation modeling (SEM) with survey data of 11,360 dogs to examine the associations and correlations between seven personality and ten unwanted behavioral traits. Personality traits included insecurity, energy, training focus, aggressiveness/dominance, human sociability, dog sociability, and perseverance. Unwanted behavioral traits included fearfulness, noise sensitivity, fear of surfaces/heights, separation anxiety, barking, stranger-directed aggression, owner-directed aggression, dog-directed aggression, hyperactivity/impulsivity, and inattention. We first fitted confirmatory factor models for the unwanted behavioral traits and the best model grouped unwanted behaviors into four latent traits: fear-related behavior, fear-aggression, aggression, and impulsivity/inattention and used this structure in the subsequent SEM model. Especially, insecurity, which resembles the human neuroticism trait, was strongly associated with unwanted behavior, paralleling the association between neuroticism and psychopathology. Similarly, training focus, resembling conscientiousness, was negatively related to impulsivity/inattention, and aggressiveness/dominance was associated with aggressive behaviors, resembling associations of conscientiousness and agreeableness with attention deficit hyperactivity disorder and aggression-related psychopathology, respectively. These results indicate that dog personality traits resemble human personality traits, suggesting that their neurological and genetic basis may also be similar and making the dog a suitable animal model for human behavior and psychiatric disorders.

Genetic testing of dogs predicts problem behaviors in clinical and nonclinical samples

BACKGROUND

Very little is known about the etiology of personality and psychiatric disorders. Because the core neurobiology of many such traits is evolutionarily conserved, dogs present a powerful model. We previously reported genome scans of breed averages of ten traits related to fear, anxiety, aggression and social behavior in multiple cohorts of pedigree dogs. As a second phase of that discovery, here we tested the ability of markers at 13 of those loci to predict canine behavior in a community sample of 397 pedigree and mixed-breed dogs with individual-level genotype and phenotype data.

RESULTS

We found support for all markers and loci. By including 122 dogs with veterinary behavioral diagnoses in our cohort, we were able to identify eight loci associated with those diagnoses. Logistic regression models showed subsets of those loci could predict behavioral diagnoses. We corroborated our previous findings that small body size is associated with many problem behaviors and large body size is associated with increased trainability. Children in the home were associated with anxiety traits; illness and other animals in the home with coprophagia; working-dog status with increased energy and separation-related problems; and competitive dogs with increased aggression directed at familiar dogs, but reduced fear directed at humans and unfamiliar dogs. Compared to other dogs, Pit Bull-type dogs were not defined by a set of our markers and were not more aggressive; but they were strongly associated with pulling on the leash. Using severity-threshold models, Pit Bull-type dogs showed reduced risk of owner-directed aggression (75th quantile) and increased risk of dog-directed fear (95th quantile).

CONCLUSIONS

Our association analysis in a community sample of pedigree and mixed-breed dogs supports the interbreed mapping. The modeling shows some markers are predictive of behavioral diagnoses. Our findings have broad utility, including for clinical and breeding purposes, but we caution that thorough understanding is necessary for their interpretation and use.

Myxomatous mitral valve disease in Miniature Schnauzers and Yorkshire Terriers: 134 cases (2007-2016)

OBJECTIVE

To characterize features of myxomatous mitral valve disease (MMVD) in Miniature Schnauzers and Yorkshire Terriers.

ANIMALS

69 Miniature Schnauzers and 65 Yorkshire Terriers, each with MMVD.

PROCEDURES

Medical record data for each dog were collected; the study period was January 2007 through December 2016. If available, radiographic data were evaluated, and a vertebral heart scale score was assigned for each dog. Statistical analysis was performed with Student t and Fisher exact tests.

RESULTS

Compared with Yorkshire Terriers, the prevalence of MMVD was significantly higher in Miniature Schnauzers and affected dogs were significantly younger at the time of diagnosis. Miniature Schnauzers were significantly more likely to have mitral valve prolapse and syncope, compared with Yorkshire Terriers. Yorkshire Terriers were significantly more likely to have coughing and have had previous or current treatment with cardiac medications, compared with Miniature Schnauzers. There was no statistical difference between breeds with regard to abnormally high vertebral heart scale scores or radiographic evidence of congestive heart failure.

CONCLUSIONS AND CLINICAL RELEVANCE

With regard to MMVD, features of the disease among Miniature Schnauzers and Yorkshire Terriers were similar, but there were also a few discernable differences between these 2 breeds and from historical findings for dogs with MMVD of other breeds. Clinical signs at the time of diagnosis differed between the 2 breeds, which may have reflected concurrent breed-specific conditions (sick sinus syndrome or airway disease [eg, tracheal collapse]). Future work should include prospective studies to provide additional information regarding the natural progression of MMVD in these dog breeds.

Heritability and Genomic Architecture of Episodic Exercise-Induced Collapse in Border Collies

An episodic nervous system disorder triggered by strenuous exercise, termed border collie collapse (BCC), exists in border collies and related breeds. The genetic basis of BCC is unknown but is believed to be a complex genetic disorder. Our goal was to estimate the heritability (h²_SNP) of BCC, define its underlying genetic architecture, and identify associated genomic loci using dense whole-genome single-nucleotide polymorphism (SNP) genotyping data. Genotype data were obtained for ~440,000 SNPs from 343 border collies (168 BCC cases and 175 controls). h²_SNP was calculated to be 49–61% depending on the estimated BCC prevalence. A total of 2407 SNPs across the genome accounted for nearly all the h2_SNP of BCC, with an estimated 2003 SNPs of small effect, 349 SNPs of moderate effect, and 56 SNPs of large effect. Genome-wide association analyses identified significantly associated loci on chromosomes 1, 6, 11, 20, and 28, which accounted for ~5% of the total BCC h²_SNP. We conclude that BCC is a moderately- to highly-heritable complex polygenetic disease resulting from contributions from hundreds to thousands of genetic variants with variable effect sizes. Understanding how much the BCC phenotype is determined by genetics and whether major gene mutations are likely to exist inform veterinarians and working/stock dog communities of the true nature of this condition.

Breed group differences in the unsolvable problem task: herding dogs prefer their owner, while solitary hunting dogs seek stranger proximity

The communicating skills of dogs are well documented and especially their contact-seeking behaviours towards humans. The aim of this study was to use the unsolvable problem paradigm to investigate differences between breed groups in their contact-seeking behaviours towards their owner and a stranger. Twenty-four dogs of ancient breeds, 58 herding dogs, and 17 solitary hunting dogs were included in the study, and their behaviour when presented with an unsolvable problem task (UPT) was recorded for 3 min. All breed groups interacted with the test apparatus the same amount of time, but the herding dogs showed a longer gaze duration towards their owner compared to the other groups and they also preferred to interact with their owner instead of a stranger. Interestingly, the solitary hunting dogs were more in stranger proximity than the other groups, and they also showed a preference to make contact with a stranger instead of their owner. Hence, we found differences in contact-seeking behaviours, reflecting the dog–human relationship, between breed groups that might not only be related to their genetic similarity to wolves, but also due to the specific breeding history of the dogs.

An aggressive CD4- CD8- T-cell neoplasm in young English bulldogs

T-cell leukemia/lymphoma accounts for roughly 30% of all types of lymphoproliferative neoplasia in dogs. Two forms of T-cell lymphoma (T-zone and peripheral T-cell lymphoma) exhibit breed-specific predilections. During the course of routine immunophenotyping, we observed a breed-specific presentation of a unique form of T-cell leukaemia in young English bulldogs. To describe the clinical presentation and outcome of a novel T-cell leukaemia in English bulldogs and determine the frequency of this neoplasm in other breeds. The Clinical Hematopathology database, containing immunophenotyping data from peripheral blood of nearly 11 900 dogs, was queried for the phenotype observed in young English bulldogs: CD45⁺ CD4^- CD8^- CD5⁺ CD3⁺ class II major histocompatibility complex (MHC)-low T-cell leukaemia. Clinical presentation, treatment, and survival data were collected for a subset of cases. Fifty-five English bulldog cases and 64 cases of other breeds were identified. No other breed was represented by >5 cases. Complete medical records were obtained for 50 bulldogs. Median age at diagnosis was 3 years and 76% of cases were male. Median lymphocyte count was 44 286 lymphocytes/μl (range, 1800-317 684/μl) and lymphocytes were described as small to intermediate-sized. Many dogs were thrombocytopenic and had liver and spleen involvement, but not lymphadenopathy. Bulldogs that received multi-agent chemotherapy had longer median survival times (83 days) compared to dogs that received no treatment (6 days) or less aggressive therapy (15 days) (p = .001). Non-bulldogs had similar outcomes. CD4^- CD8^- class II MHC-low T-cell leukaemia has an aggressive clinical course and predilection for young English bulldogs. Breed-specific presentation suggests an underlying genetic cause.

The evolution and ecology of benign tumors

Tumors are usually classified into two main categories - benign or malignant, with much more attention being devoted to the second category given that they are usually associated with more severe health issues (i.e., metastatic cancers). Here, we argue that the mechanistic distinction between benign and malignant tumors has narrowed our understanding of neoplastic processes. This review provides the first comprehensive discussion of benign tumors in the context of their evolution and ecology as well as interactions with their hosts. We compare the genetic and epigenetic profiles, cellular activities, and the involvement of viruses in benign and malignant tumors. We also address the impact of intra-tumoral cell composition and its relationship with the tumoral microenvironment. Lastly, we explore the differences in the distribution of benign and malignant neoplasia across the tree of life and provide examples on how benign tumors can also affect individual fitness and consequently the evolutionary trajectories of populations and species. Overall, our goal is to bring attention to the non-cancerous manifestations of tumors, at different scales, and to stimulate research on the evolutionary ecology of host-tumor interactions on a broader scale. Ultimately, we suggest that a better appreciation of the differences and similarities between benign and malignant tumors is fundamental to our understanding of malignancy both at mechanistic and evolutionary levels.

Panobinostat Effectively Increases Histone Acetylation and Alters Chromatin Accessibility Landscape in Canine Embryonic Fibroblasts but Does Not Enhance Cellular Reprogramming

Robust and reproducible protocols to efficiently reprogram adult canine cells to induced pluripotent stem cells are still elusive. Somatic cell reprogramming requires global chromatin remodeling that is finely orchestrated spatially and temporally. Histone acetylation and deacetylation are key regulators of chromatin condensation, mediated by histone acetyltransferases and histone deacetylases (HDACs), respectively. HDAC inhibitors have been used to increase histone acetylation, chromatin accessibility, and somatic cell reprogramming in human and mice cells. We hypothesized that inhibition of HDACs in canine fibroblasts would increase their reprogramming efficiency by altering the epigenomic landscape and enabling greater chromatin accessibility. We report that a combined treatment of panobinostat (LBH589) and vitamin C effectively inhibits HDAC function and increases histone acetylation in canine embryonic fibroblasts in vitro, with no significant cytotoxic effects. We further determined the effect of this treatment on global chromatin accessibility via Assay for Transposase-Accessible Chromatin using sequencing. Finally, the treatment did not induce any significant increase in cellular reprogramming efficiency. Although our data demonstrate that the unique epigenetic landscape of canine cells does not make them amenable to cellular reprogramming through the proposed treatment, it provides a rationale for a targeted, canine-specific, reprogramming approach by enhancing the expression of transcription factors such as CEBP.

OptM: estimating the optimal number of migration edges on population trees using Treemix

The software Treemix has become extensively used to estimate the number of migration events, or edges (m), on population trees from genome-wide allele frequency data. However, the appropriate number of edges to include remains unclear. Here, I show that an optimal value of m can be inferred from the second-order rate of change in likelihood (Δm) across incremental values of m. Repurposed from its original use to estimate the number of population clusters in the software Structure (ΔK), I show using simulated populations that Δm performs equally as well as current recommendations for Treemix. A demonstration of an empirical dataset from domestic dogs indicates that this method may be preferable in large, complex population histories and can prioritize migration events for subsequent investigation. The method has been implemented in a freely available R package called “OptM” and as a web application (https://rfitak.shinyapps.io/OptM/) to interface directly with the output files of Treemix.

Description of breed ancestry and genetic health traits in arctic sled dog breeds

Background

This study describes the presence and frequency of health traits among three populations of dogs traditionally used for sledding and explores their ancestry and breed composition as provided by the commercially available Embark dog DNA test. The three populations include the purebred Siberian Husky and the admixed populations of Alaskan sled dogs and Polar Huskies. While the Siberian Husky represents a well-established breed with extensive historical and health data, the Alaskan sled dog is less studied but has been the subject of nutritional, physiological, and genetic studies related to ancestry and performance. In contrast, the Polar Husky is a relatively obscure and rare group of dogs used for arctic exploration with very little-known information. The three populations were compared using Embark results, providing new insight into the health traits circulating within the populations and the potential ancestral linkage of the health traits between the sledding populations. Embark results are based upon 228,588 single-nucleotide polymorphisms (SNPs) spanning the canine genome, characterized using a custom-designed Illumina beadchip array.

Results

Specifically, breed composition was summarized for the two admixed populations with most of the dogs being predominantly categorized as Alaskan husky- type dog or “Supermutt”. Mitochondrial and Y chromosome haplogroups and haplotypes were found with Alaskan sled dogs carrying most of the haplogroups and types found in Siberian and Polar Huskies. Genomic principal component analysis reflected population structure corresponding to breed and substructure within the Alaskan sled dogs related to sprint or distance competition. Genetic markers associated with Alanine Aminotransferase activity, Alaskan Husky Encephalopathy, dilated cardiomyopathy, Collie eye anomaly, degenerative myelopathy, ichthyosis, and factor VII deficiency were identified in the populations of sledding breeds.

Conclusion

These results provide a preliminary description of genetic characteristics found in sledding breeds, improving the understanding and care of working sled dogs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40575-021-00108-z.

Did Dog Domestication Contribute to Language Evolution?

Different factors seemingly account for the emergence of present-day languages in our species. Human self-domestication has been recently invoked as one important force favoring language complexity mostly via a cultural mechanism. Because our self-domestication ultimately resulted from selection for less aggressive behavior and increased prosocial behavior, any evolutionary or cultural change impacting on aggression levels is expected to have fostered this process. Here, we hypothesize about a parallel domestication of humans and dogs, and more specifically, about a positive effect of our interaction with dogs on human self-domestication, and ultimately, on aspects of language evolution, through the mechanisms involved in the control of aggression. We review evidence of diverse sort (ethological mostly, but also archeological, genetic, and physiological) supporting such an effect and propose some ways of testing our hypothesis.

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.

The genetic consequences of dog breed formation-Accumulation of deleterious genetic variation and fixation of mutations associated with myxomatous mitral valve disease in cavalier King Charles spaniels

Selective breeding for desirable traits in strictly controlled populations has generated an extraordinary diversity in canine morphology and behaviour, but has also led to loss of genetic variation and random entrapment of disease alleles. As a consequence, specific diseases are now prevalent in certain breeds, but whether the recent breeding practice led to an overall increase in genetic load remains unclear. Here we generate whole genome sequencing (WGS) data from 20 dogs per breed from eight breeds and document a ~10% rise in the number of derived alleles per genome at evolutionarily conserved sites in the heavily bottlenecked cavalier King Charles spaniel breed (cKCs) relative to in most breeds studied here. Our finding represents the first clear indication of a relative increase in levels of deleterious genetic variation in a specific breed, arguing that recent breeding practices probably were associated with an accumulation of genetic load in dogs. We then use the WGS data to identify candidate risk alleles for the most common cause for veterinary care in cKCs–the heart disease myxomatous mitral valve disease (MMVD). We verify a potential link to MMVD for candidate variants near the heart specific NEBL gene in a dachshund population and show that two of the NEBL candidate variants have regulatory potential in heart-derived cell lines and are associated with reduced NEBL isoform nebulette expression in papillary muscle (but not in mitral valve, nor in left ventricular wall). Alleles linked to reduced nebulette expression may hence predispose cKCs and other breeds to MMVD via loss of papillary muscle integrity.

As a consequence of selective breeding, specific disease-causing mutations have become more frequent in certain dog breeds. Whether the breeding practice also resulted in a general increase in the overall number of disease-causing mutations per dog genome is however not clear. To address this question, we compare the amount of harmful, potentially disease-causing, mutations in dogs from eight common breeds that have experienced varying degrees of intense selective breeding. We find that individuals belonging to the breed affected by the most intense breeding—cavalier King Charles spaniel (cKCs)—carry more harmful variants than other breeds, indicating that past breeding practices may have increased the overall levels of harmful genetic variation in dogs. The most common disease in cKCs is myxomatous mitral valve disease (MMVD). To identify variants linked to this disease we next characterize mutations that are common in cKCs, but rare in other breeds, and then investigate if these mutations can predict MMVD in dachshunds. We find that variants that regulate the expression of the gene NEBL in papillary muscles may increase the risk of the disease, indicating that loss of papillary muscle integrity could contribute to the development of MMVD.

Is the Prevalence of Leishmania infantum Linked to Breeds in Dogs? Characterization of Seropositive Dogs in Ibiza

Simple Summary

Leishmaniosis is an important zoonotic protozoan disease. Leishmania infantum is a protozoan species that accounts for the majority of cases in the Mediterranean. In this study, we analyzed the prevalence of infection in different dog breeds from Ibizan Island. Our results showed that the Doberman Pinscher and Boxer breeds present a higher prevalence of infection, and the relationship between antibodies’ serum titer and staging of disease was confirmed. Differences between age and sex were not found.

Abstract

Leishmaniosis is an important zoonotic protozoan disease primarily spread to the Mediterranean region by Leishmania infantum, the predominant protozoan species, which accounts for the majority of cases. Development of disease depends on the immune response of the definitive host and, predictably, their genetic background. Recent studies have revealed breed-typical haplotypes that are susceptible to the spread of the protozoan parasite. The objective of this study was to analyze the prevalence of leishmaniosis on a Mediterranean island and determine the relationship between disease prevalence and breed. In addition, information on seropositive animals was recorded to characterize animals affected by the disease. To study the prevalence, a total of 3141 dogs were analyzed. Of these, the 149 infected animals were examined for age, sex, antibody titer, and disease stage. We observed a prevalence of 4.74%, which varied between breeds (p < 0.05). The Doberman Pinscher and Boxer breeds had the highest prevalence of leishmaniosis. Significant differences were observed between breeds with common ancestors, emphasizing the important genetic component. Finally, regarding the characterization of seropositive animals, the distribution is similar to other studies. We discovered a relationship (p < 0.05) between the number of antibody titers and the clinical disease stage, which was also present in Leishmania infantum, suggesting that the development of the disease depends on the humoral or Th2 immune response with ineffective antibodies.

Neurodevelopmental scaling is a major driver of brain-behavior differences in temperament across dog breeds

Behavioral traits like aggression, anxiety, and trainability differ significantly across dog breeds and are highly heritable. However, the neural bases of these differences are unknown. Here we analyzed structural MRI scans of 62 dogs in relation to breed-average scores for the 14 major dimensions in the Canine Behavioral Assessment and Research Questionnaire, a well-validated measure of canine temperament. Several behavior categories showed significant relationships with morphologically covarying gray matter networks and regional volume changes. Networks involved in social processing and the flight-or-fight response were associated with stranger-directed fear and aggression, putatively the main behaviors under selection pressure during wolf-to-dog domestication. Trainability was significantly associated with expansion in broad regions of cortex, while fear, aggression, and other "problem" behaviors were associated with expansion in distributed subcortical regions. These results closely overlapped with regional volume changes with total brain size, in striking correspondence with models of developmental constraint on brain evolution. This suggests that the established link between dog body size and behavior is due at least in part to disproportionate enlargement of later-developing regions in larger brained dogs. We discuss how this may explain the known correlation of increasing reactivity with decreasing body size in dogs.

Exceptional Changes in Skeletal Anatomy under Domestication: The Case of Brachycephaly

"Brachycephaly" is generally considered a phenotype in which the facial part of the head is pronouncedly shortened. While brachycephaly is characteristic for some domestic varieties and breeds (e.g., Bulldog, Persian cat, Niata cattle, Anglo-Nubian goat, Middle White pig), this phenotype can also be considered pathological. Despite the superficially similar appearance of "brachycephaly" in such varieties and breeds, closer examination reveals that "brachycephaly" includes a variety of different cranial modifications with likely different genetic and developmental underpinnings and related with specific breed histories. We review the various definitions and characteristics associated with brachycephaly in different domesticated species. We discern different types of brachycephaly ("bulldog-type," "katantognathic," and "allometric" brachycephaly) and discuss morphological conditions related to brachycephaly, including diseases (e.g., brachycephalic airway obstructive syndrome). Further, we examine the complex underlying genetic and developmental processes and the culturally and developmentally related reasons why brachycephalic varieties may or may not be prevalent in certain domesticated species. Knowledge on patterns and mechanisms associated with brachycephaly is relevant for domestication research, veterinary and human medicine, as well as evolutionary biology, and highlights the profound influence of artificial selection by humans on animal morphology, evolution, and welfare.

Oral bacterial DNA-based discrimination of human and canine saliva for the analysis of indistinct bite marks

Analyzing ambiguous bite marks using conventional morphological approaches to identify attackers is difficult; thus, applying molecular biological methods for identifying an attacker from their saliva is a possible approach in a forensic investigation. This study aimed to establish oral bacterial DNA-based human and canine saliva markers and develop a practical method for their discrimination. We considered Streptococcus oralis and Pasteurella canis as human and canine saliva marker candidates, respectively. Duplex bacterial DNA detection using melting curve analysis was designed and evaluated for forensic applicability using proof-of-concept experiments. S. oralis DNA was detected from human saliva samples from 30 out of 30 individuals, and P. canis DNA was detected from canine saliva samples from 73 out of 77 individuals (26 dog breeds). Additionally, both bacterial DNA markers were accurately detected from human blood-contaminated saliva samples and mock indistinct bite marks. Our results indicate that both bacterial DNA markers were sensitive, robust, and discriminating saliva markers. We consider that our duplex bacterial DNA examination is a simple, practical, and useful method for the detection of saliva from indistinct bite marks and discrimination between human and canine saliva.

The Emergence of SARS-CoV-2 within the Dog Population in Croatia: Host Factors and Clinical Outcome

Over a year into the COVID-19 pandemic, there is growing evidence that SARS-CoV-2 infections among dogs are more common than previously thought. In this study, the prevalence of SARS-CoV-2 antibodies was investigated in two dog populations. The first group was comprised of 1069 dogs admitted to the Veterinary Teaching Hospital for any given reason. The second group included dogs that shared households with confirmed COVID-19 cases in humans. This study group numbered 78 dogs. In COVID-19 infected households, 43.9% tested ELISA positive, and neutralising antibodies were detected in 25.64% of dogs. Those data are comparable with the secondary attack rate in the human population. With 14.69% of dogs in the general population testing ELISA positive, there was a surge of SARS-CoV-2 infections within the dog population amid the second wave of the pandemic. Noticeably seroprevalence of SARS-CoV-2 in the dog and the human population did not differ at the end of the study period. Male sex, breed and age were identified as significant risk factors. This study gives strong evidence that while acute dog infections are mostly asymptomatic, they can pose a significant risk to dog health. Due to the retrospective nature of this study, samples for viral isolation and PCR were unavailable. Still, seropositive dogs had a 1.97 times greater risk for developing central nervous symptoms.

Whole Genome Sequencing Reveals Signatures for Artificial Selection for Different Sizes in Japanese Primitive Dog Breeds

Body size is an important trait in companion animals. Recently, a primitive Japanese dog breed, the Shiba Inu, has experienced artificial selection for smaller body size, resulting in the "Mame Shiba Inu" breed. To identify loci and genes that might explain the difference in the body size of these Shiba Inu dogs, we applied whole genome sequencing of pooled samples (pool-seq) on both Shiba Inu and Mame Shiba Inu. We identified a total of 13,618,261 unique SNPs in the genomes of these two breeds of dog. Using selective sweep approaches, including F _ST, H _p and XP-CLR with sliding windows, we identified a total of 12 genomic windows that show signatures of selection that overlap with nine genes (PRDM16, ZNF382, ZNF461, ERGIC2, ENSCAFG00000033351, CCDC61, ALDH3A2, ENSCAFG00000011141, and ENSCAFG00000018533). These results provide candidate genes and specific sites that might be associated with body size in dogs. Some of these genes are associated with body size in other mammals, but 8 of the 9 genes are novel candidate genes that need further study.

Clinical features, diagnosis, and survival analysis of dogs with glioma

BACKGROUND

Gliomas in dogs remain poorly understood.

OBJECTIVES

To characterize the clinicopathologic findings, diagnostic imaging features and survival of a large sample of dogs with glioma using the Comparative Brain Tumor Consortium diagnostic classification.

ANIMALS

Ninety-one dogs with histopathological diagnosis of glioma.

METHODS

Multicentric retrospective case series. Signalment, clinicopathologic findings, diagnostic imaging characteristics, treatment, and outcome were used. Tumors were reclassified according to the new canine glioma diagnostic scheme.

RESULTS

No associations were found between clinicopathologic findings or survival and tumor type or grade. However, definitive treatments provided significantly (P = .03) improved median survival time (84 days; 95% confidence interval [CI], 45-190) compared to palliative treatment (26 days; 95% CI, 11-54). On magnetic resonance imaging (MRI), oligodendrogliomas were associated with smooth margins and T1-weighted hypointensity compared to astrocytomas (odds ratio [OR], 42.5; 95% CI, 2.42-744.97; P = .04; OR, 45.5; 95% CI, 5.78-333.33; P < .001, respectively) and undefined gliomas (OR, 84; 95% CI, 3.43-999.99; P = .02; OR, 32.3; 95% CI, 2.51-500.00; P = .008, respectively) and were more commonly in contact with the ventricles than astrocytomas (OR, 7.47; 95% CI, 1.03-53.95; P = .049). Tumor spread to neighboring brain structures was associated with high-grade glioma (OR, 6.02; 95% CI, 1.06-34.48; P = .04).

CONCLUSIONS AND CLINICAL IMPORTANCE

Dogs with gliomas have poor outcomes, but risk factors identified in survival analysis inform prognosis and the newly identified MRI characteristics could refine diagnosis of tumor type and grade.

Malassezia otitis unresponsive to primary care: outcome in 59 dogs

BACKGROUND

Otitis externa (OE) is a common disorder in dogs. Infection by the commensal yeast, Malassezia pachydermatis, may result in chronic disease that does not respond to standard primary care. Chronic infectious OE may be associated with otitis media (OM).

HYPOTHESIS/OBJECTIVE

To report medical management, clinical outcomes and frequency of middle ear involvement, in dogs with Malassezia otitis unresponsive to primary care.

ANIMALS

Fifty-nine dogs from one referral veterinary hospital from January 2007 to September 2018.

METHODS AND MATERIALS

Retrospective analysis of medical records of dogs referred with chronic otitis and treated for Malassezia otitis at a referral veterinary hospital.

RESULTS

Chronic Malassezia OE was treated successfully in 91% of ears, in 87% of these cases with one ear flush intervention. Median time-to-resolution was 27 days after ear flush intervention. Neither duration of otitis, presence of neutrophils in aural discharge nor administration of oral itraconazole affected clinical outcome. Malassezia OM occurred concurrently in 17% of ears.

CONCLUSIONS AND CLINICAL RELEVANCE

These findings assist clinicians and carers of affected dogs in decision-making, by documenting that most cases of canine Malassezia otitis that have not resolved with standard primary care, can be treated successfully with a well-staged and intense medical treatment plan. Malassezia OM should be suspected to occur concurrently in around a fifth of affected ears.

The presence, morphology and clinical significance of vertebral body malformations in an Australian population of French Bulldogs and Pugs

OBJECTIVE

To describe the incidence, morphology and clinical significance of congenital vertebral malformations (CVM) in two breeds of brachycephalic dogs presenting to a referral veterinary hospital.

DESIGN

Prospective cohort study series.

MATERIALS AND METHODS

Forty-nine French Bulldogs and Pugs were prospectively evaluated and placed in one of two groups based on whether or not they presented for neurological signs referable to spinal cord disease. A computed tomography (CT) of their entire spine was obtained and the presence and classification of CVM along with the degree of spinal kyphosis recorded for all dogs. Statistical analysis was performed to identify clinical associations between these factors (P < 0.05).

RESULTS

CVM were prevalent across both breeds with the French Bulldog having more malformations than the Pug (Kruskal-Wallis nonparametric analysis of variance, P < 0.0001). Breed associated vertebral malformation subtypes included butterfly subtype in French Bulldogs (Chi-square, P = 0.0002), and transitional subtype in Pugs (odds ratio, 22.7; P = 0.000). A new subtype, dorsal wedge, was observed in 12 cases. The presence, number and subtype of vertebral malformation were not reliable for predicting the development of neurological signs across both breeds (Chi-square, P > 0.05). However, spinal kyphosis >35° calculated via Cobb angle was associated with Pugs that had neurological deficits (Chi-square, P = 0.028).

CONCLUSIONS

Congenital vertebral malformations largely appear to be incidental findings in this population of French Bulldogs but are of more significance in the Pug breed when spinal kyphosis is >35°. French Bulldogs that have spinal cord disease and CVM are more likely to have pathology distant to CVM with intervertebral disc herniation most common.

Towards Forensic DNA Phenotyping for Predicting Visible Traits in Dogs

The popularity of dogs as human companions explains why these pets regularly come into focus in forensic cases such as bite attacks or accidents. Canine evidence, e.g., dog hairs, can also act as a link between the victim and suspect in a crime case due to the close contact between dogs and their owners. In line with human DNA identification, dog individualization from crime scene evidence is mainly based on the analysis of short tandem repeat (STR) markers. However, when the DNA profile does not match a reference, additional information regarding the appearance of the dog may provide substantial intelligence value. Key features of the dog's appearance, such as the body size and coat colour are well-recognizable and easy to describe even to non-dog experts, including most investigating officers and eyewitnesses. Therefore, it is reasonable to complement eyewitnesses' testimonies with externally visible traits predicted from associated canine DNA samples. Here, the feasibility and suitability of canine DNA phenotyping is explored from scratch in the form of a proof of concept study. To predict the overall appearance of an unknown dog from its DNA as accurately as possible, the following six traits were chosen: (1) coat colour, (2) coat pattern, (3) coat structure, (4) body size, (5) ear shape, and (6) tail length. A total of 21 genetic markers known for high predicting values for these traits were selected from previously published datasets, comprising 15 SNPs and six INDELS. Three of them belonged to SINE insertions. The experiments were designed in three phases. In the first two stages, the performance of the markers was tested on DNA samples from dogs with well-documented physical characteristics from different breeds. The final blind test, including dogs with initially withheld appearance information, showed that the majority of the selected markers allowed to develop composite sketches, providing a realistic impression of the tested dogs. We regard this study as the first attempt to evaluate the possibilities and limitations of forensic canine DNA phenotyping.

Multi-omics approach identifies germline regulatory variants associated with hematopoietic malignancies in retriever dog breeds

Histiocytic sarcoma is an aggressive hematopoietic malignancy of mature tissue histiocytes with a poorly understood etiology in humans. A histologically and clinically similar counterpart affects flat-coated retrievers (FCRs) at unusually high frequency, with 20% developing the lethal disease. The similar clinical presentation combined with the closed population structure of dogs, leading to high genetic homogeneity, makes dogs an excellent model for genetic studies of cancer susceptibility. To determine the genetic risk factors underlying histiocytic sarcoma in FCRs, we conducted multiple genome-wide association studies (GWASs), identifying two loci that confer significant risk on canine chromosomes (CFA) 5 (P_wald = 4.83x10^-9) and 19 (P_wald = 2.25x10^-7). We subsequently undertook a multi-omics approach that has been largely unexplored in the canine model to interrogate these regions, generating whole genome, transcriptome, and chromatin immunoprecipitation sequencing. These data highlight the PI3K pathway gene PIK3R6 on CFA5, and proximal candidate regulatory variants that are strongly associated with histiocytic sarcoma and predicted to impact transcription factor binding. The CFA5 association colocalizes with susceptibility loci for two hematopoietic malignancies, hemangiosarcoma and B-cell lymphoma, in the closely related golden retriever breed, revealing the risk contribution this single locus makes to multiple hematological cancers. By comparison, the CFA19 locus is unique to the FCR and harbors risk alleles associated with upregulation of TNFAIP6, which itself affects cell migration and metastasis. Together, these loci explain ~35% of disease risk, an exceptionally high value that demonstrates the advantages of domestic dogs for complex trait mapping and genetic studies of cancer susceptibility.

We have identified two regions of the canine genome that explain a striking 35% of risk for developing histiocytic sarcoma in FCRs. The disease is uniformly lethal, affects 20% of FCRs, and parallels a cancer of the same name in humans. Both regions harbor genes involved in cell migration and cancer-related pathways. The first includes variants in regulatory regions at the tumor suppressor PIK3R6 locus that are strongly associated with histiocytic sarcoma and likely confer risk for other hematopoietic cancers. FCRs with risk alleles at the second locus demonstrate increased expression of TNFAIP6, which correlates with poor prognosis in multiple human cancers. In identifying genomic differences between affected and unaffected dogs, we advance our understanding of both canine and human health biology and set the stage for the development of diagnostic and therapeutic strategies.

Genetics of canine subvalvular aortic stenosis (SAS)

Subvalvular aortic stenosis (SAS) is one of the most common congenital heart defects of dogs. The disease is characterized by obstruction of the left ventricular outflow tract, resulting in pressure overload on the left ventricle. The etiology of obstruction is a fibromuscular nodule, ridge, or ring of tissue that increases aortic outflow tract velocity. This review is focused on the prevalence, inheritance pattern, and current genetic insights of canine SAS. The prevalence of this disease was reported at 4.7 % in a large veterinary referral hospital. The mode of inheritance for this disease has also been described in breeds with a high disease prevalence such as the Bullmastiff, Bouvier des Flandres, Dogue de Bordeaux, Golden Retriever, Newfoundland, and Rottweiler. Genetic investigations seeking to identify causative mutations for SAS are lacking with only a single published variant associated with SAS in Newfoundlands.

Reliability and Validity of a Dog Personality and Unwanted Behavior Survey

Dogs have distinct, consistent personalities, but the structure of dog personality is still unclear. Dog personality and unwanted behavior are often studied with behavioral questionnaires. Even though many questionnaires are reliable and valid measures of behavior, all new questionnaire tools should be extensively validated. Here, we examined the structure of personality and six unwanted behavior questionnaire sections: noise sensitivity, fearfulness, aggression, fear of surfaces and heights, separation anxiety and impulsivity/inattention with factor analyses. Personality consisted of seven factors: Insecurity, Training focus, Energy, Aggressiveness/dominance, Human sociability, Dog sociability and Perseverance. Most unwanted behavior sections included only one factor, but the impulsivity/inattention section divided into two factors (Hyperactivity/impulsivity and Inattention) and the aggression section into four factors (Barking, Stranger directed aggression, Owner directed aggression and Dog directed aggression). We also examined the internal consistency, test-retest reliability, inter-rater reliability and convergent validity of the 17 personality and unwanted behavior traits and discovered excellent reliability and validity. Finally, we investigated the discriminant validity of the personality traits, which was good. Our findings indicate that this personality and unwanted behavior questionnaire is a reliable and valid tool that can be used to study personality and behavior extensively.

Sequencing of the Canine Cytochrome P450 CYP2C41 Gene and Genotyping of Its Polymorphic Occurrence in 36 Dog Breeds

Cytochrome P450 (CYP) drug metabolizing enzymes play an important role in efficient drug metabolism and elimination. Many CYPs are polymorphic and, thereby, drug metabolism can vary between individuals. In the case of canine CYP2C41, gene polymorphism was identified. However, as the first available canine genome sequences all were CYP2C41 negative, this polymorphism could not be clarified at the genomic level. The present study provides an exact characterization of the CYP2C41 gene deletion polymorphism at the genomic level and presents a PCR-based genotyping method that was used for CYP2C41 genotyping of 1,089 individual subjects from 36 different dog breeds. None of the Bearded Collie, Bernese Mountain, Boxer, Briard, French Bulldog or Irish Wolfhound subjects had the CYP2C41 gene in their genomes. In contrast, in the Chinese Char-Pei, Siberian Husky, Schapendoes and Kangal breeds, the CYP2C41 allele frequency was very high, with values of 67, 57, 43, and 34%, respectively. Interestingly, the site of gene deletion was identical for all CYP2C41 negative dogs, and all CYP2C41 positive dogs showed highly homologous sequence domains upstream and downstream from the CYP2C41 gene. CYP2C41 genotyping can now be routinely used in future pharmacokinetic studies in canines, in order to identify genetically-based poor or extensive drug metabolizers. This, together with more extensive in vitro drug screening for CYP2C41 substrates will help to determine the clinical relevance of CYP2C41, and to optimize drug treatment. Although the relative abundance of the CYP2C41 protein in the canine liver seems to not be very high, this CYP could substantially contribute to hepatic drug metabolism in dogs expressing CYP2C41 from both alleles and, when CYP2C41 shows higher catalytic activity to a given drug than other hepatic metabolic enzymes.

Long-term stress in dogs is related to the human-dog relationship and personality traits

Previously, we found that dogs belonging to the herding breed group, selected for human cooperation, synchronise their long-term stress levels with their owners. The aim of the current study was to investigate features that could influence long-term stress levels in ancient dog breeds, genetically closer to wolves, and dogs specifically selected to work independently of their owner. Twenty-four ancient breed dogs and 18 solitary hunting dogs were recruited and hair samples were obtained from both dogs and owners from which hair cortisol concentration (HCC) was analysed. Additionally, the owners completed lifestyle surveys, the Monash Dog Owner Relationship Scale (MDORS) on human–dog relationship, and both dog and owner personality questionnaires (Dog Personality questionnaire and Big Five Inventory survey). The results from the MDORS indicate that the subscale Perceived cost correlated to the dog HCC of tested breed groups: solitary hunting breeds (χ² = 4.95, P = 0.026, β = 0.055), ancient breeds (χ² = 2.74, P = 0.098, β = 0.027), and herding dogs included from a previous study (χ² = 6.82, P = 0.009, β = − 0.061). The HCC of the solitary hunting dogs was also related to the owner personality traits Agreeableness (χ² = 12.30, P < 0.001, β = − 0.060) and Openness (χ² = 9.56, P = 0.002, β = 0.048) suggesting a more substantial influence of the owner on the solitary hunting dog’s HCC compared to the ancient breeds. No effect of owner HCC on dog HCC was found in either ancient or in solitary hunting breeds. Hence, the long-term stress synchronisation is likely to be a trait in breeds selected for human cooperation. In conclusion, dog HCC is often related to the owners’ personality, but is primarily influenced by the owner-dog relationship.

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.

The impact of identity by descent on fitness and disease in dogs

Dogs and humans have coexisted together for thousands of years, but it was not until the Victorian Era that humans practiced selective breeding to produce the modern standards we see today. Strong artificial selection during the breed formation period has simplified the genetic architecture of complex traits and caused an enrichment of identity-by-descent (IBD) segments in the dog genome. This study demonstrates the value of IBD segments and utilizes them to infer the recent demography of canids, predict case-control status for complex traits, locate regions of the genome potentially linked to inbreeding depression, and to identify understudied breeds where there is potential to discover new disease-associated variants.

Domestic dogs have experienced population bottlenecks, recent inbreeding, and strong artificial selection. These processes have simplified the genetic architecture of complex traits, allowed deleterious variation to persist, and increased both identity-by-descent (IBD) segments and runs of homozygosity (ROH). As such, dogs provide an excellent model for examining how these evolutionary processes influence disease. We assembled a dataset containing 4,414 breed dogs, 327 village dogs, and 380 wolves genotyped at 117,288 markers and data for clinical and morphological phenotypes. Breed dogs have an enrichment of IBD and ROH, relative to both village dogs and wolves, and we use these patterns to show that breed dogs have experienced differing severities of bottlenecks in their recent past. We then found that ROH burden is associated with phenotypes in breed dogs, such as lymphoma. We next test the prediction that breeds with greater ROH have more disease alleles reported in the Online Mendelian Inheritance in Animals (OMIA). Surprisingly, the number of causal variants identified correlates with the popularity of that breed rather than the ROH or IBD burden, suggesting an ascertainment bias in OMIA. Lastly, we use the distribution of ROH across the genome to identify genes with depletions of ROH as potential hotspots for inbreeding depression and find multiple exons where ROH are never observed. Our results suggest that inbreeding has played a large role in shaping genetic and phenotypic variation in dogs and that future work on understudied breeds may reveal new disease-causing variation.

Analysis of whole-genome re-sequencing data of ducks reveals a diverse demographic history and extensive gene flow between Southeast/South Asian and Chinese populations

BACKGROUND

The most prolific duck genetic resource in the world is located in Southeast/South Asia but little is known about the domestication and complex histories of these duck populations.

RESULTS

Based on whole-genome resequencing data of 78 ducks (Anas platyrhynchos) and 31 published whole-genome duck sequences, we detected three geographic distinct genetic groups, including local Chinese, wild, and local Southeast/South Asian populations. We inferred the demographic history of these duck populations with different geographical distributions and found that the Chinese and Southeast/South Asian ducks shared similar demographic features. The Chinese domestic ducks experienced the strongest population bottleneck caused by domestication and the last glacial maximum (LGM) period, whereas the Chinese wild ducks experienced a relatively weak bottleneck caused by domestication only. Furthermore, the bottleneck was more severe in the local Southeast/South Asian populations than in the local Chinese populations, which resulted in a smaller effective population size for the former (7100-11,900). We show that extensive gene flow has occurred between the Southeast/South Asian and Chinese populations, and between the Southeast Asian and South Asian populations. Prolonged gene flow was detected between the Guangxi population from China and its neighboring Southeast/South Asian populations. In addition, based on multiple statistical approaches, we identified a genomic region that included three genes (PNPLA8, THAP5, and DNAJB9) on duck chromosome 1 with a high probability of gene flow between the Guangxi and Southeast/South Asian populations. Finally, we detected strong signatures of selection in genes that are involved in signaling pathways of the nervous system development (e.g., ADCYAP1R1 and PDC) and in genes that are associated with morphological traits such as cell growth (e.g., IGF1R).

CONCLUSIONS

Our findings provide valuable information for a better understanding of the domestication and demographic history of the duck, and of the gene flow between local duck populations from Southeast/South Asia and China.