利用微衛(wèi)星標(biāo)記指導(dǎo)紅鰭東方鲀親本選配

劉永新，劉奕，周勤，張紅濤

(1.中國水產(chǎn)科學(xué)研究院，北京 100141；2.中國水產(chǎn)科學(xué)研究院 珠江水產(chǎn)研究所，廣東 廣州 510380；3.中國水產(chǎn)科學(xué)研究院 北戴河中心實驗站，河北 秦皇島 066100)

?

利用微衛(wèi)星標(biāo)記指導(dǎo)紅鰭東方鲀親本選配

劉永新1，劉奕2，周勤3，張紅濤3

(1.中國水產(chǎn)科學(xué)研究院，北京 100141；2.中國水產(chǎn)科學(xué)研究院 珠江水產(chǎn)研究所，廣東 廣州 510380；3.中國水產(chǎn)科學(xué)研究院 北戴河中心實驗站，河北 秦皇島 066100)

摘要:為制定紅鰭東方鲀家系配組方案提供可行性指導(dǎo)，利用微衛(wèi)星標(biāo)記輔助紅鰭東方鲀Takifugu rubripes家系的建立,選擇34個微衛(wèi)星標(biāo)記對紅鰭東方鲀兩個群體進行遺傳評估。結(jié)果表明：A群體的平均等位基因數(shù)(Na)和Nei基因多樣性指數(shù)(He)分別為6.647 0和0.711 5，B群體的相應(yīng)值分別為4.647 1和0.646 1，且兩個群體之間的遺傳差異達到顯著性水平(P<0.05)；群體間的遺傳分化系數(shù)(GST)為0.050 7，基因流(Nm)為4.679 6, 表明紅鰭東方鲀?nèi)后w間存在低程度的遺傳分化和一定程度的基因交流；分子方差分析(AMOVA)表明，遺傳變異主要存在于群體內(nèi)，所占比例為78.49%，而群體之間僅占21.51%；根據(jù)群體間Nei氏遺傳距離構(gòu)建UPGMA系統(tǒng)樹，群體內(nèi)個體之間的遺傳距離為0.11~0.82，依據(jù)遺傳距離的遠近將全部個體分成多個分支，不同分支內(nèi)含有數(shù)個個體；各項遺傳參數(shù)表明，試驗群體具備一定程度的遺傳變異，不同個體之間擁有相對較遠的親緣關(guān)系，可以根據(jù)個體之間遺傳距離制定育種計劃，建立家系進行選育研究。研究表明，利用微衛(wèi)星標(biāo)記信息構(gòu)建紅鰭東方鲀家系的方法是切實可行的。

關(guān)鍵詞:紅鰭東方鲀；群體；微衛(wèi)星；遺傳距離

紅鰭東方鲀Takifugurubripes隸屬于硬骨魚綱Osteichthyes、鲀形目Tetraodontiformes、鲀亞目Tetraodontoidei、鲀總科Tetraodntoidea、鲀科Tetraodontidae、東方鲀屬Takifugu，俗稱河豚，主要分布于中國的渤海、黃海和東海，以及日本沿海和朝鮮半島[1]。紅鰭東方鲀?nèi)赓|(zhì)細嫩、味道鮮美、高蛋白、低脂肪并含有豐富的維生素和微量元素，素有“魚中之王”的美譽[2-5]。紅鰭東方鲀的養(yǎng)殖始于20世紀(jì)80年代日本，其養(yǎng)殖方式以陸上水池人工養(yǎng)殖為主，但養(yǎng)殖產(chǎn)量無法滿足日本民眾的需求，因此，每年還需從中國進口大量的商品魚。中國的紅鰭東方鲀?nèi)斯ゐB(yǎng)殖始于1991年， 但由于中國自然資源量明顯小于日本，而且紅鰭東方鲀和假睛東方鲀Takifugupseudommus兩個物種具有相近的繁殖季節(jié)和棲息海區(qū)，非常容易發(fā)生混交，從而為中國的紅鰭東方鲀養(yǎng)殖生產(chǎn)和市場穩(wěn)定帶來了不利影響[6]。國內(nèi)養(yǎng)殖的品種價格不高、種質(zhì)不純，出口受限制是養(yǎng)殖戶面臨的實際問題。中國水產(chǎn)科學(xué)研究院黃海水產(chǎn)研究所在1996年首次引進日本紅鰭東方鲀種群，在國內(nèi)確立了純種引進技術(shù)，同時形成了全人工純種育苗的大規(guī)模生產(chǎn)。迄今為止，國內(nèi)有很多養(yǎng)殖公司仍然采用從日本引種的方式進行生產(chǎn)，在繁育了1~2代后養(yǎng)殖苗種明顯出現(xiàn)生長速度緩慢、抗逆性降低的現(xiàn)象，為了解決這些問題只能再次從國外引種，因此，紅鰭東方鲀養(yǎng)殖業(yè)受國外種質(zhì)資源控制的影響非常嚴(yán)重。如何利用引進群體和國內(nèi)本土資源開展良種選育研究，盡快培育出滿足市場需求的優(yōu)良品種并提供長期穩(wěn)定的供應(yīng)，是當(dāng)前育種工作者亟待解決的根本問題。

