The DoGA consortium expression atlas of promoters and genes in 100 canine tissues.

The dog, Canis lupus familiaris, is an important model for studying human diseases. Unlike many model organisms, the dog genome has a comparatively poor functional annotation, which hampers gene discovery for development, morphology, disease, and behavior. To fill this gap, we established a comprehensive tissue biobank for both the dog and wolf samples. The biobank consists of 5485 samples representing 132 tissues from 13 dogs, 12 dog embryos, and 24 wolves. In a subset of 100 tissues from nine dogs and 12 embryos, we characterized gene expression activity for each promoter, including alternative and novel, i.e., previously not annotated, promoter regions, using the 5' targeting RNA sequencing technology STRT2-seq. We identified over 100,000 promoter region candidates in the recent canine genome assembly, CanFam4, including over 45,000 highly reproducible sites with gene expression and respective tissue enrichment levels. We provide a promoter and gene expression atlas with interactive, open data resources, including a data coordination center and genome browser track hubs. We demonstrated the applicability of Dog Genome Annotation (DoGA) data and resources using multiple examples spanning canine embryonic development, morphology and behavior, and diseases across species.

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.

Dominant Distal Myopathy 3 (MPD3) Caused by a Deletion in the HNRNPA1 Gene.

BACKGROUND AND OBJECTIVES: To determine the genetic cause of the disease in the previously reported family with adult-onset autosomal dominant distal myopathy (myopathy, distal, 3; MPD3). METHODS: Continued clinical evaluation including muscle MRI and muscle pathology. A linkage analysis with single nucleotide polymorphism arrays and genome sequencing were used to identify the genetic defect, which was verified by Sanger sequencing. RNA sequencing was used to investigate the transcriptional effects of the identified genetic defect. RESULTS: Small hand muscles (intrinsic, thenar, and hypothenar) were first involved with spread to the lower legs and later proximal muscles. Dystrophic changes with rimmed vacuoles and cytoplasmic inclusions were observed in muscle biopsies at advanced stage. A single nucleotide polymorphism array confirmed the previous microsatellite-based linkage to 8p22-q11 and 12q13-q22. Genome sequencing of three affected family members combined with structural variant calling revealed a small heterozygous deletion of 160 base pairs spanning the second last exon 10 of the heterogeneous nuclear ribonucleoprotein A1 (HNRNPA1) gene, which is in the linked region on chromosome 12. Segregation of the mutation with the disease was confirmed by Sanger sequencing. RNA sequencing showed that the mutant allele produces a shorter mutant mRNA transcript compared with the wild-type allele. Immunofluorescence studies on muscle biopsies revealed small p62 and larger TDP-43 inclusions. DISCUSSION: A small exon 10 deletion in the gene HNRNPA1 was identified as the cause of MPD3 in this family. The new HNRNPA1-related phenotype, upper limb presenting distal myopathy, was thus confirmed, and the family displays the complexities of gene identification.

Dog colour patterns explained by modular promoters of ancient canid origin.

Distinctive colour patterns in dogs are an integral component of canine diversity. Colour pattern differences are thought to have arisen from mutation and artificial selection during and after domestication from wolves but important gaps remain in understanding how these patterns evolved and are genetically controlled. In other mammals, variation at the ASIP gene controls both the temporal and spatial distribution of yellow and black pigments. Here, we identify independent regulatory modules for ventral and hair cycle ASIP expression, and we characterize their action and evolutionary origin. Structural variants define multiple alleles for each regulatory module and are combined in different ways to explain five distinctive dog colour patterns. Phylogenetic analysis reveals that the haplotype combination for one of these patterns is shared with Arctic white wolves and that its hair cycle-specific module probably originated from an extinct canid that diverged from grey wolves more than 2 million years ago. Natural selection for a lighter coat during the Pleistocene provided the genetic framework for widespread colour variation in dogs and wolves.

Missense variant in LOXHD1 is associated with canine nonsyndromic hearing loss.

