Canine Epithelial Skin Tumours: Expression of the Stem Cell Markers Lgr5, Lgr6 and Sox9 in Light of New Cancer Stem Cell Theories.

Evidence is accumulating that tumour development is driven by cancer stem cells (CSCs). In order to understand the presence and potential contribution of stem cells (SCs) as tumour-initiating cells in canine cutaneous tumours, we selected three putative SC markers (Lgr5, Lgr6 and Sox9) and investigated their expression pattern, level of protein and mRNA expression, in 43 canine hair follicle (HF) and 18 canine cutaneous epidermal tumours by immunohistochemistry and qRT-PCR, using normal skin samples as controls. Lgr5 protein expression was not detected in epidermal and HF tumours; however, Lgr5 mRNA overexpression was evident in some HF tumours. Sox9 was expressed in several tumour cases, both at the protein and mRNA level. The Lgr6 antibody tested was not suitable for formalin-fixed paraffin-embedded tissue samples, but Lgr6 gene showed higher expression in several samples of both HF and epidermal tumours compared with normal skin. Significantly higher mRNA expression levels of the three SC markers were found in trichoblastomas (TB) compared with basal cell carcinomas (BCC). The present results indicated that canine HF and epidermal tumours might have common tumour-initiating cells. The mRNA expression of the three selected SC markers, especially Lgr5, could be potentially useful in the distinction between canine TB and BCC

A KRT71 Loss-of-Function Variant Results in Inner Root Sheath Dysplasia and Recessive Congenital Hypotrichosis of Hereford Cattle

Genodermatoses, such as heritable skin disorders, mostly represent Mendelian conditions. Congenital hypotrichosis (HY) characterize a condition of being born with less hair than normal. The purpose of this study was to characterize the clinicopathological phenotype of a breed-specific non-syndromic form of HY in Hereford cattle and to identify the causative genetic variant for this recessive disorder. Affected calves showed a very short, fine, wooly, kinky and curly coat over all parts of the body, with a major expression in the ears, the inner part of the limbs, and in the thoracic-abdominal region. Histopathology showed a severely altered morphology of the inner root sheath (IRS) of the hair follicle with abnormal Huxley and Henle’s layers and severely dysplastic hair shafts. A genome-wide association study revealed an association signal on chromosome 5. Homozygosity mapping in a subset of cases refined the HY locus to a 690 kb critical interval encompassing a cluster of type II keratin encoding genes. Protein-coding exons of six positional candidate genes with known hair or hair follicle function were re-sequenced. This revealed a protein-changing variant in the KRT71 gene that encodes a type II keratin specifically expressed in the IRS of the hair follicle (c.281delTGTGCCCA; p.Met94AsnfsX14). Besides obvious phenocopies, a perfect concordance between the presence of this most likely pathogenic loss-of-function variant located in the head domain of KRT71 and the HY phenotype was found. This recessive KRT71-related form of hypotrichosis provides a novel large animal model for similar human conditions. The results have been incorporated in the Online Mendelian Inheritance in Animals (OMIA) database (OMIA 002114-9913)

Independent COL5A1 Variant 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

Genomic and Transcriptomic Characterization of Atypical Recurrent Flank Alopecia in the Cesky Fousek.

Non-inflammatory alopecia is a frequent skin problem in dogs, causing damaged coat integrity and compromised appearance of affected individuals. In this study, we examined the Cesky Fousek breed, which displays atypical recurrent flank alopecia (aRFA) at a high frequency. This type of alopecia can be quite severe and is characterized by seasonal episodes of well demarcated alopecic areas without hyperpigmentation. The genetic component responsible for aRFA remains unknown. Thus, here we aimed to identify variants involved in aRFA using a combination of histological, genomic, and transcriptomic data. We showed that aRFA is histologically similar to recurrent flank alopecia, characterized by a lack of anagen hair follicles and the presence of severely shortened telogen or kenogen hair follicles. We performed a genome-wide association study (GWAS) using 216 dogs phenotyped for aRFA and identified associations on chromosomes 19, 8, 30, 36, and 21, highlighting 144 candidate genes, which suggests a polygenic basis for aRFA. By comparing the skin cell transcription pattern of six aRFA and five control dogs, we identified 236 strongly differentially expressed genes (DEGs). We showed that the GWAS genes associated with aRFA are often predicted to interact with DEGs, suggesting their joint contribution to the development of the disease. Together, these genes affect four major metabolic pathways connected to aRFA: collagen formation, muscle structure/contraction, lipid metabolism, and the immune system

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

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

Identification of Two Independent COL5A1 Variants in Dogs with Ehlers-Danlos Syndrome.

