徐桂利，張文香，段玲欣，鞏元芳*，葛慕湘，劉謝榮，王書朋，果新苓，劉錚鑄*

(1.河北科技師范學(xué)院動(dòng)物科技學(xué)院，秦皇島 066004； 2.河北東光縣農(nóng)業(yè)局，東光 061600)

?

狐貍MC1R基因編碼區(qū)c.40A>C和c.41C>T相鄰變異研究

徐桂利1，張文香1，段玲欣1，鞏元芳1*，葛慕湘1，劉謝榮1，王書朋2，果新苓2，劉錚鑄1*

(1.河北科技師范學(xué)院動(dòng)物科技學(xué)院，秦皇島 066004； 2.河北東光縣農(nóng)業(yè)局，東光 061600)

為了檢測(cè)狐貍MC1R基因多態(tài)性及其與毛色表型的相關(guān)性，本研究采集兩個(gè)狐屬12種毛色共計(jì)163只狐貍的皮膚組織樣，利用PCR擴(kuò)增和產(chǎn)物直接測(cè)序的方法獲得狐貍MC1R基因1 054 bp長的核苷酸序列，并進(jìn)行了SNPs篩查。用PopGen32和SHEsis軟件對(duì)突變位點(diǎn)進(jìn)行了群體遺傳學(xué)分析，用PANTHER軟件評(píng)估了突變對(duì)基因產(chǎn)生的功能影響，用SPSS二元變量相關(guān)統(tǒng)計(jì)方法分析了多態(tài)位點(diǎn)與毛色表型間的相關(guān)性。結(jié)果表明，狐貍MC1R基因編碼區(qū)40(c.40A>C)和41 位點(diǎn)(c.41C>T)存在2 個(gè)相鄰錯(cuò)義突變，導(dǎo)致其編碼的第14位氨基酸發(fā)生了變異：當(dāng)?shù)?0 位點(diǎn)為A時(shí)，氨基酸由蘇氨酸(Thr)轉(zhuǎn)變?yōu)楫惲涟彼?Ile)；當(dāng)?shù)?0 位點(diǎn)為C時(shí)，氨基酸由脯氨酸(Pro)轉(zhuǎn)變?yōu)榱涟彼?Leu)。北極狐屬狐在41 位點(diǎn)基因型全部為TT型，而狐屬狐大部分個(gè)體均為CC型，不存在TT型，推測(cè)該位點(diǎn)可能是區(qū)分狐貍屬間的一個(gè)重要功能位點(diǎn)。PANTHER預(yù)測(cè)獲知第41 位點(diǎn)突變導(dǎo)致的氨基酸替換(p.Pro14Leu)對(duì)MC1R功能有顯著影響。SPSS二元變量相關(guān)分析結(jié)果表明，41 位點(diǎn)多態(tài)性與狐貍毛色表型存在顯著低度相關(guān)性，推測(cè)狐貍MC1R基因編碼區(qū)第41位點(diǎn)可能是參與其毛色形成的一個(gè)相對(duì)重要功能位點(diǎn)。

狐貍；MC1R基因；40和41位點(diǎn)；多態(tài)性；毛色

毛色是動(dòng)物表型的一種，對(duì)于狐貍來說毛色是衡量其毛皮經(jīng)濟(jì)價(jià)值的一個(gè)重要指標(biāo)。動(dòng)物的毛色是由黑色素基因控制的。黑色素主要包括兩大類：一類為真黑色素，另一類為褐黑色素，隨著黑色素的種類和含量不同，動(dòng)物毛發(fā)的顏色也會(huì)發(fā)生相應(yīng)變化，真黑色素和褐黑色素所含比例大小決定了毛色的深淺程度[1]。目前，世界上人工飼養(yǎng)的狐貍大約有40多種不同色型，分類學(xué)上歸屬于兩個(gè)屬：狐屬(Vulpes)和北極狐屬(Alopex)[2-3]。