國內(nèi)對紅鰭東方鲀遺傳的研究雖有所開展，但多以形態(tài)性狀[7-9]、群體結(jié)構(gòu)和生化遺傳變異[10-11]、分子標(biāo)記開發(fā)[12-14]和種內(nèi)雜交[15-16]等研究為主，未見有系統(tǒng)性的優(yōu)良品種選育進程報道。到目前為止，沒有紅鰭東方鲀品種通過國家原良種審定委員會審定并予以推廣養(yǎng)殖。選擇育種技術(shù)是水產(chǎn)動物品種遺傳改良過程中應(yīng)用最為廣泛的方法，雜交育種和分子標(biāo)記育種等育種技術(shù)都要以選擇育種為基礎(chǔ)來開展。即使是現(xiàn)在，部分水產(chǎn)動物品種都已完成了基因組測序，隨之而產(chǎn)生的全基因組育種技術(shù)也同樣離不開選擇育種。從已有研究基礎(chǔ)來看，紅鰭東方鲀主要經(jīng)濟性狀的遺傳改良采用常規(guī)的選擇育種技術(shù)更為適合。開展選擇育種首先需要構(gòu)建具有豐富遺傳變異的基礎(chǔ)群體；其次，要借助分子標(biāo)記技術(shù)分析群體的遺傳結(jié)構(gòu)，判斷遺傳變異度和群體內(nèi)個體之間的遺傳距離，制定相應(yīng)的親本配組方案，建立選育家系[17]。多種分子標(biāo)記都能夠進行群體遺傳結(jié)構(gòu)分析，隨機擴增片段長度多態(tài)性DNA、擴增片段長度多態(tài)性和微衛(wèi)星DNA技術(shù)先后在不同歷史時期應(yīng)用于水產(chǎn)動物遺傳研究中?，F(xiàn)階段，微衛(wèi)星DNA分子標(biāo)記是進行群體分析的首選標(biāo)記，眾多水產(chǎn)動物品種如牙鲆Paralichthysolivaceus[18]、大菱鲆Scophthalmusmaximus[19]、黃顙魚Pseudobagrusfulvidraco[20]、大黃魚Pseudosciaenacrocea[21]、吉富羅非魚Oreochromisniloticus[22]、赤點石斑魚Epinephelusakaara[23]等都利用這種標(biāo)記分析了目標(biāo)群體的遺傳基礎(chǔ)。本研究中，以從日本進口的紅鰭東方鲀受精卵發(fā)育至性成熟的親本組成育種群體，利用34個微衛(wèi)星DNA分子標(biāo)記分析該群體的遺傳結(jié)構(gòu)，了解其遺傳變異水平并進行UPGMA系統(tǒng)聚類，確定個體之間的遺傳距離，旨在為大規(guī)模建立選育家系和開展分子標(biāo)記輔助育種研究工作提供基礎(chǔ)數(shù)據(jù)和理論指導(dǎo)。

1材料與方法

1.1材料

試驗用紅鰭東方鲀親魚共計225尾，其中199尾來自河北省唐山市天河水產(chǎn)養(yǎng)殖公司(記為A群體)，26尾來自河北省昌黎縣河鲀魚良種場(記為B群體)。天河養(yǎng)殖有限公司從1995—1999年，每年由日本直接引進受精卵進行苗種孵化，并從事養(yǎng)殖，2000年時上述引進的苗種達到性成熟，將這些成熟的河鲀作為親魚；之后為防止近交衰退，每兩到三年由日本引進一次受精卵更新種魚，試驗所采集的A群體為該公司養(yǎng)殖至今保有的親魚群體。河北省昌黎縣河鲀魚良種場親魚B群體也是由日本進口的受精卵經(jīng)孵化、培育至性成熟的個體產(chǎn)生的后代，但其引種時間不同于天河養(yǎng)殖有限公司，為2005—2007年。所有親魚年齡為4～6齡，體質(zhì)量為0.8～3.5 kg。采集親魚尾部鰭條用于提取基因組DNA。

