A comprehensive biomedical variant catalogue based on whole genome sequences of 582 dogs and eight wolves

A Missense Variant in SLC39A4 in a Litter of Turkish Van Cats with Acrodermatitis Enteropathica

In a litter of Turkish Van cats, three out of six kittens developed severe signs of skin disease, diarrhea, and systemic signs of stunted growth at 6 weeks of age. Massive secondary infections of the skin lesions evolved. Histopathological examinations showed a mild to moderate hyperplastic epidermis, covered by a thick layer of laminar to compact, mostly parakeratotic keratin. The dermis was infiltrated with moderate amounts of lymphocytes and plasma cells. Due to the severity of the clinical signs, one affected kitten died and the other two had to be euthanized. We sequenced the genome of one affected kitten and compared the data to 54 control genomes. A search for private variants in the two candidate genes for the observed phenotype, MKLN1 and SLC39A4, revealed a single protein-changing variant, SLC39A4:c.1057G>C or p.Gly353Arg. The solute carrier family 39 member 4 gene (SLC39A4) encodes an intestinal zinc transporter required for the uptake of dietary zinc. The variant is predicted to change a highly conserved glycine residue within the first transmembrane domain, which most likely leads to a loss of function. The genotypes of the index family showed the expected co-segregation with the phenotype and the mutant allele was absent from 173 unrelated control cats. Together with the knowledge on the effects of SLC39A4 variants in other species, these data suggest SLC39A4:c.1057G>C as candidate causative genetic variant for the phenotype in the investigated kittens. In line with the human phenotype, we propose to designate this disease acrodermatitis enteropathica (AE).

Hereditary sensory and autonomic neuropathy in a family of mixed breed dogs associated with a novel RETREG1 variant

Background

Hereditary sensory and autonomic neuropathies (HSANs) are a group of genetic disorders affecting the peripheral nervous system. Two different associated variants have been identified in dogs: 1 in Border Collies and 1 in Spaniels and Pointers.

Objectives

Clinically and genetically characterize HSAN in a family of mixed breed dogs.

Animals

Five 7‐month‐old mixed breed dogs from 2 related litters were presented for evaluation of a 2‐month history of acral mutilation and progressive pelvic limb gait abnormalities.

Methods

Complete physical, neurological, electrodiagnostic, and histopathological evaluations were performed. Whole genome sequencing of 2 affected dogs (1 from each litter) was used to identify variants that were homozygous or heterozygous in both cases, but wild type in 217 control genomes of 100 breeds. Immunohistochemistry was used to assess protein expression.

Results

Complete physical, neurological, electrodiagnostic, and histopathological evaluations confirmed a disorder affecting sensory and autonomic nerves. Whole genome sequencing identified a missense variant in the RETREG1 (reticulophagy regulator 1) gene (c.656C > T, p.P219L). All affected dogs were homozygous for the variant, which was not detected in 1193 dogs from different breeds. Immunohistochemistry showed no expression of RETREG1 in the cerebellum of affected dogs. One of the affected dogs lived for 5 years and showed gradual progression of the clinical signs.

Conclusions and Clinical Importance

We confirmed the diagnosis of HSAN in a family of mixed breed dogs and identified a novel and possibly pathogenic RETREG1 variant. Affected dogs experienced gradual deterioration over several years.

Deletion of the SELENOP gene leads to CNS atrophy with cerebellar ataxia in dogs

We investigated a hereditary cerebellar ataxia in Belgian Shepherd dogs. Affected dogs developed uncoordinated movements and intention tremor at two weeks of age. The severity of clinical signs was highly variable. Histopathology demonstrated atrophy of the CNS, particularly in the cerebellum. Combined linkage and homozygosity mapping in a family with four affected puppies delineated a 52 Mb critical interval. The comparison of whole genome sequence data of one affected dog to 735 control genomes revealed a private homozygous structural variant in the critical interval, Chr4:66,946,539_66,963,863del17,325. This deletion includes the entire protein coding sequence of SELENOP and is predicted to result in complete absence of the encoded selenoprotein P required for selenium transport into the CNS. Genotypes at the deletion showed the expected co-segregation with the phenotype in the investigated family. Total selenium levels in the blood of homozygous mutant puppies of the investigated litter were reduced to about 30% of the value of a homozygous wildtype littermate. Genotyping >600 Belgian Shepherd dogs revealed an additional homozygous mutant dog. This dog also suffered from pronounced ataxia, but reached an age of 10 years. Selenop^-/- knock-out mice were reported to develop ataxia, but their histopathological changes were less severe than in the investigated dogs. Our results demonstrate that deletion of the SELENOP gene in dogs cause a defect in selenium transport associated with CNS atrophy and cerebellar ataxia (CACA). The affected dogs represent a valuable spontaneous animal model to gain further insights into the pathophysiological consequences of CNS selenium deficiency.

We studied a form of inherited ataxia in a family of Belgian Shepherd dogs that we termed CNS atrophy and cerebellar ataxia (CACA). Clinical signs were evident at 2 weeks of age and the affected puppies had to be euthanized at 4 weeks of age. The pedigree of the index family with 4 affected and 4 unaffected puppies suggested autosomal recessive inheritance. Using a purely positional cloning approach, we identified a complete deletion of the SELENOP gene as the most likely causative variant. SELENOP encodes selenoprotein P, a protein with multiple selenocysteine residues, which is required for the transport of selenium into the CNS. Selenium measurements in affected dogs demonstrated blood selenium levels of about 30% compared to normal control dogs. Genotyping a cohort of additional Belgian Shepherd dogs with unexplained ataxia identified another CACA case that had a relatively stable clinical condition and reached an age of 10 years. Selenop^-/- knock-out mice show a related but not identical ataxia phenotype. Our finding of a SELENOP gene deletion in CACA affected dogs identifies a spontaneous animal model to gain further insights into the pathophysiological consequences of CNS selenium deficiency.

Canine junctional epidermolysis bullosa due to a novel mutation in LAMA3 with severe upper respiratory involvement

BACKGROUND

Junctional epidermolysis bullosa (JEB) is a group of congenital blistering skin diseases characterized by clefting through the lamina lucida of the basement membrane zone.

OBJECTIVES

To characterize the clinical and morphological features of a congenital mechanobullous disease in a litter of puppies with severe upper respiratory involvement, and to identify an associated genetic variant.

ANIMALS

Five of eight puppies in an Australian cattle dog cross-bred litter showed signs of skin fragility. Three were stillborn and one died at one month of age. The two surviving puppies were presented with blistering skin disease and severe respiratory distress. Additionally, one unaffected sibling was examined and blood was obtained for genetic testing.

METHODS AND MATERIALS

Post-mortem examination, histopathological evaluation and electron microscopy were performed. Whole genome sequencing (WGS) of one affected puppy was compared to a database of 522 dogs of 55 different breeds for variant analysis. Sanger sequencing of one additional affected and one unaffected sibling confirmed the variant.

RESULTS

Clinically, severe mucocutaneous ulcers occurred in frictional areas with claw sloughing. Histopathological results revealed subepidermal clefts and electron microscopy confirmed the split in the lamina lucida. Post-mortem examination documented extensive pharyngeal and laryngeal lesions with granulation tissue and fibrinous exudate obscuring the airway. Moderate tracheal hypoplasia contributed. The WGS revealed a novel missense variant in the laminin α3-chain XP_537297.2p(Asp2867Val), with an autosomal recessive mode of inheritance.

CONCLUSIONS AND CLINICAL RELEVANCE

A novel variant in LAMA3 caused a generalized and severe phenotype of JEB with an unique clinical presentation of upper airway obstruction.

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

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

Dog10K_Boxer_Tasha_1.0: A Long-Read Assembly of the Dog Reference Genome

