趙子明，葉建生，劉美劍，沈冬冬，黃世熠

(江蘇農(nóng)牧科技職業(yè)學院，泰州 225300)

?

五個湖泊中華刺鰍(Mastacembelus sinensis)種群線粒體細胞色素b遺傳多樣性分析*

趙子明，葉建生，劉美劍，沈冬冬，黃世熠

(江蘇農(nóng)牧科技職業(yè)學院，泰州 225300)

采用特異性引物對中華刺鰍基因組DNA進行擴增和測序，獲得了江蘇省5個湖泊的中華刺鰍種群線粒體cyt b基因767bp的序列. 分析表明，cyt b基因中A+T含量略高于G+C含量，144個中華刺鰍共獲得cyt b基因單倍型28個，變異位點38個，其中單倍型Hap1分布最為廣泛，可能為祖先型單倍型. 平均單倍型多樣性和平均核苷酸多樣性指數(shù)分別為0.58751和0.00214，遺傳多樣性較低. 群體間遺傳分化系數(shù)為 0.0356，絕大多數(shù)的遺傳變異都來自群體內(nèi)部，而群體間的遺傳分化極小. 構(gòu)建的系統(tǒng)進化樹未出現(xiàn)明顯的以單個湖泊群體的聚簇，5個湖泊的種群間基因交流極其強烈. 研究結(jié)果表明，當前5個湖泊的中華刺鰍遺傳多樣性水平比較低，亟待加強保護.

中華刺鰍；線粒體DNA；cyt b；遺傳多樣性；淡水湖

中華刺鰍(Mastacembelus sinensis)屬于鱸形目(Perciformes)、刺鰍科(Mastacembelidae)、刺楸屬(Mastacembelus)，俗稱石鏟、刀鰍、鋼鰍等，是廣泛分布我國遼河、黃河、長江、錢塘江、珠江等水系的中小型淡水魚類[1]. 改革開放以來，由于人為活動和江河環(huán)境污染，特別是農(nóng)作物秸稈向水體中拋置，導致水體污染嚴重，中華刺鰍野生資受到嚴重破壞. 目前，中華刺鰍資源量不斷下降.

cyt b基因作為常用線粒體DNA分子標記，具有較強的中間解析和鑒別能力，廣泛應(yīng)用于動物的群體遺傳結(jié)構(gòu)、物種鑒別和系統(tǒng)進化的研究中[2-5]. 近年來，cyt b基因也越來越多地應(yīng)用于魚類遺傳多樣性研究. 趙凱等[6]分析了黃河上游和柴達木盆地托索湖的裸裂尻魚cyt b基因序列變異，初步認為黃河和柴達木水系托索湖的裸裂尻魚未顯著分化，兩種群表現(xiàn)出較低水平的遺傳多樣性. 閻雪嵐等[7]采用cyt b序列分析了采自新疆克孜河3個群體的塔里木裂腹魚及1個斑重唇魚群體的種群遺傳多樣性和遺傳結(jié)構(gòu)，結(jié)果顯示克孜河的塔里木裂腹魚和斑重唇魚均處于很低的遺傳多樣性水平，物種維持力較弱. 趙峰等[8]采用cyt b基因片段為遺傳標記，對黃海南部和東海的6個銀鯧群體遺傳結(jié)構(gòu)進行了分析，結(jié)果表明6個群體均呈現(xiàn)出高單倍型多樣性和低核苷酸多樣性的特點，群體間具有高度的基因交流，是一個隨機交配群體. 田輝伍等[9]則利用cyt b基因研究了長江上游的長江干流、岷江、赤水河和嘉陵江的8個群體長薄鰍的遺傳結(jié)構(gòu)，表明長薄鰍種群內(nèi)的變異大于種群間的變異，種群基因交流十分頻繁. 目前，有關(guān)刺鰍的研究很少，僅有的幾篇報道主要集中于基礎(chǔ)生物學[10-15]、病害研究[16]以及分類地位討論[17-22]. 尤其是其分類地位一直存在爭議，存在刺鰍(Mastacembelus aculeatus)和中華光蓋刺鰍(Pararhynchobdlla sinensis)等多種學名，倪勇與伍漢霖先生[1]根據(jù)其形態(tài)特征比對，最終確定了其分類地位，學名為中華刺鰍(Mastacembelus sinensis). 本研究對采自江蘇省5個湖泊中華刺鰍種群cyt b片段進行了測序與比較分析，研究5個湖泊中華刺鰍群體間的遺傳多樣性和遺傳結(jié)構(gòu)，以深入了解和掌握5個湖泊中華刺鰍種群當前的種質(zhì)資源現(xiàn)狀，為其今后的保護與開發(fā)提供基礎(chǔ)分子生物學依據(jù).