1.2方法

1.2.1微衛(wèi)星標(biāo)記的篩選與PCR擴增選擇由日本學(xué)者[24]和中國學(xué)者[10-12]開發(fā)的34個高多態(tài)性微衛(wèi)星標(biāo)記進行試驗分析。標(biāo)記名稱、引物序列和退火溫度列于表1。由上海生工生物工程服務(wù)有限公司負責(zé)引物合成。利用TIANGEN海洋動物基因組DNA提取試劑盒進行基因組DNA的提取。提取完成后，采用瓊脂糖凝膠電泳檢測DNA樣品的質(zhì)量，用核酸蛋白測定儀檢測 DNA 濃度(A260/280值在1.8～2.0 之間較好)；之后，稀釋 DNA 至 50 ng/μL，于冰箱(-20 ℃)中保存?zhèn)溆谩CR反應(yīng)體系(共25 μL)：10×Buffer 2.5 μL,Mg2+(25 mmol/L)1 μL,dNTPs(各2 mmol/L)1 μL,上、下游引物(10 μmol/L)各1 μL,模板DNA (50 ng/μL)1 μL,TaqDNA聚合酶1 U，用ddH2O補足至25 μL。PCR反應(yīng)程序：94 ℃下預(yù)變性3 min；94 ℃下變性30 s，退火30 s，72 ℃下延伸30 s，共進行30個循環(huán)；最后在72 ℃下延伸10 min。在AB9700型PCR儀上進行PCR擴增。

1.2.2PCR擴增產(chǎn)物的檢測PCR擴增產(chǎn)物經(jīng)8%非變性聚丙烯酰胺凝膠電泳檢測，上樣緩沖液由0.25%的溴酚藍和40%的蔗糖水溶液組成。電泳后用1%硝酸銀染色10 min，顯色液顯色10 min。凝膠在HP scanjet G4010掃描儀上成像，用Gel-Pro Analyzer 4.5凝膠分析軟件對電泳譜帶進行數(shù)據(jù)采集和分析。

1.3數(shù)據(jù)處理

每對微衛(wèi)星引物檢測一個位點，每條多態(tài)性帶為一個等位基因。利用Popgene 3.2軟件[25]統(tǒng)計每個微衛(wèi)星標(biāo)記的等位基因數(shù)(Number of alleles，Na)、Nei基因多樣性指數(shù)(Nei’s gene diversity index，He)和Nei遺傳距離(Genetic Distance，GD)[26]。采用Wilcoxon 配對檢驗法進行A群體和B群體等位基因數(shù)和Nei 基因多樣性指數(shù)的顯著性檢驗。根據(jù)Nei遺傳距離進行非加權(quán)類平均法(UPGMA)系統(tǒng)聚類?；诿恳粋€位點的等位基因差異，采用Arlequin 3.0軟件[27]中的分子方差分析程序(AMOVA)計算F統(tǒng)計量(F-statistics，F(xiàn)st)，分析紅鰭東方鲀?nèi)后w的遺傳分化，采用Excoffier 等[28]的非參數(shù)置換方法(3000次)進行Fst的顯著性檢驗。

2結(jié)果與分析

2.1紅鰭東方鲀?nèi)后w的遺傳參數(shù)

利用34個微衛(wèi)星標(biāo)記進行紅鰭東方鲀?nèi)后w的遺傳分析，統(tǒng)計每個位點的等位基因數(shù)和Nei基因

表1　34個微衛(wèi)星標(biāo)記的名稱、引物序列和退火溫度

