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.

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.

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

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

SOAT1 missense variant in two cats with sebaceous gland dysplasia

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

Heterozygous ATP2A2 missense variant identified in a Shih Tzu with Darier disease

Darier disease is caused by heterozygous loss of function variants in the ATP2A2 gene encoding the endoplasmic/sarcoplasmic reticulum Ca²⁺ pump ATP2A2. Defective intracellular calcium signaling in the epidermis results in a loss of desmosomal adhesion and the development of characteristic skin lesions. In this study, we investigated a Shih Tzu that developed erythematous papules on the ventrum and, over time, the dorsal neck and a nodule in the right ear canal with secondary ear infection. Histopathologic examination demonstrated discrete foci of acantholysis affecting suprabasal layers of the epidermis. Whole genome sequencing of the affected dog identified a heterozygous missense variant, p.N809H, affecting an evolutionarily conserved amino acid residue of the ATP2A2 protein. The highly characteristic clinical and histopathologic findings together with a plausible variant in the only known functional candidate gene establish the diagnosis of canine Darier disease in the studied dog and highlight the potential of genetic analyses as complementary diagnostic approach in veterinary medicine.

KRT5 missense variant in a Cardigan Welsh Corgi with epidermolysis bullosa simplex

Epidermolysis bullosa (EB) is a group of blistering disorders that includes several subtypes, classified according to their level of cleavage. Typical clinical signs are blisters and erosions resulting from minimal trauma. The disease has been described in many mammalian species and pathogenic variants in at least 18 different genes have been identified. In the present study, we investigated a Cardigan Welsh Corgi with congenital clinical signs consistent with epidermolysis bullosa. The puppy had blisters and erosions on the paw pads, and the oral mucosa. Histologic examination demonstrated the typical clefting between the dermis and epidermis and confirmed the clinical suspicion. We obtained whole genome sequencing data from the affected puppy and searched for variants in candidate genes known to cause EB. This revealed a heterozygous missense variant, KRT5:p.(E476K), affecting the highly conserved KLLEGE motif of keratin 5. The mutant allele in the affected puppy arose owing to a de novo mutation event as it was absent from both unaffected parents. Knowledge of the functional impact of KRT5 variants in other species together with the demonstration of the de novo mutation event establishes KRT5:p.(E476K) as causative variant for the observed EBS.

A de novo variant in the keratin 1 gene (KRT1) in a Chinese shar-pei dog with severe congenital cornification disorder and non-epidermolytic ichthyosis

A 3-months old Chinese shar-pei puppy with ichthyosis was investigated. The dog showed generalized scaling, alopecia and footpad lesions. Histopathological examinations demonstrated a non-epidermolytic hyperkeratosis. The parents of the affected puppy did not show any skin lesions. A trio whole genome sequencing analysis identified a heterozygous de novo 3 bp deletion in the KRT1 gene in the affected dog. This variant, NM_001003392.1:c.567_569del, is predicted to delete a single asparagine from the conserved coil 1A motif within the rod domain of KRT1, NP_001003392.1:p.(Asn190del). Immunohistochemistry demonstrated normal levels of KRT1 expression in the epidermis and follicular epithelia. This might indicate that the variant possibly interferes with keratin dimerization or another function of KRT1. Missense variants affecting the homologous asparagine residue of the human KRT1 cause epidermolytic hyperkeratosis. Histologically, the investigated Chinese shar-pei showed a non-epidermolytic ichthyosis. The finding of a de novo variant in an excellent functional candidate gene strongly suggests that KRT1:p.Asn190del caused the ichthyosis phenotype in the affected Chinese shar-pei. To the best of our knowledge, this is the first description of a KRT1-related non-epidermolytic ichthyosis in domestic animals.

A COL5A2 In-Frame Deletion in a Chihuahua with Ehlers-Danlos Syndrome

