牛志剛,王 珊,常 璐,史洪才

(新疆畜牧科學院 生物技術研究所/農(nóng)業(yè)農(nóng)村部 草食家畜繁育生物技術重點開放實驗室,新疆 烏魯木齊 830026)

多浪羊是新疆一個以肉用為主、肉脂兼用型地方綿羊品種,主要分布在新疆南部地區(qū)的葉爾羌河流域、塔里木盆地西南邊緣地區(qū)。多浪羊具有體型高大、生長發(fā)育快、早熟、產(chǎn)肉和繁殖率高、適應性強、遺傳性穩(wěn)定、飼料報酬高等特點。多浪羊母羊四季均可發(fā)情,以4～5月和9～11月發(fā)情較多,發(fā)情周期18 d,發(fā)情持續(xù)期24～48 h,妊娠期150 d,平均產(chǎn)羔率130%左右。多浪羊是新疆南疆地區(qū)非常優(yōu)良的肉羊品種,因此提高多浪羊的繁殖性能對新疆南疆畜牧業(yè)發(fā)展具有重大意義。陳曉軍等[1]、鐘發(fā)剛等[2]、史洪才等[3]對多浪羊的FecB基因進行了多態(tài)性分析,發(fā)現(xiàn)新疆多浪羊存在兩種基因型(B+、++ ),其中B+型在多浪羊體內(nèi)為優(yōu)勢基因型,說明多浪羊的多胎性能與FecB基因密切相關。本課題組近年來利用FecB基因分子標記輔助育種培育多浪羊多羔新品系,已取得了顯著的效果,多羔育種核心群的平均產(chǎn)羔率達到149%。在此基礎上,本課題組利用全基因組重測序技術,通過比較產(chǎn)單羔組與產(chǎn)雙羔組的差異,挖掘了影響多浪羊繁殖性狀的相關基因和分子標記,最終富集分析篩選了NCOA1、INHBA、ACP1、ING5、ARNT、KLHL1等一些候選基因。

核受體輔激活因子1(Nuclear receptor coactivator 1, NCOA1),又被命名為類固醇受體輔活化因子-1(steroid receptor coactivator-1, SRC-1),屬于對配體有依賴性的轉(zhuǎn)錄因子,是p160類固醇受體輔助活化因子家族成員中的一員[4]。Razeto A等[5]提出了NCOA1是生殖激素“開關”的觀點,認為生殖激素長時間打開是導致哺乳動物繁殖性能增加的原因。NCOA1蛋白的作用是通過與結(jié)合到DNA上的核受體相互作用來增強轉(zhuǎn)錄活性,NCOA1 蛋白可以調(diào)節(jié)與繁殖性狀有關基因的轉(zhuǎn)錄,從而影響這些基因的表達,進而影響哺乳動物的繁殖性能。NCOA1基因和NCOA3基因協(xié)同作用可以調(diào)節(jié)胎盤的形態(tài)以及胚胎的存活率[6]。NCOA1基因在不同組織(下丘腦、垂體、卵巢和子宮)中均有表達,并且對妊娠期起到重要的調(diào)節(jié)作用[7]。我們以新疆多浪羊為研究對象,對NCOA1基因的1個內(nèi)顯子和1個外顯子序列進行克隆,分析了Ncoa1基因不同基因型SNP位點的多態(tài)性,以期探討NCOA1基因與多浪羊產(chǎn)羔性狀的相關性,為多浪羊多胎新品系的選育奠定理論基礎。

1 材料與方法

1.1 試驗材料

在新疆喀什麥蓋提縣的多浪羊養(yǎng)殖場采集產(chǎn)單羔母羊(112只),以及產(chǎn)雙羔和三羔母羊(56只)；從頸靜脈采血,所采血樣量均為5 mL/只,用肝素鈉抗凝,于-20 ℃凍存。

1.2 主要試劑

TaqDNA聚合酶、dNTP購于北京天根生物技術有限公司,限制性內(nèi)切酶AluⅠ、DNA Marker等購自BBI公司。

1.3 DNA的提取

用北京全式金生物技術有限公司生產(chǎn)的Easy Pure Blood Genomic DNA Kit試劑盒提取DNA。

1.4 引物設計

根據(jù)GenBank上公布的綿羊NCOA1基因序列(登陸號為NC-019460.2),利用引物設計軟件Pimer 6.0進行引物設計,對NCOA1基因的1個外顯子、1個內(nèi)含子設計2對引物,引物由上海生工生物工程技術服務有限公司合成,引物信息見表1。

表1 設計的引物信息