The domestic dog has evolved to be an important biomedical model for studies regarding the genetic basis of disease, morphology and behavior. Genetic studies in the dog have relied on a draft reference genome of a purebred female boxer dog named "Tasha" initially published in 2005. Derived from a Sanger whole genome shotgun sequencing approach coupled with limited clone-based sequencing, the initial assembly and subsequent updates have served as the predominant resource for canine genetics for 15 years. While the initial assembly produced a good-quality draft, as with all assemblies produced at the time, it contained gaps, assembly errors and missing sequences, particularly in GC-rich regions, which are found at many promoters and in the first exons of protein-coding genes. Here, we present Dog10K_Boxer_Tasha_1.0, an improved chromosome-level highly contiguous genome assembly of Tasha created with long-read technologies that increases sequence contiguity >100-fold, closes >23,000 gaps of the CanFam3.1 reference assembly and improves gene annotation by identifying >1200 new protein-coding transcripts. The assembly and annotation are available at NCBI under the accession GCF_000002285.5.

Candidate Genetic Modifiers for RPGR Retinal Degeneration

Purpose

To define genetic variants associated with variable severity of X-linked progressive retinal atrophy 1 (XLPRA1) caused by a five-nucleotide deletion in canine RPGR exon ORF15.

Methods

A genome-wide association study (GWAS) was performed in XLPRA1 phenotype informative pedigree. Whole genome sequencing (WGS) was used for mutational analysis of genes within the candidate genomic region. Retinas of normal and mutant dogs were used for gene expression, gene structure, and RNA duplex analyses.

Results

GWAS followed by haplotype phasing identified an approximately 4.6 Mb candidate genomic interval on CFA31 containing seven protein-coding genes expressed in retina (ROBO1, ROBO2, RBM11, NRIP1, HSPA13, SAMSN1, and USP25). Furthermore, we identified and characterized two novel lncRNAs, ROBO1-AS and ROBO2-AS, that display overlapping gene organization with axon guidance pathway genes ROBO1 and ROBO2, respectively, producing sense-antisense gene pairs. Notably, ROBO1-AS and ROBO2-AS act in cis to form lncRNA/mRNA duplexes with ROBO1 and ROBO2, respectively, suggesting important roles for these lncRNAs in the ROBO regulatory network. A subsequent WGS identified candidate genes within the genomic region on CFA31 that might be implicated in modifying severity of XLPRA1. This approach led to discovery of genetic variants in ROBO1, ROBO1-AS, ROBO2-AS, and USP25 that are strongly associated with the XLPRA1 moderate phenotype.

Conclusions

The study provides new insights into the genetic basis of phenotypic variation in severity of RPGR^orf15-associated retinal degeneration. Our findings suggest an important role for ROBO pathways in disease progression further expanding on our previously reported changes of ROBO1 expression in XLPRA1 retinas.

Identification of a variant in NDP associated with X-linked retinal dysplasia in the English cocker spaniel dog

Purpose

Three related male English Cocker Spaniels (ECS) were reported to be congenitally blind. Examination of one of these revealed complete retinal detachment. A presumptive diagnosis of retinal dysplasia (RD) was provided and pedigree analysis was suggestive of an X-linked mode of inheritance. We sought to investigate the genetic basis of RD in this family of ECS.

Methods

Following whole genome sequencing (WGS) of the one remaining male RD-affected ECS, two distinct investigative approaches were employed: a candidate gene approach and a whole genome approach. In the candidate gene approach, COL9A2, COL9A3, NHEJ1, RS1 and NDP genes were investigated based on their known associations with RD and retinal detachment in dogs and humans. In the whole genome approach, affected WGS was compared with 814 unaffected canids to identify candidate variants, which were filtered based on appropriate segregation and predicted pathogenic effects followed by subsequent investigation of gene function. Candidate variants were tested for appropriate segregation in the ECS family and association with disease was assessed using samples from a total of 180 ECS.

Results

The same variant in NDP (c.653_654insC, p.Met114Hisfs*16) that was predicted to result in 15 aberrant amino acids before a premature stop in norrin protein, was identified independently by both approaches and was shown to segregate appropriately within the ECS family. Association of this variant with X-linked RD was significant (P = 0.0056).

Conclusions

For the first time, we report a variant associated with canine X-linked RD. NDP variants are already known to cause X-linked RD, along with other abnormalities, in human Norrie disease. Thus, the dog may serve as a useful large animal model for research.

L2HGDH Missense Variant in a Cat with L-2-Hydroxyglutaric Aciduria

A 7-month-old, spayed female, domestic longhair cat with L-2-hydroxyglutaric aciduria (L-2-HGA) was investigated. The aim of this study was to investigate the clinical signs, metabolic changes and underlying genetic defect. The owner of the cat reported a 4-month history of multiple paroxysmal seizure-like episodes, characterized by running around the house, often in circles, with abnormal behavior, bumping into obstacles, salivating and often urinating. The episodes were followed by a period of disorientation and inappetence. Neurological examination revealed an absent bilateral menace response. Routine blood work revealed mild microcytic anemia but biochemistry, ammonia, lactate and pre- and post-prandial bile acids were unremarkable. MRI of the brain identified multifocal, bilaterally symmetrical and T2-weighted hyperintensities within the prosencephalon, mesencephalon and metencephalon, primarily affecting the grey matter. Urinary organic acids identified highly increased levels of L-2-hydroxyglutaric acid. The cat was treated with the anticonvulsants levetiracetam and phenobarbitone and has been seizure-free for 16 months. We sequenced the genome of the affected cat and compared the data to 48 control genomes. L2HGDH, coding for L-2-hydroxyglutarate dehydrogenase, was investigated as the top functional candidate gene. This search revealed a single private protein-changing variant in the affected cat. The identified homozygous variant, XM_023255678.1:c.1301A>G, is predicted to result in an amino acid change in the L2HGDH protein, XP_023111446.1:p.His434Arg. The available clinical and biochemical data together with current knowledge about L2HGDH variants and their functional impact in humans and dogs allow us to classify the p.His434Arg variant as a causative variant for the observed neurological signs in this cat.

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

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

A Missense Mutation in the KLF7 Gene Is a Potential Candidate Variant for Congenital Deafness in Australian Stumpy Tail Cattle Dogs

Congenital deafness is prevalent among modern dog breeds, including Australian Stumpy Tail Cattle Dogs (ASCD). However, in ASCD, no causative gene has been identified so far. Therefore, we performed a genome-wide association study (GWAS) and whole genome sequencing (WGS) of affected and normal individuals. For GWAS, 3 bilateral deaf ASCDs, 43 herding dogs, and one unaffected ASCD were used, resulting in 13 significantly associated loci on 6 chromosomes, i.e., CFA3, 8, 17, 23, 28, and 37. CFA37 harbored a region with the most significant association (−log₁₀(9.54 × 10⁻²¹) = 20.02) as well as 7 of the 13 associated loci. For whole genome sequencing, the same three affected ASCDs and one unaffected ASCD were used. The WGS data were compared with 722 canine controls and filtered for protein coding and non-synonymous variants, resulting in four missense variants present only in the affected dogs. Using effect prediction tools, two variants remained with predicted deleterious effects within the Heart development protein with EGF like domains 1 (HEG1) gene (NC_006615.3: g.28028412G>C; XP_022269716.1: p.His531Asp) and Kruppel-like factor 7 (KLF7) gene (NC_006619.3: g.15562684G>A; XP_022270984.1: p.Leu173Phe). Due to its function as a regulator in heart and vessel formation and cardiovascular development, HEG1 was excluded as a candidate gene. On the other hand, KLF7 plays a crucial role in the nervous system, is expressed in the otic placode, and is reported to be involved in inner ear development. 55 additional ASCD samples (28 deaf and 27 normal hearing dogs) were genotyped for the KLF7 variant, and the variant remained significantly associated with deafness in ASCD (p = 0.014). Furthermore, 24 dogs with heterozygous or homozygous mutations were detected, including 18 deaf dogs. The penetrance was calculated to be 0.75, which is in agreement with previous reports. In conclusion, KLF7 is a promising candidate gene causative for ASCD deafness.

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.

