文章來(lái)源：Identification of Laying-Related SNP Markers in Geese Using RAD Sequencing
閱讀者：樊云鵬

摘要

產(chǎn)蛋性能是鵝生產(chǎn)的一個(gè)重要經(jīng)濟(jì)性狀。由于產(chǎn)蛋性能的遺傳力較低愤惰，因此制定一種標(biāo)記輔助選擇(MAS)策略具有重要意義造成。確定與目標(biāo)性狀相關(guān)的序列變異是量化MAS的前提英妓，但目前對(duì)鵝基因組的研究較少赃阀，這嚴(yán)重阻礙了鵝產(chǎn)蛋性狀遺傳標(biāo)記的識(shí)別丑瞧。最近開(kāi)發(fā)的限制性位點(diǎn)相關(guān)DNA(RAD)測(cè)序是識(shí)別大規(guī)模單核苷酸多態(tài)性(SNP)和降低基因組復(fù)雜性的一種可能方法甫恩，而不需要參考基因組信息逆济。在本研究中，我們開(kāi)發(fā)了一種用于檢測(cè)鵝產(chǎn)蛋相關(guān)SNP的混合RAD測(cè)序策略磺箕。構(gòu)建了兩個(gè)DNA庫(kù)奖慌，每個(gè)庫(kù)由10個(gè)具有高估計(jì)育種值(HEBV)和低估計(jì)育種值(LEBV)的個(gè)體的基因組DNA組成。共從42,291,356個(gè)序列中獲得139,013個(gè)SNP松靡，其中18,771,943個(gè)序列為L(zhǎng)EBV序列简僧，23,519,413序列為HEBV序列。在LEBV和HEBV兩群體中雕欺，通過(guò)個(gè)體AS-PCR基因分型進(jìn)一步驗(yàn)證了在兩個(gè)DNA池中具有不同等位基因頻率的55個(gè)SNP岛马。55個(gè)SNP中有10個(gè)在這兩個(gè)群體中表現(xiàn)出明顯的等位基因分布。這10個(gè)SNP在492只鵝群中進(jìn)一步進(jìn)行了基因分型屠列，以驗(yàn)證其與產(chǎn)蛋量的相關(guān)性啦逆。結(jié)果表明，10個(gè)SNP中有8個(gè)與產(chǎn)蛋量有關(guān)脸哀。此外蹦浦，線性回歸分析表明，SNP記錄-111407撞蜂、106975和112359涉及到影響產(chǎn)蛋率的多基因網(wǎng)絡(luò)盲镶。我們使用IPCR擴(kuò)展在候選RAD標(biāo)記的兩側(cè)未知區(qū)域。對(duì)獲得的序列進(jìn)行BLAST檢索蝌诡，得到鴨和雞的同源基因溉贿。克隆5個(gè)新的基因浦旱，這些基因含有候選的與產(chǎn)蛋相關(guān)的SNP宇色，分別是膜相關(guān)鳥(niǎo)苷酸激酶(MAGI-1)、KIAA 1462颁湖、Rho GTPase激活蛋白21(ARHGAP 21)宣蠕、酰基-CoA合成酶家族成員2(ACSF 2)甥捺、胃動(dòng)蛋白2(ASTN 2)抢蚀。綜上所述，我們的數(shù)據(jù)表明镰禾，8個(gè)SNP和5個(gè)基因可能是鵝產(chǎn)量數(shù)標(biāo)記輔助選擇的候選標(biāo)記或靶點(diǎn)皿曲。

引言

鵝的繁殖力不穩(wěn)定唱逢，孵化性能差，易受遺傳屋休、營(yíng)養(yǎng)坞古、環(huán)境、疾病等多種因素的影響劫樟。由于生殖遺傳力較低痪枫，傳統(tǒng)的選擇方法難以改善繁殖性狀。標(biāo)記輔助選擇(MAS)是改善低遺傳力性狀的有效途徑叠艳。然而听怕，探索性狀連鎖序列變異或功能成因需要發(fā)展MAS技術(shù)。單核苷酸多態(tài)性(SNP)是最豐富的遺傳標(biāo)記類型虑绵，它的遺傳穩(wěn)定性為研究基因組區(qū)的遺傳奠定了基礎(chǔ)尿瞭。然而，目前還沒(méi)有關(guān)于鵝基因組序列的數(shù)據(jù)翅睛，這極大地阻礙了該物種在分子水平上對(duì)經(jīng)濟(jì)性狀的研究声搁。

候選基因方法是識(shí)別與重要經(jīng)濟(jì)性狀相關(guān)的遺傳標(biāo)記的常用方法。Chen等人(2012年)在萬(wàn)江白鵝的催乳素(PRL)內(nèi)含子2中發(fā)現(xiàn)了30多個(gè)SNP捕发，在催乳素受體(PRLR)第10外顯子中發(fā)現(xiàn)了5個(gè)SNP疏旨。這些多態(tài)性與產(chǎn)蛋率顯著相關(guān)。趙等人(2011)發(fā)現(xiàn)兩個(gè)促性腺激素釋放激素(GnRH)和PRL的SNP分別與烏龍鵝的繁殖性狀有關(guān)扎酷。Zhang等人(2013)在紫鵝的不同階段證明了黃體生成素(LH)檐涝、PRL及其受體的基因表達(dá)，丁等人(2006)用抑制消減雜交(SSH)技術(shù)鑒定了產(chǎn)蛋鵝肝臟中許多差異表達(dá)基因法挨。這些基因包括卵黃生成素I谁榜、apoVLDL-Ⅱ、乙醇胺激酶凡纳、G蛋白γ-5亞基和白亮氨酸-tRNA合成酶窃植。最近，郭等人(2011年)采用了類似的方法荐糜，在產(chǎn)蛋期和產(chǎn)蛋期找到了幾個(gè)差異表達(dá)的基因巷怜，包括PRLR、雌激素受體1和抗穆勒激素受體II暴氏。