Hearing loss is a common sensory deficit in both humans and dogs. In canines, the genetic basis is largely unknown, as genetic variants have only been identified for a syndromic form of hearing impairment. We observed a congenital or early-onset sensorineural hearing loss in a Rottweiler litter. Assuming an autosomal recessive inheritance, we used a combined approach of homozygosity mapping and genome sequencing to dissect the genetic background of the disorder. We identified a fully segregating missense variant in LOXHD1, a gene that is known to be essential for cochlear hair cell function and associated with nonsyndromic hearing loss in humans and mice. The canine LOXHD1 variant was specific to the Rottweiler breed in our study cohorts of pure-bred dogs. However, it also was present in some mixed-breed dogs, of which the majority showed Rottweiler ancestry. Low allele frequencies in these populations, 2.6% and 0.04%, indicate a rare variant. To summarize, our study describes the first genetic variant for canine nonsyndromic hearing loss, which is clinically and genetically similar to human LOXHD1-related hearing disorder, and therefore, provides a new large animal model for hearing loss. Equally important, the affected breed will benefit from a genetic test to eradicate this LOXHD1-related hearing disorder from the population.

In-frame deletion in canine PITRM1 is associated with a severe early-onset epilepsy, mitochondrial dysfunction and neurodegeneration.

We investigated the clinical, genetic, and pathological characteristics of a previously unknown severe juvenile brain disorder in several litters of Parson Russel Terriers. The disease started with epileptic seizures at 6-12 weeks of age and progressed rapidly to status epilepticus and death or euthanasia. Histopathological changes at autopsy were restricted to the brain. There was severe acute neuronal degeneration and necrosis diffusely affecting the grey matter throughout the brain with extensive intraneuronal mitochondrial crowding and accumulation of amyloid-ß (Aß). Combined homozygosity mapping and genome sequencing revealed an in-frame 6-bp deletion in the nuclear-encoded pitrilysin metallopeptidase 1 (PITRM1) encoding for a mitochondrial protease involved in mitochondrial targeting sequence processing and degradation. The 6-bp deletion results in the loss of two amino acid residues in the N-terminal part of PITRM1, potentially affecting protein folding and function. Assessment of the mitochondrial function in the affected brain tissue showed a significant deficiency in respiratory chain function. The functional consequences of the mutation were modeled in yeast and showed impaired growth in permissive conditions and an impaired respiration capacity. Loss-of-function variants in human PITRM1 result in a childhood-onset progressive amyloidotic neurological syndrome characterized by spinocerebellar ataxia with behavioral, psychiatric and cognitive abnormalities. Homozygous Pitrm1-knockout mice are embryonic lethal, while heterozygotes show a progressive, neurodegenerative phenotype characterized by impairment in motor coordination and Aß deposits. Our study describes a novel early-onset PITRM1-related neurodegenerative canine brain disorder with mitochondrial dysfunction, Aß accumulation, and lethal epilepsy. The findings highlight the essential role of PITRM1 in neuronal survival and strengthen the connection between mitochondrial dysfunction and neurodegeneration.

Intronic variant in POU1F1 associated with canine pituitary dwarfism.

The anterior pituitary gland secretes several endocrine hormones, essential for growth, reproduction and other basic physiological functions. Abnormal development or function of the pituitary gland leads to isolated or combined pituitary hormone deficiency (CPHD). At least 30 genes have been associated with human CPHD, including many transcription factors, such as POU1F1. CPHD occurs spontaneously also in mice and dogs. Two affected breeds have been reported in dogs: German Shepherds with a splice defect in the LHX3 gene and Karelian Bear Dogs (KBD) with an unknown genetic cause. We obtained samples from five KBDs presenting dwarfism and abnormal coats. A combined analysis of genome-wide association and next-generation sequencing mapped the disease to a region in chromosome 31 and identified a homozygous intronic variant in the fourth exon of the POU1F1 gene in the affected dogs. The identified variant, c.605-3C>A, resided in the splice region and was predicted to affect splicing. The variant's screening in three new prospective cases, related breeds, and ~ 8000 dogs from 207 breeds indicated complete segregation in KBDs with a carrier frequency of 8%, and high breed-specificity as carriers were found at a low frequency only in Lapponian Herders, a related breed. Our study establishes a novel canine model for CPHD with a candidate POU1F1 defect.

