王歡歡，陳美玲，樓立峰，張 雷，張成先，陳賢惠，章學(xué)東*

(1. 杭州市農(nóng)業(yè)科學(xué)研究院，杭州 310024； 2. 浙江光大種禽業(yè)有限公司，海寧 310016)

雞GNAS基因啟動(dòng)子突變及其與膚色性狀的相關(guān)性

王歡歡1，陳美玲1，樓立峰1，張 雷1，張成先2，陳賢惠2，章學(xué)東1*

(1. 杭州市農(nóng)業(yè)科學(xué)研究院，杭州 310024； 2. 浙江光大種禽業(yè)有限公司，海寧 310016)

為研究雞GNAS基因突變情況及其對烏骨雞膚色的影響，本試驗(yàn)選擇東鄉(xiāng)綠殼蛋雞樣本，采用PCR-RFLP和DNA測序方法，篩查膚色相關(guān)的GNAS突變位點(diǎn)；采用熒光定量PCR方法，檢測EDN3、MC1R、GNAS、ADCY3基因的表達(dá)量；本試驗(yàn)還對3個(gè)雞種進(jìn)行了突變檢測。結(jié)果顯示，在家雞20號染色體的11 167 660堿基位置存在T/C突變，該T2270C突變點(diǎn)是GNAS基因啟動(dòng)子的關(guān)鍵位點(diǎn)。PCR產(chǎn)物經(jīng)DraI酶切電泳后，CC純合基因型出現(xiàn)一條342 bp的條帶；TC雜合基因型出現(xiàn)171和342 bp兩條帶；TT隱性純合基因型出現(xiàn)一條171 bp的條帶。3種基因型的膚色L值差異極顯著(P<0.01)，CC基因型的L值最小，表現(xiàn)為膚色最深；TT基因型的L值最大，表現(xiàn)為膚色淺白。EDN3基因的表達(dá)量以CC基因型最高，TT基因型最低，3種基因型間差異極顯著(P<0.01)。MC1R、GNAS和ADCY3基因的表達(dá)量以CC基因型最高，但基因型間的表達(dá)差異均不顯著(P>0.05)；GNAS與上游MC1R基因的表達(dá)量無顯著相關(guān)性，而與下游ADCY3基因的表達(dá)量極顯著相關(guān)(P<0.01)。江山烏骨雞、伊莎蛋雞、東鄉(xiāng)×伊莎測交組合雞的突變分型結(jié)果符合理論值(P>0.05)。本研究表明，雞GNAS基因啟動(dòng)子T2270C突變與膚色性狀強(qiáng)相關(guān)，突變檢測可輔助用于烏骨雞的膚色分型育種。

雞；GNAS基因；突變；膚色

鳥嘌呤核苷酸結(jié)合蛋白α亞基(Guanine nucleotide-binding protein Gs subunit alpha，GNAS)是位于細(xì)胞膜內(nèi)側(cè)的一類信號傳導(dǎo)蛋白，參與細(xì)胞的多種生命活動(dòng)[1]。在黑色素代謝通路中，GNAS與具有跨膜結(jié)構(gòu)的黑素皮質(zhì)素受體(MC1R)偶聯(lián)，激活腺苷酸環(huán)化酶(ADCY)，催化ATP生成第二信使cAMP，進(jìn)而影響胞內(nèi)黑色素的合成[2]。據(jù)家雞全基因組測序數(shù)據(jù)，雞GNAS基因位于20號染色體，全長為48 456 bp，含有10個(gè)外顯子。研究表明，烏骨雞真皮、內(nèi)臟等部位的黑色素沉積現(xiàn)象，又稱纖維化黑色變(Fibromelanosis，F(xiàn)m)，是由20號染色體上的Fm區(qū)域所決定的，烏骨雞在該區(qū)域存在著兩段拷貝數(shù)重復(fù)區(qū)(CNVs)，其中一段CNVs包含有與黑色素沉積密切相關(guān)的EDN3基因[3-4]，而GNAS等數(shù)個(gè)基因則處在兩段CNVs的中間序列。同時(shí)，醫(yī)學(xué)研究發(fā)現(xiàn)，人的GNAS基因突變導(dǎo)致黑色素瘤(Melanoma)、纖維骨性病變(Fibro-osseous)、纖維性發(fā)育不良(Fibrous dysplasia)等疾病的發(fā)生[5-7]。以上結(jié)果均提示，GNAS可能與黑色素性狀有關(guān)。

