養(yǎng)殖大鱗鲃發(fā)病死亡F3個體的微衛(wèi)星標記分析

常玉梅，魯翠云，苗建發(fā)，李崇文，梁利群

（1.中國水產科學研究院黑龍江水產研究所，黑龍江 哈爾濱 150070；2.天津市天祥水產有限責任公司，天津 310500）

養(yǎng)殖大鱗鲃發(fā)病死亡F3個體的微衛(wèi)星標記分析

常玉梅1，魯翠云1，苗建發(fā)2，李崇文2，梁利群1

（1.中國水產科學研究院黑龍江水產研究所，黑龍江 哈爾濱 150070；2.天津市天祥水產有限責任公司，天津 310500）

采用24個已報道的多態(tài)微衛(wèi)星標記對天津市天祥水產有限責任公司養(yǎng)殖的大鱗鲃Barbus capito F2親本及F3成活和病發(fā)死亡個體進行基因分型，結果有8個標記無多態(tài)，13個標記用以分析本研究群體的遺傳多樣性、親緣關系及遺傳結構。遺傳參數(shù)分析顯示，13個標記在87個大鱗鲃個體中分別檢測到2～7個等位基因，平均等位基因數(shù)（4.27～4.53）極顯著高于平均有效等位基因數(shù)（3.09～3.22）（P＜0.01）；觀察雜合度（0.58～0.64）整體低于期望雜合度（0.64～0.65），多態(tài)位點比例明顯下降（8/24），表明大鱗鲃F2及F3個體的純合度增加，遺傳多樣性下降。聚類分析顯示，F(xiàn)2親本與成活子代親緣關系更近，聚在一起；而死亡子代單獨聚為1支。遺傳結構分析顯示，親本與子代清晰地劃分成兩大類群，多數(shù)親本與多數(shù)成活個體歸為第一類群（劃分概率為51.8%～97.7%），少數(shù)親本與多數(shù)死亡個體歸為第二類群（劃分概率為52.9%～98.4%），這與親緣關系分析的結果一致，表明少數(shù)親本是F3死亡個體遺傳組成的直接貢獻者。本研究初步認定，由于幾個親本的近親交配，子代基因組純合度增加，遺傳多樣性下降，抗逆性下降，這可能是F3出現(xiàn)不明緣由病發(fā)死亡的主要原因之一。

人工養(yǎng)殖大鱗鲃；微衛(wèi)星標記；遺傳多樣性；遺傳結構；近交衰退

大鱗鲃Barbus capito隸屬于鯉科Cyprinidae，鲃亞科Barbinae，鲃屬，主要分布在咸海（Aral Sea）和里海（Caspian Sea）兩大咸水湖，具有食性廣、生長速度快、肉質鮮美、耐鹽堿、適應性強等優(yōu)良特點，是湖區(qū)主要的名貴經濟魚類之一[1,2]。我國內陸蘊藏有豐富的鹽堿水域，但大部分礦化度高，水體離子復雜等，仍處于荒蕪狀態(tài)。雖有小部分鹽堿水體可被漁業(yè)利用，但適于鹽堿水養(yǎng)殖的種類少，特別缺乏大型經濟魚類，其漁業(yè)生產力低下。2003年黑龍江水產研究所首次從烏茲別克斯坦咸海引進大鱗鲃魚種。經過十幾年的研究，基本掌握大鱗鲃的全人工繁殖技術[3]，詳細調查研究其繁殖生物學[4,5]、耐鹽堿特性[6-8]、生長和養(yǎng)殖模式[9]、營養(yǎng)需求[10]和藥物毒理[11]等，為大鱗鲃在國內的大規(guī)模推廣養(yǎng)殖積累了豐富的科研數(shù)據(jù)。目前，大鱗鲃已成功在肇東、天津、保定和青島等低洼鹽堿水域進行推廣養(yǎng)殖，經濟效益可觀，已得到市場認可[12]。