有關(guān)狐貍黑素皮質(zhì)素受體1(Melanocortin 1 receptor，MC1R)基因多態(tài)性及其與毛色表型的相關(guān)性研究很少，目前僅有少量報(bào)道[4-5]。D.I.V?ge等[4]1997年報(bào)道，MC1R基因編碼區(qū)373位點(diǎn)的突變(c.373T>C)能引起狐屬狐銀黑色毛的產(chǎn)生，該課題組2005年又發(fā)現(xiàn)MC1R基因編碼區(qū)13和839兩個(gè)位點(diǎn)的錯(cuò)義突變(c.13G>T和c.839T>G)可阻滯藍(lán)色北極狐冬季藍(lán)色被毛的表達(dá)[5]。近年來，國內(nèi)外有關(guān)人類和哺乳動(dòng)物MC1R基因多態(tài)性與毛色表型的相關(guān)性研究較多，但大多數(shù)都是針對(duì)基因單個(gè)位點(diǎn)的突變與毛色表型的相關(guān)性研究，J.S.Palmer等[6]和E.E.Bashmakova等[7]均報(bào)道，人類MC1R基因編碼蛋白151、160和294(Arg151Cys、Arg160Trp和Asp294His)的變異與紅發(fā)有關(guān)。M.Abitbol等[8]報(bào)道驢MC1R基因編碼區(qū)629位點(diǎn)的一個(gè)單堿基突變(c.629T>C)與紅色毛有很強(qiáng)的相關(guān)性。李洪濤等[9]報(bào)道哈薩克羊MC1R基因218位點(diǎn)突變(T218A)對(duì)黑色毛有顯性作用。郭多等[10]報(bào)道家犬MC1R基因編碼第306位氨基酸的密碼子存在一個(gè)由CGA到TGA的終止突變，使MC1R基因的翻譯終止，導(dǎo)致黃色毛的形成。有關(guān)兩個(gè)相鄰位點(diǎn)的變異及其與動(dòng)物毛色表型的相關(guān)性研究目前鮮有報(bào)道。

鑒于以上研究背景，本研究擬采用PCR擴(kuò)增和產(chǎn)物直接測(cè)序的方法，分析MC1R基因的多態(tài)性及其與狐貍毛色表型的相關(guān)性。研究結(jié)果將為探明MC1R基因調(diào)控狐貍毛色機(jī)理以及在生產(chǎn)上培育出滿足人類需求的更多彩色狐貍奠定分子理論基礎(chǔ)。

1　材料與方法

1.1試驗(yàn)動(dòng)物

本研究所用163只成年狐貍來自河北省秦皇島市昌黎縣金島“狐、貉、貂育種場(chǎng)”，其中赤狐18只、銀黑狐28只、白銀狐14只、巧克力狐9只、黑理石狐6只、紅理石狐5只、銀十字狐11只、金十字狐2只、琥珀狐5只、日暉狐2只、藍(lán)色北極狐47只、白色北極狐16只，其中藍(lán)色北極狐和白色北極狐歸屬于北極狐屬，其余毛色狐貍均歸屬于狐屬。

樣品采集：11月底至12月初，從受試狐貍體側(cè)無菌手術(shù)采集1 cm2左右的皮膚組織塊，去毛，投至已滅菌的1.5 mL離心管中，然后將離心管迅速投入液氮中冷凍并置于-80 ℃冰箱保存，備用。

1.2基因組DNA的提取

基因組DNA的提取采用傳統(tǒng)的酚—氯仿抽提法[11]，TE溶解，-20 ℃凍存。

1.3引物設(shè)計(jì)和PCR擴(kuò)增

