BarkBase: Epigenomic Annotation of Canine Genomes

Large-scale genomic analysis of the domestic dog informs biological discovery

Recent advances in genomics, coupled with a unique population structure and remarkable levels of variation, have propelled the domestic dog to new levels as a system for understanding fundamental principles in mammalian biology. Central to this advance are more than 350 recognized breeds, each a closed population that has undergone selection for unique features. Genetic variation in the domestic dog is particularly well characterized compared with other domestic mammals, with almost 3000 high-coverage genomes publicly available. Importantly, as the number of sequenced genomes increases, new avenues for analysis are becoming available. Herein, we discuss recent discoveries in canine genomics regarding behavior, morphology, and disease susceptibility. We explore the limitations of current data sets for variant interpretation, tradeoffs between sequencing strategies, and the burgeoning role of long-read genomes for capturing structural variants. In addition, we consider how large-scale collections of whole-genome sequence data drive rare variant discovery and assess the geographic distribution of canine diversity, which identifies Asia as a major source of missing variation. Finally, we review recent comparative genomic analyses that will facilitate annotation of the noncoding genome in dogs.

Cis-regulatory evolution of the recently expanded Ly49 gene family

Comparative genomics has revealed the rapid expansion of multiple gene families involved in immunity. Members within each gene family often evolved distinct roles in immunity. However, less is known about the evolution of their epigenome and cis-regulation. Here we systematically profile the epigenome of the recently expanded murine Ly49 gene family that mainly encode either inhibitory or activating surface receptors on natural killer cells. We identify a set of cis-regulatory elements (CREs) for activating Ly49 genes. In addition, we show that in mice, inhibitory and activating Ly49 genes are regulated by two separate sets of proximal CREs, likely resulting from lineage-specific losses of CRE activity. Furthermore, we find that some Ly49 genes are cross-regulated by the CREs of other Ly49 genes, suggesting that the Ly49 family has begun to evolve a concerted cis-regulatory mechanism. Collectively, we demonstrate the different modes of cis-regulatory evolution for a rapidly expanding gene family.

Imprinting at the KBTBD6 locus involves species-specific maternal methylation and monoallelic expression in livestock animals

BACKGROUND

The primary differentially methylated regions (DMRs) which are maternally hypermethylated serve as imprinting control regions (ICRs) that drive monoallelic gene expression, and these ICRs have been investigated due to their implications in mammalian development. Although a subset of genes has been identified as imprinted, in-depth comparative approach needs to be developed for identification of species-specific imprinted genes. Here, we examined DNA methylation status and allelic expression at the KBTBD6 locus across species and tissues and explored potential mechanisms of imprinting.

RESULTS

Using whole-genome bisulfite sequencing and RNA-sequencing on parthenogenetic and normal porcine embryos, we identified a maternally hypermethylated DMR between the embryos at the KBTBD6 promoter CpG island and paternal monoallelic expression of KBTBD6. Also, in analyzed domesticated mammals but not in humans, non-human primates and mice, the KBTBD6 promoter CpG islands were methylated in oocytes and/or allelically methylated in tissues, and monoallelic KBTBD6 expression was observed, indicating livestock-specific imprinting. Further analysis revealed that these CpG islands were embedded within transcripts in porcine and bovine oocytes which coexisted with an active transcription mark and DNA methylation, implying the presence of transcription-dependent imprinting.

CONCLUSIONS

In this study, our comparative approach revealed an imprinted expression of the KBTBD6 gene in domesticated mammals, but not in humans, non-human primates, and mice which implicates species-specific evolution of genomic imprinting.

Genome sequencing of 2000 canids by the Dog10K consortium advances the understanding of demography, genome function and architecture

BACKGROUND

The international Dog10K project aims to sequence and analyze several thousand canine genomes. Incorporating 20 × data from 1987 individuals, including 1611 dogs (321 breeds), 309 village dogs, 63 wolves, and four coyotes, we identify genomic variation across the canid family, setting the stage for detailed studies of domestication, behavior, morphology, disease susceptibility, and genome architecture and function.

RESULTS

