育成品種易捕鯉F5代保種群體的遺傳結(jié)構(gòu)研究

趙新春，賈智英，李盛文，李池陶，石連玉，2

(1 中國(guó)水產(chǎn)科學(xué)研究院 黑龍江水產(chǎn)研究所，黑龍江 哈爾濱 150070 ；2 上海海洋大學(xué) 水產(chǎn)與生命學(xué)院，上海 201306)

?

育成品種易捕鯉F5代保種群體的遺傳結(jié)構(gòu)研究

趙新春1,2，賈智英1，李盛文1,2，李池陶1，石連玉1，2

(1 中國(guó)水產(chǎn)科學(xué)研究院 黑龍江水產(chǎn)研究所，黑龍江 哈爾濱 150070 ；2 上海海洋大學(xué) 水產(chǎn)與生命學(xué)院，上海 201306)

[摘要]【目的】 研究易捕鯉的群體遺傳結(jié)構(gòu)，為其優(yōu)良性狀的穩(wěn)定保持和種質(zhì)資源的可持續(xù)利用奠定基礎(chǔ)?！痉椒ā?提取96尾易捕鯉個(gè)體的鰭條基因組DNA，以其為模板，利用30對(duì)微衛(wèi)星標(biāo)記進(jìn)行PCR擴(kuò)增，根據(jù)基因型，利用生物信息學(xué)軟件對(duì)各種遺傳參數(shù)進(jìn)行分析?！窘Y(jié)果】 在96尾易捕鯉個(gè)體中，共檢測(cè)到154個(gè)等位基因，每個(gè)位點(diǎn)的等位基因數(shù)3～9個(gè)，平均有效等位基因數(shù)為4.060 7；期望雜合度為0.596 9～0.871 0，平均為0.740 0；多態(tài)信息含量為0.628 0～0.855 0，平均為0.702 2，說明此品種具有較豐富的遺傳多樣性；哈迪-溫伯格平衡檢驗(yàn)結(jié)果表明，群體處于不平衡狀態(tài)；群體平均固定系數(shù)為-0.041 0，說明該群體存在雜合子過?，F(xiàn)象，有較大的選擇壓力；瓶頸效應(yīng)分析表明，群體經(jīng)歷了瓶頸效應(yīng)；根據(jù)連鎖不平衡方法計(jì)算有效群體大小為64.2?！窘Y(jié)論】 易捕鯉具有較高的遺傳異質(zhì)性和較大的選育空間，但由于人工選擇的作用，群體處于不平衡狀態(tài)，對(duì)品種的長(zhǎng)期發(fā)展極為不利，因而在接下來的工作中應(yīng)該注重加強(qiáng)保種策略的研究，從而保持其豐富的遺傳多樣性和優(yōu)良的經(jīng)濟(jì)性狀。

[關(guān)鍵詞]易捕鯉；微衛(wèi)星；遺傳結(jié)構(gòu)；保種

鯉(CyprinuscarpioL.)作為主要淡水養(yǎng)殖魚種之一，在全球的淡水養(yǎng)殖魚類中占有較大的比重。由于鯉是底層魚類，具有較強(qiáng)的逃網(wǎng)能力和較低的回捕率，因而其在一些湖泊、水庫(kù)中較難捕獲，從而導(dǎo)致捕獲量較低。為了提高鯉的養(yǎng)殖產(chǎn)量，黑龍江水產(chǎn)研究所利用高起捕率的大頭鯉(Cyprinuspellegrinipellegrini)、德國(guó)鏡鯉(CyprinuscarpioL.mirror)和黑龍江野鯉(Cyprinuscarpiohaematopterus)3個(gè)品種作為親本進(jìn)行雜交，從而選育出一個(gè)高起捕率的鯉新品——易捕鯉。通過在池塘、水庫(kù)、室內(nèi)水族箱中對(duì)易捕鯉起捕率進(jìn)行測(cè)定，結(jié)果表明，其兩網(wǎng)起捕率平均可達(dá)到95%以上[1]。由此可見易捕鯉不僅可以在大水面進(jìn)行養(yǎng)殖而且具有起捕率高的優(yōu)良特性。為了全面深入地了解易捕鯉的種質(zhì)特性，防止該品種優(yōu)良性狀退化，對(duì)該品種種質(zhì)資源的管理顯得極為重要。

