劉坤舉，張小輝，*，龐有志，趙淑娟，祁艷霞，王乾昆

(1.河南科技大學(xué) 動(dòng)物科技學(xué)院，河南 洛陽(yáng) 471003； 2. 洛陽(yáng)市動(dòng)物遺傳育種重點(diǎn)實(shí)驗(yàn)室，河南 洛陽(yáng) 471003)

羽色是家禽重要的質(zhì)量性狀，其遺傳效應(yīng)是近年來(lái)研究的熱點(diǎn)。黑色素的數(shù)量、性質(zhì)與分布是導(dǎo)致家禽羽色多樣性的重要原因[1]。作為羽色形成的基本物質(zhì)，黑色素可分為真黑素和褐黑素。真黑素包括黑色和褐色，褐黑素包括黃色和紅色，2種黑色素之間的相互轉(zhuǎn)化影響羽色的形成[2]。影響黑色素形成的基因有很多，鳥(niǎo)嘌呤核苷酸結(jié)合蛋白α亞基(guanine nucleotide-binding protein Gs subunit alpha，GNAS)是位于細(xì)胞膜內(nèi)側(cè)的一類信號(hào)傳導(dǎo)蛋白，參與細(xì)胞的多種生命活動(dòng)[3]。在黑色素代謝通路中，GNAS與黑素皮質(zhì)素受體1(melanocortin receptor-1，MC1R)通過(guò)跨膜結(jié)構(gòu)結(jié)合，激活腺苷酸環(huán)化酶(adenylate cyclase，ADCY)，催化ATP產(chǎn)生第二信使，影響細(xì)胞內(nèi)黑色素的合成[4]。

GNAS基因是影響烏骨雞黑色素沉積的候選基因，該基因位于20號(hào)染色體Fm區(qū)域的拷貝數(shù)重復(fù)區(qū)(copy number variations，CNVs)內(nèi)，與黑色素的沉積密切相關(guān)[5]。已有學(xué)者發(fā)現(xiàn)，GNAS基因突變可能會(huì)導(dǎo)致人類患黑色素瘤病[6-7]。綜上，GNAS基因可能與黑色素形成有關(guān)。

本研究采用實(shí)時(shí)熒光定量PCR技術(shù)探究栗羽和白羽朝鮮鵪鶉不同胚胎發(fā)育階段GNAS基因的mRNA相對(duì)表達(dá)量，并對(duì)GNAS基因單核苷酸多態(tài)位點(diǎn)(single nucleotide polymorphisms，SNPs)進(jìn)行檢測(cè)，分析SNPs與朝鮮鵪鶉羽色性狀之間的關(guān)聯(lián)性，旨在為探究朝鮮鵪鶉及其他禽類羽色形成、變化等研究提供參考。

1 材料與方法

1.1 試驗(yàn)材料

栗羽和白羽朝鮮鵪鶉受精蛋購(gòu)自市場(chǎng)。取栗羽和白羽朝鮮鵪鶉受精蛋各60枚進(jìn)行孵化，取孵化至6、8、10、12、14 d的鵪鶉胚胎翅尖組織置于液氮中保存，用于后續(xù)總RNA的提??；取栗羽和白羽朝鮮鵪鶉受精蛋各300枚，孵化至10 d，取胚胎翅尖組織置于液氮中保存，用于后續(xù)基因組DNA的提取。

1.2 主要試劑和儀器

Ezup柱式動(dòng)物基因組DNA抽提試劑盒、總RNA反轉(zhuǎn)錄試劑盒、瓊脂糖凝膠回收試劑盒、pMD-19T質(zhì)粒載體、SYBR Green II RNA染料以及DL 2000 bp Maker均購(gòu)置于北京索萊寶科技有限公司。

Nano Drop 2000C超微量分光光度計(jì)和冷凍高速離心機(jī)購(gòu)于美國(guó)Thermo Fisher Scientific公司；熒光定量PCR儀、梯度PCR儀、凝膠成像儀均購(gòu)于美國(guó)BIO-RAD公司；水平電泳槽和垂直電泳槽購(gòu)于北京六一生物科技有限公司。

1.3 總RNA提取與反轉(zhuǎn)錄

