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

Frequency of RPGRIP1 and MAP9 genetic modifiers of canine progressive retinal atrophy, in 132 breeds of dog

Variants in RPGRIP1 and MAP9, termed RPGRIP1ins44 and MAP9del respectively, are both associated with a form of canine progressive retinal atrophy referred to as RPGRIP1-CRD and have both been demonstrated to modify the development and progression of this disease. In the current study both variants were genotyped in at least 50 dogs of 132 diverse breeds and the data reveal that both segregate in multiple breeds. Individually, each variant is common within largely non-overlapping subsets of breed, and there is a negative correlation between their frequencies within breeds that segregate both variants. The frequency of both variants exceeds 0.05 in a single breed only, the Miniature Longhaired Dachshund. These data indicate that both variants are likely to be ancient and predate the development and genetic isolation of modern dog breeds. That both variants are present individually at high frequency in multiple breeds is consistent with the hypothesis that homozygosity of either variant alone is not associated with a clinically relevant phenotype, whereas the negative correlation between the two variants is consistent with the application of selective pressure, from dog breeders, against homozygosity at both loci, probably due to the more severe phenotype associated with homozygosity at both loci.

Exploring the impact of epidemiological and clinical factors on the progression of canine leishmaniosis by statistical and whole genome analyses: from breed predisposition to comorbidities

Canine leishmaniosis (CanL), caused by Leishmania infantum, is a complex disease of growing importance in Europe. Clinical manifestations result from the down-modulation of the host immune response through multiple host-parasite interactions. Although several factors might influence CanL progression, this is the first known study evaluating risk factors for its different clinical stages in a large referral hospital population (n = 35.669) from an endemic area, over a 20 year period. Genome-wide scans for selection signatures were also conducted to explore the genomic component of clinical susceptibility to L. infantum infection. The prevalence of CanL was 3.2% (16.7% stage I; 43.6% stage II; 32.1% stage III; 7.6% stage IV). Dog breed (crossbreed), bodyweight (<10 kg), living conditions (indoors), regular deworming treatment, and being vaccinated against Leishmania significantly decreased the transmission risk and the risk for developing severe clinical forms. Conversely, the detection of comorbidities was associated with advanced clinical forms, particularly chronic kidney disease, neoplasia, cryptorchidism, infectious tracheobronchitis and urate urolithiasis, although those did not impact the clinical outcome. Significant associations between an increased risk of severe clinical stages and findings in the anamnesis (renal or skin-related manifestations) and physical examination (ocular findings) were also detected, highlighting their diagnostic value in referred cases of CanL. Sixteen breeds were found to be significantly more susceptible to developing severe stages of leishmaniosis (e.g. Great Dane, Rottweiler, English Springer Spaniel, Boxer, American Staffordshire Terrier, Golden Retriever), while 20 breeds displayed a clinical resistantance phenotype and, thus, are more likely to mount an efficient immune response against L. infantum (e.g. Pointer, Samoyed, Spanish Mastiff, Spanish Greyhound, English Setter, Siberian Husky). Genomic analyses of these breeds retrieved 12 regions under selection, 63 candidate genes and pinpointed multiple biological pathways such as the IRE1 branch of the unfolded protein response, which could play a critical role in clinical susceptibility to L. infantum infection.

A Nonsynonymous Substitution of Lhx3 Leads to Changes in Body Size in Dogs and Mice

Lhx3 is a LIM-homeodomain transcription factor that affects body size in mammals by regulating the secretion of pituitary hormones. Akita, Shiba Inu, and Mame Shiba Inu dogs are Japanese native dog breeds that have different body sizes. To determine whether Lhx3 plays a role in the differing body sizes of these three dog breeds, we sequenced the Lhx3 gene in the three breeds, which led to the identification of an SNP in codon 280 (S280N) associated with body size. The allele frequency at this SNP differed significantly between the large Akita and the two kinds of smaller Shiba dogs. To validate the function of this SNP on body size, we introduced this change into the Lhx3 gene of mice. Homozygous mutant mice (S279N+/+) were found to have significantly increased body lengths and weights compared to heterozygous mutant (S279N+/-) and wild-type (S279N-/-) mice several weeks after weaning. These results demonstrate that a nonsynonymous substitution in Lhx3 plays an important role in regulating body size in mammals.

From "Husky" to "Bulldog"- behavioural correlates between castration and breed groups in the domestic dog (Canis lupus familiaris)

Neutering dogs is a widespread method and is carried out for various behavioural and husbandry reasons. This study's main objective is to investigate the behavioural correlations between neutering and the breed of male dogs. In order to possibly find breed-dependent differences in the behaviour of intact and castrated dogs, a differentiation between two clades - the "Huskies"(chow chow, shar pei, akita/shiba inu, alaskan malamute, siberian/alaskan husky) and the "Bulldogs" (german boxer, english/french bulldog, old english mastiff, boston terrier, english bull terrier, staffordshire bull terrier, american staffordshire terrier), based on Parker et al. [1], was made.Using an online questionnaire,, 31 neutered and 37 intact male dogs from the clade "Huskies" and 30 neutered and 38 intact male dogs from the clade "Bulldogs", participated in the study (N = 136).The survey included detailed questions on the dogs' personality and any associated issues as well as a behavioural anamnesis. Further questions relating to four of the "big five" personality dimensions based on the "Budapest questionnaire" by Turcsán et al. from 2011 [2] were also added.The results show, that neutered males from both breed clades more frequently displayed aggression toward humans than intact males (multinomial logistic regression, p = 0.002). When it came to aggression towards other dogs, it was the "Huskies" that differed significantly from the "Bulldogs"(multinomial logistic regression, p = 0.04) with being more aggressive. There were also significant differences in stress-related behaviour depending on castration status and breed (multinomial logistic regression, p < 0.001; Cramer's V = 0.33) and only the castration status had an impact on the significance (multinomial logistic regression, p < 0.001). The analysis also revealed significance for stress-indicating behaviour with dependence on neutering status (multinomial logistic regression, p < 0.001) and showed that stress as well as uncertainty are significantly more common in neutered dogs depending on breed and neutering status (multinomial logistic regression, p < 0.001; Cramer's V = 0.42), in that only neutered "Bulldogs" were stressed, but more "Huskies" overall.According to the Budapest questionnaire data, the "Bulldog" clade had considerably greater extraversion scores overall (ordinal regression, p < 0.001) than the "Huskies".Our findings highlight the risks and potential negative effects of neutering. Gonadectomy in no way substitutes for the dog receiving the necessary socialization, training, or bonding. Although in some circumstances it might have a favourable impact on the dog's behaviour, it should not be seen as a panacea for unwanted behaviour. Given that not all behaviours are influenced by sex hormones, every castration decision must be weighed up individually.

Urinary and household chemical exposures in pet dogs with urothelial cell carcinoma