微衛(wèi)星標(biāo)記一般由1～6 bp的核心序列重復(fù)串聯(lián)而成，具有豐富的多態(tài)性和信息量，廣泛地應(yīng)用于育成品種遺傳結(jié)構(gòu)的研究中，可為選育品種種質(zhì)資源長(zhǎng)期利用和管理提供支撐，如賈智英等[2]、石連玉等[3]利用微衛(wèi)星標(biāo)記對(duì)人工育成品種松浦鯉和松荷鯉保種群體遺傳多樣性和遺傳結(jié)構(gòu)進(jìn)行研究，結(jié)果表明，松浦鯉和松荷鯉都具有較高的遺傳多樣性，2個(gè)育成品種雖然都經(jīng)歷了多代的人工選擇，但其遺傳多樣性仍處于中上等水平，為二者的群體遺傳研究奠定了堅(jiān)實(shí)的基礎(chǔ)。為了促進(jìn)易捕鯉優(yōu)良性狀的穩(wěn)定保存和種質(zhì)資源的可持續(xù)利用，本試驗(yàn)用30對(duì)多態(tài)性較高的微衛(wèi)星標(biāo)記對(duì)易捕鯉保種群體的遺傳多樣性和遺傳結(jié)構(gòu)進(jìn)行分析，以期為易捕鯉種質(zhì)資源的有效保存和長(zhǎng)期監(jiān)測(cè)提供理論依據(jù)。

1材料與方法

1.1材料

易捕鯉(F5代)96尾，養(yǎng)殖于中國(guó)水產(chǎn)科學(xué)研究院黑龍江水產(chǎn)研究所呼蘭試驗(yàn)站，編號(hào)1-1～1-96，取每尾魚的鰭條置于體積分?jǐn)?shù)75%的酒精之中，于-20 ℃的冰箱中保存。

本試驗(yàn)所用的主要試劑有瓊脂糖、聚丙烯、硝酸銀、氫氧化鈉甲醇溶液、磁珠法組織DNA提取試劑盒(廣州欣研技術(shù)有限公司，產(chǎn)品號(hào)(MD5112-03)、2×Es Tag Master Mix(康為世紀(jì)生物公司，產(chǎn)品號(hào)CW0690) 以及滅菌雙蒸餾水。本試驗(yàn)所用主要儀器有磁珠法組織DNA提取儀、微量紫外分光光度儀、PCR儀以及瓊脂糖凝膠和聚丙烯凝膠電泳儀等。

1.2方法

1.2.1DNA提取和引物的選取采用磁珠法提取易捕鯉鰭條DNA，經(jīng)1%瓊脂糖凝膠電泳和紫外分光光度法檢測(cè)其純度和濃度后，稀釋成50 ng/μL的溶液備用。微衛(wèi)星引物參見文獻(xiàn)[4-6]，具體引物信息見表1，由上海生工生物工程有限公司合成。

表 1　30 對(duì)微衛(wèi)星的信息

續(xù)表 1　Continued table 1

1.2.2PCR擴(kuò)增PCR反應(yīng)體系為25 μL：模板DNA (50 ng/μL)0.6～1 μL，2×Es Tag Master Mix 12.5 μL，上、下游引物(10 μmol/μL)各1 μL，用滅菌雙蒸餾水補(bǔ)足體積。反應(yīng)條件為：95 ℃預(yù)變性5 min；94 ℃變性30 s，最適溫度退火30 s，72 ℃延伸30 s，23～26個(gè)循環(huán)；72 ℃總延伸5 min，4 ℃保存。

1.2.3電泳及染色PCR 產(chǎn)物用8.0%的非變性聚丙烯酰胺凝膠電泳分離。電泳結(jié)束，用硝酸銀染色后，拍照記錄電泳結(jié)果。

1.2.4數(shù)據(jù)分析利用PopGene(Version 3.2)[7]軟件分別對(duì)觀測(cè)等位基因數(shù)(Observed number of alleles，A)、有效等位基因數(shù)(Effective number of alleles，Ae)、等位基因頻率(Allele frequency)、觀測(cè)雜合度(Observed heterozygosity，Ho)、期望雜合度(Expected heterozygosity，He)、連鎖不平衡等遺傳變異參數(shù)進(jìn)行統(tǒng)計(jì)。使用Genepop version 4.0軟件的Markov chain法進(jìn)行Hardy-Weinberg平衡檢驗(yàn)[8](哈迪-溫伯格平衡檢驗(yàn)概率結(jié)果用P表示)，并對(duì)偏離平衡位點(diǎn)進(jìn)行雜合過度與缺失分析，計(jì)算固定系數(shù)(FIS)。應(yīng)用NeEstimator version 1.3軟件根據(jù)基因型連鎖不平衡原理進(jìn)行有效群體大小估算[9]。應(yīng)用kingroup v2軟件計(jì)算全同胞率和半同胞率[10]。應(yīng)用Bottleneck version 1.2.02 軟件，基于無限等位基因模型(IAM)、雙突變模型(TPM) 和逐步突變模型(SMM)的假設(shè)，分析標(biāo)記檢驗(yàn)(Sigh test)和Wilcoxon 標(biāo)記秩檢驗(yàn)(Wilcoxon sigh-rank test)結(jié)果的瓶頸效應(yīng)[11]。應(yīng)用phylip軟件計(jì)算個(gè)體間遺傳距離[12]，用MEGA 4.0軟件構(gòu)建個(gè)體聚類圖[13]。

2結(jié)果與分析

2.1易捕鯉的微衛(wèi)星擴(kuò)增結(jié)果和遺傳多樣性

試驗(yàn)選取的30對(duì)微衛(wèi)星標(biāo)記在易捕鯉群體中均擴(kuò)增出了具有特異性的清晰條帶，圖1為位點(diǎn)HLJ3410和HLJ3857的部分?jǐn)U增結(jié)果。