1.5 PCR擴增及產(chǎn)物檢測

利用設計的引物對多浪羊的基因組DNA進行PCR擴增。PCR反應體系為50 μL:2 ×Buffer(含Mg2+2.5 U/μL)25 μL、上下游引物各1.0 μL、50 ng/μL的DNA 2 μL、2.5 mmol/L dNTPs 4 μL、TaqDNA聚合酶1 μL，用滅菌雙蒸水補至50 μL。PCR反應程序為:94 ℃預變性5 min;94 ℃變性30 s,62～55 ℃退火30 s,72 ℃延伸30 s,共35個循環(huán);72 ℃延伸7 min;4 ℃保存。用2%瓊脂糖凝膠電泳檢測PCR產(chǎn)物。將上述PCR產(chǎn)物進行純化回收后,送至上海生工生物有限公司進行測序。

1.6 數(shù)據(jù)分析

用DNAStar軟件對NCOA1基因外顯子的測序結(jié)果進行序列比對和分析,計算NCOA1基因的SNP位點在新疆多浪羊群體中的等位基因頻率、雜合度(heterozygosity, He)、有效等位基因數(shù)(number of effective allele, Ne)、群體遺傳純合度(homozygosity, Ho)、多態(tài)信息含量(polymorphism information content, PIC)和基因型頻率;用SPSS 19.0軟件中的最小二乘法分析個體不同基因型與多浪羊產(chǎn)羔數(shù)的關聯(lián)效應。

2 結(jié)果與分析

2.1 NCOA1基因的PCR擴增結(jié)果

對Ncoa12和Ncoa17引物進行PCR擴增,對PCR產(chǎn)物分別用2%瓊脂糖凝膠電泳進行檢測,結(jié)果發(fā)現(xiàn),擴增片段大小與預期一致,片段長度分別為461和416 bp,特異性好，且電泳條帶清晰明亮(圖1和圖2),可進行下一步檢測。

2.2 NCOA1基因測序結(jié)果

對PCR產(chǎn)物進行測序,結(jié)果在NCOA1基因序列中內(nèi)含子12發(fā)現(xiàn)了g.166 397 A>C突變位點,檢測到兩種基因型，分別為AA和AB基因型(圖3)。在外顯子17的g.187 501 G>A突變位點,檢測到3種基因型，分別為BB、AB和AA基因型。序列分析結(jié)果(圖4)表明,在外顯子17的g.187 501 G>A突變位點發(fā)生了錯義突變,由原編碼的絲氨酸(Ser)改變?yōu)楣弱０匪?Gln)。

圖1 NCOA1基因內(nèi)含子12的PCR擴增電泳圖

2.3 基因頻率和基因型頻率

對NCOA1基因多態(tài)位點的基因頻率、基因型頻率進行統(tǒng)計分析，結(jié)果(表2)顯示：在多浪羊群體中,g.166 397 A>C位點的優(yōu)勢基因型為AA基因型,基因型頻率為0.936；優(yōu)勢基因為A等位基因,基因頻率為0.968。在多浪羊群體中, g.187 501 G>A位點的優(yōu)勢基因型為BB基因型,基因型頻率為0.500；優(yōu)勢基因為B等位基因,基因頻率為0.643。

圖3 g.166 397 A>C的測序結(jié)果

圖4 g.187 501 G>A的測序結(jié)果

2.4 NCOA1基因SNP位點基因座的遺傳多態(tài)性分析