新一代高通量DNA測(cè)序技術(shù)加速了研究動(dòng)物基因組研究的速度延塑。該技術(shù)在全基因組測(cè)序、靶再測(cè)序答渔，轉(zhuǎn)錄組測(cè)序中得到了廣泛的應(yīng)用关带。最近，Xu等人(2013)利用短讀取測(cè)序技術(shù)(Illuma)鑒定了在產(chǎn)蛋鵝和雛鵝卵巢組織中294個(gè)上調(diào)和278個(gè)下調(diào)基因的572個(gè)差異表達(dá)基因研儒。不幸的是豫缨，由此產(chǎn)生的轉(zhuǎn)錄組僅提供有限的編碼區(qū)域的限制性位點(diǎn)信息，與非編碼區(qū)域相比端朵，這些區(qū)域的核苷酸多樣性要低得多好芭。

限制性位點(diǎn)相關(guān)DNA(RAD)測(cè)序是一種新發(fā)展起來(lái)的快速大規(guī)模SNP發(fā)現(xiàn)方法，可以有效地降低基因組的復(fù)雜性冲呢。它是一種經(jīng)濟(jì)有效的SNP發(fā)現(xiàn)和基因分型方法舍败。它允許較小的研究小組，或研究尚未擁有參考基因組的生物體的小組敬拓，進(jìn)行“基因組范圍的研究”邻薯。RAD測(cè)序法已成功地應(yīng)用于多種生物體，其中包括：格皮乘凸、鮭魚(yú)厕诡、歐亞海貍、喉和虹鱒魚(yú)营勤、鱘魚(yú)和菜籽灵嫌。

在本研究中，我們應(yīng)用RAD測(cè)序葛作，在鵝基因組中發(fā)現(xiàn)了新的SNP寿羞。通過(guò)比較兩個(gè)具有最低估計(jì)育種值(LEBV)和最高估計(jì)育種值(HEBV)的DNA池之間的等位基因頻率來(lái)選擇用于選擇產(chǎn)卵性能的SNP。使用等位基因特異性PCR(AS-PCR)進(jìn)行個(gè)體基因分型赂蠢，首次在LEBV和HEBV群體中證實(shí)了候選SNP-連鎖模式绪穆，并在492只鵝群體中進(jìn)一步驗(yàn)證∈瘢克隆了含有與產(chǎn)蛋相關(guān)的SNP的新基因玖院。

材料和方法

道德聲明

所有實(shí)驗(yàn)均經(jīng)南京農(nóng)業(yè)大學(xué)動(dòng)物保護(hù)和使用委員會(huì)審核批準(zhǔn)，并按照“實(shí)驗(yàn)動(dòng)物管理?xiàng)l例”(中國(guó)第岖，1988年)進(jìn)行司恳。所有的努力都是為了盡量減少采血過(guò)程中的不適。

動(dòng)物和樣品制備

本研究以江蘇麗華畜牧有限公司養(yǎng)殖場(chǎng)的492只揚(yáng)州鵝為研究對(duì)象绍傲。在試驗(yàn)過(guò)程中扔傅，鵝盡可能被喂以添加青草或水草的稻谷。這種飼料是在白天提供的烫饼，那時(shí)鵝被放生到屋外的一個(gè)空地上猎塞。在整個(gè)研究過(guò)程中，鵝暴露在自然光照和溫度下杠纵。將產(chǎn)蛋鵝分別飼養(yǎng)在不同的籠里荠耽，以記錄整個(gè)產(chǎn)蛋期的總產(chǎn)蛋量。用含液泡的肝素鈉從翼靜脈采集血樣比藻。

產(chǎn)蛋性能與分組

biao1.png

RAD文庫(kù)的制備及序列分析