圖 1位點(diǎn)HLJ3410(A)和HLJ3857(B)對(duì)易捕鯉的部分?jǐn)U增結(jié)果

Fig.1Yibu carp partial amplified results of HLJ3410(A) and HLJ3857(B)

易捕鯉群體的遺傳變異參數(shù)見表2。由表2可知，所有微衛(wèi)星座位共檢測(cè)到154個(gè)等位基因，平均等位基因數(shù)為5.133個(gè)，單個(gè)標(biāo)記檢測(cè)到的等位基因數(shù)為3～9個(gè)；有效等位基因數(shù)為2.461 5～7.486 6，平均為4.060 7；觀測(cè)雜合度為0.260 4～0.927 1，平均為0.703 5；期望雜合度為0.596 9～0.871 0，平均為0.740 0；多態(tài)信息含量為0.628 0～0.855 0，平均為0.702 2；本試驗(yàn)所選取的所有位點(diǎn)都具有高度多態(tài)性(P≥0.5)。易捕鯉的等位基因頻率分布屬于非“L”形(圖2)。

表 2　30對(duì)微衛(wèi)星遺傳標(biāo)記基因位點(diǎn)對(duì)易捕鯉的檢測(cè)結(jié)果

注：P>0.05，位點(diǎn)處于Hardy-Weinberg平衡狀態(tài)；P<0.05，位點(diǎn)處于不平衡狀態(tài)；P<0.01，位點(diǎn)處于極不平衡狀態(tài)。

Note：P>0.05,sites in Hardy-Weinberg equilibrium;P<0.05,sites in an unbalanced state;P<0.01,sites in a very unbalanced state.

圖 2　易捕鯉的等位基因頻率分布

2.2易捕鯉的Hardy-Weinberg 平衡檢驗(yàn)和FIS值

表2表明，有7個(gè)位點(diǎn)處于Hardy-Weinberg 平衡狀態(tài)，23個(gè)位點(diǎn)處于不平衡狀態(tài)(P≤0.05)，其中有17個(gè)位點(diǎn)處于極平衡狀態(tài)(P≤0.01)；易捕鯉群體內(nèi)固定系數(shù)平均為-0.041 0。

2.3易捕鯉有效群體大小和同胞家系的檢測(cè)

據(jù)遺傳連鎖不平衡得到的有效群體大小為64.2(CI95%=59.5～69.5)。用kingoup V2軟件中的SDR方法檢測(cè)易捕鯉群體的同胞率，結(jié)果表明易捕鯉群體中共檢測(cè)到20個(gè)全同胞家系和26個(gè)半同胞家系，全同胞率和半同胞率分別為41.67%和54.16%，96尾易捕鯉群體的同胞家系(全同胞與半同胞之和)含有86個(gè)個(gè)體，同胞率為89.58%。

2.4易捕鯉的瓶頸效應(yīng)分析

