胡艷 張叁保 張瑜 高小童 江雨航 申玉建 韋英明 蔣欽楊

摘要：【目的】明確轉(zhuǎn)移抑制素基因（KISS1）多態(tài)性與努比亞山羊繁殖性狀（產(chǎn)羔數(shù)、初生重和斷奶重）的關(guān)聯(lián)，為進(jìn)一步展開山羊分子標(biāo)記輔助選擇育種提供參考依據(jù)?！痉椒ā坷肈NA混合池結(jié)合Sanger測序篩選出努比亞山羊KISS1基因SNP位點(diǎn)，再通過MTA-Seq測序?qū)?55只努比亞山羊群體進(jìn)行SNP位點(diǎn)的基因分型，并采用SAS 9.2的一般線性模型（GLM）對KISS1基因SNP位點(diǎn)與努比亞山羊的產(chǎn)羔表型性狀進(jìn)行關(guān)聯(lián)分析?！窘Y(jié)果】在努比亞山羊KISS1基因內(nèi)含子1發(fā)現(xiàn)2個(gè)SNP位點(diǎn)（g.1341600A>G和g.1341674C>G）。其中，g.1341600A>G位點(diǎn)處于中度多態(tài)（0.25<PIC<0.50），且處于Hardy-Weinberg平衡狀態(tài)（PHWE>0.05）;而g.1341674C>G位點(diǎn)處于低度多態(tài)（PIC<0.25），已偏離Hardy-Weinberg平衡狀態(tài)（PHWE<0.05）。g.1341600A>G位點(diǎn)GG和AG基因型個(gè)體的第一胎、第二胎及三胎聯(lián)合產(chǎn)羔數(shù)均高于AA基因型個(gè)體， 且在第一胎達(dá)顯著水平（P<0.05，下同）;g.1341674C>G位點(diǎn)CC基因型個(gè)體則表現(xiàn)為第一胎、第三胎及三胎聯(lián)合產(chǎn)羔數(shù)均高于GG基因型個(gè)體，同樣在第一胎達(dá)顯著水平。g.1341674C>G位點(diǎn)GG基因型個(gè)體的第一胎、第三胎及三胎聯(lián)合羔羊初生重均高于CC基因型個(gè)體，且在第三胎達(dá)顯著水平;g.1341600A>G位點(diǎn)不同基因型與努比亞山羊各胎次羔羊初生重?zé)o顯著關(guān)聯(lián)（P>0.05，下同）。在羔羊斷奶重方面，2個(gè)SNP位點(diǎn)不同基因型與努比亞山羊各胎次斷奶重也無顯著關(guān)聯(lián)。g.1341600A>G位點(diǎn)和g.1341674C>G位點(diǎn)可形成連鎖不平衡單倍型模塊，其中GC、AC、AG單倍型的頻率分別為0.645、0.265和0.090;單倍型模塊與第一胎產(chǎn)羔數(shù)顯著關(guān)聯(lián)，AACC組合型產(chǎn)羔數(shù)顯著低于其他組合型?！窘Y(jié)論】努比亞山羊KISS1基因內(nèi)含子1存在2個(gè)SNP位點(diǎn)（g.1341600A>G和g.1341674C>G）與其繁殖性狀相關(guān)，其中g(shù).1341600A>G位點(diǎn)的A等位基因和g.1341674C>G位點(diǎn)的C等位基因可作為努比亞山羊選育的潛在分子標(biāo)記，且在今后的育種過程中應(yīng)淘汰KISS1基因AACC組合型母羊。

關(guān)鍵詞：努比亞山羊;KISS1基因;SNP位點(diǎn);產(chǎn)羔數(shù);初生重;斷奶重

中圖分類號：S827.89? ? ? ? ? ? ? ? ? ? ? ? ? ?文獻(xiàn)標(biāo)志碼： A 文章編號：2095-1191（2021）10-2880-07