多樣性指數(shù)。從表2可見，不同的微衛(wèi)星位點具有較大的差異性。A群體的等位基因數(shù)為3~13，B群體的為0~12；A群體的Nei基因多樣性指數(shù)為0.456 6~0.875 0，B群體的為0~0.840 2。從表3 可知，兩個群體的等位基因數(shù)和Nei 基因多樣性指數(shù)的平均值均較高, 表明這兩個群體資源均具有較豐富的多樣性。Wilcoxon 配對測驗表明，A群體平均等位基因數(shù)(6.647 0)和Nei基因多樣性指數(shù)(0.711 5)均高于B群體(4.647 1和0.646 1)，且具有顯著性差異(P<0.05)。由此表明，A群體具有更加豐富的遺傳多樣性。

2.2紅鰭東方鲀?nèi)后w的遺傳分化

根據(jù) Nei基因多樣性指數(shù)計算兩個群體總的遺傳分化系數(shù)(GST)為0.050 7，由GST計算種群間的基因流值(Nm)為4.679 6，說明這兩個紅鰭東方鲀?nèi)后w之間存在著一定程度的基因交流。分子方差分析(AMOVA)結(jié)果表明，遺傳變異主要存在于群體內(nèi)部，所占比例為78.49%，而群體之間的遺傳變異僅占21.51%(表4)。

表2　紅鰭東方鲀兩個群體在34個微衛(wèi)星位點的等位基因數(shù)和Nei基因多樣性指數(shù)

表3紅鰭東方鲀兩個群體在34個微衛(wèi)星位點的平均等位基因數(shù)和Nei基因多樣性指數(shù)

Tab.3Mean number of alleles and average Nei’s gene diversity index at 34 microsatellite loci in two populations of redfin puffer

群體group樣本數(shù)numberofsample等位基因數(shù)Na平均值±標(biāo)準(zhǔn)差mean±S.D.P值PvalueNei基因多樣性指數(shù)He平均值±標(biāo)準(zhǔn)差mean±S.D.P值PvalueAB199266.6470±2.6501a4.6471±1.9677b0.00050.7115±0.0914a0.6461±0.1445b0.0240

注：同列中標(biāo)有不同小寫字母者表示組間有顯著性差異(P<0.05)

Note:The means with different letters within the same column are significant differences at the 0.05 probability level

2.3聚類分析

表4　紅鰭東方鲀?nèi)后w的分子方差分析

兩個群體共計225個個體的Nei遺傳距離為0.11~0.82，表明這些個體之間既有遺傳差異較大的個體，也存在遺傳差異較小的個體?；贜ei遺傳距離進行UPGMA系統(tǒng)聚類(圖1)，在Nei遺傳距離為0.82處可以明顯地劃分為兩類：A分支和B分支，A分支包含221個個體，而B分支僅包含來自B群體的4個個體(B23~B26)。在遺傳距離為0.66處同樣明顯地分成兩類，較大的分支為C分支，包括192個個體，較小的分支為D分支，包括29個個體。C分支在遺傳距離為0.65處分為E、F兩個分支，分別包含95和97個個體，其中E分支在遺傳距離為0.59處又分為I和J兩個分支。而F分支在遺傳距離為0.61處分為G和H兩個分支，其中G分支在遺傳距離為0.57處分為K和L兩個分支，而H分支在遺傳距離為0.59處分為M和N兩個分支。C分支在遺傳距離小于0.60處細化成的6個分支I、J、K、L、M和N所包含的個體數(shù)依次為48、47、26、21、48和2個，其中I、J和K分支基本全部由來自A群體的個體組成，而大多數(shù)來自B群體的個體集合在M分支內(nèi)。在遺傳距離為0.11～0.60時，I、J、K、L、M和N分支內(nèi)的個體按照遺傳距離的遠近程度聚集在不同的小分支內(nèi)。

圖1　親魚群體的親緣關(guān)系聚類圖Fig.1　Cluster dengrodogram of parental population

3討論

3.1群體遺傳學(xué)參數(shù)