基于 IAM、TPM 和 SMM 假設(shè)的標(biāo)記檢驗(yàn)和 Wilcoxon 標(biāo)記秩檢驗(yàn)的總體結(jié)果見表3。由表3可知，易捕鯉在IAM、TPM假設(shè)下極顯著偏離突變-漂移模型(P<0.000 1)，極顯著的雜合過剩，在SMM假設(shè)的Wilcoxon 標(biāo)記秩檢驗(yàn)下極顯著偏離突變-漂移模型(P<0.000 1)。從微衛(wèi)星位點(diǎn)(表4)上來看，在 IAM假設(shè)下，25個(gè)位點(diǎn)偏離突變-漂移平衡(P<0.05)，其中18個(gè)位點(diǎn)極顯著偏離(P<0.01) ；在 TPM 假設(shè)下，23個(gè)位點(diǎn)偏離突變-漂移平衡(P<0.05)，其中 12個(gè)位點(diǎn)極顯著偏離(P<0.01) ；在 SMM 假設(shè)下，14個(gè)位點(diǎn)偏離突變-漂移平衡(P<0.05)，其中 5個(gè)位點(diǎn)極顯著偏離(P<0.01)。因此，從微衛(wèi)星變異來看，該品種部分位點(diǎn)偏離突變-漂移平衡，說明群體經(jīng)歷了瓶頸效應(yīng)。

表 3　不同檢驗(yàn)方法對(duì)易捕鯉品種突變-漂移平衡的分析結(jié)果

注(Note):**.P<0.01。

表 4　易捕鯉30個(gè)微衛(wèi)星位點(diǎn)的瓶頸效應(yīng)分析

注:Heq．平均期望雜合度;P．概率。

Note:Heq．Expected heterozygosity;P.probability．

2.5易捕鯉個(gè)體間遺傳距離和個(gè)體聚類分析

由表5可知，易捕鯉個(gè)體之間遺傳距離(D)為0.2～0.8，平均值為0.572 6，主要分布于0.4～0.7，約占94%。由圖3可知,96尾易捕鯉個(gè)體形成2個(gè)分支，其中個(gè)體1-6、1-84～1-96號(hào)在一個(gè)分支上，其他個(gè)體在另一個(gè)分支上。

表 5　易捕鯉個(gè)體間遺傳距離分布

圖 3　易捕鯉個(gè)體的UPGMA聚類

3討論

3.1易捕鯉的遺傳多樣性分析

遺傳多樣性不僅是種質(zhì)評(píng)價(jià)的重要指標(biāo)，同時(shí)也是生物多樣性的核心。對(duì)于一個(gè)物種而言，豐富的遺傳多樣性可增強(qiáng)其對(duì)環(huán)境的適應(yīng)能力和生存能力，反之則可能導(dǎo)致其生存能力和適應(yīng)能力降低，最終導(dǎo)致物種退化[14]。對(duì)于人工選育品種而言，在選育過程中極易出現(xiàn)近交和遺傳漂變，因而對(duì)選育群體的遺傳多樣性進(jìn)行評(píng)估具有重要意義。有效等位基因數(shù)、期望雜合度以及多態(tài)信息含量是反映群體遺傳多樣性的重要指標(biāo)，其值越高，表明該群體的遺傳多樣性越高；反之，說明群體的遺傳一致性越高[15]。本試驗(yàn)結(jié)果表明，易捕鯉F5代平均有效等位基因數(shù)為4.060 7，平均期望雜合度為0.740 0，多態(tài)信息含量0.702 2，與池喜峰等[16]對(duì)易捕鯉F1-F4代的研究結(jié)果比較得出：隨著易捕鯉選育世代的增加，其平均有效等位基因數(shù)逐代降低，從F1到F5易捕鯉平均有效等位基因數(shù)從7.761 9下降到4.060 7。趙廣泰等[17]對(duì)大黃魚連續(xù)4代選育群體的研究表明，隨著選育的進(jìn)行，F(xiàn)1到F4平均等位基因數(shù)從5.462下降到4.308。鄭荷子等[18]對(duì)翹嘴鱖連續(xù)4代選育群體遺傳多樣性及遺傳結(jié)構(gòu)的研究表明，F(xiàn)1到F4的平均等位基因數(shù)從5.14下降到 2.57。這說明，對(duì)于人工選育品種而言，連續(xù)多代的選育所面臨的人工壓力和環(huán)境壓力勢(shì)必使其遺傳參數(shù)發(fā)生變化，并且隨著選育世代的增加，群體遺傳多樣性有下降的趨勢(shì)，表明人工選擇壓力對(duì)其遺傳多樣性造成了影響。因此在后續(xù)的選育過程中，必須要完善選育措施，增加繁育親本的數(shù)量，避免近交的影響。本試驗(yàn)多態(tài)信息含量和觀測(cè)雜合度結(jié)果顯示，易捕鯉雖然經(jīng)過連續(xù)5代的選育，但其多態(tài)信息含量和期望雜合度并沒有隨著選育世代的增加而逐代下降，與目前報(bào)道的人工育成品種相比，易捕鯉的Ae和He與松浦鯉[2]相近，高于松浦紅鏡鯉[19]，低于松荷鯉[3]和建鯉[20],說明與其他人工選育品種相比，易捕鯉保種群體具有較高的遺傳多樣性水平，選育空間較大。