國內(nèi)多數(shù)烏骨雞種的群體內(nèi)膚色往往深淺不一，例如東鄉(xiāng)綠殼蛋雞的冠、喙、皮、肉、骨、脛、趾多呈烏黑色，但還有一定比例的白膚個(gè)體[8]。這種情況表明雞種內(nèi)存在烏膚主效基因的純合與雜合，會(huì)對后續(xù)的育種應(yīng)用造成影響。本研究目的在于篩選雞的GNAS基因突變，探究GNAS突變是否與雞的烏膚性狀相關(guān)，以期為尋找有效的分子標(biāo)記輔助育種方法提供依據(jù)。

1 材料與方法

1.1 突變檢測樣本

試驗(yàn)選用浙江光大種禽業(yè)有限公司保存的東鄉(xiāng)綠殼蛋雞150只、江山烏骨雞60只、伊莎蛋雞B系60只，以及東鄉(xiāng)綠殼蛋雞的測交后代67只。由專人飼養(yǎng)管理，按程序免疫。采集所有供試個(gè)體的翅靜脈血液，加入70%乙醇中，酚/氯仿法抽提基因組DNA，于-20℃保存。

1.2 PCR-RFLP分析

根據(jù)雞GNAS基因全序列(GenBank Accession No.：NC_006107.4)，采用Primer 6.0軟件設(shè)計(jì)啟動(dòng)子區(qū)引物(表1)，交由生工生物工程(上海)股份有限公司合成。

表1 突變檢測引物信息

Table 1 The information of primer for gene mutation detection

引物名稱Primer引物序列(5′?3′)Primersequence堿基位置/bpBaseposition產(chǎn)物長度/bpProductlengthGNAS＿promo3F：ACGACAGCAGGAAACGT11167489～11167505R：TGTGGCAGGCAGGTAACTGTGT11167809～11167830342

F. 上游引物；R.下游引物。表2同

F. Forward primer; R. Reverse primer. The same as Table 2

PCR反應(yīng)體系：DNA 20 ng，2× Taq PCR Master Mix 12.5 μL，上下游引物各10 μmol·L-1，加水補(bǔ)至 25 μL。上ABI 9700型PCR儀(ABI，美國)進(jìn)行擴(kuò)增，反應(yīng)程序：94 ℃ 5 min；然后進(jìn)入PCR循環(huán)，94 ℃ 30 s，57 ℃ 30 s，72 ℃ 30 s，35個(gè)循環(huán)；最后72 ℃延伸10 min。

將PCR產(chǎn)物用限制性內(nèi)切酶DraI進(jìn)行酶切，反應(yīng)體系：PCR產(chǎn)物5 μL，10×buffer 2 μL，酶1 μL，加水補(bǔ)至20 μL。37 ℃酶切1 h。用1.5%瓊脂糖凝膠電泳檢測酶切產(chǎn)物，并送杭州擎科梓熙生物技術(shù)有限公司進(jìn)行DNA測序。

1.3 膚色和表達(dá)量檢測樣本

于126日齡時(shí)，采用美能達(dá)CR400型色差計(jì)(Konica Minolta，日本)測定150只東鄉(xiāng)綠殼蛋雞的胸側(cè)膚色L、a、b值。L為明度值(0～100)，0表示黑色，100表示白色；a為紅綠光值，正值表示偏紅，負(fù)值表示偏綠；b為黃藍(lán)光值，正值標(biāo)示偏黃，負(fù)值表示偏藍(lán)。

同一日齡，分別選取各基因型的東鄉(xiāng)綠殼蛋雞各6只，屠宰取胸側(cè)皮膚樣，液氮保存。采用TRIzol Plus RNA純化試劑盒(Invitrogen，美國)提取總RNA，SuperScript Ⅲ反轉(zhuǎn)錄試劑盒(Invitrogen，美國)合成單鏈cDNA，于-20 ℃保存。

1.4 定量PCR