使用全血DNA試劑盒(OmegaBio-TEK拍谐，Doris烛缔，USA）遵循制造商的說(shuō)明從血液中提取全基因組DNA。使用Thermo Scientific Nanodrop 2000分光光度計(jì)評(píng)估每個(gè)個(gè)體樣品的DNA濃度轩拨。所有DNA樣品的終濃度調(diào)整為100 ng/ul践瓷。A260/280和A260/230比值均在標(biāo)準(zhǔn)范圍內(nèi)。每10個(gè)個(gè)體混合等量的基因組DNA準(zhǔn)備LEBV和HEBV的兩個(gè)DNA池亡蓉。用限制性內(nèi)切酶EcoRI消化基因組DNA晕翠。共構(gòu)建了2個(gè)復(fù)用的測(cè)序文庫(kù)，其中每個(gè)DNA樣本被分配一個(gè)獨(dú)特的核苷酸序列標(biāo)識(shí)符 (MID)進(jìn)行條形碼編碼砍濒。用IlluminaHiSeq 2000進(jìn)行單端(101-BP)測(cè)序崖面。

序列分析與產(chǎn)蛋相關(guān)突變檢測(cè)

從3’端將原始序列修剪到90個(gè)核苷酸，保證了97.5%以上核苷酸的質(zhì)量值在Q30以上(等于0.1%的測(cè)序誤差)梯影。通過(guò)序列相似性巫员，利用USTACKS為每個(gè)RAD位點(diǎn)產(chǎn)生獨(dú)特的候選等位基因，將修整后的讀數(shù)聚類成讀取標(biāo)簽（以下簡(jiǎn)稱RAD-標(biāo)簽）甲棍。對(duì)于自然種群简识，在該步驟中允許兩個(gè)最大堿基對(duì)錯(cuò)配。然后感猛，使用RAD-標(biāo)簽在SNP調(diào)用的默認(rèn)參數(shù)下七扰，使用Ustack將其折疊到集群中。對(duì)每個(gè)SNP陪白，LEBV和HEBV池的等位基因頻率差異進(jìn)行比較颈走。選擇兩個(gè)池間等位基因分布顯著不同的SNP作為群體進(jìn)一步驗(yàn)證的候選位點(diǎn)。

鵝群體產(chǎn)蛋相關(guān)突變的驗(yàn)證

在LEBV和HEBV組中咱士，共選擇了55個(gè)SNP進(jìn)行進(jìn)一步的個(gè)體基因分型立由。不同等位基因分布的SNP在LEBV和HEBV中的分布在492只鵝的種群中進(jìn)行了驗(yàn)證。采用AS-PCR技術(shù)對(duì)群體進(jìn)行基因分型序厉。為了提高PCR擴(kuò)增的特異性和對(duì)等位基因的可靠識(shí)別锐膜，在3’端的第三堿基處引入了一個(gè)額外的錯(cuò)配堿基對(duì)。用Primer Premier 5軟件(美國(guó)加利福尼亞州Palo Alto弛房，Premier Bio-ft)道盏，按照Liu和Hayashi的方法設(shè)計(jì)了AS-PCR引物。

引物和PCR產(chǎn)物的長(zhǎng)度在S1表中顯示。在20μL反應(yīng)中利用兩個(gè)特異性引物進(jìn)行基因分型荷逞，其中包括約50 ng模板DNA媒咳、5μL 2X PCR Taq酶(加拿大)、1μL特異引物和普通引物 (10μmol)(中國(guó)深圳)种远。擴(kuò)增條件為：94°C預(yù)變性3 min涩澡，32次擴(kuò)增(94℃ 30 s，45°C-72°C 30 s院促，72°C 30 s)，72°C擴(kuò)增5 min斧抱。用3.0%瓊脂糖凝膠電泳分離PCR產(chǎn)物常拓。

2.png

基于產(chǎn)卵相關(guān)SNP的新基因克隆

3.png

巢式PCR擴(kuò)增RAD標(biāo)記兩側(cè)未知序列。巢式聚合酶鏈反應(yīng)(Nest PCR)在50 ml中酸舍，50ng制備的DNA机错，每種引物2μ1(10μmol)和25μl LATaq酶(TAKARA，大連)父腕。擴(kuò)增條件為：98°C預(yù)變性30 s弱匪，32次擴(kuò)增(98°C 10 s，45°C-72°C 30 s，72°C 4 min)萧诫，72°C擴(kuò)增7 min斥难。第一輪PCR后，用雙蒸餾水稀釋PCR產(chǎn)物1：100帘饶。1微升稀釋液作為第二輪擴(kuò)增的模板哑诊。引物W(外引物對(duì))和N(內(nèi)引物對(duì))分別用于第一次和第二次PCR擴(kuò)增。在1.5% Tris/硼酸/EDTA(TBE)瓊脂糖凝膠上進(jìn)行PCR反應(yīng)及刻。PCR條帶在紫外光下被去除镀裤，并按照供應(yīng)商的建議使用凝膠提取試劑盒(Omega Bio-Tek，Doraville缴饭，USA)進(jìn)行純化暑劝。根據(jù)制造商的方案，使用TA克隆試劑盒(TransgenBiotech颗搂，北京担猛，中國(guó))將純化的DNA片段直接連接到噬菌粒TA載體(pasy-T3質(zhì)粒)中，然后轉(zhuǎn)化為感受態(tài)細(xì)胞(轉(zhuǎn)基因生物技術(shù))丢氢。在含氨芐西林50 mg/ml的LB瓊脂上進(jìn)行轉(zhuǎn)化傅联。選擇菌落取樣，在2.0 Eppendorf管中懸浮于1ml LB培養(yǎng)基中疚察，在37°C下生長(zhǎng)16小時(shí)蒸走。將靶DNA測(cè)序(GeneWiz，蘇州貌嫡，中國(guó))载碌。使用DNAman軟件包進(jìn)行多序列比對(duì)(版本8.0；Lynnon Bio-soft衅枫，加拿大魁北克)嫁艇。使用BlastX（http://www.ncbi.nlm.nih.）進(jìn)行數(shù)據(jù)庫(kù)搜索GOV/BLAST）。.

統(tǒng)計(jì)和數(shù)據(jù)分析

采用獨(dú)立卡方檢驗(yàn)弦撩，檢測(cè)LEBV與HEBV DNA池RAD-標(biāo)記等位基因頻率的差異步咪。對(duì)于發(fā)現(xiàn)與產(chǎn)蛋相關(guān)的SNP，用Bonferroni校正估計(jì)了5%的I類總錯(cuò)誤率的顯著性閾值益楼，給出了

4.png

其中α_Bon是Bonferroni調(diào)整后的P值猾漫，α是未校正的P值，n_Information是SNP的個(gè)數(shù)感凤。Fisher’s精確檢驗(yàn)用統(tǒng)計(jì)語(yǔ)言R版本2.11.1比較LEBV和HEBV的等位基因頻率悯周。個(gè)體的產(chǎn)蛋期估計(jì)育種值(EBV)計(jì)算如下：

5.png

其中

6.png

用PowerMarker V3.25計(jì)算基因型頻率、等位基因頻率长酗、基因多樣性溪北、雜合性、多態(tài)信息含量(PIC)和哈代-溫伯格平衡定律擬合優(yōu)度的卡方檢驗(yàn) 夺脾。所有數(shù)據(jù)均以平均±SD表示之拨。

單因素方差分析(SPSS for Windows，20.0版咧叭；IBM-SPSS蚀乔，芝加哥，IL)比較不同基因型的平均蛋數(shù)佳簸。采用Duncan的多重范圍檢驗(yàn)(SPSS for Windows乙墙，Version20.0)評(píng)估其顯著性颖变。所有顯著性水平為p<0.05的單性狀組合均被納入進(jìn)一步的多重標(biāo)記分析生均。采用線性回歸方法(SPSS for Windows，20.0版)腥刹。采用線性回歸方法马胧，分析多標(biāo)記組合和兩種數(shù)量性狀模式(加性模式：PAa(PAA P_{Aa≈(PAA + Paa)/2})/2)和顯性模式(PAa_≈P_AA或P_aa)。

結(jié)果

RAD測(cè)序

7.png

與產(chǎn)蛋相關(guān)的SNP的發(fā)現(xiàn)

對(duì)RAD測(cè)序的138,675個(gè)SNP進(jìn)行Chisquare檢驗(yàn)筐骇，分析LEBV與HEBV池等位基因頻率的差異债鸡。經(jīng)Bonferroni調(diào)節(jié)后，467個(gè)SNP顯著(p<3.69×10^?7)铛纬。采用等位基因特異性PCR方法對(duì)所有LEBV和HEBV鵝進(jìn)行個(gè)體基因分型厌均。用該方法可獲得55個(gè)SNP的穩(wěn)定基因型(S1表)。進(jìn)一步的個(gè)體基因分型結(jié)果表明告唆，55個(gè)單核苷酸多態(tài)性中有10個(gè)等位基因頻率在LEBV和HEBV隊(duì)列中有顯著性差異 (p<0.00024~4.19×10^?8) (表3)棺弊。

8.png

實(shí)驗(yàn)鵝群體產(chǎn)蛋期相關(guān)單核苷酸多態(tài)性的驗(yàn)證

通過(guò)AS-PCR對(duì)492只鵝群體進(jìn)行了基因分型 (圖1)。利用Powermark V3.25軟件對(duì)每個(gè)SNP進(jìn)行遺傳多樣性分析擒悬。如表4所示模她，10個(gè)SNP的基因多樣性(HE)、雜合性(H0)和多態(tài)性信息內(nèi)容(PIC)分別為0.4394-0.4991懂牧、0.0830-0.533和0.3161-0.3746侈净。記錄135849的SNP具有最高的基因多樣性、雜合度和PIC僧凤。PIC常用作分子標(biāo)記遺傳多態(tài)性的指標(biāo)畜侦。在本研究中，10個(gè)SNP的PIC值范圍為0.25-0.5躯保，表明這些SNP表現(xiàn)出中間水平的多態(tài)性旋膳。7個(gè)SNP顯示與Hardy-Weinberg平衡(HWE)的顯著偏離(p<0.05)，而其它三個(gè)SNP途事，包括記錄-135849验懊、記錄-88247和記錄-135057均在HWE中(p>0.05)。

9.png

10.png

如表5所示尸变，記錄-106975的GG和GA基因型產(chǎn)蛋量顯著高于AA基因型的（p<0.01）义图。GG基因型與GA基因型在產(chǎn)蛋量上無(wú)顯著差異(p>0.05)。

11.png

SNP記錄-134172具有AA基因型鵝產(chǎn)蛋量明顯高于AT基因型鵝(p<0.01)召烂，TT基因型與AA和AT基因型無(wú)顯著性差異(p>0.05)碱工。

記錄-112359，TT基因型的產(chǎn)蛋量顯著高于GG基因型(p<0.01)骑晶。TG基因型與TT基因型和GG基因型相比痛垛，產(chǎn)蛋量差異無(wú)顯著性(p>0.05)。

記錄-106582桶蛔，AA和CA基因型的產(chǎn)蛋量顯著高于GG基因型匙头，但CA基因型與AA基因型之間差異無(wú)顯著性(p>0.05)。

記錄-111407仔雷，AA型和TA型鵝產(chǎn)蛋量顯著高于TT型鵝(p<0.01)蹂析。AA和TA基因型的產(chǎn)蛋量差異不顯著(p>0.05)舔示。

在記錄-88247中，AA基因型的產(chǎn)蛋量顯著高于GG基因型(p<0.01)电抚。AA基因型產(chǎn)蛋量顯著高于GA基因型(p<0.05)惕稻。

對(duì)于記錄135057，AA和GA基因型的產(chǎn)蛋量顯著高于GG基因型(p<0.01)蝙叛。

在記錄-130775中俺祠，AA和AG基因型的產(chǎn)蛋率高于GG基因型(p<0.01)。

記錄-135849和記錄-130652的基因型與產(chǎn)蛋量之間無(wú)顯著相關(guān)性(P>0.05)借帘。結(jié)果表明蜘渣，這8個(gè)單核苷酸多態(tài)性(紀(jì)錄-106975，記錄-134172肺然，記錄-112359蔫缸，記錄-106582，記錄-111407际起，記錄-88247拾碌，記錄-135057，記錄-130652街望，記錄-130775)與產(chǎn)蛋性狀顯著相關(guān)(p<0.01)校翔。這些新鑒定的SNP分別存放在NCBI dbSNP數(shù)據(jù)庫(kù)中，登錄號(hào)分別為1714766361它匕、1714766362展融、1714766363窖认、1714766364豫柬、1714766365、1714766367扑浸、1714766368和1714766370烧给。

產(chǎn)蛋性能多標(biāo)記回歸分析

在單標(biāo)記組合中，我們發(fā)現(xiàn)8個(gè)SNP對(duì)鵝卵數(shù)有顯著影響喝噪。采用線性回歸模型分析了多個(gè)顯著標(biāo)記對(duì)鵝產(chǎn)蛋性能的影響础嫡。8個(gè)SNP參與分析，以確定該性狀的基因組合或網(wǎng)絡(luò)(圖2)酝惧。建立了兩個(gè)網(wǎng)絡(luò)榴鼎，其中一個(gè)由兩個(gè)標(biāo)記組成，另一個(gè)由三個(gè)標(biāo)記組成晚唇。記錄-111407和記錄-106975被包括在兩個(gè)標(biāo)記網(wǎng)絡(luò)中(圖2A)巫财。預(yù)測(cè)值(在矩形中左邊)與相應(yīng)的實(shí)際值(在矩形中為右)有很高的相關(guān)性(r=0.98，r=0.81)哩陕。記錄-111407和記錄-106975分別對(duì)產(chǎn)蛋性能表現(xiàn)出加性和顯性效應(yīng)平项。用TT基因型替代GG/GA赫舒，可使平均蛋數(shù)減少9.45。記錄-111407闽瓢，A到T的轉(zhuǎn)換將導(dǎo)致卵數(shù)下降7.71接癌。三個(gè)標(biāo)記網(wǎng)絡(luò)引入了一個(gè)額外的標(biāo)記記錄-112359，它對(duì)產(chǎn)蛋率表現(xiàn)出加性效應(yīng)(圖2B)扣讼。T代G使卵數(shù)減少5.16缺猛。

12.png

13.png

攜帶產(chǎn)蛋相關(guān)SNP的新基因的鑒定

在上述8個(gè)候選基因的基礎(chǔ)上，對(duì)鵝的相關(guān)功能基因進(jìn)行了鑒定椭符。首先枯夜，利用候選RAD標(biāo)簽對(duì)NCBI公共數(shù)據(jù)庫(kù)進(jìn)行BLAST搜索，獲取相應(yīng)的序列艰山。然而湖雹，由于90-bp的RAD標(biāo)記太短，無(wú)法進(jìn)行有效的比對(duì)曙搬，因此沒(méi)有獲得明顯的匹配序列摔吏。因此，我們使用IPCR來(lái)擴(kuò)展候選RAD標(biāo)簽兩側(cè)的未知區(qū)域纵装。延伸的序列被用于進(jìn)一步的BLAST征讲。因?yàn)闆](méi)有可用的參考基因組信息對(duì)于鵝，我們主要使用檢索的鴨或雞序列用于鵝基因注釋橡娄。

