Next Article in Journal
Mitochondrial Genome Sequencing Reveals orf463a May Induce Male Sterility in NWB Cytoplasm of Radish
Next Article in Special Issue
Characterization of POU2F1 Gene and Its Potential Impact on the Expression of Genes Involved in Fur Color Formation in Rex Rabbit
Previous Article in Journal
Bone Marrow Clonogenic Myeloid Progenitors from NPM1-Mutated AML Patients Do Not Harbor the NPM1 Mutation: Implication for the Cell-Of-Origin of NPM1+ AML
 
 
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 1
10.3390/genes11010075
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
Correction published on 1 March 2021, see Genes 2021, 12(3), 357.
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Pigment Intensity in Dogs is Associated with a Copy Number Variant Upstream of KITLG

by
Kalie Weich
1,
Verena Affolter
2,
Daniel York
3,
Robert Rebhun
3,
Robert Grahn
4,
Angelica Kallenberg
4 and
Danika Bannasch
1,*
1
Department of Population Health and Reproduction, University of California-Davis, Davis, CA 95616, USA
2
Department of Pathology, Microbiology, and Immunology, University of California-Davis, Davis, CA 95616, USA
3
Department of Surgical and Radiological Sciences, University of California-Davis, Davis, CA 95616, USA
4
Veterinary Genetics Laboratory, University of California-Davis, Davis, CA 95616, USA
*
Author to whom correspondence should be addressed.
Genes 2020, 11(1), 75; https://doi.org/10.3390/genes11010075
Submission received: 10 December 2019 / Revised: 6 January 2020 / Accepted: 8 January 2020 / Published: 9 January 2020
(This article belongs to the Special Issue Coat Color Genetics)
Browse Figures
Versions Notes

Abstract

:
Dogs exhibit a wide variety of coat color types, and many genes have been identified that control pigment production, appearance, and distribution. Some breeds, such as the Nova Scotia Duck Tolling Retriever (NSDTR), exhibit variation in pheomelanin pigment intensity that is not explained by known genetic variants. A genome-wide association study comparing light red to dark red in the NSDTR identified a significantly associated region on canine chromosome 15 (CFA 15:23 Mb–38 Mb). Coverage analysis of whole genome sequence data from eight dogs identified a 6 kb copy number variant (CNV) 152 kb upstream of KITLG. Genotyping with digital droplet PCR (ddPCR) confirmed a significant association between an increased copy number with the dark-red coat color in NSDTR (p = 4.1 × 10−11). The copy number of the CNV was also significantly associated with coat color variation in eumelanin based Poodles (p = 2.1 × 10−11) and across other breeds. Moreover, the copy number correlated with pigment intensity along the hair shaft in both pheomelanin and eumelanin coats. KITLG plays an important role in melanogenesis, and variants upstream of KITLG have been associated with coat color variation in mice as well as hair color in humans consistent with its role in the domestic dog.

1. Introduction

The dog, Canis familaris, was domesticated from wolves (Canis lupus), although the exact time and place are still under scientific debate [1,2,3,4]. Evidence for morphological variation between dogs and wolves is found through the study and dating of archeological remains; however, the colors of these animals are unknown. Although there is some color variation within wolves, with the majority of their coat colors being muted to enhance camouflage [5,6]. A black coat color variant was introduced into wolves from dogs, potentially providing an advantage for hunting at night [5]. Color variants that lighten or darken the coat have been identified as early as 10,000 years ago in archeological dog samples, indicating that color variation existed early in the history of dogs [6]. In the mid-16th century, paintings depicting dogs show the richness of coat colors that had been developed within domesticated dogs by that time. One example of this is “The Hunters in the Snow” by Pieter Bruegel the Elder (1565), which shows rich red-, brown-, and black-pigmented hunting dogs. Modern dog breeds were developed in the last couple of hundred years and many breeds exhibit striking pigment intensity compared to the more muted colors of their ancestors, the wolf.
The basis of the main coat color types in domestic dogs have been determined [7,8,9,10]. Similar to other mammals, there are yellow (pheomelanin) and black pigments (eumelanin) whose production is controlled by the pigment-type switching genes encoding melanocortin 1 receptor (MC1R) and agouti signaling protein (ASIP) [7]. MC1R is a G protein-coupled receptor expressed by melanocytes. Signaling by its ligand, the α-melanocyte-stimulating hormone (MSH), promotes eumelanin synthesis [11]. Disruption of MC1R signaling occurs when the ASIP protein is expressed, which leads to pheomelanin synthesis. In wolves, a pulse of ASIP expression occurs during hair growth, leading to banded hair shafts on the dorsum [8]. Variants that alter either ASIP or MC1R can affect the pheomelanin or eumelanin production. Solid yellow or red dogs have a loss of function mutation in MC1R that results in the sole production of pheomelanin pigment [8]. Both a recessive loss-of-function ASIP variant and a dominant gain-of-function variant in the gene encoding β-defensin 103 (CBD103) can result in black coat colors in dogs [7,9].
Genetic variants have been identified that modify the basic coat colors. Breeds with exclusive pheomelanin production with a white or light cream-colored coat are enriched for a missense variant in the MFSD12 gene encoding major facilitator superfamily domain-containing 12 [12]. The dilution of pigment in some dogs can be the result of loss-of-function variants in the MLPH gene encoding melanophilin, resulting in abnormal clumping of pigment in the hair shaft and keratinocytes [13]. Brown coat color in dogs is an example of eumelanin lightening and is due to recessive alleles at TYRP1, which disrupt pigment production along the eumelanin pathway [14,15,16]. Pheomelanin appearance in dogs can vary from cream to deep red, with identical homozygous loss-of-function alleles at MC1R and MFSD12 ([12] and results presented here), indicating that there are additional loci that can affect pigment intensity in dogs.
In this study, the genetic basis for light and dark red pheomelanin phenotypes was evaluated using a genome-wide association in the Nova Scotia Duck Tolling Retriever (NSDTR). By analyzing sequence coverage from whole genome sequence data, a previously identified copy number variant (CNV) on chromosome 15 upstream of the KITLG coding sequence was found to be associated with the red color intensity in the NSDTR [17]. Phenotype genotype correlation across a spectrum of coat colors and breeds demonstrates that this locus governs both eumelanin and pheomelanin pigment intensity across dog breeds. Quantitative analysis of color along the hair shaft shows that low pigment intensity leading to lighter color is due to less pigment at the base of the hair. Conversely, high pigment intensity occurs when the hair is uniformly pigmented.

2. Materials and Methods

2.1. Sample Collection and DNA Extraction

