Sequence Characterization of the MC1R Gene in Yak (Poephagus grunniens) Breeds with Different Coat Colors

Abstract

Melanocortin 1 receptor (MC1R) gene plays a key role in determining coat color in several species, including the cattle. However, up to now there is no report regarding the MC1R gene and the potential association of its mutations with coat colors in yak (Poephagus grunniens). In this study, we sequenced the encoding region of the MC1R gene in three yak breeds with completely white (Tianzhu breed) or black coat color (Jiulong and Maiwa breeds). The predicted coding region of the yak MC1R gene resulted of 954 bp, the same to that of the wild-type cattle sequence, with >99% identity. None of the mutation events reported in cattle was found. Comparing the yak obtained sequences, five nucleotide substitutions were detected, which defined three haplotypes ([superscript]EY1[/superscript] , [superscript]EY2[/superscript] , and [superscript]EY3[/superscript] ). Of the five mutations, two, characterizing the [superscript]EY1[/superscript] haplotype, were nonsynonymous substitutions (c.340C>A and c.871G>A) causing amino acid changes located in the first extracellular loop (p.Q114K) and in the seventh transmembrane region (p.A291T). In silico prediction might indicate a functional effect of the latter substitution. However, all three haplotypes were present in the three yak breeds with relatively consistent frequency distribution, despite of their distinguished coat colors, which suggested that there was no across-breed association between haplotypes or genotypes and black/white phenotypes, at least in the investigated breeds. Other genes may be involved in affecting coat color in the analyzed yaks.

Shi-Yi Chen 1 and Yi Huang 1 and Qing Zhu 1 and Luca Fontanesi 2 and Yong-Gang Yao 3 and Yi-Ping Liu 1

Recommended by Rongjia Zhou

1, College of Animal Science and Technology, Sichuan Agriculture University, Ya'an, Sichuan 625014, China
2, DIPROVAL, Sezione di Allevamenti Zootecnici, Faculty of Agriculture, University of Bologna, Reggio Emilia 42100, Italy
3, Key Laboratory of Animal Models and Human Disease Mechanisms of Chinese Academy of Sciences & Yunnan Province, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, Yunnan 650223, China

Received 6 February 2009; Accepted 31 March 2009

1. Introduction

The coat color of domestic animals is highly variable, and has been widely used as a unique phenotype in the morphological selection that resulted in the constitution of the modern breeds. In a variety of animals, pigmentation is mainly determined by the relative amount of two basic melanins, eumelanin (black/brown) and phaeomelanin (yellow/red), which are genetically controlled by the Extension (E) and Agouti (A) loci, respectively [1, 2]. The wild alleles at the two loci would result in different coat colors, whereas dominant or recessive mutations can alter the eumelanin/phaeomelanin ratio and subsequently cause coat color variations [3]. The E locus encodes the melanocortin receptor 1 (MC1R), which is a seven-transmembrane G-protein-coupled receptor expressed in melanocytes [4]. Mutations affecting the MC1R gene fvunction determine different coat color phenotypes in a large number of vertebrates, including livestock species such as cattle [5-7], sheep [8], goat [9], horse [10], pig [11, 12], rabbit [13], and chicken [14]. Some MC1R mutations induce eumelanin production while others prompt phaeomelanin synthesis [4].

The cattle MC1R gene is located on chromosome 18 [15] and has four principal alleles (^E+ , ^ED , ^E1 , and e) , which are defined by three mutation events affecting the coding region [5, 7]. The wild type allele (^E+ ) produces a variety of coat colors. The dominant allele (^ED ), differing from the ^E+ allele by the g.296T>C substitution, is associated with black color. The recessive allele (e ) is characterized by a frameshift mutation at position 310 and encodes a truncated protein with only three putative transmembrane domains. Animals homozygous for this allele show red/white coat color phenotypes across a large number of cattle breeds [5, 6, 16, 17]. The ^E1 allele is characterized by a 12 bp duplicated fragment at position 669 that was detected in a few breeds which have fawn to brown or grey coat colors [7, 16, 18]. Other missense and frameshift mutations have been identified in the cattle MC1R gene, but their functional effect has not been completely investigated yet [16, 19].