1 材料和方法

1.1 研究材料

實驗用中華刺鰍于2015年5-7月采集自江蘇省蘇中和蘇北地區(qū)的5個湖泊水域，分別為江蘇淮安的白馬湖(B)、高郵湖(G)、洪澤湖(H)、駱馬湖(L)以及江蘇泰州的溱潼湖(Q)，剪取尾鰭，無水乙醇固定后，-20℃ 冷凍保存?zhèn)溆?

1.2 DNA提取

基因組DNA采用試劑盒提取(SangonEzupColumnAnimalGenomicDNAPurificationKit). 總DNA經(jīng)檢測后稀釋至100ng/μl，-20℃冷凍保存?zhèn)溆?

1.3 PCR擴增及測序

參照GenBank數(shù)據(jù)庫中大刺鰍(Mastacembelus armatus)的線粒體全序列(登入號：KJ184553)，設(shè)計引物擴增cyt b 基因，引物序列為F：5’-ATCCCTCCTGGGACTCTGTC-3’；R：5’-CGTGTGGAGGATTGGGACAA-3’，引物由上海生工生物公司合成.PCR反應(yīng)體系為50μl：模板DNA50ng；2×PCRMix25μl(包含Taq酶 2.5U，dNTPs10μmol，MgCl20.1mmol)，上、下游引物(10μmol/L)各1μl，其余體積用水補足. 反應(yīng)程序為：94℃預變性2min，94℃變性45s，62℃退火1min，72℃延伸1min，經(jīng)35個循環(huán)后再72℃延伸10min.PCR產(chǎn)物經(jīng)割膠純化后，直接送上海生工生物公司利用正反引物進行雙向測序.

1.4 數(shù)據(jù)分析

根據(jù)測序峰圖對雙向測序序列進行人工拼接，之后用BioEdit7.0.5[23]軟件比對獲得同源序列，MEGA6.0軟件[24]分析堿基的組成，并構(gòu)建系統(tǒng)進化樹.DnaSP軟件[25]統(tǒng)計變異位點類型與數(shù)目，計算單倍型種類、平均核苷酸差異數(shù)、單倍型多樣性指數(shù)和核苷酸多態(tài)性等，用Arlequin3.01 軟件[26]中的Kimura2-Parameters方法構(gòu)建5個種群間的遺傳距離，分子生物學方差分析(AMOVA)方法計算遺傳分化指數(shù)(Fst)及遺傳變異組成.

2 結(jié)果

2.1 序列分析

5個湖泊中華刺鰍種群144個個體的cyt b基因，經(jīng)拼接比對后，獲得長度為 767bp的同源序列.MEGA6軟件分析表明，T、C、A、G堿基組成平均為29.2%、30.1%、26.1%和14.6%，A+T含量為 55.3%，略高于G+C含量. 序列共發(fā)現(xiàn)多態(tài)性位點38個，占核苷酸總數(shù)的4.95%，其中包含18個單態(tài)突變位點和20個簡約信息位點，未檢測到堿基缺失. 144個個體中共檢出28個單倍型序列，所有單倍型序列已提交至GeneBank(登入號：KT448786～KT448813)，種群間共享單倍型共8個，其中單倍型Hap1 5個群體間都有分布且數(shù)量也大，可能是原始單倍型(表1). 5個種群的單倍型平均多樣性指數(shù)為0.58751，平均核苷酸差異數(shù)為1.63947，平均核苷酸多樣性指數(shù)為0.00214(表2). 5個湖泊的中華刺鰍群體中，高郵湖群體的單倍型多樣性指數(shù)最高，溱潼湖群體的核苷酸多樣性指數(shù)最高，而駱馬湖群體的單倍型多樣性指數(shù)和核苷酸多樣性指數(shù)均為最低.

