Characterization of the effect of Melanocortin 1 Receptor,a member of the hair color genetic locus,in alpaca (Lama pacos) fleece color differentiation

Little is known about the inheritance and influence of the fleece color gene Melanocortin 1 Receptor (MC1R). Melanocortin 1 Receptor (MC1R) is a well-known gene responsible for red versus black fleece pigmentation and is hypothesized to be a candidate gene for variation in alpaca coloration patterns. Inheritance of red versus black pigmentation in the context of genetic mutation is well understood in many domesticated mammals. We characterized the MC1R gene in a population of multi-colored alpacas in order to better understand its effect on coat color in the alpaca. Our characterization of the alpaca MC1R gene revealed 11 mutations. Of these one is a 4 bp deletion, four are silent mutations and six are single nucleotide polymorphisms (SNPs) that alter the amino acid sequence (T28V, M87V, S126G, T128I, S196F, R301C). No mutation correlated completely with fleece color in alpacas at the MC1R locus. This may be due to the epistatic relationship of MC1R with other coat color genes especially agouti signaling protein (ASIP).     

Mutations in ASIP and MC1R: dominant black and recessive black alleles segregate in native Swedish sheep populations

By studying genes associated with coat colour, we can understand the role of these genes in pigmentation but also gain insight into selection history. North European short‐tailed sheep, including Swedish breeds, have variation in their coat colour, making them good models to expand current knowledge of mutations associated with coat colour in sheep. We studied ASIP and MC1R, two genes with known roles in pigmentation, and their association with black coat colour. We did this by sequencing the coding regions of ASIP in 149 animals and MC1R in 129 animals from seven native Swedish sheep breeds in individuals with black, white or grey fleece. Previously known mutations in ASIP [recessive black allele: g.100_105del (D₅) and/or g.5172T>A] were associated with black coat colour in Klövsjö and Roslag sheep breeds and mutations in both ASIP and MC1R (dominant black allele: c.218T>A and/or c.361G>A) were associated with black coat colour in Swedish Finewool. In Gotland, Gute, Värmland and Helsinge sheep breeds, coat colour inheritance was more complex: only 11 of 16 individuals with black fleece had genotypes that could explain their black colour. These breeds have grey individuals in their populations, and grey is believed to be a result of mutations and allelic copy number variation within the ASIP duplication, which could be a possible explanation for the lack of a clear inheritance pattern in these breeds. Finally, we found a novel missense mutation in MC1R (c.452G>A) in Gotland, Gute and Värmland sheep and evidence of a duplication of MC1R in Gotland sheep.     

Nucleotide diversity of the melanocortin 1 receptor gene (MC1R) in the gayal (Bos frontalis)

The melanocortin 1 receptor gene (MC1R) plays a crucial role in determining coat colour of mammals. To investigate the relationship of polymorphism of the MC1R with coat colour in gayal, the coding sequence (CDS), and the 5'- and 3'-untranslated regions (UTR) of the MC1R were sequenced from 63 samples from the gayal and compared with the sequences of the MC1R from other ruminant species. A sequence of 1,136 bp including the whole CDS (954 bp) and parts of the 5'- and 3'-UTR (164 and 18 bp, respectively) of the gayal MC1R was obtained. A total of nine single nucleotide polymorphisms (SNPs) including four SNPs (c.-129T>C, c.-127A>C, c.-106C>T, c.-1G>A) in the 5'-UTR and five SNPs (c.201C>T, c.583C>T, c.663T>C, c.871A>G and c.876T>C) in the CDS were detected, revealing high genetic diversity. Three novel coding SNPs including c.201C>T, c.583C>T and c.876T>C, which have not been reported previously in bovid species, were retrieved. Within five coding SNPs, c.201C>T, c.663T>C and c.876T>C were silent mutations, while c.583C>T and c.871A>G were mis-sense mutations, resulting in changes in the amino acids located in the fifth (p.L195F) and seventh (p.T291A) transmembrane regions, respectively. The alignment of amino acid sequences was found to be very similar to those for other bovid species. It was demonstrated, using the functional effect prediction, that the p.T291A amino acid replacement could have an effect on MC1R protein function but not for the p.L195F substitution. Using phylogenetic analyses it was revealed that the gayal has a close genetic relationship with the yak. However, three classical bovine MC1R loci the E (D), E (+) and e were not retrieved in the gayal, indicating other genes or factors could affect coat colour in this species.     