Abstract：【Objective】Clarifying the relationship between transfer inhibin gene （KISS1） polymorphism and reproductive traits （litter size， birth weight and weaning weight） in Nubian goats， so as to provide reference basis for further molecular marker assisted selection breeding of goats. 【Method】 In this study， the SNP sites of KISS1 gene in Nubian goats were screened by DNA mixing pool and Sanger sequencing. Then 355 Nubian goats were genotyped by MTA-Seq sequencing. The association between SNP sites of KISS1 gene and lambing phenotypic traits of Nubian goats was analyzed by SAS 9.2 general linear model （GLM）.【Result】The results showed that two SNP of g.1341600A>G and g.1341674C>G were found in KISS1 gene intron. Population genetic analysis showed that the g.1341600A>G locus was moderately polymorphic （0.25<P<0.50）， and the? test showed that the g.1341600A>G locus was in Hardy-Weinberg equilibrium （PHWE>0.05）， while the g.1341674C>G site was low polymorphism （PIC<0.25）， and? test showed that the g.1341600A>G locus was not in Hardy-Weinberg equilibrium （PHWE<0.05）.Association analysis showed that at g.1341600A>G locus， the first， second and total parity litter size of GG and AG genotypes were higher than those of AA genotypes and reached significant level at the first parity litter size（P<0.05， the same below）. At g.1341674C>G locus， the first ， third and total parity litter size of CC genotypes were higher than those of GG genotypes and reached significant level at the first parity litter size.At g.1341674C>G locus， the first， third parity and total parity birth weight of GG genotypes were higher than those of CC genotypes， and reached significant level at the third parity birth weight.There was no significant difference among different genotypes of Nubian goats in diverse parities birth weight at g.1341600A>G（P>0.05， the same below）. There was no significant difference among different genotypes of Nubian goats in diverse parities weaning weight at g.1341600A>G and g.1341674C>G.The results of linkage disequilibrium and haplotype analysis showed that g.1341600 C>G locus and g.1341674C>G locus could form linkage disequilibrium module， and the frequencies of GC， AC and AG were 0.645， 0.265 and 0.090， respectively. The correlation analysis between haplotype combination and reproductive traitsshowed that there was a significant correlation between haplotype combination and the first parity litter size， and the litter size of AACC was significantly lower than those of other genotypes.【Conclusion】There are two SNP loci in KISS1 gene intron （g.1341600A>G and g.1341674C>G）， which are related to their reproductive traits in Nubian goats. The A allele of g.1341600A>G and the C allele of g.1341674C>G can be used as potential molecular markers for Nubian goat breeding， and the AACC combination ewes of KISS1 gene should be eliminated in the breeding process in the future.

Key words：Nubian goat; KISS1gene; SNP locus; litter size; birth weight; weaning weight

Foundation item：Guangxi Key Research and Development Plan （Guike AB18221067，Guike AB19245010）;Guangxi Cattle and Sheep Industry Innovation Team Construction Project （nycytxgxcxtd-2021-09）

0 引言

【研究意義】山羊繁殖性狀遺傳力較低，采用常規(guī)育種方法很難在短時(shí)間內(nèi)實(shí)現(xiàn)高效率的品種改良。轉(zhuǎn)移抑制素（Kisspeptins，KISS1）被認(rèn)為是哺乳動物生殖分子機(jī)制的關(guān)鍵調(diào)控因子，KISS1基因編碼產(chǎn)物Kisspeptin是G蛋白偶聯(lián)受體54（GPR54）內(nèi)源性受體，能與GPR54特異性結(jié)合組成KISS1/GPR54系統(tǒng)（Ohtaki et al.，2001;彭曉利，2018），參與下丘腦—垂體—性腺（HPG）軸的調(diào)控，通過促性腺激素釋放激素（GnRH）調(diào)節(jié)促卵泡生成素（FSH）和促黃體生成素（LH）的分泌水平，進(jìn)而調(diào)控哺乳動物的繁殖性狀（秦津等，2011;劉海斌等，2019;徐唯嘉和王家宏，2019）。因此，分析KISS1基因在哺乳動物季節(jié)性繁殖中的作用及其多態(tài)性，對有效篩選出繁殖性狀相關(guān)分子標(biāo)記應(yīng)用于輔助育種具有重要意義?！厩叭搜芯窟M(jìn)展】目前，針對KISS1基因的研究主要集中在動物機(jī)體初情期啟動（金慧慧等，2015）及其季節(jié)性繁殖（安雪姣等，2019）等方面。KISS1基因在機(jī)體內(nèi)廣泛表達(dá)，以在胎盤和睪丸等性腺組織中的表達(dá)量較高（Tena-Sempere，2006）。有研究表明，KISS1基因突變或缺失后人類或小鼠的性腺功能退化，而導(dǎo)致其繁殖性能喪失（de Roux et al.，2003;de Tassigny，2007）。Cao等（2010）研究發(fā)現(xiàn)，濟(jì)寧青山羊KISS1基因內(nèi)含子1存在3處突變（G269、G454和T505A），內(nèi)含子2則存在18 bp的插入缺失片段，均對產(chǎn)羔數(shù)有顯著影響。Hou等（2011）以西農(nóng)薩能、關(guān)中山羊和波爾山羊?yàn)檠芯繉ο?，發(fā)現(xiàn)KISS1基因內(nèi)含子2存在T2643C突變及8 bp的缺失片段（AGTTCCCC），且均與產(chǎn)羔數(shù)相關(guān)。An等（2013）在分析西農(nóng)薩能、關(guān)中山羊和波爾山羊的KISS1基因多態(tài)性時(shí)發(fā)現(xiàn)，有4個(gè)多態(tài)性位點(diǎn)（SNP）與其產(chǎn)仔數(shù)相關(guān)。李隱俠等（2014）在130只蘇淮山羊群體發(fā)現(xiàn)2個(gè)KISS1基因SNP位點(diǎn)，其中g(shù).2270C>;T位點(diǎn)與產(chǎn)羔數(shù)相關(guān)。Maritra等（2014）對印度山羊KISS1基因多態(tài)性進(jìn)行檢測分析，結(jié)果發(fā)現(xiàn)有3個(gè)SNP位點(diǎn)與產(chǎn)羔數(shù)相關(guān)。金慧慧（2015）將豬源KISS1基因轉(zhuǎn)移至小鼠體內(nèi)后，可在小鼠下丘腦內(nèi)表達(dá)，且轉(zhuǎn)基因小鼠窩產(chǎn)仔數(shù)高于野生型小鼠，故推測KISS1基因與小鼠產(chǎn)仔數(shù)有關(guān)。陳祥等（2018）研究黔北麻羊KISS1基因多態(tài)性時(shí)發(fā)現(xiàn)，在其內(nèi)含子1檢測到3個(gè)SNP位點(diǎn)，其中G277C位點(diǎn)和T486A位點(diǎn)與產(chǎn)羔數(shù)存在關(guān)聯(lián)性?！颈狙芯壳腥朦c(diǎn)】努比亞山羊具有性情溫馴、羊膻味低、肉質(zhì)口感好及繁殖力強(qiáng)等優(yōu)點(diǎn)，可作為雜交父本用于改良本地山羊品種。廣西是山羊傳統(tǒng)養(yǎng)殖區(qū)，年出欄肉羊超過200萬只，努比亞山羊作為引進(jìn)品種已在全區(qū)廣泛養(yǎng)殖并用于本地山羊的品種改良（蔣小剛等，2011;顏燦健等，2018），但至今鮮見有關(guān)努比亞山羊繁殖性狀的研究報(bào)道?！緮M解決的關(guān)鍵問題】以努比亞山羊?yàn)檠芯繉ο螅肈NA混合池法篩選KISS1基因SNP位點(diǎn)，并通過多重PCR測序進(jìn)行SNP分型檢測，旨在明確KISS1基因多態(tài)性與努比亞山羊繁殖性狀（產(chǎn)羔數(shù)、初生重和斷奶重）的關(guān)聯(lián)，為進(jìn)一步展開山羊分子標(biāo)記輔助選擇育種提供參考依據(jù)。

