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.

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.

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.

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.

Abnormal keratinocyte differentiation in the nasal planum of Labrador Retrievers with hereditary nasal parakeratosis (HNPK)

Hereditary nasal parakeratosis (HNPK) is an inherited disorder described in Labrador Retrievers and Greyhounds. It has been associated with breed-specific variants in the SUV39H2 gene encoding a histone 3 methyltransferase involved in epigenetic silencing. Formalin-fixed biopsies of the nasal planum of Labrador Retrievers were screened by immunofluorescence microscopy for the presence and distribution of epidermal proliferation and differentiation markers. Gene expression of these markers was further analysed using RNA sequencing (RNA-seq) and ultrastructural epidermal differences were investigated by electron microscopy. Differentiation of the nasal planum in the basal and suprabasal epidermal layers of HNPK-affected dogs (n = 6) was similar compared to control dogs (n = 6). In the upper epidermal layers, clear modifications were noticed. Loricrin protein was absent in HNPK-affected nasal planum sections in contrast to sections of the same location of control dogs. However, loricrin was present in the epidermis of paw pads and abdominal skin from HNPK dogs and healthy control dogs. The patterns of keratins K1, K10 and K14, were not markedly altered in the nasal planum of HNPK-affected dogs while the expression of the terminal differentiation marker involucrin appeared less regular. Based on RNA-seq, LOR and IVL expression levels were significantly decreased, while KRT1, KRT10 and KRT14 levels were up-regulated (log2fold-changes of 2.67, 3.19 and 1.71, respectively) in HNPK-affected nasal planum (n = 3) compared to control dogs (n = 3). Electron microscopical analysis revealed structural alterations in keratinocytes and stratum corneum, and disrupted keratinocyte adhesions and distended intercellular spaces in lesional samples (n = 3) compared to a sample of a healthy control dog (n = 1). Our findings demonstrate aberrant keratinocyte terminal differentiation of the nasal planum of HNPK-affected Labrador Retrievers and provide insights into biological consequences of this inactive SUV39H2 gene variant.

A large deletion in the GP9 gene in Cocker Spaniel dogs with Bernard-Soulier syndrome

Inherited bleeding disorders including abnormalities of platelet number and function rarely occur in a variety of dog breeds, but are probably underdiagnosed. Genetically characterized canine forms of platelet disorders provide valuable large animal models for understanding similar platelet disorders in people. Breed-specific disease associated genetic variants in only eight different genes are known to cause intrinsic platelet disorders in dogs. However, the causative genetic variant in many dog breeds has until now remained unknown. Four cases of a mild to severe bleeding disorder in Cocker Spaniel dogs are herein presented. The affected dogs showed a platelet adhesion defect characterized by macrothrombocytopenia with variable platelet counts resembling human Bernard-Soulier syndrome (BSS). Furthermore, the lack of functional GPIb-IX-V was demonstrated by immunocytochemistry. Whole genome sequencing of one affected dog and visual inspection of the candidate genes identified a deletion in the glycoprotein IX platelet (GP9) gene. The GP9 gene encodes a subunit of a platelet surface membrane glycoprotein complex; this functions as a receptor for von Willebrand factor, which initiates the maintenance of hemostasis after injury. Variants in human GP9 are associated with Bernard-Soulier syndrome, type C. The deletion spanned 2460 bp, and included a significant part of the single coding exon of the canine GP9 gene on dog chromosome 20. The variant results in a frameshift and premature stop codon which is predicted to truncate almost two-thirds of the encoded protein. PCR-based genotyping confirmed recessive inheritance. The homozygous variant genotype seen in affected dogs did not occur in 98 control Cocker Spaniels. Thus, it was concluded that the structural variant identified in the GP9 gene was most likely causative for the BSS-phenotype in the dogs examined. These findings provide the first large animal GP9 model for this group of inherited platelet disorders and greatly facilitate the diagnosis and identification of affected and/or normal carriers in Cocker Spaniels.

Two MC1R loss-of-function alleles in cream-coloured Australian Cattle Dogs and white Huskies