表2 5個湖泊的中華刺鰍種群的單倍型及遺傳多樣性參數(shù)Tab.2 Haplotypes and genetic diversity of spiny eel populations from five lakes

湖泊樣本數(shù)單倍型數(shù)多態(tài)位點數(shù)單倍型多樣性平均核苷酸差異數(shù)核苷酸多態(tài)性白馬湖30890.466670.901150.00117高郵湖2811210.714292.359790.00308洪澤湖2810110.671961.129630.00147駱馬湖29860.642860.817730.00107溱潼湖295160.421182.788180.00364合計14428380.587511.639470.00214

2.2 群體變異和遺傳結(jié)構(gòu)

將測序獲得的144個個體的cyt b基因序列輸入Arlequin軟件，使用Kimura2-Parameters方法構(gòu)建5個湖泊種群間的相對遺傳距離和Fst(表3). 5個種群的平均遺傳距離為0.00216，高郵湖種群與溱潼湖種群間的遺傳距離最遠為0.00338，遺傳距離最近的是白馬湖和駱馬湖種群之間為0.00116. 而5個湖泊刺鰍種群內(nèi)部遺傳距離從大到小依次為：溱潼湖(0.00368)>高郵湖(0.00311)>洪澤湖(0.00148)>白馬湖(0.00118)>駱馬湖 (0.00107). 5個種群的平均Fst為0.03556，5個湖泊的中華刺鰍種群間存在一定的遺傳分化，但是種群間的Fst值均較低，種群間Fst顯著性檢驗有近一半(P>0.05)，表明5個湖泊中華刺鰍種群間存在較大的基因流，遺傳分化差異較低(表4).

表3 5個湖泊的中華刺鰍種群間K 2-P遺傳距離和遺傳分化系數(shù)Tab.3 Average K 2-P distance and fixation indice between spiny eel populations from five lakes

*表示差異顯著(P<0.05)；**表示差異極顯著(P<0.01)；對角線上方數(shù)值表示種內(nèi)分化系數(shù)；對角線下方數(shù)值表示群體間遺傳距離.

表4 5個湖泊中華刺鰍種群間遺傳差異的AMOVA分析Tab.4 Analysis of AMOVA of spiny eelpopulations from five lakes

變異來源自由度平方和變異組成百分比%種群間46.5660.029363.56種群內(nèi)139110.6570.7960996.44合計143117.2220.35778

AMOVA法估算5個群體間遺傳變異結(jié)構(gòu)和來源的結(jié)果表明，種群間的變異占總變異的3.56%，而種群內(nèi)的變異對總變異的貢獻率達到96.44%(表4)，種群內(nèi)未出現(xiàn)遺傳分化. 另外，用DnaSP軟件分析單倍型數(shù)據(jù)所得遺傳分化系數(shù)(Gst)為0.01092，基因流(Nm)為22.65；而144個個體的Gst為0.05601，Nm為4.21，也表明種群間的基因交流很頻繁.

2.3 系統(tǒng)進化分析

將144個中華刺鰍個體的cyt b基因序列輸入MEGA6軟件，對5個湖泊群體進行系統(tǒng)進化分析，得到UPGMA系統(tǒng)進化樹(圖1A). 5個中華刺鰍群體沒有明顯聚類，形成了單系類群，京杭大運河水系連接的白馬湖、駱馬湖、洪澤湖和高郵湖群體首先聚類，最后與溱潼湖聚合. 而144個個體的構(gòu)建的UPGMA系統(tǒng)進化樹則表明，5個群體的個體都沒有獨立聚類，相互散布在不同的地理群體中(圖1B).

圖1 5個湖泊中華刺鰍群體間(A)及個體間(B)的UPGMA系統(tǒng)進化樹Fig.1 Dendrogram of spiny eel populations (A) and 144 individuals (B) from five lakes by UPGMA method based genetic distance

2.4 群體歷史動態(tài)

