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

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

Association of missense variants in GDF9 with litter size in Entlebucher Mountain dogs.

In the past two decades, average litter size (ALS) in Entlebucher Mountain dogs decreased by approximately 0.8 puppies. We conducted a GWAS for ALS using the single-step methodology to take advantage of 1632 pedigree records, 892 phenotypes and 372 genotypes (173 662 markers) for which only 12% of the dogs had both phenotypes and genotypes available. Our analysis revealed associations towards the growth differentiation factor 9 gene (GDF9), which is known to regulate oocyte maturation. The trait heritability was estimated at 43.1%, from which approximately 15% was accountable by the GDF9 locus alone. Therefore, markers flanking GDF9 explained approximately 6.5% of the variance in ALS. Analysis of WGSs revealed two missense substitutions in GDF9, one of which (g.11:21147009G>A) affected a highly conserved nucleotide in vertebrates. The derived allele A was validated in 111 dogs and shown to be associated with decreased ALS (-0.75 ± 0.22 puppies per litter). The variant was further predicted to cause a proline to serine substitution. The affected residue was immediately followed by a six-residue deletion that is fixed in the canine species but absent in non-canids. We further confirmed that the deletion is prevalent in the Canidae family by sequencing three species of wild canids. Since canids uniquely ovulate oocytes at the prophase stage of the first meiotic division, requiring maturation in the oviduct, we conjecture that the amino acid substitution and the six-residue deletion of GDF9 may serve as a model for insights into the dynamics of oocyte maturation in canids.

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

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

A missense variant in the NSDHL gene in a Chihuahua with a congenital cornification disorder resembling inflammatory linear verrucous epidermal nevi.

Congenital hemidysplasia with ichthyosiform nevus and limb defects syndrome in humans is a genodermatosis characterized by inflammatory linear verrucous epidermal nevi (ILVEN), often showing a striking lateralization pattern. It is caused by variants in the NSDHL gene encoding a 3ß-hydroxysteroid dehydrogenase involved in the cholesterol biosynthesis pathway. In the present study, we investigated a female Chihuahua, which showed clinical and histological signs of ILVEN. We performed a candidate gene analysis in the affected animal. This analysis revealed a single missense variant in the NSDHL gene in the affected dog (XM_014111859.2:c.700G>A). The variant is predicted to cause a non-conservative amino acid change from glycine to arginine, XP_013967334.1:p.(Gly234Arg). The mutant allele was absent from WGS data of 594 genetically diverse dogs and eight wolves. Sanger sequencing confirmed that the variant was heterozygous in the affected dog and absent from 22 control Chihuahuas. Based on the knowledge about the functional impact of NSDHL variants in dogs and other species, c.700G>A is probably pathogenic and a convincing candidate causative variant for the observed skin lesions in the affected Chihuahua.

An ABCA12 missense variant in a Shorthorn calf with ichthyosis fetalis.

Two clinical forms of ichthyosis in cattle have been reported, ichthyosis fetalis and congenital ichthyosis. Ichthyosis poses animal welfare and economic issues and the more severe form, ichthyosis fetalis, is lethal. A Shorthorn calf with ichthyosis fetalis was investigated and a likely causal missense variant on chromosome 2 in the ABCA12 gene (NM_001191294.2:c.6776T>C) was identified by whole genome sequencing. Mutations in the ABCA12 gene are known to cause ichthyosis fetalis in cattle and Harlequin ichthyosis in humans. Sanger sequencing of the affected calf and the dam confirmed the variant was homozygous in the affected calf and heterozygous in the dam. Further genotyping of 130 Shorthorn animals from the same property revealed an estimated allele frequency of 3.8%. The presented findings enable genetic testing for breeding and diagnostics.

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