1 材料與方法

1. 1 試驗(yàn)材料

選擇有產(chǎn)羔記錄（包括產(chǎn)羔數(shù)、初生重和斷奶重）的努比亞山羊355只，其中，227只來自廣西扶綏廣羊農(nóng)牧有限公司，128只來自廣西武鳴種羊場，飼養(yǎng)管理?xiàng)l件一致。從頸靜脈采集5 mL血液樣品，放入含抗凝劑的真空采血管中，低溫運(yùn)回實(shí)驗(yàn)室，-20 ℃保存?zhèn)溆?。全血基因組DNA提取試劑盒及其Marker購自天根生化科技（北京）有限公司;PCR擴(kuò)增試劑盒購自南京諾唯贊生物科技股份有限公司。主要儀器設(shè)備有核酸凝膠成像儀（美國Bio-Rad公司）、PCR儀（Applide Biosystems，美國Bio-Rad公司）及核酸濃縮儀（美國Scientist公司）等。

1. 2 DNA提取

采用全血基因組DNA提取試劑盒從努比亞山羊血液提取基因組DNA，紫外分光光度計(jì)測定260和280 nm處的OD以檢測DNA質(zhì)量，并利用0.8%瓊脂糖凝膠電泳檢測其完整性。然后將部分DNA濃度調(diào)至50 ng/μL，用于構(gòu)建DNA混合池，剩余DNA置于-20 ℃冰箱保存?zhèn)溆谩?/p>

1. 3 引物設(shè)計(jì)與合成

參照An等（2013）、李鵬程（2017）的方法設(shè)計(jì)PCR擴(kuò)增引物，預(yù)期產(chǎn)物長度為377 bp，委托生工生物工程（上海）股份有限公司合成。其中，P1引物對（F：5'-CCCGCTGTAACTAGAGAAAG-3';R：5'-CA TCCAGGGTGAGTGATACT-3'）用于DNA混合池SNP檢測;P2引物對（F：5'-GTCTCATCCAGGGTGAGT GATA-3';R：5'-GCTCTTTCTGGGTAAGGGAGGAT-3'）P2用于MTA-Seq基因分型。

1. 4 PCR擴(kuò)增

PCR反應(yīng)體系及其擴(kuò)增程序參照鄒輝等（2019）的方法。

1. 5 Sanger測序及基因分型

從355管的DNA中隨機(jī)選取150個(gè)樣組成3個(gè)混合池，50個(gè)樣/池，分別從每組50管DNA樣品中取出0.5 μL，使每組樣品分別在3個(gè)離心管中混合均勻，獲得3個(gè)混合池用于擴(kuò)增特異性片段，經(jīng)1.5%瓊脂糖凝膠電泳檢測達(dá)到預(yù)期要求后，送至生工生物工程（上海）股份有限公司進(jìn)行Sanger測序。經(jīng)序列比對分析和文獻(xiàn)查閱總結(jié)，將DNA混合池檢測到的SNP位點(diǎn)與已報(bào)道的產(chǎn)羔性狀相關(guān)基因SNP位點(diǎn)進(jìn)行匯總，根據(jù)多態(tài)性位點(diǎn)上、下游100 bp堿基序列設(shè)計(jì)多重PCR組合引物，檢測所有SNP位點(diǎn)，然后送至生工生物工程（上海）股份有限公司進(jìn)行MTA-Seq基因分型。