根據(jù)GenBank數(shù)據(jù)庫中登載的赤狐MC1R基因序列(序列號(hào)：X90844)，用Primer3(http：//www.genome.wi.mit.edu/cgibin/primer/primer3-www.cgi)設(shè)計(jì)擴(kuò)增產(chǎn)物包含整個(gè)編碼區(qū)、大小為1 055 bp的1對(duì)引物，其上游序列為：5′-GAACTG-AGCGAGACACCTGA-3′，下游序列：5′-ATCACC-ACCTCCCTTTGCCCA-3′。引物由上海生工生物工程技術(shù)服務(wù)有限公司合成。

PCR擴(kuò)增反應(yīng)體系為50 μL：基因組DNA 1 μL(75 ng·μL-1)，10×LA PCR Buffer Ⅱ(Mg2+Plus)5 μL，2.5 mmol·L-1dNTPs 8 μL，上游引物(20 μmol·L-1)1 μL，下游引物(20 μmol·L-1)1 μL，TaKaRa LATaq?聚合酶(5 U·μL-1)0.5 μL，滅菌dH2O 33.5 μL。

PCR 反應(yīng)條件：94 ℃預(yù)變性5 min；94 ℃變性45 s，64 ℃退火45 s，72 ℃延伸1 min，30個(gè)循環(huán)；72 ℃總延伸10 min。擴(kuò)增產(chǎn)物經(jīng)1.5%的瓊脂糖凝膠電泳并在凝膠自動(dòng)成像儀上檢測(cè)。

1.4PCR產(chǎn)物回收、測(cè)序

對(duì)于效果良好且量足夠的擴(kuò)增產(chǎn)物，回收純化后，由上海生工生物技術(shù)服務(wù)有限公司進(jìn)行測(cè)序。

1.5SNPs篩查及多態(tài)性分析

將163個(gè)不同毛色狐貍個(gè)體MC1R基因所得序列與網(wǎng)上登載的赤狐相應(yīng)序列(X90844)用BioEdit (v7.0.5.2)[12]和DNAMAN(5.2.2.0)軟件進(jìn)行比對(duì)及相似性分析。變異位點(diǎn)的等位基因、基因型以及單倍型及其頻率的分析，分別用PopGen32(version 1.31)[13]和SHEsis軟件[14]計(jì)算完成。

1.6功能預(yù)測(cè)

編碼區(qū)突變對(duì)基因產(chǎn)生的功能影響采用PANTHER軟件[15-16]進(jìn)行評(píng)估。該軟件根據(jù)已知功能蛋白質(zhì)家族多條序列的比對(duì)信息，分析蛋白質(zhì)家族中不同位置氨基酸發(fā)生變化時(shí)，對(duì)蛋白質(zhì)功能產(chǎn)生的影響程度?？赡墚a(chǎn)生的影響用特定位置取代進(jìn)化保守值(Substitution position-specific evolutionaryconservation score，subPSEC)來衡量，該指標(biāo)進(jìn)一步表述為異義替換SNP 對(duì)蛋白質(zhì)功能損害的概率(Pdeleterious)。subPSEC值或Pdeleterious越小(即越是趨向負(fù)值方向)表明越有可能對(duì)蛋白質(zhì)功能產(chǎn)生損害。

1.7統(tǒng)計(jì)分析

利用SPSS16.0軟件中二元變量相關(guān)分析法統(tǒng)計(jì)分析狐貍MC1R基因多態(tài)位點(diǎn)與毛色表型的相關(guān)性。

2　結(jié)　果

2.1PCR擴(kuò)增結(jié)果

利用所設(shè)計(jì)的引物，以其中1只供試銀黑狐基因組DNA為模板進(jìn)行PCR擴(kuò)增，結(jié)果見圖1。

1.PCR產(chǎn)物；M.DNA 相對(duì)分子質(zhì)量標(biāo)準(zhǔn)1.PCR product；M.DNA marker DL 2000圖1　狐貍MC1R基因的PCR擴(kuò)增Fig.1　Agarose gel electrophoretogram of fox MC1R PCR product (1 054 bp)

從圖1可以看出，PCR條帶接近于1 000 bp，且整齊、清晰、亮度也較好，初步說明擴(kuò)增產(chǎn)物就是預(yù)期的目的條帶。

2.2測(cè)序結(jié)果分析

將以上銀黑狐的擴(kuò)增產(chǎn)物回收純化后委托上海生工生物技術(shù)服務(wù)有限公司進(jìn)行雙向測(cè)序，測(cè)序結(jié)果經(jīng)BioEdit(v7.0.5.2)和DNAMAN(5.2.2.0)拼接整理，獲得1 054 bp長的一段核苷酸序列，該序列與GenBank上登載的赤狐MC1R基因相應(yīng)序列比對(duì)后，相似性高達(dá)99.53%，說明所得序列就是銀黑狐MC1R基因的序列。進(jìn)一步分析獲知該序列包括5′UTR區(qū)54 bp、編碼區(qū)954 bp和3′UTR區(qū)46 bp。

2.3SNPs篩查及分析

以上所得銀黑狐序列加上剩余162只狐貍的PCR產(chǎn)物測(cè)序后，用BioEdit(v7.0.5.2)軟件進(jìn)行比對(duì)分析，發(fā)現(xiàn)編碼區(qū)第40(c.40A>C)和41位點(diǎn)(c.41C>T)存在2個(gè)相鄰的錯(cuò)義突變，導(dǎo)致其編碼的第14位氨基酸的變異：當(dāng)?shù)?0 位點(diǎn)為A時(shí)，氨基酸由蘇氨酸(Thr)轉(zhuǎn)變?yōu)楫惲涟彼?Ile)，當(dāng)?shù)?0 位點(diǎn)為C時(shí)，氨基酸由脯氨酸(Pro)轉(zhuǎn)變?yōu)榱涟彼?Leu)。另外，對(duì)于第40位點(diǎn)，狐屬狐(包括赤狐、銀黑狐等10種毛色狐)存在3種基因型：AA、AC和CC，而北極狐屬狐(包括藍(lán)色北極狐和白色北極狐2種毛色狐)只有CC 一種基因型。對(duì)于第41位點(diǎn)，狐屬狐(包括赤狐、銀黑狐等10種毛色狐)存在2種基因型：CC和CT型，而北極狐屬狐(包括藍(lán)色北極狐和白色北極狐2種毛色狐)只有TT 一種基因型。兩位點(diǎn)基因型的分布和變異情況見表1和圖2。