14.png

對(duì)于記錄-134172，IPCR獲得了1,964bp的側(cè)翼序列長(zhǎng)度塑顺。DNA序列與鴨和雞的同源性分別為94%和79%汤求。根據(jù)記錄-134172的序列鑒定了KIAA1462基因。

記錄-112359的側(cè)翼序列長(zhǎng)度為2,164 bp严拒。DNA序列與鴨和雞的同源性分別為93%和79%扬绪。Rho GTPase激活蛋白21(ARHGAP 21)基因是根據(jù)記錄-112359序列鑒定的。

對(duì)于記錄-106582裤唠，通過(guò)IPCR獲得2,623 b長(zhǎng)度的側(cè)翼序列挤牛。DNA用鴨和雞序列分別顯示81％和65％的同源性，根據(jù)記錄-106582的序列鑒定了ACYCoA合成酶家族成員2(ACSF2)基因巧骚。

對(duì)于記錄-111407赊颠，通過(guò)IPCR獲得了1,508bp的側(cè)翼序列長(zhǎng)度格二。DNA 序列與鴨和雞序列的同源性分別為80%和72%。 根據(jù)記錄-111407序列鑒定Astrotactin2（AS TN2）基因竣蹦。

記錄-106975, 記錄-134172, 記錄-112359, 記錄-106582和記錄-111407的派生序列已存入GenBank(http://www.ncbi.nlm.nih.gov/genbank)))顶猜，登錄號(hào)分別為KP271033, KP271035, KP271036, KP271032和KP271034。

對(duì)記錄-88247痘括、記錄-135057和記錄-130775长窄，通過(guò)IPCR獲得3,100 bp、1,300bp和4,711 bp的側(cè)翼序列纲菌。3個(gè)DNA序列與鴨序列的同源性分別為90%挠日、93%和92%。記錄-88247和記錄-135057的DNA序列與雞的同源性分別為78%和82%翰舌。我們沒(méi)有從雞的全基因組shotgun contigs數(shù)據(jù)庫(kù)(分類號(hào)：9031)中找到記錄-130775相應(yīng)的序列嚣潜。

討論

基于池的RAD排序

在本研究中，我們采用了一種經(jīng)濟(jì)有效的比較RAD排序方法來(lái)發(fā)現(xiàn)與鵝產(chǎn)蛋性能有關(guān)的SNP椅贱。在沒(méi)有參考基因組的動(dòng)物中挖掘SNP有許多研究懂算。由于沒(méi)有參考基因組可供鵝使用，RAD測(cè)序提供了一種成本效益高的方法庇麦，可以在鵝基因組中發(fā)現(xiàn)大量的SNP计技。

轉(zhuǎn)換和顛換(ts/tv)的比值是DNA序列演化的一般性質(zhì)。對(duì)于已經(jīng)研究的所有基因組序列山橄，人們注意到垮媒，轉(zhuǎn)換發(fā)生的頻率在比顛換高，因?yàn)檗D(zhuǎn)換不需要改變構(gòu)象航棱。在本研究中睡雇，鵝的TS/TV比率為1.10，其符合轉(zhuǎn)換偏倚的規(guī)則丧诺。