Loss-of-function variants in the MC1R gene cause recessive red or yellow coat-colour phenotypes in many species. The canine MC1R:c.916C>T (p.Arg306Ter) variant is widespread and found in a homozygous state in many uniformly yellow- or red-coloured dogs. We investigated cream-coloured Australian Cattle Dogs whose coat colour could not be explained by this variant. A genome-wide association study with 10 cream and 123 red Australian Cattle Dogs confirmed that the cream locus indeed maps to MC1R. Whole-genome sequencing of cream dogs revealed a single nucleotide variant within the MITF binding site of the canine MC1R promoter. We propose to designate the mutant alleles at MC1R:c.916C>T as e¹ and at the new promoter variant as e² . Both alleles segregate in the Australian Cattle Dog breed. When we considered both alleles in combination, we observed perfect association between the MC1R genotypes and the cream coat colour phenotype in a cohort of 10 cases and 324 control dogs. Analysis of the MC1R transcript levels in an e¹ /e² compound heterozygous dog confirmed that the transcript levels of the e² allele were markedly reduced with respect to the e¹ allele. We further report another MC1R loss-of-function allele in Alaskan and Siberian Huskies caused by a 2-bp deletion in the coding sequence, MC1R:c.816_817delCT. We propose to term this allele e³ . Huskies that carry two copies of MC1R loss-of-function alleles have a white coat colour.

Canine NAPEPLD-associated models of human myelin disorders

Canine leukoencephalomyelopathy (LEMP) is a juvenile-onset neurodegenerative disorder of the CNS white matter currently described in Rottweiler and Leonberger dogs. Genome-wide association study (GWAS) allowed us to map LEMP in a Leonberger cohort to dog chromosome 18. Subsequent whole genome re-sequencing of a Leonberger case enabled the identification of a single private homozygous non-synonymous missense variant located in the highly conserved metallo-beta-lactamase domain of the N-acyl phosphatidylethanolamine phospholipase D (NAPEPLD) gene, encoding an enzyme of the endocannabinoid system. We then sequenced this gene in LEMP-affected Rottweilers and identified a different frameshift variant, which is predicted to replace the C-terminal metallo-beta-lactamase domain of the wild type protein. Haplotype analysis of SNP array genotypes revealed that the frameshift variant was present in diverse haplotypes in Rottweilers, and also in Great Danes, indicating an old origin of this second NAPEPLD variant. The identification of different NAPEPLD variants in dog breeds affected by leukoencephalopathies with heterogeneous pathological features, implicates the NAPEPLD enzyme as important in myelin homeostasis, and suggests a novel candidate gene for myelination disorders in people.

Genomic comparison of bovine papillomavirus 1 isolates from bovine, equine and asinine lesional tissue samples

Several attempts have been made to categorize equid- and bovid-specific bovine papillomavirus 1 (BPV1) isolates based on sequence tags. This study includes newly determined sequence information from 33 BPV1 isolates of equine, asinine and bovine origin and investigates sequence bias due to host species. Twenty of the viral genomes were sequenced over their entire length and a further thirteen were sequenced, including flanking sequences, at two specific sites, the LCR and the E5 ORF. Alignment and analyses of the sequences did not reveal statistically significant site differences between the sequences of bovine and equid origin. None of the proposed sites of divergence noted by other authors demonstrated significant species-specific characteristics. Our results suggest that BPV1 is shared between equine, asinine and bovine host species, and that viral transfer between bovines and equids is a repeated and ongoing phenomenon.

Canine Brachycephaly Is Associated with a Retrotransposon-Mediated Missplicing of SMOC2

In morphological terms, “form” is used to describe an object’s shape and size. In dogs, facial form is stunningly diverse. Facial retrusion, the proximodistal shortening of the snout and widening of the hard palate is common to brachycephalic dogs and is a welfare concern, as the incidence of respiratory distress and ocular trauma observed in this class of dogs is highly correlated with their skull form. Progress to identify the molecular underpinnings of facial retrusion is limited to association of a missense mutation in BMP3 among small brachycephalic dogs. Here, we used morphometrics of skull isosurfaces derived from 374 pedigree and mixed-breed dogs to dissect the genetics of skull form. Through deconvolution of facial forms, we identified quantitative trait loci that are responsible for canine facial shapes and sizes. Our novel insights include recognition that the FGF4 retrogene insertion, previously associated with appendicular chondrodysplasia, also reduces neurocranium size. Focusing on facial shape, we resolved a quantitative trait locus on canine chromosome 1 to a 188-kb critical interval that encompasses SMOC2. An intronic, transposable element within SMOC2 promotes the utilization of cryptic splice sites, causing its incorporation into transcripts, and drastically reduces SMOC2 gene expression in brachycephalic dogs. SMOC2 disruption affects the facial skeleton in a dose-dependent manner. The size effects of the associated SMOC2 haplotype are profound, accounting for 36% of facial length variation in the dogs we tested. Our data bring new focus to SMOC2 by highlighting its clinical implications in both human and veterinary medicine.