We report the analysis of > 48 M single-nucleotide, indel, and structural variants spanning the autosomes, X chromosome, and mitochondria. We discover more than 75% of variation for 239 sampled breeds. Allele sharing analysis indicates that 94.9% of breeds form monophyletic clusters and 25 major clades. German Shepherd Dogs and related breeds show the highest allele sharing with independent breeds from multiple clades. On average, each breed dog differs from the UU_Cfam_GSD_1.0 reference at 26,960 deletions and 14,034 insertions greater than 50 bp, with wolves having 14% more variants. Discovered variants include retrogene insertions from 926 parent genes. To aid functional prioritization, single-nucleotide variants were annotated with SnpEff and Zoonomia phyloP constraint scores. Constrained positions were negatively correlated with allele frequency. Finally, the utility of the Dog10K data as an imputation reference panel is assessed, generating high-confidence calls across varied genotyping platform densities including for breeds not included in the Dog10K collection.

CONCLUSIONS

We have developed a dense dataset of 1987 sequenced canids that reveals patterns of allele sharing, identifies likely functional variants, informs breed structure, and enables accurate imputation. Dog10K data are publicly available.

Integrative mapping of the dog epigenome: Reference annotation for comparative intertissue and cross-species studies

Dogs have become a valuable model in exploring multifaceted diseases and biology relevant to human health. Despite large-scale dog genome projects producing high-quality draft references, a comprehensive annotation of functional elements is still lacking. We addressed this through integrative next-generation sequencing of transcriptomes paired with five histone marks and DNA methylome profiling across 11 tissue types, deciphering the dog's epigenetic code by defining distinct chromatin states, super-enhancer, and methylome landscapes, and thus showed that these regions are associated with a wide range of biological functions and cell/tissue identity. In addition, we confirmed that the phenotype-associated variants are enriched in tissue-specific regulatory regions and, therefore, the tissue of origin of the variants can be traced. Ultimately, we delineated conserved and dynamic epigenomic changes at the tissue- and species-specific resolutions. Our study provides an epigenomic blueprint of the dog that can be used for comparative biology and medical research.

ICAnnoLncRNA: A Snakemake Pipeline for a Long Non-Coding-RNA Search and Annotation in Transcriptomic Sequences

Long non-coding RNAs (lncRNAs) are RNA molecules longer than 200 nucleotides that do not encode proteins. Experimental studies have shown the diversity and importance of lncRNA functions in plants. To expand knowledge about lncRNAs in other species, computational pipelines that allow for standardised data-processing steps in a mode that does not require user control up until the final result were actively developed recently. These advancements enable wider functionality for lncRNA data identification and analysis. In the present work, we propose the ICAnnoLncRNA pipeline for the automatic identification, classification and annotation of plant lncRNAs in assembled transcriptomic sequences. It uses the LncFinder software for the identification of lncRNAs and allows the adjustment of recognition parameters using genomic data for which lncRNA annotation is available. The pipeline allows the prediction of lncRNA candidates, alignment of lncRNA sequences to the reference genome, filtering of erroneous/noise transcripts and probable transposable elements, lncRNA classification by genome location, comparison with sequences from external databases and analysis of lncRNA structural features and expression. We used transcriptomic sequences from 15 maize libraries assembled by Trinity and Hisat2/StringTie to demonstrate the application of the ICAnnoLncRNA pipeline.

Identification of Genetic Risk Factors for Monogenic and Complex Canine Diseases

Advances in DNA sequencing and other technologies have greatly facilitated the identification of genetic risk factors for inherited diseases in dogs. We review recent technological developments based on selected examples from canine disease genetics. The identification of disease-causing variants in dogs with monogenic diseases may become a widely employed diagnostic approach in clinical veterinary medicine in the not-too-distant future. Diseases with complex modes of inheritance continue to pose challenges to researchers but have also become much more tangible than in the past. In addition to strategies for identifying genetic risk factors, we provide some thoughts on the interpretation of sequence variants that are largely inspired by developments in human clinical genetics.

Environmental risk factors in puppies and kittens for developing chronic disorders in adulthood: A call for research on developmental programming

Many dogs and cats are affected by chronic diseases that significantly impact their health and welfare and relationships with humans. Some of these diseases can be challenging to treat, and a better understanding of early-life risk factors for diseases occurring in adulthood is key to improving preventive veterinary care and husbandry practices. This article reviews early-life risk factors for obesity and chronic enteropathy, and for chronic behavioral problems, which can also be intractable with life-changing consequences. Aspects of early life in puppies and kittens that can impact the risk of adult disorders include maternal nutrition, establishment of the gut microbiome, maternal behavior, weaning, nutrition during growth, growth rate, socialization with conspecifics and humans, rehoming and neutering. Despite evidence in some species that the disorders reviewed here reflect the developmental origins of health and disease (DOHaD), developmental programming has rarely been studied in dogs and cats. Priorities and strategies to increase knowledge of early-life risk factors and DOHaD in dogs and cats are discussed. Critical windows of development are proposed: preconception, gestation, the suckling period, early growth pre-neutering or pre-puberty, and growth post-neutering or post-puberty to adult size, the durations of which depend upon species and breed. Challenges to DOHaD research in these species include a large number of breeds with wide genetic and phenotypic variability, and the existence of many mixed-breed individuals. Moreover, difficulties in conducting prospective lifelong cohort studies are exacerbated by discontinuity in pet husbandry between breeders and subsequent owners, and by the dispersed nature of pet ownership.