1. 6 統(tǒng)計(jì)分析

使用Excel 2013統(tǒng)計(jì)努比亞山羊KISS1基因SNP位點(diǎn)的基因型頻率、等位基因頻率、多態(tài)信息含量（PIC）、群體純合度（Ho）、群體雜合度（He）及有效等位基因數(shù)（Ne），并進(jìn)行Hardy-Weinberg平衡檢測。利用SAS 9.2的一般線性模型（GLM）對KISS1基因SNP位點(diǎn)與產(chǎn)羔表型性狀進(jìn)行關(guān)聯(lián)分析，GLM模型如下：

Yijl = μ+Gi+Lj+Pl+eijl

式中，Yijl為產(chǎn)羔數(shù)，[μ]為總體均值，Gi為基因型效應(yīng)，Lj為場次效應(yīng)，Pl為胎次效應(yīng)，eijl為隨機(jī)殘差效應(yīng)。由于山羊繁殖性狀（產(chǎn)羔數(shù)、初生重和斷奶重）是由微效多基因決定，故選用Haploview 4.2對2個(gè)SNP位點(diǎn)與繁殖性狀間進(jìn)行連鎖不平衡和單倍型分析。

2 結(jié)果與分析

2. 1 努比亞山羊KISS1基因SNP位點(diǎn)篩選結(jié)果

提取的努比亞山羊血液基因組DNA經(jīng)檢測，其濃度和純度均符合要求，條帶清晰單一，無任何雜帶（圖1），可供后續(xù)試驗(yàn)使用。以提取的努比亞山羊血液基因組DNA為模板，PCR擴(kuò)增獲得的KISS1基因片段經(jīng)1.5%瓊脂糖凝膠電泳檢測，其目的片段大?。?77 bp）與預(yù)期結(jié)果相符（圖2）。

將DNA混合池PCR擴(kuò)增產(chǎn)物送至生工生物工程（上海）股份有限公司進(jìn)行測序，經(jīng)BLAST比對分析，結(jié)果發(fā)現(xiàn)在努比亞山羊KISS1基因內(nèi)含子1存在1處突變位點(diǎn)g.1341674C>G（圖3）。此外，通過多重PCR測序在KISS 1基因內(nèi)含子1找出另一個(gè)SNP位點(diǎn)g.1341600A>G。

2. 2 努比亞山羊KISS1基因SNP位點(diǎn)的群體遺傳學(xué)分析結(jié)果

群體遺傳學(xué)分析結(jié)果（表1）表明，努比亞山羊KISS1基因在g.1341600A>G位點(diǎn)的A、G等位基因頻率分別為0.35和0.65，其中G為優(yōu)勢等位基因;AA、AG、GG基因型頻率分別為0.12、0.47和0.41，其中AG基因型為優(yōu)勢基因型。努比亞山羊KISS1基因在g.1341674C>G位點(diǎn)的C、G等位基因頻率分別為0.83和0.17，其中C為優(yōu)勢等位基因;CC、GG基因型頻率分別為0.83和0.17，以CC基因型為優(yōu)勢基因型;該位點(diǎn)未見CG基因型，可能是由于CG基因型卵泡發(fā)育受阻或卵巢內(nèi)無初級卵泡卵母細(xì)胞發(fā)育，且其卵泡周圍顆粒細(xì)胞發(fā)育異常或突變雜合子導(dǎo)致子宮和卵巢發(fā)育不全，從而導(dǎo)致機(jī)體不育。此外，g.1341600 A>G位點(diǎn)處于中度多態(tài)（0.25<PIC<0.50），χ2適應(yīng)性檢驗(yàn)結(jié)果表明該位點(diǎn)處于Hardy-Weinberg平衡狀態(tài)（PHWE>0.05）;而g.1341674C>G位點(diǎn)處于低度多態(tài)（PIC<0.25），χ2適應(yīng)性檢驗(yàn)結(jié)果表明該位點(diǎn)已偏離Hardy-Weinberg平衡狀態(tài)（PHWE<0.05）。

2. 3 KISS1基因多態(tài)性與努比亞山羊繁殖性狀間的關(guān)聯(lián)性

在產(chǎn)羔數(shù)（表2）方面，g.1341600A>G位點(diǎn)GG和AG基因型努比亞山羊個(gè)體的第一胎、第二胎及三胎聯(lián)合產(chǎn)羔數(shù)均高于AA基因型個(gè)體，且在第一胎達(dá)顯著水平（P<0.05，下同）;g.1341674C>G位點(diǎn)CC基因型個(gè)體則表現(xiàn)為第一胎、第三胎及三胎聯(lián)合產(chǎn)羔數(shù)均高于GG基因型個(gè)體，同樣在第一胎達(dá)顯著水平。在羔羊初生重（表3）方面，g.1341674C>G位點(diǎn)GG基因型個(gè)體的第一胎、第三胎及三胎聯(lián)合羔羊初生重均高于CC基因個(gè)體型，其中第三胎達(dá)顯著水平;g.1341600A>G位點(diǎn)不同基因型與努比亞山羊各胎次羔羊初生重?zé)o顯著關(guān)聯(lián)（P>0.05，下同）。在羔羊斷奶重（表4）方面，2個(gè)SNP位點(diǎn)不同基因型與努比亞山羊各胎次斷奶重也無顯著關(guān)聯(lián)。