Chromosome-length genome assembly and structural variations of the primal Basenji dog (Canis lupus familiaris) genome

BACKGROUND

Basenjis are considered an ancient dog breed of central African origins that still live and hunt with tribesmen in the African Congo. Nicknamed the barkless dog, Basenjis possess unique phylogeny, geographical origins and traits, making their genome structure of great interest. The increasing number of available canid reference genomes allows us to examine the impact the choice of reference genome makes with regard to reference genome quality and breed relatedness.

RESULTS

Here, we report two high quality de novo Basenji genome assemblies: a female, China (CanFam_Bas), and a male, Wags. We conduct pairwise comparisons and report structural variations between assembled genomes of three dog breeds: Basenji (CanFam_Bas), Boxer (CanFam3.1) and German Shepherd Dog (GSD) (CanFam_GSD). CanFam_Bas is superior to CanFam3.1 in terms of genome contiguity and comparable overall to the high quality CanFam_GSD assembly. By aligning short read data from 58 representative dog breeds to three reference genomes, we demonstrate how the choice of reference genome significantly impacts both read mapping and variant detection.

CONCLUSIONS

The growing number of high-quality canid reference genomes means the choice of reference genome is an increasingly critical decision in subsequent canid variant analyses. The basal position of the Basenji makes it suitable for variant analysis for targeted applications of specific dog breeds. However, we believe more comprehensive analyses across the entire family of canids is more suited to a pangenome approach. Collectively this work highlights the importance the choice of reference genome makes in all variation studies.

Web resource on available DNA variant tests for hereditary diseases and genetic predispositions in dogs and cats: An Update