根據(jù)NCBI數(shù)據(jù)庫中相關(guān)基因的mRNA序列，采用Primer 6.0軟件分別設(shè)計(jì)EDN3、MC1R、GNAS、ADCY3 4個(gè)目的基因和18S內(nèi)參基因的mRNA引物(表2)，并交由生工生物工程(上海)股份有限公司合成。

表2 表達(dá)量檢測引物信息

Table 2 The information of primer for gene expression detection

基因名稱GeneGenBank序列號GenBankNo．引物序列(5′?3′)Primersequence產(chǎn)物長度/bpProductlength18SAF173612．1F：CCGGACACGGACAGGATTGACAR：CAGACAAATCGCTCCACCAACTAAG94EDN3XM＿015296553．1F：CATCTGGATCAACACCCCAGR：ACTTTGATTCTGGCCCGTAG97MC1RNM＿001031462．1F：CCTGCTCTGCCTCATTGGCTTCTR：GGTAGATGGTGGGCTGCTTCTG131GNASXM＿015296547．1F：GCACCTCCAGTTGAACTAGCCAATR：CTCCGGTGGGAAATCGAAGTCT96ADCY3XM＿015285193．1F：CTGACTTTGACGCACTCCTGGATGR：CCATGTAGGTGCTGCCGATGGT78

定量PCR反應(yīng)體系：cDNA 1 μL，Power SYBR Green Master Mix 10 μL，上、下游引物各0.5 μL，加水補(bǔ)至20 μL。上CFX384多重實(shí)時(shí)熒光定量PCR儀(Bio-Rad，美國)進(jìn)行擴(kuò)增，反應(yīng)程序：95℃ 1 min；然后進(jìn)入PCR循環(huán)，95℃ 15 s，63℃ 25 s，40個(gè)循環(huán)。55～95℃繪制熔解曲線，每0.2℃讀板一次。每個(gè)樣本做3個(gè)重復(fù)，記錄Ct值，2-△△Ct法計(jì)算目的基因的相對表達(dá)量。

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

人(Homosapiens)、豬(Susscrofa)、小鼠(Musmusculus)、爪蛙(Xenopustropicalis)、雞(Gallusgallus)、鵪鶉(Coturnixjaponica)、斑馬魚(Daniorerio)的GNAS基因序列來自NCBI數(shù)據(jù)庫，分別以轉(zhuǎn)錄起始位點(diǎn)為基準(zhǔn)，取其上游5 000 bp至下游1 000 bp共6 kb序列進(jìn)行生物信息學(xué)分析。采用MEGA 5.1軟件構(gòu)建進(jìn)化樹[9]；比對雞-鵪鶉-人、雞-鵪鶉-豬、雞-鵪鶉-小鼠、雞-鵪鶉-爪蛙、雞-鵪鶉-斑馬魚共5組序列，計(jì)算各組每500 bp的保守堿基率(保守堿基數(shù)/堿基總數(shù))。

