Next Article in Journal
Structural Insights into TOR Signaling
Previous Article in Journal
Exploring the lncRNAs Related to Skeletal Muscle Fiber Types and Meat Quality Traits in Pigs
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 11
Issue 8
10.3390/genes11080884
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Transcriptome Profiling and Differential Gene Expression in Canine Microdissected Anagen and Telogen Hair Follicles and Interfollicular Epidermis

by
Dominique J. Wiener
1,*,
Kátia R. Groch
1,
Magdalena A.T. Brunner
2,3,
Tosso Leeb
2,3,
Vidhya Jagannathan
2 and
Monika M. Welle
3,4
1
Department of Veterinary Pathobiology, College of Veterinary Medicine & Biomedical Science, Texas A&M University, College Station, TX 77843, USA
2
Institute of Genetics, Vetsuisse Faculty, University of Bern, 3012 Bern, Switzerland
3
Dermfocus, Vetsuisse Faculty, University Hospital of Bern, 3010 Bern, Switzerland
4
Institute of Animal Pathology, Vetsuisse Faculty, University of Bern, 3012 Bern, Switzerland
*
Author to whom correspondence should be addressed.
Genes 2020, 11(8), 884; https://doi.org/10.3390/genes11080884
Submission received: 30 June 2020 / Revised: 31 July 2020 / Accepted: 3 August 2020 / Published: 4 August 2020
(This article belongs to the Section Animal Genetics and Genomics)
Browse Figures
Versions Notes

Abstract

:
The transcriptome profile and differential gene expression in telogen and late anagen microdissected hair follicles and the interfollicular epidermis of healthy dogs was investigated by using RNAseq. The genes with the highest expression levels in each group were identified and genes known from studies in other species to be associated with structure and function of hair follicles and epidermis were evaluated. Transcriptome profiling revealed that late anagen follicles expressed mainly keratins and telogen follicles expressed GSN and KRT15. The interfollicular epidermis expressed predominately genes encoding for proteins associated with differentiation. All sample groups express genes encoding for proteins involved in cellular growth and signal transduction. The expression pattern of skin-associated genes in dogs is similar to humans. Differences in expression compared to mice and humans include BMP2 expression mainly in telogen and high KRT17 expression in the interfollicular epidermis of dogs. Our data provide the basis for the investigation of the structure and function of canine skin or skin disease and support the use of dogs as a model for human cutaneous disease by assigning gene expression to specific tissue states.

1. Introduction

The skin is the largest organ, with various essential functions, including protection/barrier against outside influences and thermoregulation. The epidermis and its appendages form the outer layer of the skin and are the first shield of the body from external harm [1]. To fulfill the protective functions, an intact epidermis and appropriate numbers of hair follicles (HF) are necessary.
In all species, the interfollicular epidermis (IFE) is composed of a squamous, stratified, multilayered epithelium that frequently self-renews, involving epidermal stem and progenitor cells that reside in the basal layer of the epidermis [2,3]. After the basal epidermal cells detach from the basement membrane, their proliferation ceases and they undergo a terminal cell differentiation program, during which keratinocytes move upwards until they reach the cornified layer as flattened denucleated keratinocytes and are finally shed as squames at the skin surface [4,5,6].
The disturbance of this process can have a devastating impact, such as excessive water loss, invasion of pathogens and development of cancer, which may become life-threatening. The main role of HFs is the production of hair shafts to maintain the hair essential for a variety of functions, such as sensation and the dispersion of sebum and pheromones. In animals, the fur also has life-preserving functions, such as thermal regulation, visual communication and protection from solar radiation and physical trauma.
In contrast to the continuous renewal of the IFE, HFs renew periodically, going through life-long recurrent phases of growth (anagen), regression (catagen) and rest (telogen) [7,8]. The hair shaft is produced during anagen. The renewal of the IFE and the HF relies on stem cells and their progenitors. Several stem cell markers, such as LGR5, CD34, SOX9 and KRT15 (HF bulge stem cells markers) [9,10,11,12,13] or LGR-6 (expressed in the upper HF, sebaceous gland and IFE) have been described in humans and rodents [14,15] and their expression has also been investigated in dogs [10,16,17,18]. Stem cell quiescence and activation, as well as the subsequent events leading to epidermal differentiation and the maintenance of the hair cycle, are tightly regulated and involve numerous signals derived from epithelial, neuroendocrine and mesenchymal cells [19,20,21]. The reciprocal communication between signals from the immediate stem cell microenvironment within the IFE or the HF and signals from the larger macroenvironment are assembled in the stem cells. The interplay of proteins and all transcription factors implicated in several signaling pathways finally result in stem cell activation or quiescence. Once activated, stem cells give rise to rapidly proliferating transient amplifying cells, which terminally differentiate into the epidermal and the different hair follicle layers and the layers of the hair shaft. Important signaling pathways that regulate stem cell quiescence, proliferation and differentiation include the wingless-type mouse mammary tumor virus integration site (WNT)/β-catenin, transforming growth factor (TGF)-β, fibroblast growth factor (FGF)5, bone morphogenic protein (BMP), FGF18 and Notch signaling (reviewed in [22,23]). In general, BMP signals derived from dermal, adipose and epithelial tissue repress cell proliferation, while WNT and other signals promote stem cell activation and growth.
The maintenance of the protective function of the skin is not only dependent on a functional stem cell compartment but also on proteins involved in the structural and functional integrity. Examples are keratin intermediate filaments and proteins involved in cell–cell adhesion such as adherens junctions and desmosomes [24,25].
The immune function of the epidermis is maintained by the presence of specialized antigen presenting cells (Langerhans cells) and of numerous pattern recognition receptors, such as Toll-like receptors [26,27]. HFs harbor, for example, Langerhans cells and CD4+ or CD8+ T cells and produce chemokines and cytokines important in immune defense [28,29].
To date, many molecules involved in skin homeostasis have been extensively studied and their role in disease development in humans and rodents is well known [30,31,32]. Additionally, in domestic animals, numerous gene variants involved in genodermatoses were identified by whole genome sequencing and other approaches [33]. However, for some genes associated with skin disease, the function and anatomic location of the expressed proteins are unknown. Examples for this are SUV39H2 variants in nasal parakeratosis or FAM83G variants in hereditary footpad hyperkeratosis [34,35].
Gene expression analysis is a useful tool to study the molecular mechanisms underlying skin development and homeostasis. The identification of differentially expressed genes (DEGs) in numerous disease phenotypes of humans and animals is increasingly applied in research.
The transcriptome sequencing of skin biopsies in dogs to identify disease phenotypes revealed numerous deregulated genes, some of which are well known and whose expression can thus be associated with a specific location. However, many of the expressed genes have neither been associated with disease nor attributed to a specific tissue or cell type [17,36]. Studies based on the transcriptome profiling of murine epidermis and HFs [37] and human skin [38,39,40] are available. To the best of our knowledge, for domestic animals, transcriptome data for specific anatomic locations of healthy skin are only available from plugged hairs or cell cultures derived from microdissected HFs of cashmere goats [41,42,43,44]. Although dogs are increasingly recognized as an attractive model for human disease, no whole transcriptome expression data of microdissected skin tissue have been published [45,46].
The aim of our study was to identify global gene expression and differential gene expression in microdissected canine HFs in late anagen and telogen and in the IFE using transcriptome analysis by sequencing of mRNA (RNA-seq). Furthermore, we performed ontological molecular analyses for DEGs in telogen and late anagen HFs. Our data provide the basis for the investigation of the structure and function of the canine skin or canine skin disease and for the use of dogs as a model for human cutaneous disease, by assigning the expression of different genes to specific tissue types.

2. Materials and Methods

2.1. Ethics Statement

All skin samples were taken from dogs that were euthanized for reasons unrelated to this study (such as internal disease or accident). All dogs were privately owned pets that were euthanized according to ethical standards of Texas A&M University, College Station, USA. The owners signed a consent form to donate their dead animal.

2.2. Skin Biopsies

Skin samples measuring approx. 2 × 2 cm were taken from the trunk from eight freshly euthanized dogs (dead for no longer than 12 h) of various breeds and age with no skin abnormalities (Supplementary Table S1). Most of the samples were derived from female spayed dogs.

2.3. Microdissection

From the skin samples of each dog, 20 single HFs were microdissected as described before by cutting above the sebaceous glands to separate HFs from IFE and, subsequently, the sebaceous glands were removed (Supplementary Figure S1) [47,48]. The follicular stages in each sample were assigned by the length of the hair shaft and the depth of the follicle within the tissue (according to the canine hair cycle guide by Müntener et al. [49]). All anagen HFs collected were in the anagen VI stage, identified by the hair bulbs located deep in the subcutis and the hair shafts emerging through the HF ostia [49]. While, in most dogs, either anagen or telogen follicles could be identified, in both dogs 4 and 8, late anagen and telogen hair cycle stages were visible. In these dogs, 20 single HFs in telogen and 20 HFs in late anagen were collected.
The IFE was microdissected by cutting it off from the dermis and the HFs were removed with a scalpel blade. Between 20 and 30 small IFE pieces of 1–3 mm were collected. Microdissected HFs and IFE pieces were stored in RNAlater at −20 °C until RNA extraction.

2.4. Histological Analysis

Whole skin, adjacent to skin used for microdissection, was collected and fixed immediately in 10% buffered formalin. The formalin-fixed tissue was embedded in paraffin, cut as a 4 µm section and stained with hematoxylin and eosin (H&E) using standard procedures. The histological specimens were evaluated to confirm that no pathologic changes were present and to confirm the hair cycle stages determined by microdissection.

2.5. RNA Extraction and Transcriptome Sequencing (RNA-seq)