然而，由于建群數(shù)量較小，近交衰退的弊端在大鱗鲃養(yǎng)殖過程中逐漸顯現(xiàn)。如后代體色分化，養(yǎng)殖群體F2和F3出現(xiàn)白色和紅色個體，且群體生長速度下降，體色變異個體的生長速度顯著低于未變異個體[13]；病害頻發(fā)，特別是2016年春季天津市天祥水產有限責任公司養(yǎng)殖的F3個體出現(xiàn)不明原因的大量死亡，死亡個體主要表現(xiàn)為身體極度貧血，造成一定的經濟損失。

分子標記輔助育種技術是解決養(yǎng)殖魚類近交衰退的有效途徑，其中微衛(wèi)星分子標記的多態(tài)性高、共顯性分離、操作簡單等優(yōu)勢，是多數(shù)養(yǎng)殖魚類分子標記輔助育種的首選標記[14]。Geng等[15]利用磁珠富集法克隆大鱗鲃的微衛(wèi)星標記，并應用在大鱗鲃F1的繁殖中；隨后，魯翠云等[16]利用24個多態(tài)微衛(wèi)星標記檢測F1的遺傳多樣性，認為F1群體的遺傳多樣性處于中度多態(tài)水平，遺傳結構相對平衡。但后續(xù)有關F2及F3的遺傳多樣性篩查未見相關報道。本研究針對天津市天祥水產有限責任公司2016年養(yǎng)殖的大鱗鲃F3病發(fā)死亡、成活個體及其繁殖親本，采用微衛(wèi)星標記進行基因分型及數(shù)據(jù)統(tǒng)計，以檢測F2繁殖子代的遺傳多樣性，探討F3個體病發(fā)死亡的遺傳因素，為大鱗鲃的規(guī)?；庇坝N提供理論支撐。

1 材料與方法

1.1 實驗魚來源

實驗魚采自天津市天祥水產有限公司，親本（3齡）20尾，雌雄各10尾（來自黑龍江水產研究所繁殖的F2），子代為天津天祥水產有限公司自繁F3，67尾（2齡），其中不明原因病發(fā)死亡個體34尾，成活個體33尾。

1.2 基因組DNA的提取

剪取少量95%酒精保存的大鱗鲃鰭條至離心管中，加入600μL DNA提取裂解液（成分：1MTris（pH8.0）、5M NaCl、0.5M EDTA（pH8.0）、10%SDS 和200μg/mL Proteinase K）于 55℃消化 2～4h；加入等體積的酚 /氯仿混合液（1∶1）抽提 1次，12 000r/min室溫離心10min，吸取上清（大約550μL）；加入2倍體積的無水乙醇沉淀，12 000r/min室溫離心10min；加入500μL70%酒精洗滌，12 000r/min室溫離心5min；吸凈殘留酒精后，室溫干燥10min，加入50～100μL 0.1×TE溶解。提取后的DNA采用Nanodrop8000測定濃度，1%的瓊脂糖凝膠電泳檢測完整性。檢測合格的DNA稀釋至10ng/μL于4℃冰箱保存?zhèn)溆谩?/p>

1.3 微衛(wèi)星標記的PCR擴增及檢測

參照魯翠云等[16]報道的序列合成24對微衛(wèi)星引物，對F2親本、F3成活個體和F3死亡個體3組樣品進行PCR擴增。方法均參考魯翠云等[16]。

1.4 數(shù)據(jù)統(tǒng)計與分析

根據(jù)“魚類種質資源遺傳分析平臺（專利號：ZL200710144749.3）”分析軟件，將3組實驗魚的基因型數(shù)據(jù)轉換成0、1。然后根據(jù)遺傳參數(shù)計算要求轉換成PopGene（Version 1.32）軟件識別的數(shù)據(jù)格式。采用PopGene軟件統(tǒng)計3個群體在所有微衛(wèi)星標記的等位基因頻率、等位基因數(shù)（N）、有效等位基因數(shù)（Ne）、觀察雜合度（Ho）和期望雜合度（He）等遺傳參數(shù)；同時對每個群體進行哈德－溫伯格（Hardy-Weinberg）平衡檢驗（PHWE），用 Bonferroni算法進行校正。多態(tài)信息含量（PIC）由Botstein等[17]報道的公式計算。