2. 4 努比亞山羊KISS1基因SNP位點(diǎn)連鎖不平衡及單倍型分析結(jié)果

SNP位點(diǎn)連鎖不平衡及單倍型分析結(jié)果表明，g.1341600A>G位點(diǎn)和g.1341674C>G位點(diǎn)可形成1個(gè)連鎖不平衡單倍型模塊（圖4），其中GC、AC、AG單倍型的頻率分別為0.645、0.265和0.090。將單倍型模塊與努比亞山羊繁殖性狀（產(chǎn)羔數(shù)、初生重和斷奶重）進(jìn)行關(guān)聯(lián)分析，結(jié)果發(fā)現(xiàn)與第一胎產(chǎn)羔數(shù)顯著關(guān)聯(lián)，AACC組合型產(chǎn)羔數(shù)顯著低于其他組合型（表5）;而與初生重和斷奶重的關(guān)聯(lián)均未達(dá)顯著水平。

3 討論

內(nèi)含子雖然不編碼蛋白，無法對蛋白的結(jié)構(gòu)產(chǎn)生影響，但其變異對轉(zhuǎn)錄和剪切過程均產(chǎn)生影響，通過控制基因表達(dá)來調(diào)控動物表型（田志龍等，2019）。至今，有關(guān)內(nèi)含子變異影響動物繁殖性狀的研究已有較多報(bào)道。E1-Tarabany等（2017）研究發(fā)現(xiàn)，在KISS1基因內(nèi)含子1的第121位堿基處存在T→A現(xiàn)象，并證實(shí)該突變位點(diǎn)與大馬士革羊產(chǎn)羔數(shù)相關(guān)。陳祥等（2018）研究表明，黔北麻羊KISS1基因內(nèi)含子1的2個(gè)SNP位點(diǎn)致使其mRNA表達(dá)存在差異，進(jìn)而導(dǎo)致產(chǎn)羔數(shù)顯著增加。Hui等（2020）在陜北絨山羊DNMT3B基因內(nèi)含子22上發(fā)現(xiàn)1個(gè)11 bp的插入突變，并證實(shí)該突變位點(diǎn)與第一胎產(chǎn)羔數(shù)顯著相關(guān)。本研究在努比亞山羊KISS1基因內(nèi)含子1發(fā)現(xiàn)2個(gè)SNP位點(diǎn)（g.1341600A>G和g.1341674C>G），χ[2]適應(yīng)性檢驗(yàn)結(jié)果表明，g.1341600A>G位點(diǎn)處于Hardy-Weinberg平衡狀態(tài)，即該位點(diǎn)在適應(yīng)性方面具有一定優(yōu)勢;而g.1341674C>G位點(diǎn)已偏離Hardy-Weinberg平衡狀態(tài)，推測該位點(diǎn)是日常飼養(yǎng)過程中人為干預(yù)力度較大的結(jié)果，也有可能是試驗(yàn)樣本數(shù)量太少所致。

本研究的關(guān)聯(lián)分析結(jié)果表明，g.1341600A>G位點(diǎn)GG和AG基因型個(gè)體的第一胎、第二胎及三胎聯(lián)合產(chǎn)羔數(shù)均高于AA基因型個(gè)體，且在第一胎達(dá)顯著水平。值得注意的是，g.1341600A>G位點(diǎn)突變純合型產(chǎn)羔數(shù)大多高于野生純合型，從而推測優(yōu)勢等位基因（G）在一定程度上提高了努比亞山羊的產(chǎn)羔能力，與陳祥等（2018）的研究結(jié)論一致，也進(jìn)一步證實(shí)g.1341600A>G位點(diǎn)與努比亞山羊繁殖性狀顯著相關(guān)。g.1341674C>G位點(diǎn)CC基因型個(gè)體的產(chǎn)羔數(shù)在第一胎、第三胎及三胎聯(lián)合均高于GG基因型個(gè)體，同樣在第一胎達(dá)顯著水平，但g.1341674C>G位點(diǎn)野生純合型產(chǎn)羔數(shù)大多高于突變純合型，故推測優(yōu)勢等位基因（C）在一定程度上提高了努比亞山羊的產(chǎn)羔能力，與An等（2013）的研究結(jié)果一致，說明g.1341674A>G位點(diǎn)與努比亞山羊繁殖性狀顯著相關(guān)。