用Fu’Fs檢驗方法檢驗5個湖泊中華刺鰍群體的中性突變，結(jié)果表明種群的整體Fs=-23.887(P<0.05)，顯著偏離中性突變，表明5個中華刺鰍群體整體上有一定的擴張趨勢，但白馬湖、駱馬湖、高郵湖和洪澤湖群體的Fs值都小于0，而溱潼湖群體的Fs值大于0. 所有個體的線粒體cyt b基因序列進行堿基歧點分布分析，個體間堿基歧點分布的觀測曲線與期望曲線大體一致，中間只有一個較低的凸起(圖2)，表明5個湖泊的中華刺鰍歷史上種群相對比較穩(wěn)定，經(jīng)歷過一個比較小的種群擴張事件.

圖2 5個湖泊的中華刺鰍群體線粒體cyt b基因錯配堿基分布Fig.2 The mismatch distribution of the Mitochondrial DNA cyt b of spiny eel populations from five lakes

3 討論

3.1 群體遺傳多樣性

本研究獲得了5個湖泊中華刺鰍群體的線粒體 cyt b 基因片段，序列堿基組分析表明，中華刺鰍線粒體cyt b 基因的A+T含量為55.3%，略高于G+C平均含量，而G含量顯著低于其他堿基，表現(xiàn)出反G偏倚，與其他動物線粒體DNA的研究結(jié)果一致[27].

遺傳多樣性是一個物種長期生存和進化的前提，物種的遺傳多樣性越高，該物種的進化潛力也就越大，對環(huán)境變化的適應(yīng)能力也就越強[5]. 本研究獲得5個湖泊144 尾中華刺鰍cyt b基因767bp的序列，共獲得基因單倍型28個，變異位點38個，單倍型多樣性和核苷酸多樣性分別為0.58751 和0.00214，說明江蘇省5個湖泊的中華刺鰍單倍型多樣性和核苷酸多樣性都很低，遠低于黃顙魚[4](0.945、0.00419)和福建近海竹莢魚[28](0.937、0.336)等經(jīng)濟魚類的cyt b基因單倍型多樣性和核苷酸多樣性，甚至低于野生瀕危物種鱇浪白魚[5](0.996、0.00434) cyt b基因的單倍型多樣性和核苷酸多樣性，僅與瀕危物種青海湖裸鯉[29](0.783、0.00205)、黃河裸體裂尻[6](0.700、0.0026)和長薄鰍[9](0.60852、0.00089) cyt b基因的單倍型多樣性和核苷酸多樣性相當，高于怒江角魚[30](0.579、0.00070) cyt b基因單倍型多樣性和核苷酸多樣性. 如此低的遺傳多樣性，可能意味著中華刺鰍種群曾經(jīng)歷過“瓶頸效應(yīng)”的打擊[31]，說明當前該物種的進化潛力比較低，對環(huán)境變化的適應(yīng)能力很弱，存在著較大的滅絕風險.

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

AMOVA分析表明，種群內(nèi)的遺傳變異顯著超過各種群間的遺傳變異，占全部遺傳變異的96.44%，這表明中華刺鰍群體間差異不顯著，群體間未出現(xiàn)遺傳分化.

