基于線粒體DNA D-loop的都柳江鲇形目?jī)煞N經(jīng)濟(jì)魚(yú)類種群遺傳多樣性研究

劉 偉，代應(yīng)貴，袁振興，孫際佳，劉 麗

( 1. 貴州大學(xué) 動(dòng)物科學(xué)學(xué)院，貴州 貴陽(yáng) 550025； 2. 貴州大學(xué) 特種水產(chǎn)研究所，貴州 貴陽(yáng) 550025；3. 華南農(nóng)業(yè)大學(xué) 動(dòng)物科學(xué)學(xué)院，廣東 廣州 510642 )

基于線粒體DNA D-loop的都柳江鲇形目?jī)煞N經(jīng)濟(jì)魚(yú)類種群遺傳多樣性研究

劉 偉1,2，代應(yīng)貴1,2，袁振興1,2，孫際佳3，劉 麗3

( 1. 貴州大學(xué) 動(dòng)物科學(xué)學(xué)院，貴州 貴陽(yáng) 550025； 2. 貴州大學(xué) 特種水產(chǎn)研究所，貴州 貴陽(yáng) 550025；3. 華南農(nóng)業(yè)大學(xué) 動(dòng)物科學(xué)學(xué)院，廣東 廣州 510642 )

以采自貴州都柳江的鲇魚(yú)、斑鳠各39尾個(gè)體為樣品，采用PCR和DNA測(cè)序技術(shù)，研究了這兩種魚(yú)類線粒體DNA D-loop的結(jié)構(gòu)和種群遺傳多樣性。獲得了鲇魚(yú)、斑鳠D-loop 5′端長(zhǎng)度分別為544～545 bp、546～547 bp的DNA序列，并分別檢測(cè)到58、15個(gè)變異位點(diǎn)。同時(shí)，識(shí)別了該2種魚(yú)類D-loop終止序列區(qū)和中央保守區(qū)的核心序列，在中央保守區(qū)之后均有一個(gè)Poly-T結(jié)構(gòu)。都柳江鲇魚(yú)、斑鳠種群D-loop序列中，分別有51、12個(gè)多態(tài)位點(diǎn)，并分別定義了11、8種單倍型。都柳江鲇魚(yú)、斑鳠種群?jiǎn)伪缎投鄻有灾笖?shù)、核苷酸多樣性指數(shù)和平均核苷酸差異數(shù)分別為 0.772、0.02301、12.451和0.331、0.00131、0.714。目前，都柳江鲇魚(yú)種群遺傳多樣性較為豐富，具有良好的種質(zhì)價(jià)值和保護(hù)前景; 斑鳠種群遺傳多樣性貧乏，保護(hù)和恢復(fù)該種群的數(shù)量和遺傳多樣性刻不容緩。

都柳江；鲇魚(yú)；斑鳠；D-loop；遺傳多樣性

都柳江屬珠江水系，發(fā)源于貴州獨(dú)山縣，流經(jīng)貴州境內(nèi)三都縣、榕江縣和從江縣，全長(zhǎng)360 km，為柳江的源流河段[1]。鲇魚(yú)(Silurusasotus)、斑鳠(Mystusguttatus)分別隸屬于鲇形目鲇科、鲿科，其中鲇魚(yú)廣泛分布于我國(guó)各大水系(除青藏高原及新疆外)，斑鳠分布于珠江、元江及錢塘江等水系[2]。這兩種魚(yú)類均為我國(guó)重要淡水經(jīng)濟(jì)魚(yú)類，深受當(dāng)?shù)叵M(fèi)者喜愛(ài)[2-3]。近年來(lái)，由于過(guò)度的捕撈、壞境污染等因素的影響，這兩種魚(yú)類野生資源遭到了一定程度的破壞。為了更好地保護(hù)和利用其資源，開(kāi)展這兩種魚(yú)類種群遺傳多樣性研究具有重要意義。近年來(lái)，已有國(guó)內(nèi)學(xué)者分別采用RAPD[4]、Cyt b基因[5]、微衛(wèi)星DNA[6]等分子標(biāo)記研究了鲇魚(yú)不同種群的遺傳多樣性。有關(guān)斑鳠種群遺傳多樣的研究報(bào)道甚少，僅見(jiàn)陳賽等[7]利用同工酶標(biāo)記對(duì)珠江水系西江斑鳠種群遺傳多樣性的研究。迄今為止，尚無(wú)貴州都柳江鲇魚(yú)、斑鳠種群遺傳多樣性的研究報(bào)道。