Urothelial cell carcinoma (UCC) has been linked to environmental chemical exposures in people, but these risk factors are not well understood in dogs with UCC. We hypothesised that household chemical exposures contribute to the risk of UCC in pet dogs. This prospective cross-sectional case-control study included 37 dogs with UCC and 37 unaffected breed-, sex-, and age-matched controls. Dog owners completed an environmental questionnaire and household samples were collected and analysed for arsenic (in tap water and room dust) and acrolein (in room air). Urine samples from UCC dogs, control dogs, and consenting owners were analysed for inorganic arsenic species, the acrolein metabolite 3-HPMA, and the phenoxy herbicide 2,4-D. Public data on chlorination byproducts (total trihalomethanes) in municipal drinking water were also compared between case and control households. Dogs with UCC were more likely to swim in a pool (15.2%) compared with control dogs (0%) (OR 1.69, 95% CI = 1.69-∞; p = .02). Dogs with UCC also had more than 4-fold higher reported municipal water concentrations of chlorination byproducts (median 28.0 ppb) compared with controls (median 6.9 ppb; p < .0001). Dust arsenic concentrations were unexpectedly lower in case households (median 0.277 ng/cm2) compared with control households (median 0.401 ng/cm2; p = .0002). Other outcomes were not significantly different between groups. These data suggest that dog owners, especially those of breeds known to be at higher risk for UCC, consider limiting access to swimming pools and installing water filtration units that remove total trihalomethanes.

Genome-wide analyses reveals an association between invasive urothelial carcinoma in the Shetland sheepdog and NIPAL1

Naturally occurring canine invasive urinary carcinoma (iUC) closely resembles human muscle invasive bladder cancer in terms of histopathology, metastases, response to therapy, and low survival rate. The heterogeneous nature of the disease has led to the association of large numbers of risk loci in humans, however most are of small effect. There exists a need for new and accurate animal models of invasive bladder cancer. In dogs, distinct breeds show markedly different rates of iUC, thus presenting an opportunity to identify additional risk factors and overcome the locus heterogeneity encountered in human mapping studies. In the association study presented here, inclusive of 100 Shetland sheepdogs and 58 dogs of other breeds, we identify a homozygous protein altering point mutation within the NIPAL1 gene which increases risk by eight-fold (OR = 8.42, CI = 3.12-22.71), accounting for nearly 30% of iUC risk in the Shetland sheepdog. Inclusion of six additional loci accounts for most of the disease risk in the breed and explains nearly 75% of the phenotypes in this study. When combined with sequence data from tumors, we show that variation in the MAPK signaling pathway is an overarching cause of iUC susceptibility in dogs.

Structure and scaling of the middle ear in domestic dog breeds

Although domestic dogs vary considerably in both body size and skull morphology, behavioural audiograms have previously been found to be similar in breeds as distinct as a Chihuahua and a St Bernard. In this study, we created micro-CT reconstructions of the middle ears and bony labyrinths from the skulls of 17 dog breeds, including both Chihuahua and St Bernard, plus a mongrel and a wolf. From these reconstructions, we measured middle ear cavity and ossicular volumes, eardrum and stapes footplate areas and bony labyrinth volumes. All of these ear structures scaled with skull size with negative allometry and generally correlated better with condylobasal length than with maximum or interaural skull widths. Larger dogs have larger ear structures in absolute terms: the volume of the St Bernard's middle ear cavity was 14 times that of the Chihuahua. The middle and inner ears are otherwise very similar in morphology, the ossicular structure being particularly well-conserved across breeds. The expectation that larger ear structures in larger dogs would translate into hearing ranges shifted towards lower frequencies is not consistent with the existing audiogram data. Assuming that the audiograms accurately reflect the hearing of the breeds in question, oversimplifications in existing models of middle ear function or limitations imposed by other parts of the auditory system may be responsible for this paradox.

Transcript abundance of hepatic drug-metabolizing enzymes in two dog breeds compared with 14 species including humans

Drug-metabolizing enzymes are important in drug development and therapy, but have not been fully identified and characterized in many species, lines, and breeds. Liver transcriptomic data were analyzed for phase I cytochromes P450, flavin-containing monooxygenases, and carboxylesterases and phase II UDP-glucuronosyltransferases, sulfotransferases, and glutathione S-transferases. Comparisons with a variety of species (humans, rhesus macaques, African green monkeys, baboons, common marmosets, cattle, sheep, pigs, cats, dogs, rabbits, tree shrews, rats, mice, and chickens) revealed both general similarities and differences in the transcript abundances of drug-metabolizing enzymes. Similarly, Beagle and Shiba dogs were examined by next-generation sequencing (RNA-seq). Consequently, no substantial differences in transcript abundance were noted in different breeds of pigs and dogs and in different lines of mice and rats. Therefore, the expression profiles of hepatic drug-metabolizing enzyme transcripts appear to be similar in Shiba and Beagle dogs and pig breeds and the rat and mouse lines analyzed, although some differences were found in other species.

Transposons in the Williams-Beuren Syndrome Critical Region are Associated with Social Behavior in Assistance Dogs

A strong signature of selection in the domestic dog genome is found in a five-megabase region of chromosome six in which four structural variants derived from transposons have previously been associated with human-oriented social behavior, such as attentional bias to social stimuli and social interest in strangers. To explore these genetic associations in more phenotypic detail-as well as their role in training success in a specialized assistance dog program-we genotyped 1001 assistance dogs from Canine Companions for Independence®, including both successful graduates and dogs released from the training program for behaviors incompatible with their working role. We collected phenotypes on each dog using puppy-raiser questionnaires, trainer questionnaires, and both cognitive and behavioral tests. Using Bayesian mixed models, we found strong associations (95% credibility intervals excluding zero) between genotypes and certain behavioral measures, including separation-related problems, aggression when challenged or corrected, and reactivity to other dogs. Furthermore, we found moderate differences in the genotypes of dogs who graduated versus those who did not; insertions in GTF2I showed the strongest association with training success (β = 0.23, CI95% = - 0.04, 0.49), translating to an odds-ratio of 1.25 for one insertion. Our results provide insight into the role of each of these four transposons in canine sociability and may inform breeding and training practices for working dog organizations. Furthermore, the observed importance of the gene GTF2I supports the emerging consensus that variation in GTF2I genotypes and expression have important consequences for social behavior broadly.

PAPipe: A Pipeline for Comprehensive Population Genetic Analysis

Advancements in next-generation sequencing (NGS) technologies have led to a substantial increase in the availability of population genetic variant data, thus prompting the development of various population analysis tools to enhance our understanding of population structure and evolution. The tools that are currently used to analyze population genetic variant data generally require different environments, parameters, and formats of the input data, which can act as a barrier preventing the wide-spread usage of such tools by general researchers who may not be familiar with bioinformatics. To address this problem, we have developed an automated and comprehensive pipeline called PAPipe to perform nine widely used population genetic analyses using population NGS data. PAPipe seamlessly interconnects and serializes multiple steps, such as read trimming and mapping, genetic variant calling, data filtering, and format converting, along with nine population genetic analyses such as principal component analysis, phylogenetic analysis, population tree analysis, population structure analysis, linkage disequilibrium decay analysis, selective sweep analysis, population admixture analysis, sequentially Markovian coalescent analysis, and fixation index analysis. PAPipe also provides an easy-to-use web interface that allows for the parameters to be set and the analysis results to be browsed in intuitive manner. PAPipe can be used to generate extensive results that provide insights that can help enhance user convenience and data usability. PAPipe is freely available at https://github.com/jkimlab/PAPipe.

Frameshift Variant in AMPD2 in Cirneco dell'Etna Dogs with Retinopathy and Tremors