3.2易捕鯉的遺傳結(jié)構(gòu)分析

微衛(wèi)星具有高度多態(tài)性和共顯性孟德爾遺傳的特點(diǎn)，因此用其對(duì)遺傳結(jié)構(gòu)進(jìn)行研究結(jié)果更加精確，尤其是對(duì)小群體以及群體近期的變異，如瓶頸效應(yīng)、哈代-溫伯格平衡檢驗(yàn)等的檢測(cè)，效果尤為突出。

對(duì)于養(yǎng)殖群體而言，當(dāng)種群有效數(shù)量急劇減少時(shí)，群體有效等位基因丟失，會(huì)導(dǎo)致群體經(jīng)歷瓶頸效應(yīng)，近親交配幾率將大幅上升，從而導(dǎo)致種群遺傳性能衰退[21]。因此瓶頸效對(duì)保種策略的制定具有重要的影響。賈智英等[2]的研究表明，松浦鯉等位基因頻率分布于0.1～0.2的個(gè)體所占比例高于等位基因頻率分布<0.1的個(gè)體，屬于非“L”形，說明該品種經(jīng)歷了瓶頸效應(yīng)。本研究易捕鯉等位基因頻率分布也呈非“L”形，說明其也經(jīng)歷了瓶頸效應(yīng)，這與在IAM、TPM和SMM模型下的假設(shè)檢驗(yàn)結(jié)果相同。在目前已經(jīng)報(bào)道的鯉魚育成品種之中，松浦鯉[2]、松荷鯉[3]、松浦紅鏡鯉[19]等位基因頻率分布都屬于非“L”形，說明這些選育品種均經(jīng)歷了瓶頸效應(yīng)。在本研究易捕鯉等位基因頻率分布不僅屬非“L”形，而且呈現(xiàn)出接近正態(tài)分布曲線的形狀，說明易捕鯉群體顯著地經(jīng)歷了瓶頸效應(yīng)，出現(xiàn)此現(xiàn)象的原因可能是由于在選育的過程中所選取的親本數(shù)量過少，未能剔除親緣關(guān)系較近的個(gè)體，使得部分稀有等位基因丟失。所以在對(duì)易捕鯉的保種過程中，應(yīng)該充分考慮其遺傳背景，嚴(yán)格控制近親繁殖，以保持群體豐富的遺傳多樣性。

在哈迪-溫伯格平衡檢驗(yàn)中，易捕鯉大部分(77%)位點(diǎn)處于不平衡狀態(tài)。易捕鯉作為人工選育的品種，由于人工選擇、近交、遺傳漂變和親本數(shù)量受限等因素的影響，導(dǎo)致大部分位點(diǎn)顯著偏離哈迪-溫伯格平衡，這與其他學(xué)者研究得出的人工選育品種大部分位點(diǎn)(高于50%)偏離哈迪-溫伯格平衡的結(jié)論[2,22]是一致的。比較本試驗(yàn)結(jié)果與遲喜峰等[16]對(duì)易捕鯉F1-F4的研究結(jié)果發(fā)現(xiàn)，顯著偏離平衡的位點(diǎn)由F1的7個(gè)上升到本研究F5的23個(gè)，這表明易捕鯉經(jīng)過人工定向5代的選育，群體遺傳結(jié)構(gòu)發(fā)生了變化，隨選育代數(shù)的增加，群體基因頻率的變化加劇，使一些低頻等位基因丟失。本試驗(yàn)易捕鯉平均FIS<0，表明該群體存在雜合子過剩現(xiàn)象，這可能是由于人工選擇壓力所造成的。從本研究結(jié)果看，易捕鯉保種群體遺傳結(jié)構(gòu)極不平衡，因而應(yīng)當(dāng)制定適宜的保種策略，保持該品種的遺傳性能。