以雞GNAS基因的轉(zhuǎn)錄起始位點(diǎn)為基準(zhǔn)，取其上游1～3 000 bp序列進(jìn)行生物信息學(xué)分析。采用Neural Network Promoter Prediction(http://www.fruitfly.org/seq_tools/promoter.html)和FFROM(http://linux1.softberry.com/berry.phtml?topic=fprom&group=programs&subgroup=promoter)在線軟件預(yù)測啟動(dòng)子。

膚色和表達(dá)量數(shù)據(jù)采用SPSS 21.0軟件進(jìn)行方差齊性檢驗(yàn)，對齊性數(shù)據(jù)用單因素方差分析(Duncan法)，對不齊性數(shù)據(jù)用獨(dú)立樣本非參數(shù)檢驗(yàn)(Kruskal-Wallis法)。表達(dá)量數(shù)據(jù)之間的相關(guān)性采用SPSS 21.0軟件中的雙變量相關(guān)分析(Pearson相關(guān)系數(shù))。突變檢測值與理論值的符合度采用SPSS 21.0軟件中的卡方檢驗(yàn)(Fisher’s法)。上述分析均以P<0.05為差異顯著，以P<0.01為差異極顯著。

2 結(jié) 果

2.1 突變篩查和位點(diǎn)情況

7個(gè)物種的進(jìn)化樹見圖1A，雞和鵪鶉具有較高同源性，因此設(shè)定都包含有雞和鵪鶉的5組比對序列，其保守堿基率曲線見圖1B。在0～-3 000 bp范圍內(nèi)，5組保守堿基率稍低于兩側(cè)區(qū)域；同時(shí)，從雞GNAS基因的SNP統(tǒng)計(jì)(圖1C)來看，0～-3 000 bp的SNP個(gè)數(shù)稍高于兩側(cè)，這些均表明GNAS基因的轉(zhuǎn)錄起始位點(diǎn)上游3 000 bp以內(nèi)更具篩查意義，因此選擇-3 000～-1區(qū)域進(jìn)行啟動(dòng)子突變探查。

A. 7個(gè)動(dòng)物物種根據(jù)GNAS基因生成的進(jìn)化樹，采用最大似然法的Tamura-Nei模型，以斑馬魚為遍歷根。B. 5組比對序列的每500 bp的保守堿基率曲線，X軸0表示GNAS基因的轉(zhuǎn)錄起始位點(diǎn)。C. 雞GNAS基因的每500 bp的SNP個(gè)數(shù)統(tǒng)計(jì)值曲線，X軸0表示轉(zhuǎn)錄起始位點(diǎn)A. Show the phylogenetic tree of 7 species by GNAS genetic distance, based on the Tamura-Nei model of maximum likelihood method, in zebrafish traversal root. B. Show the GNAS conserved rates for each 500 bp of 5 groups, X-axis “0” represents the transcription start site. C. Show the GNAS SNP statistical curve for each 500 bp, X-axis “0” represents the transcription start site圖1 GNAS基因進(jìn)化樹、保守堿基率和SNP統(tǒng)計(jì)Fig.1 The phylogenetic tree, conserved rates and SNP statistics of GNAS

經(jīng)探查和篩選，得到位于-2 270 bp處的突變位點(diǎn)(染色體上堿基位置11 167 660，rs號314048279)。該SNP為T/C突變(簡稱T2270C)，DraI酶切電泳后有3種條帶結(jié)果：CC純合基因型出現(xiàn)一條342 bp條帶；TC雜合基因型出現(xiàn)171和342 bp兩條帶；TT隱性純合基因型出現(xiàn)一條171 bp條帶，酶切圖和測序圖分別見圖2A和圖2B。

兩種軟件的啟動(dòng)子預(yù)測結(jié)果見表3，當(dāng)T2270C堿基(即查找序列的714處)為T時(shí)，在該點(diǎn)附近均預(yù)測到啟動(dòng)子；而T2270C突變?yōu)镃時(shí)，該點(diǎn)附近的啟動(dòng)子消失，表明T2270C是啟動(dòng)子的關(guān)鍵位點(diǎn)。

A. Dra I酶切電泳結(jié)果，CC、TC、TT分別代表3種基因型，M為DNA相對分子質(zhì)量標(biāo)準(zhǔn)。B. 3種基因型的測序結(jié)果，箭頭所指堿基為突變位點(diǎn)A. The electrophoresis result after Dra I enzyme digestion, CC, TC, TT represent the 3 genotypes respectively, M mean DNA marker. B. The sequencing result of the 3 genotypes, and the bases pointed by arrows are the mutation sites圖2 突變位點(diǎn)酶切和測序圖Fig.2 The digestion and sequencing maps of the mutation site

表3GNAS基因啟動(dòng)子預(yù)測結(jié)果

Table 3 The prediction results ofGNASgene promoter

軟件名稱SoftwareT2270C為T時(shí)的預(yù)測結(jié)果PredictionresultwhenT2270CwasTT2270C為C時(shí)的預(yù)測結(jié)果PredictionresultwhenT2270CwasCNeuralNetworkStartEndScorePromoterSequence7037530．94AGGTAGTCCTTAAAAACAGTGGTGCCTGGCTTC?ACCCTGGGCATATCAG160416540．90CATGTGATAATCAAAATTCCCCCCACTCTTGGACT?GCTGCATCTGCCCTGStartEndScorePromoterSequence160416540．90CATGTGATAATCAAAATTCCCCCCACTCTTGGAC?TGCTGCATCTGCCCTGFFROM1promoter/enhancer(s)arepredictedPromoterPos：746LDF：-0．722TATAboxat712+7．427TTAAAAA0promoter/enhancer(s)arepredicted

