閆春梅, 杜曉燕, 鄭 偉，劉慧吉，肖志國，張家松，楊曉波

( 1.吉林省水產(chǎn)科學研究院，吉林 長春 130033；2.和龍市青龍漁業(yè)有限公司，吉林 延吉 133000 )

中國東北部山區(qū)細鱗魚4個群體遺傳多樣性分析

閆春梅1, 杜曉燕1, 鄭 偉1，劉慧吉1，肖志國1，張家松1，楊曉波2

( 1.吉林省水產(chǎn)科學研究院，吉林 長春 130033；
2.和龍市青龍漁業(yè)有限公司，吉林 延吉 133000 )

采用微衛(wèi)星法，利用11對微衛(wèi)星引物對黑龍江省松花江流域、吉林省圖們江支流紅旗河流域、吉林省白山市鴨綠江流域及吉林省圖們江支流密江流域隨機采集的94尾細鱗魚樣本進行遺傳多樣性分析，探討4個不同流域細鱗魚群體的遺傳多樣性差異及其親緣關系。結(jié)果顯示，4個群體共檢測到等位基因147個，各群體有效等位基因數(shù)為1.6541～2.4143，平均觀測雜合度為0.3849～0.5186，平均期望雜合度為0.3488～0.5003，4個群體平均多態(tài)信息含量為0.3320～0.4704。親緣關系分析顯示，紅旗河流域群體與密江流域群體遺傳距離最小(0.0954)，親緣關系最近。松花江流域群體與其他3個群體的遺傳距離較大，親緣關系較遠。群體間相似指數(shù)在0.7112以上，相似性較高。

細鱗魚；遺傳多樣性；微衛(wèi)星

細鱗魚(Brachymystaxlenok)，屬鮭形目、鮭科、細鱗魚屬，屬名貴冷水性經(jīng)濟魚類，營養(yǎng)價值和經(jīng)濟價值均較高，在我國主要分布于黑龍江、圖們江、鴨綠江、松花江、遼河支流太子河、河北省小灤河上游、秦嶺及額爾齊斯河等，目前僅在中國東北部有較集中分布。由于生態(tài)環(huán)境破壞和人為因素的雙重影響，細鱗魚自然繁殖受到影響，生殖親體出現(xiàn)退化趨勢，種群數(shù)量逐年減少，分布范圍日趨縮小，野生資源量急劇衰退，已列入中國瀕危動物紅皮書(魚類)，屬國家二級保護水生野生動物，瀕危種。近年來，對細鱗魚物種的保護逐漸受到重視，先后有水產(chǎn)行業(yè)專家學者對細鱗魚生存環(huán)境[1]及其生物學特性[2]、人工繁育[3]、遺傳多樣性[4-8]等方面展開研究。本研究在諸多研究基礎上，應用微衛(wèi)星分子標記方法，對分布于黑龍江省松花江流域、吉林省圖們江支流紅旗河流域、吉林省白山市鴨綠江流域及吉林省圖們江支流密江流域的細鱗魚群體進行了遺傳多樣性分析及其親緣關系比較研究，旨在進一步為細鱗魚種質(zhì)資源保護及遺傳改良和合理有效地開發(fā)利用提供分子生物學理論依據(jù)。

1 材料與方法

1.1 材料

試驗用細鱗魚組織樣品來源于中國東北部山區(qū)4個不同地域，分別為黑龍江省松花江流域，隨機采集11尾，命名為松花江群體；吉林省圖們江支流紅旗河流域，隨機采集30尾，命名為紅旗河群體；吉林省白山市鴨綠江流域，隨機采集23尾，命名為鴨綠江群體；吉林省圖們江支流密江流域，隨機采集30尾，命名為密江群體。所有鰭條組織樣品均保存于95%乙醇中。

基因組DNA試劑盒購自北京鼎國昌盛生物技術有限公司。引物合成委托上海生物工程有限公司完成。

1.2 基因組DNA提取

保存于95%乙醇中的組織樣品經(jīng)超純水清洗后晾干，應用基因組DNA提取試劑盒提取所有組織樣品基因組DNA。1%瓊脂糖凝膠電泳檢測，并用紫外分光光度計法測定其含量和純度，采用滅菌超純水稀釋至100 ng/μL，-20 ℃保存?zhèn)溆谩?/p>

1.3 引物合成及PCR擴增

NCBI查找相關細鱗魚微衛(wèi)星引物序列27對，合成普通引物。對普通引物進行可用性篩選，摸索引物條件和體系。進行PCR擴增后產(chǎn)物上2%瓊脂糖凝膠電泳檢測，確認符合目的長度的引物11對，引物序列見表1。