利用PopGen32軟件計算親本及子代間的Nei’s遺傳距離。使用STRUCTURE（Version 2.3.3）聚類軟件推測群體的遺傳結構，評估親本和死亡子代之間的遺傳結構差異。

2 結果與分析

2.1 微衛(wèi)星標記的擴增效果

24對微衛(wèi)星引物，經篩選有8對引物無多態(tài)，3對引物存在缺失。將13對引物擴增穩(wěn)定、條帶清晰、呈多態(tài)的擴增結果用于后續(xù)的遺傳參數(shù)分析。

13個微衛(wèi)星標記在87個大鱗鲃個體中均獲得穩(wěn)定、清晰的DNA條帶，在個體間表現(xiàn)出不同程度的多態(tài)性，等位基因數(shù)在2～7個之間，共檢測到57個等位基因（圖1）。

2.2 親本與子代的遺傳參數(shù)分析

親本在13個微衛(wèi)星座位的Ne為1.54（B68）～5.23（B51），平均 3.22；Ho為 0.45（B1）～0.90（B51），平均 0.60；He為 0.36（B26）～0.83（B51），平均 0.65；PIC 為 0.29（B68）～0.78（B51），平均 0.52。死亡個體的 Ne為 1.35（B68）～5.11（B51），平均 3.14；Ho為0.41（B1）～0.88（B51），平均 0.58；He為 0.326（B26）～0.82（B51），平均 0.65；PIC 為 0.22（B68）～0.78（B51），平均 0.51。成活個體的 Ne為 1.44（B1）～5.24（B51），平均 3.09；Ho為 0.21（B37）～0.88（B27），平均0.64；He為 0.31（B1）～0.82（B51），平均 0.64；PIC 為0.28（B1）～0.78（B51），平均 0.50（表 1，表 2）。親本及子代在各遺傳參數(shù)間無顯著差異；但各群體內Ne顯著低于 N（P＜0.05），雖然 Ho與 He無顯著差異，但Ho整體低于 He。

各群體Hardy-Weinberg平衡χ2檢驗的無偏估測（PHWE）顯示，親本在 3 個座位（B32、B35、B43）出現(xiàn)遺傳偏離；死亡個體出現(xiàn)偏離座位最多，有7個座位（B1、B20、B26、B32、B37、B43、B58）；成活個體只在 2個座位（B37、B58）出現(xiàn)顯著偏離（P＜0.05）。

表2 親本與子代在所有微衛(wèi)星座位的群體遺傳參數(shù)統(tǒng)計Tab.2 Summary of genetic parameters of parents and offsprings across all microsatellite loci

2.3 親本與子代的親緣關系分析

親本與子代的遺傳距離和遺傳相似度如表3所示。由表3可知，親本與子代的遺傳距離非常?。?.0234～0.0466），具有很高的遺傳相似性（0.9545～0.9769），說明子代來源于親本，是親本的直系后代。其中，親本與成活個體的遺傳距離最小，相似性最高，說明大部分親本與成活個體的遺傳組成直接相關，而只有少數(shù)親本對死亡個體具有直接的遺傳貢獻。根據(jù)各群體之間的遺傳距離，利用MEGA6.06軟件，采用UPGMA法對3個群體進行聚類分析，結果親本和成活個體聚在一起，死亡個體單獨聚為1支（圖2）。

圖1 大鱗鲃部分親本及子代在微衛(wèi)星標記B26和B35的擴增電泳圖Fig.1 Electrophoresis pattern of PCR products in some parents and their offsprings of barbell at B26 and B35 loci

2.4 親本與子代的遺傳結構分析

STRUCTURE軟件對親本和子代遺傳結構分析結果顯示（圖3），當群體k＝2時，3個群體清晰地劃分成兩大結構。有14個親本與25個成活個體及6 個死亡個體（29、36、37、39、48、50）劃分到第一類群（劃分概率為51.8%～97.7%）；有 6個親本（7、8、9、10、15、20） 與 28 個死亡個體及 8 個成活個體（61、62、64、65、70、71、74、76）劃分到第二類群（劃分概率為52.9%～98.4%）。