傳統(tǒng)的水產(chǎn)動物選擇育種在進行家系構(gòu)建時往往僅根據(jù)個體表型性狀是否優(yōu)良確定其是否作為繁殖用親本。伴隨著分子生物學(xué)技術(shù)日新月異的發(fā)展，不同種類的分子標(biāo)記已經(jīng)引入到家系選育過程中，通過判斷基礎(chǔ)群體的遺傳變異和群體內(nèi)個體之間的遺傳距離，可為親本個體的選擇和配組方案的制定提供參考依據(jù)。利用這種方法組建家系的魚類品種有建鯉Cyprinuscarpio[29]、團頭魴Megalobramaamblycephala[30]、虹鱒Oncorhynchusmykiss[31]等。遺傳多樣性是種群進化和遺傳改良的前提和基礎(chǔ)[32]，物種的遺傳多樣性越高，其獲得的遺傳進展就越大。本研究中，所選群體的平均等位基因數(shù)和Nei基因多樣性指數(shù)分別為4.647 1和0.646 1以上，與已有報道的紅鰭東方鲀?nèi)毡救后w和秦皇島群體的參數(shù)基本相同[10]，而顯著高于公司養(yǎng)殖群體和實驗室養(yǎng)殖群體[11]。需要指出的是，雖然本試驗中選擇了幾百個個體進行分析，但這些個體全部來自引進的受精卵，不敢確定一定存在遺傳漂變，同時也無法確認(rèn)是否能夠代表日本原來的自然群體。為此，本研究中還選擇了萬玉美等[10]進行紅鰭東方鲀兩個自然群體遺傳學(xué)研究時所用的微衛(wèi)星標(biāo)記，因為他們選擇和評估這些標(biāo)記，依據(jù)的是一個能代表自然群體的理想群體。除了上述兩個遺傳指標(biāo)外，筆者也計算了Hardy-Weinberg遺傳偏離指數(shù)，發(fā)現(xiàn)在當(dāng)前群體中存在平衡偏離問題，而依據(jù)標(biāo)記計算遺傳多樣性，描述群體結(jié)構(gòu)的意義并不大，因為如果放歸這些個體，這一結(jié)構(gòu)馬上會解體而不存在。對指導(dǎo)配組工作而言，試驗的主要目的不是描述多樣性和群體結(jié)構(gòu)，而是獲取個體之間遺傳差異的信息，為實際育種服務(wù)。以本研究中所選的遺傳參數(shù)指標(biāo)等位基因數(shù)為例，每個位點的等位基因數(shù)都不相同，這其中就存在一個取樣飽和度的問題。一個位點有十幾個等位基因，盡管有幾百個個體，仍會有很多等位基因只有一個，頻率極低，這時群體內(nèi)個體數(shù)量就不是30個，而需要300個甚至更多，因此，不是等位基因數(shù)愈多越好，只是在尋找個體差異的時候，等位基因數(shù)是越多越好。這里所強調(diào)的不是等位基因數(shù)量對群體研究的要求，而是等位基因數(shù)量對尋找個體差異的重要意義。

3.2群體內(nèi)個體選配