Ehlers-Danlos syndrome (EDS) is a group of heterogeneous, rare diseases affecting the connective tissues. The main clinical signs of EDS are skin hyperextensibility, joint hypermobility, and skin fragility. Currently, the classification of EDS in humans distinguishes 13 clinical subtypes associated with variants in 20 different genes, reflecting the heterogeneity of this set of diseases. At present, variants in three of these genes have also been identified in dogs affected by EDS. The purpose of this study was to characterize the clinical and histopathological phenotype of an EDS-affected Chihuahua and to identify the causative genetic variant for the disease. The clinical examination suggested a diagnosis of classical EDS. Skin histopathology revealed an abnormally thin dermis, which is compatible with classical EDS. Whole-genome sequencing identified a heterozygous de novo 27 bp deletion in the COL5A2 gene, COL5A2:c.3388_3414del. The in-frame deletion is predicted to remove 9 amino acids in the triple-helical region of COL5A2. The molecular analysis and identification of a likely pathogenic variant in COL5A2 confirmed the subtype as a form of classical EDS. This is the first report of a COL5A2-related EDS in a dog.

Independent COL5A1 Variants in Cats with Ehlers-Danlos Syndrome

We investigated four cats with similar clinical skin-related signs strongly suggestive of Ehlers-Danlos syndrome (EDS). Cases no. 1 and 4 were unrelated and the remaining two cases, no. 2 and 3, were reportedly siblings. Histopathological changes were characterized by severely altered dermal collagen fibers. Transmission electron microscopy in one case demonstrated abnormalities in the collagen fibril organization and structure. The genomes of the two unrelated affected cats and one of the affected siblings were sequenced and individually compared to 54 feline control genomes. We searched for private protein changing variants in known human EDS candidate genes and identified three independent heterozygous COL5A1 variants. COL5A1 is a well-characterized candidate gene for classical EDS. It encodes the proα1 chain of type V collagen, which is needed for correct collagen fibril formation and the integrity of the skin. The identified variants in COL5A1 are c.112_118+15del or r.spl?, c.3514A>T or p.(Lys1172*), and c.3066del or p.(Gly1023Valfs*50) for cases no. 1, 2&3, and 4, respectively. They presumably all lead to nonsense-mediated mRNA decay, which results in haploinsufficiency of COL5A1 and causes the alterations of the connective tissue. The whole genome sequencing approach used in this study enables a refinement of the diagnosis for the affected cats as classical EDS. It further illustrates the potential of such experiments as a precision medicine approach in animals with inherited diseases.

ABHD5 frameshift deletion in Golden Retrievers with ichthyosis

Ichthyoses are hereditary skin disorders characterized by the formation of scales and defects in the outermost layer of the epidermis. In dogs, at least six different breed-specific ichthyoses including a relatively common PNPLA1-related autosomal recessive ichthyosis in Golden Retrievers are known. In this study, we investigated 14 Golden Retrievers with scales that were not homozygous for the mutant PNPLA1 allele suggesting a genetically distinct new form of ichthyosis. Histopathological examinations showed lamellar, orthokeratotic hyperkeratosis, and mildly hyperplastic epidermis that led to the diagnosis of a nonepidermolytic ichthyosis. Combined linkage and homozygosity mapping in 14 cases and 30 nonaffected family members delimited a critical interval of ∼12.7 Mb on chromosome 23. Whole-genome sequencing of an affected dog revealed a single protein-changing variant within this region that was not present in 795 control genomes. The identified variant is a 14 bp deletion in the ABHD5 gene (c.1006_1019del), leading to a frameshift and altering the last 14 codons p.(Asp336Serfs*6). The genotypes at this variant showed perfect cosegregation with the ichthyosis phenotype in a large family comprising 14 cases and 72 controls. ABHD5 encodes an acyltransferase required for lipid metabolism. In humans, variants in ABHD5 cause Chanarin-Dorfman syndrome, a neutral lipid storage disease with ichthyosis. Our data in dogs together with the knowledge on the effects of ABHD5 variants in humans strongly suggest ABHD5:c.1006_1019del as candidate causative genetic variant for a new canine form of ichthyosis, which we propose to designate as Golden Retriever ichthyosis type 2 (ICH2).

Variants Affecting the C-Terminal Tail of UNC93B1 Are Not a Common Risk Factor for Systemic Lupus Erythematosus