•

A population-based genetics study of dogs that required diagnostic imaging

•

Resolution of a QTL associated with face length reduction (brachycephaly)

•

Association of brachycephaly with a retrotransposon that disrupts SMOC2 splicing

•

The SMOC2 locus explains up to 36% of face length variation in dogs

Neuronal ceroid lipofuscinosis (NCL) is caused by the entire deletion of CLN8 in the Alpenländische Dachsbracke dog

Neuronal ceroid lipofuscinoses (NCLs) are inherited lysosomal storage diseases that have been described in a variety of dog breeds, where they are caused by different mutations in different genes. However, the causative gene defect in the breed Alpenländische Dachsbracke remained unknown so far. Here we present two confirmed cases of NCL in Alpenländische Dachsbracke dogs from different litters of the same sire with a different dam harboring the same underlying novel mutation in the CLN8 gene. Case 1, a 2-year-old male Alpenländische Dachsbracke was presented with neurological signs including disorientation, character changes including anxiety states and aggressiveness, sudden blindness and reduction of food intake. Magnetic resonance imaging (MRI) scans showed cerebral atrophy with dilation of all cerebral ventricles, thinning of the intermediate mass of the thalamus and widening of the cerebral sulci. Postmortem examination of the central nervous system (CNS) showed neuronal loss in the cerebral cortex, cerebellum and spinal cord with massive intracellular deposits of ceroid pigment. Additional ceroid-lipofuscin deposits were observed in the enteric nervous system and in macrophages within spleen, lymph nodes and lung. Ultrastructural analyses confirmed NCL with the presence of osmiophilic membrane bounded lamellar-like structures. Case 2, a 1,5-year old female Alpenländische Dachsbracke was presented with progressive generalized forebrain disease including mental changes such as fearful reactions to various kinds of external stimuli and disorientation. The dog also displayed seizures, absence of menace reactions and negative cotton-ball test with normal pupillary light reactions. The clinical and post mortem examination yielded similar results in the brain as in Case 1. Whole genome sequencing of Case 1 and PCR results of both cases revealed a homozygous deletion encompassing the entire CLN8 gene as the most likely causative mutation for the NCL form observed in both cases. The deletion follows recessive inheritance since the dam and a healthy male littermate of Case 1 were tested as heterozygous carriers. This is the first detailed description of CLN8 gene associated NCL in Alpenländische Dachsbracke dogs and thus provides a novel canine CLN8 model for this lysosomal storage disease. The presence of ceroid lipofuscin in extracerebral tissues may help to confirm the diagnosis of NCL in vivo, especially in new dog breeds where the underlying mutation is not known.

Neuronal ceroid lipofuscinosis in an adult American Staffordshire Terrier

A female, 5-year-old American Staffordshire Terrier with severe progressive neurological deficits, particularly in terms of ataxia and keeping balance, was examined pathomorphologically and a genetic analysis was performed. In neurons of various localizations of the central nervous system an accumulation of a finely granular pale eosinophilic or light brown material was found. In addition, the cerebellum revealed marked degeneration and loss of Purkinje and inner granule cells. The accumulated PAS-positive, argyrophilic, autofluorescent material showed ultrastructurally a lamellar appearance suggestive of lipofuscin. Genetic analysis revealed the presence of a sequence variant in the ARSG gene encoding the lysosomal enzyme arylsulfatase G. This case report describes an adult-onset of a neuronal ceroid lipofuscinosis that shows similarities with a human disorder termed Kufs disease.

A missense change in the ATG4D gene links aberrant autophagy to a neurodegenerative vacuolar storage disease