表1 親本與子代在每個微衛(wèi)星座位的遺傳參數(shù)統(tǒng)計Tab.1 Summary of genetic parameters of parents and offsprings at each microsatellite locus

表3 親本與子代的遺傳距離及遺傳相似性指數(shù)Tab.3 Genetic distance and similarities between parents and their offsprings

圖2 親本與子代親緣關系的NJ聚類圖Fig.2 NJ dendrogram of genetic relationship between parents and their offsprings

3 討論

3.1 人工養(yǎng)殖大鱗鲃F2及F3群體的遺傳多樣性

多種遺傳參數(shù)均可衡量群體的遺傳多樣性，如雜合度、多態(tài)位點比例、等位基因數(shù)等[18]。其中，等位基因數(shù)和雜合度是微衛(wèi)星分子標記數(shù)據(jù)分析中常用的兩種多樣性指標。這兩種指標表明：本研究采用13個多態(tài)性微衛(wèi)星標記在87個大鱗鲃個體中分別檢測到2～7個等位基因，平均等位基因數(shù)（4.27～4.53）極顯著高于平均有效等位基因數(shù)（3.09～3.22）（P＜0.01），表明人工養(yǎng)殖大鱗鲃 F2及F3等位基因分布不均勻；雖然觀察雜合度（0.58～0.64）整體低于期望雜合度（0.64～0.65），但依據(jù)Botstein等[17]關于多態(tài)信息含量的劃分標準，大鱗鲃F2及F3群體處于中高度多態(tài)性水平（PIC≥0.50）。然而，遺傳多樣性水平的評估，還需考慮所用標記數(shù)量及群體大小。魯翠云等[19]和李歐等[20]通過設定標記和樣本容量梯度，分別探討?zhàn)B殖群體和野生群體遺傳多樣性評估所需的最適樣本量與標記量，認為標記數(shù)量在20～25，樣本量在40～50之間，各遺傳參數(shù)值均趨于穩(wěn)定，能夠客觀地評價一個物種的遺傳多樣性水平。本研究雖然13個標記檢測結果顯示大鱗鲃群體的遺傳多樣性水平較高，但是與魯翠云等[16]報道的多態(tài)性標記數(shù)量相比，多態(tài)位點比例明顯下降（24個標記中有8個位點無多態(tài)），說明大鱗鲃F2及F3純合度增加，遺傳多樣性實際呈下降趨勢。

圖3 每個個體分配到群體k的概率Fig.3 Probability of each individual assigned to group k

3.2 人工養(yǎng)殖大鱗鲃F3子代死亡的遺傳因素

與人工養(yǎng)殖大鱗鲃F1相比，F(xiàn)2及F3的遺傳多樣性下降。HWE平衡指數(shù)分析顯示，F(xiàn)3死亡子代有7個座位出現(xiàn)顯著偏離，占總標記數(shù)量的53.8%，而雜合子缺失及非隨機交配群體是造成偏離HWE平衡的主要原因[21]。其次，親緣關系及遺傳結構分析顯示，少數(shù)親本自交（圖2）導致大多數(shù)死亡個體與親本及成活個體的遺傳相似度降低，單獨聚為1支（圖1）。有研究學者認為，群體的遺傳多樣性每喪失10%，就會對其繁育力、存活率、生長等重要性狀產生很大的負面影響[22,23]。綜合這兩點研究結果，認為，由于幾個親本的近親交配，子代基因組純合度增加，遺傳多樣性和抗逆能力下降，這可能是F3出現(xiàn)不明緣由發(fā)病死亡的主要原因之一。

3.3 人工養(yǎng)殖大鱗鲃的一些建議