Systemic lupus erythematosus (SLE) is a heterogeneous multifactorial disease. Upregulated TLR7 signaling is a known risk factor for SLE. Recently, it was shown that specific genetic variants in UNC93B1 affect the physiological regulation of TLR7 signaling and cause characteristic autoimmune phenotypes with monogenic autosomal recessive inheritance in mutant mice and dogs. We therefore hypothesized that homologous variants in the human UNC93B1 gene might be responsible for a fraction of human SLE patients. We analyzed 536 patients of the Swiss SLE Cohort Study for the presence of genetic variants affecting the C-terminal tail of UNC93B1. None of the investigated patients carried bi-allelic UNC93B1 variants that were likely to explain their SLE phenotypes. We conclude that genetic variants affecting the C-terminal tail of UNC93B1 are not a common risk factor for SLE. It cannot be excluded that such variants might contribute to other heritable autoimmune diseases.

A COL7A1 Variant in a Litter of Neonatal Basset Hounds with Dystrophic Epidermolysis Bullosa

We investigated three neonatal Basset Hound littermates with lesions consistent with epidermolysis bullosa (EB), a group of genetic blistering diseases. A clinically normal bitch was bred to her grandfather by artificial insemination. Out of a litter of seven puppies, two affected puppies died and one was euthanized, with these puppies being submitted for diagnostic necropsy. All had multiple bullae and ulcers involving the nasal planum and paw pads, as well as sloughing claws; one puppy also had oral and esophageal ulcers. The complete genome of one affected puppy was sequenced, and 37 known EB candidate genes were assessed. We found a candidate causative variant in COL7A1, which encodes the collagen VII alpha 1 chain. The variant is a complex rearrangement involving duplication of a 107 bp region harboring a frameshift deletion of 7 bp. The variant is predicted to truncate more than 75% of the open reading frame, p.(Val677Serfs*11). Targeted genotyping of this duplication confirmed that all three affected puppies were homozygous for the duplication, whereas 12 unaffected Basset Hounds did not carry the duplication. This variant was also not seen in the genomes of more than 600 dogs of other breeds. COL7A1 variants have been identified in humans and dogs with dystrophic epidermolysis bullosa (DEB). The identified COL7A1 variant therefore most likely represents the causative variant and allows the refinement of the preliminary EB diagnosis to DEB.

LAMB3 Missense Variant in Australian Shepherd Dogs with Junctional Epidermolysis Bullosa

In a highly inbred Australian Shepherd litter, three of the five puppies developed widespread ulcers of the skin, footpads, and oral mucosa within the first weeks of life. Histopathological examinations demonstrated clefting of the epidermis from the underlying dermis within or just below the basement membrane, which led to a tentative diagnosis of junctional epidermolysis bullosa (JEB) with autosomal recessive inheritance. Endoscopy in one affected dog also demonstrated separation between the epithelium and underlying tissue in the gastrointestinal tract. As a result of the severity of the clinical signs, all three dogs had to be euthanized. We sequenced the genome of one affected puppy and compared the data to 73 control genomes. A search for private variants in 37 known candidate genes for skin fragility phenotypes revealed a single protein-changing variant, LAMB3:c.1174T>C, or p.Cys392Arg. The variant was predicted to change a conserved cysteine in the laminin β3 subunit of the heterotrimeric laminin-322, which mediates the binding of the epidermal basement membrane to the underlying dermis. Loss-of-function variants in the human LAMB3 gene lead to recessive forms of JEB. We confirmed the expected co-segregation of the genotypes in the Australian Shepherd family. The mutant allele was homozygous in two genotyped cases and heterozygous in three non-affected close relatives. It was not found in 242 other controls from the Australian Shepherd breed, nor in more than 600 other controls. These data suggest that LAMB3:c.1174T>C represents the causative variant. To the best of our knowledge, this study represents the first report of a LAMB3-related JEB in domestic animals.

Novel Brown Coat Color (Cocoa) in French Bulldogs Results from a Nonsense Variant in HPS3

Brown or chocolate coat color in many mammalian species is frequently due to variants at the B locus or TYRP1 gene. In dogs, five different TYRP1 loss-of-function alleles have been described, which explain the vast majority of dogs with brown coat color. Recently, breeders and genetic testing laboratories identified brown French Bulldogs that did not carry any of the known mutant TYRP1 alleles. We sequenced the genome of a TYRP1^+/+ brown French Bulldog and compared the data to 655 other canine genomes. A search for private variants revealed a nonsense variant in HPS3, c.2420G>A or p.(Trp807*). The brown dog was homozygous for the mutant allele at this variant. The HPS3 gene encodes a protein required for the correct biogenesis of lysosome-related organelles, including melanosomes. Variants in the human HPS3 gene cause Hermansky–Pudlak syndrome 3, which involves a mild form of oculocutaneous albinism and prolonged bleeding time. A variant in the murine Hps3 gene causes brown coat color in the cocoa mouse mutant. We genotyped a cohort of 373 French Bulldogs and found a strong association of the homozygous mutant HPS3 genotype with the brown coat color. The genotype–phenotype association and the comprehensive knowledge on HPS3 function from other species strongly suggests that HPS3:c.2420G>A is the causative variant for the observed brown coat color in French Bulldogs. In order to clearly distinguish HPS3-related from the TYRP1-related brown coat color, and in line with the murine nomenclature, we propose to designate this dog phenotype as “cocoa”, and the mutant allele as HPS3^co.