Inherited neurodegenerative disorders are debilitating diseases that occur across different species. We have performed clinical, pathological and genetic studies to characterize a novel canine neurodegenerative disease present in the Lagotto Romagnolo dog breed. Affected dogs suffer from progressive cerebellar ataxia, sometimes accompanied by episodic nystagmus and behavioral changes. Histological examination revealed unique pathological changes, including profound neuronal cytoplasmic vacuolization in the nervous system, as well as spheroid formation and cytoplasmic aggregation of vacuoles in secretory epithelial tissues and mesenchymal cells. Genetic analyses uncovered a missense change, c.1288G>A; p.A430T, in the autophagy-related ATG4D gene on canine chromosome 20 with a highly significant disease association (p = 3.8 x 10^-136) in a cohort of more than 2300 Lagotto Romagnolo dogs. ATG4D encodes a poorly characterized cysteine protease belonging to the macroautophagy pathway. Accordingly, our histological analyses indicated altered autophagic flux in affected tissues. The knockdown of the zebrafish homologue atg4da resulted in a widespread developmental disturbance and neurodegeneration in the central nervous system. Our study describes a previously unknown canine neurological disease with particular pathological features and implicates the ATG4D protein as an important autophagy mediator in neuronal homeostasis. The canine phenotype serves as a model to delineate the disease-causing pathological mechanism(s) and ATG4D function, and can also be used to explore treatment options. Furthermore, our results reveal a novel candidate gene for human neurodegeneration and enable the development of a genetic test for veterinary diagnostic and breeding purposes.

Neurodegenerative disorders affect millions of people worldwide. We describe a novel neurodegenerative disease in a canine model, characterized by progressive cerebellar ataxia and cellular vacuolization. Our genetic analyses identified a single nucleotide change in the autophagy-related ATG4D gene in affected dogs. The ATG4D gene has not been linked to inherited diseases before. The autophagy-lysosome pathway plays an important role in degrading and recycling different cellular components. Disturbed autophagy has been reported in several different diseases but mutations in core autophagy components are rare. Histological analyses of affected canine brain tissues revealed altered autophagic flux, and a knockdown of the gene in the zebrafish model caused marked neurodevelopmental alterations and neurodegeneration. Our findings identify a new disease-causing pathway and implicate the ATG4D protease as an important mediator for neuronal homeostasis. Furthermore, our study establishes a large animal model to investigate the role of ATG4D in autophagy and to test possible treatment options.

The transcriptome of equine peripheral blood mononuclear cells

Complete transcriptomic data at high resolution are available only for a few model organisms with medical importance. The gene structures of non-model organisms are mostly computationally predicted based on comparative genomics with other species. As a result, more than half of the horse gene models are known only by projection. Experimental data supporting these gene models are scarce. Moreover, most of the annotated equine genes are single-transcript genes. Utilizing RNA sequencing (RNA-seq) the experimental validation of predicted transcriptomes has become accessible at reasonable costs. To improve the horse genome annotation we performed RNA-seq on 561 samples of peripheral blood mononuclear cells (PBMCs) derived from 85 Warmblood horses. The mapped sequencing reads were used to build a new transcriptome assembly. The new assembly revealed many alternative isoforms associated to known genes or to those predicted by the Ensembl and/or Gnomon pipelines. We also identified 7,531 transcripts not associated with any horse gene annotated in public databases. Of these, 3,280 transcripts did not have a homologous match to any sequence deposited in the NCBI EST database suggesting horse specificity. The unknown transcripts were categorized as coding and noncoding based on predicted coding potential scores. Among them 230 transcripts had high coding potential score, at least 2 exons, and an open reading frame of at least 300 nt. We experimentally validated 9 new equine coding transcripts using RT-PCR and Sanger sequencing. Our results provide valuable detailed information on many transcripts yet to be annotated in the horse genome.

A deletion in the VLDLR gene in Eurasier dogs with cerebellar hypoplasia resembling a Dandy-Walker-like malformation (DWLM)