下一代測(cè)序的淺測(cè)序深度是決定從序列數(shù)據(jù)生成基因型調(diào)用的質(zhì)量和測(cè)序成本的主要因素入桂。Catchen等人(2011)模擬RAD-seq過(guò)程奄薇，以測(cè)試USTACKS識(shí)別三葉松粘背植物基因座的能力驳阎。他們證明，平均測(cè)序深度為20×和40×馁蒂，是可靠的下一代測(cè)序測(cè)序深度呵晚，有較低的錯(cuò)誤率。在本研究中沫屡，LEBV組和HEBV組的平均測(cè)序深度分別為20×和17×饵隙，表明在深度上得到了可靠的測(cè)序結(jié)果。

與產(chǎn)蛋相關(guān)的SNP的發(fā)現(xiàn)

采用兩步策略沮脖，將基于池的RAD測(cè)序與個(gè)體驗(yàn)證相結(jié)合金矛，在較大的種群中發(fā)現(xiàn)與產(chǎn)蛋相關(guān)的SNP芯急。個(gè)體池的下一代測(cè)序(NGS)在基因組范圍內(nèi)的SNP發(fā)現(xiàn)中往往更有效，即使考慮到測(cè)序錯(cuò)誤驶俊，也能提供更準(zhǔn)確的等位基因頻率估計(jì)娶耍。DNA池NGS的成本效益更高的方法在各種研究中得到了廣泛的應(yīng)用，證明了DNA池的NGS允許在單個(gè)SNP上估計(jì)等位基因頻率饼酿，具有較高的準(zhǔn)確性榕酒，但文庫(kù)的構(gòu)建和測(cè)序工作相對(duì)較少。在我們的研究中故俐，通過(guò)比較LEBV和HEBV DNA池中估計(jì)的等位基因頻率想鹰，我們確定了467個(gè)與產(chǎn)蛋量相關(guān)的SNP。在適用于AS-PCR的467個(gè)SNP中药版，有55個(gè)在LEBV和HEBV隊(duì)列中進(jìn)行了個(gè)體基因分型辑舷。10個(gè)SNP在兩組間有不同的等位基因分布，陽(yáng)性率為18.2%槽片。與以前的研究相比惩妇，Turner等人(2010年)在蛇紋石和非蛇紋石土壤中檢測(cè)到擬南芥兩個(gè)DNA庫(kù)之間的840萬(wàn)個(gè)多態(tài)性。在840萬(wàn)個(gè)多態(tài)性中筐乳，96個(gè)在土壤類型之間的等位基因頻率差異大于80%歌殃。同時(shí)，在96個(gè)多態(tài)性的基礎(chǔ)上發(fā)現(xiàn)了81個(gè)基因蝙云。目前還沒(méi)有統(tǒng)一的標(biāo)準(zhǔn)來(lái)評(píng)估基于池測(cè)序的陽(yáng)性率多態(tài)數(shù)氓皱。然而，Gautier等人(2013) 基于池和個(gè)體的RAD測(cè)序數(shù)據(jù)評(píng)估等位基因頻率估計(jì)的準(zhǔn)確性勃刨。結(jié)果表明波材，DNA池測(cè)序是估算大量SNP位點(diǎn)等位基因頻率的一種經(jīng)濟(jì)有效的方法。朱等人(2012)實(shí)驗(yàn)證明身隐，DNA池測(cè)序是一種非常強(qiáng)大和成本效益高的技術(shù)廷区，可以在全基因組范圍內(nèi)發(fā)現(xiàn)SNP。