在進(jìn)行物種遺傳分化及遺傳多樣性研究時(shí)，DNA序列分析比其他分子標(biāo)記如RELP、RAPD等更為直接和可靠[8]。線粒體DNA因其具有分子小、結(jié)構(gòu)簡(jiǎn)單、多態(tài)性高及母系遺傳等特點(diǎn)而常被用作遺傳多樣性研究的分子標(biāo)記[9]。線粒體DNA D-loop，又稱線粒體DNA控制區(qū)，是mtDNA的非編碼區(qū)，不受自然選擇壓力的影響，具有較高的突變積累，適用于進(jìn)行種內(nèi)、種群或個(gè)體間遺傳多樣性的研究[10]。本文擬對(duì)貴州境內(nèi)都柳江鲇魚(yú)、斑鳠種群mtDNA D-loop序列進(jìn)行測(cè)序分析，研究該序列的結(jié)構(gòu)特征和種群遺傳多樣性，以期為開(kāi)展都柳江鲇形目這兩種經(jīng)濟(jì)魚(yú)類種質(zhì)資源的利用及保護(hù)提供科學(xué)依據(jù)。

1 材料與方法

1.1 試驗(yàn)材料

試驗(yàn)所用鲇魚(yú)、斑鳠各39尾鮮活個(gè)體，于2014年8月至10月，分別采自貴州都柳江干流上、中、下游的獨(dú)山、三都、榕江、從江。每尾活魚(yú)取肌肉3～5 g用乙醇固定，帶回實(shí)驗(yàn)室后于冰箱-20 ℃保存，備用。

1.2 DNA模板制備

基因組DNA的提取采用北京天根生化科技有限公司提供的DNA提取試劑盒，并參照其提供的方法提取DNA，用1%的瓊脂糖凝膠電泳檢測(cè)所提取DNA的質(zhì)量和純度，于冰箱-20 ℃儲(chǔ)存、備用。

1.3 D-loop PCR擴(kuò)增及序列測(cè)定

擴(kuò)增鲇魚(yú)、斑鳠D-loop的引物由上海生物工程有限公司合成，引物序列分別為：D-loop-F: 5′-ACCCCTGGCTCCCAAAGC-3′ 和D-loop-R: 5′-ATCTTAGCATCTTCAGTG-3′[11-12]。

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

對(duì)測(cè)得序列進(jìn)行Blastn在線同源性檢測(cè)，確定序列的邊界及長(zhǎng)度，并用SeqMan軟件進(jìn)行比對(duì)，輔以人工校對(duì)。采用DnsSP(version 5.10)軟件計(jì)算遺傳多樣性參數(shù)。用MEGA 5.0軟件統(tǒng)計(jì)序列的堿基含量并用Kimura雙參數(shù)法計(jì)算遺傳距離。分別以GenBank中大口鲇(S.meridionalis)mtDNA D-loop同源序列(登錄號(hào)：NC_014866.1)、大鰭鳠(M.macropterus)mtDNA D-loop同源序列(登錄號(hào)：NC_019592.1)作為外類群，構(gòu)建都柳江鲇魚(yú)、斑鳠種群的鄰接系統(tǒng)樹(shù)。采用劉煥章等[13-15]的方法進(jìn)行鲇魚(yú)、斑鳠D-loop的結(jié)構(gòu)分析。

2 結(jié) 果

2.1 都柳江鲇魚(yú)、斑鳠種群的mtDNA D-loop序列堿基組成及多態(tài)性

PCR擴(kuò)增產(chǎn)物經(jīng)單向測(cè)序，再去除所獲序列測(cè)序峰圖中兩端波峰混亂的序列后，得到了都柳江鲇魚(yú)、斑鳠各39尾個(gè)體mtDNA D-loop 5′端長(zhǎng)度分別為544～545 bp、546～547 bp 的DNA序列，用于其各自種群遺傳多樣性的研究。結(jié)果顯示，都柳江鲇魚(yú)、斑鳠種群D-loop序列堿基組成相似，且均表現(xiàn)為A+T占?jí)A基總數(shù)的比例顯著高于G+C的比例，并均以A所占比例最高、以G所占比例最低，具有明顯的堿基偏倚性(表1)。同時(shí)，該兩個(gè)種群所獲D-loop序列末端均有一個(gè)poly-T結(jié)構(gòu)。

表1 都柳江鲇魚(yú)、斑鳠種群mtDNA D-loop堿基組成(n=39) %