Dandy-Walker-like malformation (DWLM) is the result of aberrant brain development and mainly characterized by cerebellar hypoplasia. DWLM affected dogs display a non-progressive cerebellar ataxia. Several DWLM cases were recently observed in the Eurasier dog breed, which strongly suggested a monogenic autosomal recessive inheritance in this breed. We performed a genome-wide association study (GWAS) with 9 cases and 11 controls and found the best association of DWLM with markers on chromosome 1. Subsequent homozygosity mapping confirmed that all 9 cases were homozygous for a shared haplotype in this region, which delineated a critical interval of 3.35 Mb. We sequenced the genome of an affected Eurasier and compared it with the Boxer reference genome and 47 control genomes of dogs from other breeds. This analysis revealed 4 private non-synonymous variants in the critical interval of the affected Eurasier. We genotyped these variants in additional dogs and found perfect association for only one of these variants, a single base deletion in the VLDLR gene encoding the very low density lipoprotein receptor. This variant, VLDLR:c.1713delC is predicted to cause a frameshift and premature stop codon (p.W572Gfs*10). Variants in the VLDLR gene have been shown to cause congenital cerebellar ataxia and mental retardation in human patients and Vldlr knockout mice also display an ataxia phenotype. Our combined genetic data together with the functional knowledge on the VLDLR gene from other species thus strongly suggest that VLDLR:c.1713delC is indeed causing DWLM in Eurasier dogs.

Congenital hepatic fibrosis in the Franches-Montagnes horse is associated with the polycystic kidney and hepatic disease 1 (PKHD1) gene

Congenital hepatic fibrosis has been described as a lethal disease with monogenic autosomal recessive inheritance in the Swiss Franches-Montagnes horse breed. We performed a genome-wide association study with 5 cases and 12 controls and detected an association on chromosome 20. Subsequent homozygosity mapping defined a critical interval of 952 kb harboring 10 annotated genes and loci including the polycystic kidney and hepatic disease 1 (autosomal recessive) gene (PKHD1). PKHD1 represents an excellent functional candidate as variants in this gene were identified in human patients with autosomal recessive polycystic kidney and hepatic disease (ARPKD) as well as several mouse and rat mutants. Whereas most pathogenic PKHD1 variants lead to polycystic defects in kidney and liver, a small subset of the human ARPKD patients have only liver symptoms, similar to our horses with congenital hepatic fibrosis. The PKHD1 gene is one of the largest genes in the genome with multiple alternative transcripts that have not yet been fully characterized.

We sequenced the genomes of an affected foal and 46 control horses to establish a comprehensive list of variants in the critical interval. We identified two missense variants in the PKHD1 gene which were strongly, but not perfectly associated with congenital hepatic fibrosis. We speculate that reduced penetrance and/or potential epistatic interactions with hypothetical modifier genes may explain the imperfect association of the detected PKHD1 variants. Our data thus indicate that horses with congenital hepatic fibrosis represent an interesting large animal model for the liver-restricted subtype of human ARPKD.

Hepatic fungal infection in a young beagle with unrecognised hereditary cobalamin deficiency (Imerslund-Gräsbeck syndrome)

A 12-month-old beagle presented for anorexia, pyrexia and vomiting. The dog had been treated intermittently with antibiotics and corticosteroids for inappetence and lethargy since five months of age. Previous laboratory abnormalities included macrocytosis and neutropenia. At presentation, the dog was lethargic, febrile and thin. Laboratory examination findings included anaemia, a left shift, thrombocytopenia, hypoglycaemia and hyperbilirubinaemia. Multiple, small, hypoechoic, round hepatic lesions were observed on abdominal ultrasound. Cytological examination of hepatic fine needle aspirates revealed a fungal infection and associated pyogranulomatous inflammation. The dog's general condition deteriorated despite supportive measures and treatment with fluconazole, and owners opted for euthanasia before hypocobalaminaemia was identified. Subsequent genomic analysis revealed a CUBN:c.786delC mutation in a homozygous state, confirming hereditary cobalamin malabsorption (Imerslund-Gräsbeck syndrome). Similar to human infants, dogs with Imerslund-Gräsbeck syndrome may rarely be presented for infectious diseases, distracting focus from the underlying primary disorder.

A mutation in the FAM83G gene in dogs with hereditary footpad hyperkeratosis (HFH)