Blood and saliva samples were collected from privately owned dogs through the University of California, Davis Veterinary Medical Teaching Hospital (VMTH) or submitted directly to the study from owners and veterinarians. Collection was performed under the supervision of the UC Davis Institutional Animal Care and Use Committee (protocol #20356). Genomic DNA from blood samples was extracted using the Qiagen Gentra Puregene Whole Blood Extraction Kit (QIAGEN, Valencia, CA, USA). Saliva samples were collected with Performagene Animal Collection swabs (DNA Genotek, Ottawa, ON, Canada) and extracted with the provided kit and protocol. Breed, date of birth, sex, and color were provided by owners.

2.2. Genome-Wide Association Study

Single-nucleotide variant (SNV) genotyping of 35 NSDTR was performed on the Illumina Canine HD 170 K BeadChip (Illumina, San Diego, CA, USA), which is mapped to the CanFam2.0 reference genome [18,19]. Dogs were categorized as light (n = 23) or dark (n = 13) based on visual inspection of their coat colors. The dogs were inspected in person or from photographs of the dogs in diffused lighting and graded by a single individual. Quality control and a Chi-square association analysis was performed with PLINK 1.9 [20,21]. All samples had a genotyping rate of at least 90%. In total, 63,746 SNVs were excluded due to a minor allele frequency of less than 5%, and 4238 SNVs were excluded due to a genotyping rate below 90%. Final analysis was performed with 105,678 SNVs, and the genome-wide significance was evaluated based on a Bonferroni threshold (p ≤ 4.7 × 10−7). Figures were created with R packages qqman [22], GenAbel [23], and cgmisc [24].

2.3. Whole-Genome Sequencing

Whole-genome sequencing of 100 canine genomes was performed as described by Brown et al. [25]. Thirteen dogs were NSDTR (6 light, 6 medium, and 1 dark red) and 1 was an Irish Setter (dark red). The sequencing depth was an average 8.7× coverage per sample, and reads were aligned to the canine reference genome (CanFam 3.1) [26]. Variants from CanFam3.1 chr15: 23–38 Mb were screened for homozygous segregation with a light or dark coat color. Variants were screened for conservation across human, mice, and rats with the 4-Way Multiz Alignment and Conservation track on the UCSC genome browser CanFam2 (https://genome.ucsc.edu/index.html) [27,28,29,30]. Variant prediction for biological consequences was performed with the Ensembl Variant Effect Predictor (VEP) [31,32]. Coverage analysis across the region surrounding the CNV (chr15:29,800,000–29,850,000) was performed with igvtools in sliding 2 kb windows for each sample [33]. The average coverage excluding the 12 kb of the CNV was then used to normalize each sample to the reference genome and graphed in Microsoft Excel.

2.4. Digital Droplet PCR Genotyping

Copy number quantification of samples was performed with a digital droplet PCR (ddPCR) assay on a BioRad QX200 Droplet Digital System (BioRad, Hercules, CA, USA). Primers for the KITLG CNV were designed from the CanFam3.1 reference genome (F: GAGTAGGTGTAATTTACCGGACA, R: AGCTATTTGCACAGGCTTTTT, probe: CCCATCCATCTTTACCTTCAGAAACA) with a 5′ 6-FAM (fluorescein) dye and quenched with 3′ Zen Iowa Black dye as designed by Integrated DNA Technologies (IDT, Coralville, IA, USA). A control gene assay for ETS proto-Oncogene 1 (ETS1) was designed by BioRad with HEX (hexachloro-fluorescein) dye to a region conserved across mammalian Taxa chr5:6,166,572–6,166,747 (BioRad Assay ID: 10042961). Prior to PCR, DNA was digested with New England Biosciences AvaII restriction enzyme (Product #R0153L, NEB, Ipswich, MA, USA). The copy number was calculated using the BioRad QX200 software standard protocol as a ratio of positive KITLG droplets to positive ETS1 droplets over all accepted droplets. Sample and primer concentrations were optimized to produce single templates per droplets. Statistical analysis of the results and generation of figures was performed with GraphPad Prism Version 8 (GraphPad, San Diego, CA, USA). A Mann-Whitney test was used to evaluate the differences between dark and light coat colors in NSDTRs and Poodles.

2.5. Quantitative Hair Pigment Analysis

In the Poodle, coat color information was provided by the owners (i.e., red, apricot, cream, white, black, silver, etc.). For breeds like the NSDTR, whose color variation is not recognized in the registered color of the dog, further quantification of the hair color from root to tip was performed by photographing the dogs’ coats under diffused light using a WhiBal G7 photography gray card (Michael Tapes Design, Melbourne, AUS). The color intensity was measured using Adobe photoshop in triplicate over 50 pixels at the root and tip in a single photograph. The mean intensity of root/tip was compared to the ddPCR results and evaluated for statistical significance with a linear regression test in GraphPad Prism.

3. Results

3.1. Genome-Wide Association Study in Light and Dark Red NSDTR

The NSDTR breed comes in a spectrum of coat colors ranging from a light golden red to a dark coppery red (Figure 1A and AKC NSDTR breed standard). None of the previously identified canine coat color variants explain this phenotypic variation. A GWAS was performed to identify a region in the genome that is associated with this variation in coat color using 23 light red NSDTR and 13 dark red NSDTR. All dogs were wild type for the previously identified variant at MFSD12 (orthologous to the human variant rs751,585,493 (ENST00000355415: Chr19:3,557,253 G > A (C > T in dogs); p.(Arg51Cys))) associated with lighter pheomelanin coat color [12]. One SNV was genome-wide significant with a pBonferroni of <0.05 on CanFam2 chr15: 32,383,555 (CanFam3.1 chr15: 29,371,013) (praw = 3.09 × 10−8, pBonferroni = 0.0033) (Figure 1B). Analysis of the linkage disequilibrium (LD) identified a very large critical region on CanFam 2 chromosome 15 from 23 to 38 Mb (CanFam3.1 chr15: 20–35 Mb) with r2 values >0.8 with respect to the most associated SNV (Figure 1C).

3.2. Whole Genome Sequence Variant and Coverage Analysis

In an initial screen for candidate causative variants, paired-end whole genome sequences of six light red NSDTR, one dark red NSDTR, and one dark red Irish Setter were investigated for variants in the critical region surrounding the best associated SNV. There were 81,560 SNV and 53,277 small insertion/deletion (indel) variants identified, and 112 of those variants segregated with coat color intensity (90 SNVs and 22 indels). None of these 112 variants were predicted to have protein coding effects, and none of the variants were in highly conserved regions of the genome. This region was also visually inspected on alignment files to screen for large indels or structural changes that were not identified by the variant calling software, and a literature search was performed to identify any known large structural variants in the region. Preliminary analysis of the aligned sequences of the abovementioned eight dogs plus an additional six medium red NSDTR showed that coverage over a previously identified 6 kb CNV at CanFam3.1 chr15:29,821,450–29,832,950 segregated with coat color, with dark red dogs showing an increase in the relative copy number compared to light red dogs (Figure 2A). Note that the reference genome shows two tandem copies of the CNV spanning a region of approximately 12 kb. The CNV is located in an intergenic region about 152 kb upstream from the closest gene, KITLG.

3.3. Validation of the KITLG CNV in Pheomelanin-Based Coat Colors

To validate the association of the high copy number at the KITLG CNV with the dark red coat color in the NSDTR, 28 dark red and 342 light red dogs confirmed to be homozygous for the loss of function alleles at MC1R were genotyped at the CNV with a digital droplet PCR (ddPCR) assay. Dark red NSDTR (median = 5 genomic copies) had a significantly higher copy number of the CNV compared to light red NSDTR (median = 2 genomic copies, p = 4.1 × 10−11) (Figure 2B). After verifying an association in the NSDTR, the CNV was then genotyped in additional breeds with a fixed dark red phenotype. Irish Setters (n = 48), a breed characterized by intense dark red pigment, had a median genotype of seven genomic copies. Another breed with dark red pigment, the Brittany (n = 4), had a median genotype of six genomic copies.
In order to confirm the segregation of the high copy number with dark red coats, comparisons were made of white and red-colored Poodles (Figure 3A). The comparison of white Poodles (n = 21) to red Poodles (n = 29) showed an increase in the copy number in red Poodles (median = 8 genomic copies) compared to white (median = 7 genomic copies, p = 0.003) however it was not statistically significant when white dogs with a homozygous genotype for pheomelanin dilution at MFSD12 (orthologous to human rs751,585,493: G > A) were excluded from the analysis (Figure 3B). Two other phenomenalistic breeds with non-functional MC1R that show coat color variation are the Golden Retriever (GR) and the yellow Labrador Retriever (LR). Their variation is not explained by the variant identified at MFSD12 [12]. Neither of these breeds showed a significant association with the copy number and light or dark coat color (GR, p = 0.9977, light n = 8, medium n = 19, dark n = 8; LR, p = 0.7377, light n = 8, dark n = 8). The breeds and number of dogs of each coat color are shown in Table S1.

3.4. Validation of the KITLG CNV in Eumelanin-Based Coat Colors

In the Poodle breed, pigment intensity differences occur in both pheomelanin and eumelanin-based coat colors (Figure 3A). Eumelanin-based colors range from silver or light grey to black, and these were evaluated for the CNV upstream of KITLG. Black Poodles (n = 42) had a significantly higher copy number (median = 6 genomic copies) compared to light grey (median = 2 genomic copies, n = 27) dogs (p = 2.1 × 10−11) (Figure 3B). Breeds fixed for light grey or black coat color were also genotyped, including light grey breeds like the Bearded Collie (median = 4 genomic copies, n = 5) and Old English Sheepdog (median = 2 genomic copies, n = 7), and black breeds like the Border Collie (BC, median = 3, n = 26), the Flat-Coated Retriever (FCR, median = 7, n = 20), and the Rottweiler (median = 8, n = 17). Additionally, grey wolves with unknown coat colors from North America (n = 2), Asia (n = 1), Europe (n = 1), and Africa (n = 1) were genotyped and all were found to have a very low copy number (approximately two genomic copies) (Figure 4).

3.5. Comparison of the Hair Shaft in High and Low Copy Number Dogs

Two breeds did not fit the association of the copy number with pigment intensity. The black and white BC did not have a high copy number as expected, and white Poodles had higher copy number than expected likely due to the presence of other mutations that affect their pheomelanin intensity. Eumelanin based Poodles did have a strong correlation of copy number and color. Poodles have mutations that alter their hair length, which functions by extending the anagen phase of the hair growth cycle. This led us to hypothesize that the pigment intensity might be due to changes along the length of the hair in dog breeds with more typical anagen phase growth. Changes in the pigment color along the hair shaft were evaluated quantitatively using high-resolution photographs and compared between high and low copy number dogs. Low copy number dogs appear to have less pigmentation at the root compared to high copy number dogs (Figure 5A). Linear regression analysis of 15 NSDTR with color intensity ratios from root to tip and estimated ddPCR genomic copy number identified a significant association of the high copy number with a low ratio (p = 0.0022) (Figure 5B). Additionally, a eumelanin breed with light roots (Border Collie) had significantly lower copy number compared to uniformly pigmented Flat Coated Retrievers FCR (p = 2.3 × 10−7) (Figure 5C).

4. Discussion

Segregation of pheomelanin intensity in the NSDTR enabled identification of a region on CFA 15 associated with pigment intensity in dogs. Variant analysis identified a 6 kb CNV upstream of KITLG that was highly associated with pheomelanin intensity within the NSDTR (p < 0.0001) and in other breeds. The copy number of this CNV is significantly associated eumelanin intensity in the Poodle and across breeds and to a lesser extent pheomelanin. All tested wolves had two copies of this element, indicating that chromosomes with higher copy numbers and high pigment intensity represent a derived state. Pigment intensity is due to the differing distribution of pigment along the hair shaft, with low copy number individuals having lighter hair at the root while more intensely pigmented individuals have no difference in pigment intensity from the root to hair tip. In the Poodle breed, the hair is continuously growing, such that the difference in colors associated with the different CNV alleles are particularly striking.
KITLG and its receptor, KIT, have previously been shown to function in hematopoiesis, melanogenesis, and gametogenesis [34]. Pleiotropic effects of genes responsible for coat color variation are common in mammals [35]. In this case, this ligand receptor pair is essential for stem cell populations to migrate and proliferate during development [36,37]. These stem cell populations include the primordial germ cells, hematopoietic progenitor cells, and melanoblasts. In addition to this early developmental function in melanoblast migration, KITLG is important for postnatal cutaneous melanogenesis [38]. Specifically, KITLG expression is essential in hair follicle epithelial cells for melanocyte terminal differentiation [39]. The essential role that KITLG plays in melanogenesis, both developmentally and in the hair shaft, is consistent with its identification as a pigmentation intensity locus in dogs.
Variants of KITLG and its receptor KIT cause pigmentation alternations in humans. Coding sequence variants in KITLG have been shown to cause familial progressive hyperpigmentation with or without hypopigmentation [40] as well as Waardenburg syndrome type 2 (with pigmentary abnormalities) [41]. Additionally, coding variants within the receptor KIT have been associated with classic autosomal dominant piebaldism characterized by congenital patches of skin with no melanocytes present [42,43,44,45]. However, regulatory variants have been shown to affect pigment intensity. Sulem et al. (2007) found that a variant upstream of KITLG was associated with hair color in a genome-wide association scan in Icelanders and Dutch [46]. This variant likely lies in linkage disequilibrium, with an extended haplotype upstream of KITLG showing a strong signature of selection in humans [46,47,48]. Guenther et al. (2014) found that this region drove expression exclusively in hair follicles using reporter constructs in mice [49]. KITLG has diverse important developmental functions; however, the alternation of just hair follicle expression may be a means to limit the effects to pigmentation phenotypes.
Regulatory variants of KITLG (or Kitl in the mouse) have been identified as causing pigment differences in other animal species. The original Kitl mouse mutants were called steel due to their uniformly diluted coat colors [50]. The steel panda allele is caused by an inversion that disrupts the upstream regulatory region of Kitl, which also leads to diluted pigmentation [49,51]. In both goats and mink, KITLG expression is lower in lighter-colored animals as compared to more darkly pigmented animals [52,53]. In this work, expression levels of KITLG were not evaluated due to challenges in obtaining skin samples from healthy pet dogs; however, based on the work done in other species, the CNV likely alters expression levels of KITLG in the hair follicle.
While many dog breeds tested had a pigment intensity correlated to the KITLG CNV, some breeds with pheomelanin variation did not. Notably, the Golden Retriever and Labrador Retriever had variability at the CNV, but it did not correlate to their coat colors, indicating that there are still additional pigment intensity loci or variants to be identified in dogs. In addition, white Poodles likely have multiple genes including the KITLG CNV and MFSD12 [12] involved with contributing to their light phenotype. It is not surprising that many loci affect pigmentation since dogs are under such strong artificial selection for the variety in coat colors.
The canine CNV upstream of KITLG was previously associated with susceptibility for digital squamous cell carcinoma (DSCC) in darkly colored dogs [54]. The analysis was performed in darkly pigmented cases of DSCC and control dogs, and a high copy number of this CNV was identified as the susceptibility locus; however, this effect was only found in eumelanistic cases with functional MC1R. KITLG variants have been implicated in melanoma, testicular cancer, colorectal cancer, and general cancer risk in humans [55,56,57,58]. The risk allele for DSCC in dogs had >four copies of the CNV and homozygous animals were more likely to have DSCC [54]. Black-pigmented dogs are highly overrepresented for DSCC, as 92% of samples were from black-coated dogs [59]. However, there are breeds like the Flat-Coated Retrievers that have a high copy number, black hair coat, and are not predisposed to this type of tumor [59]. It is possible that a high copy number at the KITLG CNV is another risk factor for breeds that are already more susceptible to this particular tumor type. Alternatively, some black dogs with high KITLG copy numbers may also have protective alleles.
Regulatory variants of human KITLG that affect hair color are located in a similar region to the CNV identified in dogs. Conservation tracks within the CNV are 86.1% identical to the human sequence. The human homologous region to the canine CNV is located 200 kb upstream of human KITLG. The strongest selective signature for KITLG in humans lies 218 kb upstream of the coding sequence [60] while the SNVs associated with light pigmentation are located 355 kb upstream of KITLG [47]. In dogs, there is also evidence that KITLG is under selective pressure as it was one of the top 20 regions of the genome to have a signature of selection within domestic dogs [61,62]. Interestingly, in canines, the selection by humans is for darker pigmentation over the wild-type wolf colors. Vivid and rich pigmentation within domestic dogs is one way that they may have been distinguished from wolves. Wolves have just a single copy of this 6 kb region, indicating that amplification of the region occurred within domestic dogs. It is interesting to speculate that this may have been one useful means of distinguishing partner proto-dogs from their wild ancestors.

5. Conclusions

The KITLG CNV is associated with coat pigment intensity in the domestic dog. An increase in the copy number is associated with a darker, more intense appearance to the coat and a more uniform pigment intensity across the hair shaft in many breeds. Wolves do not show copy number variation at this locus, and previous genome-wide screens have identified the KITLG CNV region as being under selection in domestic dogs. KITLG plays an important role in melanogenesis, both in development and in the adult hair follicle, and genetic variants upstream of KITLG in both mice and humans are associated with coat and hair color variation. This research suggests that the KITLG CNV is a new intensity locus in domestic dogs.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4425/11/1/75/s1, Table S1: Number of dogs genotyped for KITLG CNV copy number on ddPCR.

Author Contributions

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

Funding

This research was funded in part by the Maxine Adler Endowed Chair Funds, the Maxine Adler Graduate Fellowship, and the Center for Companion Animal Health, grant number CCAH 2019-8-F.

Acknowledgments

We would like to acknowledge and thank the owners who donated samples from their dogs for this study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Larson, G.; Karlsson, E.K.; Perri, A.; Webster, M.T.; Ho, S.Y.W.; Peters, J.; Stahl, P.W.; Piper, P.J.; Lingaas, F.; Fredholm, M.; et al. Rethinking dog domestication by integrating genetics, archeology, and biogeography. Proc. Natl. Acad. Sci. USA 2012, 109, 8878–8883. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Shannon, L.M.; Boyko, R.H.; Castelhano, M.; Corey, E.; Hayward, J.J.; McLean, C.; White, M.E.; Abi Said, M.; Anita, B.A.; Bondjengo, N.I.; et al. Genetic structure in village dogs reveals a Central Asian domestication origin. Proc. Natl. Acad. Sci. USA 2015, 112, 13639–13644. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Wang, G.-D.; Peng, M.-S.; Yang, H.-C.; Savolainen, P.; Zhang, Y.-P. Questioning the evidence for a Central Asian domestication origin of dogs. Proc. Natl. Acad. Sci. USA 2016, 113, E2554–E2555. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Botigue, L.R.; Song, S.; Scheu, A.; Gopalan, S.; Pendleton, A.L.; Oetjens, M.; Taravella, A.M.; Seregely, T.; Zeeb-Lanz, A.; Arbogast, R.M.; et al. Ancient European dog genomes reveal continuity since the Early Neolithic. Nat. Commun. 2017, 8, 16082. [Google Scholar] [CrossRef]
  5. Anderson, T.M.; vonHoldt, B.M.; Candille, S.I.; Musiani, M.; Greco, C.; Stahler, D.R.; Smith, D.W.; Padhukasahasram, B.; Randi, E.; Leonard, J.A.; et al. Molecular and evolutionary history of melanism in North American gray wolves. Science (New York N.Y.) 2009, 323, 1339–1343. [Google Scholar] [CrossRef] [Green Version]
  6. Ollivier, M.; Tresset, A.; Hitte, C.; Petit, C.; Hughes, S.; Gillet, B.; Duffraisse, M.; Pionnier-Capitan, M.; Lagoutte, L.; Arbogast, R.-M.; et al. Evidence of Coat Color Variation Sheds New Light on Ancient Canids. PLoS ONE 2013, 8, e75110. [Google Scholar] [CrossRef]
  7. Kerns, J.A.; Newton, J.; Berryere, T.G.; Rubin, E.M.; Cheng, J.F.; Schmutz, S.M.; Barsh, G.S. Characterization of the dog Agouti gene and a nonagoutimutation in German Shepherd Dogs. Mamm. Genome 2004, 15, 798–808. [Google Scholar] [CrossRef]
  8. Newton, J.M.; Wilkie, A.L.; He, L.; Jordan, S.A.; Metallinos, D.L.; Holmes, N.G.; Jackson, I.J.; Barsh, G.S. Melanocortin 1 receptor variation in the domestic dog. Mamm. Genome 2000, 11, 24–30. [Google Scholar] [CrossRef]
  9. Candille, S.I.; Kaelin, C.B.; Cattanach, B.M.; Yu, B.; Thompson, D.A.; Nix, M.A.; Kerns, J.A.; Schmutz, S.M.; Millhauser, G.L.; Barsh, G.S. A -defensin mutation causes black coat color in domestic dogs. Science (New York N.Y.) 2007, 318, 1418–1423. [Google Scholar] [CrossRef] [Green Version]
  10. Berryere, T.G.; Kerns, J.A.; Barsh, G.S.; Schmutz, S.M. Association of an Agouti allele with fawn or sable coat color in domestic dogs. Mamm. Genome 2005, 16, 262–272. [Google Scholar] [CrossRef]
  11. Robbins, L.S.; Nadeau, J.H.; Johnson, K.R.; Kelly, M.A.; Rosellirehfuss, L.; Baack, E.; Mountjoy, K.G.; Cone, R.D. Pigmentation phenotypes of variant extension locus alleles result from point mutations that alter msh receptor function. Cell 1993, 72, 827–834. [Google Scholar] [CrossRef]
  12. Hedan, B.; Cadieu, E.; Botherel, N.; Dufaure de Citres, C.; Letko, A.; Rimbault, M.; Drogemuller, C.; Jagannathan, V.; Derrien, T.; Schmutz, S.; et al. Identification of a Missense Variant in MFSD12 Involved in Dilution of Phaeomelanin Leading to White or Cream Coat Color in Dogs. Genes 2019, 10, 386. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Philipp, U.; Hamann, H.; Mecklenburg, L.; Nishino, S.; Mignot, E.; Gunzel-Apel, A.R.; Schmutz, S.M.; Leeb, T. Polymorphisms within the canine MLPH gene are associated with dilute coat color in dogs. BMC Genet. 2005, 6, 34. [Google Scholar] [CrossRef] [Green Version]
  14. Schmutz, S.M.; Berryere, T.G.; Goldfinch, A.D. TYRP1 and MC1R genotypes and their effects on coat color in dogs. Mamm. Genome 2002, 13, 380–387. [Google Scholar] [CrossRef] [PubMed]
  15. Cargill, E.J.; Famula, T.R.; Schnabel, R.D.; Strain, G.M.; Murphy, K.E. The color of a Dalmatian’s spots: Linkage evidence to support the TYRP1gene. BMC Vet. Res. 2005, 1, 1. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Hrckova Turnova, E.; Majchrakova, Z.; Bielikova, M.; Soltys, K.; Turna, J.; Dudas, A. A novel mutation in the TYRP1 gene associated with brown coat colour in the Australian Shepherd Dog Breed. Anim. Genet. 2017, 48, 626. [Google Scholar] [CrossRef]
  17. Nicholas, T.J.; Cheng, Z.; Ventura, M.; Mealey, K.; Eichler, E.E.; Akey, J.M. The genomic architecture of segmental duplications and associated copy number variants in dogs. Genome Res. 2009, 19, 491–499. [Google Scholar] [CrossRef] [Green Version]
  18. Lindblad-Toh, K.; Wade, C.M.; Mikkelsen, T.S.; Karlsson, E.K.; Jaffe, D.B.; Kamal, M.; Clamp, M.; Chang, J.L.; Kulbokas, E.J., III; Zody, M.C.; et al. Genome sequence, comparative analysis and haplotype structure of the domestic dog. Nature 2005, 438, 803–819. [Google Scholar] [CrossRef]
  19. Kim, K.S.; Lee, S.E.; Jeong, H.W.; Ha, J.H. The complete nucleotide sequence of the domestic dog (Canis familiaris) mitochondrial genome. Mol. Phylogenet. Evol. 1998, 10, 210–220. [Google Scholar] [CrossRef] [Green Version]
  20. Chang, C.C.; Chow, C.C.; Tellier, L.C.A.M.; Vattikuti, S.; Purcell, S.M.; Lee, J.J. Second-generation PLINK: Rising to the challenge of larger and richer datasets. GigaScience 2015, 4. [Google Scholar] [CrossRef]
  21. Steiß, V.; Letschert, T.; Schäfer, H.; Pahl, R. PERMORY-MPI: A program for high-speed parallel permutation testing in genome-wide association studies. Bioinformatics (Oxf. Engl.) 2012, 28, 1168–1169. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Turner, S.D. qqman: An R package for visualizing GWAS results using Q-Q and manhattan plots. bioRxiv 2014. [Google Scholar] [CrossRef]
  23. Karssen, L.C.; van Duijn, C.M.; Aulchenko, Y.S. The GenABEL Project for statistical genomics. F1000Res 2016, 5, 914. [Google Scholar] [CrossRef] [PubMed]
  24. Kierczak, M.; Jablonska, J.; Forsberg, S.K.; Bianchi, M.; Tengvall, K.; Pettersson, M.; Scholz, V.; Meadows, J.R.; Jern, P.; Carlborg, O.; et al. cgmisc: Enhanced genome-wide association analyses and visualization. Bioinformatics (Oxf. Engl. ) 2015, 31, 3830–3831. [Google Scholar] [CrossRef] [PubMed]
  25. Brown, E.A.; Dickinson, P.J.; Mansour, T.; Sturges, B.K.; Aguilar, M.; Young, A.E.; Korff, C.; Lind, J.; Ettinger, C.L.; Varon, S.; et al. FGF4 retrogene on CFA12 is responsible for chondrodystrophy and intervertebral disc disease in dogs. Proc. Natl. Acad. Sci. USA 2017, 114, 11476–11481. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Hoeppner, M.P.; Lundquist, A.; Pirun, M.; Meadows, J.R.S.; Zamani, N.; Johnson, J.; Sundström, G.; Cook, A.; FitzGerald, M.G.; Swofford, R.; et al. An Improved Canine Genome and a Comprehensive Catalogue of Coding Genes and Non-Coding Transcripts. PLoS ONE 2014, 9, e91172. [Google Scholar] [CrossRef] [PubMed]
  27. Kent, W.J.; Sugnet, C.W.; Furey, T.S.; Roskin, K.M.; Pringle, T.H.; Zahler, A.M.; Haussler, D. The human genome browser at UCSC. Genome Res. 2002, 12, 996–1006. [Google Scholar] [CrossRef] [Green Version]
  28. Kent, W.J.; Baertsch, R.; Hinrichs, A.; Miller, W.; Haussler, D. Evolution’s cauldron: Duplication, deletion, and rearrangement in the mouse and human genomes. Proc. Natl. Acad. Sci. USA 2003, 100, 11484–11489. [Google Scholar] [CrossRef] [Green Version]
  29. Chiaromonte, F.; Yap, V.B.; Miller, W. Scoring pairwise genomic sequence alignments. Pac. Symp. Biocomput. 2002, 115–126. [Google Scholar] [CrossRef]
  30. Blanchette, M.; Kent, W.J.; Riemer, C.; Elnitski, L.; Smit, A.F.; Roskin, K.M.; Baertsch, R.; Rosenbloom, K.; Clawson, H.; Green, E.D.; et al. Aligning multiple genomic sequences with the threaded blockset aligner. Genome Res. 2004, 14, 708–715. [Google Scholar] [CrossRef] [Green Version]
  31. McLaren, W.; Gil, L.; Hunt, S.E.; Riat, H.S.; Ritchie, G.R.S.; Thormann, A.; Flicek, P.; Cunningham, F. The Ensembl Variant Effect Predictor. Genome Biol. 2016, 17, 122. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. McLaren, W.; Pritchard, B.; Rios, D.; Chen, Y.; Flicek, P.; Cunningham, F. Deriving the consequences of genomic variants with the Ensembl API and SNP Effect Predictor. Bioinformatics (Oxf. Engl.) 2010, 26, 2069–2070. [Google Scholar] [CrossRef] [PubMed]
  33. Robinson, J.T.; Thorvaldsdóttir, H.; Winckler, W.; Guttman, M.; Lander, E.S.; Getz, G.; Mesirov, J.P. Integrative genomics viewer. Nat. Biotechnol. 2011, 29, 24. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Besmer, P. The kit ligand encoded at the murine Steel locus: A pleiotropic growth and differentiation factor. Curr. Opin. Cell Biol. 1991, 3, 939–946. [Google Scholar] [CrossRef]
  35. Reissmann, M.; Ludwig, A. Pleiotropic effects of coat colour-associated mutations in humans, mice and other mammals. Semin. Cell Dev. Biol. 2013, 24, 576–586. [Google Scholar] [CrossRef]
  36. Geissler, E.N.; Ryan, M.A.; Housman, D.E. The dominant-white spotting (W) locus of the mouse encodes the c-kit proto-oncogene. Cell 1988, 55, 185–192. [Google Scholar] [CrossRef]
  37. Brannan, C.I.; Lyman, S.D.; Williams, D.E.; Eisenman, J.; Anderson, D.M.; Cosman, D.; Bedell, M.A.; Jenkins, N.A.; Copeland, N.G. Steel-Dickie mutation encodes a c-kit ligand lacking transmembrane and cytoplasmic domains. Proc. Natl. Acad. Sci. USA 1991, 88, 4671–4674. [Google Scholar] [CrossRef] [Green Version]
  38. Grichnik, J.M.; Burch, J.A.; Burchette, J.; Shea, C.R. The SCF/KIT pathway plays a critical role in the control of normal human melanocyte homeostasis. J. Investig. Dermatol. 1998, 111, 233–238. [Google Scholar] [CrossRef] [Green Version]
  39. Liao, C.P.; Booker, R.C.; Morrison, S.J.; Le, L.Q. Identification of hair shaft progenitors that create a niche for hair pigmentation. Genes Dev. 2017, 31, 744–756. [Google Scholar] [CrossRef] [Green Version]
  40. Wang, Z.Q.; Si, L.; Tang, Q.; Lin, D.; Fu, Z.; Zhang, J.; Cui, B.; Zhu, Y.; Kong, X.; Deng, M.; et al. Gain-of-function mutation of KIT ligand on melanin synthesis causes familial progressive hyperpigmentation. Am. J. Hum. Genet. 2009, 84, 672–677. [Google Scholar] [CrossRef] [Green Version]
  41. Zazo Seco, C.; Serrao de Castro, L.; van Nierop, J.W.; Morin, M.; Jhangiani, S.; Verver, E.J.; Schraders, M.; Maiwald, N.; Wesdorp, M.; Venselaar, H.; et al. Allelic Mutations of KITLG, Encoding KIT Ligand, Cause Asymmetric and Unilateral Hearing Loss and Waardenburg Syndrome Type 2. Am. J. Hum. Genet. 2015, 97, 647–660. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Giebel, L.B.; Spritz, R.A. Mutation of the KIT (mast/stem cell growth factor receptor) protooncogene in human piebaldism. Proc. Natl. Acad. Sci. USA 1991, 88, 8696–8699. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Yang, Y.J.; Zhao, R.; He, X.Y.; Li, L.P.; Wang, K.W.; Zhao, L.; Tu, M.; Tang, J.S.; Xie, Z.G.; Zhu, Y.M. A novel splicing mutation of KIT results in piebaldism and auburn hair color in a Chinese family. BioMed Res. Int. 2013, 2013, 689756. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Richards, K.A.; Fukai, K.; Oiso, N.; Paller, A.S. A novel KIT mutation results in piebaldism with progressive depigmentation. J. Am. Acad. Dermatol. 2001, 44, 288–292. [Google Scholar] [CrossRef]
  45. Murakami, T.; Fukai, K.; Oiso, N.; Hosomi, N.; Kato, A.; Garganta, C.; Barnicoat, A.; Poppelaars, F.; Aquaron, R.; Paller, A.S.; et al. New KIT mutations in patients with piebaldism. J. Dermatol. Sci. 2004, 35, 29–33. [Google Scholar] [CrossRef]
  46. Sulem, P.; Gudbjartsson, D.F.; Stacey, S.N.; Helgason, A.; Rafnar, T.; Magnusson, K.P.; Manolescu, A.; Karason, A.; Palsson, A.; Thorleifsson, G.; et al. Genetic determinants of hair, eye and skin pigmentation in Europeans. Nat. Genet. 2007, 39, 1443–1452. [Google Scholar] [CrossRef]
  47. Miller, C.T.; Beleza, S.; Pollen, A.A.; Schluter, D.; Kittles, R.A.; Shriver, M.D.; Kingsley, D.M. cis-Regulatory changes in Kit ligand expression and parallel evolution of pigmentation in sticklebacks and humans. Cell 2007, 131, 1179–1189. [Google Scholar] [CrossRef] [Green Version]
  48. Yang, Z.; Shi, H.; Ma, P.; Zhao, S.; Kong, Q.; Bian, T.; Gong, C.; Zhao, Q.; Liu, Y.; Qi, X.; et al. Darwinian positive selection on the pleiotropic effects of KITLG explain skin pigmentation and winter temperature adaptation in Eurasians. Mol. Biol. Evol. 2018. [Google Scholar] [CrossRef]
  49. Guenther, C.A.; Tasic, B.; Luo, L.; Bedell, M.A.; Kingsley, D.M. A molecular basis for classic blond hair color in Europeans. Nat. Genet. 2014, 46, 748–752. [Google Scholar] [CrossRef] [Green Version]
  50. Sarvella, P.A.; Russell, L.B. Steel, a new dominant gene in the house mouse. J. Hered. 1956, 47, 123–128. [Google Scholar] [CrossRef] [Green Version]
  51. Bedell, M.A.; Brannan, C.I.; Evans, E.P.; Copeland, N.G.; Jenkins, N.A.; Donovan, P.J. DNA rearrangements located over 100 kb 5’ of the Steel (Sl)-coding region in Steel-panda and Steel-contrasted mice deregulate Sl expression and cause female sterility by disrupting ovarian follicle development. Genes Dev. 1995, 9, 455–470. [Google Scholar] [CrossRef] [Green Version]
  52. Song, X.; Xu, C.; Liu, Z.; Yue, Z.; Liu, L.; Yang, T.; Cong, B.; Yang, F. Comparative Transcriptome Analysis of Mink (Neovison vison) Skin Reveals the Key Genes Involved in the Melanogenesis of Black and White Coat Colour. Sci. Rep. 2017, 7, 12461. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Wu, S.; Li, J.; Ma, T.; Li, J.; Li, Y.; Jiang, H.; Zhang, Q. MiR-27a regulates WNT3A and KITLG expression in Cashmere goats with different coat colors. Anim. Biotechnol. 2019. [Google Scholar] [CrossRef]
  54. Karyadi, D.M.; Karlins, E.; Decker, B.; vonHoldt, B.M.; Carpintero-Ramirez, G.; Parker, H.G.; Wayne, R.K.; Ostrander, E.A. A copy number variant at the KITLG locus likely confers risk for canine squamous cell carcinoma of the digit. PLoS Genet. 2013, 9, e1003409. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Duffy, D.L.; Zhu, G.; Li, X.; Sanna, M.; Iles, M.M.; Jacobs, L.C.; Evans, D.M.; Yazar, S.; Beesley, J.; Law, M.H.; et al. Novel pleiotropic risk loci for melanoma and nevus density implicate multiple biological pathways. Nat. Commun. 2018, 9, 4774. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Landero-Huerta, D.A.; Vigueras-Villasenor, R.M.; Yokoyama-Rebollar, E.; Arechaga-Ocampo, E.; Rojas-Castaneda, J.C.; Jimenez-Trejo, F.; Chavez-Saldana, M. Epigenetic and risk factors of testicular germ cell tumors: A brief review. Front. Biosci. (Landmark. Ed.) 2017, 22, 1073–1098. [Google Scholar] [PubMed] [Green Version]
  57. Yang, S.; Li, W.S.; Dong, F.; Sun, H.M.; Wu, B.; Tan, J.; Zou, W.J.; Zhou, D.S. KITLG is a novel target of miR-34c that is associated with the inhibition of growth and invasion in colorectal cancer cells. J. Cell. Mol. Med. 2014, 18, 2092–2102. [Google Scholar] [CrossRef]
  58. Zeron-Medina, J.; Wang, X.; Repapi, E.; Campbell, M.R.; Su, D.; Castro-Giner, F.; Davies, B.; Peterse, E.F.; Sacilotto, N.; Walker, G.J.; et al. A polymorphic p53 response element in KIT ligand influences cancer risk and has undergone natural selection. Cell 2013, 155, 410–422. [Google Scholar] [CrossRef] [Green Version]
  59. Belluco, S.; Brisebard, E.; Watrelot, D.; Pillet, E.; Marchal, T.; Ponce, F. Digital squamous cell carcinoma in dogs: Epidemiological, histological, and immunohistochemical study. Vet. Pathol. 2013, 50, 1078–1082. [Google Scholar] [CrossRef] [Green Version]
  60. Williamson, S.H.; Hubisz, M.J.; Clark, A.G.; Payseur, B.A.; Bustamante, C.D.; Nielsen, R. Localizing recent adaptive evolution in the human genome. PLoS Genet. 2007, 3, e90. [Google Scholar] [CrossRef]
  61. Boyko, A.R.; Quignon, P.; Li, L.; Schoenebeck, J.J.; Degenhardt, J.D.; Lohmueller, K.E.; Zhao, K.; Brisbin, A.; Parker, H.G.; vonHoldt, B.M.; et al. A simple genetic architecture underlies morphological variation in dogs. PLoS Biol. 2010, 8, e1000451. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Vaysse, A.; Ratnakumar, A.; Derrien, T.; Axelsson, E.; Rosengren Pielberg, G.; Sigurdsson, S.; Fall, T.; Seppala, E.H.; Hansen, M.S.; Lawley, C.T.; et al. Identification of genomic regions associated with phenotypic variation between dog breeds using selection mapping. PLoS Genet. 2011, 7, e1002316. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Genes 11 00075 g001 550
Figure 1. Genome-wide association coat color in the NSDTR. (A) Representative NSDTR pictures from light red (left) and dark red (right) coat-colored dogs. (B) Quantile-quantile plot from NSDTR GWAS showing a genomic inflation, λ, of 1.03. Chromosome 15 SNVs are highlighted in red. (C) Manhattan plot. A single locus on chromosome 15 (in red) reached genome-wide significance after a Bonferroni correction (CanFam2 chr15: 32,383,555, praw = 3.09 × 10−8, pBonferroni = 0.0033). Solid line indicates genome-wide significance based on Bonferroni multiple testing correction. (D) SNVs from a 30 Mb region surrounding the most associated SNV (indicated by *). SNVs are color coded by the r2 value to show linkage disequilibrium in the region. Red dashed line indicates the Bonferroni threshold.
Figure 1. Genome-wide association coat color in the NSDTR. (A) Representative NSDTR pictures from light red (left) and dark red (right) coat-colored dogs. (B) Quantile-quantile plot from NSDTR GWAS showing a genomic inflation, λ, of 1.03. Chromosome 15 SNVs are highlighted in red. (C) Manhattan plot. A single locus on chromosome 15 (in red) reached genome-wide significance after a Bonferroni correction (CanFam2 chr15: 32,383,555, praw = 3.09 × 10−8, pBonferroni = 0.0033). Solid line indicates genome-wide significance based on Bonferroni multiple testing correction. (D) SNVs from a 30 Mb region surrounding the most associated SNV (indicated by *). SNVs are color coded by the r2 value to show linkage disequilibrium in the region. Red dashed line indicates the Bonferroni threshold.
Genes 11 00075 g001
Genes 11 00075 g002 550
Figure 2. Relative coverage across the region of the KITLG CNV. (A) Relative coverage from paired-end whole-genome sequencing is plotted on the y-axis, and CanFam3.1 chr15 base pair position is plotted on the x-axis. The region of the CNV is marked by the red bar on the x-axis. Each line represents a single individual dog, and the color of the line represents the coat color of that dog. The solid lines represent NSDTR and the dashed lines represent Brittanys and an Irish Setter. (B) NSDTR comparison between light and dark dogs using digital droplet PCR to estimate the copy number. Comparison using a Mann–Whitney test. p = 4.1 × 10−11. Dark red, n = 28. Light red, n = 34. **** = highly significant.
Figure 2. Relative coverage across the region of the KITLG CNV. (A) Relative coverage from paired-end whole-genome sequencing is plotted on the y-axis, and CanFam3.1 chr15 base pair position is plotted on the x-axis. The region of the CNV is marked by the red bar on the x-axis. Each line represents a single individual dog, and the color of the line represents the coat color of that dog. The solid lines represent NSDTR and the dashed lines represent Brittanys and an Irish Setter. (B) NSDTR comparison between light and dark dogs using digital droplet PCR to estimate the copy number. Comparison using a Mann–Whitney test. p = 4.1 × 10−11. Dark red, n = 28. Light red, n = 34. **** = highly significant.
Genes 11 00075 g002
Genes 11 00075 g003 550
Figure 3. Eumelanin and pheomelanin-based Poodle coat colors and CNV genotypes. (A) Representative Poodle hair samples from eumelanin low intensity silver (top left), eumelanin high intensity black (top right), pheomelanin low intensity cream (bottom left), and pheomelanin high intensity red (bottom right). (B) Genomic copy number measured with digital droplet PCR from high and low intensity Poodles with pheomelanin and eumelanin-based coat colors. Silver compared to black was significant (p = 1.5 × 10−8) while white compared to red was not. Sample numbers from left to right are as follows: n = 27, n = 425, n = 21, and n = 29. **** = highly significant.
Figure 3. Eumelanin and pheomelanin-based Poodle coat colors and CNV genotypes. (A) Representative Poodle hair samples from eumelanin low intensity silver (top left), eumelanin high intensity black (top right), pheomelanin low intensity cream (bottom left), and pheomelanin high intensity red (bottom right). (B) Genomic copy number measured with digital droplet PCR from high and low intensity Poodles with pheomelanin and eumelanin-based coat colors. Silver compared to black was significant (p = 1.5 × 10−8) while white compared to red was not. Sample numbers from left to right are as follows: n = 27, n = 425, n = 21, and n = 29. **** = highly significant.
Genes 11 00075 g003
Genes 11 00075 g004 550
Figure 4. Survey of the digital droplet-estimated genomic copy number in common dog breeds and wolves. Breeds known to segregate light and dark coat colors in eumelanin and/or pheomelanin-based colors were separated by color. Breeds are color coded by an approximation of the overall coat color, and sample numbers are provided in Table S1. * Breeds whose copy number estimates did not seem to match their coat colors.
Figure 4. Survey of the digital droplet-estimated genomic copy number in common dog breeds and wolves. Breeds known to segregate light and dark coat colors in eumelanin and/or pheomelanin-based colors were separated by color. Breeds are color coded by an approximation of the overall coat color, and sample numbers are provided in Table S1. * Breeds whose copy number estimates did not seem to match their coat colors.
Genes 11 00075 g004
Genes 11 00075 g005 550
Figure 5. Hair sample quantification from pheomelanin and eumelanin-based breeds. (A) Representative hair samples from (left to right) low, medium, and high copy number NSDTR, and low and high copy number eumelanin-based dogs (BC and FCR). (B) Linear regression analyses revealed a significant association between an increased copy number and lower color intensity ratio between root and tip of hair in NSDTR (r2 = 0.525, p = 0.0022, n = 15). (C) FCRs (n = 20) have a significantly higher copy number in comparison to BC (n = 26) in a ddPCR assay (p = 3.9 × 10−9). The graph color represents the approximate hair color at the root of the hair shaft. **** = highly significant.
Figure 5. Hair sample quantification from pheomelanin and eumelanin-based breeds. (A) Representative hair samples from (left to right) low, medium, and high copy number NSDTR, and low and high copy number eumelanin-based dogs (BC and FCR). (B) Linear regression analyses revealed a significant association between an increased copy number and lower color intensity ratio between root and tip of hair in NSDTR (r2 = 0.525, p = 0.0022, n = 15). (C) FCRs (n = 20) have a significantly higher copy number in comparison to BC (n = 26) in a ddPCR assay (p = 3.9 × 10−9). The graph color represents the approximate hair color at the root of the hair shaft. **** = highly significant.
Genes 11 00075 g005

Share and Cite

MDPI and ACS Style

Weich, K.; Affolter, V.; York, D.; Rebhun, R.; Grahn, R.; Kallenberg, A.; Bannasch, D. Pigment Intensity in Dogs is Associated with a Copy Number Variant Upstream of KITLG. Genes 2020, 11, 75. https://doi.org/10.3390/genes11010075

AMA Style

Weich K, Affolter V, York D, Rebhun R, Grahn R, Kallenberg A, Bannasch D. Pigment Intensity in Dogs is Associated with a Copy Number Variant Upstream of KITLG. Genes. 2020; 11(1):75. https://doi.org/10.3390/genes11010075

Chicago/Turabian Style

Weich, Kalie, Verena Affolter, Daniel York, Robert Rebhun, Robert Grahn, Angelica Kallenberg, and Danika Bannasch. 2020. "Pigment Intensity in Dogs is Associated with a Copy Number Variant Upstream of KITLG" Genes 11, no. 1: 75. https://doi.org/10.3390/genes11010075

APA Style

Weich, K., Affolter, V., York, D., Rebhun, R., Grahn, R., Kallenberg, A., & Bannasch, D. (2020). Pigment Intensity in Dogs is Associated with a Copy Number Variant Upstream of KITLG. Genes, 11(1), 75. https://doi.org/10.3390/genes11010075

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