2.4突變位點(diǎn)群體遺傳學(xué)分析

利用PopGen32(version 1.31)和SHEsis軟件對(duì)163只受試狐貍MC1R基因編碼區(qū)第40和41相鄰?fù)蛔兾稽c(diǎn)等位基因和基因型頻率進(jìn)行了計(jì)算(表2)，同時(shí)對(duì)兩位點(diǎn)進(jìn)行了單倍型的構(gòu)建(表3)。由表2可知，兩位點(diǎn)C等位基因的頻率都較高，分別為0.855和0.524。由表3可知，兩位點(diǎn)共構(gòu)建3種單倍型：H1(AC)、H2(CC)和H3(CT)，其中H2和H3的頻率較高，分別為0.380和0.476，為優(yōu)勢(shì)單倍型。

表1不同毛色狐貍MC1R基因編碼區(qū)第40和41位點(diǎn)的基因型變異

Table 1Identified genetic variations in 12 fox sub-populations with different coat color patterns

屬Genus不同毛色狐貍Differentcoatcolorfox變異Variationc.40A>Cc.41C>T狐屬Vulpes赤狐(18)RedfoxAA(3)AC(4)CC(11)CC(17)CT(1)銀黑狐(28)SilverfoxAC(10)CC(18)CC(19)CT(9)白銀狐(14)WhiteSilverfoxAA(3)AC(6)CC(5)CC(14)巧克力狐(9)ChocolatefoxCC(9)CC(4)CT(5)黑理石狐(6)BlackMarblefoxAA(2)AC(4)CC(6)紅理石狐(5)RedMarblefoxAA(2)AC(2)CC(1)CC(5)銀十字狐(11)BlackCrossfoxCC(11)CC(1)CT(10)金十字狐(2)GoldCrossfoxCC(2)CT(2)琥珀狐(5)AmberfoxCC(5)CC(2)CT(3)日暉狐(2)SunGlowfoxAA(1)CC(1)CC(2)北極狐屬Alopex藍(lán)色北極狐(47)ArcticBluefoxCC(47)TT(47)白色北極狐(16)ArcticWhitefoxCC(16)TT(16)

表2狐貍MC1R基因編碼區(qū)第40和41相鄰?fù)蛔兾稽c(diǎn)的等位基因和基因型頻率

Table 2Frequencies of genotypes and alleles at 40 and 41 sites in coding region of foxMC1Rgene

突變位點(diǎn)Mutationsite基因型頻率Genotypefrequency等位基因頻率Allelefrequency基因型Geno.頻率Freq.基因型Geno.頻率Freq.基因型Geno.頻率Freq.等位基因Allele頻率Freq.等位基因Allele頻率Freq.c.40A>CAA0.066AC0.157CC0.777A0.145C0.855c.41C>TCC0.428CT0.190TT0.379C0.524T0.476

表3狐貍MC1R基因編碼區(qū)第40和41相鄰?fù)蛔兾稽c(diǎn)單倍型頻率

Table 3Haplotype frequencies ofMC1Rgene covering the 2 polymorphic sites in foxes

單倍型Haplotype頻率Frequency突變位點(diǎn)Mutationsitec.40A>Cc.41C>TH10.145ACH20.380CCH30.476CT

