# Investigation of the premelanosome protein (PMEL or SILV) gene and identification of polymorphism excluding it as the determinant of the dilute locus in domestic rabbits (Oryctolagus cuniculus)

Luca Fontanesi, Emilio Scotti, Michela Colombo, Daniel Allain, Severine Deretz, Stefania Dall'Olio, Vincenzo Russo, Ahmad Oulmouden

► To cite this version:

Luca Fontanesi, Emilio Scotti, Michela Colombo, Daniel Allain, Severine Deretz, et al.. Investigation of the premelanosome protein (PMEL or SILV) gene and identification of polymorphism excluding it as the determinant of the dilute locus in domestic rabbits (Oryctolagus cuniculus). Archiv für Tierzucht, 2013, 56, pp.42 - 49. 10.7482/0003-9438-56-005. hal-01211974

HAL Id: hal-01211974

https://hal.science/hal-01211974

Submitted on 29 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Short communication

# Investigation of the premelanosome protein (PMEL or SILV) gene and identification of polymorphism excluding it as the determinant of the dilute locus in domestic rabbits (*Oryctolagus cuniculus*)

Luca Fontanesi<fnref n="1" />, Emilio Scotti<fnref n="1" />, Michela Colombo<fnref n="1" />, Daniel Allain², Sevérine Deretz³, Stefania Dall'Olio<fnref n="1" />, Vincenzo Russo<fnref n="1" /> and Ahmad Oulmouden⁴

## Abstract

After the rediscovery of the Mendel's laws, the domesticated European rabbit (*Oryctolagus cuniculus*) has been the objective of pioneering studies on coat colour genetics. However, despite the early role of this species in defining genetic mechanisms determining this phenotypic trait, only recently a few loci have been characterized at the molecular level analysing also in rabbits genes already shown to affect coat colour in mice. We herein investigated the rabbit premelanosome protein (PMEL) gene, also known as melanocyte protein Pmel 17 (PMEL17) or silver (SILV), as mutations in the homologous gene in mice and other species produce phenotypic effects similar to what is observed in the dilute coat colour in rabbit. The rabbit dilute locus is determined by a recessive coat colour mutation that dilutes the black to blue (grey) interacting with the basic colours influenced by the agouti and extension loci. To investigate this candidate gene, we isolated and sequenced cDNAs as well as portions of intronic and exonic regions of the PMEL gene in several rabbits with different coat colours and identified single nucleotide polymorphisms, including several missense mutations. One polymorphism, positioned in intron 7, was genotyped in a family in which there was segregation of the dilute coat colour. The results excluded PMEL as the causative gene for the dilute locus in rabbits, shortening the list of candidate genes that should be analysed to identify the mutation determining this phenotypic trait.

<footnote n="1">¹Department of Agricultural and Food Sciences, Division of Animal Sciences, University of Bologna, Bologna, Italy, ²UR631, Station d'Amélioration Génétique des Animaux, INRA, Castanet Tolosan, France, ³UR967, Génétique Expérimentale en Productions Animales, INRA, Surgères, France, ⁴INRA, UMR1061 Unité de Génétique Moléculaire Animale, Université de Limoges, Limoges Cedex, France</footnote>
**Keywords:** candidate gene, coat colour, premelanosome protein, rabbit, SNP

**Abbreviations:** cDNA: complementary DNA, DNA: deoxyribonucleic acid, KLD: Kringle-like domain, NTR: N-terminal region, PCR: polymerase chain reaction, PKD: polycystic kidney disease, PMEL: premelanosome protein gene, PMEL17: melanocyte protein Pmel 17, RNA: ribonucleic acid, SNP: single nucleotide polymorphism

## Introduction

After the rediscovery of the Mendel's laws, the domesticated European rabbit (*Orycolagus cuniculus*) has been the objective of pioneering studies on coat colour genetics based mainly on fancy breeds/lines (Castle 1905, Castle 1930, Robinson 1958).

Despite the early role of the rabbit in defining genetic mechanisms determining this phenotypic trait, only recently a few loci have been characterized at the molecular level analysing also in rabbits (Aigner et al. 2000, Fontanesi et al. 2006, 2010a, 2010b) genes already shown to affect coat colour in mice (Kwon et al. 1987, Bultman et al. 1992, Robbins et al. 1993). Among these loci, *agouti* and *extension* control the production and relative amount of eumelanin (black/brown pigments) and pheomelanin (red/yellow pigments). They encode the *agouti* signalling protein (*ASIP*) and melanocortin 1 receptor (*MC1R*) genes, respectively (Bultman et al. 1992, Robbins et al. 1993).