RNA for transcriptome sequencing was extracted from the 20 pooled microdissected HF of the same cycle stage and the IFE collected from each dog. RNA extraction, library preparation and RNA-seq were performed as described in a previous study [17,36]. In brief, total RNA was extracted using the RNeasy Micro Kit (Qiagen, Hombrechtikon, Switzerland) including proteinase K digestion according to the manufacturer’s protocol. Total RNA content was measured using the Qubit® 2.0 Fluorometer (Invitrogen, Thermo Fisher, Basel, Switzerland) and RNA quality was assessed with a Bioanalyzer (Agilent 2100; Agilent Technologies, Basel, Switzerland). High quality RNA (RIN > 9) was reverse transcribed into cDNA with the SMART-Seq v4 Ultra Low Input RNA Kit (Takara, Saint-Germain-en-Laye, France) and libraries were constructed following the TruSeq® Nano DNA Library Prep (Illumina, San Diego, CA, USA) protocol for RNA sequencing. Multiplexed total cDNA libraries were sequenced on two lanes using an Illumina HiSq3000 instrument with 100 bp single-end sequencing cycles. On average, thirty million reads were collected, converted into FASTQ file format and demultiplexed.
The data are available from the European Nucleotide Archive (ENA), study accession number PRJEB21761 and sample accessions SAMEA7016531-SAMEA7016551 [50].

2.6. Mapping to Reference Genome and Gene Expression Analysis

The reads were quality checked with fastqc [51]. Low quality bases (−q 15) and illumina adaptors were trimmed with fastp (PMID: 30423086). All reads that passed quality control were mapped with an average of 90% uniquely mapped reads to the dog reference genome (Can.Fam3.1) by STAR aligner version 2.6.0 [52] as described in a previous study [17,36]. The read abundance was calculated using HTseq and a gtf file (version 105, downloaded and modified in March 2019) obtained from NCBI Can.Fam3.1 annotation release 105 [53]. The DESeq2 R package [52] in RSTUDIO [54] was used to read the HTseq count data. Pre-filtering was applied for low-expressed/non-expressed genes, excluding those genes where the sum of raw counts across samples was below 10.
The raw counts were subjected to a regularized-logarithm transformation and a principal component analysis (PCA) was performed using R software package “ggplot2” [55] to visualize the clustering of the samples derived from the epidermis and different hair cycle stages. In addition, PCAs, including breed and age, were performed.
To take the gene length into account, raw counts were normalized using the Transcripts Per Million (TPM) method to examine the most expressed genes in each group.
Gene descriptions and functions were identified using the National Center for Biotechnology Information (NCBI) Gene database [56].

2.7. Differentially Gene Expression

DEGs between late anagen and telogen samples were analyzed by using a generalized linear model (GLM) of DESeq2 R Package [52] to the pre-filtered raw count data, assuming a negative binomial distribution. Adjusted p-values were calculated according to the Benjamini Hochberg method. Transcripts were considered to be differentially expressed with a false discovery rate (FDR) of <0.01 and log2 fold change of 1.5 as the threshold. A heatmap was drawn for the top 30 DEGs sorted based on log2 fold change. The heat map was obtained using DESeq2 normalized gene counts and “heatmap.2” function of “gplots” R package.

2.8. Pathway and Gene Enrichment Analysis

Gene Ontology (GO) analysis was applied to analyze the functional enrichment and significant pathways associated with DEGs using the Database for Annotation Visualization and Integrated Discovery (DAVID) [57,58]. The gene symbols of the DEGs were uploaded into the database and the dog genome (Canis lupus familiaris) was selected as the background parameter. The upregulated and downregulated DE gene sets comparing late anagen and telogen were used to extract functional gene sets by “functional annotation chart” analyses. Annotations whose Benjaminin corrected p-value were lower than 0.05 were considered statistically significant.

3. Results

3.1. Microdissection

In skin samples from six dogs, either anagen or telogen HF stages predominated. However, in samples from two dogs, both hair cycle stages were present in roughly equal proportions (dogs 4 and 8).

3.2. Histologic Findings

Histological examination revealed that, while in most skin samples, both anagen and telogen hair cycle stages were present, one of the stages predominated and corresponded to the stage recognized during microdissection. No pathologic changes were present in the samples included in this study and the epidermis was of normal thickness in all samples.

3.3. Gene Expression Profile in Anagen HF, Telogen HF and IFE

Transcriptome sequencing revealed the expression of 21,790 genes (genes are listed in Supplementary Table S2). Principal component analysis (PCA), based on gene expression levels, resulted in the distinct clustering of the different groups (IFE, telogen and late anagen, respectively) (Supplementary Figure S2). However, two of the telogen samples (dogs 2 and 5) clustered distant from the other telogen samples but also from late anagen samples. We therefore concluded that these two samples represent early anagen HFs and not telogen HFs. In fact, these two hair cycle stages cannot be differentiated histologically [49]. These two samples were excluded from the following evaluations in order to compare late anagen and telogen stages only. The influence of age and breed on the differential gene expression was excluded since clustering in PCA plots was not associated with these conditions (Supplementary Figures S3 and S4, respectively). As all samples, except one, were derived from female spayed dogs, the influence of gender could not be evaluated.

3.4. Highly Expressed Genes

The genes with the highest expression in anagen and telogen HFs and in the IFE were identified and are presented in a Venn diagram using GeneVenn [59] (Figure 1). A high expression level was chosen (average normalized counts > 5000) to single out the most highly expressed genes in the different groups. In late anagen HFs, genes encoding for structural proteins (KRT16, KRT25, KRT31, KRT33A, KRT33B, KRT35, KRT71, KRT83, KRT85, LOC100684920) were highly expressed. The LOC gene is predicted to encode keratin-associated protein 4-12-like In telogen HFs, gelsolin (GSN) was highly expressed, which encodes a calcium-regulated protein that functions in the assembly and disassembly of actin filaments and therefore plays a role in cellular structure. KRT15, a gene encoding for a stem cell marker described to be present in anagen and telogen hair cycle stages in dogs [10], is one of the highest expressed gene in telogen HFs, but also shows a strong expression in late anagen HFs and the IFE, but to a lesser degree (Table 1, Supplementary Table S2).
Most genes that are highly expressed in the IFE are known to be associated with differentiation of human keratinocytes or encode for proteins in the human epidermis (KRT1, KRT10, LOR and ITM2B). One of the most expressed genes in the epidermis, LOC106559288, is predicted to encode a long non-coding RNA (lncRNA).
The basal cell markers, KRT14 and KRT5, were strongly expressed in all three groups (IFE, late anagen and telogen), as well as MIR8858, ACTG1 and ANXA2, which encode for structural proteins or for proteins associated with cellular growth and signal transduction.

3.5. Gene Expression Profile of Selected Genes

A thorough literature search was performed to search for genes most commonly associated with hair cycle, stem cell markers, HF and IFE structure/function and immune function. Altogether, 131 genes were identified during this search. Of those, 92 genes had an average normalized count of >10 in at least one of the groups (IFE, late anagen, telogen) and were further analyzed. They were assigned to a specific biological function (hair cycle (n = 32), stem cell markers (n = 12), HF (n = 25) and IFE (n = 12) structure/function and immune function (n = 11), respectively) and used for further analyses. Genes that were expressed at a very low level in all three groups (average normalized count < 10) were not included in further investigations. The average expression of each gene was calculated for IFE, anagen and telogen samples (Table 1).
Genes expressed at least four-fold higher (average normalized count > 100) in one group compared to the others are mentioned in the following. Genes expressed higher in anagen than in the other groups were mainly associated with HF structure (HOXC13, KRT25, KRT31, KRT33A, KRT33B, KRT34, KRT35, KRT71, KRT83, KRT85, KRT86, KRTAP16-1, PMEL and TCHH). PMEL encodes for melanogenic proteins. As melanogenesis takes place in the inferior portion of HFs, which is only present in anagen HFs, it is not surprising that PMEL is highly expressed in late anagen HFs. The other genes expressed higher in anagen were associated with the hair cycle (FOXN1, HSD17B14 and MSX2).
In contrast, genes expressed higher in telogen were associated mostly with the hair cycle (BMP2, FGF18, MMP7), but also with stem cell marker (KRT15) or follicular structure (KRT17).
The selected genes that showed higher expression in the IFE are encoding for epidermal structure (DSC1, KRT1, KRT10 and LOR). Some of the hair cycle-associated genes (DKK3 and WNT3) are highly expressed in the IFE. Selected genes that expressed strongly in all three groups (average normalized count > 1000 counts) were associated with epidermal and HF structure and function, such as DSP, KRT5, KRT14, KRT16 and KRT17, as well as the stem cell marker KRT15.

3.6. Expression of Genes Associated with Canine Genodermatoses

The presence and the expression level of genes identified to be associated with canine genodermatoses either affecting the IFE or HFs were analyzed in our data set of normal tissue. All genes known to be associated with canine genodermatoses could be retrieved in our samples. In most cases genes associated with either primarily epidermal or alopecic disease are expressed in both HFs and IFE and it remains to be elucidated why pathology affects only one anatomical structure. The genes that were identified with specific genodermatoses are depicted in Table 2.

3.7. Differentially Expressed Genes in Anagen and Telogen Hair Cycle Stages

A total of 1236 DEGs were identified in late anagen HFs compared with telogen HF samples, including 858 upregulated genes and 378 downregulated genes (Supplementary Table S3). The 30 most significantly DEGs were further evaluated (Figure 2). All of these genes were expressed significantly more in anagen. Most genes designated with LOC encode keratin-associated proteins, therefore, 90% of the most DEGs encode keratins and keratin-associated proteins. One of the most expressed DEGs (LOC111092804), is uncharacterized and its function is yet unknown.
From the 92 preselected genes known to be associated with HF structure and function, hair cycle, stem cells and immune function, 46 genes were differentially expressed comparing late anagen versus telogen hair cycle stages (Figure 3). As presented above, several genes associated with HF structure are expressed mainly in anagen. Genes associated with the hair cycle and expressed significantly more in late anagen are FOXN1, HSD17B14 and MSX2 and in telogen FGF18 and MMP7. Immune response-associated genes are more strongly expressed in telogen HFs, such as HLA-DRB1 and CD74 (both encoding for MHC class II proteins). KRT75, IVL and LOR were strongly expressed by both hair cycle stages, but KRT75 and IVL were expressed significantly more in late anagen and LOR in telogen.