在本研究中贾铝，通過(guò)與小規(guī)模代表性個(gè)體的比較研究隙轻，進(jìn)一步驗(yàn)證了10個(gè)候選SNP在產(chǎn)蛋量較大的鵝群中的存在。8個(gè)SNP對(duì)蛋數(shù)有顯著影響垢揩，陽(yáng)性率為80%玖绿，小規(guī)模代表性比較與大規(guī)模驗(yàn)證具有較高的一致性。我們的結(jié)論是叁巨，基于池的RAD測(cè)序結(jié)合極具代表性的個(gè)體比較斑匪，是識(shí)別感興趣性狀的相關(guān)SNP的一種成本效益高的方法。

蛋產(chǎn)量的關(guān)聯(lián)分析

提高產(chǎn)蛋性能對(duì)于鵝的生產(chǎn)具有重要的意義锋勺。然而蚀瘸，由于繁殖力較低狡蝶，表型選擇對(duì)提高產(chǎn)蛋率的效率較低。對(duì)所涉及的遺傳標(biāo)記或基因進(jìn)行鑒定贮勃，有助于提高這一低遺傳力的性狀牢酵。許多研究人員一直致力于研究鵝繁殖性狀的遺傳機(jī)制。姜等（2011）檢測(cè)到PRL基因的5"側(cè)翼區(qū)SNP衙猪，以尋找影響皖西白鵝繁殖性狀的遺傳標(biāo)記馍乙。Chen等人（2012年）揭示了PRLR第10外顯子中SNP與皖江白鵝卵子性能的顯著關(guān)聯(lián)。Xu等人（2013）進(jìn)行了新轉(zhuǎn)錄體組裝和基因表達(dá)分析垫释，鑒定了大量與卵泡發(fā)育和生殖生物學(xué)相關(guān)的基因丝格，包括膽固醇側(cè)鏈裂解酶和多巴胺β羥基。Kang等人(2014)證實(shí)了烯醇ASE1(ENO1)基因的表達(dá)與孕前鵝相比棵譬，在產(chǎn)卵鵝的卵巢表達(dá)高显蝌，在四川白鵝的卵巢卵泡中鑒定出ENO1基因的表達(dá)譜。在我們的研究中订咸，我們清楚地證明了8個(gè)SNPs對(duì)鵝的產(chǎn)蛋性狀有顯著的影響曼尊。線性回歸進(jìn)一步揭示了產(chǎn)蛋數(shù)的兩個(gè)多重SNP網(wǎng)絡(luò)，涉及記錄-111407脏嚷、106975和112359骆撇。模型預(yù)測(cè)與觀測(cè)值有較好的一致性，驗(yàn)證了這些SNP對(duì)蛋數(shù)的組合效果父叙。以前的研究也報(bào)道了多個(gè)基因或標(biāo)記可以用來(lái)預(yù)測(cè)性狀神郊。Jiang等人.(2009)通過(guò)對(duì)多個(gè)標(biāo)記物的回歸分析，證實(shí)了兩個(gè)基因或三個(gè)基因網(wǎng)絡(luò)顯著地影響了肉牛的5或8個(gè)性狀趾唱。Ghazalpour等人(2006)利用微陣列和遺傳標(biāo)記數(shù)據(jù)構(gòu)建了小鼠肝臟基因共表達(dá)網(wǎng)絡(luò)涌乳，并研究了幾個(gè)基因模塊與小鼠體重的關(guān)系。因此甜癞，這8個(gè)SNP夕晓，特別是記錄-111407、106975和112359的組合可以是用于選擇鵝產(chǎn)蛋性能的有希望的分子標(biāo)記悠咱。我們使用David生物信息學(xué)資源6.7和UniHI在線工具進(jìn)一步探討了MAGI-1蒸辆，ARHGAP 21和ASTN 2的功能聯(lián)系，分別來(lái)自記錄-111407乔煞，106975和112359吁朦。分析結(jié)果表明，這些基因在任何信號(hào)通路或基因網(wǎng)絡(luò)中都沒(méi)有直接關(guān)聯(lián)渡贾。然而，MAG1-1和ARGAP21可以直接或間接地由分層(SFN)基因調(diào)節(jié)雄右。據(jù)報(bào)道空骚，在包括卵巢癌纺讲、子宮乳頭狀漿液癌、子宮平滑肌瘤囤屹、卵巢顆粒細(xì)胞瘤和類固醇細(xì)胞瘤的各種類型的人類疾病中熬甚，SFN的表達(dá)經(jīng)常丟失。Wang等人(2012)表明SFN的表達(dá)與雌激素和孕酮受體(ER和PR)呈負(fù)相關(guān)肋坚。Khongmanee等人(2013)顯示SFN在膽管癌細(xì)胞無(wú)凋亡耐藥中起重要作用乡括。大量的證據(jù)表明MAGI-1和ARHGAP 21在疾病和生殖功能中發(fā)揮作用的可能性很大。