對上一步得到有條帶的引物進行熒光引物合成，在其中1條引物的5′端加FAM/HEX/Tamra標記。對所有DNA樣品進行PCR擴增檢測，PCR反應體系為10 μL[Primers (10 mmol/L)0.4 μL，10×PCR Buffer 1 μL，dNTP 0.2 μL，50 nmol/L MgCl20.4 μL，Taq 酶 0.04 μL，DNA模板 1 μL，ddH2O 6.96 μL]。94 ℃預變性2 min；94 ℃變性30 s,退火30 s，72 ℃延伸30 s，5個循環(huán)；94 ℃ 30 s，50 ℃ 30 s，72 ℃ 30 s，30個循環(huán)；72 ℃ 5 min。PCR產(chǎn)物上3730 xl基因分析儀檢測。

表1 用于微衛(wèi)星分析引物序列

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

應用Popgen 32軟件進行數(shù)據(jù)處理，分別完成等位基因數(shù)、有效等位基因數(shù)、觀測雜合度、期望雜合度、群體間遺傳距離和相似指數(shù)、聚類分析、哈迪溫伯格平衡檢驗(P值)以及平均多態(tài)信息含量計算。

2 結(jié)果與分析

2.1 群體遺傳多樣性分析

篩選出的11個位點在4個群體中的等位基因數(shù)、有效等位基因數(shù)、觀測雜合度、期望雜合度、平均多肽信息含量和哈迪溫伯格平衡數(shù)據(jù)見表2。由表2可見，4個群體94個樣品，應用11對微衛(wèi)星引物進行DNA檢測，共檢測到等位基因147個，其中松花江群體24個、紅旗河群體50個、鴨綠江群體43個、密江群體30個。由于不同群體檢測的等位基因不一致，共獲得51個不同等位基因。其中每個位點等位基因數(shù)2～7個。OMM1008和OMM1077各檢測到7個等位基因，OMM1044和OMM1105在松花江群體和密江群體兩個群體分別只有1個等位基因，在其他兩個群體中為2～5個等位基因。

各群體有效等位基因數(shù)為1.6541～2.4143，平均觀測雜合度為0.3849～0.5186，平均期望雜合度為0.3488～0.5003。其中紅旗河群體的上述指標數(shù)值均為最高。

利用平均多態(tài)信息含量計算軟件計算4個群體平均多態(tài)信息含量為0.3320～0.4704，紅旗河群體的多肽信息含量為最高(0.4704)。除OMM1044和OMM1105為低多態(tài)性位點(平均多肽信息含量<0.25)，OMM1077和TL35為高度多態(tài)(平均多肽信息含量>0.5)，其他7個位點均為中度多態(tài)(0.25<平均多肽信息含量<0.5)。

11對微衛(wèi)星引物在4個群體中檢測的44個位點中，OMM1044和OMM1105在松花江群體和密江群體兩個群體中為單態(tài)，對其他40個多態(tài)位點進行分析，29個位點符合哈迪溫伯格平衡(P>0.05)，其他11個位點均有偏離哈迪溫伯格平衡(P<0.05)的群體，偏離哈迪溫伯格平衡的現(xiàn)象在4個群體均存在。本試驗11對引物中有4對引物(OMM1007、OMM1008、OMM1039、OMM1077)在4個群體中均符合平衡。

2.2 群體間遺傳結(jié)構(gòu)

11對引物對4個群體進行PCR擴增，產(chǎn)物條帶經(jīng)聚類分析顯示，紅旗河群體與密江群體遺傳距離最小(0.0954)，親緣關系最近。松花江群體與其他3個群體的遺傳距離較大，親緣關系較遠(表3，圖1)。

通過基因流和群體間遺傳分化系數(shù)可見(表3)，4個群體間的遺傳分化系數(shù)值為0.0580～0.1751，其中松花江群體—鴨綠江群體間的遺傳分化系數(shù)值最高(0.1751)。松花江群體—密江群體和松花江群體—紅旗河群體的群體間遺傳分化系數(shù)值僅次于松花江群體—鴨綠江群體。紅旗河群體—鴨綠江群體和紅旗河群體—密江群體的群體間遺傳分化系數(shù)值最低。4個群體間的基因流值為1.1930～4.0581。其中紅旗河群體與鴨綠江群體和密江群體兩群體間的基因流值較高，分別為4.0581和3.6423。

表2 4個細鱗魚群體的11個微衛(wèi)星位點的遺傳多樣性分析

(續(xù)表2)