Bayesian model and selection signature analyses reveal risk factors for canine atopic dermatitis

Canine atopic dermatitis is an inflammatory skin disease with clinical similarities to human atopic dermatitis. Several dog breeds are at increased risk for developing this disease but previous genetic associations are poorly defined. To identify additional genetic risk factors for canine atopic dermatitis, we here apply a Bayesian mixture model adapted for mapping complex traits and a cross-population extended haplotype test to search for disease-associated loci and selective sweeps in four dog breeds at risk for atopic dermatitis. We define 15 associated loci and eight candidate regions under selection by comparing cases with controls. One associated locus is syntenic to the major genetic risk locus (Filaggrin locus) in human atopic dermatitis. One selection signal in common type Labrador retriever cases positions across the TBC1D1 gene (body weight) and one signal of selection in working type German shepherd controls overlaps the LRP1B gene (brain), near the KYNU gene (psoriasis). In conclusion, we identify candidate genes, including genes belonging to the same biological pathways across multiple loci, with potential relevance to the pathogenesis of canine atopic dermatitis. The results show genetic similarities between dog and human atopic dermatitis, and future across-species genetic comparisons are hereby further motivated.

Identification of regenerating island-derived protein 3E in dogs

Regenerating islet-derived protein (REG) 1A (aka pancreatic stone protein) and REG3A (aka pancreatitis-associated protein) are upregulated in humans with sepsis, pancreatitis, and gastrointestinal diseases, but little is known about this protein family in dogs. Our aim was to identify REG1 and REG3 family members in dogs. REG-family genes were computationally annotated in the canine genome and proteome, with verification of gene expression using publicly available RNA-seq data. The presence of the protein in canine pancreatic tissue and plasma was investigated with Western blot and immunohistochemistry, using anti-human REG1A and REG3A antibodies. Protein identity was confirmed with mass spectrometry. Two members of the REG3 subfamily were found in the canine genome, REG3E1 and REG3E2, both encoding for the same 176 AA protein, subsequently named REG3E. Anti-human REG3A antibodies demonstrated cross-reactivity with the canine REG3E protein in pancreas homogenates. In canine plasma, a protein band of approximately 17 kDa was apparent. Mass spectrometry confirmed this protein to be the product of the two annotated REG3E genes. Strong immunoreactivity to anti-human REG3A antibodies was found in sections of canine pancreas affected with acute pancreatitis, but it was weak in healthy pancreatic tissue. Recombinant canine REG3E protein underwent a selective trypsin digestion as described in other species. No evidence for the presence of a homolog of REG1A in dogs was found in any of the investigations. In conclusion, dogs express REG3E in the pancreas, whose role as biomarker merits further investigations. Homologs to human REG1A are not likely to exist in dogs.

Canine reference genome accuracy impacts variant calling: Lessons learned from investigating embryonic lethal variants

Deficient homozygosity of a variant maintained in a population suggests that the variant may be embryonic lethal. We examined whole genome sequence data from 675 canids to investigate for variants with missing homozygosity and high predicted impact. Our analysis identified 45 variants, in 32 genes. However, further scrutiny of the sequence reads revealed that all but one of these variants were artifacts of the variant calling process when using CanFam3.1, a widely utilized canine reference genome. We demonstrate that the use of multiple, newer reference genomes could reduce artifacts and lead to more accurate variant identification.

Ancestry-inclusive dog genomics challenges popular breed stereotypes

Behavioral genetics in dogs has focused on modern breeds, which are isolated subgroups with distinctive physical and, purportedly, behavioral characteristics. We interrogated breed stereotypes by surveying owners of 18,385 purebred and mixed-breed dogs and genotyping 2155 dogs. Most behavioral traits are heritable [heritability (h²) > 25%], and admixture patterns in mixed-breed dogs reveal breed propensities. Breed explains just 9% of behavioral variation in individuals. Genome-wide association analyses identify 11 loci that are significantly associated with behavior, and characteristic breed behaviors exhibit genetic complexity. Behavioral loci are not unusually differentiated in breeds, but breed propensities align, albeit weakly, with ancestral function. We propose that behaviors perceived as characteristic of modern breeds derive from thousands of years of polygenic adaptation that predates breed formation, with modern breeds distinguished primarily by aesthetic traits.