Molecular population genetics of the melanic plumage polymorphism in Arctic skuas (Stercorarius parasiticus): evidence for divergent selection on plumage colour

The Arctic skua (Stercorarius parasiticus) is a classic example of an avian plumage polymorphism, with variation in melanin‐based ventral plumage coloration defining pale, intermediate and dark morphs in adults of both sexes. However, despite several decades of field research, there is an incomplete understanding of how the polymorphism in ventral plumage colour is maintained and the selective forces involved. Here, we investigate selection on a locus (MC1R) that is strongly associated with plumage colour variation in Arctic skuas using patterns of nucleotide variation and comparison to neutral loci (nuclear introns and mtDNA). We find that three linked nonsynonymous mutations in MC1R, including the single mutation described previously, are associated with plumage colour in the Arctic skua. The position of nonsynonymous mutations on a MC1R haplotype network implies that divergent selection drove the initial evolution of the colour morphs. Comparisons of F_STs of MC1R vs. nuclear introns among five skua populations differing in proportion of dark morphs along an approximate north–south cline reveal a signature of divergent selection on MC1R. In contrast, we find limited evidence for balancing selection on MC1R within populations, although the power is low. Our results provide strong evidence for both past and ongoing selection on MC1R, and, by implication, plumage colour in Arctic skuas. The results suggest that a fruitful avenue for future ecological studies will be analysis of selection on morphs in colonies at the extremes along the morph ratio cline.     

A missense mutation in melanocortin 1 receptor is associated with the red coat colour in donkeys

The seven donkey breeds recognised by the French studbook are characterised by few coat colours: black, bay and grey. Normand bay donkeys seldom give birth to red foals, a colour more commonly seen and recognised in American miniature donkeys. Red resembles the equine chestnut colour, previously attributed to a mutation in the melanocortin 1 receptor gene (MC1R). We used a panel of 124 donkeys to identify a recessive missense c.629T>C variant in MC1R that showed a perfect association with the red coat colour. This variant leads to a methionine to threonine substitution at position 210 in the protein. We showed that methionine 210 is highly conserved among vertebrate melanocortin receptors. Previous in silico and in vitro analyses predicted this residue to lie within a functional site. Our in vivo results emphasised the pivotal role played by this residue, the alteration of which yielded a phenotype fully compatible with a loss of function of MC1R. We thus propose to name the c.629T>C allele in donkeys the e allele, which further enlarges the panel of recessive MC1R loss‐of‐function alleles described in animals and humans.     

Association of MC3R gene polymorphisms with body weight in the red fox and comparative gene organization in four canids

There are five genes encoding melanocortin receptors. Among canids, the genes have mainly been studied in the dog (MC1R, MC2R and MC4R). The MC4R gene has also been analysed in the red fox. In this report, we present a study of chromosome localization, comparative sequence analysis and polymorphism of the MC3R gene in the dog, red fox, arctic fox and Chinese raccoon dog. The gene was localized by FISH to the following chromosome: 24q24‐25 in the dog, 14p16 in the red fox, 18q13 in the arctic fox and NPP4p15 in the Chinese raccoon dog. A high identity level of the MC3R gene sequences was observed among the species, ranging from 96.0% (red fox – Chinese raccoon dog) to 99.5% (red fox – arctic fox). Altogether, eight polymorphic sites were found in the red fox, six in the Chinese raccoon dog and two in the dog, while the arctic fox appeared to be monomorphic. In addition, association of several polymorphisms with body weight was analysed in red foxes (the number of genotyped animals ranged from 319 to 379). Two polymorphisms in the red fox, i.e. a silent substitution c.957A>C and c.*185C>T in the 3′‐flanking sequence, showed a significant association (P < 0.01) with body weight.     

A non‐synonymous SNP in exon 3 of the KIT gene is responsible for the classic grey phenotype in alpacas (Vicugna pacos)

