Historic dog furs unravel the origin and artificial selection of modern Nordic Lapphund and Elkhound dog breeds.

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

A new Finnish flavor of feline coat coloration, "salmiak," is associated with a 95-kb deletion downstream of the KIT gene.

Cats with a distinctive white hair pattern of unknown molecular cause have been discovered in the Finnish domestic cat population. Based on the unique appearance of these cats, we have named this phenotype salmiak ("salty licorice"). The use of a commercially available panel test to genotype four salmiak-colored cats revealed the absence of all known variants associated with white-haired phenotypic loci: full White (W), Spotting (Ws) and the Birman white Gloves associated (wg) allele of the KIT proto-oncogene (KIT) gene. Whole-genome sequencing on two salmiak-colored cats was conducted to search for candidate causal variants in the KIT gene. Despite a lack of coding variants, visual inspection of the short read alignments revealed a large ~95 kb deletion located ~65 kb downstream of the KIT gene in the salmiak cats. Additional PCR genotyping of 180 domestic cats and three salmiak-colored cats confirmed the homozygous derived variant genotype fully concordant with the salmiak phenotype. We suggest the newly identified variant be designated as wsal for "w salmiak".

PCYT2 deficiency in Saarlooswolfdogs with progressive retinal, central, and peripheral neurodegeneration.

We investigated a syndromic disease comprising blindness and neurodegeneration in 11 Saarlooswolfdogs. Clinical signs involved early adult onset retinal degeneration and adult-onset neurological deficits including gait abnormalities, hind limb weakness, tremors, ataxia, cognitive decline and behavioral changes such as aggression towards the owner. Histopathology in one affected dog demonstrated cataract, retinal degeneration, central and peripheral axonal degeneration, and severe astroglial hypertrophy and hyperplasia in the central nervous system. Pedigrees indicated autosomal recessive inheritance. We mapped the suspected genetic defect to a 15 Mb critical interval by combined linkage and autozygosity analysis. Whole genome sequencing revealed a private homozygous missense variant, PCYT2:c.4A>G, predicted to change the second amino acid of the encoded ethanolamine-phosphate cytidylyltransferase 2, XP_038402224.1:(p.Ile2Val). Genotyping of additional Saarlooswolfdogs confirmed the homozygous genotype in all eleven affected dogs and demonstrated an allele frequency of 9.9% in the population. This experiment also identified three additional homozygous mutant young dogs without overt clinical signs. Subsequent examination of one of these dogs revealed early-stage progressive retinal atrophy (PRA) and expansion of subarachnoid CSF spaces in MRI. Dogs homozygous for the pathogenic variant showed ether lipid accumulation, confirming a functional PCYT2 deficiency. The clinical and metabolic phenotype in affected dogs shows some parallels with human patients, in whom PCYT2 variants lead to a rare form of spastic paraplegia or axonal motor and sensory polyneuropathy. Our results demonstrate that PCYT2:c.4A>G in dogs cause PCYT2 deficiency. This canine model with histopathologically documented retinal, central, and peripheral neurodegeneration further deepens the knowledge of PCYT2 deficiency.

Idiopathic Epilepsy Risk Allele Trends in Belgian Tervuren: A Longitudinal Genetic Analysis.

Idiopathic epilepsy (IE) has been known to be inherited in the Belgian Tervuren for many decades. Risk genotypes for IE in this breed have recently been identified on Canis familiaris chromosomes (CFA) 14 and 37. In the current study, the allele frequencies of these loci were analyzed to determine whether dog breeders had employed a purposeful selection against IE, leading to a reduction in risk-associated allele frequency within the breed over time. The allele frequencies of two generational groupings of Belgian Tervuren with and without IE were compared. Allele frequencies for risk-associated alleles on CFA14 were unchanged between 1985 and 2015, whereas those on CFA37 increased during that time in the control population (p < 0.05). In contrast, dogs with IE showed a decrease (p < 0.05) in the IE risk-associated allele frequency at the CFA37 locus. Seizure prevalence in the Belgian Tervuren appears to be increasing. These results suggest that, despite awareness that IE is inherited, selection against IE has not been successful.

