A Third MLPH Variant Causing Coat Color Dilution in Dogs

Abstract

Altered melanosome transport in melanocytes, resulting from variants in the melanophilin (MLPH) gene, are associated with inherited forms of coat color dilution in many species. In dogs, the MLPH gene corresponds to the D locus and two variants, c.−22G > A (d¹) and c.705G > C (d²), leading to the dilution of coat color, as described. Here, we describe the independent investigations of dogs whose coat color dilution could not be explained by known variants, and who report a third MLPH variant, (c.667_668insC) (d³), which leads to a frameshift and premature stop codon (p.His223Profs*41). The d³ allele is found at low frequency in multiple dog breeds, as well as in wolves, wolf-dog hybrids, and indigenous dogs. Canids in which the d³ allele contributed to the grey (dilute) phenotype were d¹/d³ compound heterozygotes or d³ homozygotes, and all non-dilute related dogs had one or two D alleles, consistent with a recessive inheritance. Similar to other loci responsible for coat colors in dogs, this, alongside likely additional allelic heterogeneity at the D locus, or other loci, must be considered when performing and interpreting genetic testing.

1. Introduction Natural and artificial selection in dogs has led to a myriad of phenotypes with respect to coat coloration and its patterns. Two pigments, phaeomelanin (red/yellow) and eumelanin (black/brown), in which physical changes to the pigments themselves, as well as the timing and location of their deposition, contribute to a variety of color phenotypes. One such phenotype, dilution, is found in multiple mammalian and avian species. In dogs, the D (dilution) locus causes lightening of the coat color—i.e., from black to grey, brown to silver, or red to cream—as well as pigmentation lightening of the nose, paw pads, and eye color.

Melanosome transport, along cortical actin filaments in melanocytes, is regulated by three major proteins: Rab27A (a member of the small GTPase Rab family), the motor protein myosin Va (MYO5A), and melanophilin (MLPH) [1]. These three proteins form a functional ternary complex, in which an alteration in any of them can result in defects in melanosome transport and a dilute phenotype [1]. In 2007, Drögemüller et al. identified a transition at the splice boundary of the first untranslated exon of the canine MLPH gene, predicted to result in reduced splicing efficiency [2]. This non-coding variant, c.−22G > A (d¹), was associated with the dilution phenotype in seven breeds of dogs at the time [2]. A small percentage of dilute dogs did not possess a d¹ genotype predicted to result in the dilution phenotype, and Bauer et al. (2018) subsequently identified a second MLPH variant, c.705G > C (d²), located at the splice boundary of exon 7, and correlating to dilution in three additional breeds [3].

The effects of MLPH variants on melanosome transport and pigmentation are distinct from the recently identified variants of phaeomelanin-based pigment intensity, which also perturb the perceived depth of coloration, although resulting coat colors are varying shades of red, cream, and white [4]. Nevertheless, the previously identified d¹ and d² MLPH variants do not explain all of the dilute coat color phenotypes in canids, and a recent report suggests that additional coat color alleles are likely present at low frequency throughout the canine population [2,3,5]. Here we combined candidate gene sequencing with a query based on a large dataset of whole genome sequences, to identify a novel MLPH variant (d³), associated with dilute coat color in multiple dog breeds and wolf-dog hybrids.

2. Materials and Methods 2.1. Animals Cheek swab DNA samples were collected from 16 dilute dogs (3 Hungarian pumis, 5 Hungarian mudis, 5 wolf-dog hybrids, 1 Chihuahua, 1 Italian greyhound, and 1 Pekingese) that were previously genotyped for the d¹ and d² MLPH alleles, and incorrectly predicted to have a non-diluted coat color and 1 diluted d¹/d¹ Shih Tzu, along with 32 non-dilute dogs (11 Hungarian pumis, 15 Hungarian mudis, 2 wolf-dog hybrids, 1 Pekingese, and 3 Shih Tzu). Written consent was obtained from all dog owners. This study was approved under University of Minnesota IACUC protocol 1903 36865A. 2.2. MLPH Exon Sequencing