2.5編碼區(qū)40和41相鄰?fù)蛔兾稽c(diǎn)對(duì)MC1R功能的影響

為了推測(cè)MC1R基因編碼區(qū)第40和41位點(diǎn)突變對(duì)蛋白功能的影響，用PANTHER軟件對(duì)比了與其進(jìn)化相關(guān)的蛋白質(zhì)中特定位點(diǎn)發(fā)生突變時(shí)對(duì)蛋白質(zhì)功能的影響程度，由此得到功能損傷的似然值(表4)。subPSEC值是特定位置的野生型和突變型氨基酸的概率比值的負(fù)對(duì)數(shù)，該值在0(中性的)～-10(極可能對(duì)蛋白質(zhì)功能有害)之間連續(xù)變化[17]。

圖2　狐貍MC1R基因編碼區(qū)40和41位點(diǎn)的突變Fig.2　The mutation of 40 and 41 sites in coding region of fox MC1R gene

表4MC1R第14位氨基酸變異的subPSEC和Pdeleterious值

Table 4The subPSEC andPdeleteriousscores of the 14th amino acid substitution

核苷酸變異Nucleotidevariation40位點(diǎn)41位點(diǎn)氨基酸變異AminoacidvariationssubPSEC值subPSECscoresPdeleterious值Pdeleteriousscores40位點(diǎn)為A時(shí)c.41C>Tp.Thr14Ile-1.625940.2019740位點(diǎn)為C時(shí)c.41C>Tp.Pro14Leu-4.204460.76932

2.6統(tǒng)計(jì)分析結(jié)果

通過SPSS二元變量相關(guān)分析得出狐貍MC1R基因編碼區(qū)第40位點(diǎn)基因型與毛色表型之間的相關(guān)系數(shù)(r)為-0.050，二者不具相關(guān)性，P=0.455，無統(tǒng)計(jì)學(xué)意義。第41位點(diǎn)基因型與毛色表型之間的相關(guān)系數(shù)(r)為-0.359，二者具低度相關(guān)，P= 0.000，說明二者有顯著相關(guān)性(表5)。

表5MC1R基因第40和41 位點(diǎn)變異與毛色表型的相關(guān)分析結(jié)果

Table 5The correlation analysis between genotype at 40 and 41 sites and coat color

位點(diǎn)Site相關(guān)系數(shù)(r)CorrelationcoefficientP值Pvalue結(jié)果Result40-0.0500.455不相關(guān),無統(tǒng)計(jì)學(xué)意義41-0.359**0.000低度相關(guān)

|r|<0.3為不相關(guān)；0.3≤|r|<0.5為低度相關(guān)；0.5≤|r|<0.8為中度相關(guān)；0.8≤|r|<1為高度相關(guān)。P<0.05時(shí)，兩個(gè)變量間相關(guān)性顯著；當(dāng)P<0.01時(shí)，兩個(gè)變量間相關(guān)性非常顯著，P>0.05 時(shí)，兩個(gè)變量間沒有顯著的相關(guān)性[18]

|r|<0.3 indicate no correlation；0.3≤|r|<0.5 indicate low correlation；0.5≤|r|<0.8 indicate moderate correlation；0.8≤|r|<1 indicate high correlation.P<0.05 indicate significant correlation；P<0.01 indicate extremely significant correlation；P>0.05 indicate no significant correlation

3　討　論