Several other coat colour loci remain to be characterized at the DNA level in the rabbit. Among them, the dilute locus is determined by a recessive coat colour mutation that dilutes the black to blue (grey) and the yellow to beige interacting with the basic colours influenced by the *agouti* and *extension* mutations (Castle 1930, Robinson 1958, Searle 1968). In other species, similar phenotypic effects are determined by mutations in a few genes including the premelanosome protein (*PMEL*) gene, also known as melanocyte protein Pmel 17 (*PMEL17*) or silver (*SILV*). PMEL plays a central role in the biogenesis of the early stages of the pigment organelle, the melanosomes (Theos et al. 2005). In stage II melanosomes, processed PMEL protein aggregates to form fibrils, to which presumably the eumelanin is attached (Berson et al. 2001). Mutations within the *PMEL* gene altering the function of the coded protein are associated with the diluted-like of the basic coat colour in mouse, cattle, horse, dog and chicken (Martinez-Esparza et al. 1999, Kerje et al. 2004, Oulmouden et al. 2005, Brunberg et al. 2006, Clark et al. 2006, Gutiérrez-Gil et al. 2007, Kühn & Weikard 2007).

To evaluate the potential role of the rabbit *PMEL* gene in determining the *dilute* locus, we isolated and sequenced *PMEL* cDNAs as well as fragments of intronic and exonic regions and identified single nucleotide polymorphisms (SNPs), one of which was used to evaluate linkage with the dilute mutation in domestic rabbits.

## Materials and methods

### Animals

Twelve rabbits with different coat colours were used for resequencing of parts of the *PMEL* gene. This group of animals was made by: one Belgian Hare rabbit; 4 Vienna Blue rabbits, expected to be homozygous for the mutated dilute allele; 7 Checkered Giant rabbits, one with blue spots (buck), expected to be homozygous for the mutated dilute allele; 6 with black spots, 2 of which (does) expected to be heterozygous for the mutated dilute allele (as
they were obtained by crossing a black spotted buck with 2 other blue spotted does that were not possible to sample) and other 4 putative homozygous animals for the wild type allele at this locus, as deduced by pedigree information. In addition, two half sib-families with 5 and 7 F1 were produced crossing the sampled blue Checkered Giant buck with the two heterozygous (but black) Checkered Giant does. The full-length *PMEL* cDNA (see below) was obtained from dorsal skin of a Rex rabbit exhibiting wild type phenotype (Fontanesi *et al.* 2010a). Two additional rabbits of different breeds (Angora and Chinchilla) were used for *PMEL* cDNA isolation from dorsal skin as described below.

### DNA/RNA extraction and cDNA synthesis

Genomic DNA was extracted from hair roots as previously described (Fontanesi *et al.* 2007) or from blood (collected after slaughtering of the animals in a commercial abattoir) using a standard protocol. Dorsal skin specimens were collected from three rabbits after slaughtering, immediately frozen in liquid nitrogen and then stored at -80°C prior to RNA isolation that was carried out as already reported (Girardot *et al.* 2005). Reverse Transcription was performed from 1 µg of DNase treated total RNA in a volume of 20 µL, using 0.5 µg oligo (dT) primer and SuperScript II RNase H. Reverse Transcriptase (Invitrogen, Carlsbad, CA, USA), according to the manufacturer's instructions. The reaction was incubated for 50 min at 42 °C and then for 15 min at 70 °C.

### Rapid Amplification of cDNA Ends (RACE) and full-length cDNA cloning

To isolate the rabbit full-length *PMEL* cDNA, we performed RACE experiments with 5 µg of total RNA extracted from dorsal skin of a Rex rabbit using the SMART RACE cDNA Amplification kit (Clontech-Takara Bio Europe, Saint-Germain-en-Laye, France), according to the manufacturer's instructions. Consensus primers (Table 1) were designed based on coding sequence homology between human and bovine *PMEL* gene (GenBank acc. no. NM_006928 and Oulmouden *et al.* 2005). Nested PCR with consensus and adaptor primers was used for 5'-untranslated regions (UTR) and 3'-UTR amplification as previously described (Fontanesi *et al.* 2010a, Table 1). PCR products were cloned into the pCR2.1 vector (Invitrogen, Carlsbad, CA, USA) and sequenced as described below. The 5'- and 3'-end sequences were determined and used to design specific primers (sense and anti-sense) to amplify full-length cDNAs from the Rex rabbit and two additional animals, one Angora and one Chinchilla (Table 1).