3.8. Functional Classification of Differentially Expressed Genes

Differentially expressed genes of three GO categories (cellular component, biological process and molecular function) were either up- or down-regulated in anagen. Upregulated and downregulated DEGs enriched GO terms with a Benjamini adjusted p value < 0.05 are shown in Supplementary Figure S5. The upregulated genes in anagen are genes related to “keratin filament”, “intermediate filament” and “structural molecule activity”. In contrast, downregulated genes are involved in the “MHC class II protein complex”, “Peptide antigen binding” and “Transmembrane receptor protein tyrosine kinase signaling pathway”.

4. Discussion

We studied the gene expression profile of canine microdissected HFs in telogen and late anagen and in the IFE to gain information about the canine-specific expression pattern. It is well known that gene expression in different species varies and our data provide the basis to interpret studies involving gene expression in canine skin disorders, either inherited or acquired, in more detail and with more confidence.
A comparative analysis of the 100 most highly expressed skin-associated genes in human and mouse transcriptomes revealed a 30.2% identity and thus highlights the significant differences between these two species [85]. To investigate similarities/differences in the expression of canine skin-associated genes in relation to humans and mice, we compared the data obtained from our canine samples to the 100 most highly expressed genes in the skin of mice and humans. We found that about 60% of the genes highly expressed in human skin are also strongly expressed in dogs (>100 read counts), whereas only about 45% of the top 100 genes expressed in murine skin were strongly expressed in dogs. Eighteen out of the top 100 skin-associated genes in humans were not expressed in our tissue samples, whereas, compared to the top 100 skin-associated genes in mice, 34 genes were not expressed in our samples (Supplementary Table S4). Of note, the expression data sets derived from human and murine skin were obtained from whole skin and therefore, not surprisingly, contained some highly expressed genes encoding for the mesenchymal components of the skin. Although comparing the expression pattern of whole skin and microdissected tissue is suboptimal, we can conclude that more genes expressed in mice are not expressed in the dog samples as compared to the human skin and gene expression in canine skin is more similar to human skin than in mouse skin. Importantly it has to be mentioned that many of the genes expressed in both, murine and human skin were also strongly expressed (average normalized counts > 100 counts) in at least one of our sample groups (e.g., DSP, KRT1, KRT5, KRT10, KRT14, KRT15, KRT17 and LOR).
These results support that the dog is a good model for human skin disease, as we show in this study that the canine gene expression in IFE and HFs is much more similar to human follicular and epidermal gene expression, as compared to the mouse. Furthermore, our results emphasize the similarity between human and dog signaling associated with hair growth and highlight a more different gene expression pattern in mice. Similarity in the gene expression patterns of dogs and humans, as opposed to mice, can be seen, for example, in WNT expression. In the human anagen HFs, WNT3 seems to be the most prominent WNT ligand, whereas in the mouse, Wnt10a is one of the most prominent WNT ligands [86,87]. Similar to the human anagen HFs, WNT3 is the most prominently expressed WNT ligand in our anagen samples in dogs.
Examples to underline the difference in the gene expression pattern between dogs and mice is the expression of KRT17 and BMP2. KRT17 was highly expressed in all our samples (Table 1) but was expressed the most in telogen HFs. This is in contrast to KRT17 expression in mice, where it is expressed weakly in telogen and expression increases during early anagen and has stable expression in mid- and late anagen [88]. Furthermore, in our samples, BMP2 was mainly expressed in telogen HFs. This is in contrast to BMP2 expression in mice, where it is reported to be mainly expressed in late anagen and absent in late telogen [89]. In Cashmere goats, similar to our findings, BMP2 expression is higher in late telogen and early anagen phases and lower in late anagen [90]. To our knowledge, the specific expression of BMP2 in different hair cycle stages of human HFs has not yet been described and further studies are warranted to compare human BMP2 expression to that in other species.
We also found differences in the gene expression in the canine IFE and HF samples compared to both mice and human. One example for this is the lack of KRT17 expression in the epidermis of wild type mice and humans [91,92], whereas KRT17 is strongly expressed in canine IFE (Table 1). Another example is the expression of KRT75 in both anagen and telogen HFs in dogs, whereas KRT75 is constitutively expressed only in the anagen-stage HFs in mice and humans [93].
Gene expression profiling in canine microdissected HFs, comparing late anagen and telogen hair cycle stages, revealed a significantly higher expression of keratins and keratin-associated proteins, as well as of two other genes, namely PADI3 and LYG2, in anagen HFs. PADI3 modulates hair structural proteins, such as filaggrin and trichohyalin in the inner root sheath during HF formation, and is therefore equally associated with HF structure and function. LYG2 was enriched in the WNT signaling pathway in a study investigating DEGs in the skin of Liaoning cashmere goats [43]. The function of this gene in canine HFs needs still to be determined. The genes expressed strongest in telogen HFs are GSN and KRT15. However, both of these genes are also strongly expressed in the other tissue groups (Supplementary Table S2).
In the IFE, the genes that are expressed highest are associated with keratinocyte differentiation. One of the most expressed genes in canine IFE, LOC106559288, encodes a lncRNA, which is transcribed from the keratin gene cluster on dog chromosome 9. Further investigation is warranted whether this predicted lncRNA represents an artifact of the imperfect canine genome assembly and annotation or whether it has a real biological function. Genes expressed highly in all three tissue groups are involved in cellular structure, the regulation of cellular growth and signal transduction (KRT5, KRT14, ACTB, ANXA2, MIR8858).
Some of the hair cycle-associated genes (DKK3 and WNT3) are also highly expressed in the IFE, similar to mice and humans [94,95]. The role of these two genes in epidermal homeostasis is not yet clear but there is evidence that WNT/β-catenin signaling is involved in the maintenance of the epidermal stem cell pool [94]. DKK3, in contrast to other members of the Dickkopf family, does not interfere with the WNT signaling pathway and is proposed to be a tumor suppressor gene.
Analyzing genes in our data sets known to be associated with HF and IFE structure and function, stem cell markers and immune function, we found that, in anagen, the hair cycle-associated genes that are most upregulated are MSX2, FOXN1 and HSD17B14, and the most significantly differently expressed genes associated with HF structure (apart from keratins) is HOXC13. It has been suggested previously that FOXN1 and HOXC13 are required for hair keratin expression in mice and humans [96,97] and HOXC13 was also shown to be highly expressed in human anagen HFs [97]. FOXN1 (formerly called whn and hfh11) is known to be predominantly expressed in hair matrix keratinocytes, the hair cortex, and in outer root sheath keratinocytes and is known, together with Msx-2, to maintain anagen in mice [98,99]. The high expression of HOXC13, FOXN1 and MSX2 in canine HFs may suggest that these genes have a similar function, as compared to humans and mice.
In telogen, the most significantly DEGs associated with the hair cycle were FGF18 and MMP7. FGF18 is an important indicator of stem cell quiescence. In mice, Fgf18 was found to be expressed in the bulge and FGF18 peaked during telogen [100]. The deletion of Fgf18 in mice results in a shorter telogen phase [22,101]. Our results are in line with this and also support that, in dogs, FGF18 is a telogen marker. While not much is known about the expression of matrix metalloproteinase 7 (MMP7) in HFs, Mmp2 and Mmp9 expression was shown in all structures of murine HFs and in the sebaceous glands of different hair cycle stages, respectively [102]. In our canine study, MMP7 seems to be predominantly expressed in telogen HFs, suggesting a role in the hair cycle of dogs.
Surprisingly, most of the selected stem cell-associated genes did not show significantly different expression between anagen and telogen. These results are similar to results from a recent study comparing anagen and telogen HFs in Cashmere goats. In this study, only LHX2 and NFATC1 were differentially expressed [42]. In our study, NFATC1 was expressed in both late anagen and telogen, whereas LHX2 was expressed strongly in telogen only (Table 1). Since stem cell activity is highest in early anagen follicles, we assume that the comparison of late anagen and telogen is not optimal to detect differences in stem cell activity.
All genes playing a role in the immune function of the HFs, including β2 microglobulin (B2M) and canine MHC class I molecules (DLA88, DLA-79, DLA-12, DLA-64) were expressed at very low levels in late anagen HFs and much higher in telogen HFs and the IFE (Supplementary Table S2). This is in line with the knowledge that the inferior portion of the anagen HF is one of the immune-privileged sites of the body. Immune-privilege is achieved by, for example, the low expression of MHC class I molecules and B2M [103] and shields the HFs from autoimmune attacks, as happens in alopecia areata [104,105], a common disease in humans and domestic animals. Our data thus strongly indicate, that, like in humans, the inferior portion of the anagen HF is an immune-privileged site and identify the molecules associated with this privilege. In addition, our results indicate that dogs have a similar gene expression pattern associated with follicular immunoregulation in late anagen and, thus, the dog may be a suitable model to study the pathogenesis of alopecia areata.
Comparing our datasets with genes that have been associated with inherited canine epidermal disease and alopecia in the past, we identified the expression of all of these genes, with the exception of RSPO in the IFE, in all three tissue groups. Why genetic variants result in pathology only in the epidermis needs to be further investigated. Generally speaking, the genes identified in conjunction with alopecia are expressed at much lower levels. Interestingly, even though ichthyosis usually exclusively affects the epidermis, one of the genes associated with ichthyosis occurring in Great Danes (FATP4 [71]) had a higher expression in the anagen HF group compared to the IFE samples in our data set. This form of ichthyosis is the only one described in the literature that also has pathologic lesions in the folliculosebaceous units [106], which could be explained by the expression of FATP4 in anagen HFs in our data.

