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.

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.

Ancestry dynamics and trait selection in a designer cat breed.

The Bengal cat breed was developed from intercrosses between the Asian leopard cat, Prionailurus bengalensis, and the domestic cat, Felis catus, with a last common ancestor approximately 6 million years ago. Predicted to derive ∼94% of their genome from domestic cats, regions of the leopard cat genome are thought to account for the unique pelage traits and ornate color patterns of the Bengal breed, which are similar to those of ocelots and jaguars. We explore ancestry distribution and selection signatures in the Bengal breed by using reduced representation and whole-genome sequencing from 947 cats. The mean proportion of leopard cat DNA in the Bengal breed is 3.48%, lower than predicted from breed history, and is broadly distributed, covering 93% of the Bengal genome. Overall, leopard cat introgressions do not show strong signatures of selection across the Bengal breed. However, two popular color traits in Bengal cats, charcoal and pheomelanin intensity, are explained by selection of leopard cat genes whose expression is reduced in a domestic cat background, consistent with genetic incompatibility resulting from hybridization. We characterize several selective sweeps in the Bengal genome that harbor candidate genes for pelage and color pattern and that are associated with domestic, rather than leopard, cat haplotypes. We identify the molecular and phenotypic basis of one selective sweep as reduced expression of the Fgfr2 gene, which underlies glitter, a trait desired by breeders that affects hair texture and light reflectivity.

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.

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.

Single-cell profiling of bronchoalveolar cells reveals a Th17 signature in neutrophilic severe equine asthma.

Severe equine asthma (SEA) is a complex respiratory condition characterized by chronic airway inflammation. It shares many clinical and pathological features with human neutrophilic asthma, making it a valuable model for studying this condition. However, the immune mechanisms driving SEA have remained elusive. Although SEA has been primarily associated with a Th2 response, there have also been reports of Th1, Th17, or mixed-mediated responses. To uncover the elusive immune mechanisms driving SEA, we performed single-cell mRNA sequencing (scRNA-seq) on cryopreserved bronchoalveolar cells from 11 Warmblood horses, 5 controls and 6 with SEA. We identified six major cell types, including B cells, T cells, monocytes-macrophages, dendritic cells, neutrophils, and mast cells. All cell types exhibited significant heterogeneity, with previously identified and novel cell subtypes. Notably, we observed monocyte-lymphocyte complexes and detected a robust Th17 signature in SEA, with CXCL13 upregulation in intermediate monocytes. Asthmatic horses exhibited expansion of the B-cell population, Th17 polarization of the T-cell populations, and dysregulation of genes associated with T-cell function. Neutrophils demonstrated enhanced migratory capacity and heightened aptitude for neutrophil extracellular trap formation. These findings provide compelling evidence for a predominant Th17 immune response in neutrophilic SEA, driven by dysregulation of monocyte and T-cell genes. The dysregulated genes identified through scRNA-seq have potential as biomarkers and therapeutic targets for SEA and provide insights into human neutrophilic asthma.

Independent COL17A1 Variants in Cats with Junctional Epidermolysis Bullosa.

Epidermolysis bullosa (EB), characterized by defective adhesion of the epidermis to the dermis, is a heterogeneous disease with many subtypes in human patients and domestic animals. We investigated two unrelated cats with recurring erosions and ulcers on ear pinnae, oral mucosa, and paw pads that were suggestive of EB. Histopathology confirmed the diagnosis of EB in both cats. Case 1 was severe and had to be euthanized at 5 months of age. Case 2 had a milder course and was alive at 11 years of age at the time of writing. Whole genome sequencing of both affected cats revealed independent homozygous variants in COL17A1 encoding the collagen type XVII alpha 1 chain. Loss of function variants in COL17A1 lead to junctional epidermolysis bullosa (JEB) in human patients. The identified splice site variant in case 1, c.3019+1del, was predicted to lead to a complete deficiency in collagen type XVII. Case 2 had a splice region variant, c.769+5G>A. Assessment of the functional impact of this variant on the transcript level demonstrated partial aberrant splicing with residual expression of wildtype transcript. Thus, the molecular analyses provided a plausible explanation of the difference in clinical severity between the two cases and allowed the refinement of the diagnosis in the affected cats to JEB. This study highlights the complexity of EB in animals and contributes to a better understanding of the genotype-phenotype correlation in COL17A1-related JEB.

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.

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.

Ancient segmentally duplicated LCORL retrocopies in equids.

LINE-1 is an active transposable element encoding proteins capable of inserting host gene retrocopies, resulting in retro-copy number variants (retroCNVs) between individuals. Here, we performed retroCNV discovery using 86 equids and identified 437 retrocopy insertions. Only 5 retroCNVs were shared between horses and other equids, indicating that the majority of retroCNVs inserted after the species diverged. A large number (17-35 copies) of segmentally duplicated Ligand Dependent Nuclear Receptor Corepressor Like (LCORL) retrocopies were present in all equids but absent from other extant perissodactyls. The majority of LCORL transcripts in horses and donkeys originate from the retrocopies. The initial LCORL retrotransposition occurred 18 million years ago (17-19 95% CI), which is coincident with the increase in body size, reduction in digit number, and changes in dentition that characterized equid evolution. Evolutionary conservation of the LCORL retrocopy segmental amplification in the Equidae family, high expression levels and the ancient timeline for LCORL retrotransposition support a functional role for this structural variant.

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.

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.