The first cases of early-onset progressive polyneuropathy appeared in the Alaskan Malamute population in Norway in the late 1970s. Affected dogs were of both sexes and were ambulatory paraparetic, progressing to non-ambulatory tetraparesis. On neurologic examination, affected dogs displayed predominantly laryngeal paresis, decreased postural reactions, decreased spinal reflexes and muscle atrophy. The disease was considered eradicated through breeding programmes but recently new cases have occurred in the Nordic countries and the USA. The N-myc downstream-regulated gene (NDRG1) is implicated in neuropathies with comparable symptoms or clinical signs both in humans and in Greyhound dogs. This gene was therefore considered a candidate gene for the polyneuropathy in Alaskan Malamutes. The coding sequence of the NDRG1 gene derived from one healthy and one affected Alaskan Malamute revealed a non-synonymous G>T mutation in exon 4 in the affected dog that causes a Gly98Val amino acid substitution. This substitution was categorized to be “probably damaging” to the protein function by PolyPhen2 (score: 1.000). Subsequently, 102 Alaskan Malamutes from the Nordic countries and the USA known to be either affected (n = 22), obligate carriers (n = 7) or healthy (n = 73) were genotyped for the SNP using TaqMan. All affected dogs had the T/T genotype, the obligate carriers had the G/T genotype and the healthy dogs had the G/G genotype except for 13 who had the G/T genotype. A protein alignment showed that residue 98 is conserved in mammals and also that the entire NDRG1 protein is highly conserved (94.7%) in mammals. We conclude that the G>T substitution is most likely the mutation that causes polyneuropathy in Alaskan Malamutes. Our characterization of a novel candidate causative mutation for polyneuropathy offers a new canine model that can provide further insight into pathobiology and therapy of human polyneuropathy. Furthermore, selection against this mutation can now be used to eliminate the disease in Alaskan Malamutes.

Primary ciliary dyskinesia (PCD) is a hereditary defect of motile cilia in humans and several domestic animal species. Typical clinical findings are chronic recurrent infections of the respiratory tract and fertility problems. We analyzed an Alaskan Malamute family, in which two out of six puppies were affected by PCD. The parents were unaffected suggesting autosomal recessive inheritance. Linkage and homozygosity mapping defined critical intervals comprising ~118 Mb. Whole genome sequencing of one case and comparison to 601 control genomes identified a disease associated frameshift variant, c.43delA, in the NME5 gene encoding a sparsely characterized protein associated with ciliary function. Nme5^-/- knockout mice exhibit doming of the skull, hydrocephalus and sperm flagellar defects. The genotypes at NME5:c.43delA showed the expected co-segregation with the phenotype in the Alaskan Malamute family. An additional unrelated Alaskan Malamute with PCD and hydrocephalus that became available later in the study was also homozygous mutant at the NME5:c.43delA variant. The mutant allele was not present in more than 1000 control dogs from different breeds. Immunohistochemistry demonstrated absence of the NME5 protein from nasal epithelia of an affected dog. We therefore propose NME5:c.43delA as the most likely candidate causative variant for PCD in Alaskan Malamutes. These findings enable genetic testing to avoid the unintentional breeding of affected dogs in the future. Furthermore, the results of this study identify NME5 as a novel candidate gene for unsolved human PCD and/or hydrocephalus cases. Motile cilia are required for clearing mucous, infectious agents and inhaled dust from the airways. Primary ciliary dyskinesia (PCD) is a hereditary defect of motile cilia. Clinical findings may include recurrent airway infections, fertility problems, and sometimes hydrocephalus. We analyzed an Alaskan Malamute family, in which two out of six puppies were affected by an autosomal recessive form of PCD. Whole genome sequencing of an affected dog identified a one base pair deletion in the NME5 gene, c.43delA, leading to an early frame-shift and premature stop codon. Later in the study, we became aware of a previously published Alaskan Malamute with PCD involving respiratory infections and hydrocephalus. We observed perfect concordance of the NME5 genotypes with the PCD phenotype in all three affected Alaskan Malamutes and more than 1000 controls. The fact that the third case, which had no documented close relationship to the initial two cases, was homozygous for the same rare mutant NME5 allele, strongly supports our hypothesis that NME5:c.43delA causes the PCD phenotype. We confirmed absence of NME5 protein expression in nasal epithelium of an affected dog. Our results enable genetic testing in dogs and identify NME5 as novel candidate gene for unsolved human PCD cases.

