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

Intragenic PNPLA1 duplication in Labrador retrievers with nonepidermolytic ichthyosis

BACKGROUND

Ichthyoses represent a heterogeneous group of cornification disorders characterised by epidermal scaling.

OBJECTIVES

To describe the clinical, histopathological and genetic analysis of a Labrador retriever with nonepidermolytic ichthyosis, and the results of a population screening for a newly detected PNPLA1 genomic duplication.

ANIMALS

Two 7-year-old male littermates, 531 population samples.

MATERIALS AND METHODS

Clinical and histopathological analysis, whole genome sequencing and digital PCR-based genotyping were performed.

RESULTS

Generalised scaling and histological laminar orthokeratotic hyperkeratosis confirmed the ichthyosis diagnosis on Dog 1. Dog 2 showed mild clinical signs possibly associated with allergies and not ichthyosis. The genome of Dog 1 was sequenced and compared to 1469 genetically diverse control genomes. The analysis identified a 6099-bp duplication spanning three internal exons of the PNPLA1 gene, which is predicted to result in an altered C-terminal tail of the protein, NP_001277038.2:p.(E558Lfs*17). Dog 2 had a heterozygous genotype and carried one copy of the duplicated PNPLA1 allele. Of the screened 531 additional Labrador retrievers, 491 were homozygous wild-type, 36 were heterozygous carriers and four carried the duplication in a homozygous state.

CONCLUSIONS AND CLINICAL RELEVANCE

Previously identified PNPLA1 variants cause autosomal recessive ichthyosis in golden retrievers and humans. Given the well-established function of PNPLA1, the identified genomic duplication represents a likely candidate causal variant for the observed ichthyosis in the examined Labrador retriever. This is the first report of a new form of autosomal recessive ichthyosis in Labrador retrievers, which provides the basis for genetic testing.

Epileptic encephalopathy in a young Bengal cat caused by CAD deficiency

Developmental and epileptic encephalopathy type 50 (DEE50) in humans is a severe early-onset neurometabolic disorder caused by biallelic loss-of-function variants in the CAD gene encoding a key multi-enzymatic protein for de novo pyrimidine nucleotide synthesis. Untreated, the condition is often fatal, but patients respond to uridine supplementation, which fuels nucleotide synthesis through CAD-independent salvage pathways. Here, we report a novel variant in the feline CAD gene in a 4-month-old Bengal kitten with intractable seizures and abnormal behavior. The variant, XP_011279586.1:p.(Ser2015Asn), was predicted to affect the oligomerization of the C-terminal aspartate transcarbamylase (ATCase) domain of CAD. Genotyping of 110 unaffected Bengal cats revealed four additional carriers of the mutant allele, confirming its presence in the breed. In a CAD-knockout human cell line dependent on uridine, the recombinant expression of human wildtype CAD, but not of the Asn2015 mutant, restored cell growth without uridine, demonstrating that the p.Ser2015Asn variant disrupts CAD function and is pathogenic. This study facilitates genetic testing of carriers and affected cats and suggests that uridine supplementation could be a potential treatment. Furthermore, CAD-deficient Bengal cats might serve as a valuable spontaneous large animal model to further investigate the pathogenic mechanisms of this rare epileptic encephalopathy in humans.

Whole genome sequencing identifies novel candidate genetic variants in canine stomatocytosis

Stomatocytosis is a rare spectrum of red blood cell (RBC) disorders. In humans, stomatocytosis is typically caused by genetic changes in specific ion exchange and transport genes. Stomatocytosis has been identified in dogs, however the underlying genetic causes are unknown. Recently, stomatocytosis was reported in a Beagle and Australian Cattle Dog for the first time. Here, whole-genome sequencing (WGS) of these dogs was undertaken to identify candidate genetic variants driving or impacting stomatocytosis. Cases were compared to WGS of 119 controls of several breeds and > 1,000 dogs from public and private datasets. Candidate genes were identified, including genes linked to stomatocytosis in humans: SPTB and KCNN4. Notably, each case carried a different homozygous intronic SNP in SPTB only 24 bases apart (Beagle - chr8:39,194,923; ACD - chr8:39,194,947; CanFam3.1), which were not homozygous in other dogs. Variants with predicted deleterious impact in additional ion transport-related genes were also identified: SLC8A3, DYSF, SLC12A8, INPP5E, SLC1A1, and a novel SLC41A3 genetic change carried by the Australian Cattle Dog. Human and mouse scRNAseq and proteomics data indicate that these candidate genes are expressed in RBCs or their immature precursors. Taken together, these genetic data obtained from spontaneous stomatocytosis in a non-human species provide novel insights and candidate genes for evaluation of rare red cell disorders in humans.

Measures of Homozygosity and Relationship to Genetic Diversity in the Bearded Collie Breed

Background: Genetic diversity in closed populations, such as pedigree dogs, is of concern for maintaining the health and vitality of the population in the face of evolving challenges. Measures of genetic diversity rely upon estimates of homozygosity without consideration of whether the homozygosity is desirable or undesirable or if heterozygosity has a functional impact. Pedigree coefficients of inbreeding have been the classical approach yet they are inadequate unless based upon the entire population. Methods: Homozygosity measures based upon pedigree analyses (n = 11,898), SNP array data (n = 244), and whole genome sequencing (n = 23) were compared in the Bearded Collie, as well as a comparison of SNP array data to a pedigree cohort (n = 5042) and a mixed-breed cohort (n = 1171). Results: Molecular measures based upon DNA are more informative on an individual's homozygosity levels than pedigree analyses, although SNP coefficients of inbreeding overestimate the level of inbreeding based on the nature of SNP array methodology. Whole genome sequence (WGS) analyses revealed that the heterozygosity observed is generally in variants having neutral or low impact, which would indicate that the variability may not contribute substantially to functional diversity in the population. The majority of high-impact variants were observed in the shortest runs of homozygosity (ROH) reflecting ancestral breeding and domestication practices. As expected, mixed-breed dogs displayed higher measures of genomic diversity than either Bearded Collies or other pedigree dogs as a whole using the current paradigm algorithm models to calculate homozygosity. Conclusions: Using typical DNA-based measures reflect only a single individual and not the population thereby failing to account for regions of homozygosity that reflect ancestral breeding, domestication history, breed-defining regions, or regions positively selected for health traits. Incorporating measures of genetic diversity into dog breeding schemes is meritorious. However, until measures of diversity can distinguish between breed-defining homozygosity and homozygosity associated with positive health alleles, the measures to use as selection tools need refinement before their widespread implementation.

The RSPO2 gene is associated with bilateral anterior amelia in Chihuahuas

Bilateral anterior amelia (BAA) is the congenital absence of thoracic limbs and has been reported in the Chihuahua as an autosomal recessive disorder. In some cases, the digits of the pelvic limbs can be variably affected, but otherwise, the pelvic limbs are generally spared. A GWAS performed with nine BAA affected Chihuahuas identified a significant association on chromosome 13, and homozygosity mapping delineated a 2.1 Mb chromosomal region containing the RSPO2 gene. Loss of function variants of RSPO2 in humans and cattle has been associated with the absence of all limbs. Six affected Chihuahuas were whole genome sequenced (WGS) and aligned to the CanFam4 assembly. SNVs, small indels, and structural variants within the critical interval that fitted a recessive model were investigated. Three SNVs (NC_049234.1:g.8891861C > T; NC_049234.1:g.8974204C > T and NC_049234.1:g.9789424G > A) were homozygous in five cases and absent from 3,418 genetically diverse control genome sequences, except for one Small Poodle that was heterozygous. One SNV resided in RSPO2's second intron, while the two others were intergenic. The three candidate variants were genotyped in 7 additional cases and 100 control Chihuahuas. Twelve of 13 cases were homozygous for the mutant allele, and one case was heterozygous. Controls were either homozygous for the reference allele (97%) or heterozygous (3%). Our data should facilitate genetic testing of Chihuahuas to prevent the unintentional production of BAA affected dogs. Moreover, the identification of these variants enhances understanding of RSPO2 gene function in limb development.