Poly(A) RNA sequencing reveals age-related differences in the prefrontal cortex of dogs

Dogs may possess a unique translational potential to investigate neural aging and dementia because they are prone to age-related cognitive decline, including an Alzheimer’s disease–like pathological condition. Yet very little is known about the molecular mechanisms underlying canine cognitive decline. The goal of the current study was to explore the transcriptomic differences between young and old dogs’ frontal cortex, which is a brain region often affected by various forms of age-related dementia in humans. RNA isolates from the frontal cortical brain area of 13 pet dogs, which represented 7 different breeds and crossbreds, were analyzed. The dogs were euthanized for medical reasons, and their bodies had been donated by their owners for scientific purposes. The poly(A) tail RNA subfraction of the total transcriptome was targeted in the sequencing analysis. Cluster analyses, differential gene expression analyses, and gene ontology analyses were carried out to assess which genes and genetic regulatory mechanisms were mostly affected by aging. Age was the most prominent factor in the clustering of the animals, indicating the presence of distinct gene expression patterns related to aging in a genetically variable population. A total of 3436 genes were found to be differentially expressed between the age groups, many of which were linked to neural function, immune system, and protein synthesis. These findings are in accordance with previous human brain aging RNA sequencing studies. Some genes were found to behave more similarly to humans than to rodents, further supporting the applicability of dogs in translational aging research.

Imprinting of canine IGF2 and H19

Genomic imprinting occurs in therian mammals and is a phenomenon whereby the two alleles of a gene are differentially expressed, based on the sex of the parent from whom the alleles were inherited. The allelic differences in expression are the consequence of different epigenetic modifications that are established in the sperm or oocyte during gametogenesis and transmitted at fertilization to offspring. A small minority of genes is regulated in this way but they have important biological functions, and aberrant regulation of imprinted genes contributes to disease aetiology in humans and other animals. The factors driving the evolution of imprinted genes are also of considerable interest, as these genes appear to forego the benefits of diploidy. To broaden the phylogenetic analysis of genomic imprinting, we began a study of imprinted genes in the domestic dog, Canis familiaris. In this report, we show that canine IGF2 and H19 are imprinted, with parent-of origin-dependent monoallelic expression patterns in neonatal umbilical cord. We identify a putative imprint control region associated with the genes, and provide evidence for differential methylation of this region in a somatic tissue (umbilical cord) and for its hypermethylation in the male germline. Canis familiaris is fast becoming a highly informative system for elucidating disease processes and evolution, and the study of imprinted genes in this species may help in understanding how these genes contribute to the generation of morphological and behavioral diversity.

Multi-omics approach identifies germline regulatory variants associated with hematopoietic malignancies in retriever dog breeds

Histiocytic sarcoma is an aggressive hematopoietic malignancy of mature tissue histiocytes with a poorly understood etiology in humans. A histologically and clinically similar counterpart affects flat-coated retrievers (FCRs) at unusually high frequency, with 20% developing the lethal disease. The similar clinical presentation combined with the closed population structure of dogs, leading to high genetic homogeneity, makes dogs an excellent model for genetic studies of cancer susceptibility. To determine the genetic risk factors underlying histiocytic sarcoma in FCRs, we conducted multiple genome-wide association studies (GWASs), identifying two loci that confer significant risk on canine chromosomes (CFA) 5 (P_wald = 4.83x10^-9) and 19 (P_wald = 2.25x10^-7). We subsequently undertook a multi-omics approach that has been largely unexplored in the canine model to interrogate these regions, generating whole genome, transcriptome, and chromatin immunoprecipitation sequencing. These data highlight the PI3K pathway gene PIK3R6 on CFA5, and proximal candidate regulatory variants that are strongly associated with histiocytic sarcoma and predicted to impact transcription factor binding. The CFA5 association colocalizes with susceptibility loci for two hematopoietic malignancies, hemangiosarcoma and B-cell lymphoma, in the closely related golden retriever breed, revealing the risk contribution this single locus makes to multiple hematological cancers. By comparison, the CFA19 locus is unique to the FCR and harbors risk alleles associated with upregulation of TNFAIP6, which itself affects cell migration and metastasis. Together, these loci explain ~35% of disease risk, an exceptionally high value that demonstrates the advantages of domestic dogs for complex trait mapping and genetic studies of cancer susceptibility.