5. Conclusions

We presented the gene expression pattern of canine IFE and microdissected HFs in different hair cycle stages (late anagen and telogen). We found that the gene expression patterns of skin-associated genes in dogs are more similar to those in humans as compared to mice, further corroborating that the dog is a good model for human skin diseases. We also found gene expression that differs from what is known in mice and humans, such as BMP2 expression mainly in telogen in dogs but not mice, and KRT17 expression in the IFE in dogs, while absent in mice and humans. Our data provide the basis for research investigating the structure and function of canine skin or skin disease, but also for scientists using a canine model for human disease, by assigning genes to specific tissue types.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4425/11/8/884/s1, Table S1. Dogs used as tissue donors for the experiments. Dogs used as skin tissue donors for the experiments, including breed, sex and age, are listed. Table S2. Overall gene expression. Normalized read counts of all genes expressed in the different hair follicle stages and the epidermis. Gene expression values of individual dogs and average values for hair follicles in late anagen, early anagen, telogen and interfollicular epidermis sample groups are shown. Table S3. Results of differential expression analysis. Significantly differentially expressed genes comparing the base mean of normalized read counts of microdissected late anagen and telogen hair follicles with a false discovery rate (FDR) < 0.01 and log2 fold change of 1.5 as threshold, sorted alphabetically; BaseMean: mean of normalized read counts across all samples; log2 fold change: condition anagen versus telogen effect; LfcSE: standard error of the log2FoldChange; stat: the log2FoldChange divided by lfcSE; pvalue: Wald test effect; padj: Benjamini Hochberg adjusted p-values; FDR refers to Benjamini Hochberg adjusted p-values. Table S4. Skin-associated genes in humans and mice not expressed in dog skin. The top 100 human and mouse skin-specific genes [59] were compared to expression of these genes in dog skin. Genes highly expressed in human and mouse but not dog skin are listed. Figure S1. Microdissected hair follicles in different hair cycle stages. Hair follicles of both hair cycle stages were cut above the sebaceous glands to clearly separate hair follicles from interfollicular epidermis. A: Late anagen hair follicle. Note that the hair shaft emerges through the hair shaft ostia. B: Telogen hair follicle. Figure S2. Principal component analysis (PCA) in the two-dimensional space of the transcriptome samples. Samples are plotted across the most variable components (PC1 and PC2) and sample clustering is based on late anagen, telogen and interfollicular epidermis. Whereas samples of the late anagen and the interfollicular epidermis group cluster nicely together two telogen samples are clustering separately from the other telogen samples. This leads to the conclusion that these samples represent rather an early anagen stage. Figure S3. Principal component analysis (PCA) in the two-dimensional space of the transcriptome samples. Samples are plotted across the most variable components (PC1 and PC2) and sample clustering is based on age. Figure S4. Principal component analysis (PCA) in the two-dimensional space of the transcriptome samples. Samples are plotted across the most variable components (PC1 and PC2) and sample clustering is based on breed. Figure S5. Functional and pathway enrichment analysis of differentially expressed genes comparing late anagen and telogen samples. Gene Ontology enrichment analysis of differentially expressed genes has been performed using the Database for Annotation Visualization and Integrated Discovery. Enriched pathways and gene ontology terms that are upregulated or downregulated in anagen hair follicles are listed and their biological function is indicated in different colors.

Author Contributions