The yak (Poephagus grunniens ) belongs to the family of Bovidae and lives in an extremely harsh habitat (> 3000 m above the sea level). More than 90% of worldwide domestic yaks are distributed in the Qinghai-Tibet Plateau region in China. There are five registered yak breeds (Jiulong Yak, Qinghai Yak, Tianzhu Yak, Maiwa Yak, and Tibet Yak) in China as well as other local populations [20]. The Tianzhu yak is wellknown for the uniform white coat color even if a few subpopulations in this breed show small dark spots. The black coat is predominant in the other four yak breeds. More than 85% individuals of Jiulong Yak have uniform black coat, and the remaining individuals have a black coat with white spots/belt. The genetic basis of coat color in yaks is poorly studied compared to the cattle [1, 21]. Chen et al. [22] investigated coat color inheritance in yaks by means of classical genetics and suggested that the black and yellow phenotypes were mainly determined by two loci (B and D ), but the results of this report were not considered in the context of the comparative approach proposed by Searle [1]. In addition, no study has analyzed pigment or melanocyte distribution in yak, yet. Therefore it is not completely clear which could be the possible candidate genes involved in the black and yellow/white phenotype determination in yak. Here, in order to evaluate if mutations in the MC1R gene could potentially be involved in determining coat color variations in Poephagus grunniens , we characterized by sequencing this gene in 44 yaks belonging to three breeds with different coat colors: Tianzhu yaks (pure white animals), Jiulong yaks, and Maiwa yaks (black animals) (Figure 1).

Pictures of Tianzhu yak (white) and Maiwa yak (black).

[figure omitted; refer to PDF]

[figure omitted; refer to PDF]

2. Materials and Methods

2.1. Sampling

Ear tissue samples were collected from 44 yak individuals belonging to three pure indigenous breeds: 14 Tianzhu yaks, 16 Maiwa yaks, and 14 Jiulong yaks. The Chinese provinces and the number of herds from which the animals were sampled as well as their coat color appearances are reported in Table 1. Figure 1 shows yaks of white or black coat color types. Biological material was stored at -70 °C.

Table 1: Coat color and distribution of the MC1R haplotypes and genotypes in three Chinese yak breeds.

Breeds	Coat color	Location	No.	Genotypes	Haplotypes
				EY1 /EY1	EY1 /EY2	EY1 /EY3	EY2 /EY2	EY2 /EY3	EY3 /EY3	EY1	EY2	EY3
Tianzhu	White	Gansu	14	7	1	4	1	0	1	19	3	6
Maiwa	Black	Sichuan	16	8	7	1	0	0	0	24	7	1
Jiulong	Black	Sichuan	14	2	2	2	1	6	1	8	10	10
Total (%)	--	--	44	17 (38.6)	10 (22.7)	7 (15.9)	2 (4.5)	6 (13.6)	2 (4.5)	51 (58.0)	20 (22.7)	17 (19.3)

2.2. Amplification of the MC1R Gene and Sequencing

Genomic DNA was extracted by a standard phenol/ chloroform method. The whole coding region of MC1R gene was amplified using a primer pair (forward: 5[variant prime] -GGACCCTGAGAGCAAGCAC-3[variant prime] ; reverse 5[variant prime] -CTCACCTTCAGGGATGGTCTA-3[variant prime] ) designed using the cattle sequence (GenBank accession number NM_174108). The amplified fragment length was 1085 bp. Amplification of MC1R gene was performed in a 50 μ L reaction mixture containing 100 ng of DNA, 10 mM Tris-HCl (pH 8.3), 2.5 mM _MgCl2 , 50 mM KCl, 10 pmol of each primer, and 1 U of Taq polymerase (Takara, Dalian, China) using the following procedures: after a denaturation cycle at 94°C for 5 minutes, amplification was performed for 35 cycles at 94°C for 45 second, 61.5°C for 45 second, and 72°C for 1 minute, and ended with a final incubation at 72°C for 10 minutes for full extension.

All PCR products were purified on spin columns and were directly sequenced by using the amplification primers and the Big Dye Terminator v. 3.1 Cycle Sequencing Kit (Applied Biosystems, Calif, USA) on an ABI Prism 3100 DNA sequencer (Applied Biosystems) according to the manufacturer's manual. Sequencing electropherograms were visualized using the SeqMan program (DNAStar Inc., Madison, Wis, USA).