近年來，隨著人們對(duì)化工染料污染的擔(dān)憂和崇尚自然風(fēng)潮的興起，綠色、天然、環(huán)保、優(yōu)質(zhì)毛皮動(dòng)物產(chǎn)品已成為廣大生產(chǎn)者與消費(fèi)者的追求目標(biāo)。因此，創(chuàng)制并選育天然彩色毛皮動(dòng)物新品種，將成為現(xiàn)代畜牧業(yè)發(fā)展的緊迫任務(wù)，而探明毛皮動(dòng)物彩色絨毛分子遺傳機(jī)制是實(shí)現(xiàn)上述目標(biāo)的重要理論基礎(chǔ)。狐貍是一種重要的毛皮動(dòng)物，其與家犬均歸屬于犬科，有關(guān)家犬MC1R基因的相關(guān)研究已有大量報(bào)道，并已證實(shí)該基因與其毛色有關(guān)[19-21]，但有關(guān)狐貍MC1R基因的相關(guān)研究則相對(duì)較少。J.M.Newton等[19-20]發(fā)現(xiàn)家犬MC1R蛋白的90和306位點(diǎn)變異(S90G和R306ter)分別與其黑/棕和紅/黃毛色相關(guān)。D.L.Dreger等[21]發(fā)現(xiàn)家犬MC1R蛋白的78位點(diǎn)變異(p.Gly78Val)與Saluki和Afghan Hound家犬的毛色有關(guān)。以上研究都是針對(duì)單堿基突變與毛色的相關(guān)性研究。本研究分析了狐屬和北極狐屬12種毛色共計(jì)163只狐貍MC1R基因編碼區(qū)序列，在第40和41位點(diǎn)發(fā)現(xiàn)了2個(gè)相鄰的錯(cuò)義突變：c.40A>C和c.41C>T。當(dāng)40位點(diǎn)為A時(shí)，導(dǎo)致其編碼的第14位氨基酸發(fā)生蘇氨酸向異亮氨酸的轉(zhuǎn)變(p.Thr14Ile)；當(dāng)40位點(diǎn)為C時(shí)，導(dǎo)致第14位氨基酸發(fā)生脯氨酸向亮氨酸的轉(zhuǎn)變(p.Pro14Leu)，這一結(jié)果與J.Nowacka-Woszuk等[22]的研究結(jié)果相一致。J.Nowacka-Woszuk等[22]在分析家犬、赤狐、北極狐和貉MC1R基因編碼區(qū)的變異時(shí)，也發(fā)現(xiàn)了赤狐編碼區(qū)第40和41位點(diǎn)的變異，但由于其缺乏表型記錄沒有分析兩位點(diǎn)與毛色表型的相關(guān)性。本研究中，北極狐屬狐MC1R基因編碼區(qū)第41位點(diǎn)只有TT一種基因型，而狐屬狐則有CC和CT 2種基因型，不存在TT基因型，推測(cè)該位點(diǎn)可能是區(qū)分狐貍屬間的一個(gè)重要功能位點(diǎn)。群體遺傳分析表明，對(duì)于40位點(diǎn)，C為優(yōu)勢(shì)等位基因(0.855)，單倍型H2(CC，0.380)和H3(CT，0.476)為優(yōu)勢(shì)單倍型，說明大部分狐貍40位點(diǎn)均為C等位基因。通過PANTHER軟件功能預(yù)測(cè)，當(dāng)40位點(diǎn)為C時(shí)，p.Pro14Leu處的subPSEC和Pdeleterious值分別是-4.204 46和0.769 32，表明該處氨基酸替換對(duì)MC1R功能有重要影響。SPSS二元變量統(tǒng)計(jì)分析結(jié)果表明，狐貍MC1R基因編碼區(qū)第41位點(diǎn)多態(tài)性與毛色表型存在顯著低度相關(guān)性，推測(cè)第41位點(diǎn)可能是參與毛色形成的一個(gè)相對(duì)重要功能位點(diǎn)。

4　結(jié)　論

本研究采用PCR擴(kuò)增和產(chǎn)物直接測(cè)序的方法獲得了狐屬和北極狐屬共計(jì)163只不同毛色狐貍MC1R基因的編碼區(qū)序列，在第40和41位點(diǎn)發(fā)現(xiàn)2個(gè)相鄰錯(cuò)義突變(c.40A>C和c.41C>T)。群體遺傳分析表明，狐貍MC1R基因編碼區(qū)41位點(diǎn)可能是區(qū)分狐貍屬間的一個(gè)重要功能位點(diǎn)。PANTHER軟件功能預(yù)測(cè)41位點(diǎn)突變(c.41C>T)導(dǎo)致的氨基酸替換(p.Pro14Leu)對(duì)MC1R功能有重要影響。SPSS二元變量統(tǒng)計(jì)分析結(jié)果進(jìn)一步表明，41位點(diǎn)多態(tài)性與狐貍毛色表型存在顯著相關(guān)性。

[1]BENNETT D C，LAMOREUX M L.The color loci of mice-a genet iccentury[J].PigmentCellRes，2003，16(4)：333-344.

[2]馬澤芳，崔凱.貂狐貉實(shí)用養(yǎng)殖技術(shù)[M].北京：中國農(nóng)業(yè)出版社，2014.