大鱗鲃是引進種，經過池塘馴化培育至性成熟的親本有限，而實際參與繁殖F1的親魚數(shù)量也很有限[12]。雖然F1的遺傳多樣性較豐富，但由于奠基者效應，F(xiàn)1的遺傳組成非常相近，很容易導致近交衰退。本研究數(shù)據(jù)表明：由于累代繁殖大鱗鲃F2及F3的多態(tài)性位點比例下降，基因組純合度增加顯著，遺傳多樣性呈下降趨勢。為避免這種情況加劇，建議對現(xiàn)有繁殖親本及后備親本進行提純復壯，如利用微衛(wèi)星標記分析繁育親本輔助遺傳多樣性和遺傳結構，累代繁殖同樣可以維持較高的雜合度，增強養(yǎng)殖群體的遺傳可塑性[24]；建議在養(yǎng)殖生產中科學地進行繁殖生產，以規(guī)避更大的經濟損失。

［1］藺玉華,耿龍武,王信海,等.池塘馴養(yǎng)大鱗鲃某些生物學特性研究［J］.天津師范大學學報:自然科學版,2009,10（4）:72-75.

［2］CanliM,Ay? and KalayM.Levels ofheavymetals（Cd,Pb,Cu,Cr and Ni）in tissue of Cyprinus carpio,Barbus capito and Chondrostoma regium from the Seyhan River,Turkey［J］.Turkish JournalofZoology,1998,22（2）:149-157.

［3］耿龍武,徐偉,藺玉華,等.大鱗鲃人工繁育技術初報［J］.吉林農業(yè)大學學報,2010,32（2）:218-220.

［4］徐偉,耿龍武,李池陶,等.大鱗鲃的人工繁殖、胚胎發(fā)育和耐鹽堿測定［J］.水產學報,2011,35（2）:255-260.

［5］徐敏,楊建,姜海峰,等.NaCl鹽度和NaHCO3堿度對鯉、鯽和大鱗鲃的精子活力及其受精率的影響［J］.中國水產科學,2014,21（4）:720-728.

［6］黨云飛,徐偉,耿龍武,等.NaCl鹽度和NaHCO3堿度對大鱗鲃幼魚生長及鰓組織特征的影響［J］.中國水產科學,2013,20（3）:577-584.

［7］耿龍武,徐偉,李池陶,等.鹽堿對大鱗鲃血清滲透壓、離子含量及鰓絲Na+/K+-ATP酶活力的影響［J］.中國水產科學,2011,18（2）:458-465.

［8］耿龍武,徐偉,李池陶,等.大鱗鲃耗氧率和窒息點的測定［J］.上海海洋大學學報,2012,21（3）:363-367.

［9］徐偉,耿龍武,李池陶,等.北方寒地池養(yǎng)大鱗鲃生長和越冬成活［J］.淡水漁業(yè),2012,42（1）:68-71.

［10］許紅,王常安,徐奇友,等.大鱗鲃幼魚氨基酸需要量［J］.華中農業(yè)大學學報,2013,32（6）:126-131.

［11］耿龍武,徐偉,金國雄,等.4種常用水產藥物對大鱗鲃的急性毒性試驗［J］.江西農業(yè)大學學報,2011,33（3）:568-572.

［12］徐偉,耿龍武,苗建生,等.耐鹽堿魚類大鱗鲃的研究現(xiàn)狀［J］.天津農學院學報,2012,19（3）:62-64.

［13］丁辰龍,藺玉華,王宣朋,等.大鱗鲃的人工繁殖及體色變異分析［J］.天津師范大學學報:自然科學版,2016,36（5）:60-64.

［14］孫效文.魚類分子育種學［M］.北京:海洋出版社,2010.

［15］Geng L W,Lu C Y,Tong G X,et al.Development and characterization of twenty microsatellite markers for the endangered fish Luciobarbus capito［J］.Conservation Genetics Resources,2012,4（4）:865-867.

［16］魯翠云,耿龍武,李超,等.微衛(wèi)星標記分析大鱗鲃養(yǎng)殖群體的遺傳結構及生長性狀［J］.水產學報,2013,37（8）:1121-1128.