A Missense Variant Affecting the C-Terminal Tail of UNC93B1 in Dogs with Exfoliative Cutaneous Lupus Erythematosus (ECLE)

Cutaneous lupus erythematosus (CLE) in humans encompasses multiple subtypes that exhibit a wide array of skin lesions and, in some cases, are associated with the development of systemic lupus erythematosus (SLE). We investigated dogs with exfoliative cutaneous lupus erythematosus (ECLE), a dog-specific form of chronic CLE that is inherited as a monogenic autosomal recessive trait. A genome-wide association study (GWAS) with 14 cases and 29 controls confirmed a previously published result that the causative variant maps to chromosome 18. Autozygosity mapping refined the ECLE locus to a 493 kb critical interval. Filtering of whole genome sequence data from two cases against 654 controls revealed a single private protein-changing variant in this critical interval, UNC93B1:c.1438C>A or p.Pro480Thr. The homozygous mutant genotype was exclusively observed in 23 ECLE affected German Shorthaired Pointers and an ECLE affected Vizsla, but absent from 845 controls. UNC93B1 is a transmembrane protein located in the endoplasmic reticulum and endolysosomes, which is required for correct trafficking of several Toll-like receptors (TLRs). The p.Pro480Thr variant is predicted to affect the C-terminal tail of the UNC93B1 that has recently been shown to restrict TLR7 mediated autoimmunity via an interaction with syndecan binding protein (SDCBP). The functional knowledge on UNC93B1 strongly suggests that p.Pro480Thr is causing ECLE in dogs. These dogs therefore represent an interesting spontaneous model for human lupus erythematosus. Our results warrant further investigations of whether genetic variants affecting the C-terminus of UNC93B1 might be involved in specific subsets of CLE or SLE cases in humans and other species.

[Reproducibility of fasting serum cholesterol and triglycerides in ambulatory patients with mixed hyperlipidemia]

Intraindividual variability of serum lipid concentrations in normal volunteers and in patients with hyperlipidemia is substantial. The aim of this study was to investigate prospectively the reproducibility of fasting serum triglyceride and total cholesterol concentrations in primary health care patients with combined hyperlipidemia, i.e. under conditions of daily medical practice. Secondary forms of hyperlipidemia were excluded. 19 general medical outpatients with primary combined hyperlipidemia were studied. Serum total cholesterol and triglyceride concentrations were measured after an overnight fast at 08.00 h 4 times at weekly intervals. To study the influence of alcohol intake on serum lipid concentrations, total cholesterol and triglycerides were measured without alcohol influence and 12 hours after consumption of a mean of 100 g alcohol in the evening. In 19 patients (10 males, 9 females, mean age 55 years, body mass index 27.9 +/- 4.4 kg/m2), mean +/- SD of serum triglycerides was 3.97 +/- 1.8 mmol/l and of total cholesterol 7.9 +/- 1.8 mmol/l. The combined intraindividual and interassay coefficient of variation was 18.7 +/- 8.2% for triglycerides and 5.1 +/- 2.5% for total cholesterol. Fasting serum triglycerides (3.5 +/- 1.1 vs. 3.7 +/- 1.4 mmol/l) and total cholesterol (7.6 +/- 1.4 vs. 7.8 +/- 1.0 mmol/l) did not significantly change 12 hours after acute alcohol consumption. Patients with primary combined hyperlipidemia in a primary health care setting show small intraindividual variations of overnight fasted serum triglyceride and total cholesterol concentrations. Moderate alcohol consumption 12 hours before blood sampling does not significantly affect triglyceride and cholesterol values.