While the manifestations of many inherited retinal disorders are limited to loss of vision, others are part of a syndrome that affects multiple tissues, particularly the nervous system. Most syndromic retinal disorders are thought to be recessively inherited. Two dogs out of a litter of Cirneco dell' Etna dogs, both males, showed signs of retinal degeneration, along with tremors and signs described as either atypical seizures or paroxysmal dyskinesias, while the other two male littermates were normal. We named this oculo-neurological syndrome CONS (Cirneco oculo-neurological syndrome), and undertook homozygosity mapping and whole-genome sequencing to determine its potential genetic etiology. Notably, we detected a 1-bp deletion in chromosome 6 that was predicted to cause a frameshift and premature stop codon within the canine AMPD2 gene, which encodes adenosine monophosphate deaminase, an enzyme that converts adenosine 5'-monophosphate (AMP) to inosine 5'-monophosphate (IMP). Genotyping of the available Cirneco population suggested perfect segregation between cases and controls for the variant. Moreover, this variant was absent in canine genomic databases comprised of thousands of unaffected dogs. The AMPD2 genetic variant we identified in dogs presents with retinal manifestations, adding to the spectrum of neurological manifestations associated with AMPD2 variants in humans.

A universal panel of STR loci for the study of polymorphism of the species Canis lupus and forensic identification of dog and wolf

Commercial panels of microsatellite (STR) loci are intended for DNA analysis of the domestic dog (Canis lupus familiaris) and, therefore, when genotyping the Grey wolf (Canis lupus lupus), most markers reveal significant deviations from the Hardy-Weinberg equilibrium and have a low informative value, which complicates their use in a forensic examination. The aim of this study was to select STR markers that equally effectively reflect population polymorphism in the wolf and the dog, and to create a universal panel for the identification of individuals in forensic science. Based on the study of polymorphisms of 34 STR loci, a CPlex panel of 15 autosomal loci and two sex loci was developed, which is equally suitable for identifying wolfs and dogs. Analysis of molecular variance (AMOVA) between samples revealed significant differentiation values (FST = 0.0828, p < 0.05), which allows the panel to be used for differentiating between wolf and dog samples. For the first time in the forensic examination of objects of animal origin in the Republic of Belarus, population subdivision coefficients (θ-values) were calculated for each of the 15 STR loci of the test system being reported. It was shown that the values of the genotype frequency, when averaged over all studied animals without and with considering the θ-value, differ by three orders of magnitude (3.39 · 10-17 and 4.71 · 10-14, respectively). The use of population subdivision coefficients will provide the researcher with the most relevant results of an expert identification study. The test system was validated in accordance with the protocol of the Scientific Working Group on DNA Analysis Methods. A computational tool was developed to automate the analysis of genetic data on the wolf and dog in the forensic examination; two guides were approved for practicing forensic experts. This methodology is being successfully used in expert practice in investigating cases of illegal hunting, animal abuse and other offenses in the Republic of Belarus.

Longevity of companion dog breeds: those at risk from early death

The companion dog is one of the most phenotypically diverse species. Variability between breeds extends not only to morphology and aspects of behaviour, but also to longevity. Despite this fact, little research has been devoted to assessing variation in life expectancy between breeds or evaluating the potential for phylogenetic characterisation of longevity. Using a dataset of 584,734 unique dogs located within the UK, including 284,734 deceased, we present variation in longevity estimates within the following: parental lineage (purebred = 1 breed, crossbred ≥ 2 breeds), breed (n = 155), body size (large, medium, small), sex (male, female) and cephalic index (brachycephalic, mesocephalic, dolichocephalic). Survival estimates were then partitioned amongst phylogenetic clades: providing evidence that canine evolutionary history (via domestication and associated artificial selection) is associated with breed lifespan. This information provides evidence to inform discussions regarding pedigree health, whilst helping current/prospective owners, breeders, policy makers, funding bodies and welfare organisations improve decision making regarding canine welfare.

Pharmacogenomics of poor drug metabolism in greyhounds: Canine P450 oxidoreductase genetic variation, breed heterogeneity, and functional characterization

Greyhounds metabolize cytochrome P450 (CYP) 2B11 substrates more slowly than other dog breeds. However, CYP2B11 gene variants associated with decreased CYP2B11 expression do not fully explain reduced CYP2B11 activity in this breed. P450 oxidoreductase (POR) is an essential redox partner for all CYPs. POR protein variants can enhance or repress CYP enzyme function in a CYP isoform and substrate dependent manner. The study objectives were to identify POR protein variants in greyhounds and determine their effect on coexpressed CYP2B11 and CYP2D15 enzyme function. Gene sequencing identified two missense variants (Glu315Gln and Asp570Glu) forming four alleles, POR-H1 (reference), POR-H2 (570Glu), POR-H3 (315Gln, 570Glu) and POR-H4 (315Gln). Out of 68 dog breeds surveyed, POR-H2 was widely distributed across multiple breeds, while POR-H3 was largely restricted to greyhounds and Scottish deerhounds (35% allele frequencies), and POR-H4 was rare. Three-dimensional protein structure modelling indicated significant effects of Glu315Gln (but not Asp570Glu) on protein flexibility through loss of a salt bridge between Glu315 and Arg519. Recombinant POR-H1 (reference) and each POR variant (H2-H4) were expressed alone or with CYP2B11 or CYP2D15 in insect cells. No substantial effects on POR protein expression or enzyme activity (cytochrome c reduction) were observed for any POR variant (versus POR-H1) when expressed alone or with CYP2B11 or CYP2D15. Furthermore, there were no effects on CYP2B11 or CYP2D15 protein expression, or on CYP2D15 enzyme kinetics by any POR variant (versus POR-H1). However, Vmax values for 7-benzyloxyresorufin, propofol and bupropion oxidation by CYP2B11 were significantly reduced by coexpression with POR-H3 (by 34-37%) and POR-H4 (by 65-72%) compared with POR-H1. Km values were unaffected. Our results indicate that the Glu315Gln mutation (common to POR-H3 and POR-H4) reduces CYP2B11 enzyme function without affecting at least one other major canine hepatic P450 (CYP2D15). Additional in vivo studies are warranted to confirm these findings.

Allele and genotype frequencies for primary hereditary cataract, multifocal retinopathy 1, and degenerative myelopathy in Pyrenean Mountain dog from Italy

Pyrenean Mountain Dog (PMD) is an ancient dog breed firstly described in XIV century in the Pyrenees Region and nowadays diffused both in Europe and in the US. Hereditary Cataract (HC), defined as the inherited opacity of the lens, involves clinical signs ranging from reduced vision to glaucoma. A molecular basis of HC was firstly described in Staffordshire Bull Terriers and then reported in multiple canine breeds. The HC-associated variation is a single nucleotide deletion in HSF4 gene that introduces a premature stop codon (c.962del, p.Ala321*). Multifocal Retinopathy 1 (MR) is an ocular disorder characterized by multiple areas of retinal degeneration, caused in various dog breeds (including PMD) by a single nucleotide variant (SNV) in BEST1 gene that generates a premature stop codon (c.73G>A, p.Arg25*). Degenerative Myelopathy (DM) is an adult-onset, progressive neurodegenerative disease and it is associated to a SNV in SOD1 gene causing a change in aminoacidic sequence of the protein (c.118G>A, p.Glu40Lys). This causative variant has been described in various dog breeds, including PMD. Aim of this study was to determine the allele frequencies for the abovementioned three genetic diseases in the Italian breeding PMD population. The survey found no dogs carrying the allele (deletion) associated with HC, while three dogs (6 %) were heterozygous (G/A) for the MR-associated variant, and seven dogs (13 %) were heterozygous (G/A) for the DM-associated alteration, indicating that the variant alleles frequency were 0  %, 3 %, and 7 %, respectively. Appropriate mating management is suggested for the prevention of genetic diseases spreading in the PMD population.