Vast progress has been made in the clinical diagnosis and molecular basis of hereditary diseases and genetic predisposition in companion animals. The purpose of this report is to provide an update on the availability of DNA testing for hereditary diseases and genetic predispositions in dogs and cats utilizing the WSAVA-PennGen DNA Testing Database web resource (URL: http://research.vet.upenn.edu/WSAVA-LabSearch ). Information on hereditary diseases, DNA tests, genetic testing laboratories and afflicted breeds added to the web-based WSAVA-PennGen DNA Testing Database was gathered. Following verification through original research and clinical studies, searching various databases on hereditary diseases in dogs and cats, and contacting laboratories offering DNA tests, the data were compared to the resource reported on in 2013. The number of molecularly defined Mendelian inherited diseases and variants in companion animals listed in the WSAVA-PennGen DNA Testing Database in 2020 drastically increased by 112% and 141%, respectively. The number of DNA variant tests offered by each laboratory has also doubled for dogs and cats. While the overall number of laboratories has only slightly increased from 43 to 47, the number of larger corporate laboratories increased, while academic laboratories have declined. In addition, there are now several laboratories that are offering breed-specific or all-breed panel tests rather than single-DNA tests for dogs and cats. This unique regularly updated searchable web-based database allows veterinary clinicians, breeders and pet owners to readily find available DNA tests, laboratories performing these DNA tests worldwide, and canine and feline breeds afflicted and also serves as a valuable resource for comparative geneticists.

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.

Roan, ticked and clear coat patterns in the canine are associated with three haplotypes near usherin on CFA38

White coat patterning is a feature of many dog breeds and is known to be coded primarily by the gene micropthalmia‐associated transcription factor (MITF). This patterning in the coat can be modified by other factors to produce the attractive phenotypes termed ‘ticked’ and ‘roan’ that describe the presence of flecks of color that vary in distribution and intensity within otherwise ‘clear’ white markings. The appearance of the pigment in the white patterning caused by ticking and roaning intensifies in the weeks after birth. We applied genome‐wide association to compare English Cocker Spaniels of roan phenotype (N = 34) with parti‐color (non‐roan) English Cocker Spaniels (N = 9) and identified an associated locus on CFA 38, CFA38:11 057 040 (P_raw = 8.9 × 10⁻¹⁰, P_genome = 2.7 × 10⁻⁵). A local case–control association in English Springer Spaniels comparing 11 ticked and six clear dogs identified indicative association with a different haplotype, CFA38:11 122 467G>T (P_raw = 1.7 × 10⁻⁵) and CFA38:11 124 294A>C (P_raw = 1.7 × 10⁻⁵). We characterize three haplotypes in Spaniels according to their putative functional variant profiles at CFA38:11 111 286C>T (missense), CFA38:11 131 841–11 143 239DUP.insTTAA (using strongly linked marker CFA38:11 143 243C>T) and CFA38:11 156 425T>C (splice site). In Spaniels, the haplotypes work as an allelic series including alleles (t, recessive clear; T, dominant ticked/parti‐color; and T_R, incomplete dominant roan) to control the appearance of pigmented spots or flecks in otherwise white areas of the canine coat. In Spaniels the associated haplotypes are t (CCT), T (TCC) and T_R (TTT) for SNP markers on CFA38 at 11 111 286C>T, 11 143 243C>T and 11 156 425T>C respectively. It is likely that other alleles exist in this series and together the haplotypes result in a complex range of patterning that is only visible when dogs have white patterning resulting from the epistatic gene Micropthalmia‐associated transcription factor (the S‐locus).

A SINE Insertion in F8 Gene Leads to Severe Form of Hemophilia A in a Family of Rhodesian Ridgebacks

Hemophilia A is the most common coagulation factor disorder in humans and dogs. The disease is characterized by the lack or diminished activity of Factor VIII (FVIII), caused by variants in the F8 gene and inherited as an X chromosomal trait. Two related male Rhodesian Ridgebacks were diagnosed with Hemophilia A due to reduced FVIII activity. The purpose of the study was to determine the genetic cause and give breeding advice for the remaining family members in order to eradicate the variant. By Sanger sequencing a short interspersed nuclear element (SINE) insertion in exon 14 of the F8 gene was found. Perfect correlation of this genetic variant with clinical signs of hemophilia A in the family tree, and the lack of this genetic variant in more than 500 unrelated dogs of the same and other breeds, confirms the hypothesis of this SINE being the underlying genetic cause of Hemophilia A in this family. The identification of clinically unaffected female carriers allows subsequent exclusion of these animals from breeding, to avoid future production of clinically affected male offspring and more subclinical female carriers.

Genetics of canine diabetes mellitus part 2: Current understanding and future directions

Part 1 of this 2-part review outlined the importance of disease classification in diabetes genetic studies, as well as the ways in which genetic variants may contribute to risk of a complex disease within an individual, or within a particular group of individuals. Part 2, presented here, describes in more detail our current understanding of the genetics of canine diabetes mellitus compared to our knowledge of the human disease. Ongoing work to improve our knowledge, using new technologies, is also introduced.

Sarcoglycan A mutation in miniature dachshund dogs causes limb-girdle muscular dystrophy 2D

BACKGROUND

A cohort of related miniature dachshund dogs with exercise intolerance, stiff gait, dysphagia, myoglobinuria, and markedly elevated serum creatine kinase activities were identified.

METHODS

Muscle biopsy histopathology, immunofluorescence microscopy, and western blotting were combined to identify the specific pathologic phenotype of the myopathy, and whole genome SNP array genotype data and whole genome sequencing were combined to determine its genetic basis.

RESULTS

Muscle biopsies were dystrophic. Sarcoglycanopathy, a form of limb-girdle muscular dystrophy, was suspected based on immunostaining and western blotting, where α, β, and γ-sarcoglycan were all absent or reduced. Genetic mapping and whole genome sequencing identified a premature stop codon mutation in the sarcoglycan A subunit gene (SGCA). Affected dachshunds were confirmed on several continents.

CONCLUSIONS

This first SGCA mutation found in dogs adds to the literature of genetic bases of canine muscular dystrophies and their usefulness as comparative models of human disease.

Special Issue "Molecular Basis of Inherited Diseases in Companion Animals"

The study of inherited diseases in companion animals has exploded over the past 15 years since the publication of the first dog genome in 2005 [...].

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.

X-linked CD40 ligand deficiency in a 1-year-old male Shih Tzu with secondary Pneumocystis pneumonia

An approximately 1‐year‐old male intact Shih Tzu dog was referred to a tertiary facility with a history of progressive tachypnea, increased respiratory effort, and weight loss over a 3‐month period that failed to improve with empirical antimicrobial treatment. Upon completion of a comprehensive respiratory evaluation, the dog was diagnosed with severe Pneumocystis pneumonia and secondary pulmonary hypertension. Clinical signs resolved and disease resolution was confirmed after completion of an 8‐week course of trimethoprim‐sulfonamide, 4‐week tapering dose of prednisone to decrease an inflammatory response secondary to acute die‐off of organisms, a 2‐week course of clopidogrel to prevent clot formation, and a 2‐week course of a phosphodiesterase‐5 inhibitor to treat pulmonary hypertension. Immunodiagnostic testing and genetic sequencing were performed to evaluate for potential immunodeficiency as an underlying cause for the development Pneumocystis pneumonia, and identified an X‐linked CD40 ligand deficiency.

A COL7A1 Variant in a Litter of Neonatal Basset Hounds with Dystrophic Epidermolysis Bullosa

We investigated three neonatal Basset Hound littermates with lesions consistent with epidermolysis bullosa (EB), a group of genetic blistering diseases. A clinically normal bitch was bred to her grandfather by artificial insemination. Out of a litter of seven puppies, two affected puppies died and one was euthanized, with these puppies being submitted for diagnostic necropsy. All had multiple bullae and ulcers involving the nasal planum and paw pads, as well as sloughing claws; one puppy also had oral and esophageal ulcers. The complete genome of one affected puppy was sequenced, and 37 known EB candidate genes were assessed. We found a candidate causative variant in COL7A1, which encodes the collagen VII alpha 1 chain. The variant is a complex rearrangement involving duplication of a 107 bp region harboring a frameshift deletion of 7 bp. The variant is predicted to truncate more than 75% of the open reading frame, p.(Val677Serfs*11). Targeted genotyping of this duplication confirmed that all three affected puppies were homozygous for the duplication, whereas 12 unaffected Basset Hounds did not carry the duplication. This variant was also not seen in the genomes of more than 600 dogs of other breeds. COL7A1 variants have been identified in humans and dogs with dystrophic epidermolysis bullosa (DEB). The identified COL7A1 variant therefore most likely represents the causative variant and allows the refinement of the preliminary EB diagnosis to DEB.

CCDC66 frameshift variant associated with a new form of early-onset progressive retinal atrophy in Portuguese Water Dogs

Aberrant photoreceptor function or morphogenesis leads to blinding retinal degenerative diseases, the majority of which have a genetic aetiology. A variant in PRCD previously identified in Portuguese Water Dogs (PWDs) underlies prcd (progressive rod-cone degeneration), an autosomal recessive progressive retinal atrophy (PRA) with a late onset at 3–6 years of age or older. Herein, we have identified a new form of early-onset PRA (EOPRA) in the same breed. Pedigree analysis suggested an autosomal recessive inheritance. Four PWD full-siblings affected with EOPRA diagnosed at 2–3 years of age were genotyped (173,661 SNPs) along with 2 unaffected siblings, 2 unaffected parents, and 15 unrelated control PWDs. GWAS, linkage analysis and homozygosity mapping defined a 26-Mb candidate region in canine chromosome 20. Whole-genome sequencing in one affected dog and its obligatory carrier parents identified a 1 bp insertion (CFA20:g.33,717,704_33,717,705insT (CanFam3.1); c.2262_c.2263insA) in CCDC66 predicted to cause a frameshift and truncation (p.Val747SerfsTer8). Screening of an extended PWD population confirmed perfect co-segregation of this genetic variant with the disease. Western blot analysis of COS-1 cells transfected with recombinant mutant CCDC66 expression constructs showed the mutant transcript translated into a truncated protein. Furthermore, in vitro studies suggest that the mutant CCDC66 is mislocalized to the nucleus relative to wild type CCDC66. CCDC66 variants have been associated with inherited retinal degenerations (RDs) including canine and murine ciliopathies. As genetic variants affecting the primary cilium can cause ciliopathies in which RD may be either the sole clinical manifestation or part of a syndrome, our findings further support a role for CCDC66 in retinal function and viability, potentially through its ciliary function.

Examination of the efficacy of small genetic panels in genomic conservation of companion animal populations

In many ways, dogs are an ideal model for the study of genetic erosion and population recovery, problems of major concern in the field of conservation genetics. Genetic diversity in many dog breeds has been declining systematically since the beginning of the 1800s, when modern breeding practices came into fashion. As such, inbreeding in domestic dog breeds is substantial and widespread and has led to an increase in recessive deleterious mutations of high effect as well as general inbreeding depression. Pedigrees can in theory be used to guide breeding decisions, though are often incomplete and do not reflect the full history of inbreeding. Small microsatellite panels are also used in some cases to choose mating pairs to produce litters with low levels of inbreeding. However, the long-term impact of such practices has not been thoroughly evaluated. Here, we use forward simulation on a model of the dog genome to examine the impact of using limited marker panels to guide pairwise mating decisions on genome-wide population-level genetic diversity. Our results suggest that in unmanaged populations, where breeding decisions are made at the pairwise-rather than population-level, such panels can lead to accelerated loss of genetic diversity at genome regions unlinked to panel markers, compared to random mating. These results demonstrate the importance of genome-wide genetic panels for managing and conserving genetic diversity in dogs and other companion animals.

A CNTNAP1 Missense Variant Is Associated with Canine Laryngeal Paralysis and Polyneuropathy

Laryngeal paralysis associated with a generalized polyneuropathy (LPPN) most commonly exists in geriatric dogs from a variety of large and giant breeds. The purpose of this study was to discover the underlying genetic and molecular mechanisms in a younger-onset form of this neurodegenerative disease seen in two closely related giant dog breeds, the Leonberger and Saint Bernard. Neuropathology of an affected dog from each breed showed variable nerve fiber loss and scattered inappropriately thin myelinated fibers. Using across-breed genome-wide association, haplotype analysis, and whole-genome sequencing, we identified a missense variant in the CNTNAP1 gene (c.2810G>A; p.Gly937Glu) in which homozygotes in both studied breeds are affected. CNTNAP1 encodes a contactin-associated protein important for organization of myelinated axons. The herein described likely pathogenic CNTNAP1 variant occurs in unrelated breeds at variable frequencies. Individual homozygous mutant LPPN-affected Labrador retrievers that were on average four years younger than dogs affected by geriatric onset laryngeal paralysis polyneuropathy could be explained by this variant. Pathologic changes in a Labrador retriever nerve biopsy from a homozygous mutant dog were similar to those of the Leonberger and Saint Bernard. The impact of this variant on health in English bulldogs and Irish terriers, two breeds with higher CNTNAP1 variant allele frequencies, remains unclear. Pathogenic variants in CNTNAP1 have previously been reported in human patients with lethal congenital contracture syndrome and hypomyelinating neuropathy, including vocal cord palsy and severe respiratory distress. This is the first report of contactin-associated LPPN in dogs characterized by a deleterious variant that most likely predates modern breed establishment.

P2Y2 and P2X4 Receptors Mediate Ca2+ Mobilization in DH82 Canine Macrophage Cells

Purinergic receptors of the P2 subclass are commonly found in human and rodent macrophages where they can be activated by adenosine 5'-triphosphate (ATP) or uridine 5'-triphosphate (UTP) to mediate Ca²⁺ mobilization, resulting in downstream signalling to promote inflammation and pain. However, little is understood regarding these receptors in canine macrophages. To establish a macrophage model of canine P2 receptor signalling, the expression of these receptors in the DH82 canine macrophage cell line was determined by reverse transcription polymerase chain reaction (RT-PCR) and immunocytochemistry. P2 receptor function in DH82 cells was pharmacologically characterised using nucleotide-induced measurements of Fura-2 AM-bound intracellular Ca²⁺. RT-PCR revealed predominant expression of P2X4 receptors, while immunocytochemistry confirmed predominant expression of P2Y₂ receptors, with low levels of P2X4 receptor expression. ATP and UTP induced robust Ca²⁺ responses in the absence or presence of extracellular Ca²⁺. ATP-induced responses were only partially inhibited by the P2X4 receptor antagonists, 2',3'-O-(2,4,6-trinitrophenyl)-ATP (TNP-ATP), paroxetine and 5-BDBD, but were strongly potentiated by ivermectin. UTP-induced responses were near completely inhibited by the P2Y₂ receptor antagonists, suramin and AR-C118925. P2Y₂ receptor-mediated Ca²⁺ mobilization was inhibited by U-73122 and 2-aminoethoxydiphenyl borate (2-APB), indicating P2Y₂ receptor coupling to the phospholipase C and inositol triphosphate signal transduction pathway. Together this data demonstrates, for the first time, the expression of functional P2 receptors in DH82 canine macrophage cells and identifies a potential cell model for studying macrophage-mediated purinergic signalling in inflammation and pain in dogs.

NSDHL Frameshift Deletion in a Mixed Breed Dog with Progressive Epidermal Nevi

Loss-of-function variants in the NSDHL gene have been associated with epidermal nevi in humans with congenital hemidysplasia, ichthyosiform nevi, and limb defects (CHILD) syndrome and in companion animals. The NSDHL gene codes for the NAD(P)-dependent steroid dehydrogenase-like protein, which is involved in cholesterol biosynthesis. In this study, a female Chihuahua cross with a clinical and histological phenotype consistent with progressive epidermal nevi is presented. All exons of the NSDHL candidate gene were amplified by PCR and analyzed by Sanger sequencing. A heterozygous frameshift variant, c.718_722delGAACA, was identified in the affected dog. In lesional skin, the vast majority of NSDHL transcripts lacked the five deleted bases. The variant is predicted to produce a premature stop codon truncating 34% of the encoded protein, p.Glu240Profs*17. The mutant allele was absent from 22 additionally genotyped Chihuahuas, as well as from 647 control dogs of diverse breeds and eight wolves. The available experimental data together with current knowledge about NSDHL variants and their functional impact in humans, dogs, and other species prompted us to classify this variant as pathogenic according to the ACMG guidelines that were previously established for human sequence variants. Therefore, we propose the c.718_722delGAACA variant as causative variant for the observed skin lesions in this dog.

SLC19A3 Loss-of-Function Variant in Yorkshire Terriers with Leigh-Like Subacute Necrotizing Encephalopathy

Sporadic occurrence of juvenile-onset necrotizing encephalopathy (SNE) has been previously reported in Yorkshire terriers. However, so far, no causative genetic variant has been found for this breed-specific form of suspected mitochondrial encephalomyopathy. Affected dogs showed gait abnormalities, central visual defects, and/or seizures. Histopathological analysis revealed the presence of major characteristics of human Leigh syndrome and SNE in Alaskan huskies. The aim of this study was to characterize the genetic etiology of SNE-affected purebred Yorkshire terriers. After SNP genotyping and subsequent homozygosity mapping, we identified a single loss-of-function variant by whole-genome sequencing in the canine SLC19A3 gene situated in a 1.7 Mb region of homozygosity on chromosome 25. All ten cases were homozygous carriers of a mutant allele, an indel variant in exon 2, that is predicted to lead to a frameshift and to truncate about 86% of the wild type coding sequence. This study reports a most likely pathogenic variant in SLC19A3 causing a form of SNE in Yorkshire terriers and enables selection against this fatal neurodegenerative recessive disorder. This is the second report of a pathogenic alteration of the SLC19A3 gene in dogs with SNE.

Genomic diversity and population structure of the Leonberger dog breed

BACKGROUND

Leonberger is a giant dog breed formed in the 1850s in Germany. Its post-World War II popularity has resulted in a current global population of ~ 30,000 dogs. The breed has predispositions to neurodegenerative disorders and cancer, which is likely due in large part to limited genetic diversity. However, to date there is no scientific literature on the overall demography and genomic architecture of this breed.

RESULTS

We assessed extensive pedigree records, SNP array genotype data, and whole-genome sequences (WGS) on 142,072, 1203 and 39 Leonberger dogs, respectively. Pedigree analyses identified 22 founder animals and revealed an apparent popular sire effect. The average pedigree-based inbreeding coefficient of 0.29 and average kinship of 0.31 show a dramatic loss of genetic diversity. The observed average life span decreased over time from 9.4 years in 1989 to 7.7 years in 2004. A global health survey confirmed a high prevalence of cancer and neurological disorders. Analysis of SNP-based runs of homozygosity (ROH) identified 125,653 ROH with an average length of 5.88 Mb, and confirmed an average inbreeding coefficient of 0.28. Genome-wide filtering of the WGS data revealed 28 non-protein-changing variants that were present in all Leonberger individuals and a list of 22 potentially pathogenic variants for neurological disorders of which 50% occurred only in Leonbergers and 50% occurred rarely in other breeds. Furthermore, one of the two mtDNA haplogroups detected was present in one dog only.

CONCLUSIONS

The increasing size of the Leonberger population has been accompanied by a considerable loss of genetic diversity after the bottleneck that occurred in the 1940s due to the intensive use of popular sires resulting in high levels of inbreeding. This might explain the high prevalence of certain disorders; however, genomic data provide no evidence for fixed coding variants that explain these predispositions. The list of candidate causative variants for polyneuropathy needs to be further evaluated. Preserving the current genetic diversity is possible by increasing the number of individuals for breeding while restricting the number of litters per sire/dam. In addition, outcrossing would help optimize long-term genetic diversity and contribute to the sustainability and health of the population.

X-Linked Duchenne-Type Muscular Dystrophy in Jack Russell Terrier Associated with a Partial Deletion of the Canine DMD Gene

A 9-month old male Jack Russell Terrier started showing paraparesis of the hindlimbs after a walk. Hospitalized, the dog went into cardiac arrest, and later died. Necroscopic examination revealed a severe thickness of the diaphragm, esophagus, and base of the tongue, leading to the diagnosis of muscular dystrophy. The histology confirmed the marked size variation, regeneration, and fibrosis replacement of the skeletal muscle fibers. Immunohistochemistry demonstrated the absence of dystrophin confirming the diagnosis. Transmission electron microscopy showed disarrangement of skeletal muscle fibers. Finally, whole-genome sequencing identified a ~368kb deletion spanning 19 exons of the canine dystrophin (DMD) gene. This pathogenic loss-of-function variant most likely explains the observed disease phenotype. The X-chromosomal variant was absent in seven controls of the same breed. Most likely, this partial deletion of the DMD gene was either transmitted on the maternal path within the family of the affected dog or arose de novo. This study revealed a spontaneous partial deletion in DMD gene in a Jack Russell Terrier showing a Duchenne-type muscular dystrophy due to non-functional dystrophin.

Sequence analysis of the coding regions of the apolipoprotein C2 (APOC2) gene in Miniature Schnauzers with idiopathic hypertriglyceridemia

It has been hypothesized that idiopathic hypertriglyceridemia in Miniature Schnauzers is hereditary, but the gene responsible has yet to be identified. The objective of this study was to determine if there were coding variants in the apolipoprotein C2 (APOC2) gene in Miniature Schnauzers with idiopathic hypertriglyceridemia. Blood samples from 12 Miniature Schnauzers with idiopathic hypertriglyceridemia were analyzed. Genomic DNA was extracted from whole blood, and the three coding exons of APOC2 were amplified by PCR. The PCR amplicons were sequenced and analyzed for variants relative to the canine reference genome (CanFam3.1 assembly). A second objective was to determine the extent of variation in coding exons of APOC2 in a large and diverse canine population using the Dog Biomedical Variant Database Consortium variant catalog, comprised of whole genome sequencing variant calls from 582 dogs of 126 breeds and eight wolves. There were no variants detected in the coding exons of APOC2 for any of the 12 Miniature Schnauzers with idiopathic hypertriglyceridemia. Variants in the coding exons of APOC2 were also rare in the Dog Biomedical Variant Database Consortium variant catalog; a single synonymous variant was identified in a heterozygous state in a Tibetan Mastiff. Thus, we concluded that coding variants in APOC2 are unlikely to be a major cause of idiopathic hypertriglyceridemia in North American Miniature Schnauzers and furthermore, that such coding variants are rare in the canine population.

ARSB gene variants causing Mucopolysaccharidosis VI in Miniature Pinscher and Miniature Schnauzer dogs

Mucopolysaccharidosis (MPS) VI is a lysosomal storage disease caused by a deficiency of N-acetylgalactosamine-4-sulfatase, also called arylsulfatase B (ARSB, EC 3.1.6.12). Dogs with MPS VI show progressive predominantly oculoskeletal signs homologous to those in human and feline patients. We report herein two pathogenic ARSB gene variants in Miniature Pinscher and Miniature Schnauzer dogs with MPS VI and a genotyping survey in these breeds. All exons and adjacent regions of the ARSB gene were sequenced from three affected Miniature Pinschers and three affected Miniature Schnauzers. Allelic discrimination assays were used for genotyping. A missense variant (NM_001048133.1:c.910G>A) was found in exon 5 of MPS VI-affected Miniature Pinschers that is predicted to result in a deleterious amino acid substitution of a highly conserved glycine to arginine (NP_001041598.1:p.Gly304Arg). In MPS VI-affected Miniature Schnauzers, a 56 bp deletion (NM_001048133.1:c.-24_32del) was found at the junction of exon 1 and its upstream region, predicting no enzyme synthesis. All clinically affected Miniature Pinschers and Miniature Schnauzers were homozygous for the respective variants, and screened healthy dogs in each breed were either heterozygous or homozygous for the wt allele. Whereas the Miniature Pinscher variant seemed to occur commonly (0.133 allele frequency), the Miniature Schnauzer variant was presumed to be rare. In conclusion, two breed-specific pathogenic ARSB gene variants were identified in Miniature Pinscher and Miniature Schnauzer dogs with MPS VI, allowing for genotyping and informed breeding to prevent the production of affected offspring.

LAMB3 Missense Variant in Australian Shepherd Dogs with Junctional Epidermolysis Bullosa

In a highly inbred Australian Shepherd litter, three of the five puppies developed widespread ulcers of the skin, footpads, and oral mucosa within the first weeks of life. Histopathological examinations demonstrated clefting of the epidermis from the underlying dermis within or just below the basement membrane, which led to a tentative diagnosis of junctional epidermolysis bullosa (JEB) with autosomal recessive inheritance. Endoscopy in one affected dog also demonstrated separation between the epithelium and underlying tissue in the gastrointestinal tract. As a result of the severity of the clinical signs, all three dogs had to be euthanized. We sequenced the genome of one affected puppy and compared the data to 73 control genomes. A search for private variants in 37 known candidate genes for skin fragility phenotypes revealed a single protein-changing variant, LAMB3:c.1174T>C, or p.Cys392Arg. The variant was predicted to change a conserved cysteine in the laminin β3 subunit of the heterotrimeric laminin-322, which mediates the binding of the epidermal basement membrane to the underlying dermis. Loss-of-function variants in the human LAMB3 gene lead to recessive forms of JEB. We confirmed the expected co-segregation of the genotypes in the Australian Shepherd family. The mutant allele was homozygous in two genotyped cases and heterozygous in three non-affected close relatives. It was not found in 242 other controls from the Australian Shepherd breed, nor in more than 600 other controls. These data suggest that LAMB3:c.1174T>C represents the causative variant. To the best of our knowledge, this study represents the first report of a LAMB3-related JEB in domestic animals.

A Missense Variant in ALDH5A1 Associated with Canine Succinic Semialdehyde Dehydrogenase Deficiency (SSADHD) in the Saluki Dog

Dogs provide highly valuable models of human disease due to the similarity in phenotype presentation and the ease of genetic analysis. Seven Saluki puppies were investigated for neurological abnormalities including seizures and altered behavior. Magnetic resonance imaging showed a diffuse, marked reduction in cerebral cortical thickness, and symmetrical T2 hyperintensity in specific brain regions. Cerebral cortical atrophy with vacuolation (status spongiosus) was noted on necropsy. Genome-wide association study of 7 affected and 28 normal Salukis revealed a genome-wide significantly associated region on CFA 35. Whole-genome sequencing of three confirmed cases from three different litters revealed a homozygous missense variant within the aldehyde dehydrogenase 5 family member A1 (ALDH5A1) gene (XM_014110599.2: c.866G>A; XP_013966074.2: p.(Gly288Asp). ALDH5A1 encodes a succinic semialdehyde dehydrogenase (SSADH) enzyme critical in the gamma-aminobutyric acid neurotransmitter (GABA) metabolic pathway. Metabolic screening of affected dogs showed markedly elevated gamma-hydroxybutyric acid in serum, cerebrospinal fluid (CSF) and brain, and elevated succinate semialdehyde in urine, CSF and brain. SSADH activity in the brain of affected dogs was low. Affected Saluki dogs had striking similarities to SSADH deficiency in humans although hydroxybutyric aciduria was absent in affected dogs. ALDH5A1-related SSADH deficiency in Salukis provides a unique translational large animal model for the development of novel therapeutic strategies.

A deletion of IDUA exon 10 in a family of Golden Retriever dogs with an attenuated form of mucopolysaccharidosis type I

Background

Mucopolysaccharidosis type I (MPS‐I) is a lysosomal storage disorder caused by a deficiency of the enzyme α‐l‐iduronidase, leading to accumulation of undegraded dermatan and heparan sulfates in the cells and secondary multiorgan dysfunction. In humans, depending upon the nature of the underlying mutation(s) in the IDUA gene, the condition presents with a spectrum of clinical severity.

Objectives

To characterize the clinical and biochemical phenotypes, and the genotype of a family of Golden Retriever dogs.

Animals

Two affected siblings and 11 related dogs.

Methods

Family study. Urine metabolic screening and leucocyte lysosomal enzyme activity assays were performed for biochemical characterization. Whole genome sequencing was used to identify the causal mutation.

Results

The clinical signs shown by the proband resemble the human attenuated form of the disease, with a dysmorphic appearance, musculoskeletal, ocular and cardiac defects, and survival to adulthood. Urinary metabolic studies identified high levels of dermatan sulfate, heparan sulfate, and heparin. Lysosomal enzyme activities demonstrated deficiency in α‐l‐iduronidase activity in leucocytes. Genome sequencing revealed a novel homozygous deletion of 287 bp resulting in full deletion of exon 10 of the IDUA gene (NC_006585.3(NM_001313883.1):c.1400‐76_1521+89del). Treatment with pentosan polyphosphate improved the clinical signs until euthanasia at 4.5 years.

Conclusion and Clinical Importance

Analysis of the genotype/phenotype correlation in this dog family suggests that dogs with MPS‐I could have a less severe phenotype than humans, even in the presence of severe mutations. Treatment with pentosan polyphosphate should be considered in dogs with MPS‐I.

Multiple FGF4 Retrocopies Recently Derived within Canids

Two transcribed retrocopies of the fibroblast growth factor 4 (FGF4) gene have previously been described in the domestic dog. An FGF4 retrocopy on chr18 is associated with disproportionate dwarfism, while an FGF4 retrocopy on chr12 is associated with both disproportionate dwarfism and intervertebral disc disease (IVDD). In this study, whole-genome sequencing data were queried to identify other FGF4 retrocopies that could be contributing to phenotypic diversity in canids. Additionally, dogs with surgically confirmed IVDD were assayed for novel FGF4 retrocopies. Five additional and distinct FGF4 retrocopies were identified in canids including a copy unique to red wolves (Canis rufus). The FGF4 retrocopies identified in domestic dogs were identical to domestic dog FGF4 haplotypes, which are distinct from modern wolf FGF4 haplotypes, indicating that these retrotransposition events likely occurred after domestication. The identification of multiple, full length FGF4 retrocopies with open reading frames in canids indicates that gene retrotransposition events occur much more frequently than previously thought and provide a mechanism for continued genetic and phenotypic diversity in canids.

A nonsense variant in the KRT14 gene in a domestic shorthair cat with epidermolysis bullosa simplex

Epidermolysis bullosa simplex (EBS) is a hereditary blistering disease affecting the skin and mucous membranes. It has been reported in humans, cattle, buffaloes and dogs, but so far not in cats. In humans, EBS is most frequently caused by variants in the KRT5 or KRT14 genes. Here, we report a case of feline epidermolysis bullosa simplex and describe the causative genetic variant. An 11-month-old male domestic shorthair cat presented with a history of sloughed paw pads and ulcerations in the oral cavity and inner aspect of the pinnae, starting a few weeks after birth. Clinical and histopathological findings suggested a congenital blistering disease with a split formation within the basal cell layer of the epidermis and oral mucous epithelium. The genetic investigation revealed a homozygous nonsense variant in the KRT14 gene (c.979C>T, p.Gln327*). Immunohistochemistry showed a complete absence of keratin 14 staining in all epithelia present in the biopsy. To the best of our knowledge, this is the first report of feline EBS, and the first report of a spontaneous pathogenic KRT14 variant in a non-human species. The homozygous genotype in the affected cat suggests an autosomal recessive mode of inheritance.

Mitochondrial PCK2 Missense Variant in Shetland Sheepdogs with Paroxysmal Exercise-Induced Dyskinesia (PED)

Four female Shetland Sheepdogs with hypertonic paroxysmal dyskinesia, mainly triggered by exercise and stress, were investigated in a retrospective multi-center investigation aiming to characterize the clinical phenotype and its underlying molecular etiology. Three dogs were closely related and their pedigree suggested autosomal dominant inheritance. Laboratory diagnostic findings included mild lactic acidosis and lactaturia, mild intermittent serum creatine kinase (CK) elevation and hypoglycemia. Electrophysiological tests and magnetic resonance imaging of the brain were unremarkable. A muscle/nerve biopsy revealed a mild type II fiber predominant muscle atrophy. While treatment with phenobarbital, diazepam or levetiracetam did not alter the clinical course, treatment with a gluten-free, home-made fresh meat diet in three dogs or a tryptophan-rich, gluten-free, seafood-based diet, stress-reduction, and acetazolamide or zonisamide in the fourth dog correlated with a partial reduction in, or even a complete absence of, dystonic episodes. The genomes of two cases were sequenced and compared to 654 control genomes. The analysis revealed a case-specific missense variant, c.1658G>A or p.Arg553Gln, in the PCK2 gene encoding the mitochondrial phosphoenolpyruvate carboxykinase 2. Sanger sequencing confirmed that all four cases carried the mutant allele in a heterozygous state. The mutant allele was not found in 117 Shetland Sheepdog controls and more than 500 additionally genotyped dogs from various other breeds. The p.Arg553Gln substitution affects a highly conserved residue in close proximity to the GTP-binding site of PCK2. Taken together, we describe a new form of paroxysmal exercise-induced dyskinesia (PED) in dogs. The genetic findings suggest that PCK2:p.Arg553Gln should be further investigated as putative candidate causal variant.

A Third MLPH Variant Causing Coat Color Dilution in Dogs

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.

Pathogenic variants in COL6A3 cause Ullrich-like congenital muscular dystrophy in young Labrador Retriever dogs

The collagen VI-related muscular dystrophies in people include a broad spectrum of diseases ranging from the severe Ullrich congenital muscular dystrophy to the mild Bethlem myopathy. Clinical features are attributable to both muscle and connective tissue and include progressive muscle weakness and respiratory failure, hyperlaxity of distal joints, and progressive contracture of large joints. Here we describe two different COL6A3 pathogenic variants in Labrador Retriever dogs that result in autosomal recessive or autosomal dominant congenital myopathies with hyperlaxity of distal joints and joint contracture, similar to the condition in people.

ATP2A2 SINE Insertion in an Irish Terrier with Darier Disease and Associated Infundibular Cyst Formation

A 4-month-old female Irish Terrier presented with a well demarcated ulcerative and crusting lesion in the right ear canal. Histological analysis revealed epidermal hyperplasia with severe acantholysis affecting all suprabasal layers of the epidermis, which prompted a presumptive diagnosis of canine Darier disease. The lesion was successfully treated by repeated laser ablation of the affected epidermis. Over the course of three years, the dog additionally developed three dermal nodules of up to 4 cm in diameter that were excised and healed without complications. Histology of the excised tissue revealed multiple infundibular cysts extending from the upper dermis to the subcutis. The cysts were lined by squamous epithelium, which presented with abundant acantholysis of suprabasal keratinocytes. Infundibular cysts represent a novel finding not previously reported in Darier patients. Whole genome sequencing of the affected dog was performed, and the functional candidate genes for Darier disease (ATP2A2) and Hailey-Hailey disease (ATP2C1) were investigated. The analysis revealed a heterozygous SINE insertion into the ATP2A2 gene, at the end of intron 14, close to the boundary of exon 15. Analysis of the ATP2A2 mRNA from skin of the affected dog demonstrated a splicing defect and marked allelic imbalance, suggesting nonsense-mediated decay of the resulting aberrant transcripts. As Darier disease in humans is caused by haploinsufficiency of ATP2A2, our genetic findings are in agreement with the clinical and histopathological data and support the diagnosis of canine Darier disease.

A DSG1 Frameshift Variant in a Rottweiler Dog with Footpad Hyperkeratosis

A single male Rottweiler dog with severe footpad hyperkeratosis starting at an age of eight weeks was investigated. The hyperkeratosis was initially restricted to the footpads. The footpad lesions caused severe discomfort to the dog and had to be trimmed under anesthesia every 8–10 weeks. Histologically, the epidermis showed papillated villous projections of dense keratin in the stratum corneum. Starting at eight months of age, the patient additionally developed signs consistent with atopic dermatitis and recurrent bacterial skin and ear infections. Crusted hyperkeratotic plaques developed at sites of infection. We sequenced the genome of the affected dog and compared the data to 655 control genomes. A search for variants in 32 candidate genes associated with human palmoplantar keratoderma (PPK) revealed a single private protein-changing variant in the affected dog. This was located in the DSG1 gene encoding desmoglein 1. Heterozygous monoallelic DSG1 variants have been reported in human patients with striate palmoplantar keratoderma I (SPPK1), while biallelic DSG1 loss of function variants in humans lead to a more pronounced condition termed severe dermatitis, multiple allergies, and metabolic wasting (SAM) syndrome. The identified canine variant, DSG1:c.2541_2545delGGGCT, leads to a frameshift and truncates about 20% of the coding sequence. The affected dog was homozygous for the mutant allele. The comparative data on desmoglein 1 function in humans suggest that the identified DSG1 variant may have caused the footpad hyperkeratosis and predisposition for allergies and skin infections in the affected dog.

YARS2 Missense Variant in Belgian Shepherd Dogs with Cardiomyopathy and Juvenile Mortality

Dog puppy loss by the age of six to eight weeks after normal development is relatively uncommon. Necropsy findings in two spontaneously deceased Belgian Shepherd puppies indicated an abnormal accumulation of material in several organs. A third deceased puppy exhibited mild signs of an inflammation in the central nervous system and an enteritis. The puppies were closely related, raising the suspicion of a genetic cause. Pedigree analysis suggested a monogenic autosomal recessive inheritance. Combined linkage and homozygosity mapping assigned the most likely position of a potential genetic defect to 13 genome segments totaling 82 Mb. The genome of an affected puppy was sequenced and compared to 645 control genomes. Three private protein changing variants were found in the linked and homozygous regions. Targeted genotyping in 96 Belgian Shepherd dogs excluded two of these variants. The remaining variant, YARS2:1054G>A or p.Glu352Lys, was perfectly associated with the phenotype in a cohort of 474 Belgian Shepherd dogs. YARS2 encodes the mitochondrial tyrosyl-tRNA synthetase 2 and the predicted amino acid change replaces a negatively charged and evolutionary conserved glutamate at the surface of the tRNA binding domain of YARS2 with a positively charged lysine. Human patients with loss-of-function variants in YARS2 suffer from myopathy, lactic acidosis, and sideroblastic anemia 2, a disease with clinical similarities to the phenotype of the studied dogs. The carrier frequency was 27.2% in the tested Belgian Shepherd dogs. Our data suggest YARS2:1054G>A as the candidate causative variant for the observed juvenile mortality.

Whole Genome Sequencing Indicates Heterogeneity of Hyperostotic Disorders in Dogs

Craniomandibular osteopathy (CMO) and calvarial hyperostotic syndrome (CHS) are proliferative, non-neoplastic disorders affecting the skull bones in young dogs. Different forms of these hyperostotic disorders have been described in many dog breeds. However, an incompletely dominant causative variant for CMO affecting splicing of SLC37A2 has been reported so far only in three Terrier breeds. The purpose of this study was to identify further possible causative genetic variants associated with CHS in an American Staffordshire Terrier, as well as CMO in seven affected dogs of different breeds. We investigated their whole-genome sequences (WGS) and filtered variants using 584 unrelated genomes, which revealed no variants shared across all affected dogs. However, filtering for private variants of each case separately yielded plausible dominantly inherited candidate variants in three of the eight cases. In an Australian Terrier, a heterozygous missense variant in the COL1A1 gene (c.1786G>A; p.(Val596Ile)) was discovered. A pathogenic missense variant in COL1A1 was previously reported in humans with infantile cortical hyperostosis, or Caffey disease, resembling canine CMO. Furthermore, in a Basset Hound, a heterozygous most likely pathogenic splice site variant was found in SLC37A2 (c.1446+1G>A), predicted to lead to exon skipping as shown before in SLC37A2-associated canine CMO of Terriers. Lastly, in a Weimaraner, a heterozygous frameshift variant in SLC35D1 (c.1021_1024delTCAG; p.(Ser341ArgfsTer22)) might cause CMO due to the critical role of SLC35D1 in chondrogenesis and skeletal development. Our study indicates allelic and locus heterogeneity for canine CMO and illustrates the current possibilities and limitations of WGS-based precision medicine in dogs.

A Missense Variant Affecting the C-Terminal Tail of UNC93B1 in Dogs with Exfoliative Cutaneous Lupus Erythematosus (ECLE)

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.

Characterisation of canine KCNIP4: A novel gene for cerebellar ataxia identified by whole-genome sequencing two affected Norwegian Buhund dogs

A form of hereditary cerebellar ataxia has recently been described in the Norwegian Buhund dog breed. This study aimed to identify the genetic cause of the disease. Whole-genome sequencing of two Norwegian Buhund siblings diagnosed with progressive cerebellar ataxia was carried out, and sequences compared with 405 whole genome sequences of dogs of other breeds to filter benign common variants. Nine variants predicted to be deleterious segregated among the genomes in concordance with an autosomal recessive mode of inheritance, only one of which segregated within the breed when genotyped in additional Norwegian Buhunds. In total this variant was assessed in 802 whole genome sequences, and genotyped in an additional 505 unaffected dogs (including 146 Buhunds), and only four affected Norwegian Buhunds were homozygous for the variant. The variant identified, a T to C single nucleotide polymorphism (SNP) (NC_006585.3:g.88890674T>C), is predicted to cause a tryptophan to arginine substitution in a highly conserved region of the potassium voltage-gated channel interacting protein KCNIP4. This gene has not been implicated previously in hereditary ataxia in any species. Evaluation of KCNIP4 protein expression through western blot and immunohistochemical analysis using cerebellum tissue of affected and control dogs demonstrated that the mutation causes a dramatic reduction of KCNIP4 protein expression. The expression of alternative KCNIP4 transcripts within the canine cerebellum, and regional differences in KCNIP4 protein expression, were characterised through RT-PCR and immunohistochemistry respectively. The voltage-gated potassium channel protein KCND3 has previously been implicated in spinocerebellar ataxia, and our findings suggest that the Kv4 channel complex KCNIP accessory subunits also have an essential role in voltage-gated potassium channel function in the cerebellum and should be investigated as potential candidate genes for cerebellar ataxia in future studies in other species.

A deletion spanning the promoter and first exon of the hair cycle-specific ASIP transcript isoform in black and tan rabbits

Black and tan animals have tan-coloured ventral body surfaces separated by sharp boundaries from black-coloured dorsal body surfaces. In the a^t mouse mutant, a retroviral 6 kb insertion located in the hair cycle-specific promoter of the murine Asip gene encoding agouti signalling protein causes the black and tan phenotype. In rabbits, three ASIP alleles are thought to exist, including an a^t allele causing a black and tan coat colour that closely resembles the mouse black and tan phenotype. The goal of our study was to identify the functional genetic variant causing the rabbit a^t allele. We performed a WGS-based comparative analysis of the ASIP gene in one black and tan and three wt agouti-coloured rabbits. The analysis identified 75 a^t -associated variants including an 11 kb deletion. The deletion is located in the region of the hair cycle-specific ASIP promoter and thus in a region homologous to the site of the retroviral insertion causing the a^t allele in mice. We observed perfect association of the genotypes at this deletion with the coat colour phenotype in 49 rabbits. The comparative analysis and the previous knowledge about the regulation of ASIP expression suggest that the 11 kb deletion is the most likely causative variant for the black and tan phenotype in rabbits.