分子方差分析表明，目標(biāo)群體的大部分遺傳變異存在于群體內(nèi)部，少量的遺傳差異存在于群體之間。由于AMOVA分析不必作Hardy-Weinberg假設(shè)，避免了數(shù)學(xué)假設(shè)帶來的系統(tǒng)誤差，其結(jié)合UPGMA系統(tǒng)聚類結(jié)果，可以更好地輔助微衛(wèi)星標(biāo)記進行個體遺傳差異評估，作為建立家系的依據(jù)。魚類的親緣關(guān)系遠近程度常用遺傳距離來衡量，根據(jù)遺傳距離對親魚進行聚類分析，能夠為親本配組和雜種優(yōu)勢預(yù)測提供參考依據(jù)[32]。本研究中的UPGMA系統(tǒng)聚類結(jié)果中，A分支和B分支在遺傳距離為0.82處分開，這表明兩個分支內(nèi)的個體具有較遠的親緣關(guān)系，原則上可以作為配組對象。但是，B分支內(nèi)個體數(shù)量過少，僅包括4個個體。這在實際育種工作中就會面臨兩個問題，其一，在繁殖季節(jié)親魚一般都需要進行營養(yǎng)強化，從而保證其能夠孕育健康、優(yōu)質(zhì)的精子和卵子。親魚強化過程需要較為精湛的培育技術(shù)，在這一過程中雄性親魚容易產(chǎn)生活力不錯的精子，但對雌性親魚往往不能保證全部待產(chǎn)個體都能夠產(chǎn)生好的卵子，卵子質(zhì)量較差或者生殖腺膨大卻不能產(chǎn)出卵子的情況時有發(fā)生，從而導(dǎo)致預(yù)先設(shè)計的交配方案不能實施。因此，對于僅包含幾個個體的B分支，一旦個體性腺發(fā)育不理想則無法進行配組。其二，即使B分支內(nèi)個體性腺正常發(fā)育，但其數(shù)量過少，很容易成為一個單性別小群體，在配組時，只能做正交或反交的單向設(shè)計，且組合數(shù)也較少，產(chǎn)生性狀優(yōu)良組合的幾率也將降低。此外，AMOVA分析結(jié)果表明，主要的遺傳變異來自于群體內(nèi)部而非群體之間，因此，應(yīng)盡量考慮設(shè)計A群體或B群體內(nèi)的個體配組。鑒于上述幾點考慮，分支D、I、J、K、L、M和N之間的遺傳距離(0.57～0.65)相對較遠，這些分支內(nèi)的個體可設(shè)計交配方案。這7個分支內(nèi)還存在細化的小分支，這些小分支的遺傳距離范圍為0.11～0.57。有文獻報道，家系生長在一定范圍內(nèi)隨著親本遺傳距離的增大而逐漸增加，但是隨著遺傳距離的進一步增大，生長有下降的趨勢。這表明配組親本的遺傳距離與生長之間存在一個閾值點，一旦超過這個閾值點，即使遺傳距離再增加其生長優(yōu)勢反而會降低，這一論斷在鯉中已經(jīng)取得試驗驗證[33]，除水產(chǎn)物種外，對部分作物的研究也得出了相似的結(jié)論[34-37]。參考這些研究結(jié)果，在本試驗群體中對親本遺傳距離為0.11～0.60范圍內(nèi)的分支也可設(shè)計交配，一方面可以獲得更多的組合，有助于在大規(guī)模組合中優(yōu)中選優(yōu)，另一方面可以驗證紅鰭東方鲀這一物種是否也存在類似的“閾值”理論。一旦發(fā)現(xiàn)閾值點，在今后的育種工作中便可集中于這一區(qū)域進行配組設(shè)計，從而加速優(yōu)良組合的產(chǎn)生，用于新品種或新品系的篩選。本試驗中對紅鰭東方鲀?nèi)后w遺傳的分析結(jié)果，可以直接用于指導(dǎo)親魚配組實踐，開展大規(guī)模家系選育，進行生長、抗病、品質(zhì)、性別等多元化性狀的遺傳改良研究。育種材料是進行選育的根本所在，本研究中僅收集了日本群體的后代，如果想要獲得持續(xù)、穩(wěn)步的遺傳進展還應(yīng)收集中國不同海域的自然資源或繼續(xù)引進韓國、朝鮮等外來資源，不斷增加基礎(chǔ)群體的遺傳變異，才能為盡快培育出滿足養(yǎng)殖需求的新品種奠定堅實基礎(chǔ)。

參考文獻：

[1]孟慶聞,蘇錦祥,繆學(xué)祖.魚類分類學(xué)[M].北京:中國農(nóng)業(yè)出版社,1995：970.

[2]盧敏德,葛志亮.對我國4種東方鲀營養(yǎng)價值的評價[J].中國水產(chǎn)科學(xué),1999,6(4):90-94.

[3]Aparicio S,Chapman J,Stupka E,et al.Whole-genome shotgun assembly and analysis of the genome ofFugurubripes[J].Science,2002,297:1301-1310.

[4]劉秀云,康大海,孫效文,等.紅鰭東方鲀F1代混養(yǎng)群體的基因型鑒定及QTL分析[J].水產(chǎn)科學(xué),2008,27(12):641-644.

[5]孫緒文,王群,顏顯輝,等.土霉素在紅鰭東方鲀體內(nèi)的殘留及休藥期研究[J].漁業(yè)科學(xué)進展,2009,30(6):75-80.

[6]陸麗君.紅鰭東方鲀國內(nèi)主要養(yǎng)殖區(qū)養(yǎng)殖群體遺傳多樣性研究及與性別相關(guān)微衛(wèi)星標(biāo)記篩選[D].上海：上海海洋大學(xué),2012.

[7]王新安,馬愛軍,莊志猛,等.紅鰭東方鲀(Takifugurubripes)形態(tài)性狀對體重的影響效果[J].海洋與湖沼,2013,44(1):135-140.

[8]王新安,馬愛軍,莊志猛,等.紅鰭東方鲀體型性狀選育指標(biāo)的綜合判定[J].漁業(yè)科學(xué)進展,2012,33(6):10-15.