Fst值可以反映群體之間的遺傳分化程度，當Fst值介于 0～0.05時表示分化較弱. 本文中5個湖泊中華刺鰍群體的Fst值為0.0356，說明5個中華刺鰍種群間存遺傳分化較弱，不存在顯著的系統(tǒng)地理分化.Shaklee等[32]根據(jù)已發(fā)表的文獻認為魚類在種群、種和屬三級水平上的遺傳距離D值分別為0.05、0.30和0.9. 中華刺鰍5個湖泊群體間遺傳距離為0.00216，說明中華刺鰍群體間的遺傳分化極低，遠沒有達到種群的分化水平. 另外，與其他魚類群體間的遺傳距離相比較(0.035～0.019)[33]，中華刺鰍群體間遺傳距離(0.00116～0.00338)可以看出，中華刺鰍種群無論是群體間還是群體內(nèi)部序列差異都極小，說明其遺傳多樣性水平極低.Wright提出種群Nm>4時，群體就是一個隨機的單位[34]. 而本研究中單倍型間的Nm高達22.65，個體間的Nm為4.21，表明5個群體間基因流動的相當強烈，也進一步說明中華刺鰍種群間未出現(xiàn)分化，應(yīng)作為一個整體進行資源保護. 系統(tǒng)進化分析表明5個湖泊的中華刺鰍群體沒有明顯聚類，形成了單系類群. 以上結(jié)果都表明中華刺鰍5個湖泊群體在cyt b水平上未出現(xiàn)群體遺傳分化.Grant等依據(jù)單倍型多樣度h 和核苷酸多樣度π推測了魚類群體的4 種進化情景[35]，本研究中5個湖泊的中華刺鰍種群遺傳多樣性低，遺傳沒有分化，接近于低h、低π情形，可能因為湖泊間地理距離近，水系連通，使得基因交流的可能性很大；同時，也可能存在建立者效應(yīng)，即5個湖泊的中華刺鰍群體都源于單一建立者；此外，江河等水體環(huán)境污染和人為活動導致棲息地減少，致使中華刺鰍群體數(shù)量急劇下降，而群體變小后，瓶頸效應(yīng)和遺傳漂變可能在某種程度上導致其遺傳多樣性的降低，損失的遺傳變異也會增多.

遺傳多樣性是生命進化的基礎(chǔ)，遺傳多樣性的缺失會威脅到種群的生存. 本研究以cyt b基因為工具，研究了5個湖泊中華刺鰍群體間的遺傳多樣性和遺傳結(jié)構(gòu)，以深入了解和掌握當前中華刺鰍種質(zhì)資源現(xiàn)狀，為其今后的保護與開發(fā)提供基礎(chǔ)分子生物學依據(jù). 5個湖泊的中華刺鰍種群處于較低的遺傳多樣性水平，且未形成遺傳分化，物種維持力較弱. 本研究結(jié)果表明當前5個湖泊的中華刺鰍資源現(xiàn)狀不容樂觀，亟待加強保護. 種群一旦遭受破壞，就很難以復，這個寶貴物種也將面臨滅絕的危險. 同時，應(yīng)盡快對其他水域中華刺鰍遺傳多樣性現(xiàn)狀進行研究，全面了解該物種種質(zhì)資源現(xiàn)狀.

[1]NiYong,WuHanlineds.FishesofJiangsuProvince.Beijing:ChinaAgriculturePress, 2006: 740-742(inChinese).[倪勇, 伍漢霖. 江蘇魚類志. 北京: 中國農(nóng)業(yè)出版社, 2006: 740-742.]

[2]BrownWM.EvolutionofanimalmitochondrialDNA.In:NeiM,KoehnRKeds.Evolutionofgenesandproteins.SunderlandMA:Sinauer, 1983: 62-88.

[3]SacconeC,PesoleG,SbisaE.ThemainregulatoryregionofmammalianmitochondrialDNA:structure-functionmodelandevolutionarypattern. J Mol Evol, 1991, 33(1): 83-91.

[4]ZhongLiqiang,LiuPengpeng,PanJianlinet al.Geneticvariationanalysisofyellowcatfish(Pelteobagrus fulvidraco)fromfivelakesinthemiddleandlowerreachesoftheYangtzeRiverbasedonMitochondrialDNAcytb. J Lake Sci, 2013, 25(2): 302-308(inChinesewithEnglishabstract).DOI::10.18307/2013.0219. [鐘立強, 劉朋朋, 潘建林等. 長江中下游5個湖泊黃顙魚種群線粒體細胞色素b基因的遺傳變異分析. 湖泊科學, 2013, 25(2): 302-308.]

[5]LiuHongyan,XiongFei,YangDonget al.MitochondrialcytochromebgenesequencediversityinwildandculturedpopulationsofAnabarilius grahami. J Huazhong Agric Univ, 2011, 30(1): 94-98(inChinesewithEnglishabstract).[劉紅艷, 熊飛, 楊東等. 基于細胞色素b基因的鱇浪白魚野生群體和養(yǎng)殖群體遺傳多樣性分析. 華中農(nóng)業(yè)大學學報, 2011, 30(1): 94-98.]