基因克隆

利用IPCR擴(kuò)增RAD標(biāo)記智厌，結(jié)合鴨(TABID：8835)和雞(TABID：9031)公共數(shù)據(jù)庫(kù)的比對(duì)诲泌，共獲得5個(gè)新基因。與以往的研究相比铣鹏，我們沒(méi)有檢測(cè)到FSHβ敷扫，PRL，GnRH诚卸，LH和 PRLR等明確的生殖相關(guān)基因葵第。相反，我們發(fā)現(xiàn)三個(gè)新基因合溺，MAGI-1卒密，ARHGAP 21和KIAA 1462，可能在產(chǎn)蛋過(guò)程中起著重要的作用棠赛。實(shí)際上栅受，RAD-測(cè)序是一種通過(guò)酶切來(lái)減少基因組表達(dá)的方法。這種方法得到的SNP僅代表整個(gè)基因組的一小部分恭朗。本研究從LEBV和HEBV DNA池中分別獲得884,827和942,117個(gè)RAD標(biāo)記屏镊。估計(jì)平均覆蓋率占整個(gè)基因組的~6.96% (1.1GB，Anas platyrunchos)痰腮，因此而芥，上述已知的產(chǎn)蛋相關(guān)基因很有可能不能被列入所獲得的基因列表中。