A missense variant in IFT122 associated with a canine model of retinitis pigmentosa.

Retinitis pigmentosa (RP) is a blinding eye disease affecting nearly two million people worldwide. Dogs are affected with a similar illness termed progressive retinal atrophy (PRA). Lapponian herders (LHs) are affected with several types of inherited retinal dystrophies, and variants in PRCD and BEST1 genes have been associated with generalized PRA and canine multifocal retinopathy 3 (cmr3), respectively. However, all retinal dystrophy cases in LHs are not explained by these variants, indicating additional genetic causes of disease in the breed. We collected DNA samples from 10 PRA affected LHs, with known PRCD and BEST1 variants excluded, and 34 unaffected LHs. A genome-wide association study identified a locus on CFA20 (praw = 2.4 × 10-7, pBonf = 0.035), and subsequent whole-genome sequencing of an affected LH revealed a missense variant, c.3176G>A, in the intraflagellar transport 122 (IFT122) gene. The variant was also found in Finnish Lapphunds, in which its clinical relevancy needs to be studied further. The variant interrupts a highly conserved residue, p.(R1059H), in IFT122 and likely impairs its function. Variants in IFT122 have not been associated with retinal degeneration in mammals, but the loss of ift122 in zebrafish larvae impaired opsin transport and resulted in progressive photoreceptor degeneration. Our study establishes a new spontaneous dog model to study the role of IFT122 in RP biology, while the affected breed will benefit from a genetic test for a recessive condition.

webGQT: A Shiny Server for Genotype Query Tools for Model-Based Variant Filtering.

SUMMARY: Genotype Query Tools (GQT) were developed to discover disease-causing variations from billions of genotypes and millions of genomes, processes data at substantially higher speed over other existing methods. While GQT has been available to a wide audience as command-line software, the difficulty of constructing queries among non-IT or non-bioinformatics researchers has limited its applicability. To overcome this limitation, we developed webGQT, an easy-to-use tool with a graphical user interface. With pre-built queries across three modules, webGQT allows for pedigree analysis, case-control studies, and population frequency studies. As a package, webGQT allows researchers with less or no applied bioinformatics/IT experience to mine potential disease-causing variants from billions. RESULTS: webGQT offers a flexible and easy-to-use interface for model-based candidate variant filtering for Mendelian diseases from thousands to millions of genomes at a reduced computation time. Additionally, webGQT provides adjustable parameters to reduce false positives and rescue missing genotypes across all modules. Using a case study, we demonstrate the applicability of webGQT to query non-human genomes. In addition, we demonstrate the scalability of webGQT on large data sets by implementing complex population-specific queries on the 1000 Genomes Project Phase 3 data set, which includes 8.4 billion variants from 2504 individuals across 26 different populations. Furthermore, webGQT supports filtering single-nucleotide variants, short insertions/deletions, copy number or any other variant genotypes supported by the VCF specification. Our results show that webGQT can be used as an online web service, or deployed on personal computers or local servers within research groups. AVAILABILITY: webGQT is made available to the users in three forms: 1) as a webserver available at https://vm1138.kaj.pouta.csc.fi/webgqt/, 2) as an R package to install on personal computers, and 3) as part of the same R package to configure on the user's own servers. The application is available for installation at https://github.com/arumds/webgqt.

A putative silencer variant in a spontaneous canine model of retinitis pigmentosa.