We have identified two regions of the canine genome that explain a striking 35% of risk for developing histiocytic sarcoma in FCRs. The disease is uniformly lethal, affects 20% of FCRs, and parallels a cancer of the same name in humans. Both regions harbor genes involved in cell migration and cancer-related pathways. The first includes variants in regulatory regions at the tumor suppressor PIK3R6 locus that are strongly associated with histiocytic sarcoma and likely confer risk for other hematopoietic cancers. FCRs with risk alleles at the second locus demonstrate increased expression of TNFAIP6, which correlates with poor prognosis in multiple human cancers. In identifying genomic differences between affected and unaffected dogs, we advance our understanding of both canine and human health biology and set the stage for the development of diagnostic and therapeutic strategies.

Man's best friend in life and death: scientific perspectives and challenges of dog brain banking

Biobanking refers to the systematic collection, storage, and distribution of pre- or post-mortem biological samples derived from volunteer donors. The demand for high-quality human specimens is clearly demonstrated by the number of newly emerging biobanking facilities and large international collaborative networks. Several animal species are relevant today in medical research; therefore, similar initiatives in comparative physiology could be fruitful. Dogs, in particular, are gaining increasing attention in translational research on complex phenomena, like aging, cancer, and neurodegenerative diseases. Therefore, biobanks gathering and storing dog biological materials together with related data could play a vital role in translational and veterinary research projects. To achieve these aims, a canine biobank should meet the same standards in sample quality and data management as human biobanks and should rely on well-designed collaborative networks between different professionals and dog owners. While efforts to create dog biobanks could face similar financial and technical challenges as their human counterparts, they can widen the spectrum of successful collaborative initiatives towards a better picture of dogs’ physiology, disease, evolution, and translational potential. In this review, we provide an overview about the current state of dog biobanking and introduce the “Canine Brain and Tissue Bank” (CBTB)—a new, large-scale collaborative endeavor in the field.

Variants That Differentiate Wolf and Dog Populations Are Enriched in Regulatory Elements

Research on the genetics of domestication most often focuses on the protein-coding exons. However, exons cover only a minor part (1–2%) of the canine genome, whereas functional mutations may be located also in regions beyond the exome, in regulatory regions. Therefore, a large proportion of phenotypical differences between dogs and wolves may remain genetically unexplained. In this study, we identified variants that have high allelic frequency differences (i.e., highly differentiated variants) between wolves and dogs across the canine genome and investigated the potential functionality. We found that the enrichment of highly differentiated variants was substantially higher in promoters than in exons and that such variants were enriched also in enhancers. Several enriched pathways were identified including oxytocin signaling, carbohydrate digestion and absorption, cancer risk, and facial and body features, many of which reflect phenotypes of potential importance during domestication, including phenotypes of the domestication syndrome. The results highlight the importance of regulatory mutations during dog domestication and motivate the functional annotation of the noncoding part of the canine genome.

A novel canine reference genome resolves genomic architecture and uncovers transcript complexity

We present GSD_1.0, a high-quality domestic dog reference genome with chromosome length scaffolds and contiguity increased 55-fold over CanFam3.1. Annotation with generated and existing long and short read RNA-seq, miRNA-seq and ATAC-seq, revealed that 32.1% of lifted over CanFam3.1 gaps harboured previously hidden functional elements, including promoters, genes and miRNAs in GSD_1.0. A catalogue of canine "dark" regions was made to facilitate mapping rescue. Alignment in these regions is difficult, but we demonstrate that they harbour trait-associated variation. Key genomic regions were completed, including the Dog Leucocyte Antigen (DLA), T Cell Receptor (TCR) and 366 COSMIC cancer genes. 10x linked-read sequencing of 27 dogs (19 breeds) uncovered 22.1 million SNPs, indels and larger structural variants. Subsequent intersection with protein coding genes showed that 1.4% of these could directly influence gene products, and so provide a source of normal or aberrant phenotypic modifications.

Improving human cancer therapy through the evaluation of pet dogs