位點群體等位基因數(shù)有效等位基因數(shù)觀測雜合度期望雜合度多態(tài)信息含量哈迪溫伯格平衡檢驗OMM1044松花江群體11.00000.00000.00000.00000.000000**紅旗河群體41.31960.06670.24630.23210.000536**鴨綠江群體42.55560.34780.62220.54930.047738*密江群體11.00000.00000.00000.00000.000000**OMM1077松花江群體32.14160.81820.55840.43160.073314紅旗河群體64.11900.66670.77010.71680.080230鴨綠江群體75.53140.77270.83830.79530.225080密江群體42.15050.70000.54410.45840.508318OMM1105松花江群體11.00000.00000.00000.00000.000000**紅旗河群體51.32060.26670.24690.23350.999754鴨綠江群體21.18880.17390.16230.14620.592661密江群體11.00000.00000.00000.00000.000000**TL24松花江群體21.11720.11110.11110.09951.000000紅旗河群體41.93940.54170.49470.40210.932953鴨綠江群體53.56230.47830.73530.67510.045465*密江群體21.18030.16670.15540.14160.546181TL35松花江群體32.39600.45450.61040.50750.467053紅旗河群體52.82570.53330.65710.60430.359580鴨綠江群體42.76960.60870.65310.56800.043662*密江群體21.99120.60000.50620.37390.299619平均值松花江群體2.18181.66740.38490.34920.3320紅旗河群體4.54552.41430.51860.50030.4704鴨綠江群體3.90912.37820.42600.49420.4378密江群體2.72731.65410.42100.34880.3500

注:*:P<0.05; **:P<0.01.

表3 4個群體的遺傳分化系數(shù)及基因流

注：對角線以下為群體間遺傳分化系數(shù)，對角線以上為基因流.

表4 4個群體的遺傳距離和相似性指數(shù)

注：對角線以下為遺傳距離，對角線以上為遺傳相似指數(shù).

圖1 4個群體的UPGMA聚類圖

3 討 論

微衛(wèi)星法是檢測共顯性遺傳特征的分辨率很高的較新遺傳學研究方法，目前已成為種群遺傳結(jié)構(gòu)與進化分析的有效分子標記。本試驗利用微衛(wèi)星法分析了4個不同流域細鱗魚群體的遺傳多樣性和遺傳結(jié)構(gòu)。

3.1 遺傳多樣性

遺傳多樣性是衡量一個種群種質(zhì)資源質(zhì)量的重要指標。遺傳多樣性越高，進化潛力越大[9]。本研究主要從等位基因數(shù)、雜合度、多態(tài)信息含量、哈迪溫伯格平衡等幾個方面對4個不同流域細鱗魚群體的遺傳多樣性進行分析。平均雜合度是被檢測位點上群體的雜合子頻率，用于檢測群體內(nèi)遺傳變異，是群體雜合程度的度量單位。雜合度高低與遺傳多樣性豐度和群體對環(huán)境的適應能力成正相關[10]。本試驗中，4個群體平均有效等位基因數(shù)為1.6541～2.4143，平均觀測雜合度0.3849～0.5186，平均期望雜合度0.3488～0.5003。表明4個群體遺傳多樣性均較豐富，尤其是紅旗河群體，群體平均觀測雜合度0.5186，平均期望雜合度0.5003，遺傳多樣性高于其他3個群體，對環(huán)境適應能力較強，具有較豐富的育種和遺傳改良潛力，應采取相應措施予以保護。

除雜合度以外，還可以通過平均多肽信息含量來檢測多態(tài)性，從而來度量遺傳多樣性和基因豐富程度。Botstein等[11]認為，當平均多肽信息含量<0.25時為低度多態(tài)位點，0.25<平均多肽信息含量<0.5時為中度多態(tài)位點，平均多肽信息含量>0.5時為高度多態(tài)位點。本研究中11個微衛(wèi)星位點中包括2個低度多態(tài)性位點、7個中度多態(tài)性位點和2個高度多態(tài)性位點，位點多態(tài)性較豐富。對4個群體分析結(jié)果表明，4個群體的平均多態(tài)信息含量為0.3192～0.4704，介于0.25～0.5之間，證明群體間遺傳多樣性中等，4個群體均屬于中度多態(tài)性。

哈迪溫伯格平衡是指在一個群體無限大，且又具備隨機交配、沒有突變、沒有選擇、沒有遺傳漂變的情況下，群體內(nèi)一個位點上的基因型頻率和基因頻率將保持世代不變，處于遺傳平衡狀態(tài)[12]。本研究中11個微衛(wèi)星引物在4個群體中均觀測到顯著偏離平衡情況，40個多態(tài)位點中11個位點偏離平衡(P<0.05)。這種現(xiàn)象可能與樣本量及無效等位基因的存在有關[13]。