都柳江鲇魚(yú)種群39尾個(gè)體mtDNA D-loop序列長(zhǎng)度為544～545 bp，有變異位點(diǎn)58個(gè)，變異位點(diǎn)數(shù)占分析序列長(zhǎng)度的10.64%～10.66%(表2)。其中，有51個(gè)簡(jiǎn)約信息位點(diǎn)(含2個(gè)簡(jiǎn)約信息位點(diǎn)同時(shí)存在缺失：206 nt位點(diǎn)為T→C轉(zhuǎn)換同時(shí)存在缺失T，536 nt位點(diǎn)為A→T顛換同時(shí)存在缺失A)、2個(gè)單一變異位點(diǎn)、4個(gè)插入位點(diǎn)(135 nt、213 nt、246 nt位點(diǎn)插入A，202 nt位點(diǎn)插入C)和3個(gè)缺失位點(diǎn)(177 nt位點(diǎn)缺失A，206 nt位點(diǎn)缺失T，536 nt位點(diǎn)缺失A)?？梢?jiàn)，該序列共有51個(gè)多態(tài)位點(diǎn)。同時(shí)，該序列含有47個(gè)轉(zhuǎn)換位點(diǎn)和7個(gè)顛換位點(diǎn)， 其中1 個(gè)位點(diǎn)同時(shí)存在轉(zhuǎn)換與顛換(244 nt，簡(jiǎn)約信息位點(diǎn)：C→T、C→A)，平均轉(zhuǎn)顛換比為6.71。

都柳江斑鳠種群39個(gè)個(gè)體mtDNA D-loop基因序列長(zhǎng)度為546～547 bp，共有15個(gè)變異位點(diǎn)，變異位點(diǎn)數(shù)占分析序列長(zhǎng)度的2.74%～2.75%(表3)。其中含1個(gè)簡(jiǎn)約信息位點(diǎn)、11個(gè)單一變異位點(diǎn)、2個(gè)插入位點(diǎn)(20 nt位點(diǎn)插入A、44 nt位點(diǎn)插入T)和1個(gè)缺失位點(diǎn)(52 nt位點(diǎn)缺失C)，共有12個(gè)多態(tài)位點(diǎn)。該序列含有11個(gè)轉(zhuǎn)換位點(diǎn)和2個(gè)顛換位點(diǎn)，其中一個(gè)位點(diǎn)同時(shí)存在轉(zhuǎn)換和顛換(111 nt，單一變異位點(diǎn)：G→A、G→C)，平均轉(zhuǎn)顛換比約為5.5。

不計(jì)插入或缺失位點(diǎn)，都柳江鲇魚(yú)、斑鳠種群分別定義了11、8個(gè)單倍型(表2、3)。在鲇魚(yú)種群11個(gè)單倍型中，單倍型1 (Hap_1) 個(gè)體數(shù)最多、達(dá)18個(gè)，為優(yōu)勢(shì)單倍型。在斑鳠種群8個(gè)單倍型中，單倍型2 (Hap_2)占個(gè)體數(shù)最多、達(dá)31個(gè)，其余單倍型均為一個(gè)個(gè)體。

表2 都柳江鲇魚(yú)種群mtDNA D-loop的多態(tài)位點(diǎn)及單倍型

注: 圓點(diǎn)(.)表示與 Hap_1 相同的堿基.表3同.

表3 都柳江斑鳠種群mtDNA D-loop的多態(tài)位點(diǎn)及單倍型

2.2 都柳江鲇魚(yú)、斑鳠mtDNA D-loop的結(jié)構(gòu)

采用劉煥章等[13-15]的方法，識(shí)別了都柳江鲇魚(yú)、斑鳠D-loop終止序列區(qū)和中央保守區(qū)，但未能識(shí)別保守序列區(qū)。

2.2.1 鲇魚(yú)mtDNA D-loop的結(jié)構(gòu)

在都柳江鲇魚(yú)種群mtDNA D-loop序列中，識(shí)別了終止序列區(qū)2個(gè)終止相關(guān)序列，分別為ETAS1(TACATAATATGCATAA-CTGG，“-”表示發(fā)生變異的堿基，下同)和ETAS2(-ACATATTATGTATAATAT)(圖1)。在終止序列區(qū)共有7個(gè)核心序列(TACAT)和1個(gè)與核心序列互補(bǔ)的序列(ATGTA)。終止序列區(qū)長(zhǎng)度為250～251 bp，含有38個(gè)變異位點(diǎn)，變異率為15.14%～15.20%。同時(shí)，還識(shí)別了中央保守區(qū)的CSB-F、CSB-E、CSB-D等核心序列，分別為ATGTTTAGTAGTGAGAGACCATCAACC、TCAGGGACAAAACTTGTGGGGGT、TATTACTGGCATC