Identification of novel genetic risk factors of dilated cardiomyopathy: from canine to human.

BACKGROUND: Dilated cardiomyopathy (DCM) is a life-threatening heart disease and a common cause of heart failure due to systolic dysfunction and subsequent left or biventricular dilatation. A significant number of cases have a genetic etiology; however, as a complex disease, the exact genetic risk factors are largely unknown, and many patients remain without a molecular diagnosis. METHODS: We performed GWAS followed by whole-genome, transcriptome, and immunohistochemical analyses in a spontaneously occurring canine model of DCM. Canine gene discovery was followed up in three human DCM cohorts. RESULTS: Our results revealed two independent additive loci associated with the typical DCM phenotype comprising left ventricular systolic dysfunction and dilatation. We highlight two novel candidate genes, RNF207 and PRKAA2, known for their involvement in cardiac action potentials, energy homeostasis, and morphology. We further illustrate the distinct genetic etiologies underlying the typical DCM phenotype and ventricular premature contractions. Finally, we followed up on the canine discoveries in human DCM patients and discovered candidate variants in our two novel genes. CONCLUSIONS: Collectively, our study yields insight into the molecular pathophysiology of DCM and provides a large animal model for preclinical studies.

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.

A frameshift deletion in F8 associated with hemophilia A in Labrador Retriever dogs.

Hemophilia A is the most common inherited coagulation factor disorder in dogs. It manifests as excessive bleeding resulting from pathogenic variants in the X-chromosomal F8 gene encoding coagulation factor VIII (FVIII) protein. In this study, we performed careful clinical phenotyping to confirm hemophilia A in two distinct Labrador Retriever (LR) pedigrees. Whole-genome sequencing on an affected dog from litter 1 identified a case-specific frameshift deletion variant in F8 predicted to cause a premature stop codon (c.2923_2924del, p.(E975Kfs*8)). This variant was hemizygous in all the affected males from litter 1 (n = 3), while all the unaffected LRs in the pedigree were heterozygous or wild-type (n = 22). Additionally, screened samples from 199 LRs were all found to be wild-type. As a result of this study, a gene test can now be developed to screen dogs before breeding to prevent further cases. However, it is important to note that the affected LR with decreased FVIII activity from litter 2 was wild-type for the identified deletion variant, and no segregating F8 variants were detected when this dog's DNA sample was whole-genome sequenced. Thus, the cause of decreased FVIII activity in this dog remains to be unraveled in future studies.

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

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

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.

Validation of a Chromosome 14 Risk Haplotype for Idiopathic Epilepsy in the Belgian Shepherd Dog Found to Be Associated with an Insertion in the RAPGEF5 Gene.

An idiopathic epilepsy (IE) risk haplotype on canine chromosome (CFA) 14 has been reported to interact with the CFA37 common risk haplotype in the Belgian shepherd (BS). Additional IE cases and control dogs were genotyped for the risk haplotypes to validate these previous findings. In the new cohort, the interaction between the two regions significantly elevated IE risk. When the haplotypes were analyzed individually, particular haplotypes on both CFA14 (ACTG) and 37 (GG) were associated with elevated IE risk, though only the CFA37 AA was significantly associated (p < 0.003) with reduced risk in the new cohort. However, the CFA14 ACTG risk was statistically significant when the new and previous cohort data were combined. The frequency of the ACTG haplotype was four-fold higher in BS dogs than in other breeds. Whole genome sequence analysis revealed that a 3-base pair predicted disruptive insertion in the RAPGEF5 gene, which is adjacent to the CFA14 risk haplotype. RAPGEF5 is involved in the Wnt-ß-catenin signaling pathway that is crucial for normal brain function. Although this risk variant does not fully predict the likelihood of a BS developing IE, the association with a variant in a candidate gene may provide insight into the genetic control of canine IE.