初生重能放映羊羔胚胎時(shí)期的營養(yǎng)水平，代表羊羔的先天素質(zhì);斷奶重則體現(xiàn)出羊羔在哺乳期的營養(yǎng)水平、生長發(fā)育情況及今后的成長潛力?？梢?，除產(chǎn)羔數(shù)外，初生重和斷奶重也是影響山羊養(yǎng)殖業(yè)收益的重要繁殖性狀。杜建文等（2020）研究表明，母羊每胎產(chǎn)羔數(shù)越多，其羔羊的初生重越輕。本研究中，g.1341674C>G位點(diǎn)GG基因型個(gè)體的羔羊初生重在第一胎、第三胎及三胎聯(lián)合均高于CC基因型個(gè)體，且在第三胎達(dá)顯著水平，可能是由于GG基因型產(chǎn)羔數(shù)少，胎兒對母體營養(yǎng)物質(zhì)的競爭小，平均每只羔羊從母體內(nèi)吸收的營養(yǎng)較多，與張磊等（2020）的研究結(jié)論一致;而g.1341600A>G位點(diǎn)不同基因型與努比亞山羊各胎次羔羊初生重?zé)o顯著關(guān)聯(lián)。g.1341600A>G位點(diǎn)和g.1341674C>G位點(diǎn)的不同基因型與努比亞山羊各胎次斷奶重也無顯著關(guān)聯(lián)，即這2個(gè)SNP位點(diǎn)不適合作為努比亞山羊斷奶重改良的分子標(biāo)記。此外，單倍型模塊與努比亞山羊繁殖性狀（產(chǎn)羔數(shù)、初生重和斷奶重）的關(guān)聯(lián)分析發(fā)現(xiàn)，單倍型模塊與第一胎產(chǎn)羔數(shù)顯著關(guān)聯(lián)，AACC組合型產(chǎn)羔數(shù)顯著低于其他組合型，因此在今后的育種過程中應(yīng)淘汰AACC組合型母羊。

4 結(jié)論

努比亞山羊KISS1基因肉含子1存在2個(gè)SNP位點(diǎn)（g.1341600A>G和g.1341674C>G）與其繁殖性狀相關(guān)，其中g(shù).1341600A>G位點(diǎn)的A等位基因和g.1341674C>G位點(diǎn)的C等位基因可作為努比亞山羊選育的潛在分子標(biāo)記，且在今后的育種過程中應(yīng)淘汰KISS1基因AACC組合型母羊。

參考文獻(xiàn)：

安雪姣，趙生國，潘章源，文禹粱，田志龍，蔡原，儲明星. 2019. 綿羊KISS1基因組織表達(dá)及其多態(tài)性與產(chǎn)羔數(shù)關(guān)聯(lián)分析[J]. 中國畜牧雜志，55（6）：45-50. [An X J，Zhao S G，Pan Z Y，Wen Y L，Tian Z L，Cai Y，Chu M X. 2019. Tissue expression and polymorphism of sheep KISS1 gene and their association with litter size[J]. Chinese Journal of Animal Science，55（6）：45-50.] doi：10.19556/j. 0258-7033.2019-06-045.

陳祥，李鵬程，孫振梅，許厚強(qiáng)，陳偉. 2018. 黔北麻羊KISS1基因多態(tài)性及組織差異表達(dá)研究[J]. 農(nóng)業(yè)生物技術(shù)學(xué)報(bào)，26（7）：1213-1218. [Chen X，Li P C，Sun Z M，Xu H Q，Chen W. 2018. Study on polymorphism and tissue diffe-rential expression of KISS1 gene in Qianbei Ma goats （Capra hircus）[J]. Journal of Agricultural Biotechnology，26（7）：1213-1218.] doi：10.3969/j.issn.1674-7968.2018. 07.013.

杜建文，譚世新，田月珍，付雪峰，徐新明，李紅燕，許建標(biāo)，胥磊，楊雪梅，黃錫霞，田可川. 2020. 影響不同品種羔羊初生重的因素分析[J]. 家畜生態(tài)學(xué)報(bào)，41（8）：33-36. [Du J W，Tan S X，Tian Y Z，F(xiàn)u X F，Xu X M，Li H Y，Xu J B，Xu L，Yang X M，Huang X X，Tian K C. 2020. Factors influencing the birth weight of different breeds of lambs[J]. Acta Ecologae Animalis Domastici，41（8）：33-36.] doi：10.3969/j.issn.1673-1182.2020.08.006.

蔣小剛，韋錦益，磨勤，周恒，蔣玉秀. 2011. 努比亞山羊引種飼養(yǎng)觀察報(bào)告[J]. 廣西畜牧獸醫(yī)，27（5）：300-301. [Jiang X G，Wei J Y，Mo Q，Zhou H，Jiang Y X. 2011. Observation report on introduction and breeding of Nubian goat in 2011[J]. Guangxi Journal of Animal Husbandry & Veterinary Medicine，27（5）：300-301.] doi：10.3969/j.issn.1002-5235.2011.05.023.