### Identification and analysis of polymorphisms in rabbit families

Four PCR primer pairs were designed on the rabbit *PMEL* gene sequence (Ensembl no. ENSOCUG00000000999) to amplify the coding regions of 4 exons and parts of downstream or upstream intronic sequences (Table 1 and Figure 1). Genomic DNA used for PCR amplification and sequencing was from the 12 rabbits already described above. PCR was carried out using a PTC-100 thermal cycler (MJ Research, Watertown, MA, USA) in a 20 µL reaction volume containing ~50 ng genomic DNA, 1 U DNA EuroTaq DNA polymerase (EuroClone Ltd., Paington, Devon, UK), 1X Euro Taq PCR buffer, 2.5 mM dNTPs, 10 pmol of each primer and 2.0 mM of MgCl₂. PCR profile was as follows: 5 min at 95 °C; 35 amplification cycles of 30 sec at 95 °C, 30 s at 56-62 °C, 30 s at 72 °C; 5 min at 72 °C (Table 1). The amplified fragments
<table><caption>Table 1
PCR primer pairs and PCR conditions</caption><thead><tr><th>Primer name and amplified region</th><th>Primer forward and reverse sequences (5'-3')</th><th>PCR<sup>2</sup></th><th>Use</th></tr></thead><tbody><tr><td>5'-UTR first primer</td><td>CAATGGGAGCCAGGTGTGGGGAGGAC</td><td>-</td><td>5'-RACE</td></tr><tr><td>5'-UTR second primer</td><td>AATGATGGGCCTACACTGATTGGKGC</td><td>-</td><td>5'-RACE</td></tr><tr><td>5'-UTR first primer</td><td>YMCCCAGCCCAGCCTGCCAGCTGGTT</td><td>-</td><td>3'-RACE</td></tr><tr><td>5'-UTR second primer</td><td>GTCAGCACCCAGCTTRTCATGCCTGG</td><td>-</td><td>3'-RACE</td></tr><tr><td>Specific cDNA</td><td>GGCCTTTGGTTGCTGGAGAGAAGAACCTTCCTTCAGAGGCAGCCATGACCCG</td><td>2091/61/1.5</td><td>Full length CDS amplification of PMEL cDNA</td></tr><tr><td>Exon 1 region<sup>1</sup></td><td>TGTAGCAGGGTGGGGAAGCCAAACCAAATACATGGGAGA</td><td>399/59/2.0</td><td>Re-sequencing genomic region</td></tr><tr><td>Exon 2 region<sup>1</sup></td><td>CCATGGGAGACCAGGAGAAAGACTCAAAGCGGCACTCAG</td><td>396/55/1.0</td><td>Re-sequencing genomic region</td></tr><tr><td>Exon 5 region<sup>1</sup></td><td>GGCCCCAGACTTCCGAATGCTTCTCCTAGCCCTGTGAC</td><td>354/58/2.0</td><td>Re-sequencing genomic region</td></tr><tr><td>Exon 7 region<sup>1</sup></td><td>GCTCTGGGTCTACCTGGAAA GAGGTGCAAGGTTAAAGACTGG</td><td>446/59/2.0</td><td>Re-sequencing genomic region, PCR-RFLP with S<sub>sil</sub></td></tr></tbody></table>

<figure><img src="image_2.png" /><figcaption>Figure 1</figcaption></figure>

Structure of the rabbit PMEL gene (Ensembl no. ENSOCUG00000000999), re-sequenced regions from genomic DNA (indicated below the schematic representation of the gene), and position of the two intronic SNPs. Size of the re-sequenced regions does not include primers.

were sequenced after treatment with 1 µL of ExoSAP-IT® (USB Corporation, Cleveland, Ohio, USA) for 15 min at 37 °C. Cycle sequencing of the treated PCR products was produced using the same PCR primers and the Big Dye v3.1 kit (Applied Biosystems, Foster City, CA, USA). Sequencing reactions, after a few purification steps using EDTA, Ethanol 100% and Ethanol 70%, were loaded on an ABI3100 Avant capillary sequencer (Applied Biosystems, Foster City, CA, USA). Sequences were edited and aligned with the help of the CodonCode Aligner software (CodonCode Corporation, Dedham, MA, USA).