[6]ZhaoKai,YangGongshe,LiJunbinget al.PhylogeneticstructureofSchizopygopsis pylzovipopulationsfromMitochondrialcytochromebgenesequencevariations. Acta Hydrob Sinica, 2006, 30(2): 129-133(inChinesewithEnglishabstract).[趙凱, 楊公社, 李俊兵等. 黃河裸裂尻魚群體遺傳結(jié)構(gòu)和cytb 序列變異. 水生生物學報, 2006, 30(2): 129-133.]

[7]YanXuelan,YangJinquan,TangWenqiaoet al.GeneticdiversityoftwoCyprinidfishes: Schizothorax biddulphiandDiptychus maculatesinKezilesuRiver,XinjiangUygurAutonomousRegionbasedonCytbgenesequencesvariation. Chinese J Zool, 2009, 44(5): 8-13(inChinesewithEnglishabstract).[閻雪嵐, 楊金權(quán), 唐文喬等. 基于線粒體Cytb基因序列變異的克孜河塔里木裂腹魚和斑重唇魚遺傳多樣性. 動物學雜志, 2009, 44(5): 8-13.]

[8]ZhaoFeng,ZhuangPing,ZhangLongzhenet al.PopulationgeneticstructureofPampus argenteusinthesouthYellowSeaandEastChinaSeabasedonthemitochondrialCytbsequence. Acta Hydrob Sinica, 2011, 35(5): 745-752(inChinesewithEnglishabstract).[趙峰, 莊平, 章龍珍等. 基于線粒體Cytb基因的黃海南部和東海銀鯧群體遺傳結(jié)構(gòu)分析. 水生生物學報, 2011, 35(5): 745-752.]

[9]TianHuiwu,DuanXinbin,WangDengqianget al.SequencevariabilityofCytochromebandgeneticstructureofLeptobotia elongataintheupperYangtzeRiver. Freshw Fish, 2013, 43(6): 13-18(inChinesewithEnglishabstract).[田輝伍, 段辛斌, 汪登強等. 長江上游長薄鰍Cytb基因的序列變異與遺傳結(jié)構(gòu)分析. 淡水漁業(yè), 2013, 43(6): 13-18.]

[10]WuLv.CellularcompositionoftheisletsinMastacembelus. Science, 1950, 32(3): 91(inChinese).[伍律. 刺鰍胰島之細胞組成. 科學, 1950, 32(3): 91.]

[11]WuLv.Thepancreasandthepancreaticisletsofateleost, Mastacembelus aculeatus. Acta Hydrob Sinica, 1956, 2(2): 272-278(inChinesewithEnglishabstract).[伍律. 刺鰍的胰臟和胰島. 水生生物學集刊, 1956, 2(2): 272-278.]

[12]ZhouRongjia,YuQixing,ChengHanhuaet al.PCRamplificationofSRYboxgeneinMonopterus albusandMastacembelus aculeatus. Chinese Sci Bull, 1996, 41(7): 640-642(inChinesewithEnglishabstract).[周榮家, 余其興, 程漢華等.PCR擴增黃鱔和刺鰍SRY盒基因. 科學通報, 1996, 41(7): 640-642.]

[13]LiuJiangdong,HuangXiao,YuQixinget al.ThecytogeneticevidencesforconfirmingtheexistenceofsexchromosomesinMastacembelus aculeatus (Osteichthyes,Perciformes). J Wuhan Univ (Nat Sci Ed), 1999, 45(2): 185-190(inChinesewithEnglishabstract).[劉江東, 黃曉, 余其興等. 刺鰍性染色體的細胞遺傳學確定證據(jù). 武漢大學學報: 自然科學版, 1999, 45(2): 185-190.]

[14]ChenJi,ZhaoGang,ZangYatinget al.TheconstructionofXchromosomelibraryofspinyeel(Mastacembelus aculeatus). Acta Hydrob Sinica, 2009, 33(4): 571-576(inChinesewithEnglishabstract).[陳戟, 趙剛, 臧亞婷等. 刺鰍X染色體DNA文庫的構(gòu)建. 水生生物學報, 2009, 33(4): 571-576.]