金慧慧，劉秋月，狄冉，胡文萍，王翔宇，儲明星. 2015. Kiss1基因調(diào)控哺乳動物繁殖性能的研究進(jìn)展[J]. 家畜生態(tài)學(xué)報(bào)，36（8）：1-6. [Jin H H，Liu Q Y，Di R，Hu W P，Wang X Y，Chu M X. 2015. Research progress on reproduction regulation of Kiss1 gene in mammals[J]. Acta Ecologae Animalis Domastici，36（8）：1-6.] doi：10.3969/j.issn.1673-1182.2015.08.001.

金慧慧. 2015. 轉(zhuǎn)豬FSHR和KISS1基因小鼠繁殖性能分析[D]. 重慶：西南大學(xué). [Jin H H. 2015. Evaluation the reproductive performance of porcine KISS1 and FSHR transgenic mice[D]. Chongqing：Southwest University.]

李鵬程. 2017. ESR、KISS-1、PGR基因mRNA表達(dá)及多態(tài)性與黔北麻羊產(chǎn)羔性狀相關(guān)性分析[D]. 貴陽：貴州大學(xué). [Li P C. 2017. Correlation analysis between polymorphism of ESR，KISS-1，PGR and litter size of Qianbei Ma goat and their mRNA expressio level testing[D]. Guiyang：Guizhou University.]

李隱俠，張俊，錢勇，李靜心，孟春花，王慧利，鐘聲，曹少先. 2014. Kiss-1基因多態(tài)性與蘇淮山羊產(chǎn)羔數(shù)的關(guān)聯(lián)分析[J]. 畜牧獸醫(yī)學(xué)報(bào)，45（12）：1917-1923. [Li Y X，Zhang J，Qian Y，Li J X，Meng C H，Wang H L，Zhong S，Cao S X. 2014. Kiss-1 gene polymorphisms and its association with litter sizes in Suhuai goats[J]. Acta Veterinaria et Zootechnica Sinica，45（12）：1917-1923.] doi：10.11843/j.issn.0366-6964.2014.12.002.

劉海斌，吳占福，張立永，官麗輝，陳寶江，胡滿. 2019. 光照強(qiáng)度對Kisspeptin及其受體GPR54在雞腦表達(dá)的影響[J]. 中國獸醫(yī)學(xué)報(bào)，39（3）：497-506. [Liu H B，Wu Z F，Zhang L Y，Guan L H，Chen B J，Hu M. 2019. Effect of light intensity on expression of Kisspeptin and GPR54 in brain of hens[J]. Chinese Journal of Veterinary Science，39（3）：497-506.] doi：10.16303/j.cnki.1005-4545.2019. 03.21.

彭曉利. 2018. KISS1/GPR54系統(tǒng)在山羊睪丸間質(zhì)細(xì)胞自分泌的功能研究[D]. 重慶：西南大學(xué). [Peng X L. 2018. The self-secretory function of KISS1/GPR54 in goat leying cells[D]. Chongqing：Southwest University.]

秦津，李瑞明，凌澤繼，謝瑩雪，宋小白，楊柄壯，許典新. 2011. 促性腺激素釋放激素在不同品種水牛卵巢中的免疫組化定位[J]. 南方農(nóng)業(yè)學(xué)報(bào)，42（1）：90-93. [Qin J，Li R M，Ling Z J，Xie Y X，Song X B，Yang B Z，Xu D X. 2011. Immuno-histochemical localization of gonadotropin-releasing hormone in ovaries of different buffalo varieties[J]. Journal of Southern Agriculture，42（1）：90-93.] doi：10.3969/j.issn.2095-1191.2011.01.021.

田志龍，湯繼順，孫慶，王玉琴，張效生，張金龍，儲明星. 2019. 綿羊SMAD1基因組織表達(dá)及其多態(tài)性與產(chǎn)羔數(shù)關(guān)聯(lián)分析[J]. 中國農(nóng)業(yè)科學(xué)，52（4）：755-766. [Tian Z L，Tang J S，Sun Q，Wang Y Q，Zhang X S，Zhang J L，Chu M X. 2019. Tissue expression and polymorphism of sheep SMAD1 gene and their association with litter size[J]. Scientia Agricultura Sinica，52（4）：755-766.] doi：10.3864/j.issn.0578-1752.2019.04.015.

徐唯嘉，王家宏. 2019. KISS-1基因相關(guān)研究進(jìn)展[J]. 系統(tǒng)醫(yī)學(xué)，4（1）：192-195. [Xu W J，Wang J H. 2019. Research progress of KISS-1 gene related study[J]. Systems Medicine，4（1）：192-195.] doi：10.3969/j.issn.2096-1782. 2019.01.072.

顏燦健，張偉，黃蘭珍，朱志旅，陳子良，楊膺白. 2018. 廣西山羊養(yǎng)殖模式及效果探討[J]. 廣西畜牧獸醫(yī)，34（2）：96-98. [Yan C J，Zhang W，Huang L Z，Zhu Z L，Chen Z L，Yang Y B. 2018. Discussion on goat breeding mode and effect in Guangxi[J]. Guangxi Journal of Animal Husbandry & Veterinary Medicine，34（2）：96-98.] doi：10. 3969/j.issn.1002-5235.2018.02.015.