TGGTTCCTATTTCA，其中CSB-E包含GTGGG-box。中央保守區(qū)序列長(zhǎng)251～252 bp，含17個(gè)變異位點(diǎn)，變異率為6.75%～6.77%。在CSB-D后有一個(gè)poly-T結(jié)構(gòu)(TATTTTTTTTTTT)。

2.2.2 斑鳠mtDNA D-loop的結(jié)構(gòu)

在都柳江斑鳠mtDNA D-loop序列中，識(shí)別了終止序列區(qū)(ETAS)1個(gè)終止相關(guān)序列：TACATAATATGTATAATAT(圖2)。在終止序列區(qū)，有6個(gè)核心序列(TACAT)和1個(gè)與其反向互補(bǔ)的序列(ATGTA)。終止序列區(qū)長(zhǎng)232～233 bp，共有9個(gè)變異位點(diǎn)，變異率3.86%～3.88%。此外，還識(shí)別了中央保守區(qū)3個(gè)核心序列，即CSB-F(ATGTAGTAAGAAACCAGCAACCCTCAATA)、CSB-E(TCAGGGACAA-AATTGTGAGGGTT)、CSB-D(TATTACTGGCATCTGGTTCCTAT)。中央保守區(qū)序列長(zhǎng)度為252 bp，含有6個(gè)變異位點(diǎn)，變異率2.38%。在CSB-D后也存在一個(gè)poly-T結(jié)構(gòu)(TATTTTTTTTTTTT)。

圖1 都柳江鲇魚(yú)mtDNA D-loop部分序列

圖2 都柳江斑鳠mtDNA D-loop部分序列

2.3 都柳江鲇魚(yú)、斑鳠種群的遺傳多樣性及種群遺傳結(jié)構(gòu)

用DnaSP 5.10軟件對(duì)都柳江鲇魚(yú)、斑鳠種群mtDNA D-loop的遺傳多樣性進(jìn)行了分析。結(jié)果表明，鲇魚(yú)種群各項(xiàng)遺傳多樣性指數(shù)均明顯高于斑鳠種群(表4)。

表4 都柳江鲇魚(yú)、斑鳠種群mtDNA遺傳多樣性指數(shù)

應(yīng)用MEGA 5.0 軟件的Kimura 雙參數(shù)模型計(jì)算都柳江鲇魚(yú)、斑鳠種群?jiǎn)伪缎椭g的遺傳距離。結(jié)果顯示，鲇魚(yú)種群11個(gè)單倍型之間的遺傳距離為0.0019～0.0837，平均遺傳距離為0.0343；斑鳠種群8個(gè)單倍型之間的遺傳距離為0.0018～0.0167，平均遺傳距離為0.0063?？梢?jiàn)，鲇魚(yú)種群?jiǎn)伪缎椭g的平均遺傳距離明顯大于斑鳠種群。

構(gòu)建的鄰接系統(tǒng)樹(shù)表明，鲇魚(yú)11個(gè)單倍型聚為兩分支(Gp1和Gp2)，其中Gp1包括9個(gè)單倍型，置信值為99，代表34個(gè)個(gè)體，為優(yōu)勢(shì)分支;Gp2包括2個(gè)單倍型，代表5個(gè)個(gè)體，該分支置信值高達(dá)100(圖3)。2個(gè)分支間平均遺傳距離為0.0745。斑鳠8個(gè)單倍型分為兩支，其中一支由7個(gè)單倍型組成，代表38個(gè)個(gè)體，置信值為82；另一支僅1個(gè)單倍型Hap_8，代表一個(gè)個(gè)體(圖4)。

圖3 基于都柳江鲇魚(yú)種群mtDNA單倍型間遺傳距離構(gòu)建的鄰接系統(tǒng)樹(shù)注：分支處的數(shù)字為Bootstrap值，重復(fù)次數(shù)為1000次；Hap_1～Hap_11代表不同的單倍型.

圖4 基于都柳江斑鳠種群mtDNA單倍型間遺傳距離構(gòu)建的鄰接系統(tǒng)樹(shù)注：分支處的數(shù)字為Bootstrap值，重復(fù)次數(shù)為1000次；Hap_1～Hap_8代表不同的單倍型.

3 討 論