GTPBP2 in-frame deletion in canine model with non-syndromic progressive retinal atrophy

Progressive retinal atrophy (PRA), caused by aberrant functioning of rod/cone photoreceptors, leads to blindness affecting mammals, including dogs. We identified a litter of three Labrador retrievers affected by non-syndromic PRA; the parents and three other siblings were unaffected. Homozygosity mapping and whole-genome sequencing detected a homozygous 3-bp deletion in the coding region of GTPBP2, located in CFA12 (NC_049233.1:12,264,348_12,264,350del, c.1606_1608del, p.Ala536del). The variant was absent from the online European Variation Archive (EVA) database, the Dog Biomedical Variants Database Consortium, and the Dog10k database. We tested 91 non-affected dogs from the same kennel and found 75 wild-type (WT) and 16 carriers, all clinically normal, and 569 Labradors from the general population (USA), all WT. GTPBP2 is associated with Jaberi-Elahi syndrome (JES) in Homo sapiens, and splice variants in Mus musculus are associated with neurodegeneration; in both cases photoreceptor degeneration may be included in its manifestation. Heterologous cellular systems were transfected with cDNA encoding WT or A536del mutant GTPBP2 protein and immunoblot analysis of total cell lysate with anti-GTPBP2 antibodies showed that the expression level of the GTPBP2 mutant protein A536del is slightly but not significantly reduced compared to WT. Immunofluorescent methods and confocal analysis of cells transfected with WT or A536del GTPBP2 protein revealed that the WT form is diffuse throughout the cytosol, while the mutant form resulted in the formation of cytoplasmic aggregates in ~70-80% of cells. The deleted amino acid falls within a conserved interval outside the GTP domain of GTPBP2, suggesting a potentially novel role of the sequence on cellular localization of the protein.

Novel ryanodine receptor 1 (RYR1) missense gene variants in two pet dogs with fatal malignant hyperthermia identified by next-generation sequencing

OBJECTIVE

Evaluate a precision medicine approach to confirm a tentative diagnosis of fatal malignant hyperthermia (MH) in isoflurane-anesthetized pet dogs by identifying novel risk variants in known MH susceptibility genes.

STUDY DESIGN

Retrospective case series.

ANIMALS