The effect of body size and inbreeding on cancer mortality in breeds of the domestic dog: a test of the multi-stage model of carcinogenesis

Cancer is a leading cause of death in domestic dogs. Deaths due to cancer vary widely among breeds, providing an opportunity for testing the multi-stage model of carcinogenesis. This model underpins evolutionary and basic studies of cancer suppression and predicts a linear increase in cancer with breed size, an expectation complicated by bigger breeds having a shorter lifespan (decreasing risk). Using three independent datasets, the weight and lifespan of breeds provided a good fit of lifetime cancer mortality to the multi-stage model, the fit suggesting many canine cancers are initiated by four driver mutations. Of 85 breeds in more than one dataset, only flat-coated retriever showed significantly elevated cancer mortality, with Scottish terrier, Bernese mountain dog and bullmastiff also showing notable risk (greater than 50% over expected). Analysis of breed clades suggested terriers experience elevated cancer mortality. There was no evidence that the lower mass-specific metabolic rate of larger breeds reduced cancer risk. Residuals indicated increased breed inbreeding shortened expected lifespan, but had no overall effect on cancer mortality. The results provide a baseline for identifying increased breed risk for specific cancers and demonstrate that, unless selection promotes increased cancer suppression, the evolution of larger longer-lived animals leads to a predictable increased cancer risk.

Characterization of CD3+/CD20+ canine large-cell lymphoma

Immunophenotyping of canine large-cell lymphoma (LCL) for B-cell and T-cell surface antigens is commonly performed to better predict the clinical outcome. Expression of surface antigen CD3 is associated with T-cell malignancies; surface antigen CD20 is expressed on B cells. However, a small subset of canine LCLs expresses both CD3 and CD20 (CD3+/CD20+); this form of lymphoma remains poorly defined at the molecular level. In a retrospective study, we aimed to better characterize immunophenotypic properties and antigen receptor clonality of CD3+/CD20+ LCL. We selected formalin-fixed, paraffin-embedded tissues from 10 cases of CD3+/CD20+ LCL and breed-matched controls of peripheral large T-cell lymphoma (PTCL) and diffuse large B-cell lymphoma (DLBCL). Using PCR for antigen receptor rearrangement (PARR), we identified monoclonal T-cell receptor gamma (TCRγ) rearrangements in all CD3+/CD20+ cases. Three of 10 cases had monoclonal rearrangements in the immunoglobulin heavy chain (IgH), supportive of cross-lineage rearrangement. There was no significant difference in the frequency of antigen receptor rearrangement between CD3+/CD20+ and PTCL cases. In comparison with DLBCL, CD3+/CD20+ LCL had TCRγ rearrangement more frequently and IgH rearrangement less frequently, respectively. Immunolabeling of the B-cell marker PAX5 occurred less frequently in all CD3+/CD20+ LCL cases compared to the DLBCL controls. Immunolabeling for BCL-2 was robust, regardless of immunophenotype. Nuclear Ki67 positivity was variable in CD3+/CD20+ cases, indicating a heterogeneity in proliferation. Overall, cases of canine CD3+/CD20+ LCL had properties similar to PTCL, suggesting a similar histogenesis of these 2 subsets.

Global Frequency Analyses of Canine Progressive Rod-Cone Degeneration-Progressive Retinal Atrophy and Collie Eye Anomaly Using Commercial Genetic Testing Data

Hundreds of genetic variants associated with canine traits and disorders have been identified, with commercial tests offered. However, the geographic distributions and changes in allele and genotype frequencies over prolonged, continuous periods of time are lacking. This study utilized a large set of genotypes from dogs tested for the progressive rod-cone degeneration-progressive retinal atrophy (prcd-PRA) G>A missense PRCD variant (n = 86,667) and the collie eye anomaly (CEA)-associated NHEJ1 deletion (n = 33,834) provided by the commercial genetic testing company (Optigen/Wisdom Panel, Mars Petcare Science & Diagnostics). These data were analyzed using the chi-square goodness-of-fit test, time-trend graphical analysis, and regression modeling in order to evaluate how test results changed over time. The results span fifteen years, representing 82 countries and 67 breeds/breed mixes. Both diseases exhibited significant differences in genotype frequencies (p = 2.7 × 10-152 for prcd-PRA and 0.023 for CEA) with opposing graphical trends. Regression modeling showed time progression to significantly affect the odds of a dog being homozygous or heterozygous for either disease, as do variables including breed and breed popularity. This study shows that genetic testing informed breeding decisions to produce fewer affected dogs. However, the presence of dogs homozygous for the disease variant, especially for prcd-PRA, was still observed fourteen years after test availability, potentially due to crosses of unknown carriers. This suggests that genetic testing of dog populations should continue.

Chemotherapy for the treatment of intracranial glioma in dogs

Gliomas are the second most common primary brain tumor in dogs and although they are associated with a poor prognosis, limited data are available relating to the efficacy of standard therapeutic options such as surgery, radiation and chemotherapy. Additionally, canine glioma is gaining relevance as a naturally occurring animal model that recapitulates human disease with fidelity. There is an intense comparative research drive to test new therapeutic approaches in dogs and assess if results translate efficiently into human clinical trials to improve the poor outcomes associated with the current standard-of-care. However, the paucity of data and controversy around most appropriate treatment for intracranial gliomas in dogs make comparisons among modalities troublesome. To further inform therapeutic decision-making, client discussion, and future studies evaluating treatment responses, the outcomes of 127 dogs with intracranial glioma, either presumed (n = 49) or histologically confirmed (n = 78), that received chemotherapy as leading or adjuvant treatment are reviewed here. This review highlights the status of current chemotherapeutic approaches to intracranial gliomas in dogs, most notably temozolomide and lomustine; areas of novel treatment currently in development, and difficulties to consensuate and compare different study observations. Finally, suggestions are made to facilitate evidence-based research in the field of canine glioma therapeutics.

CanISO: a database of genomic and transcriptomic variations in domestic dog (Canis lupus familiaris)

BACKGROUND

The domestic dog, Canis lupus familiaris, is a companion animal for humans as well as an animal model in cancer research due to similar spontaneous occurrence of cancers as humans. Despite the social and biological importance of dogs, the catalogue of genomic variations and transcripts for dogs is relatively incomplete.

RESULTS

We developed CanISO, a new database to hold a large collection of transcriptome profiles and genomic variations for domestic dogs. CanISO provides 87,692 novel transcript isoforms and 60,992 known isoforms from whole transcriptome sequencing of canine tumors (N = 157) and their matched normal tissues (N = 64). CanISO also provides genomic variation information for 210,444 unique germline single nucleotide polymorphisms (SNPs) from the whole exome sequencing of 183 dogs, with a query system that searches gene- and transcript-level information as well as covered SNPs. Transcriptome profiles can be compared with corresponding human transcript isoforms at a tissue level, or between sample groups to identify tumor-specific gene expression and alternative splicing patterns.