The domestic dog serves as an excellent model to investigate the genetic basis of disease. More than 400 heritable traits analogous to human diseases have been described in dogs. To further canine medical genetics research, we established the Dog Biomedical Variant Database Consortium (DBVDC) and present a comprehensive list of functionally annotated genome variants that were identified with whole genome sequencing of 582 dogs from 126 breeds and eight wolves. The genomes used in the study have a minimum coverage of 10× and an average coverage of ~24×. In total, we identified 23 133 692 single-nucleotide variants (SNVs) and 10 048 038 short indels, including 93% undescribed variants. On average, each individual dog genome carried ∼4.1 million single-nucleotide and ~1.4 million short-indel variants with respect to the reference genome assembly. About 2% of the variants were located in coding regions of annotated genes and loci. Variant effect classification showed 247 141 SNVs and 99 562 short indels having moderate or high impact on 11 267 protein-coding genes. On average, each genome contained heterozygous loss-of-function variants in 30 potentially embryonic lethal genes and 97 genes associated with developmental disorders. More than 50 inherited disorders and traits have been unravelled using the DBVDC variant catalogue, enabling genetic testing for breeding and diagnostics. This resource of annotated variants and their corresponding genotype frequencies constitutes a highly useful tool for the identification of potential variants causative for rare inherited disorders in dogs.

A novel KIT deletion variant in a German Riding Pony with white-spotting coat colour phenotype.

White spotting phenotypes in horses may be caused by developmental alterations impairing melanoblast differentiation, survival, migration and/or proliferation. Candidate genes for white-spotting phenotypes in horses include EDNRB, KIT, MITF, PAX3 and TRPM1. We investigated a German Riding Pony with a sabino-like phenotype involving extensive white spots on the body together with large white markings on the head and almost completely white legs. We obtained whole genome sequence data from this horse. The analysis revealed a heterozygous 1273-bp deletion spanning parts of intron 2 and exon 3 of the equine KIT gene (Chr3: 79 579 925-79 581 197). We confirmed the breakpoints of the deletion by PCR and Sanger sequencing. Knowledge of the functional impact of similar KIT variants in horses and other species suggests that this deletion represents a plausible candidate causative variant for the white-spotting phenotype. We propose the designation W28 for the mutant allele.

Compound heterozygosity for TNXB genetic variants in a mixed-breed dog with Ehlers-Danlos syndrome.

The Ehlers-Danlos syndromes (EDSs) are a heterogeneous group of inherited connective tissue disorders characterized by skin hyperextensibility, joint hypermobility and tissue fragility. Inherited disorders similar to human EDS have been reported in different mammalian species. In the present study, we investigated a female mixed-breed dog with clinical signs of EDS. Whole-genome sequencing of the affected dog revealed two missense variants in the TNXB gene, encoding the extracellular matrix protein tenascin XB. In humans, TNXB genetic variants cause classical-like EDS or the milder hypermobile EDS. The affected dog was heterozygous at both identified variants. Each variant allele was transmitted from one of the case's parents, consistent with compound heterozygosity. Although one of the variant alleles, XM_003431680.3:c.2012G>A, p.(Ser671Asn), was private to the family of the affected dog and absent from whole-genome sequencing data of 599 control dogs, the second variant allele, XM_003431680.3:c.2900G>A, p.(Gly967Asp), is present at a low frequency in the Chihuahua and Poodle population. Given that TNXB is a functional candidate gene for EDS, we suggest that compound heterozygosity for the identified TNXB variants may have caused the EDS-like phenotype in the affected dog. Chihuahuas and Poodles should be monitored for EDS cases, which might confirm the hypothesized pathogenic effect of the segregating TNXB variant.

A complex structural variant at the KIT locus in cattle with the Pinzgauer spotting pattern.

A specific white spotting phenotype, termed finching or line-backed spotting, is known for all Pinzgauer cattle and occurs occasionally in Tux-Zillertaler cattle, two Austrian breeds. The so-called Pinzgauer spotting is inherited as an autosomal incompletely dominant trait. A genome-wide association study using 27 white spotted and 16 solid-coloured Tux-Zillertaler cattle, based on 777k SNP data, revealed a strong signal on chromosome 6 at the KIT locus. Haplotype analyses defined a critical interval of 122 kb downstream of the KIT coding region. Whole-genome sequencing of a Pinzgauer cattle and comparison to 338 control genomes revealed a complex structural variant consisting of a 9.4-kb deletion and an inversely inserted duplication of 1.5 kb fused to a 310-kb duplicated segment from chromosome 4. A diagnostic PCR was developed for straightforward genotyping of carriers for this structural variant (KITPINZ ) and confirmed that the variant allele was present in all Pinzgauer and most of the white spotted Tux-Zillertaler cattle. In addition, we detected the variant in all Slovenian Cika, British Gloucester and Spanish Berrenda en negro cattle with similar spotting patterns. Interestingly, the KITPINZ variant occurs in some white spotted animals of the Swiss breeds Evolèner and Eringer. The introgression of the KITPINZ variant confirms admixture and the reported historical relationship of these short-headed breeds with Austrian Tux-Zillertaler and suggests a mutation event, occurring before breed formation.