3.2 群體間遺傳結(jié)構(gòu)

本研究的目的是檢驗4個細鱗魚群體的遺傳關系，因此，遺傳距離和遺傳相似度是重要的衡量指標。在群體的進化研究中，可通過計算遺傳距離來鑒定種間遺傳結(jié)構(gòu)及分化關系。微衛(wèi)星研究方法認為遺傳距離越遠，分化時間越長，雜種優(yōu)勢越大，即遺傳距離與分化關系呈正相關。1982年Thorp綜合一系列研究認為，遺傳距離大于0.05同時小于0.3為同種不同群體；遺傳距離大于0.3同時小于0.9為同屬不同種；遺傳距離大于0.9為不同屬[14]。本試驗4個群體間遺傳距離均大于0.0954，認為4個群體間發(fā)生了遺傳分化。其中松花江群體與鴨綠江群體遺傳距離為0.3408，稍高于0.3，認為種群間開始發(fā)生種的分化，其他3個群體遺傳距離均小于0.3，認為種群間已發(fā)生群體分化，為同種不同群體。綜合來看，4個群體間雖有分化，但分化程度不高。聚類分析結(jié)果也證實紅旗河群體與密江群體首先聚為一支，松花江群體單獨分為一支。出現(xiàn)這一結(jié)果的原因可能是紅旗河群體與密江群體均屬吉林省圖們江支流，地理位置相對較近，因此聚為一支；而鴨綠江群體屬吉林省鴨綠江流域，與紅旗河群體與密江群體同屬吉林省但流域不同，因此遺傳距離稍遠；松花江群體則屬于黑龍江省松花江流域，地理位置最遠，所以單獨分為一支。

此外，基因流和群體間遺傳分化系數(shù)也可用來分析群體間遺傳分化程度?；蛄?1說明群體分化程度較大，1<基因流<4則是分化程度一般，基因流>4說明群體間該位點的遺傳分化程度較小[15]。本研究中4個群體之間的基因流值均大于1.1930，表明群體間發(fā)生了一般程度的基因交流，而紅旗河群體與鴨綠江群體和密江群體兩群體間的遺傳分化程度更小，基因流值分別為4.0581和3.6423。遺傳分化系數(shù)是反映種群間的遺傳分化程度的重要指標。當 0<遺傳分化系數(shù)≤0.05時, 表明群體遺傳分化較弱;0.05<遺傳分化系數(shù)≤0.15時, 表明群體遺傳分化中等; 0.15<遺傳分化系數(shù)≤0.25時,表明群體遺傳分化較大; 遺傳分化系數(shù)>0.25時,表明群體遺傳分化很大[16]。本研究中松花江群體與另外3個群體種群間的遺傳分化系數(shù)分別為0.1402、0.1751、0.1733，遺傳分化程度較大，紅旗河群體與鴨綠江群體和密江群體兩群體間的種群間遺傳分化系數(shù)分別為0.0580、0.0642，鴨綠江群體和密江群體的種群間遺傳分化系數(shù)為0.1391。證明其遺傳分化程度中等。此結(jié)論與群體間聚類分析結(jié)果基本一致。本研究分析的4個中國東北部山區(qū)細鱗魚群體遺傳多樣性可見，4個群體遺傳多樣性均較高，特別是紅旗河群體，具有一定的遺傳潛力，可作為選育群體。建議人工繁殖時選取遺傳距離較遠的個體作為親本，同時對同一繁殖群體可進行不同地理區(qū)域放流，以增加群體遺傳多樣性，保證種質(zhì)資源永續(xù)利用。

[1] 高燕.鴨綠江上游支流細鱗魚(Brachymystaxlenok)的生態(tài)環(huán)境調(diào)查及人工增殖研究[D].長春:東北師范大學, 2004.

[2] 白慶利, 牟振波, 于洪賢, 等. 野生細鱗魚的生物學特性及繁殖力[J].水產(chǎn)學雜志, 2007，20(1):69-73.

[3] 牟振波, 劉洋, 徐革鋒, 等.細鱗魚攝食和生長最適水溫的研究[J].水產(chǎn)學雜志,2011，24(4):6-8，24.

[4] 王荻, 徐革鋒, 劉洋, 等.細鱗魚三個野生種群的遺傳多樣性AFLP分析[J].農(nóng)業(yè)生物技術學報,2009，17(4):638-644.