帶方框淺灰底色字符為突變位點(diǎn)

The character with border and light-gray shading are the mutation sites

2.2 突變與膚色的相關(guān)性

3種突變基因型的膚色L、a、b值見圖3(n=150)。組間L值差異極顯著(P<0.01)，CC基因型的L值最小，表現(xiàn)為膚色最深；TT基因型的L值最大，表現(xiàn)為膚色淺白。a為正值，b為負(fù)值，但3組間均無顯著差異(P>0.05)。

2.3 不同突變基因型的相關(guān)基因表達(dá)量

3種突變基因型的EDN3、MC1R、GNAS和ADCY3基因表達(dá)量見圖4(n=3×6)。EDN3基因表達(dá)量的組間差異極顯著(P<0.01)，以CC基因型的為最高，TT基因型的為最低。MC1R、GNAS和ADCY3基因表達(dá)量的組間差異均不顯著(P>0.05)，但以CC基因型的表達(dá)量最高。GNAS基因與上游MC1R基因的表達(dá)量之間的相關(guān)性不顯著，而與下游ADCY3基因表達(dá)量之間的相關(guān)性極顯著(P<0.01)。

柱形圖上方數(shù)據(jù)為各基因型的膚色平均值，肩標(biāo)小寫字母不同表示差異顯著(P<0.05)；大寫字母不同表示差異極顯著(P<0.01)。柱形圖的灰度根據(jù)各基因型的膚色L值模擬生成The data upon each column is the average value of skin color of each genotype, and the different superscript lowercase letters represent significant difference (P<0.05), meanwhile the different superscript uppercase letters represent extremely significant difference (P<0.01). The gray level is simulated according to the skin color L value of each genotype圖3 突變基因型的膚色值Fig.3 The skin color values of the mutation genotypes

圖中的散點(diǎn)表示各基因型的基因表達(dá)量，橫豎線分別表示平均值和標(biāo)準(zhǔn)差。EDN3基因上方的“**”表示基因型間的差異極顯著(P<0.01)。MC1R和GNAS、GNAS和ADCY3基因之間的數(shù)值為表達(dá)量的相關(guān)性值，“**”表示相關(guān)性極顯著(P<0.01)The scatters in each column represent the gene expression level of each genotype, and the horizontal and vertical line represent the mean and standard deviation, respectively. The “**” upon EDN3 gene represent extremely significant difference among 3 genotypes (P<0.01). The value between MC1R and GNAS, GNAS and ADCY3 are the correlation values of gene expression, and the “**”represent extremely significant difference (P<0.01)圖4 突變基因型的相關(guān)基因表達(dá)量Fig.4 The expression levels of the related genes of the mutation genotypes

2.4 不同雞種的突變檢測情況

3個(gè)雞種的GNAS突變檢測結(jié)果見表4。江山烏骨雞膚色全黑，經(jīng)檢測基本為CC基因型。伊莎蛋雞B系的膚色白，檢測全為TT基因型。將東鄉(xiāng)綠殼蛋雞公雞與伊莎蛋雞B系母雞進(jìn)行膚色測交，檢測無白膚后代的家系群體(有3個(gè)家系，3♂×18♀的后代共67只，均無白膚)基本為TC基因型?？ǚ綑z驗(yàn)表明，以上檢測結(jié)果與理論值相符(P>0.05)。

表4 雞種的GNAS突變檢測結(jié)果

Table 4 The detection results ofGNASgene mutation in chicken breeds

雞種Chickenbreed只數(shù)/只Number基因型GenotypeCCTCTTP值Pvalue江山烏骨雞JiangshanBlack?bonechicken6056400．119伊莎蛋雞B系ISAlayerBline600060-東鄉(xiāng)綠殼蛋雞×伊莎蛋雞B系DongxiangBlue?eggshellcock×ISAhen6706250．058

3 討 論