張磊，喬為盛，白前前，宋宇軒，周占琴. 2020. 湖羊胎產(chǎn)羔數(shù)對羔羊生長發(fā)育的影響[J]. 家畜生態(tài)學(xué)報(bào)，41（4）：43-46. [Zhang L，Qiao W S，Bai Q Q，Song Y X，Zhou Z Q. 2020. Effects of lambing number on the growth and development of lambs in Hu sheep[J]. Acta Ecologae Animalis Domastici，41（4）：43-46.] doi：10.3969/j.issn.1673-1182.2020.04.008.

鄒輝，瞿秋紅，夏琴，崔悅悅，江雨航，韋玲靜，嚴(yán)雪瑜，蔣欽楊，韋英明. 2019. 努比亞山羊BMPR-IB基因多態(tài)性與其產(chǎn)羔性狀的關(guān)聯(lián)分析[J]. 南方農(nóng)業(yè)學(xué)報(bào)，50（4）：860-866. [Zou H，Qu Q H，Xia Q，Cui Y Y，Jiang Y H，Wei L J，Yan X Y，Jiang Q Y，Wei Y M. 2019. Correlation ana-lysis between Nubian goats BMPR-IB gene polymorphism and its lambing traits[J]. Journal of Southern Agriculture，50（4）：860-866.] doi：10.3969/j.issn.2095-1191. 2019.04.25

An X P，Han P，Hou J X，Zhao H B，Yan Y，Ma T，F(xiàn)ang F，Meng F X，Song Y X，Wang J G，Cao B Y. 2013. Mole-cular cloning and characterization of KISS1 promoter and effect of KISS1 gene mutations on litter size in the goat[J]. Genetics and Molecular Research，12（4）：4308-4316. doi：10.4238/2013.February.28.13.

Cao G L，Chu M X，F(xiàn)ang L，Di R，F(xiàn)eng T，Li N. 2010. Analysis on DNA sequence of KiSS-1 gene and its association with litter size in goats[J]. Molecular Biology Reports，37（8）：3921-3929. doi：10.1007/s11033-010-0049-7.

de Roux N，Genin E，Carel J C，Matsuda F，Chaussain J L，Milgrom E. 2003. Hypogonadotropic hypogonadism due to loss of function of the KiSS1-derived peptide receptor GPR54[J]. Proceedings of the National Academy of Scien-ces of the United States of America，100（19）：10972-10976. doi：10.1073/pnas.1834399100.

de Tassigny X A，F(xiàn)agg L A，Dixon J P C，Day K，Leitch H G，Hendrick A G，Zahn D，F(xiàn)ranceschini I，Caraty A，Carlton M B L，Aparicio S A J R，Colledge W H. 2007. Hypogonadotropic hypogonadism in mice lacking a functional Kiss1 gene[J]. Proceedings of the National Academy of Sciences of the United States of America，104（25）：10714-10719. doi：10.1073/pnas.0704114104.

El-Tarabany M S，Zaglool A W，El-Tarabany A A，Awad A. 2017. Association analysis of polymorphism in KiSS1 gene with reproductive traits in goats[J]. Animal Reproduction Science，180：92-99. doi：10.1016/j.anireprosci. 2017.03.006.

Hou J X，An X P，Wang J G，Song Y X，Cui Y H，Wang Y F，Chen Q J，Cao B Y. 2011. New genetic polymorphisms of KiSS-1 gene and their association with litter size in goats[J]. Small Ruminant Research，96（2-3）：106-110. doi：10.1016/j.smallrumres.2010.11.013.

Hui Y Q，Zhang Y H，Wang K，Pan C Y，Chen H，Qu L，Song X Y，Lan X Y. 2020. Goat DNMT3B：An indel mutation detection，association analysis with litter size and mRNA expression in gonads[J]. Theriogenology，147：108-115. doi：10.1016/j.theriogenology.2020.02.025.

Maitra A，Sharma R，Ahlawat S，Tantia M S，Roy M，Prakash V. 2014. Association analysis of polymorphisms in caprine KiSS1 gene with reproductive traits[J]. Animal Reproduction Science，151（1-2）：71-77. doi：10.1016/j.anireprosci.2014.09.013.

Ohtaki T，Shintani Y，Honda S，Matsumoto H，Hori A，Kanehashi K，Terao Y，Kumano S，Takatsu Y，Masuda Y，Ishibashi Y，Watanabe T，Asada M，Yamada T，Suenaga M，Kitada C，Usuki S，Kurokawa T，Onda H，Nishimura O，F(xiàn)ujino M. 2001. Metastasis suppressor gene KiSS-1 encodes peptide ligand of a G-protein-coupled receptor[J]. Nature，411（6837）：613-617. doi：10.1038/35079135.

Tena-Sempere M. 2006. The roles of kisspeptins and G protein-coupled receptor-54 in pubertal development[J]. Current Opinion in Pediatrics，18（4）：442-447. doi：10. 1097/01.mop.0000236396.79580.cc.

（責(zé)任編輯 蘭宗寶）