Comparative oncology clinical trials play an important and growing role in cancer research and drug development efforts. These trials, typically conducted in companion (pet) dogs, allow assessment of novel anticancer agents and combination therapies in a veterinary clinical setting that supports serial biologic sample collections and exploration of dose, schedule and corresponding pharmacokinetic/pharmacodynamic relationships. Further, an intact immune system and natural co-evolution of tumour and microenvironment support exploration of novel immunotherapeutic strategies. Substantial improvements in our collective understanding of the molecular landscape of canine cancers have occurred in the past 10 years, facilitating translational research and supporting the inclusion of comparative studies in drug development. The value of the approach is demonstrated in various clinical trial settings, including single-agent or combination response rates, inhibition of metastatic progression and randomized comparison of multiple agents in a head-to-head fashion. Such comparative oncology studies have been purposefully included in the developmental plan for several US FDA-approved and up-and-coming anticancer drugs. Challenges for this field include keeping pace with technology and data dissemination/harmonization, improving annotation of the canine genome and immune system, and generation of canine-specific validated reagents to support integration of correlative biology within clinical trial efforts.

Chromatin accessibility in canine stromal cells and its implications for canine somatic cell reprogramming

Naturally occurring disease in pet dogs is an untapped and unique resource for stem cell-based regenerative medicine translational research, given the many similarities and complexity such disease shares with their human counterparts. Canine-specific regulators of somatic cell reprogramming and pluripotency maintenance are poorly understood. While retroviral delivery of the four Yamanaka factors successfully reprogrammed canine embryonic fibroblasts, adult stromal cells remained resistant to reprogramming in spite of effective viral transduction and transgene expression. We hypothesized that adult stromal cells fail to reprogram due to an epigenetic barrier. Here, we performed assay for transposase-accessible chromatin using sequencing (ATAC-seq) on canine stromal and pluripotent stem cells, analyzing 51 samples in total, and establishing the global landscape of chromatin accessibility before and after reprogramming to induced pluripotent stem cells (iPSC). We also studied adult stromal cells that do not yield iPSC colonies to identify potential reprogramming barriers. ATAC-seq analysis identified distinct cell type clustering patterns and chromatin remodeling during embryonic fibroblast reprogramming. Compared with embryonic fibroblasts, adult stromal cells had a chromatin accessibility landscape that reflects phenotypic differentiation and somatic cell-fate stability. We ultimately identified 76 candidate genes and several transcription factor binding motifs that may be impeding somatic cell reprogramming to iPSC, and could be targeted for inhibition or activation, in order to improve the process in canines. These results provide a vast resource for better understanding of pluripotency regulators in dogs and provide an unbiased rationale for novel canine-specific reprogramming approaches.

Epigenetic Mechanisms in Canine Cancer

A plethora of data has highlighted the role of epigenetics in the development of cancer. Initiation and progression of different cancer types are associated with a variety of changes of epigenetic mechanisms, including aberrant DNA methylation, histone modifications, and miRNA expression. At the same time, advances in the available epigenetic tools allow to investigate and reverse these epigenetic changes and form the basis for the development of anticancer drugs in human oncology. Although human and canine cancer shares several common features, only recently that studies emerged investigating the epigenetic landscape in canine cancer and applying epigenetic modulators to canine cancer. This review focuses on the existing studies involving epigenetic changes in different types of canine cancer and the use of small-molecule inhibitors in canine cancer cells.

Transcriptomes from German shepherd dogs reveal differences in immune activity between atopic dermatitis affected and control skin

Canine atopic dermatitis (CAD) is an inflammatory and pruritic allergic skin disease with both genetic and environmental risk factors described. We performed mRNA sequencing of non-lesional axillary skin biopsies from nine German shepherd dogs. Obtained RNA sequences were mapped to the dog genome (CanFam3.1) and a high-quality skin transcriptome was generated with 23,510 expressed gene transcripts. Differentially expressed genes (DEGs) were defined by comparing three controls to five treated CAD cases. Using a leave-one-out analysis, we identified seven DEGs: five known to encode proteins with functions related to an activated immune system (CD209, CLEC4G, LOC102156842 (lipopolysaccharide-binding protein-like), LOC480601 (regakine-1-like), LOC479668 (haptoglobin-like)), one (OBP) encoding an odorant-binding protein potentially connected to rhinitis, and the last (LOC607095) encoding a novel long non-coding RNA. Furthermore, high mRNA expression of inflammatory genes was found in axillary skin from an untreated mild CAD case compared with healthy skin. In conclusion, we define genes with different expression patterns in CAD case skin helping us understand post-treatment atopic skin. Further studies in larger sample sets are warranted to confirm and to transfer these results into clinical practice.