據(jù)報(bào)道，和GNAS同屬鳥嘌呤核苷酸結(jié)合蛋白(GNA)基因家族的GNAQ、GNA11基因的突變型與惡性黑色素瘤病人存在相關(guān)性[10-11]，而人GNAS基因第5外顯子T393C多態(tài)與黑色素瘤、膀胱癌、肺癌等惡性腫瘤生存期顯著相關(guān)[5,12-13]；第8外顯子R201H、R201C突變在纖維性發(fā)育不良病中的發(fā)生率達(dá)86%[7]。本研究篩查到了雞GNAS基因的啟動(dòng)子區(qū)突變，其不同基因型間的膚色明度值差異顯著，但在針對雞相應(yīng)于人GNAS的外顯子區(qū)的突變篩查中沒有發(fā)現(xiàn)類似位點(diǎn)。雖然人的黑色素瘤、纖維性發(fā)育不良和雞的纖維化黑色變(Fm)都牽涉機(jī)體的黑色素、纖維組織等成分，但與這些表象相關(guān)聯(lián)的GNAS突變并不一樣，這應(yīng)是因?yàn)槲锓N以及具體性狀不同所致。本研究同時(shí)檢測了與黑色素沉積密切相關(guān)的EDN3基因表達(dá)量，發(fā)現(xiàn)在GNAS基因T2270C突變的CC純合基因型中，其EDN3表達(dá)量要高于TT隱性基因型10倍以上，這一結(jié)果與B. Dorshorst等關(guān)于絲羽烏骨雞(烏膚純合)、新漢普夏雞(白膚對照)的檢測結(jié)論一致[3]，也反映出該GNAS突變確實(shí)與黑色素以及烏膚性狀有關(guān)。

就GNAS的上下游通路基因而言，MC1R已被證實(shí)是控制雞黑色素性狀的重要基因[14]，但主要關(guān)聯(lián)羽色，本研究檢測到MC1R在皮膚中的表達(dá)量低，且與GNAS表達(dá)量相關(guān)性不強(qiáng)；ADCY3作為GNAS下游的結(jié)合效應(yīng)物，直接受其調(diào)節(jié)，檢測發(fā)現(xiàn)二者的基因表達(dá)量呈現(xiàn)出強(qiáng)相關(guān)性。據(jù)此推測，GNAS啟動(dòng)子突變可能是通過影響本底以及下游基因的表達(dá)，最終反映到膚色差異，但更明晰的分子調(diào)控機(jī)制還需要深入研究發(fā)掘。啟動(dòng)子突變一般通過轉(zhuǎn)錄因子結(jié)合而調(diào)控靶基因表達(dá)，比如TERT基因啟動(dòng)子的C228T、C250T突變與人黑色素瘤強(qiáng)相關(guān)，機(jī)制是突變產(chǎn)生ETS轉(zhuǎn)錄因子的共有結(jié)合位點(diǎn)，從而增強(qiáng)了基因的轉(zhuǎn)錄活性[15-16]；在對雞的Piwil1基因NF-Y轉(zhuǎn)錄因子、L-FABP基因C/EBPα轉(zhuǎn)錄因子實(shí)行定點(diǎn)突變后，啟動(dòng)子活性呈顯著下降或升高[17-18]。因此，后續(xù)可重點(diǎn)分析該GNAS突變附近的轉(zhuǎn)錄因子，研究突變引起的啟動(dòng)子活性變化情況。

當(dāng)前，家禽育種領(lǐng)域內(nèi)已有一些基因突變的研究應(yīng)用，例如K. L. Wells等通過全基因組SNP掃描發(fā)現(xiàn)無毛雞品系的FGF20基因存在R179X突變，并開發(fā)出dCAPS標(biāo)記方法進(jìn)行檢測[19]；Q. Chu等通過FMO3基因T329S突變的TaqMan探針法檢測，在北京油雞群體中逐步剔除含魚腥味缺陷的個(gè)體[20]。本研究的GNAS基因突變在東鄉(xiāng)綠殼蛋雞中得到性狀關(guān)聯(lián)驗(yàn)證，即CC純合型的烏膚深色度要顯著高于TC雜合型。因此在早期選擇時(shí)，可結(jié)合GNAS分子標(biāo)記，對東鄉(xiāng)綠殼蛋雞進(jìn)行膚色分型，提高選種效率。進(jìn)一步地，可擴(kuò)展研究并應(yīng)用在其他烏骨雞種中，實(shí)現(xiàn)膚色性狀的定向培育。

[1] WEINSTEIN L S, LIU J, SAKAMOTO A, et al. Minireview: GNAS: normal and abnormal functions [J].Endocrinology, 2004, 145(12): 5459-5464.