PCR-RFLP (Table 1) was applied to genotype the g.39510357G>A SNP (nomenclature of intronic SNPs is based on system coordinates of rabbit chromosome 4 in the oryCun2.0

<footnote n="1">¹Including intronic regions as described in *Results and discussion*. ²Information is related to the different primer pairs: amplified product size (bp)/annealing temperature (°C)/[MgCl₂] (mM). For the primers used in the RACE experiments, only the specific PMEL primers are reported (PCR conditions are the same as described in Fontanesi *et al.* [2010a]).</footnote>
genome version: http://www.ensembl.org/Oryctolagus_cuniculus/Info/) in the rabbit F1 families. Briefly, the amplified fragment of 446 bp obtained using the same exon 7 primers used for re-sequencing was digested with SsI. Digestion reaction was carried out overnight at 37°C in a 25 µL reaction volume including 5 µL of PCR product, 1× restriction enzyme buffer and 2 U of SsI (MBI Fermentas, Vilnius, Lithuania). The resulting DNA fragments were electrophoresed on 10% 29:1 polyacrylamide:bis acrylamide gel and visualized with 1× GelRed Nucleic Acid Gel Stain (Biotium Inc., Hayward, CA, USA). Restriction patterns were as follow: two fragments (353+93 bp) when allele A was present; one fragment of 446 bp if allele G was present. The PCR product obtained from 3 F1 rabbits was also sequenced as described before in order to confirm the PCR-RFLP genotyping results.

# Results and discussion

## Rabbit PMEL cDNA sequences

The rabbit PMEL cDNA sequence obtained from the Rex rabbit included 29 bp of the 5'-UTR, 1974 bp of the coding sequence (CDS) and 131 of the 3'-UTR (with the polyA tail), for a total of 2135 bp (GenBank/EMBL acc. no.: FR849710). The coding sequence encoded for a predicted 657 amino acid sequence, that was 79%, 82%, 80% and 78% identical to the human (661 amino acids, P40967), horse (662 amino acids, NP_001157361), pig (636 amino acids, BAE03309) and cattle (649 amino acids, Q06154) PMEL amino acid sequences, respectively. As expected, the predicted rabbit PMEL amino acid sequence has the characteristic signature for mammalian PMEL protein (Theos *et al.* 2005). A short N-terminal signal sequence of 24 residues is followed by a N-terminal region (NTR) of 200 amino acids that contains three putative N-terminal glycosylation sites and three cysteine residues. Downstream of the NTR is a polycystic kidney disease (PKD) domain followed by a conserved cysteine residue that separates a region rich in acidic residues, prolines, serines and threonines consisting primarily of imperfect direct repeated sequences (repeat domain, RPT). Following the RPT region is a conserved dibasic sequence that serves as a site for proteolytic cleavage that is critical for PMEL function. Downstream of the cleavage site is a conserved region known as the Kringle-like domain (KLD), characterized by the presence of highly conserved cysteines. Following a highly conserved cysteine and N-glycosylation site at the C-terminus of the KLD is a poorly conserved, short proline-rich region, the length of which varies among the two major splice forms in humans that separates the KLD from the single 24-residue membrane spanning domain.

Aligning the Rex rabbit cDNA sequence with that obtained from Angora and Chinchilla rabbits (GenBank/EMBL acc. no.: FR849711 and FR849709) and the rabbit cDNA sequence XM_002711060 (retrieved from GenBank and derived by automated computational analysis with gene prediction methods of sequencing data of the rabbit genome version oryCun2.0, that, in the meantime, has been completed at the level of ~7X by the Broad Institute) we identified, on the whole, 10 SNPs, 8 of which are missense polymorphisms (c.557T>G, p.V186G; c.584A>G, p.Q195R; c.603T>C; c.854C>T, p.P285L; c.878C>A, p.P293Q; c.1016A>T, p.Q339L; c.1256C>T, p.T419M; c.1346G>T, p.S449I; c.1501G>A, p.G501R; c.1770T>C), and 2 are synonymous SNPs (c.603T>C; c.1770T>C).
## *Re-sequencing of parts of the rabbit PMEL gene and identification of intronic polymorphisms*

To identify polymorphisms that could be useful to evaluate their possible co-segregation with the black-blue spotted phenotype in the F1 families, fragments of the rabbit *PMEL* gene including 4 out of 11 exons were re-sequenced in rabbits with different putative alleles at the dilute locus (Figure 1). These fragments encompassed exon 1 and upstream flanking and downstream intronic regions (360 bp), exon 2 with part of preceding intron 1 and downstream intron 2 regions (357 bp), exon 5 and parts of intron 4 and parts of intron 5 (316 bp), exon 7 and portions of upstream and downstream introns 6 and 7 (404 bp). Within a total of 1437 bp re-sequenced fragments, only 2 single nucleotide polymorphisms were identified in intronic regions: one was in intron 4 (g.39510357G>A); the second was in intron 7 (g.39508302G>A). Sequences have been submitted to GenBank/EMBL with acc. no. FR849712 and FR849713. These SNPs are not in any apparently functional or regulatory motif.

## *Exclusion of the rabbit PMEL gene as the determinant of the dilute locus*

The two identified SNPs were in linkage disequilibrium in the sequenced rabbits. Therefore we selected one of them (g.39510357G>A) for the analysis of the rabbit families in which there was segregation of blue and black coat colours. The genotyping of the animals was carried out by PCR-RFLP and sequencing of some of the analysed samples confirmed the genotyping results. Figure 2 shows pedigree information including the genotypes and the coat colour of the animals. The results clearly indicated that there is no co-segregation between the g.39510357G>A SNP and the coat colour of the F1 rabbits. Black and blue rabbits share the same genotypes. Therefore, the *PMEL* gene can be excluded as a possible determinant of the dilute coat colour locus in rabbit. It remains to be investigated if polymorphisms in this gene can have other effects on pigmentation, like the silvering in Champagne d'Argent breed, or can modify the intensity of the blue colour in rabbits showing this phenotype. Other identified cDNA polymorphisms will be analysed to this purpose.

We have already excluded another strong candidate gene [myosin VA (heavy chain 12, myoxin), *MYO5A*] for the dilute locus in rabbits (Fontanesi *et al.* 2012) that in mouse is responsible for the corresponding dilute locus (Mercer *et al.* 1991). Therefore, the exclusion of the *PMEL* gene further shorten the list of other possible candidate genes that have been already shown to determine hypopigmentation or dilution of coat colour in other species, like melanophilin (*MLPH*), or RAB27A, member RAS oncogene family (RAB27A) (Wilson *et al.* 2000, Matesic *et al.* 2001).

The characterization of this locus could add basic information to the biology of pigmentation in mammals as well as potential applications for fancy breeders, that mainly in countries like Germany, France, Italy and USA, contribute to maintain and develop rabbit breeds and lines.

## Acknowledgements

This work was funded by University of Bologna RFO and FaGenomicH projects and by INRA/
University of Limoges, and was carried out in the framework of the COST Action TD1101
»A Collaborative European Network on Rabbit Genome Biology – RGB-Net«.
<figure><img src="image_3.png" /><figcaption>Figure 2</figcaption></figure>

Genotype of the g.39510357G>A SNP of animals of the F1 families in which there is segregation of the black-blue spotted coat colours.

## References

* Aigner B, Besenfelder U, Müller M, Brem G (2000) Tyrosinase gene variants in different rabbit strains. Mamm Genome 11, 700-702

* Berson JF, Harper DC, Tenza D, Raposo G, Marks MS (2001) Pmel17 initiates premelanosome morphogenesis within multivesicular bodies. Mol Biol Cell 12, 3451-3464

* Brunberg E, Andersson L, Cothran G, Sandberg K, Mikko S, Lindgren G (2006) A missense mutation in PMEL17 is associated with the Silver coat color in the horse. BMC Genet 7, 46

* Bultman SJ, Michaud EJ, Woychik RP (1992) Molecular characterization of the mouse agouti locus. Cell 71, 1195-1204

* Castle WE (1905) Heredity of coat characters in guinea-pigs and rabbits. Carnegie Institution, Publication 23, Washington, USA

* Castle WE (1930) The Genetics of Domestic Rabbit. Harvard University Press, Cambridge, MA, USA

* Clark LA, Wahl JM, Rees CA, Murphy KE (2006) Retrotransposon insertion in SILV is responsible for merle patterning of the domestic dog. Proc Natl Acad Sci USA 103, 1376-1381

* Fontanesi L, Tazzoli M, Beretti F, Russo V (2006) Mutations in the melanocortin 1 receptor (MC1R) gene are associated with coat colours in the domestic rabbit (*Oryctolagus cuniculus*). Anim Genet 37, 489-493

* Fontanesi L, Tazzoli M, Russo V (2007) Non-invasive and simple methods for sampling rabbit DNA for PCR analysis of melanocortin 1 receptor (MC1R) gene mutations: a technical note. World Rabbit Sci 15, 121-126

* Fontanesi L, Forestier L, Allain D, Scotti E, Beretti F, Deretz-Picoulet S, Pecchioli E, Vernesi C, Robinson TJ, Malaney JL, Russo V, Oulmouden A (2010a) Characterization of the rabbit agouti signaling protein (ASIP) gene: transcripts and phylogenetic analyses and identification of the causative mutation of the nonagouti black coat colour. Genomics 95, 166-175

* Fontanesi L, Scotti E, Colombo M, Beretti F, Forestier L, Dall'Olio S, Deretz S, Russo V, Allain D, Oulmouden A (2010b) A composite six bp in-frame deletion in the melanocortin 1 receptor (MC1R) gene is associated with the Japanese brindling coat colour in rabbits (*Oryctolagus cuniculus*). BMC Genet 11, 59
Fontanesi L, Scotti E, Dall'Olio S, Oulmouden A, Russo V (2012) Identification and analysis of single nucleotide polymorphisms in the myosin VA (MYO5A) gene and its exclusion as the causative gene of the dilute coat colour locus in rabbit. World Rabbit Sci 20, 35-41

Girardot M, Martin J, Guibert S, Leveziel H, Julien R, Oulmouden A (2005) Widespread expression of the bovine *Agouti* gene results from at least three alternative promoters. Pigment Cell Res 18, 34-41

Gutiérrez-Gil B, Wiener P, Williams JL (2007) Genetic effects on coat colour in cattle: dilution of eumelanin and phaeomelanin pigments in an F2-Backcross Charolais × Holstein population. BMC Genet 8, 56

Kerje S, Sharma P, Gunnarsson U, Kim H, Bagchi S, Fredriksson R, Schütz K, Jensen P, von Heijne G, Okimoto R, Andersson L (2004) The *Dominant white*, *Dun* and *Smoky* color variants in chicken are associated with insertion/deletion polymorphisms in the *PMEL17* gene. Genetics 168, 1507-1518

Kühn Ch, Weikard R (2007) An investigation into the genetic background of coat colour dilution in a Charolais × German Holstein F2 resource population. Anim Genet 38, 109-113

Kwon BS, Haq AK, Pomerantz SH, Halaban R (1987) Isolation and sequence of a cDNA clone for human tyrosinase that maps at the mouse c-albino locus. Proc Natl Acad Sci USA 84, 7473-7477

Martinez-Esparza M, Jiménez-Cervantes C, Bennett DC, Lozano JA, Solano F, Garcia-Borrón JC (1999) The mouse *silver* locus encodes a single transcript truncated by the *silver* mutation. Mamm Genome 10, 1168-1171

Matesic LE, Yip R, Reuss AE, Swing DA, O'Sullivan TN, Fletcher CF, Copeland NG, Jenkins NA (2001) Mutations in *Mlph*, encoding a member of the Rab effector family, cause the melanosome transport defects observed in *leaden* mice. Proc Natl Acad Sci USA 98, 10238-10243

Mercer JA, Seperack PK, Strobel MC, Copeland NG, Jenkins NA (1991) Novel myosin heavy chain encoded by murine *dilute* (d) coat colour locus. Nature 349, 709-713

Oulmouden A, Julien R, Laforet JM, Leveziel H (2005) Use of silver gene for authentication of the racial origin of animal populations, and of the derivative products thereof. Patent Publication WO2005/019473

Robbins LS, Nadeau JH, Johnson KR, Kelly MA, Roselli-Rehfuss L, Baack E, Mountjoy KG, Cone RD (1993) Pigmentation phenotypes of variant *extension* locus alleles result from point mutations that alter MSH receptor function. Cell 72, 827-834

Robinson R (1958) Genetic studies of the rabbit. Bibl Genet 17, 229-558

Searle AG (1968) Comparative genetics of coat colour in mammals. Logos Press, London, UK

Theos AC, Truschel ST, Raposo G, Marks MS (2005) The *Silver* locus product Pmel17/gp100/Silv/ME20: controversial in name and in function. Pigment Cell Res 18, 322-336

Wilson SM, Yip R, Swing DA, O'Sullivan TN, Zhang Y, Novak EK, Swank RT, Russell LB, Copeland NG, Jenkins NA (2000) A mutation in *Rab27a* causes the vesicle transport defects observed in *ashen* mice. Proc Natl Acad Sci USA 97, 7933-7938