[15]WuYuanan,LiangZhiqiang,LiChuanwuet al.AnalysisandevaluationofnutritionalcomponentsinthemuscleoftwokindsofMastacembelus. Acta Nutrimenta Sinica, 2010, 32(5): 499-502(inChinesewithEnglishabstract).[伍遠安, 梁志強, 李傳武等. 兩種刺鰍肌肉營養(yǎng)成分分析及評價. 營養(yǎng)學報, 2010, 32(5): 499-502.]

[16]HanJiugao.InfectionandcontrolofIchthyophthiriasisspp.onMastacembelus aculeatus. J Anhui Agri Sci, 2007, 35(16): 4826-4834(inChinesewithEnglishabstract).[韓九皋. 刺鰍感染小瓜蟲病及防治研究. 安徽農(nóng)業(yè)科學, 2007, 35(16): 4826-4834.]

[17]TraversRA.AreviewoftheMastacembeloidei,asuborderofsynbranchiformteleostfishes.PartII:hylogeneticanalysis. Bull Br Mus Nat Hist Zool, 1984, 47(2): 83-150.

[18]KottelatM.NotesonthetaxonomyanddistributionofsomewesternIndonesianfreshwaterfishes,withdiagnosesofanewgenusandsixnewspecies(Pisces:Cyprinidae,Belonidae,andChaudhuriidae). Ichthyol Explor Fres, 1991, 2(3): 273-287.

[19]JohnsonGD,PattersonC.Percomorphphylogeny:AsurveyofAcanthomorphsandanewproposal. Bull Mar Sci, 1993, 52(1): 554-626.

[20]KottelatM,LimKP.DiagnosisoftwonewgeneraandthreenewspeciesofearthwormeelsfromtheMalayPeninsulaBorneo(Teleostei:Chaudhuriidae). Ichthyol Explor Fres, 1994, 5(2):181-190.

[21]BritzR.OntogenyofethmoidalregionandhyopalatinearchinMacrognathus pancalus(Percomorpha,Mastacembeloidei),withcriticalremarksonmastacembeloidinter-andintrarelationships. Am Mus Novit, 1996, 3181: 1-18.

[22]VrevenEJ.Mastacembelidae(Teleostei;Synbranchiformes)subfamilydivisionandAfricangenericdivision:Anevaluation. J Nat Hist, 2005, 39(4):351-370.

[23]HallTA.BioEdit:auser-friendlybiologicalsequencealignmenteditorandanalysisprogramforWindows95/98/NT. Nucl Acids Symp Ser, 1999, 41: 95-98.

[24]TamuraK,StecherG,PetersonDet al.MEGA6:Molecularevolutionarygeneticsanalysisversion6.0. Mol Biol Evol, 2013, 30(4): 2725-2729.

[25]RozasJ,Anchez-DelBarrioJC,EsseguerXMet al.DnaSP,DNApolymorphismanalysesbythecoalescentandothermethods. Bioinformatics, 2003, 19(18):2496-2497.

[26]ExcoffierL,LavalG,SchneiderS.Arlequinver3.01:Anintegratedsoftwarepackageforpopulationgeneticsdataanalysis.Berne:ComputationalandMolecularPopulationGeneticsLaboratory(CMPG),Switzerland:UniversityofBerne, 2006.

[27]WolstenholmeDR.AnimalmitochondrialDNA:Structureandevolutionin:Mitochondrialgenomes.SanDiego:AcademicPress， 1992:173-372.

[28]NiuSufang,SuYongquan,WangJunet al.GeneticpolymorphismofmitochondrialcontrolregionandcytbinTrachurus japonicusfromFujiancoastalwaters. J Fish Sci China, 2011, 18(1): 66-74(inChinesewithEnglishabstract).[牛素芳, 蘇永全, 王軍等. 福建近海竹莢魚線粒體DNA控制區(qū)和細胞色素b遺傳多態(tài)性. 中國水產(chǎn)科學, 2011, 18(1): 66-74.]