Many inherited polyneuropathies (PN) observed in dogs have clinical similarities to the genetically heterogeneous group of Charcot-Marie-Tooth (CMT) peripheral neuropathies in humans. The canine disorders collectively show a variable expression of progressive clinical signs and ages of onset, and different breed prevalences. Previously in the Leonberger breed, a variant highly associated with a juvenile-onset PN was identified in the canine orthologue of a CMT-associated gene. As this deletion in ARHGEF10 (termed LPN1) does not explain all cases, PN in this breed may encompass variants in several genes with similar clinical and histopathological features. A genome-wide comparison of 173 k SNP genotypes of 176 cases, excluding dogs homozygous for the ARHGEF10 variant, and 138 controls, was carried out to detect further PN-associated variants. A single suggestive significant association signal on CFA15 was found. The genome of a PN-affected Leonberger homozygous for the associated haplotype was sequenced and variants in the 7.7 Mb sized critical interval were identified. These variants were filtered against a database of variants observed in 202 genomes of various dog breeds and 3 wolves, and 6 private variants in protein-coding genes, all in complete linkage disequilibrium, plus 92 non-coding variants were revealed. Five of the coding variants were predicted to have low or moderate effect on the encoded protein, whereas a 2 bp deletion in GJA9 results in a frameshift of high impact. GJA9 encodes connexin 59, a connexin gap junction family protein, and belongs to a group of CMT-associated genes that have emerged as important components of peripheral myelinated nerve fibers. The association between the GJA9 variant and PN was confirmed in an independent cohort of 296 cases and 312 controls. Population studies showed a dominant mode of inheritance, an average age of onset of approximately 6 years, and incomplete penetrance. This GJA9 variant represents a highly probable candidate variant for another form of PN in Leonberger dogs, which we have designated LPN2, and a new candidate gene for CMT disease. To date, approximately every third PN-diagnosed Leonberger dog can be explained by the ARHGEF10 or GJA9 variants, and we assume that additional genetic heterogeneity in this condition exists in the breed.

Cutaneous lupus erythematosus (CLE) in humans encompasses multiple subtypes that exhibit a wide array of skin lesions and, in some cases, are associated with the development of systemic lupus erythematosus (SLE). We investigated dogs with exfoliative cutaneous lupus erythematosus (ECLE), a dog-specific form of chronic CLE that is inherited as a monogenic autosomal recessive trait. A genome-wide association study (GWAS) with 14 cases and 29 controls confirmed a previously published result that the causative variant maps to chromosome 18. Autozygosity mapping refined the ECLE locus to a 493 kb critical interval. Filtering of whole genome sequence data from two cases against 654 controls revealed a single private protein-changing variant in this critical interval, UNC93B1:c.1438C>A or p.Pro480Thr. The homozygous mutant genotype was exclusively observed in 23 ECLE affected German Shorthaired Pointers and an ECLE affected Vizsla, but absent from 845 controls. UNC93B1 is a transmembrane protein located in the endoplasmic reticulum and endolysosomes, which is required for correct trafficking of several Toll-like receptors (TLRs). The p.Pro480Thr variant is predicted to affect the C-terminal tail of the UNC93B1 that has recently been shown to restrict TLR7 mediated autoimmunity via an interaction with syndecan binding protein (SDCBP). The functional knowledge on UNC93B1 strongly suggests that p.Pro480Thr is causing ECLE in dogs. These dogs therefore represent an interesting spontaneous model for human lupus erythematosus. Our results warrant further investigations of whether genetic variants affecting the C-terminus of UNC93B1 might be involved in specific subsets of CLE or SLE cases in humans and other species.

Altered melanosome transport in melanocytes, resulting from variants in the melanophilin (MLPH) gene, are associated with inherited forms of coat color dilution in many species. In dogs, the MLPH gene corresponds to the D locus and two variants, c.−22G > A (d¹) and c.705G > C (d²), leading to the dilution of coat color, as described. Here, we describe the independent investigations of dogs whose coat color dilution could not be explained by known variants, and who report a third MLPH variant, (c.667_668insC) (d³), which leads to a frameshift and premature stop codon (p.His223Profs*41). The d³ allele is found at low frequency in multiple dog breeds, as well as in wolves, wolf-dog hybrids, and indigenous dogs. Canids in which the d³ allele contributed to the grey (dilute) phenotype were d¹/d³ compound heterozygotes or d³ homozygotes, and all non-dilute related dogs had one or two D alleles, consistent with a recessive inheritance. Similar to other loci responsible for coat colors in dogs, this, alongside likely additional allelic heterogeneity at the D locus, or other loci, must be considered when performing and interpreting genetic testing.

Dogs represent a unique spontaneous cancer model. Osteosarcoma (OSA) is the most common primary bone tumor in dogs (OMIA 001441-9615), and strongly resembles human forms of OSA. Several large- to giant-sized dog breeds, including the Leonberger, have a greatly increased risk of developing OSA. We performed genome-wide association analysis with high-density imputed SNP genotype data from 273 Leonberger cases with a median age of 8.1 [3.1–13.5] years and 365 controls older than eight years. This analysis revealed significant associations at the CDKN2A/B gene locus on canine chromosome 11, mirroring previous findings in other dog breeds, such as the greyhound, that also show an elevated risk for OSA. Heritability (h²_SNP) was determined to be 20.6% (SE = 0.08; p-value = 5.7 × 10⁻⁴) based on a breed prevalence of 20%. The 2563 SNPs across the genome accounted for nearly all the h²_SNP of OSA, with 2183 SNPs of small effect, 316 SNPs of moderate effect, and 64 SNPs of large effect. As with many other cancers it is likely that regulatory, non-coding variants underlie the increased risk for cancer development. Our findings confirm a complex genetic basis of OSA, moderate heritability, and the crucial role of the CDKN2A/B locus leading to strong cancer predisposition in dogs. It will ultimately be interesting to study and compare the known genetic loci associated with canine OSA in human OSA.