對(duì)于保種群體的研究，主要分析的參數(shù)是個(gè)體的遺傳距離和個(gè)體聚類，根據(jù)遺傳距離和聚類圖分析結(jié)果，在后續(xù)的保種和繁殖工作中去掉親緣關(guān)系較近的個(gè)體，以降低其遺傳相似性。Barker[23]認(rèn)為，遺傳距離可以作為決策保種計(jì)劃時(shí)研究群體結(jié)構(gòu)的重要指標(biāo)。遺傳結(jié)構(gòu)不僅是研究群體遺傳多樣性的基礎(chǔ)，而且可以用來描述群體的系統(tǒng)分化。目前對(duì)于人工選育品種而言，松浦鯉[2]的遺傳距離主要分布在0.5～0.8，松浦紅鏡鯉[19]和松荷鯉[3]的遺傳距離主要分布0.4～0.6，本試驗(yàn)易捕鯉的遺傳距離主要分布在0.4～0.7。一般認(rèn)為群體分化的時(shí)間越短，遺傳距離越小，從以上的研究結(jié)果來看，目前對(duì)于人工選育的鯉魚品種而言，其遺傳距離處于中等水平，群體分化時(shí)間較長(zhǎng);另外對(duì)易捕鯉同胞率的計(jì)算結(jié)果顯示,其同胞率為89.58%，說明群體內(nèi)有較高的近交壓力。這提示在接下來的選育過程中，應(yīng)盡量選擇個(gè)體間遺傳距離較大的個(gè)體進(jìn)行配組，從而最大限度地避免近親交配，以保持群體豐富的遺傳多樣性。

3.3易捕鯉的有效群體分析

有效群體分析是種質(zhì)資源保存研究中不可或缺的內(nèi)容[24-25]，保持足夠大小的群體既可以降低近交幾率，又可以保持群體的遺傳多樣性。本試驗(yàn)采用England等[26]研究的連鎖不平衡方法對(duì)易捕鯉有效群體大小進(jìn)行評(píng)估，結(jié)果顯示其有效群體大小為64.2尾，本研究所采集樣本尾數(shù)為96尾，滿足連鎖不平衡分析的樣本數(shù)量，估算的有效群體大小能夠準(zhǔn)確說明該品種的有效群體數(shù)量。Nei等[27]研究表明，群體數(shù)量為有效群體的4～10倍時(shí)才能保持群體遺傳多樣性的穩(wěn)定性，因此保持易捕鯉種質(zhì)穩(wěn)定的群體數(shù)量是256.8～642尾，在該范圍內(nèi)數(shù)量越大越有利于種質(zhì)穩(wěn)定保存。

[參考文獻(xiàn)]

[1]李池陶,石連玉,杜海清,等.易捕鯉新品種起捕率的測(cè)定：2010年中國(guó)水產(chǎn)學(xué)會(huì)學(xué)術(shù)年會(huì)論文摘要集 [M].西安:中國(guó)水產(chǎn)學(xué)會(huì),2011:57.

Li C T,Shi L Y,Du H Q,et al.Determination of the harvesting rate for the Yibu common carp：2010 China Fisheries Society Annual Conference Abstract Book [M].Xi’an:China Society of Fisheries,2011:57.(in Chinese)

[2]賈智英,李飛,李池陶,等.松浦鯉保種群體遺傳結(jié)構(gòu)研究 [J].華北農(nóng)學(xué)報(bào),2011,26(增刊):1-8.

Jia Z Y,Li F,Li C T,et al.Study on the genetic diversity of the conservation population in Song pu common carp [J].Acta Agriculturae Boreali-Sinica,2011,26(Supple):1-8.(in Chinese)

[3]石連玉,李飛,賈智英,等.選育品種松荷鯉遺傳結(jié)構(gòu)研究 [J].海洋湖沼通報(bào),2012,34(4):104-112.

Shi L Y,Li F,Jia Z Y,et al.Study on genetic diversity and structure in songhe carp [J].Transactions of Oceanology and Limnology,2012,34(4):104-112.(in Chinese)

[4]Zhang Y,Xu P,Lu C,et al.Genetic linkage mapping and analysis of muscle fiber-related QTLs in common carp (CyprinuscarpioL.) [J].Marine Biotechnology,2011,13(3):376-392.

[5]Crooijmans R,Poel J J,Groenen M,et al.Microsatellite markers in common carp (CyprinuscarpioL.) [J].Animal Genetics,1997,28(2):129-134.

[6]Tanck M W T,Palstra A P,Weerd M,et al.Segregation of microsatellite alleles and residual heterozygosity at single loci in homozygous androgenetic common carp (CyprinuscarpioL.) [J].Genome,2001,44(5):743-751.

[7]Yeh F C,Yang R C,Boyle T.Popgene version 1.31 quick user guide [M].Canada:University of Alberta and Center for International Forestry Research,1999.