The Ehlers-Danlos syndromes (EDS) are a heterogeneous group of heritable disorders affecting connective tissues. The mutations causing the various forms of EDS in humans are well characterized, but the genetic mutations causing EDS-like clinical pathology in dogs are not known, thus hampering accurate clinical diagnosis. Clinical analysis of two independent cases of skin hyperextensibility and fragility, one with pronounced joint hypermobility was suggestive of EDS. Whole-genome sequencing revealed de novo mutations of COL5A1 in both cases, confirming the diagnosis of the classical form of EDS. The heterozygous COL5A1 p.Gly1013ValfsTer260 mutation characterized in case 1 introduced a premature termination codon and would be expected to result in α1(V) mRNA nonsense-mediated mRNA decay and collagen V haploinsufficiency. While mRNA was not available from this dog, ultrastructural analysis of the dermis demonstrated variability in collagen fibril diameter and the presence of collagen aggregates, termed 'collagen cauliflowers', consistent with COL5A1 mutations underlying classical EDS. In the second case, DNA sequencing demonstrated a p.Gly1571Arg missense variant in the COL5A1 gene. While samples were not available for further analysis, such a glycine substitution would be expected to destabilize the strict molecular structure of the collagen V triple helix and thus affect protein stability and/or integration of the mutant collagen into the collagen V/collagen I heterotypic dermal fibrils. This is the first report of genetic variants in the COL5A1 gene causing the clinical presentation of EDS in dogs. These data provided further evidence of the important role of collagen V in dermal collagen fibrillogenesis. Importantly, from the clinical perspective, we showed the utility of DNA sequencing, combined with the established clinical criteria, in the accurate diagnosis of EDS in dogs

A Mutation in the SUV39H2 Gene in Labrador Retrievers with Hereditary Nasal Parakeratosis (HNPK) Provides Insights into the Epigenetics of Keratinocyte Differentiation.

Peer reviewe

An Intronic MBTPS2 Variant Results in a Splicing Defect in Horses with Brindle Coat Texture

We investigated a family of horses exhibiting irregular vertical stripes in their hair coat texture along the neck, back, hindquarters, and upper legs. This phenotype is termed "brindle" by horse breeders. We propose the term "brindle 1 (BR1)" for this specific form of brindle. In some BR1 horses the stripes were also differentially pigmented. Pedigree analyses were suggestive of a monogenic X-chromosomal semi-dominant mode of inheritance. Haplotype analyses identified a 5 Mb candidate region on chromosome X. Whole genome sequencing of 4 BR1 and 60 non-brindle horses identified 61 private variants in the critical interval, none of them located in an exon of an annotated gene. However, one of the private variants was close to an exon/intron boundary in intron 10 of the MBPTS2 gene encoding the membrane bound transcription factor peptidase, site 2 (c.1437+4T>C). Different coding variants in this gene lead to three related genodermatoses in human patients. We therefore analyzed MBPTS2 transcripts in skin and identified an aberrant transcript in a BR1 horse, which lacked the entire exon 10 and parts of exon 11. The MBPTS2:c1437+4T>C variant showed perfect co-segregation with the brindle phenotype in the investigated family and was absent from 457 control horses of diverse breeds. Altogether, our genetic data and the previous knowledge on MBTPS2 function in the skin suggest that the identified MBTPS2 intronic variant leads to partial exon skipping and causes the BR1 phenotype in horses

MKLN1 splicing defect in dogs with lethal acrodermatitis

Lethal acrodermatitis (LAD) is a genodermatosis with monogenic autosomal recessive inheritance in Bull Terriers and Miniature Bull Terriers. The LAD phenotype is characterized by poor growth, immune deficiency, and skin lesions, especially at the paws. Utilizing a combination of genome wide association study and haplotype analysis, we mapped the LAD locus to a critical interval of similar to 1.11 Mb on chromosome 14. Whole genome sequencing of an LAD affected dog revealed a splice region variant in the MKLN1 gene that was not present in 191 control genomes (chr14:5,731,405T>G or MKLN/:c.400+3A>C). This variant showed perfect association in a larger combined Bull Terrier/Miniature Bull Terrier cohort of 46 cases and 294 controls. The variant was absent from 462 genetically diverse control dogs of 62 other dog breeds. RT-PCR analysis of skin RNA from an affected and a control dog demonstrated skipping of exon 4 in the MKLN1 transcripts of the LAD affected dog, which leads to a shift in the MKLN1 reading frame. MKLN1 encodes the widely expressed intracellular protein muskelin 1, for which diverse functions in cell adhesion, morphology, spreading, and intracellular transport processes are discussed. While the pathogenesis of LAD remains unclear, our data facilitate genetic testing of Bull Terriers and Miniature Bull Terriers to prevent the unintentional production of LAD affected dogs. This study may provide a starting point to further clarify the elusive physiological role of muskelin 1 in vivo.Peer reviewe