[3] DORSHORST B, MOLIN A M, RUBIN C J, et al. A complex genomic rearrangement involving the endothelin 3 locus causes dermal hyperpigmentation in the chicken[J].PLoSGenet,2011,7(12):e1002412.

[4] SHINOMIYA A, KAYASHIMA Y, KINOSHITA K, et al. Gene duplication of endothelin 3 is closely correlated with the hyperpigmentation of the internal organs (Fibromelanosis) in Silky chickens [J].Genetics, 2012, 190(2): 627-638.

[5] FREY U H, FRITZ A, ROTTERDAM S, et al. GNAS1 T393C polymorphism and disease progression in patients with malignant melanoma [J].EurJMedRes, 2010, 15(10): 422-427.

[6] TABAREAU-DELALANDE F, COLLIN C, GOMEZ-BROUCHET A, et al. Diagnostic value of investigatingGNASmutations in fibro-osseous lesions: a retrospective of 91 cases of fibrous dysplasia and 40 other fibro-osseous lesions [J].ModPathol, 2013, 26(7): 911-921.

[7] SHI R R, LI X F, ZHANG R, et al. GNAS mutational analysis in differentiating fibrous dysplasia and ossifying fibroma of the jaw [J].ModPathol, 2013, 26(8): 1023-1031.

[8] 國家畜禽遺傳資源委員會(huì).中國畜禽遺傳資源志·家禽志 [M]. 北京：中國農(nóng)業(yè)出版社, 2011.

CHINA NATIONAL COMMISION OF ANIMAL RESOURCES. Animal genetic resources in China·chicken [M].Beijing: China Agriculture Press, 2011. (in Chinese)

[9] TAMURA K, PETERSON D, PETERSON N, et al. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods [J].MolBiolEvol, 2011, 28(10): 2731-2739.

[10] VAN RAAMSDONK C D, BEZROOKOVE V, GREEN G, et al. Frequent somatic mutations of GNAQ in uveal melanoma and blue naevi [J].Nature, 2009, 457(7229): 599-602.

[11] VAN RAAMSDONK C D, GRIEWANK K G, CROSBY M B, et al. Mutations in GNA11 in uveal melanoma [J].NEnglJMed, 2010, 363(23): 2191-2199.

[12] FREY U H, EISENHARDT A, LüMMEN G, et al. The T393C polymorphism of the Gαs gene (GNAS1) is a novel prognostic marker in bladder cancer [J].CancerEpidemiolBiomarkersPrev, 2005, 14(4): 871-877.

[13] 龔虹云.GNAS1 T393C單核苷酸多態(tài)性與肺癌生存及放化療副反應(yīng)相關(guān)性研究 [D]. 武漢: 武漢大學(xué), 2012.

GONG H Y. Association of the GNAS1 T393C polymorphism with survival in advanced lung cancer and side effects in chemotherapy and radiotherapy [D]. Wuhan: Wuhan University, 2012. (in Chinese)

[14] 楊永升, 鄧學(xué)梅, 李 寧, 等.MC1R是控制雞黑色素形成的候選主效基因 [J]. 生物化學(xué)與生物物理進(jìn)展, 2004, 31(6): 500-505.

YANG Y S, DENG X M, LI N, et al. MC1R is the candidate gene regulating melanin synthesis in chicken[J].ProgressinBiochemistryandBiophysics, 2004, 31(6): 500-505.(in Chinese)

[15] HUANG F W, HODIS E, XU M J, et al. Highly recurrent TERT promoter mutations in human melanoma [J].Science,2013, 339(6122): 957-959.

[16] HORN S, FIGL A, RACHAKONDA P S, et al. TERT promoter mutations in familial and sporadic melanoma [J].Science, 2013, 339(6122): 959-961.

[17] 夏明秀, 常國斌, 陳 蓉, 等. 雞Piwil1基因啟動(dòng)子區(qū)轉(zhuǎn)錄調(diào)控元件的初步分析 [J]. 畜牧獸醫(yī)學(xué)報(bào), 2013, 44(9): 1349-1354.