[8]Rousset F.Genepop’ 007:A complete re-implementation of the genepop software for Windows and Linux [J].Molecular Ecology Resources,2008,8(1):103-106.

[9]Peel D,Ovenden J R,Peel S L.Ne estimator:software for estimating effective population size [M].Version 1.3.Queensland Government:Department of Primary Industries and Fisheries,2004.

[10]Konovalov D A,Manning C,Henshaw M T.KINGROUP:A program for pedigree relationship reconstruction and kin group assignments using genetic markers [J].Molecular Ecology Notes,2004,4(4):779-782.

[11]Cornuet J M,Luikart G.Description and power analysis of two tests for detecting recent population bottlenecks from allele frequency data [J].Genetics,1996,144(4):2001-2014.

[12]Felsenstein J.{PHYLIP}:Phylogenetic inference package [M].Version 3.5 c.Washington:The Department of Genome Science and the Department of Biology at the University of Washington,1995.

[13]Tamura K,Dudley J,Nei M,et al.MEGA4:Molecular evolutionary genetics analysis (MEGA) software version 4.0 [J].Molecular Biology and Evolution,2007,24(8):1596-1599.

[14]Na-Nakorn U,Kamonrat W,Ngamsiri T.Genetic diversity of walking catfish,Clarias macrocephalus in Thailand and evidence of genetic introgression from introduced farmedC.gariepinus[J].Aquaculture,2004,240(1):145-163.

[15]蒙子寧,莊志猛,金顯仕,等.黃海和東海小黃魚遺傳多樣性的 RAPD 分析 [J].生物多樣性,2003(3):197-203.

Meng Z N,Zhuang Z M,Jin X S,et al.Genetic diversity in small yellow croaker (Pseudosciaenapolyactis) by PAPD analysis [J].Biodiversity Science,2003(3):197-203.(in Chinese)

[16]池喜峰,賈智英,李池陶,等.鯉易捕性狀選育群體不同世代微衛(wèi)星分析 [J].上海海洋大學(xué)學(xué)報(bào),2010,19(3):308-313.

Chi X F,Jia Z Y,Li C T,et al.SSR analysis on selected high-harvesting populations of common carp (Cyprinuscarpio) [J].Journal of Shanghai Ocean University,2010,19(3):308-313.(in Chinese)

[17]趙廣泰,劉賢德,王志勇,等.大黃魚連續(xù)4代選育群體遺傳多樣性與遺傳結(jié)構(gòu)的微衛(wèi)星分析 [J].水產(chǎn)學(xué)報(bào),2010,34(4):500-508.

Zhao G T,Liu X D,Wang Z Y,et al.Genetic structure and genetic diversity analysis of four consecutive breeding generations of large yellow croaker(Pseudosciaenacrocea) using microsatellite markers [J].Journal of Fisheries of China,2010,34(4):500-508.(in Chinese)

[18]鄭荷子,易提林,梁旭方,等.翹嘴鱖連續(xù)4代選育群體遺傳多樣性及遺傳結(jié)構(gòu)分析 [J].淡水漁業(yè),2013,43(6):8-12.

Zheng H Z,Yi T L,Liang X F,et al.Genetic structure and genetic diversity analysis of four consecutive breeding generations ofSinipercachuatsi[J].Freshwater Fisheries,2013,43(6):8-12.(in Chinese)

[19]李盛文,賈智英,柏盈盈,等.松浦紅鏡鯉保種群體的遺傳結(jié)構(gòu) [J].中國(guó)水產(chǎn)科學(xué),2014,21(1):67-74.

Li S W,Jia Z Y,Bai Y Y,et al.The genetic diversity of a conservation population of red mirror carp,Cyprinuscarpiovar.Songpu [J].Journal of Fishery Sciences of China,2014,21(1):67-74.(in Chinese)

[20]李建林,唐永凱,李紅霞,等.利用微衛(wèi)星標(biāo)記分析建鯉種質(zhì)資源的研究 [J].中國(guó)農(nóng)學(xué)通報(bào),2012,28(32):107-112.

Li J L,Tang Y K,Li H X,et al.Germplasm resources analysis inCyprinuscarpiovar.Jan by microsatellite markers [J].Chinese Agricultural Science Bulletin,2012,28(32):107-112.(in Chinese)

[21]孫莉,楊國(guó)梁,王軍毅,等.利用微衛(wèi)星標(biāo)記分析長(zhǎng)鰭鯉、錦鯉和龍鳳鯉的遺傳多樣性 [J].揚(yáng)州大學(xué)學(xué)報(bào):農(nóng)業(yè)與生命科學(xué)版,2009,30(3):35-40.