A male Pit Bull mix aged 7 years (case #1), a male Golden Retriever aged 12 months (case #2) and the dam and sire of case #2.

METHODS

Available case histories and medical records were reviewed. Missense variants in MH susceptibility genes RYR-1, CACNA1S and STAC3 (case #2 only) were identified by next-generation sequencing of DNA from each case and the parents of case #2 with confirmation by Sanger sequencing. The pathogenicity of variants was evaluated by multiple in silico approaches.

RESULTS

Both cases demonstrated clinical signs during isoflurane anesthesia consistent with volatile anesthetic-induced MH, including tachypnea, tachycardia, severe hyperthermia and muscle rigidity. Despite whole body cooling and other treatments, both dogs died after cardiac arrest within 15 minutes of detecting hyperthermia. Gene sequencing identified novel missense RYR-1 variants in case #1 (p.Gly2375Arg) and case #2 (p.Pro152Leu). Both variants were likely pathogenic based on multiple criteria, including gene location, amino acid alteration and population allele frequency. The case #1 variant was identical to a known human diagnostic MH variant (p.Gly2375Arg). Neither parent of case #2 had the case #2 variant, indicating this variant was not inherited, but arose de novo in a germ cell of either parent or early in embryogenesis. Whole genome sequence analysis confirmed parentage. Two missense variants were identified in CACNA1S. Both variants were considered nonpathogenic. No variants were identified in STAC3.

CONCLUSIONS AND CLINICAL RELEVANCE

Like humans, MH susceptibility in dogs is associated with different rare variants located in pathogenic hotspots in the RYR-1 gene. Next-generation sequencing is a useful tool to assist in the definitive diagnosis of MH in dogs.

Intragenic duplication disrupting the reading frame of MFSD8 in Small Swiss Hounds with neuronal ceroid lipofuscinosis

Neuronal ceroid lipofuscinosis (NCL) represents a heterogenous group of lysosomal storage diseases resulting in progressive neurodegeneration. We investigated two Small Swiss Hound littermates that showed progressive ataxia and loss of cognitive functions and vision starting around the age of 12 months. Both dogs had to be euthanized a few months after the onset of disease owing to the severity of their clinical signs. Pathological investigation of one affected dog revealed cerebral and cerebellar atrophy with cytoplasmic accumulation of autofluorescent material in degenerating neurons. The clinical signs in combination with the characteristic histopathology led to a tentative diagnosis of NCL. In the subsequent genetic investigation, the genome of one affected dog was sequenced. This revealed a duplication of 18 819 bp within the MFSD8 gene. The duplication breakpoints were located in intron 3 and exon 12 of the gene and were predicted to disrupt the reading frame. Both affected dogs carried the duplication in a homozygous state and there was perfect cosegregation of the genotypes with the phenotype in a large pedigree, consistent with autosomal recessive inheritance. MFSD8 loss-of-function variants are a known cause of NCL7 in human patients, dogs and other mammalian species. The existing knowledge on MFSD8 together with the experimental data strongly suggests that the identified intragenic MFSD8 duplication caused the disease in the Small Swiss Hounds. These results allow their diagnosis to be refined to NCL7 and enable genetic testing in the breed to avoid further unintentional carrier × carrier matings.

Intragenic dystrophin (DMD) duplication variant in Entlebucher Mountain Dogs with Duchenne muscular dystrophy

Muscular dystrophies represent a group of disorders characterized by progressive muscle degeneration and weakness. An important subgroup are the dystrophin-related muscular dystrophies caused by variants in the DMD gene. They can be divided into the more severe Duchenne muscular dystrophy and the milder Becker muscular dystrophy. Here, we characterize the clinical, histopathological and molecular genetic aspects of two male Entlebucher Mountain Dogs with clinical signs of muscular dystrophy. The two dogs presented with marked dysphagia starting at the age of several weeks and in the later course recognizable exercise intolerance with highly increased serum creatine kinase levels. Histopathological signs of a dystrophic myopathy represented by degeneration of muscle fibers and signs of regeneration were present. Whole genome sequencing of one affected dog identified an intragenic 8.6 kb duplication in the X-chromosomal DMD gene, c.7528-4048_7645 + 4450dup. No other protein-changing variants in candidate genes for muscular dystrophy were identified. The duplication includes exon 52 of DMD and is predicted to lead to a frameshift and truncation of 30% of the wild-type open reading frame. Genotyping of the whole family confirmed the presence of the mutant allele in both affected dogs and the unaffected dam. The correct co-segregation of the mutant allele in the affected family as well as knowledge from humans and other species suggest the identified DMD variant as the most likely candidate variant for the muscular dystrophy phenotype in the two investigated dogs.

A map of canine sequence variation relative to a Greenland wolf outgroup

For over 15 years, canine genetics research relied on a reference assembly from a Boxer breed dog named Tasha (i.e., canFam3.1). Recent advances in long-read sequencing and genome assembly have led to the development of numerous high-quality assemblies from diverse canines. These assemblies represent notable improvements in completeness, contiguity, and the representation of gene promoters and gene models. Although genome graph and pan-genome approaches have promise, most genetic analyses in canines rely upon the mapping of Illumina sequencing reads to a single reference. The Dog10K consortium, and others, have generated deep catalogs of genetic variation through an alignment of Illumina sequencing reads to a reference genome obtained from a German Shepherd Dog named Mischka (i.e., canFam4, UU_Cfam_GSD_1.0). However, alignment to a breed-derived genome may introduce bias in genotype calling across samples. Since the use of an outgroup reference genome may remove this effect, we have reprocessed 1929 samples analyzed by the Dog10K consortium using a Greenland wolf (mCanLor1.2) as the reference. We efficiently performed remapping and variant calling using a GPU-implementation of common analysis tools. The resulting call set removes the variability in genetic differences seen across samples and breed relationships revealed by principal component analysis are not affected by the choice of reference genome. Using this sequence data, we inferred the history of population sizes and found that village dog populations experienced a 9-13 fold reduction in historic effective population size relative to wolves.

A RETREG1 variant is associated with hereditary sensory and autonomic neuropathy with acral self-mutilation in purebred German Spitz

Hereditary sensory and autonomic neuropathies (HSAN) represent a group of genetic diseases affecting the peripheral nervous system. In humans, at least 16 loci have been associated with the disorder but do not explain the disease origin of all patients. In dogs, similar conditions have been documented for decades in various breeds with a severe impact on life quality and are often referred to as acral mutilation syndrome (AMS). Causal variants in three genes have been identified to date, suggesting larger genetic heterogeneity in the dog population. Our aim was to explain the genetic etiology of an early-onset HSAN/AMS in a purebred German Spitz. The affected dog showed progressive loss of pain sensation in the distal extremities, which led to intense licking, biting, and self-mutilation of digits and paw pads. Whole-genome sequencing identified a single candidate causal variant on chromosome 4 in the RETREG1 gene (c.656C>T, p.Pro219Leu). This missense variant was previously recognized as deleterious in a mixed breed dog family with similar clinical signs. Haplotype analyses and targeted genotyping revealed a likely German Spitz ancestry of these mixed breed dogs. Further screening of an extensive cohort of ~900 000 dogs of various breeds hinted at the variant allele origin in the German Spitz breed. Disruption of RETREG1 inhibits endoplasmic reticulum turnover and leads to neuron degeneration. Our findings provide evidence that this variant underlies the recessive form of HSAN/AMS in the German Spitz and support the use of whole-genome sequencing-based veterinary precision medicine for early diagnosis and prevention via a genetic test.

Canine RNF170 Single Base Deletion in a Naturally Occurring Model for Human Neuroaxonal Dystrophy

BACKGROUND

Neuroaxonal dystrophy (NAD) is a group of inherited neurodegenerative disorders characterized primarily by the presence of spheroids (swollen axons) throughout the central nervous system. In humans, NAD is heterogeneous, both clinically and genetically. NAD has also been described to naturally occur in large animal models, such as dogs. A newly recognized disorder in Miniature American Shepherd dogs (MAS), consisting of a slowly progressive neurodegenerative syndrome, was diagnosed as NAD via histopathology.

OBJECTIVES

To describe the clinical and pathological phenotype together with the identification of the underlying genetic cause.

METHODS

Clinical and postmortem evaluations, together with a genome-wide association study and autozygosity mapping approach, followed by whole-genome sequencing.

RESULTS

Affected dogs were typically young adults and displayed an abnormal gait characterized by pelvic limb weakness and ataxia. The underlying genetic cause was identified as a 1-bp (base pair) deletion in RNF170 encoding ring finger protein 170, which perfectly segregates in an autosomal recessive pattern. This deletion is predicted to create a frameshift (XM_038559916.1:c.367delG) and early truncation of the RNF170 protein (XP_038415844.1:(p.Ala123Glnfs*11)). The age of this canine RNF170 variant was estimated at ~30 years, before the reproductive isolation of the MAS breed.

CONCLUSIONS

RNF170 variants were previously identified in human patients with autosomal recessive spastic paraplegia-85 (SPG85); this clinical phenotype shows similarities to the dogs described herein. We therefore propose that this novel MAS NAD could serve as an excellent large animal model for equivalent human diseases, particularly since affected dogs demonstrate a relatively long lifespan, which represents an opportunity for therapeutic trials. © 2024 The Author(s). Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

Heterozygous DSP in-frame deletion in a poodle with syndromic ichthyosis involving additional hair and tooth abnormalities

Ichthyoses comprise a large heterogeneous group of skin disorders, characterized by generalized scaly and hyperkeratotic skin. We investigated a miniature poodle with early onset generalized scaling, dry and irregularly thickened skin, paw pad hyperkeratosis and abnormalities in hair and teeth. The clinical signs of ichthyosis were confirmed by histopathological examination, which revealed mild epidermal hyperplasia and lamellar orthokeratotic hyperkeratosis. A hereditary condition was suspected and a genetic investigation was initiated. We sequenced the whole genome of the affected dog and searched for potentially causative variants in functional candidate genes for the observed phenotype. The analysis revealed a heterozygous in-frame deletion in DSP, NC_049256.1:g.8804542_8804544del resulting from a de novo mutation event as evidenced by genotyping leukocyte DNA from both parents. The 3 bp deletion is predicted to remove one aspartic acid without disrupting the open reading frame (XM_038584124.1:c.1821_1823del, XP_038440052.1:p.(Asp608del)). The DSP gene encodes desmoplakin, a desmosomal plaque protein, responsible for cell-cell adhesion to provide resistance to mechanical stress in epidermal and cardiac tissues. We hypothesize that the deletion of one amino acid in the N-terminal globular head domain acts in a dominant negative manner and thus impairs the proper connection with other proteins. Several variants in DSP in humans and cattle have been described to result in different phenotypes associated with hair and skin abnormalities, sometimes in combination with variable cardiac and/or dental manifestations. In conclusion, we characterized a new syndromic ichthyosis phenotype in a dog and identified a de novo 3 bp deletion in the DSP gene as causal variant.

A Variant-Centric Analysis of Allele Sharing in Dogs and Wolves

Canines are an important model system for genetics and evolution. Recent advances in sequencing technologies have enabled the creation of large databases of genetic variation in canines, but analyses of allele sharing among canine groups have been limited. We applied GeoVar, an approach originally developed to study the sharing of single nucleotide polymorphisms across human populations, to assess the sharing of genetic variation among groups of wolves, village dogs, and breed dogs. Our analysis shows that wolves differ from each other at an average of approximately 2.3 million sites while dogs from the same breed differ at nearly 1 million sites. We found that 22% of the variants are common across wolves, village dogs, and breed dogs, that ~16% of variable sites are common across breed dogs, and that nearly half of the differences between two dogs of different breeds are due to sites that are common in all clades. These analyses represent a succinct summary of allele sharing across canines and illustrate the effects of canine history on the apportionment of genetic variation.

Identification of nonsense variants in the genomes of 15 Murciano-Granadina bucks and analysis of their segregation in parent-offspring trios

Nonsense variants can inactivate gene function by causing the synthesis of truncated proteins or by inducing nonsense mediated decay of messenger RNAs. The occurrence of such variants in the genomes of livestock species is modulated by multiple demographic and selective factors. Even though nonsense variants can have causal effects on embryo lethality, abortions, and disease, their genomic distribution and segregation in domestic goats have not been characterized in depth yet. In this work, we have sequenced the genomes of 15 Murciano-Granadina bucks with an average coverage of 32.92 × ± 1.45 × . Bioinformatic analysis revealed 947 nonsense variants consistently detected with SnpEff and Ensembl-VEP. These variants were especially abundant in the 3'end of the protein-coding regions. Genes related to olfactory perception, ATPase activity coupled to transmembrane movement of substances, defense to virus, hormonal response, and sensory perception of taste were particularly enriched in nonsense variants. Seventeen nonsense variants expected to have harmful effects on fitness were genotyped in parent-offspring trios. We observed that several nonsense variants predicted to be lethal based on mouse knockout data did not have such effect, a finding that could be explained by the existence of multiple mechanisms counteracting lethality. These findings demonstrate that predicting the effects of putative nonsense variants on fitness is extremely challenging. As a matter of fact, such a goal could only be achieved by generating a high quality telomere-to-telomere goat reference genome combined with carefully curated annotation and functional testing of promising candidate variants.

Neuronal ceroid lipofuscinosis in a Schapendoes dog is caused by a missense variant in CLN6

Neuronal ceroid lipofuscinosis (NCL) is a group of neurodegenerative disorders that occur in humans, dogs, and several other species. NCL is characterised clinically by progressive deterioration of cognitive and motor function, epileptic seizures, and visual impairment. Most forms present early in life and eventually lead to premature death. Typical pathological changes include neuronal accumulation of autofluorescent, periodic acid-Schiff- and Sudan black B-positive lipopigments, as well as marked loss of neurons in the central nervous system. Here, we describe a 19-month-old Schapendoes dog, where clinical signs were indicative of lysosomal storage disease, which was corroborated by pathological findings consistent with NCL. Whole genome sequencing of the affected dog and both parents, followed by variant calling and visual inspection of known NCL genes, identified a missense variant in CLN6 (c.386T>C). The variant is located in a highly conserved region of the gene and predicted to be harmful, which supports a causal relationship. The identification of this novel CLN6 variant enables pre-breeding DNA-testing to prevent future cases of NCL6 in the Schapendoes breed, and presents a potential natural model for NCL6 in humans.

KRT5 in-frame deletion in a family of German Shepherd dogs with split paw pad disease resembling localized epidermolysis bullosa simplex in human patients

Split paw pad disease is a scarcely defined phenotype characterized by skin lesions on the paw pads of dogs. We studied a family of German Shepherd dogs, in which four dogs developed intermittent paw pad lesions and lameness. The paw pads of two of the affected dogs were biopsied and demonstrated cleft formation in the stratum spinosum and stratum corneum, the outermost layers of the epidermis. Whole genome sequencing data from an affected dog revealed a private heterozygous 18 bp in frame deletion in the KRT5 gene. The deletion NM_001346035.1:c.988_1005del or NP_001332964.1:p.(Asn330_Asp335del) is predicted to lead to a loss of six amino acids in the L12 linker domain of the encoded keratin 5. KRT5 variants in human patients lead to various subtypes of epidermolysis bullosa simplex (EBS). Localized EBS is the mildest of the KRT5-related human diseases and may be caused by variants affecting the L12 linker domain of keratin 5. We therefore think that the detected KRT5 deletion in dogs represents a candidate causal variant for the observed skin lesions in dogs. However, while the clinical phenotype of KRT5-mutant dogs of this study closely resembles human patients with localized EBS, there are differences in the histopathology. EBS is defined by cleft formation within the basal layer of the epidermis while the cleft formation in the dogs described herein occurred in the outermost layers, a hallmark of split paw pad disease. Our study provides a basis for further studies into the exact relation of split paw pad disease and EBS.

Intragenic MFSD8 duplication and histopathological findings in a rabbit with neuronal ceroid lipofuscinosis

Neuronal ceroid lipofuscinoses (NCL) are among the most prevalent neurodegenerative disorders of early life in humans. Disease-causing variants have been described for 13 different NCL genes. In this study, a refined pathological characterization of a female rabbit with progressive neurological signs reminiscent of NCL was performed. Cytoplasmic pigment present in neurons was weakly positive with Sudan black B and autofluorescent. Immunohistology revealed astrogliosis, microgliosis and axonal degeneration. During the subsequent genetic investigation, the genome of the affected rabbit was sequenced and examined for private variants in NCL candidate genes. The analysis revealed a homozygous ~10.7 kb genomic duplication on chromosome 15 comprising parts of the MFSD8 gene, NC_013683.1:g.103,727,963_103,738,667dup. The duplication harbors two internal protein coding exons and is predicted to introduce a premature stop codon into the transcript, truncating ~50% of the wild-type MFSD8 open reading frame encoding the major facilitator superfamily domain containing protein 8, XP_002717309.2:p.(Glu235Leufs*23). Biallelic loss-of-function variants in MFSD8 have been described to cause NCL7 in human patients, dogs and a single cat. The available clinical and pathological data, together with current knowledge about MFSD8 variants and their functional impact in other species, point to the MFSD8 duplication as a likely causative defect for the observed phenotype in the affected rabbit.

The Association between the Abundance of Homozygous Deleterious Variants and the Morbidity of Dog Breeds

It is well known that highly inbred dogs are more prone to diseases than less inbred or outbred dogs. This is because inbreeding increases the load of recessive deleterious variants. Using the genomes of 392 dogs belonging to 83 breeds, we investigated the association between the abundance of homozygous deleterious variants and dog health. We used the number of non-routine veterinary care events for each breed to assess the level of morbidity. Our results revealed a highly significant positive relationship between the number of homozygous deleterious variants located within the runs of homozygosity (RoH) tracts of the breeds and the level of morbidity. The dog breeds with low morbidity had a mean of 87 deleterious SNVs within the RoH, but those with very high morbidity had 187 SNVs. A highly significant correlation was also observed for the loss-of-function (LoF) SNVs within RoH tracts. The dog breeds that required more veterinary care had 2.3 times more homozygous LoF SNVs than those that required less veterinary care (112 vs. 50). The results of this study could be useful for understanding the disease burden on breed dogs and as a guide for dog breeding programs.

EDA Missense Variant in a Cat with X-Linked Hypohidrotic Ectodermal Dysplasia

Hypohidrotic ectodermal dysplasia is a developmental defect characterized by sparse or absent hair, missing or malformed teeth and defects in eccrine glands. Loss-of-function variants in the X-chromosomal EDA gene have been reported to cause hypohidrotic ectodermal dysplasia in humans, mice, dogs and cattle. We investigated a male cat exhibiting diffuse truncal alopecia with a completely absent undercoat. The cat lacked several teeth, and the remaining teeth had an abnormal conical shape. Whole-genome sequencing revealed a hemizygous missense variant in the EDA gene, XM_011291781.3:c.1042G>A or XP_011290083.1:p.(Ala348Thr). The predicted amino acid exchange is located in the C-terminal TNF signaling domain of the encoded ectodysplasin. The corresponding missense variant in the human EDA gene, p.Ala349Thr, has been reported as a recurring pathogenic variant in several human patients with X-linked hypohidrotic ectodermal dysplasia. The identified feline variant therefore represents the likely cause of the hypohidrotic ectodermal dysplasia in the investigated cat, and the genetic investigation confirmed the suspected clinical diagnosis. This is the first report of an EDA-related hypohidrotic ectodermal dysplasia in cats.

A naturally occurring canine model of syndromic congenital microphthalmia

In humans, the prevalence of congenital microphthalmia is estimated to be 0.2-3.0 for every 10,000 individuals, with nonocular involvement reported in ∼80% of cases. Inherited eye diseases have been widely and descriptively characterized in dogs, and canine models of ocular diseases have played an essential role in unraveling the pathophysiology and development of new therapies. A naturally occurring canine model of a syndromic disorder characterized by microphthalmia was discovered in the Portuguese water dog. As nonocular findings included tooth enamel malformations, stunted growth, anemia, and thrombocytopenia, we hence termed this disorder Canine Congenital Microphthalmos with Hematopoietic Defects. Genome-wide association study and homozygosity mapping detected a 2 Mb candidate region on canine chromosome 4. Whole-genome sequencing and mapping against the Canfam4 reference revealed a Short interspersed element insertion in exon 2 of the DNAJC1 gene (g.74,274,883ins[T70]TGCTGCTTGGATT). Subsequent real-time PCR-based mass genotyping of a larger Portuguese water dog population found that the homozygous mutant genotype was perfectly associated with the Canine Congenital Microphthalmos with Hematopoietic Defects phenotype. Biallelic variants in DNAJC21 are mostly found to be associated with bone marrow failure syndrome type 3, with a phenotype that has a certain degree of overlap with Fanconi anemia, dyskeratosis congenita, Shwachman-Diamond syndrome, Diamond-Blackfan anemia, and reports of individuals showing thrombocytopenia, microdontia, and microphthalmia. We, therefore, propose Canine Congenital Microphthalmos with Hematopoietic Defects as a naturally occurring model for DNAJC21-associated syndromes.

Increased susceptibility to Mycobacterium avium complex infection in miniature Schnauzer dogs caused by a codon deletion in CARD9

Mammals are generally resistant to Mycobacterium avium complex (MAC) infections. We report here on a primary immunodeficiency disorder causing increased susceptibility to MAC infections in a canine breed. Adult Miniature Schnauzers developing progressive systemic MAC infections were related to a common founder, and pedigree analysis was consistent with an autosomal recessive trait. A genome-wide association study and homozygosity mapping using 8 infected, 9 non-infected relatives, and 160 control Miniature Schnauzers detected an associated region on chromosome 9. Whole genome sequencing of 2 MAC-infected dogs identified a codon deletion in the CARD9 gene (c.493_495del; p.Lys165del). Genotyping of Miniature Schnauzers revealed the presence of this mutant CARD9 allele worldwide, and all tested MAC-infected dogs were homozygous mutants. Peripheral blood mononuclear cells from a dog homozygous for the CARD9 variant exhibited a dysfunctional CARD9 protein with impaired TNF-α production upon stimulation with the fungal polysaccharide β-glucan that activates the CARD9-coupled C-type lectin receptor, Dectin-1. While CARD9-deficient knockout mice are susceptible to experimental challenges by fungi and mycobacteria, Miniature Schnauzer dogs with systemic MAC susceptibility represent the first spontaneous animal model of CARD9 deficiency, which will help to further elucidate host defense mechanisms against mycobacteria and fungi and assess potential therapies for animals and humans.

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.

NDUFS7 variant in dogs with Leigh syndrome and its functional validation in a Drosophila melanogaster model

Two Jack-Russell Terrier × Chihuahua mixed-breed littermates with Leigh syndrome were investigated. The dogs presented with progressive ataxia, dystonia, and increased lactate levels. Brain MRI showed characteristic bilateral symmetrical T2 hyperintense lesions, histologically representing encephalomalacia. Muscle histopathology revealed accumulation of mitochondria. Whole genome sequencing identified a missense variant in a gene associated with human Leigh syndrome, NDUFS7:c.535G > A or p.(Val179Met). The genotypes at the variant co-segregated with the phenotype in the investigated litter as expected for a monogenic autosomal recessive mode of inheritance. We investigated the functional consequences of the missense variant in a Drosophila melanogaster model by expressing recombinant wildtype or mutant canine NDUFS7 in a ubiquitous knockdown model of the fly ortholog ND-20. Neither of the investigated overexpression lines completely rescued the lethality upon knockdown of the endogenous ND-20. However, a partial rescue was found upon overexpression of wildtype NDUFS7, where pupal lethality was moved to later developmental stages, which was not seen upon canine mutant overexpression, thus providing additional evidence for the pathogenicity of the identified variant. Our results show the potential of the fruit fly as a model for canine disease allele validation and establish NDUFS7:p.(Val179Met) as causative variant for the investigated canine Leigh syndrome.

Polioencephalopathy in Eurasier dogs

BACKGROUND

Polioencephalopathies secondary to inborn errors of metabolism have been described in dogs, but few genetically characterized.

OBJECTIVES

Clinically and genetically characterize polioencephalopathy in a family of Eurasier dogs.

ANIMALS

Three Eurasier dogs (littermates) presented with early onset movement disorders (9 weeks in 2, 4-6 months in 1). Progressive gait abnormalities were detected in 2 of the dogs, persistent divergent strabismus in 1, whereas consciousness and behavior remained intact in all dogs. One dog was euthanized at 25 months.

METHODS

Video footage was assessed in all dogs, and Dogs 1 and 2 had examinations and investigations performed. Whole genome sequencing of Dog 1 and further genetic analyses in the family were performed. A cohort of 115 Eurasier controls was genotyped for specific variants.

RESULTS

Episodes were characterized by generalized ataxia, as well as a hypermetric thoracic limb gait, dystonia, and irregular flexion and extension movements of the thoracic limbs. Magnetic resonance imaging of the brain in Dogs 1 and 2 identified symmetrical, bilateral T2 and fluid attenuated inversion recovery hyperintense, T1 hypo to isointense, nonenhancing lesions of the caudate nucleus, lateral and medial geniculate nuclei, thalamus, hippocampus, rostral colliculus and mild generalized brain atrophy. Genetic analyses identified a homozygous mitochondrial trans-2-enoyl-CoA reductase (MECR) missense variant in all 3 dogs, and a homozygous autophagy-related gene 4D (ATG4D) missense variant in Dogs 1 and 2.

CONCLUSIONS AND CLINICAL IMPORTANCE

We describe a presumed hereditary and progressive polioencephalopathy in a family of Eurasier dogs. Further research is needed to establish the role of the MECR gene in dogs and the pathogenic effects of the detected variants.

Multiocular defect in the Old English Sheepdog: A canine form of Stickler syndrome type II associated with a missense variant in the collagen-type gene COL11A1

Multiocular defect has been described in different canine breeds, including the Old English Sheepdog. Affected dogs typically present with multiple and various ocular abnormalities. We carried out whole genome sequencing on an Old English Sheepdog that had been diagnosed with hereditary cataracts at the age of five and then referred to a board-certified veterinary ophthalmologist due to owner-reported visual deterioration. An ophthalmic assessment revealed that there was bilateral vitreal degeneration, macrophthalmos, and spherophakia in addition to cataracts. Follow-up consultations revealed cataract progression, retinal detachment, uveitis and secondary glaucoma. Whole genome sequence filtered variants private to the case, shared with another Old English Sheepdog genome and predicted to be deleterious were genotyped in an initial cohort of six Old English Sheepdogs (three affected by multiocular defect and three control dogs without evidence of inherited eye disease). Only one of the twenty-two variants segregated correctly with multiocular defect. The variant is a single nucleotide substitution, located in the collagen-type gene COL11A1, c.1775T>C, that causes an amino acid change, p.Phe1592Ser. Genotyping of an additional 14 Old English Sheepdogs affected by multiocular defect revealed a dominant mode of inheritance with four cases heterozygous for the variant. Further genotyping of hereditary cataract-affected Old English Sheepdogs revealed segregation of the variant in eight out of nine dogs. In humans, variants in the COL11A1 gene are associated with Stickler syndrome type II, also dominantly inherited.

Hemophilia A in a litter of Border Collies caused by a one base pair deletion in the F8 gene

In dogs, hemophilia A is known to affect different breeds. This is a case report describing hemophilia A in a litter of Border Collies. A privately owned bitch and her puppies (n = 7) were presented to the referring veterinarian after acute hematoma formation in the male offspring (n = 3) following microchip implantation. Global coagulation testing, as well as determination of factor VIII and IX activity, were carried out. Based on the results, factor VIII deficiency was suspected. Two of the affected male puppies were euthanized within a few days. Genetic testing of the mother and the surviving male puppy resulted in the description of a deletion in exon 14 of the F8 gene. This c.3206delA variant leads to a frameshift in amino acid sequence and a premature stop codon (p.Asn1069IlefsTer7). The detection of the mutation and consequent testing of related dogs revealed that the deletion most likely had occurred spontaneously in the mother and had been transmitted to several of her offspring in different litters. Identified carriers were taken out of the breeding scheme. It is concluded that genetic testing in the context of suspected genetic disease can lead to preventive measures, including timely exclusion of carriers from breeding.

A SACS deletion variant in Great Pyrenees dogs causes autosomal recessive neuronal degeneration

ARSACS (autosomal recessive spastic ataxia of Charlevoix-Saguenay) is a human neurological disorder characterized by progressive cerebellar ataxia and peripheral neuropathy. A recently recognized disorder in Great Pyrenees dogs is similarly characterized by widespread central nervous system degeneration leading to progressive cerebellar ataxia and spasticity, combined with peripheral neuropathy. Onset of clinical signs occurred in puppies as young as 4 months of age, with slow progression over several years. A multi-generation pedigree suggested an autosomal recessive mode of inheritance. Histopathology revealed consistent cerebellar Purkinje cell degeneration, neuronal degeneration in brainstem nuclei, widespread spinal cord white matter degeneration, ganglion cell degeneration, inappropriately thin myelin sheaths or fully demyelinated peripheral nerve fibers, and normal or only mild patterns of denervation atrophy in skeletal muscles. Genome-wide single nucleotide polymorphism (SNP) genotype data was collected from 6 cases and 26 controls, where homozygosity mapping identified a 3.3 Mb region on CFA25 in which all cases were homozygous and all controls were either heterozygous or homozygous for alternate haplotypes. This region tagged the SACS gene where variants are known to cause ARSACS. Sanger sequencing of SACS in affected dogs identified a 4 bp deletion that causes a frame shift and truncates 343 amino acids from the C terminus of the encoded sacsin protein (p.Val4244AlafsTer32). Our clinical and histopathological descriptions of this canine disorder contribute to the description of human ARSACS and represents the first naturally occurring large animal model of this disorder.

A frameshift-deletion mutation in Reelin causes cerebellar hypoplasia in White Swiss Shepherd dogs

Cerebellar hypoplasia is a heterogeneous neurological condition in which the cerebellum is smaller than usual or not completely developed. The condition can have genetic origins, with Mendelian-effect mutations described in several mammalian species. Here, we describe a genetic investigation of cerebellar hypoplasia in White Swiss Shepherd dogs, where two affected puppies were identified from a litter with a recent common ancestor on both sides of their pedigree. Whole genome sequencing was conducted for 10 dogs in this family, and filtering of these data based on a recessive transmission hypothesis highlighted five protein-altering candidate variants - including a frameshift-deletion of the Reelin (RELN) gene (p.Val947*). Given the status of RELN as a gene responsible for cerebellar hypoplasia in humans, sheep and mice, these data strongly suggest the loss-of-function variant as underlying these effects. This variant has not been found in other dog breeds nor in a cohort of European White Swiss Shepherds, suggesting a recent mutation event. This finding will support the genotyping of a more diverse sample of dogs, and should aid future management of the harmful allele through optimised mating schemes.

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

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

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.

Identification of an ADAMTS2 frameshift variant in a cat family with Ehlers-Danlos syndrome

We investigated 4 European domestic shorthair kittens with skin lesions consistent with the dermatosparaxis type of the Ehlers-Danlos syndrome, a connective tissue disorder. The kittens were sired by the same tomcat but were born by 3 different mothers. The kittens had easily torn skin resulting in nonhealing skin wounds. Both clinically and histologically, the skin showed thin epidermis in addition to inflammatory changes. Changes in collagen fibers were visible in electron micrographs. The complete genome of an affected kitten was sequenced. A one base pair duplication leading to a frameshift in the candidate gene ADAMTS2 was identified, p.(Ser235fs*3). All 4 affected cats carried the frameshift duplication in a homozygous state. Genotypes at this variant showed perfect cosegregation with the autosomal recessive Ehlers-Danlos syndrome phenotype in the available family. The mutant allele did not occur in 48 unrelated control cats. ADAMTS2 loss-of-function variants cause autosomal recessive forms of Ehlers-Danlos syndrome in humans, mice, dogs, cattle, and sheep. The available evidence from our investigation together with the functional knowledge on ADAMTS2 in other species allows to classify the identified ADAMTS2 variant as pathogenic and most likely causative variant for the observed Ehlers-Danlos syndrome.

SGCD Missense Variant in a Lagotto Romagnolo Dog with Autosomal Recessively Inherited Limb-Girdle Muscular Dystrophy

An 8-month-old female Lagotto Romagnolo dog was presented for a 1-month history of an initial severe reluctance to move, rapidly progressing to a marked stiff gait and progressive muscular weakness and evolving to tetraparesis, which persuaded the owner to request euthanasia. A primary muscle pathology was supported by necropsy and histopathological findings. Macroscopically, the muscles were moderately atrophic, except for the diaphragm and the neck muscles, which were markedly thickened. Histologically, all the skeletal muscles examined showed atrophy, hypertrophy, necrosis with calcification of the fibers, and mild fibrosis and inflammation. On immunohistochemistry, all three dystrophin domains and sarcoglycan proteins were absent. On Western blot analysis, no band was present for delta sarcoglycan. We sequenced the genome of the affected dog and compared the data to more than 900 control genomes of different dog breeds. Genetic analysis revealed a homozygous private protein-changing variant in the SGCD gene encoding delta- sarcoglycan in the affected dog. The variant was predicted to induce a SGCD:p.(Leu242Pro) change in the protein. In silico tools predicted the change to be deleterious. Other 770 Lagotto Romagnolo dogs were genotyped for the variant and all found to be homozygous wild type. Based on current knowledge of gene function in other mammalian species, including humans, hamsters, and dogs, we propose the SGCD missense variant as the causative variant of the observed form of muscular dystrophy in the index case. The absence of the variant allele in the Lagotto Romagnolo breeding population indicates a rare allele that has appeared recently.

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.

Whole Animal Genome Sequencing: user-friendly, rapid, containerized pipelines for processing, variant discovery, and annotation of short-read whole genome sequencing data

Advancements in massively parallel short-read sequencing technologies and the associated decreasing costs have led to large and diverse variant discovery efforts across species. However, processing high-throughput short-read sequencing data can be challenging with potential pitfalls and bioinformatics bottlenecks in generating reproducible results. Although a number of pipelines exist that address these challenges, these are often geared toward human or traditional model organism species and can be difficult to configure across institutions. Whole Animal Genome Sequencing (WAGS) is an open-source set of user-friendly, containerized pipelines designed to simplify the process of identifying germline short (SNP and indel) and structural variants (SVs) geared toward the veterinary community but adaptable to any species with a suitable reference genome. We present a description of the pipelines [adapted from the best practices of the Genome Analysis Toolkit (GATK)], along with benchmarking data from both the preprocessing and joint genotyping steps, consistent with a typical user workflow.

RALGAPA1 Deletion in Belgian Shepherd Dogs with Cerebellar Ataxia

Several genetically distinct forms of cerebellar ataxia exist in Belgian shepherd dogs. We investigated a litter in which two puppies developed cerebellar ataxia. The clinical signs stabilized at around six weeks of age, but remained visible into adulthood. Combined linkage and homozygosity mapping delineated a 5.5 Mb critical interval. The comparison of whole-genome sequence data of one affected dog to 929 control genomes revealed a private homozygous ~4.8 kb deletion in the critical interval, Chr8:14,468,376_14,473,136del4761. The deletion comprises exon 35 of the RALGAPA1 gene, XM_038544497.1:c.6080-2893_6944+1003del. It is predicted to introduce a premature stop codon into the transcript, truncating ~23% of the wild-type open reading frame of the encoded Ral GTPase-activating protein catalytic subunit α 1, XP_038400425.1:(p.Val2027Glnfs*7). Genotypes at the deletion showed the expected co-segregation with the phenotype in the family. Genotyping additional ataxic Belgian shepherd dogs revealed three additional homozygous mutant dogs from a single litter, which had been euthanized at five weeks of age due to their severe clinical phenotype. Histopathology revealed cytoplasmic accumulation of granular material within cerebellar Purkinje cells. Genotyping a cohort of almost 900 Belgian shepherd dogs showed the expected genotype-phenotype association and a carrier frequency of 5% in the population. Human patients with loss-of-function variants in RALGAPA1 develop psychomotor disability and early-onset epilepsy. The available clinical and histopathological data, together with current knowledge about RALGAPA1 variants and their functional impact in other species, suggest the RALGAPA1 deletion is the likely causative defect for the observed phenotype in the affected dogs.

Natural disease history of a canine model of oligogenic RPGRIP1-cone-rod dystrophy establishes variable effects of previously and newly mapped modifier loci

Canine RPGRIP1-cone-rod dystrophy (CRD), a model for human inherited retinal diseases (IRDs), was originally identified as autosomal recessive early-onset blindness. However, later studies revealed extensive phenotypic variability among RPGRIP1 mutants. This led to the identification of a homozygous MAP9 variant as a modifier associated with early-onset disease. Based on further phenotypic variation affecting cone photoreceptor function, we report mapping of L3 as an additional modifier locus, within a 4.1-Mb locus on canine chromosome 30. We establish the natural disease history of RPGRIP1-CRD based on up to 9-year long-term functional and structural retinal data from 58 dogs including 44 RPGRIP1 mutants grouped according to the modifier status. RPGRIP1 mutants affected by both MAP9 and L3 modifiers exhibited the most severe phenotypes with rapid disease progression. MAP9 alone was found to act as an overall accelerator of rod and cone diseases, while L3 had a cone-specific effect. Ultrastructural analysis of photoreceptors revealed varying degrees of rod and cone damage, while the connecting cilia appeared structurally preserved in all groups. We conclude that RPGRIP1-CRD is an oligogenic disease with at least three loci contributing to the pathogenesis. While the RPGRIP1 variant is required for developing the disease, MAP9 and L3 modifiers exacerbate the phenotype, individually and cumulatively. Oligogenic canine RPGRIP1-CRD illustrates the impact of multiple genetic modifiers on disease phenotype and thus has the potential to reveal new targets for broad-spectrum therapies for oligogenic or polygenic forms of human IRDs.

Comprehensive analysis of geographic and breed-purpose influences on genetic diversity and inherited disease risk in the Doberman dog breed

BACKGROUND

Publicly available phenotype data and genotyping array data from two citizen science projects: "Doberman Health Surveys" and "The Doberman Diversity Project" were analyzed to explore relative homozygosity, diversity, and disorder risk according to geographical locale and breeding purpose in the Doberman.

RESULTS

From the phenotypic data cohort, life expectancy of a Doberman at birth is 9.1 years. The leading causes of death were heart disease (accounting for 28% of deaths) and cancers (collectively accounting for 14% of deaths). By genotyping, the world Doberman population exists as four major cohorts (European exhibition-bred, Americas exhibition-bred, European work, Americas pet/informal). Considering the entire Doberman population, four genomic regions longer than 500 Kb are fixed in 90% or more of 3,226 dogs included in this study. The four fixed regions reside on two autosomal chromosomes: CFA3:0.8-2.3 Mb (1.55 Mb); CFA3: 57.9-59.8 Mb (1.8 Mb); CFA31:0-1.2 Mb (1.2 Mb); and CFA31:4.80-6.47 Mb (1.67 Mb). Using public variant call files including variants for eight Doberman pinschers, we observed 30 potentially functional alternate variants that were evolutionarily diverged relative to the wider sequenced dog population within the four strongly homozygous chromosomal regions. Effective population size (Ne) is a statistical measure of breed diversity at the time of sampling that approximates the number of unique individuals. The major identified sub-populations of Dobermans demonstrated Ne in the range 70-236. The mean level of inbreeding in the Doberman breed is 40% as calculated by the number of array variants in runs of homozygosity divided by the assayed genome size (excluding the X chromosome). The lowest observed level of inbreeding in the Dobermans assayed was 15% in animals that were first generation mixes of European and USA bred Dobermans. Array variant analysis shows that inter-crossing between European and USA-bred Dobermans has capacity to re-introduce variation at many loci that are strongly homozygous.

CONCLUSIONS

We conclude that efforts to improve breed diversity first should focus on regions with the highest fixation levels, but managers must ensure that mutation loads are not worsened by increasing the frequencies of rarer haplotypes in the identified regions. The analysis of global data identified regions of strong fixation that might impact known disorder risks in the breed. Plausible gene candidates for future analysis of the genetic basis of cardiac disease and cancer were identified in the analysis.

PKD1 Nonsense Variant in a Lagotto Romagnolo Family with Polycystic Kidney Disease

A female Lagotto Romagnolo dog with polycystic kidney disease (PKD) and her progeny, including PKD-affected offspring, were studied. All affected dogs appeared clinically inconspicuous, while sonography revealed the presence of renal cysts. The PKD-affected index female was used for breeding and produced two litters with six affected offspring of both sexes and seven unaffected offspring. The pedigrees suggested an autosomal dominant mode of inheritance of the trait. A trio whole genome sequencing analysis of the index female and her unaffected parents identified a de novo heterozygous nonsense variant in the coding region of the PKD1 gene. This variant, NM_001006650.1:c.7195G>T, is predicted to truncate 44% of the open reading frame of the wild-type PKD1 protein, NP_001006651.1:p.(Glu2399*). The finding of a de novo variant in an excellent functional candidate gene strongly suggests that the PKD1 nonsense variant caused the observed phenotype in the affected dogs. Perfect co-segregation of the mutant allele with the PKD phenotype in two litters supports the hypothesized causality. To the best of our knowledge, this is the second description of a PKD1-related canine form of autosomal dominant PKD that may serve as an animal model for similar hepatorenal fibrocystic disorders in humans.

A TNR Frameshift Variant in Weimaraner Dogs with an Exercise-Induced Paroxysmal Movement Disorder

BACKGROUND

Some paroxysmal movement disorders remain without an identified genetic cause.

OBJECTIVES

The aim was to identify the causal genetic variant for a paroxysmal dystonia-ataxia syndrome in Weimaraner dogs.

METHODS

Clinical and diagnostic investigations were performed. Whole genome sequencing of one affected dog was used to identify private homozygous variants against 921 control genomes.

RESULTS

Four Weimaraners were presented for episodes of abnormal gait. Results of examinations and diagnostic investigations were unremarkable. Whole genome sequencing revealed a private frameshift variant in the TNR (tenascin-R) gene in an affected dog, XM_038542431.1:c.831dupC, which is predicted to truncate more than 75% of the open read frame. Genotypes in a cohort of 4 affected and 70 unaffected Weimaraners showed perfect association with the disease phenotype.

CONCLUSIONS

We report the association of a TNR variant with a paroxysmal dystonia-ataxia syndrome in Weimaraners. It might be relevant to include sequencing of this gene in diagnosing humans with unexplained paroxysmal movement disorders. © 2023 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

Case report: Sacral agenesis in two boxer dogs: clinical presentation, diagnostic investigations, and outcome

Two boxer dogs from the same litter were presented at 3 months of age for urinary and fecal incontinence. Both dogs had an abnormal tail consisting of a small stump, an atonic anal sphincter, and absent perineal reflex and sensation. Neurological evaluation was indicative of a lesion of the cauda equina or sacral spinal cord. Radiology and CT scan of the spine displayed similar findings in the two dogs that were indicative of sacral agenesis. Indeed, they had 6 lumbar vertebrae followed by a lumbosacral transitional vertebra, lacking a complete spinous process, and a hypoplastic vertebra carrying 2 hypoplastic sacral transverse processes as the only remnant of the sacral bone. Caudal vertebrae were absent in one of the dogs. On MRI, one dog had a dural sac occupying the entire spinal canal and ending in a subfascial fat structure. In the other dog, the dural sac finished in an extracanalar, subfascial, well-defined cystic structure, communicating with the subarachnoid space, and consistent with a meningocele. Sacral agenesis-that is the partial or complete absence of the sacral bones-is a neural tube defect occasionally reported in humans with spina bifida occulta. Sacral agenesis has been described in human and veterinary medicine in association with conditions such as caudal regression syndrome, perosomus elumbis, and Currarino syndrome. These neural tube defects are caused by genetic and/or environmental factors. Despite thorough genetic investigation, no candidate variants in genes with known functional impact on bone development or sacral development could be found in the affected dogs. To the best of the authors' knowledge, this is the first report describing similar sacral agenesis in two related boxer dogs.

An ABCC9 Missense Variant Is Associated with Sudden Cardiac Death and Dilated Cardiomyopathy in Juvenile Dogs

Sudden cardiac death in the young (SCDY) is a devastating event that often has an underlying genetic basis. Manchester Terrier dogs offer a naturally occurring model of SCDY, with sudden death of puppies as the manifestation of an inherited dilated cardiomyopathy (DCM). We performed a genome-wide association study for SCDY/DCM in Manchester Terrier dogs and identified a susceptibility locus harboring the cardiac ATP-sensitive potassium channel gene ABCC9. Sanger sequencing revealed an ABCC9 p.R1186Q variant present in a homozygous state in all SCDY/DCM-affected dogs (n = 26). None of the controls genotyped (n = 398) were homozygous for the variant, but 69 were heterozygous carriers, consistent with autosomal recessive inheritance with complete penetrance (p = 4 × 10-42 for the association of homozygosity for ABCC9 p.R1186Q with SCDY/DCM). This variant exists at low frequency in human populations (rs776973456) with clinical significance previously deemed uncertain. The results of this study further the evidence that ABCC9 is a susceptibility gene for SCDY/DCM and highlight the potential application of dog models to predict the clinical significance of human variants.

SOAT1 missense variant in two cats with sebaceous gland dysplasia

Spontaneously arisen hereditary diseases in domestic animals provide an excellent opportunity to study the physiological functions of the altered genes. We investigated two 4-month-old sibling domestic short haired kittens with dry dark debris around the eyes, nose, and ears, dark crusting on the legs and a thin poor hair coat. Skin biopsies revealed abnormal sebaceous gland morphology with lack of normal sebocyte arrangement and differentiation. Hair follicles had a distorted silhouette, interpreted as a change secondary to the observed sebaceous gland dysplasia. Whole genome sequencing on both affected kittens and 65 genetically diverse feline genomes was performed. Filtering for variants that were present in both kittens but absent from the control genomes revealed a homozygous missense variant in SOAT1, encoding sterol O-acyltransferase 1. The protein is localized in the endoplasmic reticulum and catalyzes the formation of cholesteryl esters, an essential component of sebum and meibum. The identified SOAT1:c.1531G > A variant is predicted to change a highly conserved glycine residue within the last transmembrane domain of SOAT1, p.Gly511Arg. In mice, variants in Soat1 or complete knockout of the gene lead to the "hair interior defect" (hid) or abnormal Meibomian glands, respectively. SOAT1:c.1531G > A represents a plausible candidate variant for the observed sebaceous gland dysplasia in both kittens of this study. The variant was not present in 10 additional cats with a similar clinical and histopathological phenotype suggesting genetic heterogeneity. SOAT1 variants should be considered as potential cause in hereditary sebaceous gland dysplasias of humans and domestic animals.

A clue to the etiology of disorders of sex development from identity-by-descent analysis in dogs with cryptic relatedness

Disorders of sex development (DSDs) are discrepancies between sex chromosomes and phenotypical sex. Quite common forms of DSD in canine populations include testicular and ovotesticular XX DSDs with a normal set of sex chromosomes. The objective of this study was to identify genes and putative harmful variants for canine XX DSDs. I have reanalyzed data from the whole-genome sequencing of 11 XX DSD French Bulldogs and six XX DSD American Staffordshire Terriers. Identity-by-descent analysis revealed cryptic relatedness in affected French Bulldogs. Causative genes were sought in chromosomal segments shared identical-by-descent by close relatives. In French Bulldogs, the reanalysis identified 19 regions of importance with a total length of just 65.9 Mb. Variant filtering within the regions implicated AKAP2, PIWIL1, POLR3A and SH2D4B as genes that may be involved in individual cases of testicular and ovotesticular XX DSD in French Bulldogs and American Staffordshire Terriers.