The alpaca classic grey phenotype is of particular interest to the industry. Until now, there were only indirect data suggesting that the KIT gene was involved in the classic grey phenotype. All exons of KIT in three black and three classic silvergrey alpacas were sequenced. Five non‐synonymous SNPs were observed. There was only one SNP found that was present only in the silvergrey alpacas, and this was also the only SNP predicted to be damaging. This variant results in a change of a glycine (Gly) to an arginine (Arg) at amino acid position 126 (c.376G>A), occurring in the second Ig‐like domain of the extracellular domain of KIT. Basic protein modelling predicted that this variant is likely destabilising. Therefore, an additional 488 alpacas were genotyped for this SNP using the tetra‐primer amplification refractory mutation system PCR (Tetra‐primer ARMS‐PCR). All classic grey alpacas were observed to be heterozygous, and 99.3% of non‐grey dark base colour alpacas were found to be homozygous for the wildtype allele in this position. These results confirm that the classic grey phenotype in alpacas is the result of a c.376G>A (p.Gly126Arg) SNP in exon 3 of KIT. These data also support the hypothesis that the grey phenotype is autosomal dominant and that the mutation is most likely homozygous lethal.     

High follicle density does not decrease sweat gland density in Huacaya alpacas

When exposed to high ambient temperatures, mammals lose heat evaporatively by either sweating from glands in the skin or by respiratory panting. Like other camelids, alpacas are thought to evaporate more water by sweating than panting, despite a thick fleece, unlike sheep which mostly pant in response to heat stress. Alpacas were brought to Australia to develop an alternative fibre industry to sheep wool. In Australia, alpacas can be exposed to ambient temperatures higher than in their native South America. As a young industry there is a great deal of variation in the quality and quantity of the fleece produced in the national flock. There is selection pressure towards animals with finer and denser fleeces. Because the fibre from secondary follicles is finer than that from primary follicles, selecting for finer fibres might alter the ratio of primary and secondary follicles. In turn the selection might alter sweat gland density because the sweat glands are associated with the primary follicle. Skin biopsy and fibre samples were obtained from the mid-section of 33 Huacaya alpacas and the skin sections were processed into horizontal sections at the sebaceous gland level. Total, primary, and secondary follicles and the number of sweat gland ducts were quantified. Fibre samples from each alpaca were further analysed for mean fibre diameter. The finer-fibred animals had a higher total follicle density (P<0.001) and more sweat glands (P<0.001) than the thicker-fibred animals. The fibre diameter and total follicle density were negatively correlated (R²=0.56, P<0.001). Given that the finer-fibred animals had higher follicle density and more sweat glands than animals with thicker fibres, we conclude that alpacas with high follicle density should not be limited for potential sweating ability.     

Increased constraints on MC4R during primate and human evolution

The melanocortin 4 receptor (MC4R) is routinely investigated for the role it plays in human obesity, as mutations in MC4R are the most common dominantly inherited form of the disease. As little is known about the evolutionary history of this locus, we investigated patterns of variation at MC4R in a worldwide sample of 1,015 humans from 51 populations, and in 8 central chimpanzees. There is a significant paucity of diversity at MC4R in humans, but not in chimpanzees. The spectrum of mutations in humans, combined with the overall low level of diversity, suggests that most (if not all) of the observed non-synonymous polymorphisms are likely to be transient deleterious mutations. The MC4R coding region was resequenced in 12 primate species and sequences from an additional 29 vertebrates were included in molecular evolutionary analyses. MC4R is highly conserved throughout vertebrate evolution, and has apparently been subject to high levels of continuous purifying selection that increased approximately threefold during primate evolution. Furthermore, the strong selection extends to codon usage bias, where most silent mutations are expected to be either quickly fixed or removed from the population, which may help explain the unusually low levels of silent polymorphisms in humans. Finally, there is a significant tendency for non-synonymous mutations that impact MC4R function to occur preferentially at sites that are identified by evolutionary analyses as being subject to very strong purifying selection. The information from this study should help inform future epidemiological investigations of MC4R. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Loss‐of‐Function Mutations in MC4R Are Very Rare in the Greek Severely Obese Adult Population

Melanocortin‐4 receptor (MC4R) loss‐of‐function mutations are the commonest genetic cause of human monogenic obesity, so far. The contribution of MC4R coding mutations to severe obesity in the high‐obesity prone Greek population has not been investigated to date. We determined the MC4R coding sequence of 510 obese and 469 lean control subjects of Greek origin, and we estimated the prevalence and the penetrance on obesity of MC4R loss‐of‐function mutations. The functional impact of novel nonsynonymous variants detected was investigated in vitro. We found two novel synonymous mutations (L23L and I102I), four nonsynonymous mutations (T112M, S127L, N274S, and S295L), and two polymorphisms (V103I and I251L) previously described in literature. We also detected a novel mutation (L207V) in a severely obese 69‐year‐old female patient, although the mutation did not cosegregate with obesity in the corresponding pedigree and had no functional consequences on MC4R protein function. Loss‐of‐function mutations represented 75% of all nonsynonymous rare mutations identified among lean carriers and only 25% among obese subjects (P = 0.0001). The prevalence of loss‐of‐function mutations was lower in the obese group than in lean control subjects (0.20 vs. 0.64%) but this difference was not significant. Therefore, the estimated penetrance of deleterious MC4R mutations was very low (6.3%) in heterozygous Greek carriers of MC4R loss‐of‐function mutations. Our data suggest that MC4R loss‐of‐function mutations are rare in the Greek population. MC4R genetic deficiency is unlikely to explain the high propensity to develop severe obesity in this specific population.     