［17］Botstein D,White R L,Skolnich M,et al.Construction ofa genetic linkage map in man using restriction fragment length polymorphisms［J］.American Journal of Human.Genetics,1980,32（3）:314-331.

［18］Richards C and Leberg P L.Temporal changes in allele frequencies and a population’s history of severe bottlenecks［J］.Conservation Biology,1996,10（3）:832-839.

［19］魯翠云,金萬昆,孫效文,等.樣本容量對養(yǎng)殖群體內主要遺傳結構分析參數(shù)的影響［J］.水產學報,2008,32（5）:674-683.

［20］李鷗,趙瑩瑩,郭娜,等.草魚種群SSR分析中樣本量及標記數(shù)量對遺傳多度的影響［J］.動物學研究,2009,30（2）:121-130.

［21］Reed D H and Frankham R.Correlation between fitness and genetic diversity［J］.Conservation Biology,2003,17（1）:230-237.

［22］ChangY M,LiangLQ,Ma H T,et al.Microsatellite analysis of genetic diversity and population structure of Chinese mitten crab（Eriocheir sinensis）［J］.Journal of Genetics and Genomics,2008,35（3）:171-176.

［23］Allendorf F W and Ryman N.Genetic management of hatchery stock.Population genetics and fishery management［M］.Seattle：University of Washington Press,1987:141-159.

［24］常玉梅,王文文,徐萬土,等.人工繁育大黃魚（Pseudosciaena crocea）群體 F2及 F3遺傳差異分析［J］.海洋與湖沼,2009,40（4）:414-422.

Microsatellite Marker Analysis of Disease Mortality of F3Individuals in Cultured Barbel Barbus capito

CHANG Yu-mei1,LU Cui-yun1,MIAO Jian-fa2,LI Chong-wen2,LIANG Li-qun1
(1.Heilongjiang River Fisheries Research Institute,Chinese Academy of Fishery Sciences,Harbin 150070,China;2.Tianjin Tianxiang Aquatic Company Limited,Tianjin 301500,China)

The genotypes were analyzed in F2(parents)and F3(offspring)of cultured barbel Barbus capito sampled from Tianjin Tianxiang Aquatic Company Limited by 24 reported microsatellite markers.The results showed that there were eight markers presenting no polymorphisms,and that genetic diversity,genetic relationship and genetic structure between F2and F3individuals were ultimately calculated by the genotyping data of 13 markers.The genetic parameters revealed that 2～7 alleles were detected in a total of 87 individuals at each microsatellite locus;there was higher number of average alleles (4.27～4.53)than that of the average effective alleles(3.09～3.22);the values of observed heterozogosity (0.58～0.64)lower than the values of expected heterozogosity(0.64～0.65)across all individuals and microsatellite loci,as well as the decrease of the ratio of polymorphic loci(8/24),indicating the trends of increase in homozogosity and decrease in genetic diversity in the genomes of F2and F3.Phylogenetic analysis showed that the genetic similarities of parents were more close to the surviving individuals,and clustered together,whereas,the dead individuals clustered singly.Genetic structure analysis revealed that most parents and most of the surviving individuals were classified as the first group with probability of 52.9%～98.4%,while a few parents and most of the dead individuals were classified into second groups with probability of 51.8%～97.7%,which consistent with the results of phylogenetic analysis,demonstrating that a few parents are the direct contributors to the genetic composition of the dead individuals.The findings indicate that the genetic purity of offspring increased and the genetic diversity decreased due to the inbreeding of several parents,which may be one of the main factors leading to disease m ortalityin F3individuals.

cultured barbel Barbus capito;m icrosatellite marker;geneticdiversity;genetic structure;inbreeding depression

S917

A

1005-3832（2017）05-0001-06

2017-04-17

中國水產科學研究院中央級公益性科研院所基本科研業(yè)務費專項（2014A11JC04）；科技基礎性工作專項（2013FY110700）.

常玉梅（1978-），女，博士，從事魚類抗逆育種研究.E-mail:changyumei@hrfri.ac.cn

梁利群，研究員.E-mail:liangliqun@hrfri.ac.cn