- Research article
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 Geneticsvolume 11, Article number: 59 (2010)
In the domestic rabbit (Oryctolagus cuniculus), classical genetic studies have identified five alleles at the Extension locus: ED(dominant black), ES(steel, weaker version of ED), E (wild type, normal extension of black), eJ(Japanese brindling, mosaic distribution of black and yellow) and e (non-extension of black, yellow/red with white belly). Sequencing almost the complete coding sequence (CDS) of the rabbit MC1R gene, we recently identified two in-frame deletions associated with dominant black (c.280_285del6; alleles EDor ES) and recessive red (c.304_333del30; allele e) coat colours. It remained to characterize the eJallele whose phenotypic effect is similar to the Orange and Sex-linked yellow loci of cat and Syrian hamster.
We sequenced the whole CDS in 25 rabbits of different coat colours including 10 Japanese and 10 Rhinelander (tricolour) rabbits and identified another 6 bp-in frame deletion flanked by a G > A transition in 5' (c.[124G>A;125_130del6]) that was present in all animals with Japanese brindling coat colour and pattern. These mutations eliminate two amino acids in the first transmembrane domain and, in addition, cause an amino acid substitution at position 44 of the wild type sequence. Genotyping 371 rabbits of 31 breeds with different coat colour this allele (eJ) was present in homozygous state in Japanese, Rhinelander and Dutch tricolour rabbits only (except one albino rabbit). Rabbits with eJ/eJgenotype were non fixed at the non-agouti mutation we previously identified in the ASIP gene. Segregation in F1 and F2 families confirmed the order of dominance already determined by classical genetic experiments with a possible dose effect evident comparing eJ/eJand eJ/e animals. MC1R mRNA was expressed in black hair skin regions only.
The c.[124A;125_130del6] allele may be responsible for a MC1R variant determining eumelanin production in the black areas. However, the mechanism determining the presence of both red and black hairs in the same animal seems more complex. Expression analyses of the c.[124A;125_130del6] allele suggest that MC1R transcription may be regulated epigenetically in rabbits with the Japanese brindling phenotype. Further studies are needed to clarify this issue.
Coat colour in mammals is determined by the presence, distribution and biochemical activity of the melanocytes, which are specialized cells where eumelanins (black/brown pigments) and pheomelanins (yellow/red pigments) are synthesized. More than 300 loci have been shown to affect coat colour in mice regulating or altering melanocyte development and migration during embryogenesis, melanocyte morphology and functions, its components and its enzymatic machinery [1, 2]. However, Extension and Agouti are the main loci that affect the production and relative amount of these two melanin types in the melanocytes . These loci show epistatic interactions and usually wild type Extension is required for Agouti expression. Dominant alleles at the Extension locus produce black pigmentation, whereas recessive alleles extend the production of pheomelanins, causing yellow/red/pale pigmentation. On the contrary, dominant Agouti alleles determine pheomelanic phenotypes whereas recessive alleles cause black coat colour with a few exceptions.
The Extension locus encodes the melanocortin 1 receptor (MC1R) . This protein belongs to the seven transmembrane G protein coupled receptors that binds the α melanocyte-stimulating hormone (αMSH) inducing eumelanin synthesis. The Agouti locus encodes the agouti signaling protein (ASIP) that is a small paracrine signalling protein (131-135 amino acids in different mammals) . ASIP affects pigmentation blocking the αMSH-MC1R interaction which, in turn, causes a switch in pigment type from eumelanins to phaeomelanins [6, 7].
Mutations of the MC1R gene associated with different coat colours have been described in several species (mice , humans , guinea pigs , cattle [10–12], pigs [13, 14], horses , sheep , goats , dogs [18, 19], foxes [20, 21], bears , felids , pocket and beach mice [24, 25], squirrels , chickens , Japanese quails , bananaquits , guinea fowl  and reptiles ) in which gain of function mutations produce black/dark coat colour, whereas loss of function mutations usually cause yellow/red coat colour.
In the domestic rabbit (Oryctolagus cuniculus), classical genetic studies involving crossbreeding experiments among breeds with different coat colours have identified five alleles at the Extension locus: ED(dominant black), ES(steel, weaker version of ED), E (wild type, normal grey or normal extension of black), eJ(Japanese brindling, mosaic distribution of black and yellow) and e (non-extension of black, yellow/red with white belly) [3, 32–34]. The order of dominance is the following: ED> ES> E > eJ> e, with possible partial dominance of EDover ESand of ESover E[3, 32–34]. Sequencing almost the complete coding sequence (CDS) of the rabbit MC1R gene, we recently identified two in-frame deletions associated with dominant black (c.280_285del6; alleles EDand/or ES) and recessive red (c.304_333del30; allele e) coat colours, respectively . The dominant black deletion eliminates two amino acids in the second transmembrane domain (Figure 1). We could not clarify if this deletion corresponds to the EDor ESallele as sequencing and genotyping in a large number of rabbits across 30 breeds/strains did not identify any other mutation associated with dominant black (and/or steel). Therefore it could be hypothesized the existence of only one dominant black allele and the presence of modifier genes might be the reason of different intensities of black (steel) of the two putative alleles identified by classical genetic studies. The recessive red deletion eliminates 10 amino acids of the first extracellular loop (Figure 1) and was in homozygous state in all yellow/red rabbit breeds and strains (Burgundy Fawn, Gold Saxony, New Zealand Red, and Thuringian breeds and yellow/red animals of other breeds not fixed for any coat colour) . It remained to characterize the eJ(Japanese brindling) allele. The Japanese brindling coat colour variety probably appeared in France within the second half of the 19th century and subsequently was introduced in England and other countries [36, 37]. The Japanese brindling pattern can be defined as a yellow coat mottled with black, therefore this coat colour seems determined by the presence of two different types of melanocytes in different skin areas, one producing eumelanin and another one producing pheomelanin [3, 33, 36, 37]. The eJallele should be fixed in the Japanese and Rhinelander breeds [33, 34]. The former breed (named after its coat colour) has the classic mottled phenotype (Figure 2A) that is similar to the tortoiseshell pattern observed in female cats and female Syrian hamsters. The coat colour in these two species is caused by the Orange and Sex-linked yellow loci, respectively, that are two non-homologous X chromosome loci [38, 39]. The coat colour phenotype determined by the rabbit eJallele is also similar to that determined in guinea-pig by the ep(tortoiseshell) allele at the Extension locus described by classical genetic studies [40, 41] but not characterized at the molecular level, yet . The Rhinelander rabbit breed has a tricolour pattern with black, yellow and white areas (Figure 2B). The white regions are caused by the absence of mature melanocytes due to a dominant English spotting locus allele [3, 32, 33, 35, 42]. The coloured patches are much larger than those of the tortoiseshell-like pattern probably because a reduced number of melanoblasts in the skin allows for a spatial expansion of the melanoblast-derived clones. This tricolour pattern is similar to the calico coat colour in cat determined by the epistatic interaction between the Spotting and Orange loci  and in guinea-pigs to the effect of the white-spotting factor s in ep/ epanimals . An Extension allele with similar (but not overlapping) phenotypic effect to the rabbit eJallele has been described in pig . In this species the Epallele, that is caused by a frameshift mutation (a 2 bp insertion) in the MC1R coding sequence, determines red coat colour in Hereford, Linderöd and Tamworth pigs but frequent somatic reversions re-establish the correct reading frame and produce black coat colour in skin areas determining a spotted phenotype in a red background. As a matter of fact, epistatic effects of the Dominant white/KIT locus on this Extension allele might determine the white and black spotted phenotype or the complete white colour of some pig breeds [14, 43].
Here we investigated the MC1R gene in Japanese and Rhinelander rabbits and identified a novel in-frame deletion associated with the eJcoat colour. In addition we analysed in these rabbits the MC1R gene expression in skin specimens of different coloured regions.
Rabbit MC1R gene sequences
In our previous study we did not completely sequence the CDS of the rabbit MC1R gene and for some breeds the Extension allele status was obtained by genotyping the two in-frame deletions identified and not by sequencing . Here we amplified and sequenced 1348 bp (considering a wild type sequence) of the MC1R gene (the whole CDS of 954 bp, 14 bp of the 5'-untranslated region and 380 bp of the 3'-untranslated region) in 25 rabbits of different breeds including 10 animals of the Japanese and 10 animals of the Rhinelander breeds. The sequences obtained from Burgundy Fawn and Thuringian rabbits confirmed the presence of the 30 bp in-frame deletion (c.304_333del30) as the determinant of the e allele (EMBL accession number FN658678). The c.280_285del6 was identified in a Checkered Giant rabbit carrying the ED(or ES) allele confirming what we previously reported  (EMBL accession number FN658677). Surprisingly, another 6 bp deletion flanked by a G>A transition in 5' (c.[124G>A;125_130del6]), was identified in all Japanese and Rhinelander sequenced rabbits (EMBL accession number FN658679). The 6 bp deletion included 2 nucleotides of codon 42, three nucleotide of codon 43 and one nucleotide of codon 44. Together with the G>A transition this 6 bp deletion eliminates two amino acids in the first transmembrane domain (D and G at positions 42 and 43 of the wild type sequences) and, in addition, causes an amino acid substitution at position 44 (p.L44T), considering the wild type sequences (Figure 1). EMBL accession numbers of the wild type sequences (obtained from Giant Grey and Blue Vienna rabbits; see ) are FN658675 and FN658676.
Genotyping in different rabbit breeds
In order to confirm the putative role of this new allele in determining the Japanese and tricolour phenotypes we genotyped by fragment analysis the c.125_130del6 mutation in a total of 371 rabbits belonging to 31 breeds. Genotyping data are reported in Table 1 (see Figures 3E and 3F as examples of the genotyping results). These animals were also genotyped for the other two in-frame deletions associated with the e and ED(or ES) alleles . All rabbits of the Japanese, Rhinelander and Dutch tricolour breeds were homozygous for the new composite mutation. This mutation was also identified: 1) in heterozygous condition with the c.280_285del6 allele, in two black and white Checkered Giant rabbits belonging to a colony in which there was segregation of the tricolour and black and white phenotypes; 2) in heterozygous condition with a wild type Extension allele in one Giant Grey rabbit (with wild type coat colour); 3) homozygous in one Angora albino rabbit. The presence of this mutation in other animals that do not have the Japanese or tricolour phenotype demonstrates the dominance of the c.280_285del6 and of the wild type Extension alleles over the c.[124A;125_130del6] allele, confirming the results of classical genetic studies [1, 33, 36, 37]. In addition, albino rabbits do not express any other coat colour locus allele due to its status caused by a mutation in the tyrosinase (TYR) gene . All rabbits were also genotyped for the insertion in the ASIP gene (c.5_6insA) we recently identified to be the causative mutation of the nonagouti black coat colour . Interestingly, all three possible genotypes were identified in the Japanese and Rhinelander rabbits whereas 7 Dutch tricolour rabbits were homozygous for the c.5_6insA mutation and one was heterozygous. To exclude the presence of additional ASIP mutations we sequenced the three coding exons and parts of the intronic regions in two Japanese or Rhinelander animals for each c.5_6insA genotype. No other polymorphisms were identified apart those already reported .
Family based analyses and colour segregation
Three F1 families were created crossing animals of different breeds (Additional file 1). All parental animals were homozygous for the c.5_6insA ASIP mutation. The first family was obtained by crossing a Rhinelander buck with a Thuringian doe homozygous for the c.304_333del30 mutation (e allele). Of the 7 obtained F1 rabbits, 2 had the Japanese coat colour pattern (but with predominance of the yellow/red colour) and 5 were similar to the buck (tricolour) even if, again, they showed a predominance of the yellow/red spots over the black spots. This confirms the dominance of the c.[124A;125_130del6] allele over the e allele and that there was a dose effect of the new mutation that was in heterozygous condition in the F1 animals. A second F1 family was created by crossing the same Rhinelander buck with a Japanese doe obtaining 9 F1 rabbits. Four of them had the same coat colour pattern of the buck and 5 were like the doe. In this family, the tricolour or Japanese brindling phenotypes were determined by the segregation of a dominant English spotting allele (heterozygous in the Rhinelander father), like in the previous family [3, 32, 33, 35, 42]. The third F1 family was constituted by 10 F1 rabbits obtained by crossing a Checkered Giant buck (homozygous for the c.280_285del6 allele) and the Japanese doe of the previous family. All F1 animals were completely black (no. = 3) or black and white (no. = 7). Again, the black and white F1 rabbits were determined by a dominant English spotting locus allele, heterozygous in the Checkered Giant buck [3, 32, 33, 35, 42]. This family confirms the dominance of the c.280_285del6 allele over the c.[124A;125_130del6] allele. F2 animals were obtained by crossing two F1 rabbits of this family. The F1 animals were carriers of the c.[124A;125_130del6] MC1R allele and of the English spotting locus allele. Of the nine F2 rabbits, 4 were completely black (two were homozygous for the c.280_285del6 allele and two had the heterozygous genotype c.280_285del6/c.[124A;125_130del6]), 3 were black and white (one was homozygous for the c.280_285del6 allele and two had the heterozygous genotype c.280_285del6/c.[124A;125_130del6]), 2 had a tricolour phenotype, confirming the segregation of the Japanese brindling coat colour together with homozygosity for the c.[124A;125_130del6] allele in these latter rabbits (data not shown).
MC1R gene expression in different skin regions
Several skin specimens were sampled from two rabbits (one Rhinelander rabbit and one Japanese solid rabbit) just after moulting. Skin samples were collected from regions with different hair colours (white, black and red). Total RNA was isolated, retrotranscribed and amplified to evaluate (as a simple fragment size analysis) the presence of MC1R transcripts including the c.125_130del6 mutation (see Methods). In all skin samples with black hair a transcript with expected size was identified whereas no amplification was obtained from skins of both red or white hair regions (Figures 3A and 3B). A control glyceraldehyde 3-phosphate dehydrogenase (GAPDH) cDNA fragment was amplified from all skin samples (Additional file 2). Amplification of DNA extracted from the same skin samples used for RNA analysis showed the expected MC1R gene fragment (with the c.125_130del6 mutation) (Figure 3).
According to the results we obtained, the c.[124G>A;125_130del6] composite mutation is associated with the eJallele at the rabbit Extension series. This is the third in-frame MC1R mutation associated with a coat colour in rabbit, making this species quite unique. One MC1R allele determined by a nonframeshift deletion associated only with melanism has been observed in each of other four species (jaguar, jaguarundi, squirrel and guinea fowl) [23, 26, 30].
The c.280_285del6 rabbit allele (EDor ES), that eliminates two amino acids in the second transmembrane domain), may determine a constitutive activation of the MC1R protein, that in turn, drives eumelanin production. It seems that the eJin-frame mutation that eliminates two amino acids in the first transmembrane domain (including an amino acid substitution at position 44) could cause the melanic phenotype of the black skin areas by either determining constitutive activation of the receptor or increasing the affinity to αMSH. Otherwise, the presence of black spots or black hairs in Rhinelander and Japanese rabbits could not be explained. The black of the Japanese and tricolour rabbits has been also suggested to be dominant by the pioneering studies of Castle  and Punnett . The nonagouti insertion  is not the causative mutation of the black areas in these animals for the following reasons: i) in our population survey, all three genotypes for the c.5_6insA ASIP mutation  were identified in Japanese and Rhinelander rabbits, but their coat colour did not apparently differ; ii) as the Extension locus should be epistatic over the Agouti locus, the presence of at least a wild type allele is needed to show the allelic status at the Agouti locus, and the eJallele seems clearly a non wild-type allele, at least in our crossbreeding experiments; iii) classical genetic studies indicated that the presence of any Agouti allele does not alter the black in the Japanese rabbits [36, 37]. It could be also possible that (an)other not yet identified mutation(s) in the promoter region of the ASIP gene can determine the black phenotype in rabbits not homozygous for the c.5_6insA ASIP mutation but this hypothesis is quite remote as we did not have any evidence of the presence of other non-agouti mutations in a survey that included 31 rabbit breeds . According to classical genetic studies [36, 37] and to our results, the eJallele is dominant (or partially dominant) over the e allele but recessive over the other alleles of the Extension series. Another example of recessive melanism determined by an in-frame deletion in the MC1R gene has been recently reported in guinea fowl . However, to obtain a direct experimental proof of the effects of the eJin-frame composite deletion on the receptor activities and interactions in rabbit pharmacological studies should be carried out.
It also remains to be explained the presence of red coat colour areas associated with this mutation. The effect of the rabbit eJallele resembles the mosaicism of the sex linked Orange locus in female cat and that of Sex linked yellow in female Syrian hamster [38, 39]. The Lyon hypothesis for random X inactivation  can explain the two sex linked phenotypes in heterozygous animals carrying both a wild-type allele and a mutated allele. However the rabbit Extension is an autosomal locus and the mosaic phenotype is present in homozygous eJ/eJanimals (or in eJ/e rabbits, in which the eJallele is dominant). Komai [47, 48] suggested that the Japanese brindling coat colour in rabbit could be determined by both EDand e alleles that are present on the same chromosome as a result of unequal crossing-over. Searle  re-interpreted Komai hypothesis [47, 48] from the point of view of dosage compensation as Lyon  looked at the situation of X-linked loci. Of the two alleles (EDand e) on the same chromosome only one can be active in each chromosome so as to achieve the correct gene dosage. As the EDis dominant, all cell types in which at least a copy of EDwould be active (in one or the other chromosome) produce black hair, whereas the cells in which all EDare not active (that phenotypically appear e/e) would produce yellow/red hair. The result would be a brindled effect with black predominating over red. In case of heterozygous eJ/e this hypothesis would predict that half the cells give black hair and half red hair, but as inactivation is random there could be eJ/e rabbits that are completely or almost completely red, as reported by Castle  and Punnett . According to our results, the Komai-Searle hypothesis cannot be confirmed. The eJallele is different from both the EDand e alleles that are caused by distinct nonframeshift deletions.
Another hypothesis on the presence of autosomal variegation was suggested by Chase  who indicated that this phenotype in ep/epguinea-pigs could be the result of unstable genes. Somatic mutations to black could occur in yellow mutants. In pigs, a high frequency of somatic reversion of a frameshift mutation in the MC1R gene cause black spotting in a red background . This reversion re-establishes the correct reading frame in a MC1R transcript that includes another missense mutation (p.D121N) associated with the dominant black coat colour . However, the evidence obtained for the rabbit Japanese brindling coat colour may indicate that the mechanism determining this phenotype is different from what was observed in pigs. In pigs, the 2 bp mutation determining the Epallele causes a frameshift mutation in the CDS sequence of the pig MC1R gene producing a transcript that includes a premature stop at codon 56 whereas the rabbit eJallele is determined by a nonframeshift deletion that could be predicted to maintain the MC1R action. In pigs, gene expression evaluation in red skin spots showed the presence of the mutated MC1R transcript only (including the 2 bp insertion) whereas two MC1R transcripts, one with the insertion of 2 bp and one reversed (without this insertion), were detected in black skin spots . This was not unexpected since partial expression of the non mutated transcript should be sufficient to restore the dominant black coat colour in pigs . In Japanese and Rhinelander rabbits, only one MC1R transcript (having the 6 bp deletion of the eJallele) was observed in black hair skin regions whereas red regions did not express any MC1R transcript (at least at the level of sensitivity of our qualitative RT-PCR). However, the MC1R gene (including the 6 bp deletion of the eJallele) was amplified from DNA isolated from black and red skin regions indicating that somatic disrupting mutations did not occur in pheomelanic regions. Even if these results were obtained in a few animals, based on these evidences it could be possible to suppose that a regulatory mechanism, driven by another mutation in linkage disequilibrium with the c.[124G>A;125_130del6] composite mutation or by the same composite mutation, may determine the Japanese brindling phenotype in rabbits. To identify a promoter polymorphism we sequenced about 1 kb of the MC1R gene 5'-flanking region in two Japanese and in two Checkered Giant rabbits but we did not find any sequence difference (data not shown). Regulation of MC1R expression might be epigenetically determined as already reported for other genes [50, 51]. Constitutive methylation of the MC1R promoter could inhibits MC1R expression causing red coat colour. On the other hand, DNA hypomethylation in some cellular skin clones could lead MC1R expression that might determine black coat colour. In this case, as also similarly hypothesized in pigs, since the expression of one of the two copies is sufficient to restore a dominant black coat colour , this mechanism could act by chance in one or the other gene promoter of the diploid somatic skin cells. A dose-effect mechanism that could confirm this hypothesis was evident in the eJ/e tricolour rabbits (in which the e allele produces a non functional MC1R receptor ) obtained crossing a Rhinelander buck with a Thuringian doe (Additional file 3). In these rabbits there is a preponderance of the red spots over the black spots compared to tricolour eJ/eJrabbits. A similar observation has been reported by Castle  and Punnett  who even showed that, as an extreme phenotype, some eJ/e rabbits did not have black hair. It is worthwhile to mention that the data in our hands cannot formally exclude the possibility that the Japanese coat colour is caused by a frequent somatic mutation in a regulatory region that re-establishes the MC1R gene expression. However, this mechanism cannot completely explain the Japanese brindling in solid coloured rabbits. These animals could have striped-like patterns of black and red as well as brindle-like phenotypes in specific body regions for which a directional epigenetic activation or silencing could be a plausible explanation instead of a random somatic mutation reverting the phenotype in clone-derived cell lines, as in the case of the Epallele in spotted pigs.
We have identified a novel composite mutation in the rabbit MC1R gene associated with the Japanese brindling coat colour phenotype that is a mixture of yellow/red and black hairs. In the case of the epistatic interaction with another locus (English spotting), it assumes the tricolour pattern with well defined red and black spots. This mutation might be responsible for a dominant (over allele e) or recessive (over other Extension alleles) eumelanin production in the black hair skin regions. However, the mechanism(s) that cause(s) the presence of both red and black hairs in the same animal seems more complex. Results we obtained about the expression of the eJallele indicate that a regulatory mechanism, driven by another mutation in linkage disequilibrium or by the same identified composite mutation, could determine the brindling phenotype. Epigenetics could be involved in determining this coat colour phenotype. Other studies are needed to evaluate this hypothesis.
Animals and samples
Blood samples were collected from 10 unrelated Japanese, 10 unrelated Rhinelander and one Burgundy Fawn, one Thuringian, one Checkered Giant, one Giant Grey and one Blue Vienna rabbits registered to the breed herdbook (all these animals were used for sequencing as reported below). Hair root samples were also obtained from: 22 Japanese, 1 Rhinelander and 8 Dutch tricolour rabbits. Three F1 families were constructed crossing 1) a Rhinelander buck with a Thuringian doe (obtaining 7 F1 rabbits), 2) the same Rhinelander buck with a Japanese doe (obtaining 9 F1 rabbits) and 3) a Checkered Giant buck (with blue spots) with the same Japanese doe (obtaining 10 F1 rabbits). Two F1 rabbits of the latter family were crossed producing one F2 family (9 rabbits). Blood samples were collected from all these animals after slaughtering at slaughterhouse and stored at -20°C. Pictures were taken from all F1, F2 and parental rabbits. Skin samples were collected at the slaughterhouse from one Japanese and two Rhinelander rabbits slaughtered just after moulting. These samples were obtained from completely red, black and, in Rhinelander rabbits, also from completely white hair skin regions. In the Japanese rabbit used for skin sampling it was possible to identify completely red and completely black hair regions. Skin specimens were immediately snap frozen in liquid nitrogen and stored at -80°C.
Hair root samples were also collected from additional 315 rabbits belonging to other 29 different breeds (including the Dutch breed, for which animals with other coat colours were sampled) (Table 1). All these rabbits were registered to their breed herdbooks.
Sequencing of the MC1R and ASIP genes
DNA was extracted from blood and skin samples using the Wizard® Genomic DNA Purification kit (Promega Corporation, Madison, WI) or using rapid extraction methods for hair roots [35, 52]. Primers for rabbit MC1R amplification and sequencing (Additional file 4) were already reported  or newly designed using PRIMER3 (Whitehead Institute for Biomedical Research, Cambridge, MA) on consensus gene regions obtained aligning the MC1R gene in different species or obtained using rabbit trace record sequences belonging to the MC1R gene identified by BLASTN analysis. Primers for ASIP sequencing (the three coding exons and part of intronic regions) were previously reported . MC1R gene sequences were obtained from the 10 Japanese and 10 Rhinelander sampled rabbits and from one Burgundy Fawn, one Thuringian, one Giant Grey, one Vienna Blue and one Checkered Giant rabbits already genotyped for the previously detected in-frame deletions [35, 52]. ASIP gene sequences were obtained for two Japanese or Rhinelander rabbits for each c.5_6insA ASIP genotype, as previously determined . PCR was performed using a TGradient thermal cycler (Biometra, Goettingen, Germany) or a PT-100 thermal cycler (MJ Research, Watertown, MA, USA) in a volume of 20 μL containing 10-100 ng DNA template, 1 U DNA EuroTaq DNA polymerase (EuroClone Ltd., Paington, Devon, UK), 1× PCR Buffer, 2.5 mM dNTPs, 10 pmol of each primer and optimised MgCl2 concentrations (from 2.0 to 2.5 mM). PCR profile was as follows: 5 min at 95°C; 35 amplification cycles of 30 s at 95°C, 30 s at 60/65°C, 30 s at 72°C; 5 min at 72°C. For sequencing, 3-5 μL of PCR product was treated with 2 μL of ExoSAP-IT® (USB Corporation, Cleveland, Ohio, USA) following the manufacturer's protocol. Cycle sequencing of the PCR products was obtained with the Big Dye v3.1 kit (Applied Biosystems, Foster City, CA, USA) and sequencing reactions, after a few purification steps using EDTA 0.125 M, Ethanol 100% and Ethanol 70%, were loaded on an ABI3100 Avant sequencer (Applied Biosystems). All sequences were visually inspected, edited, assembled, and aligned with the help of the BioEdit software v. 18.104.22.168 http://www.mbio.ncsu.edu/BioEdit/bioedit.html and the CodonCode Aligner software http://www.codoncode.com/aligner.
Analysis of the novel 6 bp deletion was obtained by amplification of a portion of the MC1R gene (Additional file 4) and by capillary electrophoresis for fragment analysis. Genotyping of this mutation was carried out using DNA samples extracted from blood or hair root samples collected from the 371 rabbits of different breeds and from the rabbits of the three F1 and F2 families. Genotypes were also obtained from DNA extracted from skin samples. All these animals were also genotyped for the other two in-frame MC1R deletions associated with the e and ED(or ES) alleles already identified [35, 52] and by the ASIP insertion associated with the black nonagouti coat colour we recently reported . To analyse the novel 6 bp deletion, PCR profile was as described above and as reported in Additional file 4. Amplification was obtained in a final volume of 10 μL with the forward primer labelled with 6-FAM at 5' (Additional file 4). Amplified products were loaded on a ABI3100 Avant capillary sequencer (Applied Biosystems). For this analysis 1-2 μL of reaction product was diluted in 10 μL of Hi-Di formamide (Applied Biosystems) and added to 0.1 μL of Rox labelled DNA ladders (500HD Rox, Applied Biosystems). Labelled DNA fragments were sized using GeneScan v. 3.7 and Genotyper v. 3.7 software (Applied Biosystems).
RNA extraction and RT-PCR
Isolation of total RNA from rabbit skin specimens (about 100 mg) was carried out using the RNeasy® Lipid Tissue kit (Qiagen) following the manufacturer's instructions. After RNA extraction, about 1 μg of total RNA was treated by RNase-Free DNase set (Qiagen) and 40 ng were retrotranscribed with Improm II Reverse Transcription system (Promega) using oligo(dT) primers and following the manufacturer's protocol. cDNA was amplified with the same primers used for genotyping the 6 bp novel MC1R deletion as reported above and in Additional file 4. Two-three μL of amplified fragments were loaded in a ABI3100 Avant capillary sequencer (Applied Biosystems) as described above. GAPDH cDNA amplification was used as reference (Additional files 2 and 4). GAPDH fragments were electrophoresed in 10% polyacryamide:bisacrylamide 29:1 TBE 1× gels and visualized with 1× GelRed Nucleid Acid Gel Stain (Biotium Inc., Hayward, CA, USA). PCR fragments of the MC1R gene obtained with non labelled primers were also confirmed by sequencing as described above.
Bennett DC, Lamoreux ML: The color loci of mice - a genetic century. Pigment Cell Res. 2003, 16: 333-344. 10.1034/j.1600-0749.2003.00067.x.
Montoliu L, Oetting WS, Bennett DC: Color Genes. 2009, European Society for Pigment Cell Research, [http://www.espcr.org/micemut]
Searle AG: Comparative Genetics of Coat Colour in Mammals. 1968, London: Logos Press
Robbins LS, Nadeau JH, Johnson KR, Kelly MA, Roselli-Rehfuss L, Baack E, Mountjoy KG, Cone RD: Pigmentation phenotypes of variant extension locus alleles result from point mutations that alter MSH receptor function. Cell. 1993, 72: 827-834. 10.1016/0092-8674(93)90572-8.
Bultman SJ, Michaud EJ, Woychik RP: Molecular characterization of the mouse agouti locus. Cell. 1992, 71: 1195-1204. 10.1016/S0092-8674(05)80067-4.
Lu D, Willard D, Patel IR, Kadwell S, Overton L, Kost T, Luther M, Chen W, Woychik RP, Wilkison WO, Cone RD: Agouti protein is an antagonist of the melanocyte-stimulating-hormone receptor. Nature. 1994, 371: 799-802. 10.1038/371799a0.
Ollmann MM, Lamoreux ML, Wilson BD, Barsh GS: Interaction of Agouti protein with the melanocortin 1 receptor in vitro and in vivo. Genes Dev. 1998, 12: 316-330. 10.1101/gad.12.3.316.
Valverde P, Healy E, Jackson I, Rees JL, Thody AJ: Variants of the melanocyte-stimulating hormone receptor gene are associated with red hair and fair skin in humans. Nat Genet. 1995, 11: 328-330. 10.1038/ng1195-328.
Cone RD, Lu D, Koppula S, Våge DI, Klungland H, Boston B, Chen W, Orth DN, Pouton C, Kesterson RA: The melanocortin receptors: agonists, antagonists, and the hormonal control of pigmentation. Recent Prog Horm Res. 1996, 51: 287-317.
Klungland H, Våge DI, Gomez-Raya L, Adalsteinsson S, Lien S: The role of melanocyte-stimulating hormone (MSH) receptor in bovine coat color determination. Mamm Genome. 1995, 6: 636-639. 10.1007/BF00352371.
Joerg H, Fries HR, Meijerink E, Stranzinger GF: Red coat color in Holstein cattle is associated with a deletion in the MSHR gene. Mamm Genome. 1996, 7: 317-318. 10.1007/s003359900090.
Rouzaud F, Martin J, Gallet PF, Delourme D, Goulemot-Leger V, Amigues Y, Ménissier F, Levéziel H, Julien R, Oulmouden A: A first genotyping assay of French cattle breeds based on a new allele of the extension gene encoding the melanocortin-1 receptor (Mc1r). Genet Sel Evol. 2000, 32: 511-520. 10.1186/1297-9686-32-5-511.
Kijas JMH, Wales R, Törnsten A, Chardon P, Moller M, Andersson L: Melanocortin receptor 1 (MC1R) mutations and coat color in pigs. Genetics. 1998, 150: 1177-1185.
Kijas JM, Moller M, Plastow G, Andersson L: A frameshift mutation in MC1R and a high frequency of somatic reversions cause black spotting in pigs. Genetics. 2001, 158: 779-785.
Marklund L, Johansson Moller M, Sandberg K, Andersson L: A missense mutation in the gene for melanocyte-stimulating hormone receptor (MC1R) is associated with the chestnut coat color in horses. Mamm Genome. 1996, 7: 895-899. 10.1007/s003359900264.
Våge DI, Klungland H, Lu D, Cone RD: Molecular and pharmacological characterization of dominant black coat color in sheep. Mamm Genome. 1999, 10: 39-43. 10.1007/s003359900939.
Fontanesi L, Beretti F, Riggio V, Dall'Olio S, Gómez González E, Finocchiaro R, Davoli R, Russo V, Portolano B: Missense and nonsense mutations in melanocortin 1 receptor (MC1R) gene of different goat breeds: association with red and black coat colour phenotypes but with unexpected evidences. BMC Genet. 2009, 10: 47-10.1186/1471-2156-10-47.
Newton JM, Wilkie AL, He L, Jordan SA, Metallinos DL, Holmes NG, Jackson IJ, Barsh GS: Melanocortin 1 receptor variation in the domestic dog. Mamm Genome. 2000, 11: 24-30. 10.1007/s003350010005.
Everts RE, Rothuizen J, van Oost BA: 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. 10.1046/j.1365-2052.2000.00639.x.
Våge DI, Lu D, Klungland H, Lien S, Adalsteinsson S, Cone RD: A non-epistatic interaction of agouti and extension in the fox, Vulpes vulpes. Nat Genet. 1997, 15: 311-315. 10.1038/ng0397-311.
Våge DI, Fuglei E, Snipstad K, Beheim J, Landsem VM, Klungland H: Two cysteine substitutions in the MC1R generate the blue variant of the Arctic fox (Alopex lagopus) and prevent expression of the white winter coat. Peptides. 2005, 26: 1814-1817. 10.1016/j.peptides.2004.11.040.
Ritland K, Newton C, Marshall HD: Inheritance and population structure of the white-phased "Kermode" black bear. Curr Biol. 2001, 11: 1468-1472. 10.1016/S0960-9822(01)00448-1.
Eizirik E, Yuhki N, Johnson WE, Menotti-Raymond M, Hannah SS, O'Brien SJ: Molecular genetics and evolution of melanism in the cat family. Curr Biol. 2003, 13: 448-453. 10.1016/S0960-9822(03)00128-3.
Nachman MW, Hoekstra HE, D'Agostino SL: The genetic basis of adaptive melanism in pocket mice. Proc Natl Acad Sci USA. 2003, 100: 5268-5273. 10.1073/pnas.0431157100.
Hoekstra HE, Hirschmann RJ, Bundey RA, Insel PA, Crossland JP: A single amino acid mutation contributes to adaptive beach mouse color pattern. Science. 2006, 313: 101-104. 10.1126/science.1126121.
McRobie H, Thomas A, Kelly J: The genetic basis of melanism in the gray squirrel (Sciurus carolinensis). J Hered. 2009, 100: 709-714. 10.1093/jhered/esp059.
Kerje S, Lind J, Schütz K, Jensen P, Andersson L: Melanocortin 1-receptor (MC1R) mutations are associated with plumage colour in chicken. Anim Genet. 2003, 34: 241-248. 10.1046/j.1365-2052.2003.00991.x.
Nadeau NJ, Minvielle F, Mundy NI: Association of a Glu92Lys substitution in MC1R with extended brown in Japanese quail (Coturnix japonica). Anim Genet. 2006, 37: 287-289. 10.1111/j.1365-2052.2006.01442.x.
Theron E, Hawkins K, Bermingham E, Ricklefs RE, Mundy NI: The molecular basis of an avian plumage polymorphism in the wild: a melanocortin-1-receptor point mutation is perfectly associated with the melanic plumage morph of the bananaquit, Coereba flaveola. Curr Biol. 2001, 11: 550-557. 10.1016/S0960-9822(01)00158-0.
Vidal O, Araguas RM, Fernández R, Heras S, Pla C: Melanism in guinea fowl (Numida meleagris) is associated with a deletion of Phenilalanine-256 in the MC1R gene. Anim Genet. doi:10.1111/j.1365-2052.2010.02056.x
Rosenblum EB, Hoekstra HE, Nachman MW: Adaptive reptile color variation and the evolution of the Mc1r gene. Evolution. 2004, 58: 1794-1808.
Castle WE: The Genetics of Domestic Rabbit. 1930, London, Cambridge Harvard University Press
Robinson R: Genetic studies of the rabbit. Bibliogr Genet. 1958, 17: 229-558.
Fox RR: Taxonomy and Genetics. The Biology of the Laboratory Rabbit. Edited by: Manning PJ, Ringler DH, Newcomer CE. 1994, Academic Press, San Diego, CA, 1-26. Second
Fontanesi L, Tazzoli M, Beretti F, Russo V: Mutations in the melanocortin 1 receptor (MC1R) gene are associated with coat colours in the domestic rabbit (Oryctolagus cuniculus). Anim Genet. 2006, 37: 489-493. 10.1111/j.1365-2052.2006.01494.x.
Castle WE: Genetics of the Japanese rabbit. J Genet. 1924, 14: 225-229. 10.1007/BF02983096.
Punnett RC: On the "Japanese" rabbit. J Genet. 1924, 14: 230-240. 10.1007/BF02983097.
Schmidt-Küntzel A, Nelson G, David VA, Schäffer AA, Eizirik E, Roelke ME, Kehler JS, Hannah SS, O'Brien SJ, Menotti-Raymond M: A domestic cat X chromosome linkage map and the sex-linked orange locus: mapping of orange, multiple origins and epistasis over nonagouti. Genetics. 2009, 181: 1415-1425. 10.1534/genetics.108.095240.
Alizadeh A, Hong LZ, Kaelin CB, Raudsepp T, Manuel H, Barsh GS: Genetics of Sex-linked yellow in the Syrian hamster. Genetics. 2009, 181: 1427-1436. 10.1534/genetics.108.095018.
Wright S: Colour inheritance in mammals - V. The guinea pig. J Hered. 1917, 8: 476-480.
Wright S: The relation between piebald and tortoiseshell color patterns in guinea pigs. Anat Rec. 1923, 26: 393-
Fontanesi L, Vargiolu M, Scotti E, Mazzoni M, Clavenzani P, De Giorgio R, Romeo G, Russo V: Endothelin receptor B (EDNRB) is not the causative gene of the English spotting locus in the domestic rabbit (Oryctolagus cuniculus). Anim Genet. doi:10.1111/j.1365-2052.2010.02084.x
Marklund S, Kijas J, Rodriguez-Martinez H, Rönnstrand L, Funa K, Moller M, Lange D, Edfors-Lilja I, Andersson L: Molecular basis for the dominant white phenotype in the domestic pig. Genome Res. 1998, 8: 826-833.
Aigner B, Besenfelder U, Müller M, Brem G: Tyrosinase gene variants in different rabbit strains. Mamm Genome. 2000, 11: 700-702. 10.1007/s003350010120.
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: 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. 2010, 95: 166-175. 10.1016/j.ygeno.2009.11.003.
Lyon MF: Gene action in the X-chromosome of the mouse (Mus musculus L.). Nature. 1961, 190: 372-373. 10.1038/190372a0.
Komai T: Semi-allelic genes. Amer Nat. 1950, 84: 381-392. 10.1086/281636.
Komai T: On mosaic inheritance in mammals. Amer Nat. 1951, 85: 333-334. 10.1086/281685.
Chase HB: Studies on the tricolor pattern of the guinea-pig - II. The distribution of black and yellow as affected by white spotting and by imperfect dominance in the tortoiseshell series of alleles. Genetics. 1939, 24: 610-643.
Kass SU, Pruss D, Wolffe AP: How does DNA methylation repress transcription?. Trends Genet. 1997, 13: 444-449. 10.1016/S0168-9525(97)01268-7.
Shen L, Kondo Y, Guo Y, Zhang J, Zhang L, Ahmed S, Shu J, Chen X, Waterland RA, Issa JP: Genome-wide profiling of DNA methylation reveals a class of normally methylated CpG island promoters. PLoS Genet. 2007, 3: 2023-2036. 10.1371/journal.pgen.0030181.
Fontanesi L, Tazzoli M, Russo V: Non-invasive and simple methods for sampling DNA for PCR analysis of melanocortin 1 receptor (MC1R) gene mutations: a technical note. World Rabbit Sci. 2007, 15: 121-126.
We thank people who helped in the sampling: a particular thank is for Renate Regitz (Germany) who provided many Japanese and the Dutch tricolour samples. Associazione Nazionale Coniglicoltori Italiani (ANCI) (Italy) collaborated during the collection of samples. This work was financed by University of Bologna RFO and FAGenomicH funds to LuF and was supported by University of Limoges and INRA.
LuF conceived the study, analysed sequences and obtained data, contributed to the sampling and to the construction of rabbit families, coordinated and organized the laboratory work and drafted the manuscript. ES reared the rabbit families, genotyped and collaborated in the RT-PCR analyses. MC performed the RT-PCR experiment. FB and LiF sequenced and genotyped. SD, DA and SDe collected samples. DA, VR supervised the work and were involved in the design of the study. AO co-conceived the study, contributed to the sampling and to drafting the manuscript. All authors reviewed the manuscript and accepted the final version.