對(duì)于MAG1-1基因膀值，Kranjecetal(2014)證明在HPV-陽(yáng)性細(xì)胞中可促進(jìn)細(xì)胞-細(xì)胞接觸棍丐，具有抑制細(xì)胞增殖和促進(jìn)的功能凋亡。ARHGAP21優(yōu)先用作GTPase活化蛋白(GAP)沧踏，CDC42調(diào)節(jié)ARP2/3復(fù)合物歌逢。它位于細(xì)胞核、細(xì)胞質(zhì)或核周區(qū)域翘狱，參與細(xì)胞-細(xì)胞粘附形成和細(xì)胞遷移秘案。KIAA1462是一種位于細(xì)胞核、胞漿和質(zhì)膜上的蛋白編碼基因。與KIAA1462相關(guān)的疾病包括動(dòng)脈疾病和冠狀動(dòng)脈疾病阱高。Akashi等人(2011)鑒定KIAA 1462為一種新的蛋白赚导，位于細(xì)胞與細(xì)胞的連接處，并得出結(jié)論：KIAA 1462在內(nèi)皮細(xì)胞與細(xì)胞的連接中的積累依賴于VE-cadherin介導(dǎo)的細(xì)胞-細(xì)胞粘附。卵母細(xì)胞的生長(zhǎng)是由膜細(xì)胞和顆粒細(xì)胞支持的，它們?cè)诼涯讣?xì)胞周圍建立了動(dòng)態(tài)和高度組織的細(xì)胞層捡硅。卵母細(xì)胞和顆粒細(xì)胞之間的縫隙連接是復(fù)雜的，并且在卵母細(xì)胞生長(zhǎng)的支持圈暗、減數(shù)分裂停滯的維持以及整個(gè)卵泡上皮的信號(hào)轉(zhuǎn)導(dǎo)中起主要作用。一個(gè)眾所周知的與細(xì)胞-細(xì)胞連接的建立相關(guān)的效應(yīng)是抑制細(xì)胞增殖裕膀。以上證據(jù)表明员串，這些基因(MAGI-1、ARHGAP 21和KIAA 1462)極有可能影響顆粒細(xì)胞的增殖和凋亡魂角，進(jìn)而干擾卵母細(xì)胞的生長(zhǎng)昵济。此外，KIAA1462在減數(shù)分裂重組中起著非常重要的作用野揪。Chowdhuryetal(2009)發(fā)現(xiàn)访忿，在人的重組率中KIAA1462是與變異相關(guān)的六個(gè)基因位點(diǎn)中的一個(gè)。減數(shù)分裂中的失敗或錯(cuò)誤可能導(dǎo)致不育斯稳、流產(chǎn)或出生缺陷海铆。

另外兩個(gè)克隆基因包括ACSF 2和ASTN2.ACSF 2是Acyl-CoA合成酶(ACS)家族成員，參與脂肪酸合成和三羧酸循環(huán)挣惰。ACSF 2是一種線粒體基質(zhì)酶卧斟，位于線粒體基質(zhì)中。ACSF 2基因在細(xì)胞代謝過(guò)程中的特征與線粒體功能有關(guān)憎茂，暗示ACSF 2基因可能在生殖過(guò)程中起一定作用珍语。ASTN 2在移行性小腦顆粒神經(jīng)元中高表達(dá)。它在神經(jīng)元的功能中起著重要的作用竖幔。Lesch等人(2008)鑒定ASTN 2基因參與細(xì)胞粘附和神經(jīng)細(xì)胞通訊板乙。Ahn等人(2010)發(fā)現(xiàn)一種來(lái)自ASTN 2基因內(nèi)含子的新的微RNA，并優(yōu)先在性腺中表達(dá)拳氢。

結(jié)論

我們將基于池的RAD測(cè)序策略應(yīng)用于鵝的SNP發(fā)現(xiàn)募逞。通過(guò)個(gè)體關(guān)聯(lián)分析驗(yàn)證了8個(gè)與產(chǎn)蛋相關(guān)的SNP。以產(chǎn)蛋相關(guān)的SNP為基礎(chǔ)馋评，克隆了5個(gè)鵝的新基因放接。我們的數(shù)據(jù)表明，這些SNP或基因可能是標(biāo)記輔助選擇動(dòng)物的候選標(biāo)記或靶標(biāo)留特。這些方法可以在其他生產(chǎn)動(dòng)物中進(jìn)行纠脾，可以提供給人類更高效玛瘸、更好的動(dòng)物食用/使用。我們的研究還表明乳乌，分子方法可以有目的的使用(而不是)簡(jiǎn)單地確定一些生理級(jí)聯(lián)的分子機(jī)制捧韵。事實(shí)上市咆，這方面的更多研究將有助于各種動(dòng)物的生產(chǎn)汉操。