Genomic DNA was extracted from the cheek swab samples using Qiagen’s Puregene tissue kit or Qiagen’s DNEasy DNA extraction kit. The entire coding sequence (except exon 11), as well as flanking intronic boundaries, were Sanger sequenced in the Italian greyhound (Supplementary Table S1). For the pumis, mudis, Chihuahua and wolf-dog hybrids, an amplicon encompassing exon 7 of the MLPH gene, and extending into the flanking introns (including alleles d² and the 1 bp insertion reported here), was obtained (Supplementary Table S1). For the mudis and wolf-dog hybrids not previously genotyped for the d¹ allele, a 268 bp amplicon was obtained with d¹ PCR primers F and R (Supplementary Table S1). Genotyping was performed with direct Sanger sequencing of the amplicon at Eurofins Genomics (Louisville, KY, USA), and visualized with Sequencher software (Gene Codes Corporation, Ann Arbor, MI, USA), or sequenced by Laboklin and analyzed with SeqScanner Software, respectively.

2.3. Whole-Genome Database Query

The coding sequence and surrounding splice site regions of the canine MLPH gene, were inspected for potential loss-of-function mutations in a variant call file (VCF). This contained 804 dog (from >140 breeds) and 9 wolf genomes, available in an updated version of the Dog Biomedical Variant Database Consortium (DBVDC), utilizing Golden Helix GenomeBrowse, with public annotations from Ensembl Genes 89 and RefSeq Genes 104 [6,7].

3. Results 3.1. MLPH Exon Sequencing