Retinitis pigmentosa (RP) is the leading cause of blindness with nearly two million people affected worldwide. Many genes have been implicated in RP, yet in 30-80% of the RP patients the genetic cause remains unknown. A similar phenotype, progressive retinal atrophy (PRA), affects many dog breeds including the Miniature Schnauzer. We performed clinical, genetic and functional experiments to identify the genetic cause of PRA in the breed. The age of onset and pattern of disease progression suggested that at least two forms of PRA, types 1 and 2 respectively, affect the breed, which was confirmed by genome-wide association study that implicated two distinct genomic loci in chromosomes 15 and X, respectively. Whole-genome sequencing revealed a fully segregating recessive regulatory variant in type 1 PRA. The associated variant has a very recent origin based on haplotype analysis and lies within a regulatory site with the predicted binding site of HAND1::TCF3 transcription factor complex. Luciferase assays suggested that mutated regulatory sequence increases expression. Case-control retinal expression comparison of six best HAND1::TCF3 target genes were analyzed with quantitative reverse-transcriptase PCR assay and indicated overexpression of EDN2 and COL9A2 in the affected retina. Defects in both EDN2 and COL9A2 have been previously associated with retinal degeneration. In summary, our study describes two genetically different forms of PRA and identifies a fully penetrant variant in type 1 form with a possible regulatory effect. This would be among the first reports of a regulatory variant in retinal degeneration in any species, and establishes a new spontaneous dog model to improve our understanding of retinal biology and gene regulation while the affected breed will benefit from a reliable genetic testing.

Recessive missense LAMP3 variant associated with defect in lamellar body biogenesis and fatal neonatal interstitial lung disease in dogs.

Neonatal interstitial lung diseases due to abnormal surfactant biogenesis are rare in humans and have never been reported as a spontaneous disorder in animals. We describe here a novel lung disorder in Airedale Terrier (AT) dogs with clinical symptoms and pathology similar to the most severe neonatal forms of human surfactant deficiency. Lethal hypoxic respiratory distress and failure occurred within the first days or weeks of life in the affected puppies. Transmission electron microscopy of the affected lungs revealed maturation arrest in the formation of lamellar bodies (LBs) in the alveolar epithelial type II (AECII) cells. The secretory organelles were small and contained fewer lamellae, often in combination with small vesicles surrounded by an occasionally disrupted common limiting membrane. A combined approach of genome-wide association study and whole exome sequencing identified a recessive variant, c.1159G>A, p.(E387K), in LAMP3, a limiting membrane protein of the cytoplasmic surfactant organelles in AECII cells. The substitution resides in the LAMP domain adjacent to a conserved disulfide bond. In summary, this study describes a novel interstitial lung disease in dogs, identifies a new candidate gene for human surfactant dysfunction and brings important insights into the essential role of LAMP3 in the process of the LB formation.

Improved Criteria for the Classification of Titin Variants in Inherited Skeletal Myopathies.

BACKGROUND: Extensive genetic screening results in the identification of thousands of rare variants that are difficult to interpret. Because of its sheer size, rare variants in the titin gene (TTN) are detected frequently in any individual. Unambiguous interpretation of molecular findings is almost impossible in many patients with myopathies or cardiomyopathies. OBJECTIVE: To refine the current classification framework for TTN-associated skeletal muscle disorders and standardize the interpretation of TTN variants. METHODS: We used the guidelines issued by the American College of Medical Genetics and Genomics (ACMG) and the Association for Molecular Pathology (AMP) to re-analyze TTN genetic findings from our patient cohort. RESULTS: We identified in the classification guidelines three rules that are not applicable to titin-related skeletal muscle disorders; six rules that require disease-/gene-specific adjustments and four rules requiring quantitative thresholds for a proper use. In three cases, the rule strength need to be modified. CONCLUSIONS: We suggest adjustments are made to the guidelines. We provide frequency thresholds to facilitate filtering of candidate causative variants and guidance for the use and interpretation of functional data and co-segregation evidence. We expect that the variant classification framework for TTN-related skeletal muscle disorders will be further improved along with a better understanding of these diseases.

Improving Copy Number Variant Detection from Sequencing Data with a Combination of Programs and a Predictive Model.