MA Z F，CUI K.Practical breeding techniques of mink，fox and raccoon dog[M].Beijing：China Agriculture Press，2014.(in Chinese)

[3]佟煜仁，張志明.毛皮動(dòng)物毛色遺傳及繁育新技術(shù)[M].北京：金盾出版社，2009.

TONG Y R，ZHANG Z M.New technology for coat color genetics and breeding of fur animal[M].Beijing：Jindun Press，2009.(in Chinese)

[4]V?GE D I，LU D，KLUNGLAND H，et al.A non-epistatic interaction of agouti and extension in the fox，Vulpesvulpes[J].NatGenet，1997，15(3)：311-315.

[5]V?GE D I，F(xiàn)UGLEI E，SNIPSTAD K，et al.Two cysteine substitutions in the MC1R generate the blue variant of the Arctic fox (Alopexlagopus) and prevent expression of the white winter coat[J].Peptides，2005，26(10)：1814-1817.

[6]PALMER J S，DUFFY D L，BOX N F，et al.Melanocortin-1 receptor polymorphisms and risk of melanoma：is the association explained solely by pigmentation phenotype?[J].AmJHumGenet，2000，66(1)：176-186.

[7]BASHMAKOVA E E，KRASITSKAYA V V，BONDAR A A，et al.Bioluminescent assay to detect melanocortin-1 receptor (MC1R) polymorphisms (R160W，R151C，and D294H)[J].MolBiol(Mosk)，2015，49(6)：953-958.

[8]ABITBOL M，LEGRAND R，TIRET L.A missense mutation in melanocortin 1 receptor is associated with the red coat colour in donkeys[J].AnimGenet，2014，45(6)：878-880.

[9]李洪濤，曾獻(xiàn)存，張文祥，等.哈薩克綿羊MC1R和ASIP基因多態(tài)性及表達(dá)量與被毛顏色表型相關(guān)性的研究[J].畜牧獸醫(yī)學(xué)報(bào)，2013，44(3)：366-375.

LI H T，ZENG X C，ZHANG W X，et al.Analysis of the relationship between polymorphism/mRNA expression levels ofMC1RandASIPand coat color phenotype in Kazakh sheep[J].ActaVeterinariaetZootechnicaSinica，2013，44(3)：366-375.(in Chinese)

[10]郭多，巴彩鳳，蘇玉虹，等.犬MC1R基因R306ter與毛色性狀相關(guān)性研究[J].中國實(shí)驗(yàn)動(dòng)物學(xué)報(bào)，2004，12(2)：100-102.

GUO D，BA C F，SU Y H，et al.Study on the relationship between r306ter loci in mc1r gene and coat color trait in dogs[J].ActaLaboratoriumAnimalisScientiaSinica，2004，12(2)：100-102.(in Chinese)

[11]SAMBROOK J，F(xiàn)RITSCH E F，MANIATIS T.Molecular cloning[M].New York：Cold Spring Habor Laboratory Press，1989.

[12]HALL T A.BioEdit：A user-friendly biological sequences alignment editor and analysis program for Windows 95/98/NT[J].NucleicAcidsSympSer，1999，41：95-98.

[13]YEH F C，YANG R C，BOYLE T.POPGENE version 1.31：Microsoft Windows-based freeware for population genetic analysis-Quick User Guide[M].University of Alberta and Centre for International Forestry Research.Edmonton，AB Canada，1999：1-28.

[14]SHI Y Y，HE L.SHEsis，a powerful software platform for analyses of linkage disequilibrium，haplotype construction，and genetic association at polymorphism loci[J].CellRes，2005，15(2)：97-98.

[15]BRUNHAM L R，SINGARAJA R R，PAPE T D，et al.Accurate prediction of the functional significance of single nucleotide polymorphisms and mutations in the ABCA1 gene[J].PLoSGenet，2005，1(6)：e83.

[16]THOMAS P D，KEJARIWAL A，GUO N，et al.Applications for protein sequence-function evolution data：mRNA/protein expression analysis and coding SNP scoring tools[J].NucleicAcidsRes，34 (Web Server issue)：W645-650.

[17]苗永旺，吳桂生，王磊，等.黑色素皮質(zhì)素受體1 基因?qū)λＣz傳的影響[J].中國科學(xué)C輯：生命科學(xué)，2009，39(12)：1155-1161.