Investigation of a genome wide association signal for obesity: synthetic association and haplotype analyses at the melanocortin 4 receptor gene locus

Background

Independent genome-wide association studies (GWAS) showed an obesogenic effect of two single nucleotide polymorphisms (SNP; rs12970134 and rs17782313) more than 150 kb downstream of the melanocortin 4 receptor gene (MC4R). It is unclear if the SNPs directly influence MC4R function or expression, or if the SNPs are on a haplotype that predisposes to obesity or includes functionally relevant genetic variation (synthetic association). As both exist, functionally relevant mutations and polymorphisms in the MC4R coding region and a robust association downstream of the gene, MC4R is an ideal model to explore synthetic association.

Methodology/Principal Findings

We analyzed a genomic region (364.9 kb) encompassing the MC4R in GWAS data of 424 obesity trios (extremely obese child/adolescent and both parents). SNP rs12970134 showed the lowest p-value (p = 0.004; relative risk for the obesity effect allele: 1.37); conditional analyses on this SNP revealed that 7 of 78 analyzed SNPs provided independent signals (p≤0.05). These 8 SNPs were used to derive two-marker haplotypes. The three best (according to p-value) haplotype combinations were chosen for confirmation in 363 independent obesity trios. The confirmed obesity effect haplotype includes SNPs 3′ and 5′ of the MC4R. Including MC4R coding variants in a joint model had almost no impact on the effect size estimators expected under synthetic association.

Conclusions/Significance

A haplotype reaching from a region 5′ of the MC4R to a region at least 150 kb from the 3′ end of the gene showed a stronger association to obesity than single SNPs. Synthetic association analyses revealed that MC4R coding variants had almost no impact on the association signal. Carriers of the haplotype should be enriched for relevant mutations outside the MC4R coding region and could thus be used for re-sequencing approaches. Our data also underscore the problems underlying the identification of relevant mutations depicted by GWAS derived SNPs.     

Prediction of alpaca fibre quality by near-infrared reflectance spectroscopy

Rapid and efficient methods to evaluate variables associated with fibre quality are essential in animal breeding programs and fibre trade. Near-infrared reflectance spectroscopy (NIRS) combined with multivariate analysis was evaluated to predict textile quality attributes of alpaca fibre. Raw samples of fibres taken from male and female Huacaya alpacas (n = 291) of different ages and colours were scanned and their visible–near-infrared (NIR; 400 to 2500 nm) reflectance spectra were collected and analysed. Reference analysis of the samples included mean fibre diameter (MFD), standard deviation of fibre diameter (SDFD), coefficient of variation of fibre diameter (CVFD), mean fibre curvature (MFC), standard deviation of fibre curvature (SDFC), comfort factor (CF), spinning fineness (SF) and staple length (SL). Patterns of spectral variation (loadings) were explored by principal component analysis (PCA), where the first four PC's explained 99.97% and the first PC alone 95.58% of spectral variability. Calibration models were developed by modified partial least squares regression, testing different mathematical treatments (derivative order, subtraction gap, smoothing segment) of the spectra, with or without applying spectral correction algorithms (standard normal variate and detrend). Equations were selected through one-out cross-validation according to the proportion of explained variance (R²CV), root mean square error in cross-validation (RMSECV) and the residual predictive deviation (RPD), which relates the standard deviation of the reference data to RMSECV. The best calibration models were accomplished when using the NIR region (1100 to 2500 nm) for the prediction of MFD and SF, with R²CV = 0.90 and 0.87; RMSECV = 1.01 and 1.08 μm and RPD = 3.13 and 2.73, respectively. Models for SDFD, CVFD, MFC, SDFC, CF and SL had lower predictive quality with R²CV < 0.65 and RPD < 1.5. External validation performed for MFD and SF on 91 samples was slightly poorer than cross-validation, with R² of 0.86 and 0.82, and standard error of prediction of 1.21 and 1.33 μm, for MFD and SF, respectively. It is concluded that NIRS can be used as an effective technique to select alpacas according to some important textile quality traits such as MFD and SF.     