In silico and in vitro analysis of genetic variants of the equine CYP3A94, CYP3A95 and CYP3A97 isoenzymes.

Cytochrome P450 enzymes (CYPs) of the equine CYP3A subfamily are predominantly involved in drug metabolism. In this study, genetic variants of the equine CYP3A94, CYP3A95, and CYP3A97 were identified and characterized using in silico modeling and in vitro enzyme kinetics. The genomes of 81 horses were sequenced to obtain the genetic variants. Structural CYP modifications of the most frequent variants were analyzed in silico using the 3D-structures predicted by homology modeling. Enzyme kinetic analyses were performed using testosterone as substrate. Twenty genetic variants were found including five missense variants (CYP3A94:p.Asp217Asn, CYP3A95:p.Asp214His, CYP3A95:p.Ser392Thr, CYP3A97:p.Ile119Thr, CYP3A97:p.Met500Val) with a higher percentage of minor allele frequency (MAF) (range 0.2-0.4). A splice-site variant (c.798 + 1G > A) in CYP3A94, likely to generate a truncated protein, was found in 50% of the horses. CYP3A94:p.Asp217Asn and CYP3A95:p.Asp214His were localized on the CYP F-α-helix, an important region for the substrate interactions in the human CYP3A4. Testosterone 2ß-hydroxylation was diminished in CYP3A94217Asn and CYP3A95392Thr. Ketoconazole inhibited 2ß-hydroxylation differently in the five variants with the most pronounced inhibition obtained for CYP3A95392Thr. In vitro and in silico analyses of genetic variants allow unraveling structural features in equine CYPs that correlate with changes in the CYP activity.

Whole-genome sequencing reveals a large deletion in the MITF gene in horses with white spotted coat colour and increased risk of deafness.

White spotting phenotypes in horses are highly valued in some breeds. They are quite variable and may range from the common white markings up to completely white horses. EDNRB, KIT, MITF, PAX3 and TRPM1 represent known candidate genes for white spotting phenotypes in horses. For the present study, we investigated an American Paint Horse family segregating a phenotype involving white spotting and blue eyes. Six of eight horses with the white-spotting phenotype were deaf. We obtained whole-genome sequence data from an affected horse and specifically searched for structural variants in the known candidate genes. This analysis revealed a heterozygous ~63-kb deletion spanning exons 6-9 of the MITF gene (chr16:21 503 211-21 566 617). We confirmed the breakpoints of the deletion by PCR and Sanger sequencing. PCR-based genotyping revealed that all eight available affected horses from the family carried the deletion. The finding of an MITF variant fits well with the syndromic phenotype involving both depigmentation and an increased risk for deafness and corresponds to human Waardenburg syndrome type 2A. Our findings will enable more precise genetic testing for depigmentation phenotypes in horses.

A second KRT71 allele in curly coated dogs.