CONCLUSIONS

CanISO is expected to increase understanding of the dog genome and transcriptome, as well as its functional associations with humans, such as shared/distinct mechanisms of cancer. CanISO is publicly available at https://www.kobic.re.kr/caniso/ .

Lymphoma in Border Collies: Genome-Wide Association and Pedigree Analysis

There has been considerable interest in studying cancer in dogs and its potential as a model system for humans. One area of research has been the search for genetic risk variants in canine lymphoma, which is amongst the most common canine cancers. Previous studies have focused on a limited number of breeds, but none have included Border Collies. The aims of this study were to identify relationships between Border Collie lymphoma cases through an extensive pedigree investigation and to utilise relationship information to conduct genome-wide association study (GWAS) analyses to identify risk regions associated with lymphoma. The expanded pedigree analysis included 83,000 Border Collies, with 71 identified lymphoma cases. The analysis identified affected close relatives, and a common ancestor was identified for 54 cases. For the genomic study, a GWAS was designed to incorporate lymphoma cases, putative "carriers", and controls. A case-control GWAS was also conducted as a comparison. Both analyses showed significant SNPs in regions on chromosomes 18 and 27. Putative top candidate genes from these regions included DLA-79, WNT10B, LMBR1L, KMT2D, and CCNT1.

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.

Multiple Origins and Genomic Basis of Complex Traits in Sighthounds

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

The Shepherd and the Hunter: A Genomic Comparison of Italian Dog Breeds

Shepherd and hunting dogs have undergone divergent selection for specific tasks, resulting in distinct phenotypic and behavioural differences. Italy is home to numerous recognized and unrecognized breeds of both types, providing an opportunity to compare them genomically. In this study, we analysed SNP data obtained from the CanineHD BeadChip, encompassing 116 hunting dogs (representing 6 breeds) and 158 shepherd dogs (representing 9 breeds). We explored the population structure, genomic background, and phylogenetic relationships among the breeds. To compare the two groups, we employed three complementary methods for selection signature detection: FST, XP-EHH, and ROH. Our results reveal a clear differentiation between shepherd and hunting dogs as well as between gun dogs vs. hounds and guardian vs. herding shepherd dogs. The genomic regions distinguishing these groups harbour several genes associated with domestication and behavioural traits, including gregariousness (WBSRC17) and aggressiveness (CDH12 and HTT). Additionally, genes related to morphology, such as size and coat colour (ASIP and TYRP1) and texture (RSPO2), were identified. This comparative genomic analysis sheds light on the genetic underpinnings of the phenotypic and behavioural variations observed in Italian hunting and shepherd dogs.

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

INTRODUCTION

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Pain sensitivity differs between dog breeds but not in the way veterinarians believe

Background

Veterinarians hold distinct breed-specific pain sensitivity beliefs that differ from the general public but are highly consistent with one another. This is remarkable as there is no current scientific evidence for biological differences in pain sensitivity across dog breeds. Therefore, the present study evaluated whether pain sensitivity thresholds differ across a set of dog breeds and, if so, whether veterinarians' pain sensitivity ratings explain these differences or whether these ratings are attributed to behavioral characteristics.

Methods

Pain sensitivity thresholds [using quantitative sensory testing (QST) methods] and canine behaviors (using owner questionnaires and emotional reactivity tests) were prospectively measured across selected dog breeds. Adult, healthy dogs from 10 dog breeds/breed types were recruited, representing breeds subjectively rated by veterinarians as high (chihuahua, German shepherd, Maltese, Siberian husky), average (border collie, Boston terrier, Jack Russell terrier), or low (golden retriever, pitbull, Labrador retriever) pain sensitivity. A final sample of 149 dogs was included in statistical analyses.

Results

Veterinarians' pain sensitivity ratings provided a minimal explanation for pain sensitivity thresholds measured using QST in dogs; however, dog breeds did differ in their pain sensitivity thresholds across the QST methods evaluated. Breed differences were observed for some aspects of emotional reactivity tests; however, these behavioral differences did not explain the differences in pain sensitivity thresholds found. Veterinarians' pain sensitivity ratings were positively associated with dog approach scores for the disgruntled stranger test suggesting that the way dogs greet strangers may be a factor influencing veterinarians' ratings of pain sensitivity across dog breeds.

Conclusions and clinical relevance

Overall, these findings highlight a need to investigate biological mechanisms that may explain breed differences in pain sensitivity because this may inform pain management recommendations. Further, future research should focus on when and how these breed-specific pain sensitivity beliefs developed in veterinarians, as veterinarians' beliefs could impact the recognition and treatment of pain for canine patients.

A 4-lineage Statistical Suite to Evaluate the Support of Large-Scale Retrotransposon Insertion Data to Reconstruct Evolutionary Trees