Hereditary footpad hyperkeratosis (HFH) represents a palmoplantar hyperkeratosis, which is inherited as a monogenic autosomal recessive trait in several dog breeds, such as e.g. Kromfohrländer and Irish Terriers. We performed genome-wide association studies (GWAS) in both breeds. In Kromfohrländer we obtained a single strong association signal on chromosome 5 (p_raw = 1.0×10⁻¹³) using 13 HFH cases and 29 controls. The association signal replicated in an independent cohort of Irish Terriers with 10 cases and 21 controls (p_raw = 6.9×10⁻¹⁰). The analysis of shared haplotypes among the combined Kromfohrländer and Irish Terrier cases defined a critical interval of 611 kb with 13 predicted genes. We re-sequenced the genome of one affected Kromfohrländer at 23.5× coverage. The comparison of the sequence data with 46 genomes of non-affected dogs from other breeds revealed a single private non-synonymous variant in the critical interval with respect to the reference genome assembly. The variant is a missense variant (c.155G>C) in the FAM83G gene encoding a protein with largely unknown function. It is predicted to change an evolutionary conserved arginine into a proline residue (p.R52P). We genotyped this variant in a larger cohort of dogs and found perfect association with the HFH phenotype. We further studied the clinical and histopathological alterations in the epidermis in vivo. Affected dogs show a moderate to severe orthokeratotic hyperplasia of the palmoplantar epidermis. Thus, our data provide the first evidence that FAM83G has an essential role for maintaining the integrity of the palmoplantar epidermis.

The palms and soles of mammals are covered by the palmoplantar epidermis, which has to bear immense mechanical forces and has therefore a special composition in comparison to the epidermis on regular skin. We studied a Mendelian disease in dogs, termed hereditary footpad hyperkeratosis (HFH). HFH affected dogs develop deep fissures in the paw pads, which are the consequence of a pathological thickening of the outermost layer of the epidermis. We mapped the disease causing genetic variant in the Kromfohrländer and Irish Terrier breeds to a 611 kb interval on chromosome 5. HFH affected Kromfohrländer and Irish Terriers shared the same haplotype indicating descent from a common founder. We re-sequenced the genome of an affected dog and compared it to genome sequences of 46 control dogs. The HFH affected dog had only one private non-synonymous variant in the critical interval, a missense variant of the FAM83G gene. We genotyped this variant in more than 500 dogs and found perfect association with the HFH phenotype. Our data very strongly suggest that the FAM83G variant is causative for HFH. FAM83G is a protein with unknown biochemical function. Our study thus provides the first link between this protein and the palmoplantar epidermis.

A mutation in the SUV39H2 gene in Labrador Retrievers with hereditary nasal parakeratosis (HNPK) provides insights into the epigenetics of keratinocyte differentiation

Hereditary nasal parakeratosis (HNPK), an inherited monogenic autosomal recessive skin disorder, leads to crusts and fissures on the nasal planum of Labrador Retrievers. We performed a genome-wide association study (GWAS) using 13 HNPK cases and 23 controls. We obtained a single strong association signal on chromosome 2 (p(raw) = 4.4×10⁻¹⁴). The analysis of shared haplotypes among the 13 cases defined a critical interval of 1.6 Mb with 25 predicted genes. We re-sequenced the genome of one case at 38× coverage and detected 3 non-synonymous variants in the critical interval with respect to the reference genome assembly. We genotyped these variants in larger cohorts of dogs and only one was perfectly associated with the HNPK phenotype in a cohort of more than 500 dogs. This candidate causative variant is a missense variant in the SUV39H2 gene encoding a histone 3 lysine 9 (H3K9) methyltransferase, which mediates chromatin silencing. The variant c.972T>G is predicted to change an evolutionary conserved asparagine into a lysine in the catalytically active domain of the enzyme (p.N324K). We further studied the histopathological alterations in the epidermis in vivo. Our data suggest that the HNPK phenotype is not caused by hyperproliferation, but rather delayed terminal differentiation of keratinocytes. Thus, our data provide evidence that SUV39H2 is involved in the epigenetic regulation of keratinocyte differentiation ensuring proper stratification and tight sealing of the mammalian epidermis.

IL26 gene inactivation in Equidae

Interleukin-26 (IL26) is a member of the IL10 cytokine family. The IL26 gene is located between two other well-known cytokines genes of this family encoding interferon-gamma (IFNG) and IL22 in an evolutionary conserved gene cluster. In contrast to humans and most other mammals, mice lack a functional Il26 gene. We analyzed the genome sequences of other vertebrates for the presence or absence of functional IL26 orthologs and found that the IL26 gene has also become inactivated in several equid species. We detected a one-base pair frameshift deletion in exon 2 of the IL26 gene in the domestic horse (Equus caballus), Przewalski horse (Equus przewalskii) and donkey (Equus asinus). The remnant IL26 gene in the horse is still transcribed and gives rise to at least five alternative transcripts. None of these transcripts share a conserved open reading frame with the human IL26 gene. A comparative analysis across diverse vertebrates revealed that the IL26 gene has also independently been inactivated in a few other mammals, including the African elephant and the European hedgehog. The IL26 gene thus appears to be highly variable, and the conserved open reading frame has been lost several times during mammalian evolution.