Major characteristics of coat variation in dogs can be explained by variants in only a few genes. Until now, only one missense variant in the KRT71 gene, p.Arg151Trp, has been reported to cause curly hair in dogs. However, this variant does not explain the curly coat in all breeds as the mutant 151 Trp allele, for example, is absent in Curly Coated Retrievers. We sequenced the genome of a Curly Coated Retriever at 22× coverage and searched for variants in the KRT71 gene. Only one protein-changing variant was present in a homozygous state in the Curly Coated Retriever and absent or present in a heterozygous state in 221 control dogs from different dog breeds. This variant, NM_001197029.1:c.1266_1273delinsACA, was an indel variant in exon 7 that caused a frameshift and an altered and probably extended C-terminus of the KRT71 protein NP_001183958.1:p.(Ser422ArgfsTer?). Using Sanger sequencing, we found that the variant was fixed in a cohort of 125 Curly Coated Retrievers and segregating in five of 14 additionally tested breeds with a curly or wavy coat. KRT71 variants cause curly hair in humans, mice, rats, cats and dogs. Specific KRT71 variants were further shown to cause alopecia. Based on this knowledge from other species and the predicted molecular consequence of the newly identified canine KRT71 variant, it is a compelling candidate causing a second curly hair allele in dogs. It might cause a slightly different coat phenotype than the previously published p.Arg151Trp variant and could potentially be associated with follicular dysplasia in dogs.

Comprehensive characterization of horse genome variation by whole-genome sequencing of 88 horses.

Whole-genome sequencing studies are vital to gain a thorough understanding of genomic variation. Here, we summarize the results of a whole-genome sequencing study comprising 88 horses and ponies from diverse breeds at 19.1× average coverage. The paired-end reads were mapped to the current EquCab3.0 horse reference genome assembly, and we identified approximately 23.5 million single nucleotide variants and 2.3 million short indel variants. Our dataset included at least 7 million variants that were not previously reported. On average, each individual horse genome carried ∼5.7 million single nucleotides and 0.8 million small indel variants with respect to the reference genome assembly. The variants were functionally annotated. We provide two examples for potentially deleterious recessive alleles that were identified in a heterozygous state in individual genome sequences. Appropriate management of such deleterious recessive alleles in horse breeding programs should help to improve fertility and reduce the prevalence of heritable diseases. This comprehensive dataset has been made publicly available, will represent a valuable resource for future horse genetic studies and supports the goal of accelerating the rates of genetic gain in domestic horse.

A non-coding regulatory variant in the 5'-region of the MITF gene is associated with white-spotted coat in Brown Swiss cattle.

Recently, the Swiss breeding association reported an increasing number of white-spotted cattle in the Brown Swiss breed, which is normally solid brown coloured. A total of 60 Brown Swiss cattle with variably sized white abdominal spots, facial markings and depigmented claws were collected for this study. A genome-wide association study using 40k SNP genotypes of 20 cases and 1619 controls enabled us to identify an associated genome region on chromosome 22 containing the MITF gene, encoding the melanogenesis associated transcription factor. Variants at the MITF locus have been reported before to be associated with white or white-spotted phenotypes in other species such as horses, dogs and mice. Whole-genome sequencing of a single white-spotted cow and subsequent genotyping of 172 Brown Swiss cattle revealed two significantly associated completely linked single nucleotide variants (rs722765315 and rs719139527). Both variants are located in the 5'-regulatory region of the bovine MITF gene, and comparative sequence analysis showed that the variant rs722765315, located 139 kb upstream of the transcription start site of the bovine melanocyte-specific MITF transcript, is situated in a multi-species conserved sequence element which is supposed to be regulatory important. Therefore, we hypothesize that rs722765315 represents the most likely causative variant for the white-spotting phenotype observed in Brown Swiss cattle. Presence of the mutant allele in a heterozygous or homozygous state supports a dominant mode of inheritance with incomplete penetrance and results in a variable extent of coat colour depigmentation.

A frameshift variant in the EDA gene in Dachshunds with X-linked hypohidrotic ectodermal dysplasia.

X-linked hypohidrotic ectodermal dysplasia (XLHED) is a genetic disease characterized by hypoplasia or absence of hair, teeth and sweat glands. The EDA gene, located on the X chromosome, encodes the type II transmembrane protein ectodysplasin A. Variants in the EDA gene can lead to XLHED in humans, mice, cattle and dogs. In the present study, we investigated a litter of Dachshund puppies, of which four male puppies showed clinical signs of XLHED. We performed a candidate gene analysis in one affected puppy and several non-affected relatives. This analysis revealed a single base-pair deletion in the coding sequence of the EDA gene in the affected puppy (NM_001014770.2:c.842delT). The deletion is predicted to cause a frameshift, NP_001014770.1:p.(Leu281HisfsTer22), leading to a premature stop codon which truncates more than one quarter of the EDA protein. Sanger sequencing results confirmed that this variant was inherited from the dam. Based on knowledge about the functional impact of EDA variants in dogs and other species, c.842delT is a convincing candidate causative variant for the observed XLHED in the male puppies.