采用軟件進行卡方適合性檢驗，結(jié)果表明:新疆多浪羊g.166 397 A>C和g.187 501 G>A位點均處于Hardy-Weinberg平衡狀態(tài),說明多浪羊NCOA1基因突變位點的基因型頻率沒有受到選擇、突變或遷移等因素的影響。對多浪羊群體在NCOA1基因的2個SNP位點的純合度、雜合度、有效等位基因數(shù)以及多態(tài)信息含量進行統(tǒng)計分析,計算結(jié)果(表3)表明,多浪羊群體中NCOA1基因的g.166 397 A>C和g.187 501 G>A位點的多態(tài)信息含量分別屬于低度多態(tài)(PIC<0.25)和中度多態(tài)(0.25

2.5 多浪羊群體NCOA1基因的SNP位點與產(chǎn)羔數(shù)的相關性分析

由表4可知：多浪羊NCOA1基因的g.166 397 A>C突變位點的AA基因型個體產(chǎn)羔數(shù)顯著高于BB基因型個體的(P<0.05);多浪羊NCOA1基因的g.187 501 G>A突變位點的BB基因型個體的產(chǎn)羔數(shù)顯著高于AA和AB基因型個體的(P<0.05),AA、AB基因型個體之間差異不顯著。

表2 基因頻率和基因型頻率

表3 多浪羊群體多態(tài)位點的多態(tài)信息含量、純合度、雜合度和等位基因數(shù)

注:PIC>0.50為高度多態(tài)；0.25

3 討論

研究發(fā)現(xiàn)核受體輔激活蛋白1(NCOA1)基因是影響豬繁殖性狀的一個候選基因并定位于豬3號染色體上。人的NCOA1基因總共有183.4 kb,含有22個外顯子,其cDNA長度為4721 bp。豬的NCOA1基因總共有97681 bp,含有19個外顯子,其cDNA長度為4451 bp。綿羊的NCOA1基因的編碼區(qū)全長為218477 bp,共有9種可變剪切體,含有21或22個外顯子。NCOA1蛋白通過與結(jié)合到DNA上的雌激素受體和雄激素受體相互作用來增強轉(zhuǎn)錄活性,刺激特定的雌激素和雄激素反應基因的轉(zhuǎn)錄并調(diào)節(jié)隨后的生理反應,作用于動物的生長發(fā)育、體內(nèi)微生態(tài)平衡以及繁殖等生理過程[8-10]?；糇噪p等[11]對紅眼白水貂的NCOA1基因多態(tài)性及其與繁殖性狀的相關性進行了分析，發(fā)現(xiàn)在NCOA1基因上存在1個多態(tài)位點,為第6外顯子處g.151 536T>C突變,且AB基因型個體的總產(chǎn)仔數(shù)和產(chǎn)活仔數(shù)均極顯著高于AA基因型個體的(P<0.01)。胡慧艷等[12]在英系和美系大白豬NCOA1基因多態(tài)位點均檢測到AA、AB兩種基因型。武子寅等[13]在邵伯雞母系NCOA1多態(tài)性與產(chǎn)蛋性能的相關研究中發(fā)現(xiàn),該基因在第10外顯子多態(tài)位點,AA型的產(chǎn)蛋數(shù)極顯著高于BB型的,且顯著高于AB型的。李濤等[14]對大白母豬的NCOA1基因進行了多態(tài)性與繁殖性能的關聯(lián)分析,發(fā)現(xiàn)經(jīng)產(chǎn)母豬NCOA1基因中AB型個體的總產(chǎn)仔數(shù)顯著高于AA型個體的。李永輝研究發(fā)現(xiàn)長白豬與大白豬均存在3種基因型,且在不同基因型之間總產(chǎn)仔數(shù)與產(chǎn)活仔數(shù)均差異極顯著[15],說明NCOA1可能是影響豬產(chǎn)仔數(shù)的一個候選基因。NCOA1基因是否調(diào)控綿羊的繁殖性能,尚未見相關報道。本研究將NCOA1基因作為影響多浪羊產(chǎn)羔性狀的候選基因,利用PCR直接測序技術分析了NCOA1基因的多態(tài)性及其與多浪羊繁殖性狀的相關性，發(fā)現(xiàn)在NCOA1內(nèi)含子12處的g.166 397 A>C突變位點產(chǎn)生了AA、AB兩種基因型,AA基因型為優(yōu)勢基因型，且在AA與AB基因型個體間產(chǎn)羔數(shù)存在顯著差異；在NCOA1外顯子17處的g.187 501 G>A突變位點產(chǎn)生了AA、AB、BB三種基因型,BB基因型為優(yōu)勢基因型,在該突變點發(fā)生了錯義突變,使原編碼的絲氨酸(Ser)改變?yōu)楣弱０匪?Gln)，且在AA與BB基因型、AB與BB基因型個體間產(chǎn)羔數(shù)均存在顯著差異，說明NCOA1基因?qū)Χ嗬搜虍a(chǎn)羔性狀的效應是顯著的,表明NCOA1基因可以作為多浪羊產(chǎn)羔性狀的主要候選基因。關于新疆多浪羊的繁殖性能仍需進行深入研究,比如可以將NCOA1基因作為分子遺傳標記,對新疆多浪羊進行標記輔助選種選配,從而提高多浪羊的繁殖性能，這將對新疆南疆畜牧業(yè)的發(fā)展具有重大意義。

表4 NCOA1基因SNP位點與多浪羊產(chǎn)羔數(shù)的相關性

注:在不同基因型之間進行比較,數(shù)據(jù)后有任一字母不同者為差異顯著(P<0.05)。