3.1 都柳江鲇魚(yú)、斑鳠種群mtDNA D-loop的序列變異

線粒體DNA D-loop區(qū)為非編碼區(qū)，不受自然選擇壓力的影響，存在較多的突變?nèi)鐗A基替換、缺失、插入及串聯(lián)重復(fù)等，是mtDNA中變異最多的區(qū)域[16]。本研究測(cè)定了都柳江鲇魚(yú)種群、斑鳠種群mtDNA D-loop部分序列，分別存在58、15個(gè)變異位點(diǎn)，變異類型包括堿基替換、插入和缺失，并存在同一位點(diǎn)同時(shí)發(fā)生轉(zhuǎn)換和顛換、缺失和轉(zhuǎn)換或缺失和顛換的現(xiàn)象。都柳江鲇魚(yú)、斑鳠種群D-loop堿基轉(zhuǎn)顛換比值分別為6.71、5.5，轉(zhuǎn)換頻率均大于顛換頻率，這與脊椎動(dòng)物DNA中發(fā)生轉(zhuǎn)換的頻率通常高于顛換的規(guī)律[17]相符。

3.2 都柳江鲇魚(yú)、斑鳠mtDNA D-loop的結(jié)構(gòu)

魚(yú)類mtDNA D-loop區(qū)結(jié)構(gòu)較為復(fù)雜，一般可分為3個(gè)區(qū)域：終止序列區(qū)、中央保守區(qū)和保守序列區(qū)[14]。本研究識(shí)別了都柳江鲇魚(yú)、斑鳠種群D-loop的終止序列區(qū)和中央保守區(qū)的3個(gè)核心序列(CSB-F、CSB-E、CSB-D)。與GenBank中鲇魚(yú)(登錄號(hào)：JX256247.1)、斑鳠(登錄號(hào)：NC_023976.1)mtDNA D-loop全序列的比對(duì)表明，本研究中鲇魚(yú)、斑鳠D-loop保守序列均應(yīng)在其各自所獲序列末端poly-T結(jié)構(gòu)之后而已被去除，因而未能識(shí)別出其各自的保守序列。

都柳江鲇魚(yú)、斑鳠種群D-loop終止序列區(qū)堿基變異率均高于其各自中央保守區(qū)的堿基變異率，可見(jiàn)終止序列區(qū)表現(xiàn)了較高的變異性。此外，本研究還識(shí)別了都柳江鲇魚(yú)種群D-loop終止序列區(qū)的2個(gè)終止相關(guān)序列——ETAS1、ETAS2。同樣，朱世華等[15, 18]在鲹科魚(yú)類和斑點(diǎn)叉尾(Ictaluruspunctatus)控制區(qū)中也發(fā)現(xiàn)了2個(gè)ETAS序列。當(dāng)魚(yú)類mtDNA D-loop中存在多個(gè)終止相關(guān)序列時(shí)，劉煥章等[13]認(rèn)為其中只有1個(gè)ETAS可行使功能，而其他ETAS只是復(fù)制的結(jié)果，不行使功能。都柳江斑鳠種群D-loop僅識(shí)別了1個(gè)ETAS，其核心識(shí)別序列與鳑鲏亞科魚(yú)類[13]、鲿科魚(yú)類[14]略有不同，與同屬的大鰭鳠[19]則差異較大。此外，還識(shí)別了鲇魚(yú)、斑鳠中央保守區(qū)3個(gè)核心序列(CSB-F、CSB-E、CSB-D)。其中，除鲇魚(yú)的CSB-F與鲿科魚(yú)類差異較大外，其他核心序列則與鲿科魚(yú)類基本一致[13]。本研究中，在鲇魚(yú)、斑鳠D-loop區(qū)CSB-D后均存在一個(gè)poly-T結(jié)構(gòu)(TATTTTTTTTTTT)，在中國(guó)野生香魚(yú)(Plecoglossusaltivelis)[20]、大鰭鳠[21]和大口(Psettodeserumel)[22]中也發(fā)現(xiàn)了類似的結(jié)構(gòu)，其功能尚不清楚，還有待研究。

3.3 都柳江鲇魚(yú)、斑鳠種群的遺傳多樣性