The albinism of the feral Asinara white donkeys (Equus asinus) is determined by a missense mutation in a highly conserved position of the tyrosinase (TYR) gene deduced protein

A feral donkey population (Equus asinus), living in the Asinara National Park (an island north‐west of Sardinia, Italy), includes a unique white albino donkey subpopulation or colour morph that is a major attraction of this park. Disrupting mutations in the tyrosinase (TYR) gene are known to cause recessive albinisms in humans (oculocutaneous albinism Type 1; OCA1) and other species. In this study, we analysed the donkey TYR gene as a strong candidate to identify the causative mutation of the albinism of these donkeys. The TYR gene was sequenced from 13 donkeys (seven Asinara white albino and six coloured animals). Seven single nucleotide polymorphisms were identified. A missense mutation (c.604C>G; p.His202Asp) in a highly conserved amino acid position (even across kingdoms), which disrupts the first copper‐binding site (CuA) of functional protein, was identified in the homozygous condition (G/G or D/D) in all Asinara white albino donkeys and in the albino son of a trio (the grey parents had genotype C/G or H/D), supporting the recessive mode of inheritance of this mutation. Genotyping 82 donkeys confirmed that Asinara albino donkeys had genotype G/G whereas all other coloured donkeys had genotype C/C or C/G. Across‐population association between the c.604C>G genotypes and the albino coat colour was highly significant (P = 6.17E?18). The identification of the causative mutation of the albinism in the Asinara white donkeys might open new perspectives to study the dynamics of this putative deleterious allele in a feral population and to manage this interesting animal genetic resource.     

Evolution of dark colour in toucans (Ramphastidae): a case of molecular adaptation?

In the last decades, researchers have been able to determine the molecular basis of some phenotypes, to test for evidence of natural selection upon them, and to demonstrate that the same genes or genetic pathways can be associated with convergent traits. Colour traits are often subject to natural selection because even small changes in these traits can have a large effect on fitness via camouflage, sexual selection or other mechanisms. The melanocortin‐1 receptor locus (MC1R) is frequently associated with intraspecific coat colour variation in vertebrates, but it has been far harder to demonstrate that this locus is involved in adaptive interspecific colour differences. Here, we investigate the contribution of the MC1R gene to the colour diversity found in toucans (Ramphastidae). We found divergent selection on MC1R in the clade represented by the genus Ramphastos and that this coincided with the evolution of darker plumage in members of this genus. Using phylogenetically corrected correlations, we show significant and specific relationships between the rate of nonsynonymous change in MC1R (dN) and plumage darkness across Ramphastidae, and also between the rate of functionally significant amino acid changes in MC1R and plumage darkness. Furthermore, three of the seven amino acid changes in MC1R that occurred in the ancestral Ramphastos branch are associated with melanism in other birds. Taken together, our results suggest that the dark colour of Ramphastos toucans was related to nonsynonymous substitutions in MC1R that may have been subject to positive selection or to a relaxation of selective pressure. These results also demonstrate a quantitative relationship between gene and phenotype evolution, representing an example of how MC1R molecular evolution may affect macroevolution of plumage phenotypes.     

Effect of the MC1R gene on sexual dimorphism in melanin‐based colorations

Variants of the melanocortin‐1 receptor (MC1R) gene result in abrupt, naturally selected colour morphs. These genetic variants may differentially affect sexual dimorphism if one morph is naturally selected in the two sexes but another morph is naturally or sexually selected only in one of the two sexes (e.g. to confer camouflage in reproductive females or confer mating advantage in males). Therefore, the balance between natural and sexual selections can differ between MC1R variants, as suggest studies showing interspecific correlations between sexual dimorphism and the rate of nonsynonymous vs. synonymous amino acid substitutions at the MC1R. Surprisingly, how MC1R is related to within‐species sexual dimorphism, and thereby to sex‐specific selection, has not yet been investigated. We tackled this issue in the barn owl (Tyto alba), a species showing pronounced variation in the degree of reddish pheomelanin‐based coloration and in the number and size of black feather spots. We found that a valine (V)‐to‐isoleucine (I) substitution at position 126 explains up to 30% of the variation in the three melanin‐based colour traits and in feather melanin content. Interestingly, MC1R genotypes also differed in the degree of sexual colour dimorphism, with individuals homozygous for the II MC1R variant being 2 times redder and 2.5 times less sexually dimorphic than homozygous individuals for the VV MC1R variant. These findings support that MC1R interacts with the expression of sexual dimorphism and suggest that a gene with major phenotypic effects and weakly influenced by variation in body condition can participate in sex‐specific selection processes.     