[29]ZhaoKai,HeShunping,PengZuoganget al.MitochondrialdiversityandpopulationstructureofGymnocypris przewalskii przewalskii (Teleostei:Cyprinidae). J Qinghai Univ (Nat Sci), 2006, 24(4): 1-4(inChinesewithEnglishabstract).[趙凱, 何舜平, 彭作剛等. 青海湖裸鯉的種群結(jié)構(gòu)和線粒體DNA變異. 青海大學學報, 2006, 24(4): 1-4.]

[30]ZhangDongya,ChenYong,LiuShaopinget al.PopulationgeneticstructureoftheendangeredspeciesEpalzeorhynchus bicornisintheNujiangRiver. Freshw Fish, 2009, 39(2): 3-7(inChinesewithEnglishabstract).[張東亞, 陳勇, 劉紹平等. 怒江瀕危魚類角魚種群遺傳結(jié)構(gòu)研究. 淡水漁業(yè), 2009, 39(2): 3-7.]

[31]BillingtonN,HebertPD.MitochondrialDNAdiversityinfishesanditsimplicationsforintroductions. Canadian Journal of Fisheries and Aquatic Sciences, 1991, 48(suppl.1): 8-94.

[32]ShakleeJB,TamaruCS,WaplesRS.Speciationandevolutionofmarinefishesstudiedbytheelectrophoreticanalysisofproteins. Pac Sci, 1982, 36(2): 141-157.

[33]ChenDi.AnalysisofgeneticstructureinmitochondrialcytochromebofCoilia ectenesinYangtzeRiver[Dissertation].Guangzhou:JinanUniversity, 2006 (inChinesewithEnglishabstract).[陳迪. 長江刀鱭遺傳多樣性的細胞色素b基因序列分析[學位論文]. 廣州: 暨南大學, 2006.]

[34]LüBaozhong,ZhongYang,GaoLipinget al.Molecularevolutionandphylogenetics.Beijing:HighEducationPress, 2002.[呂寶忠, 鐘揚, 高莉萍等譯. 分子進化與系統(tǒng)發(fā)育. 北京: 高等教育出版社, 2002.]

[35]GrantWS,BowenBW.Shallowpopulationhistoriesindeepevolutionarylineagesofmarinefishes:insightsfromsardinesandanchoviesandlessonsforconservation. Journal of Heredity, 1998, 89: 415-426.

Genetic diversity of spiny eel (Mastacembelus sinensis) from five lakes based on Mitochondrial DNA cyt b gene

ZHAO Ziming， YE Jiansheng， LIU Meijian， SHEN Dongdong & HUANG Shiyi

(Jiangsu Agri-animal Husbandry Vocational College, Taizhou 225300, P.R.China)

ThemitochondrialDNAcyt bgene(767bp)ofspinyeelfromfivelakesinJiangsuProvincewasamplifiedandsequenced.TheaveragecontentsofA+TwerealittlehigherthanthoseofG+Cwhen38variablesitesweredetected,and28haplotypeswererecovered.Theaveragehaplotypediversity(0.58751)andnucleotidediversity(0.00214)demonstratedverylowgeneticdiversityofthespinyeelpopulationsinthefivelakes.Thefixationindiceofanalysisofmolecularvariancewas0.0356,whichshowednosignificantgeneticdifferenceamongpopulations.TheUPGMAphylogenetictreeofcyt bgenebuiltwithMEGA6.0showedthatfivepopulationswerenotclusteredintocladesbasedongeneticdistance.TherewasanextensivegeneralflowamongspinyeelpopulationsfromthefivelakesinJiangsuProvince.Theresultshowedthatspinyeelfromthefivelakeshadlowgeneticdiversityandgeneticallyhomogeneouspopulationstructure,indicatingtheirweakviabilitytomaintainpopulation.Thespinyeelneedtobewellprotected.

Mastacembelus sinensis;MitochondrialDNA; cyt b;geneticdiversity;freshwaterlakes

*江蘇農(nóng)牧科技職業(yè)學院課題(NSFRC1304)資助. 2015-08-29收稿; 2015-11-13收修改稿. 趙子明(1962～), 男, 副教授;E-mail:zzm0282678@163.com.