[5] 徐革鋒, 牟振波, 薛淑群, 等.不同流域細鱗魚染色體遺傳多態(tài)性分析[J].水生生物學報,2009，33(5):975-979。

[6] 原居林, 朱俊杰, 王高學.秦嶺細鱗鮭黑河種群和湑水河種群的遺傳多樣性分析[J]. 南方水產(chǎn), 2009，5(3):63-66.

[7] 張艷萍, 王太, 杜巖巖, 等.秦嶺細鱗鮭人工繁育群體與野生群體遺傳變異分析[J].水生生物學報, 2014, 38(5):828-833.

[8] 馬波, 姜作發(fā).烏蘇里江2種細鱗鮭種群遺傳多樣性及親緣關系的微衛(wèi)星分析[J]. 中國水產(chǎn)科學, 2007, 14(1):39-45.

[9] 程雪艷, 滕爽爽, 肖國強, 等.基于微衛(wèi)星的泥蚶5個地理群體遺傳多樣性分析[J]. 海洋科學, 2015, 39(11):116-125.

[10] 鄒杰, 馬愛軍, 趙艷飛, 等.利用微衛(wèi)星標記分析3個暗紋東方鲀(Takifuguobscurus)養(yǎng)殖群體的遺傳多樣性[J].漁業(yè)科學進展,2015,36(2):55-62.

[11] Botstein D, White R L, Skolnick M,et al. Construction of a genetic linkage map in man using restriction fragment length polymorphisms[J]. Am J Hum Genet, 1980, 32(3):314-331.

[12] 黃代新, 楊慶恩.卡方檢驗和精確檢驗在HWE檢驗中的應用[J]. 法醫(yī)學雜志, 2004, 20(2):116-119.

[13] 王鶴, 王田田, 胡麗萍, 等.仿刺參中國群體與韓國群體雜交子代微衛(wèi)星標記分析[J].水產(chǎn)科學,2016, 35 (1):60-66.

[14] Shaklee J B, Tamaru C S, Waples R S. Speciation and evolution of marine fishes studied by the electrophretic analysis of proteins[J]. Pac Sci, 1982, 36(2):141-157.

[15] 劉博, 匡友誼, 佟廣香, 等.微衛(wèi)星分析9個哲羅魚野生群體的遺傳多樣性[J].動物學研究,2011, 32(6):597-604.

[16] Weight S. Evolution and the genetics of population, Volume 4: variability within and among natural population Chicago[M]. Chicago:University of Chicago Press, 1987.

GeneticDiversityofFourPopuationsofLenokBrachymystaxlenokinMountainAreaofNortheastChina

YAN Chunmei1, DU Xiaoyan1, ZHENG Wei1, LIU Huiji1, XIAO Zhiguo1, ZHANG Jiasong1, YANG Xiaobo2

( 1.Jilin Fisheries Research Institute, Changchun 130033，China; 2.Helong Qinglong Fisheries Limited Company, Yanji 133000，China )

The genetic diversities and relationship were investigated in ninety-four individuals of lenok(Brachymystoxlenok) collected from Songhua river basin in heilongjiang province, Jilin Tumen river tributary of Red Flag river basin, Jilin Baishan Yalu river basin and Jilin Tumen river tributary of Mijiang river basin by 11 pairs of microsate llite primers and by SSR. The results showed that 147 alleles were detected in the four populations, with effective allele number (Ne) of 1.6541-2.4143, average observed heterozygosity of 0.3849-0.5186, and the average expected heterozygosity of 0.3488-0.5003. Four populations of lenok had the average polymorphism information content (PIC) of 0.3320-0.4704. The phylogenetic relationship analysis indicated that the Red Flag river group and the Mijiang river group had the minimal genetic distance(0.0954), with the close tgenetic relationship. The Songhuajiang group showed larger genetic distance than the other 3 groups did, with similarity index of over 0.7112 between groups.

Brachymystaxlenok; genetic diversity; microsatellite

10.16378/j.cnki.1003-1111.2017.06.015

S917

A

1003-1111(2017)06-0778-06

2016-10-24；

2016-12-29.

吉林省省級公益性科研院所行業(yè)科研項目；吉林省現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術體系建設項目(2015-2016).

閆春梅(1981-)，女，工程師，碩士；研究方向：預防獸醫(yī)學、水產(chǎn)養(yǎng)殖. E-mail：yanchunmei199950@126.com.通訊作者：杜曉燕(1961-)，女，研究員；研究方向：水產(chǎn)生物. E-mail: dxypa@tom.com.