Magnetic resonance imaging and genetic investigation of a case of Rottweiler leukoencephalomyelopathy

Background

Leukoencephalomyelopathy is an inherited neurodegenerative disorder that affects the white matter of the spinal cord and brain and is known to occur in the Rottweiler breed. Due to the lack of a genetic test for this disorder, post mortem neuropathological examinations are required to confirm the diagnosis. Leukoencephalopathy with brain stem and spinal cord involvement and elevated lactate levels is a rare, autosomal recessive disorder in humans that was recently described to have clinical features and magnetic resonance imaging (MRI) findings that are similar to the histopathologic lesions that define leukoencephalomyelopathy in Rottweilers. Leukoencephalopathy with brain stem and spinal cord involvement is caused by mutations in the DARS2 gene, which encodes a mitochondrial aspartyl-tRNA synthetase. The objective of this case report is to present the results of MRI and candidate gene analysis of a case of Rottweiler leukoencephalomyelopathy to investigate the hypothesis that leukoencephalomyelopathy in Rottweilers could serve as an animal model of human leukoencephalopathy with brain stem and spinal cord involvement.

Case presentation

A two-and-a-half-year-old male purebred Rottweiler was evaluated for generalised progressive ataxia with hypermetria that was most evident in the thoracic limbs. MRI (T2-weighted) demonstrated well-circumscribed hyperintense signals within both lateral funiculi that extended from the level of the first to the sixth cervical vertebral body. A neurodegenerative disorder was suspected based on the progressive clinical course and MRI findings, and Rottweiler leukoencephalomyelopathy was subsequently confirmed via histopathology. The DARS2 gene was investigated as a causative candidate, but a sequence analysis failed to identify any disease-associated variants in the DNA sequence.

Conclusion

It was concluded that MRI may aid in the pre-mortem diagnosis of suspected cases of leukoencephalomyelopathy. Genes other than DARS2 may be involved in Rottweiler leukoencephalomyelopathy and may also be relevant in human leukoencephalopathy with brain stem and spinal cord involvement.

A locus on chromosome 5 is associated with dilated cardiomyopathy in Doberman Pinschers

Dilated cardiomyopathy (DCM) is a heterogeneous group of heart diseases with a strong genetic background. Currently, many human DCM cases exist where no causative mutation can be identified. DCM also occurs with high prevalence in several large dog breeds. In the Doberman Pinscher a specific DCM form characterized by arrhythmias and/or echocardiographic changes has been intensively studied by veterinary cardiologists. We performed a genome-wide association study in Doberman Pinschers. Using 71 cases and 70 controls collected in Germany we identified a genome-wide significant association to DCM on chromosome 5. We validated the association in an independent cohort collected in the United Kingdom. There is no known DCM candidate gene under the association signal. Therefore, DCM in Doberman Pinschers offers the chance of identifying a novel DCM gene that might also be relevant for human health.

Polymorphisms in the ABCB1 gene in phenobarbital responsive and resistant idiopathic epileptic Border Collies

BACKGROUND

Variation in the ABCB1 gene is believed to play a role in drug resistance in epilepsy.

HYPOTHESIS/OBJECTIVES

Variation in the ABCB1 gene encoding the permeability-glycoprotein could have an influence on phenobarbital (PB) resistance, which occurs with high frequency in idiopathic epileptic Border Collies (BCs).

ANIMALS

Two hundred and thirty-six client-owned BCs from Switzerland and Germany including 25 with idiopathic epilepsy, of which 13 were resistant to PB treatment.

METHODS

Prospective and retrospective case-control study. Data were collected retrospectively regarding disease status, antiepileptic drug (AED) therapy, and drug responsiveness. The frequency of a known mutation in the ABCB1 gene (4 base-pair deletion in the ABCB1 gene [c.296_299del]) was determined in all BCs. Additionally, the ABCB1 coding exons and flanking sequences were completely sequenced to search for additional variation in 41 BCs. Association analyses were performed in 2 case-control studies: idiopathic epileptic and control BCs and PB-responsive and resistant idiopathic epileptic BCs.