本研究基于mtDNA D-loop序列，選取多態(tài)位點(diǎn)數(shù)、單倍型數(shù)、單倍型多樣性指數(shù)、核苷酸多樣性指數(shù)和平均核苷酸差異數(shù)等5個(gè)參數(shù)研究了都柳江鲇魚(yú)、斑鳠種群的遺傳多樣性。其中，單倍型多樣性指數(shù)和核苷酸多樣性指數(shù)是衡量一個(gè)種群遺傳多樣性高低的2個(gè)重要指標(biāo)。由于核苷酸多樣性指數(shù)考慮了各單倍型在群體中所占的比例，因而較單倍型多樣性指數(shù)更為可靠[23]。本研究測(cè)定了都柳江鲇魚(yú)、斑鳠各39個(gè)個(gè)體mtDNA的D-loop序列，分別定義了11和8個(gè)單倍型。都柳江鲇魚(yú)種群?jiǎn)伪缎投鄻有灾笖?shù)、核苷酸多樣性指數(shù)分別為0.772、0.02301，而斑鳠種群?jiǎn)伪缎投鄻有灾笖?shù)、核苷酸多樣性指數(shù)分別為0.331、0.00131。可見(jiàn)，都柳江鲇魚(yú)種群遺傳多樣性高于斑鳠種群。

都柳江鲇魚(yú)種群?jiǎn)伪缎投鄻有灾笖?shù)略低于長(zhǎng)江水系、珠江水系和黑龍江水系8個(gè)草魚(yú)(Ctenopharyngodonidella)野生群體[24]、長(zhǎng)江中游5個(gè)青魚(yú)(Mylopharyngodonpiceus)群體[25]、長(zhǎng)江中游湖泊黃顙魚(yú)(Pelteobagrusfulvidraco)群體[26]，但其核苷酸多樣性指數(shù)均高于上述3種魚(yú)類群體。而都柳江斑鳠種群?jiǎn)伪缎投鄻有灾笖?shù)、核苷酸多樣性指數(shù)均低于上述3種魚(yú)類群體。此外，都柳江鲇魚(yú)種群?jiǎn)伪缎投鄻有灾笖?shù)、核苷酸多樣性指數(shù)均顯著高于同分布于都柳江的易危魚(yú)類小口白甲魚(yú)(Onychostomalini)種群[27]，而都柳江斑鳠種群?jiǎn)伪缎投鄻有灾笖?shù)略高于都柳江小口白甲魚(yú)種群但核苷酸多樣性指數(shù)則低于都柳江小口白甲魚(yú)種群。因此，目前都柳江鲇魚(yú)種群尚保存著較豐富的遺傳多樣性，而斑鳠種群遺傳多樣性則較為貧乏。鲇魚(yú)、斑鳠均為我國(guó)的名優(yōu)淡水經(jīng)濟(jì)魚(yú)類，其中鲇魚(yú)還為我國(guó)重要的養(yǎng)殖魚(yú)類，斑鳠為極具馴養(yǎng)開(kāi)發(fā)價(jià)值的野生魚(yú)類。都柳江鲇魚(yú)種群較高的遺傳多樣性，表明該種群具有良好的種質(zhì)價(jià)值和較好的保護(hù)前景。都柳江斑鳠種群目前遺傳多樣性貧乏，保護(hù)和恢復(fù)該種群的數(shù)量和遺傳多樣性刻不容緩。

Grant等[28]認(rèn)為，可以根據(jù)種群?jiǎn)伪缎投鄻有灾笖?shù)和核苷酸多樣性指數(shù)值推測(cè)該種群的歷史，即單倍型多樣性指數(shù)和核苷酸多樣性指數(shù)較高的種群可能由一個(gè)大而穩(wěn)定的種群經(jīng)過(guò)長(zhǎng)時(shí)間演化所產(chǎn)生或由兩個(gè)不同系群的種群二次接觸所形成，而單倍型多樣性指數(shù)和核苷酸多樣性指數(shù)較低的種群最近可能發(fā)生了瓶頸效應(yīng)或奠基者效應(yīng)。都柳江鲇魚(yú)種群較高的單倍型多樣性指數(shù)和核苷酸多樣性指數(shù)，這表明該種群可能源于大而穩(wěn)定的種群并且目前仍維持著較高的遺傳多樣性。斑鳠為西江“四大名魚(yú)”之一[3]，經(jīng)濟(jì)價(jià)值高，近年來(lái)因過(guò)度捕撈導(dǎo)致其種群數(shù)量急劇減少，進(jìn)而降低了都柳江斑鳠種群?jiǎn)伪缎投鄻有灾笖?shù)和核苷酸多樣性指數(shù)，推測(cè)該種群近期可能發(fā)生了瓶頸效應(yīng)。