Sun L,Yang G L,Wang J Y,et al.Study on genetic diversity of Long fin carp,Koi carp and Longfeng carp by microsatellite markers [J].Journal of Yangzhou University:Agricultural and Life Science Edition,2009,30(3):35-40.(in Chinese)

[22]Norris A T,Bradley D G,Cunningham E P.Microsatellite genetic variation between and within farmed and wild Atlantic salmon (Salmosalar) populations [J].Aquaculture,1999,180(3):247-264.

[23]Barker J S F.A global protocol for determining genetic distances among domestic livestock breeds [M].Canada:Proceedings of the 5th World Congress on Genetics Applied to Livestock Production,Guelph and Ontario,1994,21:501-508.

[24]Mace G M,Lande R.Assessing extinction threats:Toward a reevaluation of IUCN threatened species categories [J].Conservation Biology,1991,5(2):148-157.

[25]Frankham R.Effective population size/adult population size ratios in wildlife:A review [J].Genetical Research,1995,66(2):95-107.

[26]England P R,Cornuet J M,Berthier P,et al.Estimating effective population size from linkage disequilibrium:Severe bias in small samples [J].Conservation Genetics,2006,7(2):303-308.

[27]Nei M,Chakraborty R,Fuerst P A.Infinite allele model with varying mutation rate [J].Proceedings of the National Academy of Sciences,1976,73(11):4164-4168.

Genetic diversity of conservation population in Yibu F5common carp

ZHAO Xin-chun1,2,JIA Zhi-ying1,LI Sheng-wen1,2,LI Chi-tao1,SHI Lian-yu1，2

(1HeilongjiangRiverFisheriesResearchInstitute,ChineseAcademyofFisherySciences,Harbin,Heilongjiang150070,China;2CollegeofFisheriesandLife,ShanghaiOceanUniversity,Shanghai201306,China)

Abstract:【Objective】 Yibu common carp is a new breeding variety.In order to maintain the genetic diversity,its genetic structure was studied．【Method】 The sequences of 30 microsatellite loci were amplified by PCR using genetic DNA of 96 Yibu common carps as templates.According to the genotype,biological software was used for analysis of genetic diversity. 【Result】 From the 96 sampled individuals,154 alleles were detected with 3-9 alleles at each locus and the average number of 4.060 7.The expected heterozygosity was 0.596 9-0.871 0 with the average of 0.740 0.The polymorphic information content (PIC) was 0.628 0-0.855 0 with the average of 0.702 2.The results indicated rich polymorphism information content and large genetic diversity in the population.Hardy-weinberg equilibrium analysis showed that the population was in disequilibrium mode,indicating the artificial selection had a great effect on population.The average fixation index (FIS) was -0.041 0,showing that heterozygote excess existed in the population.Bottleneck analysis illustrated that the population had experienced bottleneck effect recently.The effective population size was 64.2 based on linkage disequilibrium method.【Conclusion】 Yibu common carp had rich genetic diversity,but the population structure was disequilibrium and bottleneck effect happened.Therefore,it is necessary to strengthen protection measure to maintain its rich diversity and economic traits.

Key words:Yibu common carp;microsatellite markers;genetic structure;reservation

DOI：網(wǎng)絡(luò)出版時(shí)間:2016-02-0209:3710.13207/j.cnki.jnwafu.2016.03.005

[收稿日期]2014-07-16

[基金項(xiàng)目]國(guó)家“十二五”科技支撐計(jì)劃項(xiàng)目“重要水產(chǎn)種質(zhì)資源群體構(gòu)建及遺傳評(píng)價(jià)方法研究”(2012BAD26B01)；國(guó)家大宗淡水魚類產(chǎn)業(yè)技術(shù)體系北方鯉魚種質(zhì)資源與育種項(xiàng)目(CARS-46-02)

[作者簡(jiǎn)介]趙新春(1988-)，女，內(nèi)蒙古赤峰人，在讀碩士，主要從事水產(chǎn)遺傳育研究。E-mail：zhaoxc0716@163.com[通信作者]石連玉(1960-)，男，遼寧西豐人，研究員，碩士生導(dǎo)師，主要從事水產(chǎn)遺傳育種研究。E-mail:sly2552@aliyun.com

[中圖分類號(hào)]S917；Q959.46+8

[文獻(xiàn)標(biāo)志碼]A

[文章編號(hào)]1671-9387(2016)03-0028-09