RESULTS

One of 236 BCs (0.4%) was heterozygous for the mutation in the ABCB1 gene (c.296_299del). A total of 23 variations were identified in the ABCB1 gene: 4 in exons and 19 in introns. The G-allele of the c.-6-180T > G variation in intron 1 was significantly more frequent in epileptic BCs resistant to PB treatment than in epileptic BCs responsive to PB treatment (P(raw) = .0025).

CONCLUSIONS AND CLINICAL IMPORTANCE

A variation in intron 1 of the ABCB1 gene is associated with drug responsiveness in BCs. This might indicate that regulatory mutations affecting the expression level of ABCB1 could exist, which may influence the reaction of a dog to AEDs.

A whole-genome scan for recurrent airway obstruction in Warmblood sport horses indicates two positional candidate regions

Recurrent airway obstruction (RAO), or heaves, is a naturally occurring asthma-like disease that is related to sensitisation and exposure to mouldy hay and has a familial basis with a complex mode of inheritance. A genome-wide scanning approach using two half-sibling families was taken in order to locate the chromosome regions that contribute to the inherited component of this condition in these families. Initially, a panel of 250 microsatellite markers, which were chosen as a well-spaced, polymorphic selection covering the 31 equine autosomes, was used to genotype the two half-sibling families, which comprised in total 239 Warmblood horses. Subsequently, supplementary markers were added for a total of 315 genotyped markers. Each half-sibling family is focused around a severely RAO-affected stallion, and the phenotype of each individual was assessed for RAO and related signs, namely, breathing effort at rest, breathing effort at work, coughing, and nasal discharge, using an owner-based questionnaire. Analysis using a regression method for half-sibling family structures was performed using RAO and each of the composite clinical signs separately; two chromosome regions (on ECA13 and ECA15) showed a genome-wide significant association with RAO at P < 0.05. An additional 11 chromosome regions showed a more modest association. This is the first publication that describes the mapping of genetic loci involved in RAO. Several candidate genes are located in these regions, a number of which are interleukins. These are important signalling molecules that are intricately involved in the control of the immune response and are therefore good positional candidates.

A comparative radiation hybrid map of sheep chromosome 10

Comparative radiation hybrid (RH) maps of individual ovine chromosomes are essential to identify genes governing traits of economic importance in sheep, a livestock species for which whole genome sequence data are not yet available. The USUoRH5000 radiation hybrid panel was used to generate a RH map of sheep chromosome 10 (OAR10) with 59 markers that span 1,422 cR over an estimated 92 Mb of the chromosome, thus providing markers every 2 Mb (equivalent to every 24 cR). The markers were derived from 46 BAC end sequences (BESs), a single EST, and 12 microsatellites. Comparative analysis showed that OAR10 shares remarkable conservation of gene order along the entire length of cattle chromosome 12 and that OAR10 contains four major homologous synteny blocks, each related to segments of the homologous human chromosome 13. Extending the comparison to the horse, dog, mouse, and chicken genome showed that these blocks share conserved synteny across species.

Microchip capillary electrophoresis-based genetic comparison of closely related cyathostomin nematode parasites of horses using randomly amplified polymorphic DNA polymerase chain reaction

The microchip-based capillary electrophoresis technology represents a valuable recent development for the analysis of complex DNA banding patterns. We have used this technology for the differentiation of the closely related cyathostomin species Cylicocyclus elongatus and C. insigne from the horse. We found that the Agilent 2100 bioanalyser in combination with the DNA 7500 Lab Chip were suited to perform a phylogenetic DNA fingerprinting analysis of the parasite species studied. The analysis of the electrophoretic data was optimised and it was possible to resolve a phylogenetic tree where all 12 individual worms of the two Cylicocyclus species studied were assigned to their species as determined by microscopic identification based on morphological traits. Thus, our data indicated that the procedure described here provides an additional powerful tool that can be employed for species delineation of closely related strains or species, such as the two taxa of Cylicocyclus investigated in the present study. Furthermore, by determining the second internal transcribed spacer region of three and nine individual worms for C. elongatus and C. insigne, respectively, low intraspecific variations of only up to 0.3% were demonstrated.