分別取栗羽和白羽朝鮮鵪鶉不同胚胎時(shí)期的翅尖組織，將同羽色、同孵化期的3個(gè)樣品混合，用電動(dòng)研磨器磨碎，采用Trizol法提取總RNA。取1.5 μL樣品滴于超微量分光光度計(jì)上檢測(cè)，合格后將樣品稀釋至500 ng·μL-1。取2 μL稀釋后RNA樣品混合1 μL 6×DNA Loading Buffer，進(jìn)行瓊脂糖凝膠電泳，檢測(cè)總RNA的質(zhì)量。

使用反轉(zhuǎn)錄試劑盒合成cDNA第一鏈。反應(yīng)條件是：總RNA 3 μL、gDNA remover 1 μL、10×g DNA remover Buffer 1 μL，補(bǔ)RNase-free H2O至10 μL，42 ℃孵育2 min，60℃孵育5 min。第二步的合成體系及反應(yīng)條件是：dNTP Mix 1μL、GoldenstarTMOligo(dT)171 μL、Randomer 1 μL、5×GoldenstarTMBuffer 4 μL、DTT 1 μL、GoldenstarTMRT6 1 μL、RNase-free H2O 2 μL，25 ℃孵育10 min，55 ℃孵育40 min，85 ℃孵育5 min，產(chǎn)物置于-20 ℃保存。

1.4 基因組DNA提取與檢測(cè)

分別取栗羽和白羽朝鮮鵪鶉胚胎孵化至10 d時(shí)的翅尖組織，使用基因組DNA抽提試劑盒進(jìn)行提取。取1.5 μL樣品滴在超微量分光光度計(jì)上檢測(cè)，取4 μL合格的DNA樣品混合1 μL 6×DNA Loading Buffer進(jìn)行電泳，檢測(cè)基因組DNA的質(zhì)量。

1.5 GNAS基因mRNA表達(dá)檢測(cè)

采用克隆獲得的朝鮮鵪鶉GNAS基因序列設(shè)計(jì)引物mRNA序列信息，用Beacon Designer 8.0軟件設(shè)計(jì)引物，交由北京擎科新業(yè)生物技術(shù)有限公司合成，引物序列見(jiàn)表1。反應(yīng)體系為：2×TaqPlus Master Mix 10 μL，上下游引物各0.8 μL，SYBR Green Ⅰ 1 μL，ddH2O 6.4 μL，cDNA模板1 μL。反應(yīng)條件為：94 ℃預(yù)變性4 min；94 ℃變性30 s，60 ℃退火30 s，72 ℃延伸30 s，40個(gè)循環(huán)；延伸10 min；熔解曲線從65 ℃到95 ℃，每5 s升高0.5 ℃。

1.6 GNAS基因CDS區(qū)克隆與分析

根據(jù)NCBI數(shù)據(jù)庫(kù)日本鵪鶉和原雞GNAS基因的保守序列設(shè)計(jì)2對(duì)引物(表2)。反應(yīng)體系為：1.1×T3 Super PCR Mix 17.6 μL，上下游引物各0.7 μL，模板cDNA 1 μL。反應(yīng)條件為：94 ℃預(yù)變性4 min；94 ℃變性15 s，退火15 s，72 ℃延伸15 s，36個(gè)循環(huán)；最后72 ℃延伸10 min。得到的PCR產(chǎn)物經(jīng)過(guò)1.5%的瓊脂糖凝膠電泳檢測(cè)后進(jìn)行膠回收，經(jīng)過(guò)膠回收的產(chǎn)物測(cè)序后運(yùn)用SeqMan軟件對(duì)序列信息進(jìn)行拼接獲得完整的CDS區(qū)序列?；厥蘸蟮漠a(chǎn)物連接到pMD-19T質(zhì)粒載體后，轉(zhuǎn)化到大腸埃希菌感受態(tài)細(xì)胞中，于恒溫培養(yǎng)箱37 ℃條件下培養(yǎng)8 h，藍(lán)白斑篩選后，選取陽(yáng)性菌落進(jìn)行PCR鑒定檢測(cè)，將鑒定為陽(yáng)性的菌落送北京擎科新業(yè)生物技術(shù)有限公司測(cè)序。將測(cè)序結(jié)果分別與NCBI的基因組進(jìn)行BLAST分析，用ExPASy預(yù)測(cè)蛋白序列編碼氨基酸、等電點(diǎn)與分子質(zhì)量等。使用MEGA 7.0軟件對(duì)朝鮮鵪鶉進(jìn)行生物種屬間比對(duì)，生成系統(tǒng)進(jìn)化樹(shù)。

表1 實(shí)時(shí)熒光定量PCR引物信息

表2 GNAS基因CDS區(qū)引物設(shè)計(jì)

1.7 朝鮮鵪鶉GNAS基因多態(tài)性檢測(cè)

1.7.1GNAS基因的不對(duì)稱PCR擴(kuò)增

采用克隆獲得的朝鮮鵪鶉GNAS基因序列設(shè)計(jì)外顯子12部分引物，參考NCBI數(shù)據(jù)庫(kù)中的日本鵪鶉GNAS基因序列信息設(shè)計(jì)3’非編碼區(qū)部分引物，用Primer 5.0軟件設(shè)計(jì)出2對(duì)引物，由北京擎科新業(yè)生物技術(shù)有限公司合成(表3)。

表3 朝鮮鵪鶉GNAS基因引物序列信息

采用不對(duì)稱PCR進(jìn)行擴(kuò)增，以提取的基因組DNA為模板，等濃度引物進(jìn)行PCR擴(kuò)增，擴(kuò)增出目的片段的雙鏈，檢測(cè)產(chǎn)物質(zhì)量；再以PCR產(chǎn)物作為模板，利用其中的任意1條引物進(jìn)行第2次PCR擴(kuò)增，擴(kuò)增出目的片段的單鏈。外顯子12部分片段的第1次PCR擴(kuò)增反應(yīng)體系為：基因組DNA 1 μL，1.1×T3 Super PCR Mix 17.6 μL，上下游引物各0.7 μL；PCR擴(kuò)增反應(yīng)條件：95 ℃預(yù)變性4 min；95 ℃變性15 s，59 ℃退火15 s，72 ℃延伸15 s，36個(gè)循環(huán)；72 ℃延伸10 min。1.8%瓊脂糖凝膠電泳檢測(cè)產(chǎn)物，得到目的條帶后進(jìn)行第2次PCR擴(kuò)增。第2次PCR擴(kuò)增反應(yīng)體系為：第1次PCR擴(kuò)增產(chǎn)物1 μL，1.1×T3 Super PCR Mix18.3 μL，上游引物0.7 μL，反應(yīng)條件同第1次PCR，循環(huán)數(shù)調(diào)整至20個(gè)。3′UTR部分?jǐn)U增的反應(yīng)體系和反應(yīng)條件除退火溫度為59 ℃外，其他與外顯子12部分一致。

1.7.2 聚丙烯酰胺凝膠電泳與分析

取1 μL SYBR Green II RNA染料加入1 mL滅菌雙蒸水中混合，再加入1 mL 6×Loading buffer緩沖液與其混勻，取5 μL上述工作液加入到20 μL電泳樣品中，94 ℃變性使兩者充分結(jié)合。取8 μL混合液，于常溫下濃度為12%的聚丙烯酰胺凝膠進(jìn)行電泳，185 V電泳10 min后轉(zhuǎn)換為85 V繼續(xù)電泳8 h，結(jié)束后在凝膠成像儀上進(jìn)行拍照保存。通過(guò)電泳結(jié)果對(duì)單鏈構(gòu)象多態(tài)性(single strand conformation polymorphism，SSCP)分析，相同基因型對(duì)應(yīng)個(gè)體各選取3份PCR擴(kuò)增產(chǎn)物送至北京擎科新業(yè)生物技術(shù)有限公司測(cè)序。測(cè)序結(jié)果使用DNAMAN軟件和Ensembl網(wǎng)站進(jìn)行對(duì)比分析，以確定突變位點(diǎn)。根據(jù)基因型統(tǒng)計(jì)結(jié)果，計(jì)算基因型頻率、基因頻率和多態(tài)信息含量(polymorphic information content，PIC)，采用卡方檢驗(yàn)檢測(cè)朝鮮鵪鶉群體是否處于Hardy-Weinberg平衡。

1.8 數(shù)據(jù)的統(tǒng)計(jì)學(xué)分析

用2-ΔΔCt法計(jì)算相對(duì)表達(dá)量，運(yùn)用SPSS 20.0進(jìn)行樣本t檢驗(yàn)，Graphpad Prism 8.0軟件作圖，用SPSS 20.0的χ2檢驗(yàn)SNP位點(diǎn)與羽色性狀的關(guān)聯(lián)性和群體Hardy-Weinberg平衡。

2 結(jié)果與分析

2.1 總RNA和基因組DNA質(zhì)量檢測(cè)

經(jīng)超微量分光光度計(jì)檢測(cè)，提取的總RNA的濃度在2 000～3 500 ng·μL-1，A260/A280均在2.0左右，A260/A230均在2.2左右。1.2%瓊脂糖凝膠檢測(cè)如圖1，可清晰看到28S、18S、5S 3條帶，表明所提取的總RNA質(zhì)量良好。

1～5：RNA驗(yàn)證結(jié)果。1-5：RNA verification results.圖1 總RNA凝膠成像Fig.1 Gel imaging of total RNA

對(duì)提取的基因組DNA檢測(cè)，A260/A280均在1.9左右，A260/A230均在2.2左右。0.8%瓊脂糖凝膠檢測(cè)如圖2，條帶清晰明亮且長(zhǎng)度符合DNA要求，表明所提取的DNA質(zhì)量良好。

1～5：DNA驗(yàn)證結(jié)果。1-5：DNA verification results.圖2 基因組DNA凝膠成像Fig.2 Gel imaging of genomic DNA

2.2 GNAS基因mRNA表達(dá)水平

朝鮮鵪鶉不同發(fā)育階段胚胎中GNAS基因mRNA相對(duì)表達(dá)水平見(jiàn)圖3。如圖所示，GNAS基因在朝鮮鵪鶉胚胎發(fā)育至6 d時(shí)的表達(dá)水平低于發(fā)育到8～14 d；當(dāng)發(fā)育至8 d時(shí)，表達(dá)量達(dá)到最高；GNAS基因在栗羽鵪鶉胚胎組織中的表達(dá)水平整體高于白羽朝鮮鵪鶉，其中胚胎發(fā)育到8～14 d時(shí)，GNAS基因在栗羽鵪鶉胚胎組織中的表達(dá)水平極顯著(P<0.01)高于白羽鵪鶉。

**表示相同發(fā)育階段不同羽色鵪鶉胚胎間差異極顯著(P<0.01)。** means great significant difference(P<0.01) between different plumage color Korean quail embryos at same developmental stage.圖3 朝鮮鵪鶉不同發(fā)育階段胚胎中GNAS基因mRNA相對(duì)表達(dá)水平Fig.3 Relative expression of GNAS gene mRNA in different developmental stages of Korean quail embryos

2.3 GNAS基因CDS區(qū)克隆及序列分析

GNAS基因CDS區(qū)的PCR擴(kuò)增結(jié)果如圖4所示，其片段大小分別為1 001、439 bp，產(chǎn)物條帶清晰，可以用于DNA測(cè)序分析。測(cè)序結(jié)果經(jīng)BLAST比對(duì)，成功獲得GNAS基因CDS區(qū)完整序列(GenBank No. MT649634)，其中開(kāi)放閱讀框(open reading frame，ORF)長(zhǎng)1 140 bp。GNAS基因位于鵪鶉20號(hào)染色體上，有14個(gè)外顯子，編碼的蛋白含379個(gè)氨基酸，等電點(diǎn)6.43，分子質(zhì)量為44 399.47 Da，帶正電荷氨基酸殘基數(shù)為54，帶負(fù)電荷氨基酸殘基數(shù)為56，GRAVY值為-0.593，該蛋白為穩(wěn)定親水性蛋白。

M：DNA marker 2 000；1～2：CDS區(qū)擴(kuò)增產(chǎn)物。M: DNA marker 2 000；1-2: CDS region amplification products.圖4 朝鮮鵪鶉GNAS基因CDS區(qū)產(chǎn)物的電泳檢測(cè)Fig.4 Electrophoresis detection of PCR product of GNAS CDS in Korean quail

2.4 GNAS基因編碼蛋白的系統(tǒng)進(jìn)化分析

由圖5可知，朝鮮鵪鶉GNAS基因編碼蛋白與日本鵪鶉和原雞同源性最高，表明GNAS基因在物種進(jìn)化過(guò)程中保持了較高的保守性。

圖5 基于GNAS蛋白的朝鮮鵪鶉種屬間系統(tǒng)進(jìn)化樹(shù)Fig.5 Interspecific phylogenetic tree of Korean quail based on the protein of GNAS

2.5 GNAS基因特定序列PCR擴(kuò)增產(chǎn)物檢測(cè)

圖6為不同羽色朝鮮鵪鶉GNAS基因外顯子12和3′UTR的PCR擴(kuò)增結(jié)果。外顯子12片段擴(kuò)增長(zhǎng)度254 bp，3′UTR的擴(kuò)增產(chǎn)物長(zhǎng)度為190 bp，產(chǎn)物條帶清晰且無(wú)雜帶，可以用于回收測(cè)序。

2.6 GNAS基因型分析

經(jīng)PCR-SSCP分析，GNAS基因外顯子12和3′UTR均有3種基因型(圖7)。圖中，左側(cè)是外顯子12基因型，右側(cè)為3′UTR基因型。

分別選取外顯子12和3′UTR的3種基因型個(gè)體PCR產(chǎn)物進(jìn)行測(cè)序。如圖8所示，外顯子12區(qū)域和3′UTR區(qū)域均有1個(gè)突變位點(diǎn)，分別是g.119221T>A和g.121181A>G。

2.7 突變位點(diǎn)遺傳信息分析

如表4所示，朝鮮鵪鶉群體中GNAS基因的g.119221T>A的優(yōu)勢(shì)等位基因?yàn)锽，優(yōu)勢(shì)等位基因頻率為0.75；g.121181A>G突變位點(diǎn)的優(yōu)勢(shì)等位基因是A，頻率為0.71。通過(guò)對(duì)朝鮮鵪鶉群體多態(tài)信息含量計(jì)算可知，GNAS基因的2個(gè)突變位點(diǎn)均表現(xiàn)為中度多態(tài)(0.250.05)。

1～4、9～12：栗羽鵪鶉PCR擴(kuò)增產(chǎn)物；5～8、13～17：白羽鵪鶉PCR擴(kuò)增產(chǎn)物；M：DNA marker 2 000。1-4， 9-12: PCR amplification products of maroon plumage quail; 5-8， 13-17: PCR amplification product of white plumage quail;M: DNA marker 2 000.圖6 朝鮮鵪鶉GNAS基因外顯子12和3′UTR PCR產(chǎn)物電泳Fig.6 Electrophoresis detection of PCR product of exon 12 and 3′UTR in GNAS gene of Korean quail

圖7 朝鮮鵪鶉GNAS基因擴(kuò)增產(chǎn)物不對(duì)稱PCR-SSCP分析Fig.7 Genotypes of GNAS gene in Korean quail by ASPCR-SSCP analysis

圖8 朝鮮鵪鶉GNAS基因的突變檢測(cè)Fig.8 Detections of GNAS mutation in Korean quail

2.8 GNAS基因多態(tài)性與羽色性狀的關(guān)聯(lián)分析

如表5所示，在g.119221T>A位點(diǎn)，栗羽朝鮮鵪鶉3種基因型的頻率分布與白羽朝鮮鵪鶉差異不顯著(P>0.05)；在g.121181A>G突變位點(diǎn)，栗羽朝鮮鵪鶉的3種基因型的頻率分布與白羽朝鮮鵪鶉差異極顯著(P<0.01)，等位基因A能顯著提高羽色沉積度，其中AA基因型個(gè)體>AB基因型個(gè)體>BB基因型個(gè)體。由此可以推斷，GNAS基因3′UTR的g.121181A>G突變位點(diǎn)與朝鮮鵪鶉羽色性狀具有關(guān)聯(lián)性。

表4 不同基因型群體分析

表5 朝鮮鵪鶉GNAS基因多態(tài)性與羽色性狀的相關(guān)性分析

3 討論

羽色的遺傳效應(yīng)是近年來(lái)許多學(xué)者研究的重點(diǎn)。類胡蘿卜素類色素和黑色素類色素是沉積在皮膚中的2種主要色素，其相互影響形成了動(dòng)物界豐富多彩的顏色。類胡蘿卜素類色素主要沉積在皮膚和卵巢中[8]；真黑色素和褐黑素統(tǒng)稱黑色素，主要存在于黑色素細(xì)胞中，黑色素細(xì)胞主要分布于結(jié)締組織、骨骼以及骨膜等各種組織中[9]，黑色素細(xì)胞在真皮層和表皮層沉積的位置不同對(duì)膚色影響也不盡相同[10]。

在黑色素皮質(zhì)素代謝途徑中，MC1R與配體結(jié)合后被激活，催化產(chǎn)生第二信使cAMP，促進(jìn)真黑素的形成[4]。當(dāng)黑色素細(xì)胞中的MC1R含量不足或者活動(dòng)受阻時(shí)，黑色素的含量就會(huì)減少，動(dòng)物的皮膚等顏色就表現(xiàn)為淺色[11-12]。本研究結(jié)果表明，GNAS基因在栗羽朝鮮鵪鶉組織中表達(dá)水平高于白羽朝鮮鵪鶉，與趙新[14]提出的GNAS基因mRNA在黑色綿羊皮膚中表達(dá)水平極顯著高于白色綿羊皮膚組織的觀點(diǎn)一致。GNAS基因在朝鮮鵪鶉胚胎發(fā)育至8 d時(shí)表達(dá)量達(dá)到最大，此時(shí)該基因處于活躍時(shí)期，催化激活第二信使并促進(jìn)真黑素的合成，當(dāng)色素的合成和沉積達(dá)到一定水平后會(huì)影響到動(dòng)物的膚色表型。朝鮮鵪鶉GNAS基因的CDS區(qū)序列與其他物種同源性對(duì)比，發(fā)現(xiàn)GNAS基因CDS區(qū)序列與日本鵪鶉和原雞等禽類具有非常高的同源性，在系統(tǒng)發(fā)育樹(shù)上先聚合，緊接著與哺乳類以及其他物種再聚合，反映出GNAS基因編碼蛋白存在著物種間的差異但同時(shí)保持了高度保守性。

據(jù)報(bào)道，與GNAS同屬一個(gè)家族的GNAQ和GNA11基因的變異與惡性黑色毒瘤的發(fā)生有密切的相關(guān)性[14-15]，說(shuō)明基因在體內(nèi)黑色素細(xì)胞的代謝通路中，發(fā)揮著重要作用。Weinstein等[3]報(bào)道，McCune-Albright綜合征患者的GNAS基因第8外顯子產(chǎn)生了錯(cuò)義突變，造成皮膚黑色素增多并出現(xiàn)咖啡樣色素斑。同時(shí)人的GNAS基因的第5外顯子上的突變位點(diǎn)與黑色素瘤、膀胱癌以及肺癌等惡性腫瘤患者的生存期顯著相關(guān)[16-17]。通過(guò)以上報(bào)道可以得知GNAS基因與黑色素的形成有密切的關(guān)系。本次試驗(yàn)在朝鮮鵪鶉GNAS基因上檢測(cè)到了2個(gè)SNP位點(diǎn)，2個(gè)突變位點(diǎn)均表現(xiàn)為中度多態(tài)，且都處于Hardy-Weinberg平衡狀態(tài)。由于體內(nèi)的黑色素形成和作用機(jī)理異常復(fù)雜，因此，在這一過(guò)程中黑色素的合成和沉積，以及GNAS基因的表達(dá)和變異都會(huì)對(duì)羽色多樣性產(chǎn)生影響。經(jīng)對(duì)比可看出這2個(gè)SNP位點(diǎn)均未能引起氨基酸的改變，但外顯子12區(qū)域的g.119221T>A突變?yōu)橥x突變，可能會(huì)引起RNA結(jié)構(gòu)發(fā)生變化進(jìn)而導(dǎo)致立體結(jié)構(gòu)改變而影響其理化性質(zhì)，從而可能會(huì)影響后續(xù)的轉(zhuǎn)錄和翻譯過(guò)程[18]；目前對(duì)于MicroRNA的研究集中在3′UTR區(qū)域[19]，本試驗(yàn)的3′UTR區(qū)域的g.121181A>G的突變位點(diǎn)與MicroRNA結(jié)合有關(guān)，進(jìn)而會(huì)影響基因轉(zhuǎn)錄活性等變化。在g.121181A>G突變位點(diǎn)，栗羽朝鮮鵪鶉3種基因型的頻率分布與白羽朝鮮鵪鶉差異極顯著(P<0.01)，等位基因A能顯著提高羽色沉積度，說(shuō)明該突變位點(diǎn)與朝鮮鵪鶉羽色性狀具有一定的關(guān)聯(lián)性。