Dog colour patterns explained by modular promoters of ancient canid origin.

Distinctive colour patterns in dogs are an integral component of canine diversity. Colour pattern differences are thought to have arisen from mutation and artificial selection during and after domestication from wolves but important gaps remain in understanding how these patterns evolved and are genetically controlled. In other mammals, variation at the ASIP gene controls both the temporal and spatial distribution of yellow and black pigments. Here, we identify independent regulatory modules for ventral and hair cycle ASIP expression, and we characterize their action and evolutionary origin. Structural variants define multiple alleles for each regulatory module and are combined in different ways to explain five distinctive dog colour patterns. Phylogenetic analysis reveals that the haplotype combination for one of these patterns is shared with Arctic white wolves and that its hair cycle-specific module probably originated from an extinct canid that diverged from grey wolves more than 2 million years ago. Natural selection for a lighter coat during the Pleistocene provided the genetic framework for widespread colour variation in dogs and wolves.

Whole Genome Sequencing Reveals Multiple Linked Genetic Variants on Canine Chromosome 12 Associated with Risk for Symmetrical Lupoid Onychodystrophy (SLO) in the Bearded Collie.

In dogs, symmetrical lupoid onychodystrophy (SLO) results in nail loss and an abnormal regrowth of the claws. In Bearded Collies, an autoimmune nature has been suggested because certain dog leukocyte antigen (DLA) class II haplotypes are associated with the condition. A genome-wide association study of the Bearded Collie revealed two regions of association that conferred risk for disease: one on canine chromosome (CFA) 12 that encompasses the DLA genes, and one on CFA17. Case-control association was employed on whole genome sequencing data to uncover putative causative variants in SLO within the CFA12 and CFA17 associated regions. Genotype imputation was then employed to refine variants of interest. Although no SLO-associated protein-coding variants were identified on CFA17, multiple variants, many with predicted damaging effects, were identified within potential candidate genes on CFA12. Furthermore, many potentially damaging alleles were fully correlated with the presence of DLA class II risk haplotypes for SLO, suggesting that the variants may reflect DLA class II haplotype association with disease or vice versa. Strong linkage disequilibrium in the region precluded the ability to isolate and assess the individual or combined effect of variants on disease development. Nonetheless, all were predictive of risk for SLO and, with judicious assessment, their application in selective breeding may prove useful to reduce the incidence of SLO in the breed.

Reliable wolf-dog hybrid detection in Europe using a reduced SNP panel developed for non-invasively collected samples.

BACKGROUND: Understanding the processes that lead to hybridization of wolves and dogs is of scientific and management importance, particularly over large geographical scales, as wolves can disperse great distances. However, a method to efficiently detect hybrids in routine wolf monitoring is lacking. Microsatellites offer only limited resolution due to the low number of markers showing distinctive allele frequencies between wolves and dogs. Moreover, calibration across laboratories is time-consuming and costly. In this study, we selected a panel of 96 ancestry informative markers for wolves and dogs, derived from the Illumina CanineHD Whole-Genome BeadChip (174 K). We designed very short amplicons for genotyping on a microfluidic array, thus making the method suitable also for non-invasively collected samples. RESULTS: Genotypes based on 93 SNPs from wolves sampled throughout Europe, purebred and non-pedigree dogs, and suspected hybrids showed that the new panel accurately identifies parental individuals, first-generation hybrids and first-generation backcrosses to wolves, while second- and third-generation backcrosses to wolves were identified as advanced hybrids in almost all cases. Our results support the hybrid identity of suspect individuals and the non-hybrid status of individuals regarded as wolves. We also show the adequacy of these markers to assess hybridization at a European-wide scale and the importance of including samples from reference populations. CONCLUSIONS: We showed that the proposed SNP panel is an efficient tool for detecting hybrids up to the third-generation backcrosses to wolves across Europe. Notably, the proposed genotyping method is suitable for a variety of samples, including non-invasive and museum samples, making this panel useful for wolf-dog hybrid assessments and wolf monitoring at both continental and different temporal scales.