[1] 伍律. 貴州魚(yú)類志[M]. 貴陽(yáng): 貴州人民出版社,1989.

[2] 褚新洛,鄭葆珊,戴定遠(yuǎn),等. 中國(guó)動(dòng)物志: 硬骨魚(yú)綱鲇形目[M]. 北京:科學(xué)出版社,1999：77-87.

[3] 吳小明,周立斌,彭飛影. 斑鳠的含肉率及肌肉營(yíng)養(yǎng)成分分析[J]. 水產(chǎn)科學(xué),2010,29(7):425-428.

[4] 王朝明,鄒桂偉,鄭蓓蓓,等. 三種鲇遺傳多樣性的RAPD分析[J]. 淡水漁業(yè),2005,35(4):14-17.

[5] 肖明松,崔峰,康健,等. 淮河野生鲇魚(yú)線粒體cyt b基因的序列變異與遺傳結(jié)構(gòu)分析[J]. 動(dòng)物學(xué)雜志,2013,48(1):75-86.

[6] 徐丹丹. 基于微衛(wèi)星標(biāo)記和線粒體cyt b基因序列的鲇遺傳多樣性研究[D]. 重慶：西南大學(xué),2013.

[7] 陳賽,焦宗垚,朱新平,等. 西江野生斑鳠乳酸脫氫酶和酯酶同工酶的電泳研究[J]. 湛江海洋大學(xué)學(xué)報(bào):自然科學(xué)版,2006,26(6):58-62.

[8] Buonnacorsi V P, Mcdowell J R, Graves J E. Reconciling patterns of interocean molecular variance from four classes of molecular markers in blue marlin (Makairanigricans) [J]. Mol Eco,2001,10(5):1179-1196.

[9] 郭新紅,劉少軍,劉巧,等. 魚(yú)類線粒體DNA研究新進(jìn)展[J]. 遺傳學(xué)報(bào),2004,31(9):983-1000.

[10] Thacker C E. Molecular phylogeny of the gobioid fishes (Teleostei:Perciformes:Gobioidei)[J]. Mol Phyl Evol,2003,26(3):354-368.

[11] Zeng Q,Wang Z,Peng Z. Mitochondrial genome ofSilurusasotus(Teleostei:Siluriformes)[J]. Mitochondrial DNA, 2011,22(5/6):162-164.

[12] Tian H,Que Y,Zhao N,et al. The complete mitochondrial genome of the spotted longbarbel catfish,Hemibagrusguttatus(Siluriformes, Bagridae)[J]. Mitochondrial DNA, 2016,27(1):467-468.

[13] 劉煥章. 魚(yú)類線粒體控制區(qū)的結(jié)構(gòu)和進(jìn)化:以鳑鲏?mèng)~類為例[J]. 自然科學(xué)進(jìn)展,2002,12(3):266-270.

[14] 張燕,張鶚,何舜平. 中國(guó)鲿科魚(yú)類線粒體DNA 控制區(qū)結(jié)構(gòu)及其系統(tǒng)發(fā)育分析[J]. 水生生物學(xué)報(bào),2003,27(5):463-467.

[15] 朱世華,鄭文娟,鄒記興,等. 鲹科魚(yú)類線粒體DNA控制區(qū)結(jié)構(gòu)及系統(tǒng)發(fā)育關(guān)系[J]. 動(dòng)物學(xué)研究,2007,28(6):606-614.

[16] Xiao Y S, Zhang Y, Gao T X, et al. Genetic diversity in the mtDNA control region and population structure in the small yellow croakerLarimichthyspolyactis[J]. African Journal of Biotechnology, 2009,85(4):303-314.

[17] Loftus R T, Machugh D E, Bradley D G, et al. Evidence for two independent domestications of cattle[J]. Proceedings of the National Academy of Sciences of the United States of America, 1994, 91(7):2575-2761.

[18] 秦欽,許志強(qiáng),邊文冀,等. 斑點(diǎn)叉尾線粒體DNA控制區(qū)結(jié)構(gòu)和群體遺傳多樣性分析[J].淡水漁業(yè),2011,41(4):50-54.

[19] 魏剛,周麗,曾令江,等. 大鰭鳠線粒體控制區(qū)結(jié)構(gòu)及其比較分析[J]. 西南大學(xué)學(xué)報(bào)：自然科學(xué)版,2015,37(1):46-52.

[20] 樂(lè)小亮. 中國(guó)野生香魚(yú)遺傳多樣性分析[D]. 廣州：暨南大學(xué),2010.