[9]姜志強,何旭穎,王冰,等.不同家系紅鰭東方鲀形態(tài)差異研究[J].現(xiàn)代農(nóng)業(yè)科技,2012,17(1):255-257.

[10]萬玉美,王蕾,譚照君,等.紅鰭東方鲀兩個群體的遺傳結(jié)構(gòu)及與經(jīng)濟性狀的相關(guān)性分析[J].淡水漁業(yè),2011,41(5):9-16.

[11]陳艷翠,高天翔,陳四清,等.紅鰭東方鲀養(yǎng)殖群體生化遺傳變異初步研究[J].海洋水產(chǎn)研究,2007,28(4):14-20.

[12]郝君,孫效文,孟雪松.紅鰭東方鲀微衛(wèi)星DNA多態(tài)性初步分析[J].上海水產(chǎn)大學(xué)學(xué)報,2006,15(1):21-24.

[13]崔建洲,申雪艷,楊官品,等.紅鰭東方鲀基因組微衛(wèi)星特征分析[J].中國海洋大學(xué)學(xué)報,2006,36(2):249-254.

[14]郝君,孫效文,孟雪松.紅鰭東方鲀BAC數(shù)據(jù)庫和ESTs數(shù)據(jù)庫中微衛(wèi)星的篩選與應(yīng)用[J].大連水產(chǎn)學(xué)院學(xué)報,2007,22(2):97-101.

[15]萬玉美，趙海濤，張福崇.紅鰭東方鲀、菊黃東方鲀及其雜交F1幼魚耗氧率與窒息點研究[J].水產(chǎn)科學(xué),2013,32(1):21-25.

[16]范文濤,劉海金,趙文江,等.菊黃東方鲀♀× 紅鰭東方鲀♂雜交后代早期形態(tài)特征及生長速度的比較[J].水產(chǎn)學(xué)報,2011,35(7):1065-1071.

[17]谷晶晶,朱迎軍,孟雪松,等.微衛(wèi)星標(biāo)記對紅鰭東方鲀繁殖的指導(dǎo)應(yīng)用及遺傳分析[J].水產(chǎn)科學(xué),2010,29(9):527-531.

[18]馬曉冰,王桂興,劉海金,等.秦皇島海域野生牙鲆群體遺傳多樣性分析[J].中國水產(chǎn)科學(xué),2012,19(6):963-969.

[19]侯仕營,馬愛軍,王新安,等.大菱鲆4個引進地理群體遺傳多樣性的微衛(wèi)星分析[J].漁業(yè)科學(xué)進展,2011,32(1):16-23.

[20]王婷婷,宋學(xué)宏,許愛國,等.應(yīng)用微衛(wèi)星標(biāo)記分析4個黃顙魚群體的遺傳多樣性[J].江蘇農(nóng)業(yè)科學(xué),2012,40(4):41-45.

[21]陳淑吟,徐士霞,張志勇,等.大黃魚野生群體與養(yǎng)殖群體遺傳多樣性研究[J].海洋科學(xué),2011,35(12):82-87.

[22]劉海情,郭昱嵩,王中鐸,等.利用微衛(wèi)星分析吉富羅非魚群體的遺傳多樣性[J].南京農(nóng)業(yè)學(xué)報,2012,43(1):94-98.

[23]陳省平,丁少雄,陳嘉慧,等.赤點石斑魚群體遺傳結(jié)構(gòu)的微衛(wèi)星分析[J].中山大學(xué)學(xué)報,2012,51(3):83-89.

[24]古川聰史.紅鰭東方鲀與高生長性狀有關(guān)的遺傳學(xué)及分子生物學(xué)的研究[D].東京:東京大學(xué),2009.

[25]Raymond M,Rousset F.GENEPOP (version 1.2):population genetics software for exact tests and ecumenicism[J].Journal of Heredity,1995,86(3):248.

[26]Nei M.Analysis of gene diversity in subdivided populations[J].Proceedings of the National Academy of Sciences of the United States of America,1973,70(12):3321-3323.

[27]Excoffier L,Laval L G,Schneider S.Arlequin ver 3.0:an integrated software package for population genetics data analysis[J/OL].Evolutionary Bioinformatics Online,2005(1):47-50.

[28]Excoffier L,Smouse P,Quattro J.Analysis of molecular variance inferred from metric distances among DNA haplotypes:application to human mitochondrial DNA restriction data[J].Genetics,1992,131(2):479-491.