MC1R variation and melanoma risk in the Swedish population in relation to clinical and pathological parameters

The genetic background of cutaneous malignant melanoma (CMM) includes both germ line aberrations in high‐penetrance genes, like CDKN2A, and allelic variation in low‐penetrance genes like the melanocortin‐1 receptor gene, MC1R. Red‐hair colour associated MC1R alleles (RHC) have been associated with red hair, fair skin and risk of CMM. We investigated MC1R and CDKN2A variation in relation to phenotype, clinical factors and CMM risk in the Swedish population. The study cohort consisted of sporadic primary melanoma patients, familial melanoma patients and a control group. An allele‐dose dependent increase in melanoma risk for carriers of variant MC1R alleles (after adjusting for phenotype), with an elevated risk among familial CMM patients, was observed. This elevated risk was found to be significantly associated with an increased frequency of dysplastic nevi (DN) among familial patients compared to sporadic patients. MC1R variation was found to be less frequent among acral lentiginous melanomas (ALM) and dependent on tumour localisation. No association was found between CDKN2A gene variants and general melanoma risk. Two new variants in the POMC gene were identified in red haired individuals without RHC alleles.     

Differential expression of MC1R gene in Liaoning Cashmere goats with different coat colors

Melanocyte stimulating hormone receptor (MC1R) has been known as a regulator of eumelanin and phaeomelanin production in the melanocytes, and MC1R mutations causing coat color changes are known in many vertebrates; however, there are no research reports about the differentially expression of MC1R gene and its coding protein in Cashmere goats with different coat color. We examined the presence of MC1R distribution and MC1R protein and gene expression in the white Cashmere goats and black Cashmere goats, respectively; q-PCR, Western blot and immunhistochemical analysis showed that the expression of the MC1R gene in the black Cashmere goats was 3.39 fold more than the white ones (p?<?0.01), and Cashmere goats with black genotype had significantly higher (2.03, p?<?0.01) MC1R protein expression than white genotype in the all investigated samples. Moreover, all Cashmere goats with different coat color available for immunhistochemical analysis showed either lower (white Cashmere goats) or higher (black Cashmere goats) expression of the MC1R protein; these findings suggested that it had a relationship between the MC1R and the coat color of Cashmere goats. That could lay the foundation for the further research of the MC1R and coat color controllability regulation of the Cashmere goats.     

Not another type of potato: MC1R and the russet coloration of Burmese cats

The Burmese is a breed of domestic cat that originated in Southeast Asia and was further developed in the United States. Variants in melanocortin 1 receptor (MC1R) causes common coat colour phenotypes in a variety of mammalian species but only limited colour variation in the domestic cat. Known as the extension (E) locus, melanocortin 1 receptor (MC1R) interacts with the agouti locus to produce the eumelanin and pheomelanin pigments. Recently, a novel reddish coloration, which is termed russet, was identified in the Burmese cat breed. Because this russet Burmese coloration changes with aging, MC1R was suggested as candidate gene. The similar colouration in specific lineages of Norwegian Forest cat known as amber (e) (c.250G>A; p.Asp84Asn) was excluded for this Burmese phenotype. The complete 954‐bp coding region of MC1R was directly sequenced in russet Burmese and suspected carriers. A 3‐bp deletion (c.439_441del) associated with the deletion of a phenyalanine (p.Phe146del) in the protein sequence was identified. All russet coloured cats were homozygous for the variant, and all obligate carriers were heterozygous, confirming that the deletion segregated concordantly with colouring in Burmese cats from the New Zealand foundation lineage. The variant was not identified in 442 cats from 26 different breeds and random‐bred cats. Twenty‐six Burmese from the USA did not have the variant. This MC1R variant defines a unique coloration and the second breed‐specific MC1R variant in cats. The interactions of the two recessive feline MC1R alleles (E  >  e, e^r) is unknown.