Retrophylogenomics makes use of genome-wide retrotransposon presence/absence insertion patterns to resolve questions in phylogeny and population genetics. In the genomics era, evaluating high-throughput data requires the associated development of appropriately powerful statistical tools. The currently used KKSC 3-lineage statistical test for estimating the significance of retrophylogenomic data is limited by the number of possible tree topologies it can assess in one step. To improve on this, we have extended the analysis to simultaneously compare four lineages, enabling us to evaluate ten distinct presence/absence insertion patterns for 26 possible tree topologies plus 129 trees with different incidences of hybridization or introgression. The new tool provides statistics for cases involving multiple ancestral hybridizations/introgressions, ancestral incomplete lineage sorting, bifurcation, and polytomy. The test is embedded in a user-friendly web R application (http://retrogenomics.uni-muenster.de:3838/hammlet/) and is available for use by the scientific community. [ancestral hybridization/introgression; ancestral incomplete lineage sorting (ILS); empirical distribution; KKSC-statistics; 4-lineage (4-LIN) insertion polymorphism; polytomy; retrophylogenomics.].

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

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

Would You Detour with Me? Association between Functional Breed Selection and Social Learning in Dogs Sheds Light on Elements of Dog-Human Cooperation

Interspecific social learning is a main synchronizing mechanism that enables dogs to adapt to the anthropogenic niche. It is not known whether dogs in general possess the capacity of learning from humans or whether more recent selective events have affected their ability to learn from humans. We hypothesized that cooperative and independent working dog breeds may behave differently in a social learning task. Dogs (N = 78 from 16 cooperative and 18 independent breeds) had to detour a transparent, V-shaped wire mesh fence. The experiment consisted of three one-minute-long trials. The control condition did not include a demonstration. In the demonstration condition, the experimenter placed a reward in the inside corner by walking around the fence. Cooperative dogs reached the target significantly faster, while independent dogs did not detour faster in trials 2 and 3 after the human demonstration. Detour latencies were not associated with the keeping conditions and training level of the subjects. As we assembled both test groups from several genetically distantly related breeds, we can exclude the explanation that higher cooperativity emerged only in particular clades of dogs; instead, functional selection for particular working tasks could enhance capacities that affect a wide range of socio-cognitive traits in dogs.

Systematically identifying genetic signatures including novel SNP-clusters, nonsense variants, frame-shift INDELs, and long STR expansions that potentially link to unknown phenotypes existing in dog breeds

BACKGROUND

In light of previous studies that profiled breed-specific traits or used genome-wide association studies to refine loci associated with characteristic morphological features in dogs, the field has gained tremendous genetic insights for known dog traits observed among breeds. Here we aim to address the question from a reserve perspective: whether there are breed-specific genotypes that may underlie currently unknown phenotypes. This study provides a complete set of breed-specific genetic signatures (BSGS). Several novel BSGS with significant protein-altering effects were highlighted and validated.

RESULTS

Using the next generation whole-genome sequencing technology coupled with unsupervised machine learning for pattern recognitions, we constructed and analyzed a high-resolution sequence map for 76 breeds of 412 dogs. Genomic structures including novel single nucleotide polymorphisms (SNPs), SNP clusters, insertions, deletions (INDELs) and short tandem repeats (STRs) were uncovered mutually exclusively among breeds. We also partially validated some novel nonsense variants by Sanger sequencing with additional dogs. Four novel nonsense BSGS were found in the Bernese Mountain Dog, Samoyed, Bull Terrier, and Basset Hound, respectively. Four INDELs resulting in either frame-shift or codon disruptions were found in the Norwich Terrier, Airedale Terrier, Chow Chow and Bernese Mountain Dog, respectively. A total of 15 genomic regions containing three types of BSGS (SNP-clusters, INDELs and STRs) were identified in the Akita, Alaskan Malamute, Chow Chow, Field Spaniel, Keeshond, Shetland Sheepdog and Sussex Spaniel, in which Keeshond and Sussex Spaniel each carried one amino-acid changing BSGS in such regions.

CONCLUSION

Given the strong relationship between human and dog breed-specific traits, this study might be of considerable interest to researchers and all. Novel genetic signatures that can differentiate dog breeds were uncovered. Several functional genetic signatures might indicate potentially breed-specific unknown phenotypic traits or disease predispositions. These results open the door for further investigations. Importantly, the computational tools we developed can be applied to any dog breeds as well as other species. This study will stimulate new thinking, as the results of breed-specific genetic signatures may offer an overarching relevance of the animal models to human health and disease.

Psychological Stress Is Associated with Increased Cancer Risk in Dogs

Although there is evidence that psychological stress may be associated with increased cancer risk, the effect of stress on cancer risk is difficult to study, both in humans, due to socioeconomic factors, and in animal models, due to questionable biological relevance. Here, we test whether heritable canine temperament that increases psychological stress is associated with cancer risk. The study data are breed-specific averages of incidences of multiple cancer types and of temperament classes. The latter are derived from a latent class analysis of behavioral questionnaires completed by owners (C-BARQ). We thus classified the dogs according to whether they are calm vs. reactive within and across breeds. Using meta-analysis approaches, we modeled the risk of multiple cancer types in calm vs. reactive dogs. We adjusted for breed averages of body mass and lifespan, which are common confounders that impact cancer. Our study confirms that body size has a significant effect of on risk of multiple types of cancers in dogs and shows for the first time that temperament also has a moderate effect. These findings suggest dog models of heritable psychological stress are suitable for molecular epidemiological and translational studies on its effects on cancer risk.

Breed, age, and social environment are associated with personality traits in dogs

Personality in pets and other domesticated animals is important for their well-being and it can also influence human-animal relationships. Genetic and environmental factors influencing unwanted behavior in dogs are somewhat well known, but the factors influencing dog personality remain understudied. Here we examined environmental and demographic factors associated with seven broad personality traits in a survey of over 11,000 dogs. We utilized linear models and extensive model validation to examine the factors that have the most significant influences on personality and calculated effect sizes to assess the importance of these variables. Breed and age had the strongest associations with dog personality traits. Some environmental factors, especially puppyhood socialization, were also associated with personality. All factors had small effect sizes, highlighting that a lot of variation in personality remains unexplained. Our results indicate that personality traits are complex and strikingly similar in dogs, humans, and other nonhuman animals.

Population genomics of American mink using genotype data

Understanding the genetic structure of the target population is critically important to develop an efficient genomic selection program in domestic animals. In this study, 2,973 American mink of six color types from two farms (Canadian Centre for Fur Animal Research (CCFAR), Truro, NS and Millbank Fur Farm (MFF), Rockwood, ON) were genotyped with the Affymetrix Mink 70K panel to compute their linkage disequilibrium (LD) patterns, effective population size (Ne), genetic diversity, genetic distances, and population differentiation and structure. The LD pattern represented by average r ², decreased to <0.2 when the inter-marker interval reached larger than 350 kb and 650 kb for CCFAR and MFF, respectively, and suggested at least 7,700 and 4,200 single nucleotide polymorphisms (SNPs) be used to obtain adequate accuracy for genomic selection programs in CCFAR and MFF respectively. The Ne for five generations ago was estimated to be 76 and 91 respectively. Our results from genetic distance and diversity analyses showed that American mink of the various color types had a close genetic relationship and low genetic diversity, with most of the genetic variation occurring within rather than between color types. Three ancestral genetic groups was considered the most appropriate number to delineate the genetic structure of these populations. Black (in both CCFAR and MFF) and pastel color types had their own ancestral clusters, while demi, mahogany, and stardust color types were admixed with the three ancestral genetic groups. This study provided essential information to utilize the first medium-density SNP panel for American mink in their genomic studies.

The Histopathological Characteristic of Gastric Carcinoma in the Belgian Tervueren and Groenendael Dog: A Comparison of Two Classification Methods

Gastric carcinoma is generally considered to be a rare disease in dogs, carrying a grave prognosis. However, in the Tervueren and Groenendael varieties of the Belgian Shepherd dog breed, the disease is highly prevalent. While histopathology is the gold standard for diagnosing gastric carcinoma, there is no general consensus on the methods for histological classification in these cases. Biopsies of a group of 61 dogs with confirmed gastric carcinoma (45 Tervueren and 16 Groenendael) were examined and classified according to World Health Organization (WHO) and Laurén classifications. Kaplan-Meier curves were used to compare survival between the different subtypes and simple and multiple linear regression were used to analyse the association between age of onset and breed variant, sex, neuter status, location of the tumour, inflammation score, and Laurén and WHO classifications. Mean age at diagnosis was significantly different in Groenendael (10.1 ± 2.01) and Tervueren dogs (8.5 ± 1.90). The Laurén classification resulted in 29 (48%) diffuse- and 32 (52%) intestinal-type tumours. Applying the WHO classification resulted in 30 (49%) tubular carcinoma growth patterns and 31 (51%) others. Median survival time was significantly reduced for the diffuse type as compared to the intestinal type according to the Laurén classification, with the same median survival time results for tubular compared to non-tubular subtypes according to the WHO classification (median survival time of 61 vs. 182 days, respectively). Using the WHO and Lauren classification on tumour biopsies may help the practising clinician in the prognostication of gastric carcinoma in Tervueren and Groenendael dogs.

A Descriptive Methodology for Studying the Ontogeny of Object Play and Breed Differences in Dogs (Canis lupus familiaris)

Play behavior is a prominent aspect of juvenile behavior for many animals, yet early development, especially play with objects, has received little attention. Our previous study on object play introduced our general methods, focusing on litter differences in the developmental trajectory of object play and toy preferences. Here, we present a detailed ethogram of more than 30 observed object play behaviors. We focus on breed differences in the development of play in the three following breeds: Welsh Terriers, Vizslas, and standard Poodles. Puppies were video recorded from 3 to 7 weeks of age at half-week intervals upon the introduction of a standard set of five toys into their home environments. Ten minutes of video from each session for each puppy were analyzed using the Noldus Observer XT program. Aside from analyzing individual behaviors, they were also grouped into three behavioral categories. These were behaviors that occurred only in a solitary context, only in a social context, or in both contexts. Solitary object play developed first, and social object play developed later across breeds. There was a significant three-way interaction between breed, developmental age, and the context in which play occurred. Pairwise comparisons within each breed, age, and context are discussed, but a prominent result is that the onset of many behaviors occurred later in Welsh Terriers compared to the other breeds.

Identification of Genomic Signatures in Bullmastiff Dogs Using Composite Selection Signals Analysis of 23 Purebred Clades

Dog breeds represent canine sub-populations with distinctive phenotypic features and limited genetic diversity. We have established a resource to study breed-specific genetic diversity. Utilising genetic resources within our laboratory biobank, public domain genotype data and the phylogenetic framework of 23 breed clades, the primary objective for this study was to identify genomic regions that differentiate the Bullmastiff breed. Through application of a composite index analysis (CSS), genomic signatures were identified in Bullmastiffs when compared to the formative breeds, Mastiffs and Bulldogs, and to 22 other breed groups. Significant regions were identified on 15 chromosomes, with the most differentiated regions found on CFA1, CFA9, and CFA18. These regions may reflect genetic drift following establishment of the breed or the effects of selective breeding during development of the modern Bullmastiff. This was supported by analysis of genes from the identified genomic regions, including 458 genes from the multi-clade analysis, which revealed enriched pathways that may be related to characteristic traits and distinct morphology of the breed. The study demonstrates the utility of the CSS method in breed-specific genome analysis and advances our understanding of genetic diversity in Bullmastiff dogs.

Large-Scale Polymorphism Analysis of Dog Leukocyte Antigen Class I and Class II Genes (DLA-88, DLA-12/88L and DLA-DRB1) and Comparison of the Haplotype Diversity between Breeds in Japan

Polymorphisms of canine leukocyte antigen (DLA) class I (DLA-88 and DLA-12/88L) and class II (DLA-DRB1) genes are important for disease susceptibility studies, but information on the genetic diversity among dog breeds is still lacking. To better elucidate the polymorphism and genetic diversity between breeds, we genotyped DLA-88, DLA-12/88L, and DLA-DRB1 loci using 829 dogs of 59 breeds in Japan. Genotyping by Sanger sequencing identified 89, 43, and 61 alleles in DLA-88, DLA-12/88L, and DLA-DRB1 loci, respectively, and a total of 131 DLA-88-DLA-12/88L-DLA-DRB1 haplotypes (88-12/88L-DRB1) were detected more than once. Of the 829 dogs, 198 were homozygotes for one of the 52 different 88-12/88L-DRB1 haplotypes (homozygosity rate: 23.8%). Statistical modeling suggests that 90% of the DLA homozygotes or heterozygotes with one or other of the 52 different 88-12/88L-DRB1 haplotypes within somatic stem cell lines would benefit graft outcome after 88-12/88L-DRB1-matched transplantation. As previously reported for DLA class II haplotypes, the diversity of 88-12/88L-DRB1 haplotypes varied remarkably between breeds but was relatively conserved within most breeds. Therefore, the genetic characteristics of high DLA homozygosity rate and poor DLA diversity within a breed are useful for transplantation therapy, but they may affect biological fitness as homozygosity progresses.

The dogs of Chernobyl: Demographic insights into populations inhabiting the nuclear exclusion zone

The 1986 Chernobyl nuclear disaster initiated a series of catastrophic events resulting in long-term and widespread environmental contamination. We characterize the genetic structure of 302 dogs representing three free-roaming dog populations living within the power plant itself, as well as those 15 to 45 kilometers from the disaster site. Genome-wide profiles from Chernobyl, purebred and free-breeding dogs, worldwide reveal that the individuals from the power plant and Chernobyl City are genetically distinct, with the former displaying increased intrapopulation genetic similarity and differentiation. Analysis of shared ancestral genome segments highlights differences in the extent and timing of western breed introgression. Kinship analysis reveals 15 families, with the largest spanning all collection sites within the radioactive exclusion zone, reflecting migration of dogs between the power plant and Chernobyl City. This study presents the first characterization of a domestic species in Chernobyl, establishing their importance for genetic studies into the effects of exposure to long-term, low-dose ionizing radiation.

Genomic Analysis of the Endangered Fonni's Dog Breed: A Comparison of Genomic and Phenotypic Evaluation Scores

The Fonni's dog is an ancient Sardinian breed for livestock and property guarding. In recent years, the number of new registrations to the breeding book has slumped and, thus, this breed risks being lost forever. This work refocuses attention to the Fonni's dog, analysing its genomic makeup and comparing different phenotypical and genetic evaluation scores. Thirty Fonni's dogs were ranked by their general accordance to the breed typicality (typicality score) and to the provisional standard by official judges (judges' score). They were genotyped with a 230K SNP BeadChip and compared with 379 dogs of 24 breeds. Genomically, the Fonni's dogs placed themselves near shepherd dogs and showed a unique genetic signature, which was used to create the genomic score. This score better correlated with typicality (ρ = 0.69, p < 0.0001) than the judges' score (ρ = 0.63, p = 0.0004), which showed little variability among the included dogs. Hair texture or colour were significantly associated in the three scores. The Fonni's dog is confirmed as a well-distinguished breed, despite being selected mainly for its work abilities. The evaluation criteria used during dog expositions can be improved to increase their variability and include elements typical of the breed. The recovery of the Fonni's dog would be possible only with a shared vision between the Italian kennel club and breeders, and the support of regional programs.