When two or more heterozygous positions were recognized in the same sequencing read, in order to clarify the haplotypes, the PCR products obtained from one individual for each heterozygous pattern were cloned after separation on 1.5% agarose gel followed by purification with the Gel DNA Extraction Kit (Tiangen, Beijing, China). The target fragments were ligated to the PMD20-T vector (Takara, Dalian, China) and were transformed into E. coli strain JM109 competent cells. Seven or eight positive clones were further selected, and plasmids were isolated for insert sequencing as described above.

2.3. Sequence Data Analysis

After editing of the sequences, the noncoding regions were truncated, and the predicted coding region was aligned with the corresponding region of the cattle MC1R wild-type allele (^E+ ) (GenBank accession numbers Y13957 and Y19103; [23, 24]). The relationship among the yak alleles, defined according to their haplotypes, was evaluated using program Network 4.1 (http://www.fluxus-engineering.com/sharenet.htm), along with the ^E+ , ^ED , and e cattle MC1R alleles. The putative transmembrane structure of the MC1R protein, conducted by the TMpred program [25], was drawn to show the position of the amino acid substitutions.

In silico functional analysis of missense mutations was obtained using PANTHER [26]. PANTHER estimates the likelihood of a particular amino acid substitution to cause a functional impact on the protein. It calculates the substitution position-specific evolutionary conservation (subPSEC) score based on an alignment of evolutionarily related proteins [26-28]. The probability that a given variant will cause a deleterious effect on protein function is estimated by _Pdeleterious , such that a subPSEC score of -3 corresponds to a _Pdeleterious of 0.5 [29]. The subPSEC score is the negative logarithm of the probability ratio of the wild-type and mutant amino acids at a particular position. PANTHER subPSEC scores are continuous values from 0 (neutral) to about -10 (most likely to be deleterious).

3. Results and Discussion

Very limited information is available on coat color genetics in yak despite the presence of coat color variability among and, in some cases, within breeds. In addition, to our knowledge, before our study, no information about any coat color gene was available in yak. Therefore, in order to add data in this field and try to link coat color phenotype variations in this species with polymorphisms at the DNA level, we first characterized by sequencing the MC1R gene, whose mutations in cattle and several other species are responsible for eumelanism and pheomelanism.

The predicted coding region of the yak MC1R gene resulted of 954 bp, the same to that of the wild type cattle sequence [5-7], with identity > 99%. Comparing the 44 yak sequences, a total of 5 single nucleotide polymorphisms (SNPs) were identified. These SNPs defined three haplotypes that were indicated as ^EY1 , ^EY2 , and ^EY3 . Their sequences were submitted to GenBank database with the accession numbers FJ624478, FJ624479, and FJ624480, respectively. Considering as reference the wild type cattle sequence, the ^EY1 haplotype was characterized by the c.340C>A and c.871G>A substitutions, whereas the ^EY3 haplotype included the c.375C>T and c.663T>C mutations (Figure 2). The ^EY2 haplotype differed from the ^EY3 haplotype by c.642C>G substitution (Figure 2). None of the mutation events reported in cattle was found in the sequenced yaks (Figure 2 and data not shown). It is worth to note that in river buffalo (Bubalus bubalis ) Klungland et al. [24] identified the MC1R e allele with quite high frequency (0.20), suggesting a recurrent origin of the deletion causing this allele or a mutation event that occurred before the radiation of the Bovinae subfamily.

Figure 2: Sequence variations of the three MC1R haplotypes in yaks. Dots (*) denote identity with the cattle MC1R allele ^E+ , which was used as a reference sequence (GenBank accession numbers Y13957 and Y19103). Amino acid substitutions (AA Subst.) were listed below the nucleotide information. The three known cattle MC1R alleles were also included to demonstrate the difference between the yak and cattle sequences. Numbering was scored relative to the first codon of the cattle MC1R allele ^E+ .

[figure omitted; refer to PDF]

The relationship of the three major cattle MC1R alleles with the yak haplotypes is shown in Figure 3. The ^EY1 yak haplotype has four- and five-mutation distance with the ^EY2 and ^EY3 haplotypes, respectively, passing through the differentiation node constituted by the wild-type cattle sequence. Therefore, it could be expected that another haplotype, having the same sequence of the ^E+ cattle allele, might be present (or was present) in some other yak populations. This network distribution of the yak haplotypes resembles the mtDNA data that suggested that the Chinese yak was domesticated from two distinct matrilineal sources or from a heterogeneous pool containing two divergent lineages, with occasional gene introgression from cattle [30].

Figure 3: Network profile of the yak and cattle MC1R haplotypes. The links are labeled by nucleotide positions to designate transition; transversions are further specified in bold. The MC1R haplotype frequency in yak was proportional to the circle area. Yaks with white coat color and black coat color were marked by white and black, respectively. The three cattle alleles (^E+ ,^ED , and e ) were denoted by dots.

[figure omitted; refer to PDF]

Of the five identified yak SNPs, three (c.375C>T, c.642C>G, and c.663T>C) were predicted as synonymous mutations while the two polymorphisms characterizing the E^Y1 haplotype (c.340C>A and c.871G>A) were deduced as nonsynonymous substitutions causing a p.Q114K and a p.A291T amino acid change, respectively (Figures 2 and 4). These amino acid changes resided in the first extracellular loop region (p.Q114K) and in the seventh transmembrane domain (p.A291T) (Figure 4). In these regions, several amino acid substitutions, a few of which with effects on pigmentation, have been detected in other species, but none of the positions that are polymorphic in yak has been associated with any coat color variation. However, it is interesting to note that in the same transmembrane domain harboring the yak p.A291T amino acid change, a substitution at position 298 (p.Y298C) in black bear is associated with a white coat color variant [31]. Moreover, at position 294, a p.D294H substitution is associated with red hair and fair skin in human [32]. In silico functional analysis of the p.Q114K substitution did not support a putative functional role of this mutation (subSPEC=-2.60994; _Pdeleterious =0.4037) while a functional effects was more likely for the p.A291T amino acid change (subSPEC=-4.64874; _Pdeleterious =0.83872) . Therefore, it could be important to confirm the putative predicted functional effect of the latter mutation with specific experiments evaluating the pharmacological activity of this variant.

Figure 4: Hypothetical structure of the yak MC1R protein with the two amino acid changes identified in this study. Amino acid changes are scored relative to the ^EY1 haplotype. Transmembrane domains were predicted using the TMpred program [25], with reference to the reported cattle MC1R structure.

[figure omitted; refer to PDF]

The distribution of the MC1R haplotypes among the three investigated breeds is shown in Table 1. In Tianzhu yaks and Maiwa yaks, ^EY1 haplotype was predominant (0.68 and 0.75, respectively), followed by haplotype ^EY3 in Tianzhu animals (0.21) or by haplotype ^EY2 in Maiwa breed (0.22), whereas in Jiulong there was a relatively equal distribution of the three haplotypes. All six possible genotypes (derived by the combination of the three haplotypes) were identified in Jiulong with relatively even frequency distribution, while five and three genotypes were found in Tianzhu and Maiwa yaks, respectively. The ^EY2 /^EY3 genotype, which was absent from Tianzhu yaks and Maiwa yaks, had the highest frequency in Jiulong yaks (0.43). When we tentatively grouped the yak samples according to coat color, white yaks had a higher frequency of the haplotype ^EY1 (0.68 versus 0.53) and a lower frequency of haplotype ^EY2 (0.11 versus 0.28) than black yaks. From these data it seems evident that there is no association between across-breed coat color phenotype and any haplotype or haplotype combination (genotype). However, the presence of possible regulatory mutation(s) in noncoding regions, affecting MC1R gene expression and, ultimately, yak coat colors (black versus white), was not directly excluded in this study, as we investigated only the coding region and partial ^5' - and ^3' - untranslated regions. Nonetheless, to account for possible regulatory mutation(s) we should hypothesize its/their occurrence only in all three haplotypes identified in one coat color type but not in the same three haplotypes observed in the other coat color type, making this possibility improbable. Therefore, based on the available data, the Extension locus in yak might not play an active role in black versus white coat color variation across breeds, at least considering the three investigated breeds.

Following these results, it remains to be clarified the nature of the white coat color in Tianzhu yaks as at present there is no study that could clarify if this phenotype is due to pheomelanism or absence of melanocytes or other mechanisms. Detailed hair and skin phenotypical evaluations and histological analyses as well as cross-breeding experiments might help in evaluating this issue.

Coat color can be considered a complex trait for which multiple genes, pathways, or networks are involved. Therefore, many genetic factors, other than the MC1R variants identified in the present study, might be more directly linked to coat color phenotypes in yak. More studies are essential to draw a complete picture of coat color determination in this species.

Other information could be extracted from the reported results. A larger level of variability seems present in the Jiulong yaks compared to the other sampled breeds. Moreover, the Maiwa yak and Jiulong yak, all with black coat color phenotype, had a significantly different distribution (P<.05 , chi square test) of the ^EY1 haplotype. One possible explanation for this pattern is that the Jiulong yaks might have a different population structure and/or origin even if mtDNA data did not support any conclusion in this direction [30]. Further studies might be needed to confirm the differences of nuclear variability and, eventually, clarify the reasons of these results.

Acknowledgments

The authors thank Dr. Songjia Lai for his help in collecting samples and the two anonymous reviewers for helpful comments on the early version of the manuscript. This work was financially supported by the Chinese national key technologies R & D program (Grant no. 2008BADB2B01), the project of science and technology bureau of Sichuan Province (Grant no. 2008HH0017), and the National 863 Project of China (no. 2008AA101001). Y.-G.Yao. was supported by the Hundreds-Talent Program of the Chinese Academy of Sciences.

References

[1] A. G. Searle Comparative Genetics of Coat Colour in Mammals , Logos Press, London, UK, 1968.

[2] I. J. Jackson, "Molecular and developmental genetics of mouse coat color," Annual Review of Genetics , vol. 28, no. 1, pp. 189-217, 1994.

[3] I. J. Jackson, "Homologous pigmentation mutations in human, mouse and other model organisms," Human Molecular Genetics , vol. 6, no. 10, pp. 1613-1624, 1997.

[4] L. S. Robbins, J. H. Nadeau, K. R. Johnson, "Pigmentation phenotypes of variant extension locus alleles result from point mutations that alter MSH receptor function," Cell , vol. 72, no. 6, pp. 827-834, 1993.

[5] H. Klungland, D. I. Vage, L. Gomez-Raya, S. Adalsteinsson, S. Lien, "The role of melanocyte-stimulating hormone (MSH ) receptor in bovine coat color determination," Mammalian Genome , vol. 6, no. 9, pp. 636-639, 1995.

[6] H. Joerg, H. R. Fries, E. Meijerink, G. F. Stranzinger, "Red coat color in Holstein cattle is associated with a deletion in the MSHR gene," Mammalian Genome , vol. 7, no. 4, pp. 317-318, 1996.

[7] F. Rouzaud, J. Martin, P. F. Gallet, "A first genotyping assay of French cattle breeds based on a new allele of the extension gene encoding the melanocortin-1 receptor (MC1R )," Genetics Selection Evolution , vol. 32, no. 5, pp. 511-520, 2000.

[8] D. I. Våge, H. Klungland, D. Lu, R. D. Cone, "Molecular and pharmacological characterization of dominant black coat color in sheep," Mammalian Genome , vol. 10, no. 1, pp. 39-43, 1999.

[9] L. Fontanesi, F. Beretti, V. Riggio, "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 Genetics . In press

[10] L. Marklund, M. J. Moller, K. Sandberg, L. Andersson, "A missense mutation in the gene for melanocyte-stimulating hormone receptor (MC1R ) is associated with the chestnut coat color in horses," Mammalian Genome , vol. 7, no. 12, pp. 895-899, 1996.

[11] J. M. Kijas, R. Wales, A. Tornsten, "Melanocortin receptor 1 (MC1R ) mutations and coat color in pigs," Genetics , vol. 150, no. 3, pp. 1177-1185, 1998.

[12] J. M. H. Kijas, M. Moller, G. Plastow, L. Andersson, "A frameshift mutation in MC1R and a high frequency of somatic reversions cause black spotting in pigs," Genetics , vol. 158, no. 2, pp. 779-785, 2001.

[13] L. Fontanesi, M. Tazzoli, F. Beretti, V. Russo, "Mutations in the melanocortin 1 receptor (MC1R ) gene are associated with coat colours in the domestic rabbit (Oryctolagus cuniculus )," Animal Genetics , vol. 37, no. 5, pp. 489-493, 2006.

[14] S. Kerje, J. Lind, K. Schütz, P. Jensen, L. Andersson, "Melanocortin 1-receptor (MC1R ) mutations are associated with plumage colour in chicken," Animal Genetics , vol. 34, no. 4, pp. 241-248, 2003.

[15] L. A. Werth, G. A. Hawkins, A. Eggen, "Rapid communication: melanocyte stimulating hormone receptor (MC1R ) maps to bovine chromosome 18," Journal of Animal Science , vol. 74, no. 1, pp. 262, 1996.

[16] C. Maudet, P. Taberlet, "Holstein's milk detection in cheeses inferred from melanocortin receptor 1 (MC1R ) gene polymorphism," Journal of Dairy Science , vol. 85, no. 4, pp. 707-715, 2002.

[17] T. G. Berryere, S. M. Schmutz, R. J. Schimpf, C. M. Cowan, J. Potter, " TYRP1 is associated with dun coat colour in Dexter cattle or how now brown cow?," Animal Genetics , vol. 34, no. 3, pp. 169-175, 2003.

[18] V. Russo, L. Fontanesi, E. Scotti, M. Tazzoli, S. Dall'Olio, R. Davoli, "Analysis of melanocortin 1 receptor (MC1R ) gene polymorphisms in some cattle breeds: their usefulness and application for breed traceability and authentication of Parmigiano Reggiano cheese," Italian Journal of Animal Science , vol. 6, no. 3, pp. 257-272, 2007.

[19] D. Graphodatskaya, H. Joerg, G. Stranzinger, "Molecular and pharmacological characterisation of the MSH-R alleles in Swiss cattle breeds," Journal of Receptor and Signal Transduction Research , vol. 22, no. 1-4, pp. 421-430, 2002.

[20] H. Qiu, Z. R. Qing, Y. C. Chen, D. A. Wang Bovine Breeds in China , Shanghai Scientific & Technical Publishers, Shanghai, China, 1988.

[21] T. A. Olson, R. Fries, A. Ruvinsky, "Genetics of colour variation," The Genetics of Cattle , pp. 33-53, CAB International, Cambridge, UK, 1999.

[22] Y. K. Chen, H. Z. Huang, B. Liu, "Study on color inheritance of yaks," China Yak , vol. 2, no. 2, pp. 33-36, 1994.

[23] B. Kriegesmann, B. Dierkes, T. Leeb, S. Jansen, B. Brenig, "Two breed-specific bovine MC1-R alleles in Brown Swiss and Saler breeds," Journal of Dairy Science , vol. 84, no. 7, pp. 1768-1771, 2001.

[24] H. Klungland, H. G. Olsen, M. S. Hassanane, K. Mahrous, D. I. Våge, "Coat colour genes in diversity studies," Journal of Animal Breeding and Genetics , vol. 117, no. 4, pp. 217-224, 2000., [email protected]

[25] K. Hofmann, W. Stoffel, "TMbase--a database of membrane spanning proteins segments," Biological Chemistry Hoppe-Seyler , vol. 374, pp. 166-168, 1993.

[26] P. D. Thomas, M. J. Campbell, A. Kejariwal, "PANTHER: a library of protein families and subfamilies indexed by function," Genome Research , vol. 13, no. 9, pp. 2129-2141, 2003.

[27] P. D. Thomas, A. Kejariwal, "Coding single-nucleotide polymorphisms associated with complex vs. Mendelian disease: evolutionary evidence for differences in molecular effects," Proceedings of the National Academy of Sciences of the United States of America , vol. 101, no. 43, pp. 15398-15403, 2004., [email protected]

[28] P. D. Thomas, A. Kejariwal, N. Guo, "Applications for protein sequence-function evolution data: mRNA/protein expression analysis and coding SNP scoring tools," Nucleic Acids Research , vol. 34, web server issue, pp. W645-W650, 2006., [email protected]

[29] L. R. Brunham, R. R. Singaraja, T. D. Pape, A. Kejariwal, P. D. Thomas, M. R. Hayden, "Accurate prediction of the functional significance of single nucleotide polymorphisms and mutations in the ABCA1 gene," PLoS Genetics , vol. 1, no. 6, article e83, pp. 739-747, 2005.

[30] S.-J. Lai, S.-Y. Chen, Y.-P. Liu, Y.-G. Yao, "Mitochondrial DNA sequence diversity and origin of Chinese domestic yak," Animal Genetics , vol. 38, no. 1, pp. 77-80, 2007.

[31] K. Ritland, C. Newton, H. D. Marshall, "Inheritance and population structure of the white-phased "Kermode" black bear," Current Biology , vol. 11, no. 18, pp. 1468-1472, 2001.

[32] P. Valverde, E. Healy, I. Jackson, J. L. Rees, A. J. Thody, "Variants of the melanocyte-stimulating hormone receptor gene are associated with red hair and fair skin in humans," Nature Genetics , vol. 11, no. 3, pp. 328-330, 1995.