XIA M X, CHANG G B, CHEN R, et al. Preliminary analysis of transcription regulation elements in the promoter region ofPiwil1 gene in chicken[J].ActaVeterinariaetZootechnicaSinica, 2013, 44(9): 1349-1354. (in Chinese)

[18] 賀 綦, 孫嬰寧, 李 輝, 等. 雞L-FABP啟動(dòng)子區(qū)C/EBPα結(jié)合位點(diǎn)的定點(diǎn)突變分析 [J]. 畜牧獸醫(yī)學(xué)報(bào), 2015, 46(9): 1496-1501.

HE Q, SUN Y N, LI H, et al. Site directed mutagenesis for C/EBPα binding sites in chickenL-FABPpromoter[J].ActaVeterinariaetZootechnicaSinica, 2015, 46(9): 1496-1501. (in Chinese)

[19] WELLS K L, HADAD Y, BEN-AVRAHAM D, et al. Genome-wide SNP scan of pooled DNA reveals nonsense mutation in FGF20 in the scaleless line of featherless chickens [J].BMCGenomics, 2012, 13:257.

[20] CHU Q, ZHANG J, ZHU S, et al. The detection and elimination of flavin-containing monooxygenase 3 gene T329S mutation in the Beijing You chicken [J].PoultSci, 2013, 92(12): 3109-3112.

(編輯 郭云雁)

GNAS Gene Promoter Mutation in Chicken and the Correlation with Skin Color Traits

WANG Huan-huan1，CHEN Mei-ling1，LOU Li-feng1，ZHANG Lei1，ZHANG Cheng-xian2，CHEN Xian-hui2，ZHANG Xue-dong1*

(1. Hangzhou Academy of Agricultural Sciences，Hangzhou 310024，China；2.ZhejiangGuangdaPoultryBreedingCompany，Haining310016，China)

The study aimed to investigate the mutation ofGNASgene and its effect on skin color traits in chicken. Dongxiang Blue-eggshell chicken (Gallusgallusdomestica) was chose and theGNASmutation sites related to skin color were screened using PCR-restriction length polymorphism (RFLP) and DNA sequencing in this study. We detected expression of theEDN3,MC1R,GNASandADCY3 genes using fluorescence quantitative PCR, and the mutation sites among 3 chicken breeds. The results showed that there was a T/C mutation at base 11 167 660 on chicken chromosome 20 and the mutation (T2270C) was crucial for theGNASgene promoter. AfterDraI enzyme digestion for PCR products, the CC homozygous genotype, the TC heterozygous genotype and the TT recessive homozygous genotype presented one band at 342 bp, two bands at 171 and 342 bp, and one band at 171 bp, respectively. The L values of skin color were statistically significantly different among 3 genotypes (P<0.01). The individuals with CC genotype had the lowest L value, which presented the darkest skin, and the individuals with TT genotype had the highest L value, which presented the lightest skin. Expression of theEDN3 gene was highest and lowest in individuals with the CC and TT genotype, respectively, and the difference among 3 genotypes was extremely significant (P<0.01). Expression of theMC1R,GNASandADCY3 genes was highest in individuals with CC genotype, but the difference among genotypes was not significant (P>0.05). The correlation ofGNASwithMC1Rwas non-significant in expression, and the correlation withADCY3 was extremely significant (P<0.01). The typing results according to mutation detection in Jiangshan Black-bone chicken, ISA layer, test-crossing offspring of Dongxiang Blue-eggshell cock and ISA hen, conformed to the theoretical values (P>0.05). The study indicate that the T2270C mutation inGNASgene promoter in chicken is correlated strongly with the skin color traits, and detection of the mutation can be used in the breeding of black-bone chicken.

chicken；GNASgene；mutation；skin color

10.11843/j.issn.0366-6964.2016.12.004

2016-05-16

國家肉雞產(chǎn)業(yè)技術(shù)體系綜合試驗(yàn)站(CARS-42-Z07)；浙江省公益技術(shù)研究農(nóng)業(yè)項(xiàng)目(2015C32106)

王歡歡(1986-)，女，黑龍江佳木斯人，畜牧師，碩士，主要從事動(dòng)物遺傳育種工作，E-mail: wennuan2009@163.com

*通信作者：章學(xué)東，研究員，E-mail：bigzhengliang@hotmail.com

S831.2

A

0366-6964(2016)12-2354-08