Identification by whole genome sequencing of genes associated with delayed postoperative hemorrhage in Scottish deerhounds

BACKGROUND

Delayed postoperative hemorrhage (DEPOH) is an important health concern for Scottish deerhounds.

HYPOTHESIS/OBJECTIVES

Identify genes associated with DEPOH in Scottish deerhounds.

ANIMALS

Two hundred sixty-nine privately owned Scottish deerhounds.

METHODS

Retrospective case-control study. DEPOH cases and controls were identified through an owner health survey. Genome-wide association analysis was performed using whole genome sequences from 8 cases and 17 controls. All cases and controls were genotyped for selected variants.

RESULTS

Of 269 dogs, 10 met inclusion and exclusion criteria for DEPOH, while 62 controls had undergone similar surgical procedures without DEPOH. Genome-wide association analysis identified a single locus on chromosome 9 spanning 40 genes. One of these genes (SERPINF2 encoding alpha-2 antiplasmin) was directly linked to the pathophysiology of DEPOH. The entire cohort was genotyped for a missense SERPINF2 variant (c.605 C>T; p.A202V). Compared to dogs with the reference C/C genotype, the likelihood of DEPOH was significantly higher for dogs with the T/T genotype (odds ratio [OR] = 1235; 95% confidence interval [CI] = 23-6752; P = 0.0005) and with the C/T genotype (OR = 28; 95% CI = 1.4-542; P = 0.03).