Clinical and Genetic Findings in 28 American Cocker Spaniels with Aural Ceruminous Gland Hyperplasia and Ectasia.

American Cocker Spaniels (ACSs) develop aural ceruminous gland hyperplasia and ectasia more often than dogs of other breeds. Data on the cause and development of these breed characteristic histopathological changes are lacking. We performed video-otoscopic examinations and dermatological work-up on 28 ACSs, obtained aural biopsies from each dog and assessed the statistical associations between the presence of ceruminous gland hyperplasia and ectasia and disease history, clinical or microbiological findings and underlying cause of otitis externa (OE). Histological lesions of ceruminous gland hyperplasia and ectasia were observed in aural biopsies from 6/13 clinically healthy ears and 13/15 ears with OE from 19/28 examined dogs. Nine of 28 dogs had histologically normal ceruminous glands (odds ratio [OR] 6.2, 95% confidence interval [CI] 1.1-36.6). Bacterial growth in microbiological culture of aural exudate (OR 14.1, 95% CI 2.1-95.3) was associated with ceruminous glandular changes, whereas previous history of OE, cutaneous findings or underlying allergies were not. Pedigree analysis and a genome-wide association study (GWAS) were performed on 18 affected and eight unaffected dogs based on histopathological diagnosis. While the GWAS indicated a tentative, but not statistically significant, association of ceruminous gland hyperplasia and ectasia with chromosome 31, a larger cohort is needed to confirm this preliminary result. Based on our results, ceruminous gland hyperplasia and ectasia may also precede clinical signs of OE in ACSs and a genetic aetiological component is likely Further studies with larger cohorts are warranted to verify our preliminary results.

Missense variant in LOXHD1 is associated with canine nonsyndromic hearing loss.

Hearing loss is a common sensory deficit in both humans and dogs. In canines, the genetic basis is largely unknown, as genetic variants have only been identified for a syndromic form of hearing impairment. We observed a congenital or early-onset sensorineural hearing loss in a Rottweiler litter. Assuming an autosomal recessive inheritance, we used a combined approach of homozygosity mapping and genome sequencing to dissect the genetic background of the disorder. We identified a fully segregating missense variant in LOXHD1, a gene that is known to be essential for cochlear hair cell function and associated with nonsyndromic hearing loss in humans and mice. The canine LOXHD1 variant was specific to the Rottweiler breed in our study cohorts of pure-bred dogs. However, it also was present in some mixed-breed dogs, of which the majority showed Rottweiler ancestry. Low allele frequencies in these populations, 2.6% and 0.04%, indicate a rare variant. To summarize, our study describes the first genetic variant for canine nonsyndromic hearing loss, which is clinically and genetically similar to human LOXHD1-related hearing disorder, and therefore, provides a new large animal model for hearing loss. Equally important, the affected breed will benefit from a genetic test to eradicate this LOXHD1-related hearing disorder from the population.

In-frame deletion in canine PITRM1 is associated with a severe early-onset epilepsy, mitochondrial dysfunction and neurodegeneration.