Genome-wide association study and heritability estimate for ectopic ureters in Entlebucher mountain dogs.

An ectopic ureter is a congenital anomaly which may lead to urinary incontinence and without a surgical intervention even to end-stage kidney disease. A genetic component contributes to the development of this anomaly in Entlebucher mountain dogs (EMD); however, its nature remains unclear. Using the Illumina CanineHD bead chip, a case-control genome-wide association study was performed to identify SNPs associated with the trait. Six loci on canine chromosomes 3, 17, 27 and 30 were identified with 16 significantly associated SNPs. There was no single outstanding SNP associated with the phenotype, and the association signals were not close to known genes involved in human congenital anomalies of the kidney or lower urinary tract. Additional research will be necessary to elucidate the potential role of the associated genes in the development of ectopic ureters in the EMD breed.

Generation of an equine biobank to be used for Functional Annotation of Animal Genomes project.

The Functional Annotation of Animal Genomes (FAANG) project aims to identify genomic regulatory elements in both sexes across multiple stages of development in domesticated animals. This study represents the first stage of the FAANG project for the horse, Equus caballus. A biobank of 80 tissue samples, two cell lines and six body fluids was created from two adult Thoroughbred mares. Ante-mortem assessments included full physical examinations, lameness, ophthalmologic and neurologic evaluations. Complete blood counts and serum biochemistries were also performed. At necropsy, in addition to tissue samples, aliquots of serum, ethylenediaminetetraacetic acid (EDTA) plasma, heparinized plasma, cerebrospinal fluid, synovial fluid, urine and microbiome samples from all regions of the gastrointestinal and urogenital tracts were collected. Epidermal keratinocytes and dermal fibroblasts were cultured from skin samples. All tissues were grossly and histologically evaluated by a board-certified veterinary pathologist. The results of the clinical and pathological evaluations identified subclinical eosinophilic and lymphocytic infiltration throughout the length of the gastrointestinal tract as well as a mild clinical lameness in both animals. Each sample was cryo-preserved in multiple ways, and nuclei were extracted from selected tissues. These samples represent the first published systemically healthy equine-specific biobank with extensive clinical phenotyping ante- and post-mortem. The tissues in the biobank are intended for community-wide use in the functional annotation of the equine genome. The use of the biobank will improve the quality of the reference annotation and allow all equine researchers to elucidate unknown genomic and epigenomic causes of disease.

A frameshift variant in the COL5A1 gene in a cat with Ehlers-Danlos syndrome.

Ehlers-Danlos syndrome (EDS) is a group of heritable connective tissue disorders caused by defective collagen synthesis or incorrect assembly of the collagen triple helical structure. EDS is characterised by joint hypermobility, skin hyperextensibility, abnormal scarring, poor wound healing and tissue friability. Human EDS may be caused by variants in several different genes including COL5A1, which encodes the collagen type V alpha 1 chain. For the present study we investigated a 1.5-year-old, spayed female, domestic shorthair cat with EDS. The affected cat showed multiple recurrent skin tears, hyperextensibility of the skin and joint abnormalities. We obtained whole genome sequencing data from the affected cat and searched for variants in candidate genes known to cause EDS. We detected a heterozygous single base-pair deletion in exon 43 of the COL5A1 gene, namely c.3420delG. The deletion was predicted to result in a frameshift and premature stop codon: p.(Leu1141SerfsTer134). Sanger sequencing confirmed that the variant was present in the affected cat and absent from 103 unaffected cats from different breeds. The variant was also absent from a Burmese cat with EDS. Based on knowledge about the functional impact of COL5A1 variants in other species, COL5A1:c.3420delG represents a compelling candidate causative variant for the observed EDS in the affected cat.