CONCLUSIONS AND CLINICAL IMPORTANCE

SERPINF2 is associated with DEPOH in Scottish deerhounds. Genetic testing might be able to identify dogs that are susceptible to DEPOH.

Genetic distance from wolves affects family dogs' reactions towards howls

Domestication dramatically changes behaviour, including communication, as seen in the case of dogs (Canis familiaris) and wolves (Canis lupus). We tested the hypothesis that domestication may affect an ancient, shared communication form of canids, the howling which seems to have higher individual variation in dogs: the perception and usage of howls may be affected by the genetic relatedness of the breeds to their last common ancestor with wolves ('root distance') and by other individual features like age, sex, and reproductive status. We exposed 68 purebred dogs to wolf howl playbacks and recorded their responses. We identified an interaction between root distance and age on the dogs' vocal and behavioural responses: older dogs from more ancient breeds responded longer with howls and showed more stress behaviours. Our results suggest that domestication impacts vocal behaviour significantly: disintegrating howling, a central, species-specific communication form of canids and gradually eradicating it from dogs' repertoire.

That brachycephalic look: Infant-like facial appearance in short-muzzled dog breeds

Brachycephalic dog breeds are highly popular, yet their conformation-related disorders represent a major welfare concern. It has been suggested that the current popularity of such breeds can be explained by their cute, infant-like facial appearances. The concept of 'kindchenschema' refers to the observation that certain physical features of infant humans and other animals can automatically stimulate positive and nurturant feelings in adult observers. But the proposal that brachycephalic dogs possess heightened 'kindchenschema' facial features, even into adulthood, has never been formally investigated. Here, we hypothesised that relative muzzle shortening across a range of breeds would be associated with known 'kindchenschema' facial features, including a relatively larger forehead, larger eyes and smaller nose. Relative fronto-facial feature sizes in exemplar photographs of adult dogs from 42 popular breeds were measured and associated with existing data on the relative muzzle length and height-at-withers of the same breeds. Our results show that, in adulthood, shorter-muzzled breeds have relatively larger (taller) foreheads and relatively larger eyes (i.e. area of exposed eyeball relative to overall face area) than longer-muzzled breeds, and that this effect is independent of breed size. In sum, brachycephalic dog breeds do show exaggeration of some, but not all, known fronto-facial 'kindchenschema' features, and this may well contribute to their apparently cute appearance and to their current popularity as companion animals. We conclude that the challenge of addressing conformation-related disorders in companion dogs needs to take account of the cute, 'kindchenschema' looks that many owners are likely to be attracted to.

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.

Genomic insight into the influence of selection, crossbreeding, and geography on population structure in poultry

BACKGROUND

In poultry, the population structure of local breeds is usually complex mainly due to unrecorded breeding. Local chicken breeds offer an interesting proxy to understand the complexity of population structure in the context of human-mediated development of diverse morphologies and varieties. We studied 37 traditional Dutch chicken breeds to investigate population structure and the corresponding genomic impact using whole-genome sequence data.

RESULTS

Looking at the genetic differences between breeds, the Dutch chicken breeds demonstrated a complex and admixed subdivided structure. The dissection of this complexity highlighted the influence of selection adhering to management purposes, as well as the role of geographic distance within subdivided breed clusters. Identification of signatures of genetic differentiation revealed genomic regions that are associated with diversifying phenotypic selection between breeds, including dwarf size (bantam) and feather color. In addition, with a case study of a recently developed bantam breed developed by crossbreeding, we provide a genomic perspective on the effect of crossbreeding.

CONCLUSIONS

This study demonstrates the complex population structure of local traditional Dutch chicken, and provides insight into the genomic basis and the factors involved in the formation of this complexity.

A Kmer-based paired-end read de novo assembler and genotyper for canine MHC class I genotyping

The major histocompatibility complex class I (MHC-I) genes are highly polymorphic. MHC-I genotyping is required for determining the peptide epitopes available to an individual's T-cell repertoire. Current genotyping software tools do not work for the dog, due to very limited known canine alleles. To address this, we developed a Kmer-based paired-end read (KPR) de novo assembler and genotyper, which assemble paired-end RNA-seq reads from MHC-I regions into contigs, and then genotype each contig and estimate its expression level. KPR tools outperform other popular software examined in typing new alleles. We used KPR tools to successfully genotype152 dogs from a published dataset. The study discovers 33 putative new alleles, finds dominant alleles in 4 dog breeds, and builds allele diversity and expression landscapes among the 152 dogs. Our software meets a significant need in biomedical research.

Kazakh national dog breed Tazy: What do we know?

The Tazy or Kazakh National sighthound has been officially recognized as the national heritage of Kazakhstan. Comprehensive genetic studies of genetic diversity and population structure that could be used for selection and conservation of this unique dog breed have not been conducted so far. The aim of this study was to determine the genetic structure of the Tazy using microsatellite and SNP markers and to place the breed in the context of the world sighthound breeds. Our results showed that all 19 microsatellite loci examined were polymorphic. The observed number of alleles in the Tazy population varied from 6 (INU030 locus) to 12 (AHT137, REN169D01, AHTh260, AHT121, and FH2054 loci) with a mean of 9.778 alleles per locus. The mean number of effective alleles was 4.869 and ranged from 3.349 f to 4.841. All markers were highly informative (PIC values greater than 0.5) and ranged from 0.543 (REN247M23 locus) to 0.865 (AHT121 locus). The observed and expected heterozygosities in a total population were 0.748 and 0.769 and ranged from 0.746 to 0.750 and 0.656 to 0.769, respectively. Overall, the results confirmed that the Tazy breed has a high level of genetic diversity, no significant inbreeding, and a specific genetic structure. Three gene pools underlie the genetic diversity of the Tazy breed. SNP analysis using the CanineHD SNP array, which contains more than 170,000 SNP markers, showed that the Tazy breed is distinct from other sighthound breeds and genetically related to ancient eastern sighthound breeds sharing the same branch with the Afghan Hound and the Saluki. The results, together with archeological findings, confirm the ancient origin of the breed. The findings can be used for the conservation and international registration of the Tazy dog breed.

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.

Evaluation of Genetic Diversity in Dog Breeds Using Pedigree and Molecular Analysis: A Review

Domestic dogs are important for many economic and social reasons, and they have become a well-known model species for human disease. According to research, dog breeds exhibit significant levels of inbreeding and genetic diversity loss, decreasing the population’s ability to adapt in certain conditions, and indicating the need of conservation strategies. Before the development of molecular markers, pedigree information was used for genetic diversity management. In recent years, genomic tools are frequently applied for accurate estimation of genetic diversity and improved genetic conservation due to incomplete pedigrees and pedigree errors. The most frequently used molecular markers include PCR-based microsatellite markers (STRs) and DNA sequencing-based single-nucleotide polymorphism markers (SNP). The aim of this review was to highlight genetic diversity studies on dog breeds conducted using pedigree and molecular markers, as well as the importance of genetic diversity conservation in increasing the adaptability and survival of dog breed populations.