[29]李建林,唐永凱,李紅霞,等.6個建鯉家系的遺傳結(jié)構(gòu)及不同親緣關(guān)系個體間的遺傳差異分析[J].大連海洋大學(xué)學(xué)報,2013,28(2):166-170.

[30]曾聰,張耀,曹小娟,等.團頭魴3個地理種群雜交效果的配合力和微衛(wèi)星標(biāo)記預(yù)測[J].水產(chǎn)學(xué)報,2012,36(6):809-814.

[31]王炳謙,谷偉,高會江,等.利用配合力和微衛(wèi)星標(biāo)記預(yù)測虹鱒品系間的雜交優(yōu)勢[J].中國水產(chǎn)科學(xué),2009,16(2):206-213.

[32]陳靈芝.中國的生物多樣性現(xiàn)狀及其保護對策[M].北京:科學(xué)出版社,1993:93-113.

[33]魯翠云,曹頂臣,孫效文,等.微衛(wèi)星分子標(biāo)記輔助鏡鯉家系構(gòu)建[J].中國水產(chǎn)科學(xué),2008,15(6):893-901.

[34]李明爽,洪脫,龔永生,等.雜種優(yōu)勢預(yù)測研究進展[J].中國農(nóng)學(xué)通報,2008,24(1):117-122.

[35]陳平華，張輝，陳如凱.利用RAPD 技術(shù)進行甘蔗育種親本輔助選擇[J].分子植物育種,2004,2(5):675-681.

[36]高世斌，榮廷昭，李晚忱，等.19 個玉米自交系的數(shù)量性狀和AFLP遺傳差異比較研究[J].華北農(nóng)學(xué)報,2004,19(2):24-27.

[37]李伯群，余國東，石有明，等.兩系雜交小麥雜種優(yōu)勢、配合力及遺傳距離分析[J].麥類作物學(xué)報,2003,23(3):17-21.

Application of microsatellite markers in parental apolegamy

of redfin pufferTakifugurubripes

LIU Yong-xin1, LIU Yi2, ZHOU Qin3, ZHANG Hong-tao3

(1.Chinese Academy of Fishery Sciences, Beijing 100141, China; 2.Pearl River Fisheries Research Institute, Chinese Academy of Fishery Sciences, Guangzhou 510380, China; 3.Beidaihe Central Experiment Station, Chinese Academy of Fishery Sciences, Qinhuangdao 066100, China)

Abstract：A family mating plan of redfin puffer Takifugu rubripes was established by parental apolegamy assisted by microsatellite markers and by analysis of genetic variations and distances of parental populations, and genetic evaluation was carried out in two populations of the redfin puffer by a set of 34 microsatellite markers. The results showed that there were significantly higher allele number (6.647 0) and Nei’s gene diversity index (0.711 5) in population A than those (allele number 4.647 1 and Nei’s gene diversity index 0.646 1) in population B, significant genetic differences between the two populations (P<0.05). The genetic differentiation coefficient GSTof 0.050 7 and Nmof 4.679 6 indicated that there were inferior differentiation and some degree of gene exchange among the redfin puffer populations. AMOVA revealed that there was 78.49% variation within the population while there was 21.51% variation between the two populations. The UPGMA tree established by genetic distance indicated that the genetic distance (GD) among individuals within the population was varied from 0.11 to 0.82. All individuals were divided into different branches on basis of genetic distance magnitude. Each genetic parameter demonstrated that there were genetic variations to some extent in experimental population and distant relatives among individuals. The breeding schedule can be made according to the genetic distances among individuals, and then families can be produced for selective breeding. The findings indicate that the method is feasible for establishment of families in redfin puffer with microsatellite marker information.

Key words：Takifugu rubripes; population; microsatellite; genetic distance

通信作者：周勤(1963—)，高級工程師。E-mail:laozhou529@hotmail.com

作者簡介：劉永新(1979—)，博士，副研究員。E-mail:liuyx@cafs.ac.cn

基金項目：國家“十二五”科技支撐計劃項目(2012BAD26B01)；中國水產(chǎn)科學(xué)研究院基本科研業(yè)務(wù)費項目(2013A0505)

收稿日期：2014-06-27

中圖分類號：S917

文獻標(biāo)志碼：A

文章編號：2095-1388(2015)02-0113-07

DOI：10.3969/J.ISSN.2095-1388.2015.02.001