MIAO Y W，WU G S，WANG L，et al.The role of MC1R gene in buffalo coat color[J].ScienceinChina(SeriesC：LifeSciences)，2009，39(12)：1155-1161.(in Chinese)

[18]孫逸敏.利用SPSS軟件分析變量間的相關(guān)性[J].新疆教育學(xué)院學(xué)報(bào)，2007，23(2)：120-123.

SUN Y M.Using SPSS software to analyze the correlation between variables[J].JournalofXinjiangEducationInstitute，2007，23(2)：120-123.(in Chinese)

[19]NEWTON J M，WILKIE A L，HE L，et al.Melanocortin 1 receptor variation in the domestic dog[J].MammGenome，2000，11(1)：24-30.

[20]WANG G D，CHENG L G，F(xiàn)AN R X，et al.Signature of balancing selection at the MC1R gene in Kunming dog populations[J].PLoSOne，2013，8(2)：e55469.

[21]DREGER D L，SCHMUTZ S M.A new mutation in MC1R explains a coat color phenotype in 2 "old" breeds：Saluki and Afghan hound[J].JHered，2010，101(5)：644-649.

[22]NOWACKA-WOSZUK J，SALAMON S，GORNA A，et al.Missense polymorphisms in the MC1R gene of dog，red fox，arctic fox and Chinese raccon dog[J].JAnimBreedGenet，2013，130(2)：136-141.

(編輯郭云雁)

Study on Adjacent Variation of c.40A>C and c.41C>T of FoxMC1RGene in Different Populations

XU Gui-li1，ZHANG Wen-xiang1，DUAN Ling-xin1，GONG Yuan-fang1*，GE Mu-xiang1，LIU Xie-rong1，WANG Shu-peng2，GUO Xin-ling2，LIU Zheng-zhu1*

(1.CollegeofAnimalScienceandTechnology，HebeiNormalUniversityofScience&Technology，Qinhuangdao066004，China；2.AgricultureBureauofDongguangCountyinHebeiProvince，Dongguang061600，China)

In order to detect the relationship between polymorphism ofMC1Rgene and coat color in fox，a total of 163 skin samples of 12 coat color foxes were collected.The nucleotide sequence (1 054 bp) of foxMC1Rgene were obtained by the method of PCR and direct sequencing，and the polymorphism were analyzed.The population genetics were analyzed using PopGen32 and SHEsis softwares.The effect of mutations on the function ofMC1Rgene was evaluated using PANTHER software.The relationships between the variable sites and coat color were analyzed by the statistical methods of SPSS bivarate correlation analysis.Two adjacent missense mutations (c.40A>C and c.41C>T) were found in the coding region of foxMC1Rgene，which resulted in codon change of p.Thr14Ile or p.Pro14Leu.When 40 site was A，it led to the substitution between threonine (Thr) and isoleucine (Ile).When 40 site was C，it led to the substitution between proline (Pro) and leucine (Leu).All the genotypes of two coat color foxes belonging toAlopexwere TT.However，the genotypes of most coat color foxes belonging toVulpeswere CC.It was supposed that 41 site was important in distinguishingAlopexandVulpes.The in silico functional analysis showed that the amino acid substitution at p.Pro14Leu had significant impact on the function of MC1R.The statistical analysis showed the polymorphism of 41 site had significant low correlation with the fox coat color.The results indicate that SNP c.41C>T in the coding region of theMC1Rgene is probably associated with the coat color in fox.

fox；MC1Rgene；40 and 41 sites；polymorphism；coat color

10.11843/j.issn.0366-6964.2016.10.010

2016-02-24

國家自然科學(xué)基金(31272412)；河北省自然科學(xué)基金(C2013407101；C2016407114)；河北省高等學(xué)校創(chuàng)新團(tuán)隊(duì)領(lǐng)軍人才培育計(jì)劃(LJRC004)；河北科技師范學(xué)院科學(xué)研究基金(1)

徐桂利(1986-)，男，山東臨沂人，碩士生，主要從事動(dòng)物分子遺傳育種研究，E-mail：296881485@163.com

鞏元芳，教授，E-mail：gyfkeyan@163.com；劉錚鑄，教授，E-mail：liuzhengzhu@163.com

S829.9;S813.3

A

0366-6964(2016)10-2020-07