Bioinformatics tools for analyzing copy number variants (CNVs) from massively parallel sequencing (MPS) data are less well developed compared with other variant types. We present an efficient bioinformatics pipeline for CNV detection from gene panel MPS data in neuromuscular disorders. CNVs were generated in silico into samples sequenced with a previously published MPS gene panel. The in silico CNVs from these samples were analyzed with four programs having complementary CNV detection ranges: CoNIFER, XHMM, ExomeDepth, and CODEX. A logistic regression model was trained with the obtained set of in silico CNV detections to predict true-positive CNV detections among all CNV detections from samples. This model was validated using 66 control samples with a verified true-positive (n = 58) or false-positive (n = 8) CNV detection. Applying all four programs together provided more sensitive detection results with in silico CNVs than other program combinations or any program alone. Furthermore, a model with CNV detection-specific scores from all four programs as variables performed overall best in the validation. No single program could detect all CNV sizes and types equally or with enough accuracy. Therefore, a combination of carefully selected programs should be used to maximize detection accuracy. In addition, the detected CNVs should be reviewed with a statistical model to streamline and standardize the filtering of the detections for annotation.

Novel protective and risk loci in hip dysplasia in German Shepherds.

Canine hip dysplasia is a common, non-congenital, complex and hereditary disorder. It can inflict severe pain via secondary osteoarthritis and lead to euthanasia. An analogous disorder exists in humans. The genetic background of hip dysplasia in both species has remained ambiguous despite rigorous studies. We aimed to investigate the genetic causes of this disorder in one of the high-risk breeds, the German Shepherd. We performed genetic analyses with carefully phenotyped case-control cohorts comprising 525 German Shepherds. In our genome-wide association studies we identified four suggestive loci on chromosomes 1 and 9. Targeted resequencing of the two loci on chromosome 9 from 24 affected and 24 control German Shepherds revealed deletions of variable sizes in a putative enhancer element of the NOG gene. NOG encodes for noggin, a well-described bone morphogenetic protein inhibitor affecting multiple developmental processes, including joint development. The deletion was associated with the healthy controls and mildly dysplastic dogs suggesting a protective role against canine hip dysplasia. Two enhancer variants displayed a decreased activity in a dual luciferase reporter assay. Our study identifies novel loci and candidate genes for canine hip dysplasia, with potential regulatory variants in the NOG gene. Further research is warranted to elucidate how the identified variants affect the expression of noggin in canine hips, and what the potential effects of the other identified loci are.

Canine models of human amelogenesis imperfecta: identification of novel recessive ENAM and ACP4 variants.

Amelogenesis imperfecta (AI) refers to a genetically and clinically heterogeneous group of inherited disorders affecting the structure, composition, and quantity of tooth enamel. Both non-syndromic and syndromic forms of AI have been described and several genes affecting various aspects of the enamel physiology have been reported. Genetically modified murine models of various genes have provided insights into the complex regulation of proper amelogenesis. Non-syndromic AI occurs spontaneously also in dogs with known recessive variants in ENAM and SLC24A4 genes. Unlike rodents with a reduced dentition and continuously erupting incisors, canine models are valuable for human AI due to similarity in the dental anatomy including deciduous and permanent teeth. We have performed a series of clinical and genetic analyses to investigate AI in several breeds of dogs and describe here two novel recessive variants in the ENAM and ACP4 genes. A fully segregating missense variant (c.716C>T) in exon 8 of ENAM substitutes a well-conserved proline to leucine, p.(Pro239Leu), resulting in a clinical hypomineralization of teeth. A 1-bp insertion in ACP4 (c.1189dupG) is predicted to lead to a frameshift, p.(Ala397Glyfs), resulting in an abnormal C-terminal part of the protein, and hypoplastic AI. The ENAM variant was specific for Parson Russell Terriers with a carrier frequency of 9%. The ACP4 variant was found in two breeds, Akita and American Akita with a carrier frequency of 22%. These genetic findings establish novel canine models of human AI with a particular interest in the case of the ACP4-deficient model, since ACP4 physiology is poorly characterized in human AI. The affected dogs could also serve as preclinical models for novel treatments while the breeds would benefit from genetic tests devised here for veterinary diagnostics and breeding programs.

A homozygous missense variant in the alkaline phosphatase gene ALPL is associated with a severe form of canine hypophosphatasia.