Examples of dogs with dilute phenotypes included in our study are provided in Figure 1. Our initial investigations focused on the Hungarian pumi and mudi dogs, as well as an Italian greyhound, whose dilute phenotypes were not explained by the d¹ and d² variants genotyped by a commercial laboratory; i.e., these dilute dogs were not homozygous for the d¹ (c.−22G > A) or d² (c.705G > C) variants, and they were not compound heterozygotes (d¹/d²). The MLPH coding sequence, as well as flanking intronic boundaries, were sequenced in the Italian greyhound, and exons 1 and 7 were sequenced in the pumis and mudis. We confirmed the commercial laboratory results, where the dilute Italian greyhound and dilute pumis were heterozygous for d¹ and homozygous wild type for d². At the same time, we identified a cytosine insertion subsequent to a run of 8 cytosine bases, (NM_001103219.2: c.667_668insC or chr25: g.48150749_50insC (CanFam 3.1 assembly), hereafter referred to as d³. No other variants altering the protein coding sequence were found.

The canine MLPH gene contains 16 exons and 581 amino acid residues. As expected, considerable cross species homology exists across the full length of the protein; however, several regions of poor conservation are also present (Figure 2). The (c.667_668insC, p.His223Profs*41) insertion is predicted to recode 41 amino acids, before a stop codon is generated that truncates 54% of the wild type protein.

3.2. Variant Database Search

Inspection of the DBVDC VCF containing 804 dog and 9 wolf genomes, also revealed the previously unrecognized 1 bp insertion (d³) within exon 7, at an allele frequency of 0.0045. The d³ variant was present in 7 genomes in the heterozygous state (4 indigenous dogs, 1 Tibetan mastiff, 1 Yorkshire terrier, and 1 wolf), and no homozygous d³ individuals were found. Additionally, 1 heterozygous Shetland sheepdog was identified within the University of Minnesota’s private database. The frequencies of the d¹ and d² alleles in the DBVDC are 0.054 and 0.0025, respectively.

3.3. MLPH Genotyping

Sanger sequencing enabled us to resolve all three d³ variant genotypes, with the d³ variant homozygous in the dilute mudi (Figure 3). In addition to the d³ variant in the dilute wolf-dog hybrid, a previously recognized synonymous variant, c.669C > T (dbSNP rs852237859), was also observed. This SNP is relatively common in the DBVDC dogs at an allele frequency of 0.062. None of the dogs with the d³ allele in DBVDC have this SNP, while the wolf from DBVDC with the d³ allele does. The d² variant (also in exon 7) was not observed in our Sanger-sequenced study cohort.

3.4. Genotype-Phenotype Association

Table 1 and Supplementary Table S2 present the MLPH variant genotypes from the available phenotyped dogs, including 14 Hungarian pumis, 20 Hungarian mudis, 7 wolf-dog hybrids, 1 Chihuahua, 1 Italian greyhound, 2 Pekingese, and 4 Shih Tzu.

A perfect association between the d³ variant and the dilute coat color was found in three pumi families: all three grey pumis were compound heterozygotes (d¹/d³), whereas all non-dilute pumis had at least one D allele. The three families—P1–4, P5–11 and P12–14—were also investigated for allelic segregation within the phenotype. No sire was available for families P1–4 and P12–14. The available female P1 (mother of P2–P4) and P12 (mother of P13 and P14) carried one d¹ allele, which was transmitted to one of their respective offspring. Both pups were phenotypically grey and had d¹/d³ genotypes, and thus inherited the d³ allele form the sire. In the third family (P5–P11), both parents were heterozygotes: the mother for d³ and the father for d¹. Four descendants showed a black coat color being homozygous D/D or heterozygous D/d¹ and D/d³, and one pup was compound heterozygous d¹/d³ and showed a grey coat.

The D locus genotypes for the 20 phenotyped mudi are also shown in Table 1. The dilute mudi were either d¹/d³ compound heterozygotes or d³/d³ homozygotes, while all non-dilute dogs had at least one D allele. Transmission of the d³ allele was also demonstrated in a pedigree that produced both dilute (d³/d³) and non-dilute (D/d³) offspring.

Similarly, the grey wolf-dog hybrids were discordant at their d¹ locus genotype. These dogs were evaluated at the d² and d³ loci, and the genotypes d¹/d³ or d³ /d³ were confirmed, respectively. The mother and sibling of the d¹/d³ heterozygote were subsequently genotyped. The d¹ allele was inherited from the mother (AW1), and the d³ allele was likely inherited from the father, or represented a spontaneous mutation event. In three additional unrelated dilute wolf-dog hybrids (Figure 3), the previously recognized synonymous variant, c.669C > T (dbSNP rs852237859), was observed in addition to the d³ frameshift variant. The d² variant was not observed in any of the sequenced samples. Additionally, the DBVDC database contained 9 wolves: 7 were wt/wt at both d¹ and d³, 1 was heterozygous at d¹ only, and 1 was heterozygous at d³ only.

Additional examples of the association between d¹ and d³ variants, with dilute coat color in smaller cohorts of other breeds, included the Shih Tzu, Chihuahua, Pekingese, and Italian greyhound (Table 1). Examples of D/d³ heterozygotes (4 indigenous dogs, 1 Tibetan mastiff, 1 Yorkshire terrier, 1 wolf and 1 Shetland sheepdog) were identified in the aforementioned WGS variant databases. A single homozygous d³ Shih Tzu has been identified (ShT3) (Table 1, Table S2), but this individual was also homozygous for the d¹ allele. The non-dilute offspring (ShT5 and 6) from this individual, and a male wild type at all three known variants (ShT4), were heterozygous for both the d¹ and d³ alleles. However, from knowledge of the parents′ genotypes, it can be deduced that the pups received a chromosome containing both the d¹ and d³ variants from the dam, and a D chromosome from the sire, so the offspring still had one functional D allele. In this case, the d¹/d³ haplotype likely represents either a spontaneous mutation at the d³ locus on the d¹ background that is identical by state, or a recombination event.

Figure 4 provides additional images of the dogs genotyped in our study, and the impacts of additional coat color loci. Both mudis in (B), the mudi on the right in (C), and the mudis in (E) and (F) also have the merle mutation [8]. The merle mutation is inherited as autosomal dominant, and often results in uneven patches of pigment diluted to light grey, light brown, or white [8,9]. However, it is not phenotypically apparent in phaeomelanin. Additionally, the dog on the left in (B) has a TYRP1 recessive brown mutation [10,11,12]. Dogs homozygous recessive at both TYRP1 and MLPH, are diluted at brown instead of black, and have a resulting pigment color referred to in many breeds as “silver”, “isabella”, or “lilac”. In mudis, this color is referred to as “grey-brown”. The mudis in (E) and (F) have two copies of the recessive e allele of the MC1R gene, that affects phaeomelanin and results in varying shades of cream [4,10,13,14,15,16]. Additionally, the mudi in (E) is homozygous for d³, thus showing grey pigmentation on the nose and eye rims, and a lightened eye color.

4. Discussion

MLPH gene variants cause coat color dilution in dogs and other species [17,18,19,20,21]. However, not all dogs that are phenotypically dilute are explained by the two published MLPH dilution variants [2,3]. The d³ variant was identified at the UC Davis Veterinary Genetics Laboratory, in exon 7 of a phenotypically dilute Italian greyhound, with one copy of d¹ and no copies of d². At the same time, the d³ variant was detected at Laboklin in dilute pumis and wolf-dog hybrids, where testing for d¹ and d² confirmed the dogs to be compound heterozygous d¹/d³. Independently, a scan of the DBVDC for other MLPH coding sequence variants by the University of Minnesota, identified the same frameshift mutation in exon 7, after being presented with dilute mudis and wolf-dog hybrids, who were tested previously as D/D or D/d¹. All three independent identifications confirmed the mutation through Sanger sequencing.

Exon 7 harbors two of the three known canine MLPH mutations (Figure 5). The region is not highly conserved across mammals (Figure 2); however, it is in the vicinity of a predicted myosin Va interaction site, necessary for facilitating melanosome transport. The ablation of MLPH-myosin Va interaction, therefore, readily explains both the defective melanosome transport and the dilute grey phenotype. Taken together, the above arguments strongly suggest that the variant, c.667_668insC, contributes to the dilution of coat color in multiple breeds, and represents another loss-of-function allele similar to d¹ and d². We observed the association of the d³ allele with dilute coat color, both in its homozygous form and as a compound d¹/d³ heterozygote. That compound heterozygosity can result in the dilute coat color phenotype, has previously been shown with the d¹ and d² variants [3].

The d³ allele was found in a wide spectrum of breeds and other canids, minimally including the Hungarian pumi, Hungarian mudi, Chihuahua, Pekingese, Italian greyhound, Shih Tzu, Tibetan mastiff, Yorkshire terrier, and Shetland sheepdog, as well as indigenous dogs, wolves, and wolf-dog hybrids. The wolf-dog hybrids in this study that carry the d³ variant have recent ancestry behind them, including wolf, Alaskan Malamute, German shepherd dog, Siberian husky, and collie. The presence of the d³ variant in village dogs and wolves may indicate that this variant represents an ancient mutation, persisting at a low frequency across some breeds, while being lost in others. However, given that this variant is a single nucleotide insertion in an 8 bp poly C stretch, the possibility that d³ is identical by state, at least in some breeds, must be considered. The presence of the d³ variant on the d¹ background in the Shih Tzu, and the c.669C > T variant 2 bp from the wolf d³ allele, would suggest that this type of spontaneous mutation can occur. Alternatively, this could represent a recombination event between the d¹ and d³ haplotypes. Regardless, the identified Shih Tzu had two copies of the d¹/d³ haplotype, suggesting that this specific allele is present, at least in this breed. Despite a third dilution variant being discovered, not all phenotype/genotype discrepancies are explained. Overall, 99.74% of > 10,000 genotyped phenotypically dilute French bulldogs, as reported by their owners, are resolved by the known d¹ and d² variants, and the d³ variant has not yet been detected in this breed at all (private lab observation, data not shown). Since every discordant dog is heterozygous at d¹, and homozygous wild type at both d² and d³, it is possible that at least one additional MLPH variant exists, or a variant in a different gene altogether. This novel d³ variant adds to our understanding of the basis for coat color variation, and the interaction between variants in dogs and other species. Additionally, known coat color variants are useful to breeders in understanding the potential outcomes of matings. As this newly identified c.667_668insC d³ variant does not resolve all MLPH genotype/phenotype inconsistencies, additional variants are likely. This fact must be taken under consideration when performing genetic testing and counseling.

Figure 1. Examples of dogs with a dilute phenotype. A, a dilute Hungarian mudi puppy flanked by its black siblings; B, a dilute Hungarian pumi.

Figure 2. Amino acid sequence alignment of mammalian MLPH genes. The transcript IDs for the alignments are as follows: Human, NM_024101.7; cat, ENSFCAT00000039067.3; American mink, ENSNVIT00000032578.1; rabbit, ENSOCUT00000016497.4; and dog, NM_001103219.2. The position of the d3 dog variant is indicated with an *, and the recoded peptide preceding the premature stop codon follows in red.

Figure 3. Sanger sequencing traces from MLPH exon 7. The red box highlights the location of the d3 variant. A, a D/D mudi; B, a D/d3 mudi; C, a d3/d3 mudi; D, a d3/ d3 wolf-dog hybrid that also has a synonymous SNP (c.669C > T) 2 bp from the d3 mutation.

Figure 4. Hungarian mudi dogs and wolf-dog hybrids with a dilute coat color. A, a homozygous wild type mudi dog at d1, d2, and d3; B, two mudi dogs that are both compound heterozygotes for d1 and d3; C, two mudi dogs homozygous for d3; D, a d3 homozygous wolf-dog hybrid; E, an e/e, d3 homozygous mudi dog; F, an e/e, D/d3 mudi dog.

Figure 5. Locations of the d1, d2 and d3 variants within the MLPH gene. Exons are boxes, protein coding regions are blue and non-coding regions are tan. Protein interaction domains were adapted from a previous analysis in rabbit and American mink [18,21].

Breed	Color	D/D	D/d1	D/d3	d1/d3	d1/d1	d3/d3
Hungarian pumi	Black	5	4	2	0	0	0
Hungarian pumi	Grey	0	0	0	3	0	0
Hungarian mudi	Black *	5	3	6	0	0	0
Hungarian mudi	Grey *	0	0	0	2	0	2
Hungarian mudi	Cream *	0	0	1	0	0	1 #
Wolf-dog hybrid	Black	0	1	0	0	0	0
Wolf-dog hybrid	Grey	0	0	0	1	1	3
Wolf-dog hybrid	Other	1	0	0	0	0	0
Chihuahua	Grey	0	0	0	1	0	0
Italian greyhound	Grey	0	0	0	1	0	0
Pekingese	White	0	0	0	1 #	0	0
Pekingese	Brown	0	1	0	0	0	0
Shih Tzu	White	1	0	0	0	0	0
Shih Tzu	Grey	0	0	0	0	1 **
Shih Tzu	Non-dilute	0	2 **		0	0	0

Dogs were genotyped by Sanger sequencing as described in Methods. *, mudi with merle patterning are included; **, Shih Tzu chromosome containing both d1 and d3 variants; #, cream/white dogs with lightened nose, paw pad, and eye color. Genotypes capable of producing the dilute coat color are shaded in grey.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4425/11/6/639/s1, Table S1: Primers used for the PCR amplification and sequencing of MLPH exons, Table S2: Genotypes of all dogs in this study.

Author Contributions

Conceptualization, S.L.V.B., R.A.G. and A.K.; Data curation, K.M.M.; Formal analysis, S.L.V.B., K.M.M., R.A.G. and A.K.; Investigation, S.L.V.B., K.M.M., J.C.G. and J.R.C.; Methodology, K.M.M., J.M., J.R.M. and A.K.; Project administration, S.L.V.B.; Resources, K.M.M., R.A.G., J.R.M. and E.M.; Supervision, K.M.M. and R.A.G.; Validation, K.M.M., R.A.G., J.C.G. and A.K.; Visualization, S.L.V.B. and K.M.M.; Writing-original draft, S.L.V.B., K.M.M., R.A.G. and A.K.; Writing, review & editing, S.L.V.B., K.M.M., R.A.G., J.R.M., J.M., P.K. and A.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The authors have no external funding sources to acknowledge. We would first like to thank the owners who submitted samples, Terese Jonsson for providing the top left photo in Figure 4, and to all other owners for allowing us to use their images. Also, sincere gratitude to members of the DBVDC, the University of Bern, as well as Steve Friedenberg for the construction and maintenance of sequence and variant databases.

Conflicts of Interest

A.K., E.M. and P.K. are employees of Laboklin, a commercial laboratory offering genetic testing in dogs. R.A.G., J.C.G., J.M., and J.R.C. are employees at the Veterinary Genetics Laboratory, a commercial laboratory offering genetic testing in dogs.

References

1. Jo, C.; Park, H.I.; Jung, H.J.; Park, J.I.; Lee, J.E.; Myung, C.H.; Hwang, J.S. A novel function of Prohibitin on melanosome transport in melanocytes. Theranostics 2020, 10, 3880-3891.

2. Drögemüller, C.; Philipp, U.; Haase, B.; Günzel-Apel, A.R.; Leeb, T. A noncoding melanophilin gene (MLPH) SNP at the splice donor of exon 1 represents a candidate causal mutation for coat color dilution in dogs. J. Hered. 2007, 98, 468-473.

3. Bauer, A.; Kehl, A.; Jagannathan, V.; Leeb, T. A novel MLPH variant in dogs with coat colour dilution. Anim. Genet. 2018, 49, 94-97.

4. 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.

5. Dreger, D.L.; Hooser, B.N.; Hughes, A.M.; Ganesan, B.; Donner, J.; Anderson, H.; Holtvoigt, L.; Ekenstedt, K.J. True Colors: Commercially-acquired morphological genotypes reveal hidden allele variation among dog breeds, informing both trait ancestry and breed potential. PLoS ONE 2019, 14, e0223995.

6. Jagannathan, V.; Drögemüller, C.; Leeb, T. Dog Biomedical Variant Database Consortium (DBVDC). A comprehensive biomedical variant catalogue based on whole genome sequences of 582 dogs and eight wolves. Anim. Genet. 2019, 50, 695-704.

7. Golden Helix GenomeBrowse® visualization tool (Version 2.x). Available online: http://www.goldenhelix.com (accessed on 11 April 2020).

8. Clark, L.A.; Wahl, J.M.; Rees, C.A. Murphy KE. Retrotransposon insertion in SILV is responsible for merle patterning of the domestic dog. Proc. Natl. Acad. Sci. USA 2006, 103, 1376-1381.

9. Langevin, M.; Synkova, H.; Jancuskova, T.; Pekova, S. Merle phenotypes in dogs-SILV SINE insertions from Mc to Mh. PLoS ONE 2018, 13, e0198536.

10. 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.

11. Hrckova, T.E.; Majchrakovam, Z.; Bielikovam, M.; Soltysm, K.; Turnam, J.; Dudasm, A. A novel mutation in the TYRP1 gene associated with brown coat colour in the Australian Shepherd Dog Breed. Anim. Genet. 2017, 48, 626.

12. Wright, H.E.; Schofield, E.; Mellersh, C.S.; Burmeister, L.M. A novel TYRP1 variant is associated with liver and tan coat colour in Lancashire Heelers. Anim. Genet. 2019, 50, 783.

13. Everts, R.E.; Rothuizen, J.; van Oost, B.A. Identification of a premature stop codon in the melanocyte-stimulating hormone receptor gene (MC1R) in Labrador and Golden retrievers with yellow coat colour. Anim. Genet. 2000, 31, 194-199.

14. Schmutz, S.M.; Berryere, T.G.; Ellinwood, N.M.; Kerns, J.A.; Barsh, G.S. MC1R studies in dogs with melanistic mask or brindle patterns. J. Hered. 2003, 94, 69-73.

15. Dreger, D.L.; Schmutz, S.M. A new mutation in MC1R explains a coat color phenotype in 2 "old" breeds: Saluki and Afghan hound. J. Hered. 2010, 10, 644-649.

16. Dürig, N.; Letko, A.; Lepori, V.; Hadji Rasouliha, S.; Loechel, R.; Kehl, A.; Hytönen, M.K.; Lohi, H.; Mauri, N.; Dietrich, J.; et al. Two MC1R loss-of-function alleles in cream-coloured Australian Cattle Dogs and white Huskies. Anim. Genet. 2018, 49, 284-290.

17. Ishida, Y.; David, V.A.; Eizirik, E.; Schäffer, A.A.; Neelam, B.A.; Roelke, M.E.; Hannah, S.S.; O'Brien, S.J. Menotti-Raymond, M. A homozygous single-base deletion in MLPH causes the dilute coat color phenotype in the domestic cat. Genomics 2006, 88, 698-705.

18. Lehner, S.; Gähle, M.; Dierks, C.; Stelter, R.; Gerber, J.; Brehm, R.; Distl, O. Two-exon skipping within MLPH is associated with coat color dilution in rabbits. PLoS ONE 2013, 8, e84525.

19. Li, W.; Sartelet, A.; Tamma, N.; Coppieters, W.; Georges, M.; Charlier, C. Reverse genetic screen for loss-of-function mutations uncovers a frameshifting deletion in the melanophilin gene accountable for a distinctive coat color in Belgian Blue cattle. Anim. Genet. 2016, 47, 110-113.

20. Vaez, M.; Follett, S.A.; Bed'hom, B.; Gourichon, D.; Tixier-Boichard, M.; Burke, T. A single point-mutation within the melanophilin gene causes the lavender plumage colour dilution phenotype in the chicken. BMC Genet. 2008, 9.

21. Cirera, S.; Markakis, M.N.; Christensen, K.; Anistoroaei, R. New insights into the melanophilin (MLPH) gene controlling coat color phenotypes in American mink. Gene 2013, 527, 48-54.