[21] 周麗. 大鰭鳠遺傳多樣性及其線粒體控制區(qū)結(jié)構(gòu)的研究[D]. 重慶：西南大學(xué),2008.

[22] 張艷春,孔曉瑜,王忠明. 大口鳒線粒體DNA控制區(qū)結(jié)構(gòu)和鰈形目魚(yú)類的系統(tǒng)進(jìn)化初步研究[J]. 熱帶海洋學(xué)報(bào),2010,29(6):71-78.

[23] Zhou H,Li D Q,Zhang Y G,et al. Study on mitochondrial DNA genetic diversity of Tibetan antelope[J]. Yi Chuan,2006,28(3):299-305.

[24] 傅建軍,王榮泉,沈玉幫,等. 我國(guó)草魚(yú)野生群體D-Loop序列遺傳變異分析[J]. 水生生物學(xué)報(bào),2015,39(2):349-357.

[25] 楊宗英,汪登強(qiáng),陳大慶,等. 基于mtDNA序列分析青魚(yú)群體遺傳結(jié)構(gòu)[J]. 淡水漁業(yè),2015,45(2):3-7.

[26] 劉朋朋,鐘立強(qiáng),潘建林,等. 基于線粒體DNA D-loop分析黃顙魚(yú)五個(gè)淡水湖泊群體的遺傳多樣性[J]. 海洋與湖沼,2013,44(3):728-733.

[27] 代應(yīng)貴,韓雪,張曉杰. 小口白甲魚(yú)都柳江種群mtDNA D環(huán)的序列變異及遺傳多樣性[J].動(dòng)物學(xué)雜志,2010,45(2):115-120.

[28] Grant W S,Bowen B W. Shallow population histories in deep evolutionary lineages of marine fishes: insights from sardines and anchovies and lessons for conservation[J]. The Journal of Heredity,1998,89(5):415-426.

GeneticDiversityof2EconomicallyImportantFishesinSiluriformesfromDuliuRiverBasedonmtDNAD-loopSequences

LIU Wei1, 2, DAI Yinggui1, 2, YUAN Zhenxing1, 2,SUN Jijia3, LIU Li3

( 1. College of Animal Sciences, Guizhou University, Guiyang 550025, China;2. Special Fisheries Research Institute, Guizhou University, Guiyang 550025,China;3. College of Animal Sciences, Agricultural University of South China, Guangzhou 510642, China )

The composition of mtDNA D-loop and genetic diversity were studied in 39 individuals of catfishesSilurusasotusandMystusguttatuscollected in the Guizhou section of Duliu River by the methods of PCR and DNA sequencing on the basis of sequences of mtDNA D-loop.There were 544—545 bp and 546—547 bp of D-loop sequences near the 5′ end of D-loop, including 58 variable loci inS.asotusand 15 variable loci inM.guttatus. The characteristic sequences were identified in the extended termination associated sequence and the central conserved domain and a poly-T next to the 3′ end of central conserved domain in the D-loop sequence in each of the two species.There were 51 and 12 polymorphic loci and 11 and 8 haplotypes in the populations ofS.asotusandM.guttatus, respectively.S.Asotushad haploidtype diversity of 0.772,nucleotide diversity of 0.02301,and average nucleotide differences of 12.451, andM.guttatushad haploidtype diversity of 0.331,nucleotide diversity of 0.00131,and verage nucleotide differences of 0.714 .The population ofS.asotusfrom the Duliu River showed abundant genetic diversity,good germplasm value and protection outlook. While the population ofM.guttatushad much low genetic diversity in the Duliu River. Therefore,it is urgent to protect and recover the genetic diversity ofM.guttatuspopulation from the Duliu River.

Duliu River;Silurusasotus;Mystusguttatus; D-loop; genetic diversity

S917

A

1003-1111(2016)04-0386-07

10.16378/j.cnki.1003-1111.2016.04.013

2015-09-11；

2016-01-12.

國(guó)家公益性行業(yè)(農(nóng)業(yè))科研專項(xiàng)(201303048);貴州省留學(xué)人員科技活動(dòng)項(xiàng)目[黔人項(xiàng)目資助合同(2013)10].

劉偉(1992—)，男，碩士研究生；研究方向：水產(chǎn)動(dòng)物遺傳育種與種質(zhì)資源學(xué).E-mail: 892913900@qq.com.通訊作者：代應(yīng)貴(1968—)，男，教授；研究方向：魚(yú)類學(xué)、漁業(yè)資源與環(huán)境.E-mail: daiygui@163.com.