Inherited skeletal disorders affect both humans and animals. In the current study, we have performed series of clinical, pathological and genetic examinations to characterize a previously unreported skeletal disease in the Karelian Bear Dog (KBD) breed. The disease was recognized in seven KBD puppies with a variable presentation of skeletal hypomineralization, growth retardation, seizures and movement difficulties. Exome sequencing of one affected dog revealed a homozygous missense variant (c.1301T > G; p.V434G) in the tissue non-specific alkaline phosphatase gene, ALPL. The identified recessive variant showed full segregation with the disease in a cohort of 509 KBDs with a carrier frequency of 0.17 and was absent from 303 dogs from control breeds. In humans, recessive and dominant ALPL mutations cause hypophosphatasia (HPP), a metabolic bone disease with highly heterogeneous clinical manifestations, ranging from lethal perinatal hypomineralization to a relatively mild dental disease. Our study reports the first naturally occurring HPP in animals, resembling the human infantile form. The canine HPP model may serve as a preclinical model while a genetic test will assist in breeding programs.

Author Correction: Assembly and Analysis of Unmapped Genome Sequence Reads Reveal Novel Sequence and Variation in Dogs.

A correction to this article has been published and is linked from the HTML and PDF versions of this paper. The error has not been fixed in the paper.

Assembly and Analysis of Unmapped Genome Sequence Reads Reveal Novel Sequence and Variation in Dogs.

Dogs are excellent animal models for human disease. They have extensive veterinary histories, pedigrees, and a unique genetic system due to breeding practices. Despite these advantages, one factor limiting their usefulness is the canine genome reference (CGR) which was assembled using a single purebred Boxer. Although a common practice, this results in many high-quality reads remaining unmapped. To address this whole-genome sequence data from three breeds, Border Collie (n = 26), Bearded Collie (n = 7), and Entlebucher Sennenhund (n = 8), were analyzed to identify novel, non-CGR genomic contigs using the previously validated pseudo-de novo assembly pipeline. We identified 256,957 novel contigs and paired-end relationships together with BLAT scores provided 126,555 (49%) high-quality contigs with genomic coordinates containing 4.6 Mb of novel sequence absent from the CGR. These contigs close 12,503 known gaps, including 2.4 Mb containing partially missing sequences for 11.5% of Ensembl, 16.4% of RefSeq and 12.2% of canFam3.1+ CGR annotated genes and 1,748 unmapped contigs containing 2,366 novel gene variants. Examples for six disease-associated genes (SCARF2, RD3, COL9A3, FAM161A, RASGRP1 and DLX6) containing gaps or alternate splice variants missing from the CGR are also presented. These findings from non-reference breeds support the need for improvement of the current Boxer-only CGR to avoid missing important biological information. The inclusion of the missing gene sequences into the CGR will facilitate identification of putative disease mutations across diverse breeds and phenotypes.

Nonsense variant in COL7A1 causes recessive dystrophic epidermolysis bullosa in Central Asian Shepherd dogs.

A rare hereditary mechanobullous disorder called epidermolysis bullosa (EB) causes blistering in the skin and the mucosal membranes. To date, nineteen EB-related genes have been discovered in human and other species. We describe here a novel EB variant in dogs. Two newborn littermates of Central Asian Shepherd dogs with severe signs of skin blistering were brought to a veterinary clinic and euthanized due to poor prognosis. In post-mortem examination, the puppies were shown to have findings in the skin and the mucosal membranes characteristic of EB. A whole-genome sequencing of one of the affected puppies was performed to identify the genetic cause. The resequencing data were filtered under a recessive model against variants from 31 other dog genomes, revealing a homozygous case-specific nonsense variant in one of the known EB-causing genes, COL7A1 (c.4579C>T, p.R1527*). The variant results in a premature stop codon and likely absence of the functional protein in the basement membrane of the skin in the affected dogs. This was confirmed by immunohistochemistry using a COL7A1 antibody. Additional screening of the variant indicated full penetrance and breed specificity at ~28% carrier frequency. In summary, this study reveals a novel COL7A1 variant causing recessive dystrophic EB and provides a genetic test for the eradication of the disease from the breed.