We investigated the clinical, genetic, and pathological characteristics of a previously unknown severe juvenile brain disorder in several litters of Parson Russel Terriers. The disease started with epileptic seizures at 6-12 weeks of age and progressed rapidly to status epilepticus and death or euthanasia. Histopathological changes at autopsy were restricted to the brain. There was severe acute neuronal degeneration and necrosis diffusely affecting the grey matter throughout the brain with extensive intraneuronal mitochondrial crowding and accumulation of amyloid-ß (Aß). Combined homozygosity mapping and genome sequencing revealed an in-frame 6-bp deletion in the nuclear-encoded pitrilysin metallopeptidase 1 (PITRM1) encoding for a mitochondrial protease involved in mitochondrial targeting sequence processing and degradation. The 6-bp deletion results in the loss of two amino acid residues in the N-terminal part of PITRM1, potentially affecting protein folding and function. Assessment of the mitochondrial function in the affected brain tissue showed a significant deficiency in respiratory chain function. The functional consequences of the mutation were modeled in yeast and showed impaired growth in permissive conditions and an impaired respiration capacity. Loss-of-function variants in human PITRM1 result in a childhood-onset progressive amyloidotic neurological syndrome characterized by spinocerebellar ataxia with behavioral, psychiatric and cognitive abnormalities. Homozygous Pitrm1-knockout mice are embryonic lethal, while heterozygotes show a progressive, neurodegenerative phenotype characterized by impairment in motor coordination and Aß deposits. Our study describes a novel early-onset PITRM1-related neurodegenerative canine brain disorder with mitochondrial dysfunction, Aß accumulation, and lethal epilepsy. The findings highlight the essential role of PITRM1 in neuronal survival and strengthen the connection between mitochondrial dysfunction and neurodegeneration.

A hypomyelinating leukodystrophy in German Shepherd dogs.

BACKGROUND: Shaking puppy syndrome is commonly attributed to abnormal myelination of the central nervous system. HYPOTHESIS/OBJECTIVES: To report the long-term clinical course and the imaging characteristics of hypomyelinating leukodystrophy in German Shepherd dogs. ANIMALS AND METHODS: Three related litters with 11 affected dogs. RESULTS: The 11 affected dogs experienced coarse, side-to-side tremors of the head and trunk, which interfered with normal goal-oriented movements and disappeared at rest. Signs were noticed shortly after birth. Nine dogs were euthanized, 3 dogs underwent pathological examination, and 2 littermates were raised by their breeder. Tremors improved gradually until 6 to 7 months of age. Adult dogs walked with severe residual pelvic limb ataxia. One dog developed epilepsy with tonic-clonic seizures at 15 months of age. Conventional magnetic resonance imaging (MRI) disclosed homogenous hyperintense signal of the entire subcortical white matter in 3 affected 7-week-old dogs and a hypointense signal in a presumably unaffected littermate. Subcortical white matter appeared isointense to gray matter at 15 and 27 weeks of age on repeated MRI. Abnormal white matter signal with failure to display normal gray-white matter contrast persisted into adulthood. Cerebellar arbor vitae was not visible at any time point. Clinical signs, MRI findings, and pathological examinations were indicative of a hypomyelinating leukodystrophy. All parents of the affected litters shared a common ancestor and relatedness of the puppies suggested an autosomal recessive mode of inheritance. CONCLUSION: We describe a novel hypomyelinating leukodystrophy in German Shepherd dogs with a suspected inherited origin.

Intronic variant in POU1F1 associated with canine pituitary dwarfism.

The anterior pituitary gland secretes several endocrine hormones, essential for growth, reproduction and other basic physiological functions. Abnormal development or function of the pituitary gland leads to isolated or combined pituitary hormone deficiency (CPHD). At least 30 genes have been associated with human CPHD, including many transcription factors, such as POU1F1. CPHD occurs spontaneously also in mice and dogs. Two affected breeds have been reported in dogs: German Shepherds with a splice defect in the LHX3 gene and Karelian Bear Dogs (KBD) with an unknown genetic cause. We obtained samples from five KBDs presenting dwarfism and abnormal coats. A combined analysis of genome-wide association and next-generation sequencing mapped the disease to a region in chromosome 31 and identified a homozygous intronic variant in the fourth exon of the POU1F1 gene in the affected dogs. The identified variant, c.605-3C>A, resided in the splice region and was predicted to affect splicing. The variant's screening in three new prospective cases, related breeds, and ~ 8000 dogs from 207 breeds indicated complete segregation in KBDs with a carrier frequency of 8%, and high breed-specificity as carriers were found at a low frequency only in Lapponian Herders, a related breed. Our study establishes a novel canine model for CPHD with a candidate POU1F1 defect.