Conceptualization: M.M.W.; Methodology: M.A.T.B., K.R.G., T.L. and D.J.W.; Software: K.R.G., V.J. and T.L.; Validation: D.J.W., K.R.G. and V.J.; Formal analysis: M.A.T.B., K.R.G. and V.J; Investigation: D.J.W., M.A.T.B. and K.R.G.; Resources: D.J.W. and M.M.W.; Data curation: M.A.T.B., K.R.G., D.J.W. and V.J.; Writing—original draft: D.J.W.; Writing—review and editing: M.A.T.B., K.R.G., V.J, T.L. and M.M.W.; Visualization: D.J.W., K.R.G and M.M.W.; Supervision: D.J.W. and M.M.W.; Project administration: D.J.W. and M.M.W.; Funding acquisition: D.J.W. and M.M.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The authors would like to thank Alen Merzdo and Joshua Freeman for helping with the skin sampling. We thank Caitlin Older, Alex Myers and Kranti Konganti for helping with data analysis. We thank Muriel Fragnière, Sabrina Schenk and Nathalie Besuchet-Schmutz of the Next Generation Sequencing Platform of the University of Bern for performing sequencing experiments and the Interfaculty Bioinformatics Unit and the HPC cluster UBELIX (http://www.id.unibe.ch/hpc) at the University of Bern for providing computational infrastructure.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Fuchs, E. Scratching the surface of skin development. Nature 2007, 445, 834–842. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Suter, M.M.; Schulze, K.; Bergman, W.; Welle, M.; Roosje, P.; Muller, E.J. The keratinocyte in epidermal renewal defence. Vet. Dermatol. 2009, 20, 515–532. [Google Scholar] [CrossRef] [PubMed]
  3. Halprin, K.M. Epidermal turnover time—A re-examination. Br. J. Dermatol. 1972, 86, 14–19. [Google Scholar] [CrossRef]
  4. Blanpain, C.; Fuchs, E. Epidermal stem cells of the skin. Annu. Rev. Cell Dev. Biol. 2006, 22, 339–373. [Google Scholar] [CrossRef] [Green Version]
  5. Watt, F.M. The stem cell compartment in human interfollicular epidermis. J. Dermatol. Sci. 2002, 28, 173–180. [Google Scholar] [CrossRef]
  6. Watt, F.M.; Fujiwara, H. Cell-Extracellular matrix interactions in normal and diseased skin. Cold Spring Harb. Perspect. Biol. 2011, 3, a005124. [Google Scholar] [CrossRef] [Green Version]
  7. Hsu, Y.C.; Pasolli, H.A.; Fuchs, E. Dynamics between stem cells, niche, and progeny in the hair follicle. Cell 2011, 144, 92–105. [Google Scholar] [CrossRef] [Green Version]
  8. Zhang, Y.V.; Cheong, J.; Ciapurin, N.; McDermitt, D.J.; Tumbar, T. Distinct self-renewal and differentiation phases in the niche of infrequently dividing hair follicle stem cells. Cell Stem Cell 2009, 5, 267–278. [Google Scholar] [CrossRef] [Green Version]
  9. Trempus, C.S.; Morris, R.J.; Bortner, C.D.; Cotsarelis, G.; Faircloth, R.S.; Reece, J.M.; Tennant, R.W. Enrichment for living murine keratinocytes from the hair follicle bulge with the cell surface marker CD34. J. Investig. Dermatol. 2003, 120, 501–511. [Google Scholar]
  10. Gerhards, N.M.; Sayar, B.S.; Origgi, F.C.; Galichet, A.; Müller, E.J.; Welle, M.M.; Wiener, D.J. Stem Cell-Associated Marker Expression in Canine Hair Follicles. J. Histochem. Cytochem. 2016, 64, 190–204. [Google Scholar] [CrossRef] [Green Version]
  11. Liu, Y.; Lyle, S.; Yang, Z.; Cotsarelis, G. Keratin 15 promoter targets putative epithelial stem cells in the hair follicle bulge. J. Investig. Dermatol. 2003, 121, 963–968. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Jaks, V.; Barker, N.; Kasper, M.; van Es, J.H.; Snippert, H.J.; Clevers, H.; Toftgård, R. Lgr5 marks cycling, yet long-lived, hair follicle stem cells. Nat. Genet. 2008, 40, 1291–1299. [Google Scholar] [CrossRef] [PubMed]
  13. Nowak, J.A.; Polak, L.; Pasolli, H.A.; Fuchs, E. Hair follicle stem cells are specified and function in early skin morphogenesis. Cell Stem Cell 2008, 3, 33–43. [Google Scholar] [CrossRef] [Green Version]
  14. Snippert, H.J.; Haegebarth, A.; Kasper, M.; Jaks, V.; van Es, J.H.; Barker, N.; van de Wetering, M.; van den Born, M.; Begthel, H.; Vries, R.G.; et al. Lgr6 marks stem cells in the hair follicle that generate all cell lineages of the skin. Science 2010, 327, 1385–1389. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Jang, B.G.; Lee, C.; Kim, H.S.; Shin, M.S.; Cheon, M.S.; Kim, J.W.; Kim, W.H. Distinct expression profile of stem cell markers, LGR5 and LGR6, in basaloid skin tumors. Virchows Arch. 2017, 470, 301–310. [Google Scholar] [CrossRef] [PubMed]
  16. Kobayashi, T.; Iwasaki, T.; Amagai, M.; Ohyama, M. Canine follicle stem cell candidates reside in the bulge and share characteristic features with human bulge cells. J. Investig. Dermatol. 2010, 130, 1988–1995. [Google Scholar] [CrossRef] [Green Version]
  17. Brunner, M.A.T.; Jagannathan, V.; Waluk, D.P.; Roosje, P.; Linek, M.; Panakova, L.; Leeb, T.; Wiener, D.J.; Welle, M.M. Novel insights into the pathways regulating the canine hair cycle and their deregulation in alopecia X. PLoS ONE 2017, 12, e0186469. [Google Scholar] [CrossRef] [Green Version]
  18. de Castro, R.V.G.; Tavares, M.R.; Bressan, F.F.; Pieri, N.C.G.; Hill, A.B.T.; Souza, A.F.; Cruz, N.d.R.N.; Martins, D.S.; Ambrosio, C.E.; Meirelles, F.V.; et al. In vitro identification of a stem cell population from canine hair follicle bulge region. Tissue Cell 2018, 50, 43–50. [Google Scholar] [CrossRef] [Green Version]
  19. Jahoda, C.A.; Christiano, A.M. Niche crosstalk: Intercellular signals at the hair follicle. Cell 2011, 146, 678–681. [Google Scholar] [CrossRef] [Green Version]
  20. Plikus, M.V.; Chuong, C.M. Complex hair cycle domain patterns and regenerative hair waves in living rodents. J. Investig. Dermatol. 2008, 128, 1071–1080. [Google Scholar] [CrossRef] [Green Version]
  21. Festa, E.; Fretz, J.; Berry, R.; Schmidt, B.; Rodeheffer, M.; Horowitz, M.; Horsley, V. Adipocyte lineage cells contribute to the skin stem cell niche to drive hair cycling. Cell 2011, 146, 761–771. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Lee, J.; Tumbar, T. Hairy tale of signaling in hair follicle development and cycling. Semin. Cell Dev. Biol. 2012, 23, 906–916. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Rompolas, P.; Greco, V. Stem cell dynamics in the hair follicle niche. Semin. Cell Dev. Biol. 2014, 25, 34–42. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Sumigray, K.D.; Lechler, T. Cell adhesion in epidermal development and barrier formation. Curr. Top. Dev. Biol. 2015, 112, 383–414. [Google Scholar] [PubMed] [Green Version]
  25. Coulombe, P.A.; Kopan, R.; Fuchs, E. Expression of keratin K14 in the epidermis and hair follicle: Insights into complex programs of differentiation. J. Cell Biol. 1989, 109, 2295–2312. [Google Scholar] [CrossRef]
  26. Romani, N.; Ebner, S.; Tripp, C.H.; Flacher, V.; Koch, F.; Stoitzner, P. Epidermal Langerhans Cells—Changing views on their function in vivo. Immunol. Lett. 2006, 106, 119–125. [Google Scholar] [CrossRef]
  27. Klicznik, M.M.; Szenes-Nagy, A.B.; Campbell, D.J.; Gratz, I.K. Taking the lead-how keratinocytes orchestrate skin T cell immunity. Immunol. Lett. 2018, 200, 43–51. [Google Scholar] [CrossRef]
  28. Adachi, T.; Kobayashi, T.; Sugihara, E.; Yamada, T.; Ikuta, K.; Pittaluga, S.; Saya, H.; Amagai, M.; Nagao, K. Hair follicle-derived IL-7 and IL-15 mediate skin-resident memory T cell homeostasis and lymphoma. Nat. Med. 2015, 21, 1272–1279. [Google Scholar] [CrossRef]
  29. Nagao, K.; Kobayashi, T.; Moro, K.; Ohyama, M.; Adachi, T.; Kitashima, D.Y.; Ueha, S.; Horiuchi, K.; Tanizaki, H.; Kabashima, K.; et al. Stress-induced production of chemokines by hair follicles regulates the trafficking of dendritic cells in skin. Nat. Immunol. 2012, 13, 744–752. [Google Scholar] [CrossRef] [Green Version]
  30. Vahlquist, A.; Fischer, J.; Torma, H. Inherited Nonsyndromic Ichthyoses: An Update on Pathophysiology, Diagnosis and Treatment. Am. J. Clin. Dermatol. 2018, 19, 51–66. [Google Scholar] [CrossRef]
  31. Duverger, O.; Morasso, M.I. To grow or not to grow: Hair morphogenesis and human genetic hair disorders. Semin. Cell Dev. Biol. 2014, 25–26, 22–33. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Ahmed, A.; Almohanna, H.; Griggs, J.; Tosti, A. Genetic Hair Disorders: A Review. Dermatol. Ther. (Heidelb.) 2019, 9, 421–448. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Leeb, T.; Müller, E.J.; Roosje, P.; Welle, M. Genetic testing in veterinary dermatology. Vet. Dermatol. 2017, 28, 4-e1. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Bauer, A.; Nimmo, J.; Newman, R.; Brunner, M.; Welle, M.M.; Jagannathan, V.; Leeb, T. A splice site variant in the SUV39H2 gene in Greyhounds with nasal parakeratosis. Anim. Genet. 2018, 49, 137–140. [Google Scholar] [CrossRef] [Green Version]
  35. Drogemuller, M.; Jagannathan, V.; Becker, D.; Drogemuller, C.; Schelling, C.; Plassais, J.; Kaerle, C.; de Citres, C.D.; Thomas, A.; Muller, E.J.; et al. A mutation in the FAM83G gene in dogs with hereditary footpad hyperkeratosis (HFH). PLoS Genet. 2014, 10, e1004370. [Google Scholar] [CrossRef]
  36. Brunner, M.A.T.; Rufenacht, S.; Bauer, A.; Erpel, S.; Buchs, N.; Braga-Lagache, S.; Heller, M.; Leeb, T.; Jagannathan, V.; Wiener, D.J.; et al. Bald thigh syndrome in sighthounds-Revisiting the cause of a well-known disease. PLoS ONE 2019, 14, e0212645. [Google Scholar] [CrossRef]
  37. Joost, S.; Zeisel, A.; Jacob, T.; Sun, X.; La Manno, G.; Lonnerberg, P.; Linnarsson, S.; Kasper, M. Single-Cell Transcriptomics Reveals that Differentiation and Spatial Signatures Shape Epidermal and Hair Follicle Heterogeneity. Cell Syst. 2016, 3, 221–237.e9. [Google Scholar] [CrossRef] [Green Version]
  38. Acosta, J.R.; Joost, S.; Karlsson, K.; Ehrlund, A.; Li, X.; Aouadi, M.; Kasper, M.; Arner, P.; Ryden, M.; Laurencikiene, J. Single cell transcriptomics suggest that human adipocyte progenitor cells constitute a homogeneous cell population. Stem Cell Res. Ther. 2017, 8, 250. [Google Scholar] [CrossRef] [Green Version]
  39. Edqvist, P.H.; Fagerberg, L.; Hallstrom, B.M.; Danielsson, A.; Edlund, K.; Uhlen, M.; Ponten, F. Expression of human skin-specific genes defined by transcriptomics and antibody-based profiling. J. Histochem. Cytochem. 2015, 63, 129–141. [Google Scholar] [CrossRef] [Green Version]
  40. Takahashi, R.; Grzenda, A.; Allison, T.F.; Rawnsley, J.; Balin, S.J.; Sabri, S.; Plath, K.; Lowry, W.E. Defining Transcriptional Signatures of Human Hair Follicle Cell States. J. Investig. Dermatol. 2020, 140, 764–773.e4. [Google Scholar] [CrossRef]
  41. Gao, Y.; Wang, X.; Yan, H.; Zeng, J.; Ma, S.; Niu, Y.; Zhou, G.; Jiang, Y.; Chen, Y. Comparative Transcriptome Analysis of Fetal Skin Reveals Key Genes Related to Hair Follicle Morphogenesis in Cashmere Goats. PLoS ONE 2016, 11, e0151118. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. He, N.; Su, R.; Wang, Z.; Zhang, Y.; Li, J. Exploring differentially expressed genes between anagen and telogen secondary hair follicle stem cells from the Cashmere goat (Capra hircus) by RNA-Seq. PLoS ONE 2020, 15, e0231376. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Wang, Y.; Zheng, Y.; Guo, D.; Zhang, X.; Guo, S.; Hui, T.; Yue, C.; Sun, J.; Guo, S.; Bai, Z.; et al. m6A Methylation Analysis of Differentially Expressed Genes in Skin Tissues of Coarse and Fine Type Liaoning Cashmere Goats. Front. Genet. 2019, 10, 1318. [Google Scholar] [CrossRef] [PubMed]
  44. Su, R.; Fan, Y.; Qiao, X.; Li, X.; Zhang, L.; Li, C.; Li, J. Transcriptomic analysis reveals critical genes for the hair follicle of Inner Mongolia cashmere goat from catagen to telogen. PLoS ONE 2018, 13, e0204404. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Hoffman, J.M.; Creevy, K.E.; Franks, A.; O’Neill, D.G.; Promislow, D.E.L. The companion dog as a model for human aging and mortality. Aging Cell 2018, 17, e12737. [Google Scholar] [CrossRef]
  46. Martel, B.C.; Lovato, P.; Baumer, W.; Olivry, T. Translational Animal Models of Atopic Dermatitis for Preclinical Studies. Yale J. Biol. Med. 2017, 90, 389–402. [Google Scholar]
  47. Kobayashi, T.; Shimizu, A.; Nishifuji, K.; Amagai, M.; Iwasaki, T.; Ohyama, M. Canine hair-follicle keratinocytes enriched with bulge cells have the highly proliferative characteristic of stem cells. Vet. Dermatol. 2009, 20, 338–346. [Google Scholar] [CrossRef]
  48. Wiener, D.J.; Doherr, M.G.; Müller, E.J.; Welle, M.M. Spatial Distribution of Stem Cell-Like Keratinocytes in Dissected Compound Hair Follicles of the Dog. PLoS ONE 2016, 11, e0146937. [Google Scholar] [CrossRef] [Green Version]
  49. Müntener, T.; Doherr, M.G.; Guscetti, F.; Suter, M.M.; Welle, M.M. The canine hair cycle—A guide for the assessment of morphological and immunohistochemical criteria. Vet. Dermatol. 2011, 22, 383–395. [Google Scholar] [CrossRef] [Green Version]
  50. European Nucleotide Archive (ENA). Accession Number PRJEB21761. Available online: http://www.ebi.ac.uk/ena/data/view/PRJEB21761 (accessed on 26 June 2020).
  51. Andrews, S. FastQC: A Quality Control Tool for High Throughput Sequence Data. 2010. Available online: http://www.bioinformatics.babraham.ac.uk/projects/fastqc/ (accessed on 19 July 2020).
  52. Love, M.I.; Huber, W.; Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014, 15, 550. [Google Scholar] [CrossRef] [Green Version]
  53. Anders, S.; Pyl, P.T.; Huber, W. HTSeq—A Python framework to work with high-throughput sequencing data. Bioinformatics 2015, 31, 166–169. [Google Scholar] [CrossRef] [PubMed]
  54. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2019. [Google Scholar]
  55. Wickham, H. ggplot2—Elegant Graphics for Data Analysis, 2nd ed.; Springer International Publishing: New York, NY, USA, 2016. [Google Scholar]
  56. National Center for Biotechnology Information (NCBI). Gene. Available online: https://www.ncbi.nlm.nih.gov/gene (accessed on 19 July 2020).
  57. da Huang, W.; Sherman, B.T.; Lempicki, R.A. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat. Protoc. 2009, 4, 44–57. [Google Scholar] [CrossRef] [PubMed]
  58. Database for Annotation Visualization and Integrated Discovery (DAVID). Available online: https://david.ncifcrf.gov/ (accessed on 19 June 2020).
  59. Pirooznia, M.; Nagarajan, V.; Deng, Y. GeneVenn—A web application for comparing gene lists using Venn diagrams. Bioinformation 2007, 1, 420–422. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  60. Bauer, A.; De Lucia, M.; Jagannathan, V.; Mezzalira, G.; Casal, M.L.; Welle, M.M.; Leeb, T. A Large Deletion in the NSDHL Gene in Labrador Retrievers with a Congenital Cornification Disorder. G3 (Bethesda) 2017, 7, 3115–3121. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  61. Leuthard, F.; Lehner, G.; Jagannathan, V.; Leeb, T.; Welle, M. A missense variant in the NSDHL gene in a Chihuahua with a congenital cornification disorder resembling inflammatory linear verrucous epidermal nevi. Anim. Genet. 2019, 50, 768–771. [Google Scholar] [CrossRef]
  62. Linek, M.; Doelle, M.; Leeb, T.; Bauer, A.; Leuthard, F.; Henkel, J.; Bannasch, D.; Jagannathan, V.; Welle, M.M. ATP2A2 SINE Insertion in an Irish Terrier with Darier Disease and Associated Infundibular Cyst Formation. Genes 2020, 11, 481. [Google Scholar] [CrossRef]
  63. Baldeschi, C.; Gache, Y.; Rattenholl, A.; Bouillé, P.; Danos, O.; Ortonne, J.P.; Bruckner-Tuderman, L.; Meneguzzi, G. Genetic correction of canine dystrophic epidermolysis bullosa mediated by retroviral vectors. Hum. Mol. Genet. 2003, 12, 1897–1905. [Google Scholar] [CrossRef] [Green Version]
  64. Niskanen, J.; Dillard, K.; Arumilli, M.; Salmela, E.; Anttila, M.; Lohi, H.; Hytönen, M.K. Nonsense variant in COL7A1 causes recessive dystrophic epidermolysis bullosa in Central Asian Shepherd dogs. PLoS ONE 2017, 12, e0177527. [Google Scholar] [CrossRef]
  65. Capt, A.; Spirito, F.; Guaguere, E.; Spadafora, A.; Ortonne, J.P.; Meneguzzi, G. Inherited junctional epidermolysis bullosa in the German Pointer: Establishment of a large animal model. J. Investig. Dermatol. 2005, 124, 530–535. [Google Scholar] [CrossRef] [Green Version]
  66. Olivry, T.; Linder, K.E.; Wang, P.; Bizikova, P.; Bernstein, J.A.; Dunston, S.M.; Paps, J.S.; Casal, M.L. Deficient plakophilin-1 expression due to a mutation in PKP1 causes ectodermal dysplasia-skin fragility syndrome in Chesapeake Bay retriever dogs. PLoS ONE 2012, 7, e32072. [Google Scholar] [CrossRef] [Green Version]
  67. Mauldin, E.A.; Wang, P.; Olivry, T.; Henthorn, P.S.; Casal, M.L. Epidermolysis bullosa simplex in sibling Eurasier dogs is caused by a PLEC non-sense variant. Vet. Dermatol. 2017, 28, 10–e13. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  68. Jagannathan, V.; Bannoehr, J.; Plattet, P.; Hauswirth, R.; Drogemuller, C.; Drogemuller, M.; Wiener, D.J.; Doherr, M.; Owczarek-Lipska, M.; Galichet, A.; et al. A mutation in the SUV39H2 gene in Labrador Retrievers with hereditary nasal parakeratosis (HNPK) provides insights into the epigenetics of keratinocyte differentiation. PLoS Genet. 2013, 9, e1003848. [Google Scholar] [CrossRef] [PubMed]
  69. Bannoehr, J.; Balmer, P.; Stoffel, M.H.; Jagannathan, V.; Gaschen, V.; Kühni, K.; Sayar, B.; Drögemüller, M.; Howald, D.; Wiener, D.J.; et al. Abnormal keratinocyte differentiation in the nasal planum of Labrador Retrievers with hereditary nasal parakeratosis (HNPK). PLoS ONE 2020, 15, e0225901. [Google Scholar] [CrossRef] [PubMed]
  70. Bauer, A.; Waluk, D.P.; Galichet, A.; Timm, K.; Jagannathan, V.; Sayar, B.S.; Wiener, D.J.; Dietschi, E.; Muller, E.J.; Roosje, P.; et al. A de novo variant in the ASPRV1 gene in a dog with ichthyosis. PLoS Genet. 2017, 13, e1006651. [Google Scholar] [CrossRef]
  71. Metzger, J.; Wöhlke, A.; Mischke, R.; Hoffmann, A.; Hewicker-Trautwein, M.; Küch, E.M.; Naim, H.Y.; Distl, O. A Novel SLC27A4 Splice Acceptor Site Mutation in Great Danes with Ichthyosis. PLoS ONE 2015, 10, e0141514. [Google Scholar] [CrossRef]
  72. Credille, K.M.; Barnhart, K.F.; Minor, J.S.; Dunstan, R.W. Mild recessive epidermolytic hyperkeratosis associated with a novel keratin 10 donor splice-site mutation in a family of Norfolk terrier dogs. Br. J. Dermatol. 2005, 153, 51–58. [Google Scholar] [CrossRef] [Green Version]
  73. Mauldin, E.A.; Wang, P.; Evans, E.; Cantner, C.A.; Ferracone, J.D.; Credille, K.M.; Casal, M.L. Autosomal Recessive Congenital Ichthyosis in American Bulldogs Is Associated With NIPAL4 (ICHTHYIN) Deficiency. Vet. Pathol. 2015, 52, 654–662. [Google Scholar] [CrossRef] [Green Version]
  74. Casal, M.L.; Wang, P.; Mauldin, E.A.; Lin, G.; Henthorn, P.S. A Defect in NIPAL4 Is Associated with Autosomal Recessive Congenital Ichthyosis in American Bulldogs. PLoS ONE 2017, 12, e0170708. [Google Scholar] [CrossRef]
  75. Briand, A.; Cochet-Faivre, N.; Reyes-Gomez, E.; Jaraud-Darnault, A.; Tiret, L.; Chevallier, L. NIPAL4 deletion identified in an American Bully with autosomal recessive congenital ichthyosis and response to topical therapy. Vet. Med. Sci. 2019, 5, 112–117. [Google Scholar] [CrossRef] [Green Version]
  76. Grall, A.; Guaguère, E.; Planchais, S.; Grond, S.; Bourrat, E.; Hausser, I.; Hitte, C.; Le Gallo, M.; Derbois, C.; Kim, G.J.; et al. PNPLA1 mutations cause autosomal recessive congenital ichthyosis in golden retriever dogs and humans. Nat. Genet. 2012, 44, 140–147. [Google Scholar] [CrossRef]
  77. Credille, K.M.; Minor, J.S.; Barnhart, K.F.; Lee, E.; Cox, M.L.; Tucker, K.A.; Diegel, K.L.; Venta, P.J.; Hohl, D.; Huber, M.; et al. Transglutaminase 1-deficient recessive lamellar ichthyosis associated with a LINE-1 insertion in Jack Russell terrier dogs. Br. J. Dermatol. 2009, 161, 265–272. [Google Scholar] [CrossRef] [PubMed]
  78. Drogemuller, C.; Karlsson, E.K.; Hytonen, M.K.; Perloski, M.; Dolf, G.; Sainio, K.; Lohi, H.; Lindblad-Toh, K.; Leeb, T. A mutation in hairless dogs implicates FOXI3 in ectodermal development. Science 2008, 321, 1462. [Google Scholar] [CrossRef]
  79. Casal, M.L.; Scheidt, J.L.; Rhodes, J.L.; Henthorn, P.S.; Werner, P. Mutation identification in a canine model of X-linked ectodermal dysplasia. Mamm. Genome 2005, 16, 524–531. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  80. Waluk, D.P.; Zur, G.; Kaufmann, R.; Welle, M.M.; Jagannathan, V.; Drogemuller, C.; Muller, E.J.; Leeb, T.; Galichet, A. A Splice Defect in the EDA Gene in Dogs with an X-Linked Hypohidrotic Ectodermal Dysplasia (XLHED) Phenotype. G3 2016, 6, 2949–2954. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  81. Rasouliha, S.H.; Bauer, A.; Dettwiler, M.; Welle, M.M.; Leeb, T. A frameshift variant in the EDA gene in Dachshunds with X-linked hypohidrotic ectodermal dysplasia. Anim. Genet. 2018, 49, 651–654. [Google Scholar] [CrossRef]
  82. Parker, H.G.; Whitaker, D.T.; Harris, A.C.; Ostrander, E.A. Whole Genome Analysis of a Single Scottish Deerhound Dog Family Provides Independent Corroboration That a SGK3 Coding Variant Leads to Hairlessness. G3 (Bethesda) 2020, 10, 293–297. [Google Scholar] [CrossRef] [Green Version]
  83. Hytönen, M.K.; Lohi, H. A frameshift insertion in SGK3 leads to recessive hairlessness in Scottish Deerhounds: A candidate gene for human alopecia conditions. Hum. Genet. 2019, 138, 535–539. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  84. Parker, H.G.; Chase, K.; Cadieu, E.; Lark, K.G.; Ostrander, E.A. An insertion in the RSPO2 gene correlates with improper coat in the Portuguese water dog. J. Hered. 2010, 101, 612–617. [Google Scholar] [CrossRef] [Green Version]
  85. Gerber, P.A.; Buhren, B.A.; Schrumpf, H.; Homey, B.; Zlotnik, A.; Hevezi, P. The top skin-associated genes: A comparative analysis of human and mouse skin transcriptomes. Biol. Chem. 2014, 395, 577–591. [Google Scholar] [CrossRef]
  86. Hawkshaw, N.J.; Hardman, J.A.; Alam, M.; Jimenez, F.; Paus, R. Deciphering the molecular morphology of the human hair cycle: Wnt signalling during the telogen-anagen transformation. Br. J. Dermatol. 2019, 182, 1184–1193. [Google Scholar] [CrossRef]
  87. Reddy, S.; Andl, T.; Bagasra, A.; Lu, M.M.; Epstein, D.J.; Morrisey, E.E.; Millar, S.E. Characterization of Wnt gene expression in developing and postnatal hair follicles and identification of Wnt5a as a target of Sonic hedgehog in hair follicle morphogenesis. Mech. Dev. 2001, 107, 69–82. [Google Scholar] [CrossRef]
  88. Panteleyev, A.A.; Paus, R.; Wanner, R.; Nurnberg, W.; Eichmuller, S.; Thiel, R.; Zhang, J.; Henz, B.M.; Rosenbach, T. Keratin 17 gene expression during the murine hair cycle. J. Investig. Dermatol. 1997, 108, 324–329. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  89. Plikus, M.V.; Mayer, J.A.; de la Cruz, D.; Baker, R.E.; Maini, P.K.; Maxson, R.; Chuong, C.M. Cyclic dermal BMP signalling regulates stem cell activation during hair regeneration. Nature 2008, 451, 340–344. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  90. Su, R.; Zhang, W.; Sharma, R.; Chang, Z.; Yin, J.; Li, J. Characterization of BMP2 gene expression in embryonic and adult Inner Mongolia Cashmere goat (Capra hircus) hair follicles. Can. J. Anim. Sci. 2009, 89, 457–462. [Google Scholar] [CrossRef]
  91. McGowan, K.M.; Tong, X.; Colucci-Guyon, E.; Langa, F.; Babinet, C.; Coulombe, P.A. Keratin 17 null mice exhibit age- and strain-dependent alopecia. Genes Dev. 2002, 16, 1412–1422. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  92. Yang, L.; Zhang, S.; Wang, G. Keratin 17 in disease pathogenesis: From cancer to dermatoses. J. Pathol. 2019, 247, 158–165. [Google Scholar] [CrossRef]
  93. Wang, Z.; Wong, P.; Langbein, L.; Schweizer, J.; Coulombe, P.A. Type II epithelial keratin 6hf (K6hf) is expressed in the companion layer, matrix, and medulla in anagen-stage hair follicles. J. Investig. Dermatol. 2003, 121, 1276–1282. [Google Scholar] [CrossRef] [Green Version]
  94. Lim, X.; Nusse, R. Wnt signaling in skin development, homeostasis, and disease. Cold Spring Harb. Perspect. Biol. 2013, 5, a008029. [Google Scholar] [CrossRef] [Green Version]
  95. Du, G.; Kataoka, K.; Sakaguchi, M.; Abarzua, F.; Than, S.S.; Sonegawa, H.; Makino, T.; Shimizu, T.; Huh, N.H. Expression of REIC/Dkk-3 in normal and hyperproliferative epidermis. Exp. Dermatol. 2011, 20, 273–277. [Google Scholar] [CrossRef]
  96. Prowse, D.M.; Lee, D.; Weiner, L.; Jiang, N.; Magro, C.M.; Baden, H.P.; Brissette, J.L. Ectopic expression of the nude gene induces hyperproliferation and defects in differentiation: Implications for the self-renewal of cutaneous epithelia. Dev. Biol. 1999, 212, 54–67. [Google Scholar] [CrossRef] [Green Version]
  97. Jave-Suarez, L.F.; Winter, H.; Langbein, L.; Rogers, M.A.; Schweizer, J. HOXC13 is involved in the regulation of human hair keratin gene expression. J. Biol. Chem. 2002, 277, 3718–3726. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  98. Mecklenburg, L.; Tychsen, B.; Paus, R. Learning from nudity: Lessons from the nude phenotype. Exp. Dermatol. 2005, 14, 797–810. [Google Scholar] [CrossRef] [PubMed]
  99. Meier, N.; Dear, T.N.; Boehm, T. Whn mHa3 are components of the genetic hierarchy controlling hair follicle differentiation. Mech. Dev. 1999, 89, 215–221. [Google Scholar] [CrossRef]
  100. Kawano, M.; Komi-Kuramochi, A.; Asada, M.; Suzuki, M.; Oki, J.; Jiang, J.; Imamura, T. Comprehensive analysis of FGF and FGFR expression in skin: FGF18 is highly expressed in hair follicles and capable of inducing anagen from telogen stage hair follicles. J. Investig. Dermatol. 2005, 124, 877–885. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  101. Plikus, M.V.; Baker, R.E.; Chen, C.C.; Fare, C.; de la Cruz, D.; Andl, T.; Maini, P.K.; Millar, S.E.; Widelitz, R.; Chuong, C.M. Self-organizing and stochastic behaviors during the regeneration of hair stem cells. Science 2011, 332, 586–589. [Google Scholar] [CrossRef] [Green Version]
  102. Hou, C.; Miao, Y.; Wang, X.; Chen, C.; Lin, B.; Hu, Z. Expression of matrix metalloproteinases and tissue inhibitor of matrix metalloproteinases in the hair cycle. Exp. Ther. Med. 2016, 12, 231–237. [Google Scholar] [CrossRef] [Green Version]
  103. Paus, R.; Nickoloff, B.J.; Ito, T. A hairy privilege. Trends Immunol. 2005, 26, 32–40. [Google Scholar] [CrossRef]
  104. Paus, R.; Bertolini, M. The role of hair follicle immune privilege collapse in alopecia areata: Status and perspectives. J. Investig. Dermatol. Symp. Proc. 2013, 16, S25–S27. [Google Scholar] [CrossRef]
  105. Paus, R.; Slominski, A.; Czarnetzki, B.M. Is alopecia areata an autoimmune-response against melanogenesis-related proteins, exposed by abnormal MHC class I expression in the anagen hair bulb? Yale J. Biol. Med. 1993, 66, 541–554. [Google Scholar]
  106. Hoffmann, A.; Metzger, J.; Wöhlke, A.; Peters, M.; Junginger, J.; Mischke, R.; Distl, O.; Hewicker-Trautwein, M. Congenital Ichthyosis in 14 Great Dane Puppies With a New Presentation. Vet. Pathol. 2016, 53, 614–620. [Google Scholar] [CrossRef] [Green Version]
Genes 11 00884 g001 550
Figure 1. Venn Diagram depicting the most strongly expressed genes (average normalized counts > 5000) in telogen and late anagen canine hair follicles and the interfollicular epidermis, respectively. Genes with the highest expression in each group (late anagen, telogen and interfollicular epidermis) and those which are expressed in more than one group are outlined.
Figure 1. Venn Diagram depicting the most strongly expressed genes (average normalized counts > 5000) in telogen and late anagen canine hair follicles and the interfollicular epidermis, respectively. Genes with the highest expression in each group (late anagen, telogen and interfollicular epidermis) and those which are expressed in more than one group are outlined.
Genes 11 00884 g001
Genes 11 00884 g002 550
Figure 2. Heat map presenting hierarchical clustering of the 30 most significantly differentially expressed genes comparing late anagen and telogen hair cycle stages. High gene expression is shown in red, low expression in blue, intermediate DESeq2 normalized values in white and different shades of red and blue.
Figure 2. Heat map presenting hierarchical clustering of the 30 most significantly differentially expressed genes comparing late anagen and telogen hair cycle stages. High gene expression is shown in red, low expression in blue, intermediate DESeq2 normalized values in white and different shades of red and blue.
Genes 11 00884 g002
Genes 11 00884 g003 550
Figure 3. Heat map presenting hierarchical clustering of a preselected subset of 46 genes, based on their known function associated with hair cycle, hair follicle structure and function and stem cells, comparing late anagen and telogen hair cycle stages. High expression is shown in red, low expression in blue, intermediate DESeq2 normalized values in white and different shades of red and blue.
Figure 3. Heat map presenting hierarchical clustering of a preselected subset of 46 genes, based on their known function associated with hair cycle, hair follicle structure and function and stem cells, comparing late anagen and telogen hair cycle stages. High expression is shown in red, low expression in blue, intermediate DESeq2 normalized values in white and different shades of red and blue.
Genes 11 00884 g003
Table
Table 1. Expression levels of 92 selected genes. Preselected hair cycle-associated, hair follicle structure-associated, stem cell-associated, immune function-associated and epidermal structure-associated genes are listed, and the expression levels based on the average TPM normalized values in telogen, anagen and interfollicular epidermis samples are depicted. Darker colors represent higher expression levels (average normalized counts ≥ 400), lighter shades of blue represent intermediate and low expression levels (average normalized count of ≥100 and <400 and ≥10 and <100, respectively). White boxes represent very low expression of genes (average normalized counts < 10).
Table 1. Expression levels of 92 selected genes. Preselected hair cycle-associated, hair follicle structure-associated, stem cell-associated, immune function-associated and epidermal structure-associated genes are listed, and the expression levels based on the average TPM normalized values in telogen, anagen and interfollicular epidermis samples are depicted. Darker colors represent higher expression levels (average normalized counts ≥ 400), lighter shades of blue represent intermediate and low expression levels (average normalized count of ≥100 and <400 and ≥10 and <100, respectively). White boxes represent very low expression of genes (average normalized counts < 10).
Hair Cycle-Associated Genes
Gene SymbolFull Name of GeneAnagenTelogenIFE
ADCY7Activin A receptor type 118.643.775.48
BMP2Bone morphogenetic protein 244.37281.7966.23
BMP4Bone morphogenetic protein 499.167.2912.64
BMP7Bone morphogenetic protein 725.657.9722.11
BNC1Basonuclin 170.4759.6883.80
BNC2Basonuclin 248.3440.384.74
CCNA2Cyclin A229.046.119.87
CYP1A1Cytochrome P450 family 1 subfamily A member 117.3943.1874.82
CYP1B1Cytochrome P450 family 1 subfamily B member 196.94214.59129.88
CYP24A1Cytochrome P450 family 24 subfamily A member 118.890.120.00
CYP27B1Cytochrome P450 family 27 subfamily B member 157.641.450.96
DKK3Dickkopf 3102.28129.90218.67
FGF18Fibroblast growth factor 188.70133.544.28
FGF5Fibroblast growth factor 553.664.651.78
FOXN1Forkhead box protein N1182.115.5824.39
FZD10Frizzled 1015.663.343.48
FZD3Frizzled 316.469.322.64
HSD17B14Hydroxysteroid 17-beta dehydrogenase 14396.5817.2025.22
JAG1Jagged 130.0522.3613.30
KISS1KISS-1 metastasis-suppressor60.790.170.00
MMP7Matrix metalloproteinase 75.44218.7312.82
MSX2Homeobox protein MSX2771.0949.3864.14
SFRP1Secreted frizzled related protein 1107.4037.3118.18
SFRP2Secreted frizzled related protein 239.5314.070.26
SOX10Sex determining region Y-box 1054.9359.2712.94
TCF7/TCF1Transcription factor 155.8446.934.60
TGFB2Transforming growth factor beta 22.5311.325.62
THBS2Thrombospondin 267.86137.4962.01
VDRVitamin D receptor105.5855.5822.71
WNT3Wnt family member 322.821.9397.02
WNT5aWnt family member 5A36.115.1512.76
WNT10AWnt family member 10A27.3218.1220.25
Hair Follicle Structure-Associated Genes
Gene SymbolFull Name of GeneAnagenTelogenIFE
CD133/PROM1Prominin15.863.440.86
CD71/TFRCTransferrin receptor61.6131.8823.85
CDH3P-Cadherin 3196.09225.99463.31
GATA6GATA binding protein 614.8059.764.07
HOXC13Homeobox C3374.2241.3913.79
KRT16Keratin 165931.8210,187.551790.56
KRT17Keratin 1715,476.3777,099.968158.26
KRT25Keratin 2510,456.2733.582.57
KRT31Keratin 315422.802.051.89
KRT33AKeratin 33A10,416.912.247.39
KRT33BKeratin 33B10,528.671.772.25
KRT34Keratin 343966.390.550.81
KRT35Keratin 355365.721.110.36
KRT6AKeratin 6A1018.091513.1937.94
KRT6BKeratin 6B3692.268621.6280.91
KRT71Keratin 7117,621.0411.187.17
KRT75Keratin 52585.702895.722.59
KRT79Keratin 79279.523752.501216.86
KRT83Keratin 835671.732.302.07
KRT85Keratin 8510,133.202.732.81
KRT86Keratin 863068.601.201.06
KRTAP16-1Keratin associated protein 16-12984.340.330.72
LEF1Lymphoid enhancer binding factor 187.2718.8421.65
PMELPremelanosome protein1907.9021.59135.56
TCHHTrichohyalin3997.425.248.35
Stem Cell-Associated Genes
Gene SymbolFull Name of GeneAnagenTelogenIFE
CD200CD200 Molecule23.6933.446.71
CD34CD34 Molecule111.79131.5547.56
ITGA6Integrin subunit alpha 670.85189.30150.07
ITGB1Integrin subunit beta 1219.44341.60265.09
KRT15Keratin 152485.1217,673.044104.93
KRT19Keratin 1959.14180.624.04
LGR5Leucin rich repeat containing G protein-coupled receptor 59.4310.240.02
LHX2LIM homeobox 2105.49200.130.34
NFATC1Nuclear factor of activated T cells19.1676.5421.90
SOX9SRY-box transcription factor 9563.03817.93104.80
TCF3Transcription factor 3127.36105.9582.50
TCF4Transcription factor 442.90104.5340.37
Epidermal Structure-Associated Genes
Gene SymbolFull Name of GeneAnagenTelogenIFE
DSC1Desmocollin 1328.6771.791523.99
DSG1Desmoglein 1505.83479.42984.35
DSPDesmoplakin1406.581492.091610.92
FLGFilaggrin12.3571.7060.74
IVLInvolucrin2244.92480.5585.91
KRT1Keratin 1203.64988.4421,712.07
KRT5Keratin 513,982.2718,257.9412,733.06
KRT10Keratin 10339.132042.7533,847.68
KRT14Keratin 1423,040.6519,320.4929,472.32
LORLoricrin65.36399.387229.10
PLECPlectin55.7887.5695.43
TP63Tumor protein p63137.68179.75230.67
Immune Function-Associated Genes
Gene SymbolFull Name of GeneAnagenTelogenIFE
CCL20C-C motif chemokine ligand 200.783.9414.63
CD74CD74 molecule113.38737.231498.77
DLA88MHC class I DLA-88116.10439.94584.87
DLA-DQA1MHC class II DQ alpha 120.24115.09257.38
HLA-DRB1MHC class II DR beta 134.90211.64455.74
IL17BInterleukin 17B15.0111.494.35
IL18Interleukin 1818.5562.5194.44
IL1AInterleukin 1 alpha19.316.5913.99
IL33Interleukin 3310.4028.9117.75
TLR2Toll like receptor 21.7915.282.43
TLR4Toll like receptor 47.5517.1412.41
Table
Table 2. Expression of genes associated with canine genodermatoses. Genes that are associated with a genetic variant leading to epidermal or follicular disease or phenotype are listed and the expression levels based on the average TPM normalized values in telogen, anagen and interfollicular epidermis samples are depicted.
Table 2. Expression of genes associated with canine genodermatoses. Genes that are associated with a genetic variant leading to epidermal or follicular disease or phenotype are listed and the expression levels based on the average TPM normalized values in telogen, anagen and interfollicular epidermis samples are depicted.
Tissue AffectedDiseaseDog BreedGene MutationAnagenTelogenIFERef.
IFECHILD’s like syndromeLabrador retriever, ChihuahuaNSDHL41.236.521.2[60,61]
Darier diseaseIrish TerrierATP2A214.319.318.9[62]
Epidermolysis bullosa, dystrophicGolden retrieverCOL7A19.658.661.4[63,64]
Epidermolysis bullosa simplex, junctionalGerman shorthaired pointerLAMA320.658.641.4[65]
Epidermolysis bullosa simplex, suprabasalChesapeake Bay retrieverPKP10.10.20.3[66]
Epidermolysis bullosa simplex, basalEurasierPLEC55.887.695.4[67]
Hereditary nasal parakeratosisLabrador retriever, greyhoundSUV39H29.813.912.5[34,68,69]
Hereditary footpad hyperkeratosisIrish terrier, KromfohrländerFAM83G174.2190.6158.6[35]
IchthyosisGerman shepherd dogASPRV1252.9111.1698.8[70]
Great DaneSLC27A4 (FATP4)178.9125.9147.6[71]
Norfolk terrierKRT10339.12042.833,847.7[72]
American bulldogNIPAL460.968.2250.1[73,74,75]
Golden retrieverPNPLA158.046.890.3[76]
Jack Russell terrierTGM183.134.2138.0[77]
HFCanine ectodermal dysplasiaChinese crested, Mexican hairless dogs, Peruvian hairless dogFOXI30.47.90.1[78]
German shepherd dog, dachshundEDA2.47.84.8[79,80,81]
Chesapeake Bay retrieverPKP10.10.20.3[66]
HairlessnessAmerican terrier, Scottish deerhoundSGK374.138.147.9[82,83]
Improper coatPortuguese water dogRSPO1.20.60.0[84]

Share and Cite

MDPI and ACS Style

Wiener, D.J.; Groch, K.R.; Brunner, M.A.T.; Leeb, T.; Jagannathan, V.; Welle, M.M. Transcriptome Profiling and Differential Gene Expression in Canine Microdissected Anagen and Telogen Hair Follicles and Interfollicular Epidermis. Genes 2020, 11, 884. https://doi.org/10.3390/genes11080884

AMA Style

Wiener DJ, Groch KR, Brunner MAT, Leeb T, Jagannathan V, Welle MM. Transcriptome Profiling and Differential Gene Expression in Canine Microdissected Anagen and Telogen Hair Follicles and Interfollicular Epidermis. Genes. 2020; 11(8):884. https://doi.org/10.3390/genes11080884

Chicago/Turabian Style

Wiener, Dominique J., Kátia R. Groch, Magdalena A.T. Brunner, Tosso Leeb, Vidhya Jagannathan, and Monika M. Welle. 2020. "Transcriptome Profiling and Differential Gene Expression in Canine Microdissected Anagen and Telogen Hair Follicles and Interfollicular Epidermis" Genes 11, no. 8: 884. https://doi.org/10.3390/genes11080884

APA Style

Wiener, D. J., Groch, K. R., Brunner, M. A. T., Leeb, T., Jagannathan, V., & Welle, M. M. (2020). Transcriptome Profiling and Differential Gene Expression in Canine Microdissected Anagen and Telogen Hair Follicles and Interfollicular Epidermis. Genes, 11(8), 884. https://doi.org/10.3390/genes11080884

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop