季 敏，劉學(xué)洪＊，余長(zhǎng)林，余選富，劉定江，史憲偉

（1.云南農(nóng)業(yè)大學(xué)動(dòng)物科學(xué)技術(shù)學(xué)院，云南 昆明 650201；2.云南省騰沖縣畜牧工作站，云南 騰沖 679100）

1 引言

檳榔江水牛是我國(guó)迄今為止發(fā)現(xiàn)的唯一的河流型水牛遺傳資源，分布在云南省騰沖縣及周邊。檳榔江水牛有著獨(dú)特的遺傳特性，尤其是在產(chǎn)奶性能比較優(yōu)秀。有研究結(jié)果顯示檳榔江水牛平均產(chǎn)奶量是德宏水牛的2倍多，比摩拉水牛高15%，但個(gè)體間生產(chǎn)性能差異較大［1］。章純熙等［2］對(duì)101頭檳榔江水牛的105個(gè)泌乳期泌乳量結(jié)果進(jìn)行測(cè)定，結(jié)果發(fā)現(xiàn)檳榔江水牛平均泌乳天數(shù)為270d，比其他河流型和沼澤型水牛平均泌乳天數(shù)短，而平均泌乳量為2 450kg，比摩拉和尼里拉菲等河流型水牛高8.4%～15.0%，比廣東水牛和溫州水牛等沼澤型水牛高124.4%～140.3%。另外通過(guò)分析騰沖水牛場(chǎng)提供的31頭母牛產(chǎn)奶數(shù)據(jù)資料，一個(gè)產(chǎn)奶周期平均產(chǎn)奶量為2 531.32±641.45kg，產(chǎn)奶量最高可達(dá)3 615.7 kg，最低的為1 543.4kg［3］。闡明檳榔江水牛產(chǎn)奶性狀的分子基礎(chǔ)，為今后開(kāi)展產(chǎn)奶性能標(biāo)記輔助選擇，提高檳榔江水牛的產(chǎn)奶性能具有重要的意義。

信號(hào)轉(zhuǎn)導(dǎo)及轉(zhuǎn)錄激活因子（STAT5A）主要參與哺乳動(dòng)物乳腺組織中細(xì)胞因子的信號(hào)轉(zhuǎn)導(dǎo)，對(duì)泌乳活動(dòng)具有重要的調(diào)控作用。STAT5A也被稱作乳腺因子，因?yàn)槠渲饕{(diào)控靶基因中催乳素和生長(zhǎng)激素的活動(dòng)［4，5］，并且已經(jīng)有研究證明乳蛋白基因的表達(dá)受催乳素和腎上腺皮質(zhì)激素的協(xié)同調(diào)控［6，7］，對(duì)奶牛注射生長(zhǎng)激素對(duì)奶牛的產(chǎn)奶量有提高作用［8］。許多研究也證明了STAT5A蛋白與奶牛的產(chǎn)奶性能有著密切的聯(lián)系。Maria Selvaggi在意大利布朗牛STAT5A基因中發(fā)現(xiàn)多態(tài)性，并且CC和CT基因型與產(chǎn)奶量和乳蛋白含量有非常明顯的相關(guān)性。CC基因型比CT基因型產(chǎn)奶量高，CC基因型乳蛋白質(zhì)含量也高于CT基因型，但乳脂含量沒(méi)有明顯的差異［9］。Pawe Brym在荷斯坦和新澤西州牛的內(nèi)含子9發(fā)現(xiàn)一個(gè)SNP（A／G）。對(duì)荷斯坦牛來(lái)說(shuō)，不同的基因型與泌乳性狀之間沒(méi)有明顯的相關(guān)性，但在新澤西牛中，不同基因型與頭胎次和二胎次泌乳量，乳蛋白率和乳脂率之間存在著明顯的相關(guān)性，GG型牛有高產(chǎn)奶量，而AA和AG型牛蛋白質(zhì)產(chǎn)量高［10］。我們前期的研究發(fā)現(xiàn)STAT5A基因在檳榔江水牛上有高度多態(tài)性，為本研究繼續(xù)開(kāi)展多態(tài)性與產(chǎn)奶性能的關(guān)聯(lián)研究奠定了基礎(chǔ)。

2 材料與方法

2.1 試驗(yàn)材料

83頭檳榔江水牛耳組織樣本采自云南省騰沖縣巴福樂(lè)檳榔江水牛良種繁育公司。浸泡于無(wú)水酒精中于－20℃保存?zhèn)溆?。鑒于檳榔江水牛第三胎產(chǎn)奶量最高，本研究采用檳榔江水牛第三胎次產(chǎn)奶量、乳蛋白率和乳脂率數(shù)據(jù)為產(chǎn)奶性狀資料。

2.2 試驗(yàn)方法

2.2.1 基因組DNA的提取 采用傳統(tǒng)的苯酚／氯仿抽提法。

2.2.2 引物合成 根據(jù)NCBI上提供的牛（Bostaurus）的STAT5A全基因組序列（Gene bank登錄號(hào)：NW＿003104496.1）利用primer3在線設(shè)計(jì)3對(duì)引物，引物由上海生工生物工程公司合成。引物序列、PCR產(chǎn)物大小、PCR體系與條件見(jiàn)文獻(xiàn)［11］。

2.2.3 序列分析及數(shù)據(jù)分析 采用DNAStar軟件包對(duì)測(cè)序結(jié)果進(jìn)行分析，采用人工校正。選用適當(dāng)?shù)南拗菩詢?nèi)切酶建立PCR－RFLP檢測(cè)方法。等位基因頻率、基因型頻率及Hardy－Weinberg遺傳平衡分析采用SAS9.0進(jìn)行。

2.2.4 PCR－RFLP 通過(guò)對(duì)引物 P3（5＇－GGCTGGAACTACACCTTCTGG－3＇，5＇－GGGTGTTCTCGTTCTTGAGC－3＇）擴(kuò)增產(chǎn)物片段測(cè)序結(jié)果進(jìn)行分析，結(jié)果發(fā)現(xiàn)在內(nèi)含子15的一個(gè)C／T堿基替換產(chǎn)生了一個(gè)MspI限制性內(nèi)切酶位點(diǎn)，因此本研究建立了PCR－RFLP分子標(biāo)記技術(shù)。酶切反應(yīng)體系為20 μL：10×酶解緩沖液2μL，MspI（20I／μL）0.5μL，4 μL PCR產(chǎn)物，加水至20μL，37℃水浴4h以上。MspI酶切位點(diǎn)由C和G兩個(gè)等位基因控制，經(jīng)酶切得到三種基因型，分別為CC，CG，GG。基因型為CC的片段被切成452bp和330bp兩條帶；雜合子CG為782bp，452bp和330bp三條帶；基因型為GG的片段只有一條782bp的條帶，未被切開(kāi)。

2.2.5 產(chǎn)奶性狀分析 本研究中的產(chǎn)奶量、乳脂率、乳蛋白率均為牛場(chǎng)DHI生產(chǎn)記錄數(shù)據(jù)。根據(jù)最小二乘法建立下列模型：

式中Yij為表型值；μ為群體平均值；Gi為第i種基因型效應(yīng)；eij為隨機(jī)誤差。根據(jù)以上線性模型，利用SAS9.0對(duì)產(chǎn)奶數(shù)據(jù)進(jìn)行分析。

3 結(jié)果與分析

檳榔江水牛耳組織基因組DNA的提取及STAT5A基因PCR擴(kuò)增結(jié)果良好，可以用于測(cè)序和酶切反應(yīng)。為了檢測(cè)檳榔江水牛STAT5A基因多態(tài)性，選擇8個(gè)具有特征特性的水牛個(gè)體的DNA用三對(duì)引物進(jìn)行擴(kuò)增，然后測(cè)序。通過(guò)對(duì)所有測(cè)序結(jié)果進(jìn)行分析，共發(fā)現(xiàn)38個(gè)變異位點(diǎn)。其中引物P1擴(kuò)增產(chǎn)物中發(fā)現(xiàn)24個(gè)變異位點(diǎn)；引物P2擴(kuò)增產(chǎn)物中發(fā)現(xiàn)11個(gè)變異位點(diǎn)；引物P3擴(kuò)增產(chǎn)物中發(fā)現(xiàn)3個(gè)變異位點(diǎn)。所檢測(cè)到的所有變異位點(diǎn)均為單核苷酸變異（SNP），所有變異位點(diǎn)主要以轉(zhuǎn)換、顛換為主，未發(fā)現(xiàn)堿基插入或缺失。其中外顯子8產(chǎn)生2處C924T，C975T突變，氨基酸分析表明，沒(méi)有導(dǎo)致308位氨基酸（脯氨酸）和325位氨基酸（絲氨酸）的變化；外顯子13產(chǎn)生2處G1 482A，C1 548T突變，沒(méi)有導(dǎo)致494位氨基酸（脯氨酸）和516位氨基酸（丙氨酸）的變化。具體結(jié)果見(jiàn)文獻(xiàn)［11］。

3.1 檳榔江水牛STAT 5A基因序列測(cè)定及基因分型

本研究對(duì)所有個(gè)體進(jìn)行了基因分型，選取了STAT5A基因變異中4個(gè)外顯子變異位點(diǎn)（924／E8，975／E8，1482／E13，1548／E13）進(jìn)行分析。各位點(diǎn)測(cè)序結(jié)果見(jiàn)圖1～4。

圖1 924／E8變異位點(diǎn)序列

圖2 975／E8變異位點(diǎn)序列

圖3 1482／E13變異位點(diǎn)序列

圖4 1548／E13變異位點(diǎn)序列

3.2 檳榔江水牛STAT5A基因 MspI標(biāo)記建立

用MspI限制性內(nèi)切酶對(duì)全部83個(gè)樣品進(jìn)行酶切反應(yīng)，得到的結(jié)果如下。圖5為擴(kuò)增產(chǎn)物酶切電泳結(jié)果，表1為對(duì)MspI酶切位點(diǎn)統(tǒng)計(jì)的基因型頻率和基因頻率。

表1 引物P3擴(kuò)增產(chǎn)物PCR－RFLP檢測(cè)到的基因型頻率和等位基因頻率

通過(guò)計(jì)算發(fā)現(xiàn)，CG為優(yōu)勢(shì)基因型，C為優(yōu)勢(shì)基因，等位基因C的頻率明顯高于等位基因G的頻率。

3.3 檳榔江水牛STAT5A基因多態(tài)信息分析

對(duì)4個(gè)外顯子變異及MspI酶切位點(diǎn)進(jìn)行多態(tài)信息分析，利用遺傳學(xué)理論對(duì)每個(gè)位點(diǎn)的多態(tài)信息含量、雜合度和有效等位基因數(shù)進(jìn)行分析并進(jìn)行χ2適合性檢驗(yàn)，結(jié)果見(jiàn)表2。

由統(tǒng)計(jì)結(jié)果可以看出，在924／E8變異位點(diǎn)，群體多態(tài)信息含量（PIC）為0.336，達(dá)到中度多態(tài)性（0.25＜PIC＜0.5）。有效等位基因數(shù)接近觀察到的實(shí)際等位基因數(shù)。在975／E8變異位點(diǎn)中，PIC為0.100，屬于低度多態(tài)（PIC＜0.25），雜合度和有效基因數(shù)的數(shù)值也反映該位點(diǎn)遺傳變異較小。在1482／E13變異位點(diǎn)，PIC為0.346，達(dá)到中度多態(tài)，雜合度和有效等位基因數(shù)的數(shù)值也反映相同的趨勢(shì)。在外顯子1548／E13變異位點(diǎn)，PIC為0.240，屬于低度多態(tài)，遺傳變異較小，雜合度和有效等位基因數(shù)也反映了這一趨勢(shì)。MspI酶切位點(diǎn)的PIC為0.349，達(dá)到中度多態(tài)性。χ2適合性檢驗(yàn)結(jié)果顯示，5個(gè)變異位點(diǎn)差異水平均不顯著（P＞0.05），處于Hardy－Weinberg平衡狀態(tài)。

圖5 引物P3擴(kuò)增產(chǎn)物MspI酶切結(jié)果M：DGL 2000DNA Marker

3.4 檳榔江水牛STAT5A基因多態(tài)性與檳榔江水牛產(chǎn)奶性狀的關(guān)聯(lián)性分析

選取STAT5A基因變異中4個(gè)外顯子變異位點(diǎn)（924／E8，975／E8，1482／E13，1548／E13）和 MspI酶切位點(diǎn)分別與檳榔江水牛第三胎次的產(chǎn)奶量、乳蛋白率和乳脂率進(jìn)行關(guān)聯(lián)性分析，結(jié)果如下：

表2 檳榔江水牛STAT5A基因變異位點(diǎn)的遺傳特性

表3 924／E8變異位點(diǎn)多態(tài)性與檳榔江水牛產(chǎn)奶性狀的關(guān)聯(lián)性分析

表4 975／E8變異位點(diǎn)多態(tài)性與檳榔江水牛產(chǎn)奶性狀的關(guān)聯(lián)性分析

分別將STAT5A基因外顯子8的兩個(gè)變異位點(diǎn)的不同基因型與檳榔江水牛第三胎產(chǎn)奶性狀進(jìn)行關(guān)聯(lián)分析，結(jié)果表明：在924／E8位點(diǎn)，CT基因型個(gè)體的產(chǎn)奶量極顯著高于TT基因型個(gè)體（P＜0.01），三種基因型對(duì)乳蛋白率及乳脂率均沒(méi)有顯著影響（P＞0.05）；在975／E8位點(diǎn)，CT基因型個(gè)體的乳脂率顯著高于CC基因型個(gè)體（P＜0.05），產(chǎn)奶量和乳蛋白率在兩種基因型之間沒(méi)有顯著差異（P＞0.05）。分析結(jié)果見(jiàn)表3和4。

分別將STAT5A基因外顯子13的兩個(gè)變異位點(diǎn)的不同基因型與檳榔江水牛第三胎產(chǎn)奶性狀進(jìn)行關(guān)聯(lián)分析，結(jié)果表明：在1482／E13位點(diǎn)，產(chǎn)奶量，乳蛋白率和乳脂率在兩種基因型之間沒(méi)有顯著差異；在1548／E13位點(diǎn)，產(chǎn)奶量，乳蛋白率和乳脂率在兩種基因型之間也沒(méi)有顯著差異。分析結(jié)果見(jiàn)表5和6。

表5 1482／E13變異位點(diǎn)多態(tài)性與檳榔江水牛產(chǎn)奶性狀的關(guān)聯(lián)性分析

表6 1548／E13變異位點(diǎn)多態(tài)性與檳榔江水牛產(chǎn)奶性狀的關(guān)聯(lián)性分析

表7 STAT5A基因MspI酶切位點(diǎn)基因多態(tài)性與檳榔江水牛產(chǎn)奶性狀的關(guān)聯(lián)性分析

將STAT5A基因的MspI位點(diǎn)與檳榔江水牛的產(chǎn)奶性狀進(jìn)行關(guān)聯(lián)分析，得到的結(jié)果如下：CG基因型個(gè)體產(chǎn)奶量和乳脂率極顯著高于CC基因型個(gè)體（P＜0.01），乳蛋白率在三種基因型個(gè)體之間沒(méi)有顯著差異（P＞0.05）。結(jié)果見(jiàn)表7。

4 討論

4.1 檳榔江水牛STAT5A基因序列顯示高度多態(tài)

檢測(cè)到的38個(gè)變異位點(diǎn)中，發(fā)生在內(nèi)含子的變異有34個(gè)，外顯子變異有4個(gè)。本實(shí)驗(yàn)是首次以檳榔江水牛為樣本對(duì)STAT5A作檢測(cè)，他人未作過(guò)相關(guān)報(bào)道，因而無(wú)從比較。但是發(fā)生在檳榔江水牛STAT5A基因上的變異與有關(guān)報(bào)道發(fā)生在其他牛STAT5A基因上的變異相比較還是不同的。Maria Selvaggi運(yùn)用PCR－RFLP（AvaI酶切）技術(shù)在意大利布朗牛發(fā)現(xiàn)多態(tài)性，它是位于外顯子7第6 853位置的一個(gè)C／T堿基替換［9］，F(xiàn)lisikowski等對(duì)夏洛萊，紅安格斯，利木贊，海福特，西門塔爾及荷斯坦公牛進(jìn)行了STAT5A外顯子7單核苷酸多態(tài)性的研究，結(jié)果表明在STAT5A基因外顯子7存在一個(gè)SSCP位點(diǎn)［12］；Khatib等人在奶牛STAT5A基因外顯子8中發(fā)現(xiàn)一個(gè)SNP位點(diǎn)是G／C堿基替換［13］；McCracken等人［14］在STAT5A 基因內(nèi)含子12中發(fā)現(xiàn)了不同長(zhǎng)度的TG重復(fù)；Flisikowski等人在牛內(nèi)含子15中發(fā)現(xiàn)CCT缺失，而在外顯子16中發(fā)現(xiàn)一個(gè) T／C突變［15，16］；2005年Jang等報(bào)道了在韓國(guó)奶牛STAT5A基因內(nèi)含子15檢測(cè)到1個(gè)SNP變異位點(diǎn)［17］。通過(guò)比較可以看出，STAT5A基因在檳榔江水牛上的變異與在其他牛上檢測(cè)到的變異不同并且遠(yuǎn)遠(yuǎn)高于在其他牛上的變異，這也反映了檳榔江水牛STAT5A基因高度的多態(tài)性。此外Flisikowski等人還在STAT5A基因的啟動(dòng)子區(qū)發(fā)現(xiàn)兩個(gè)SNP位點(diǎn)，分別是－247A／G替換和－226G／A替換［18］。本研究未對(duì)此區(qū)域作檢測(cè)，因而無(wú)法比較。出現(xiàn)這種情況可能與檳榔江水牛所處的地理環(huán)境有關(guān)。騰沖縣位于云南省保山市西南部，緊鄰孟加拉灣和東南亞地區(qū)，西部與緬甸毗鄰，歷史上曾是古西南絲綢之路的要沖，騰沖縣的地理位置決定了檳榔江水牛與亞洲其他水牛種群之間的交流十分頻繁，而LEI等的研究也揭示該地區(qū)屬于沼澤型水牛和河流型水牛的混雜之地［19］，這可能是導(dǎo)致檳榔江水牛的遺傳變異十分豐富的原因。

本研究對(duì)5個(gè)變異位點(diǎn)作了多態(tài)信息含量、有效等位基因數(shù)和遺傳雜合度的分析，其中有3個(gè)位點(diǎn)達(dá)到了中度多態(tài)，一個(gè)位點(diǎn)非常接近中度多態(tài)，只有一個(gè)位點(diǎn)屬于低度多態(tài)，有效等位基因數(shù)和雜合度也反映了相同的趨勢(shì)，這說(shuō)明檳榔江水牛群體遺傳變異程度較高，遺傳多樣性較豐富。χ2檢驗(yàn)結(jié)果顯示所有位點(diǎn)都達(dá)到了 Hardy－Weinberg平衡，這表明在此位點(diǎn)的基因型頻率和基因頻率將穩(wěn)定遺傳，說(shuō)明了檳榔江水牛STAT5A基因在有豐富的遺傳變異的同時(shí)也具有一定的遺傳穩(wěn)定性。

4.2 STAT5A基因多態(tài)性與檳榔江水牛產(chǎn)奶性狀間有顯著關(guān)聯(lián)性

STAT5A基因作為產(chǎn)奶性能的候選基因之一，一直是遺傳學(xué)研究的熱點(diǎn)內(nèi)容。本研究選擇了STAT5A基因的5個(gè)變異位點(diǎn)與檳榔江水牛第三胎產(chǎn)奶量、乳蛋白率、乳脂率等產(chǎn)奶性狀進(jìn)行了關(guān)聯(lián)性分析，以期尋找與數(shù)量性狀位點(diǎn)相關(guān)的遺傳標(biāo)記。分析發(fā)現(xiàn)，5個(gè)變異位點(diǎn)中有3個(gè)位點(diǎn)與產(chǎn)奶性能有顯著性相關(guān)。在975／E8位點(diǎn)，雖然產(chǎn)奶量和乳蛋白率在兩種基因型之間沒(méi)有顯著差異（P＞0.05），但呈現(xiàn)CC＞CT的趨勢(shì)。在MspI酶切位點(diǎn)，CG基因型個(gè)體產(chǎn)奶量和乳脂率顯著高于CC基因型個(gè)體（P＜0.01），乳蛋白率在三種基因型個(gè)體之間沒(méi)有顯著差異（P＞0.05），但也呈現(xiàn)CG＞CC＞GG的趨勢(shì)，由此可以看出CG基因型對(duì)檳榔江水牛的產(chǎn)奶性能起到了正面的效應(yīng)。

本研究結(jié)果與其他關(guān)于STAT5A基因與產(chǎn)奶性能關(guān)聯(lián)的研究結(jié)果一致。Maria Selvaggi運(yùn)用PCR－RFLP（AvaI酶切）技術(shù)在意大利布朗牛上作的STAT5A基因與產(chǎn)奶性能關(guān)聯(lián)性分析結(jié)果發(fā)現(xiàn)，外顯子7第6 853位置的一個(gè)C／T堿基替換，并且CC和CT基因型與產(chǎn)奶量，乳蛋白含量有非常明顯的關(guān)聯(lián)性。CC基因型比CT基因型產(chǎn)奶量高，CC基因型乳蛋白質(zhì)含量也高于CT基因型，但乳脂含量沒(méi)有明顯的差異［9］。Brym等在2004年報(bào)道了STAT5A基因內(nèi)含子9的一個(gè)SNP（A／G），在186頭荷斯坦牛中，不同的基因型與泌乳性狀之間沒(méi)有明顯的相關(guān)性，而在138頭新澤西牛中，不同基因型與頭胎次和二胎次泌乳量，乳蛋白率和乳脂率之間存在著明顯的相關(guān)性，GG型牛有高產(chǎn)奶量，而AA和AG型牛蛋白質(zhì)產(chǎn)量高［10］。何峰等在2007年對(duì)758頭中國(guó)荷斯坦奶牛進(jìn)行了STAT5A基因多態(tài)性檢測(cè)，并將其與5個(gè)產(chǎn)奶性狀進(jìn)行了關(guān)聯(lián)分析。結(jié)果發(fā)現(xiàn)，在內(nèi)含子9的一個(gè)SNP（A／G）對(duì)乳蛋白率有顯著影響（P＜0.05），12 440位T→C和12 550位的CCT插入／缺失，此位點(diǎn)對(duì)產(chǎn)奶量有極顯著影響（P＜0.01），對(duì)乳脂量、乳蛋白量有顯著影響（P＜0.05）［20］。

通過(guò)本研究和國(guó)內(nèi)外其他的研究可以看出，STAT5A基因多態(tài)性在多種牛上檢測(cè)到了其與產(chǎn)奶性能的顯著關(guān)聯(lián)性，因此本研究認(rèn)為STAT5A基因可以作為與產(chǎn)奶性能有關(guān)的候選基因。前人的研究發(fā)現(xiàn)檳榔江水牛產(chǎn)奶性比云南當(dāng)?shù)仫曫B(yǎng)的其他河流型水牛和沼澤性水牛產(chǎn)奶天數(shù)短，產(chǎn)奶量卻高，但個(gè)體間產(chǎn)奶性能差異很大［1］。本研究檢測(cè)了5個(gè)變異位點(diǎn)與產(chǎn)奶性能的關(guān)系，這為提高檳榔江水牛的產(chǎn)奶性能作了很好的鋪墊，特別是與產(chǎn)奶性能有相關(guān)性的變異位點(diǎn)可以作為分子標(biāo)記而應(yīng)用，如逐步固定對(duì)產(chǎn)奶性能有正面效用的基因型，以提高生產(chǎn)性能。同時(shí)這些變異位點(diǎn)可以整合到檳榔江水牛選育方案中，建立產(chǎn)奶性能標(biāo)記輔助選擇。

4.3 檳榔江水牛STAT5A基因MspI酶切位點(diǎn)可作為一個(gè)有效分子標(biāo)記

在MspI酶切位點(diǎn)χ2檢驗(yàn)達(dá)到了Hardy－Weinberg平衡，這說(shuō)明在此位點(diǎn)的基因型頻率和基因頻率將穩(wěn)定遺傳，處于遺傳平衡狀態(tài)。PIC為0.349，達(dá)到中度多態(tài)，這說(shuō)明該遺傳標(biāo)記可以提供比較豐富和合理的遺傳信息。酶切結(jié)果顯示，此酶切位點(diǎn)CG基因型頻率高于CC和GG基因型頻率，說(shuō)明在此位點(diǎn)雜合子是優(yōu)勢(shì)基因型，在特定的條件下雜合子個(gè)體可能比正常純合子個(gè)體更有利于生存和繁殖后代。PCR－RFLP技術(shù)操作簡(jiǎn)便，分型時(shí)間短，特異性強(qiáng)，并且MspI限制性內(nèi)切酶價(jià)格便宜，性能穩(wěn)定，是十分常用的一種內(nèi)切酶。因此在大規(guī)模群體檢測(cè)中此位點(diǎn)可以作為一個(gè)快速簡(jiǎn)便并且經(jīng)濟(jì)的有效分子遺傳標(biāo)記。

［1］ 屈在久，沈雪鷹，楊加用，等.騰沖檳榔江水牛的起源及初步研究［J］.云南畜牧獸醫(yī)，2008，（增刊）：34－37.

［2］ 章純熙，苗永旺，李大林.我國(guó)首例本土河流型水牛－檳榔江水牛的種質(zhì)特征［J］.家畜生態(tài)學(xué)報(bào)，2011，32（6）：39－45.

［3］ 屈在久，李大林，苗永旺.檳榔江水牛種質(zhì)資源調(diào)查與評(píng)價(jià)［J］.云南農(nóng)業(yè)大學(xué)學(xué)報(bào)，2008，23（2）：265－269.

［4］ Wakao H，Gouilleux F，Groener B，et al.Mammary gland factor（MGF）is a novel member of the cytokine－regulated transcription factor gene family and confers the prolactin response［J］.EMBO 13（9）：2182－2191.

［5］ Argetsinger L S，Carter－Su C.Growth hormone signalling mechanisms：involvement of the tyrosine kinase JAK2［J］.Hormone Research，1996，4522－4524.

［6］ Ball R K，F(xiàn)riis R R，Schoenenberger C A，et al.Prolactin regulation ofβ－casein gene expression and of a cytosolic 120KD protein in a cloned mouse mammary epithelial cell line［J］.EMBO，1988，7：2089－2095.

［7］ Lee C S，Kim K，Yu D Y，et al.Pretreatment with glucocorticoid is essential for lactogenic induction of the bovine beta－casein／CAT expression in HC11cells［J］.Endocrine Research，1998，24：65－77.

［8］ Bauman D E，Vernon R G.Effects of exogenous bovine somatotropin on lactation［J］.Annual review of nutrition，1993，13：437－461.

［9］ Maria Selvaggi，Cataldo Dario，Giovanni Normanno，et al.Genetic polymorphism of STAT5Aprotein：relationships with production traits and milk composition in Italian Brown cattle［J］.Journal of Dairy Research，2009，76（4）：441－445.

［10］ Pawe Brym，Stanisaw Kam?ski，Anna Ruceae.New SSCP polymorphism within bovine STAT5Agene and its associations with milk performance traits in Black－and－White and Jersey cattle［J］.Appl Genet，2004，45（4）：445－452.

［11］ 季 敏，余長(zhǎng)林，余選富，等.檳榔江水牛STAT5A基因多態(tài)性及 MspI PCR－RFLP遺傳標(biāo)記的建立［J］.中國(guó)牛業(yè)科學(xué)，2012，38（6）：12－18.

［12］ Flisikowski K，Zwierzchowski L.Single strand conformation polymorphism with exon 7of the bovine STAT5Agene［J］.Animal Science Papers and Reports，2002，2（20）：133－137.

［13］ Khatib H，Monson R L，Schutzkus V，et al.Mutations in the STAT5Agene are associated with embryonic survival and milk composition in cattle［J］.Journal of Dairy Science，2008，91：784－793.

［14］ McCracken J Y，Molenaar A J，Snell R J，et al.A polymorphic TG repeat present within the bovine STAT5Agene［J］.Animal Genetics.1997，28，453－464.

［15］ Flisikowski K，Szymanowska M，Zwierzchowski L.The DNA－binding capacity of genetic variants of the bovine STAT5A transcription factor［J］.Celluar ＆ Molecular Biology Letters，2003，8：831－840.

［16］ Flisikowski K，Zwierzchowski L.Polymerase chain reactionheteroduplex （PCR－HD）polymorphism within the bovine STAT5Agene［J］.Appl Genet，2003，44：185－189.

［17］ Jang G W，Cho K H，Kim T H，et al.Association of Candidate genes with production traits in Korean dairy Proven and Young bulls［J］.Asian－Aust J Anim Sci，2005，18（2）：165－169.

［18］ Flisikowski K，Hiendleder S，Zwierzchowski L.Nucleotide Sequence Variation in the Transcription Factor STAT5Agene 5＇－noncoding Region in Bostaurus and Bos indicus Cattle［J］.Biochemical Genetics，2005，43：459－464.

［19］ Lei C Z，Zhang W，Chen H，et al.Independent maternal origion of Chinese swamp buffalo（Bubalus bubalis）［J］.AnimGenet，2007，38（2）：97－102.

［20］ 何 峰，孫東曉，俞 英，等.荷斯坦奶牛STAT5A基因的SNPs檢測(cè)及其與產(chǎn)奶性狀的關(guān)聯(lián)分析［J］.畜牧獸醫(yī)學(xué)報(bào)，2007，38（4）：326－331.