張鶴千，楊子拓，李桂峰，肖詩斌，孫際佳，楊慧榮，趙會(huì)宏，謝堂暉，韓興鵬，劉 麗

(1.華南農(nóng)業(yè)大學(xué)動(dòng)物科學(xué)學(xué)院，廣東廣州510642;2.中山大學(xué)生命科學(xué)大學(xué)院，廣東廣州510275;3.華南農(nóng)業(yè)大學(xué)動(dòng)物科學(xué)學(xué)院//廣東省農(nóng)業(yè)動(dòng)物基因組學(xué)與分子育種重點(diǎn)實(shí)驗(yàn)室，廣東廣州510642)

黃顙魚Pelteobagrus fulvidraco在分類學(xué)上屬于硬骨魚綱、鲇形目Silurforwes、鲿科Siluridae、黃顙魚屬Pelteobagrus，俗名叫黃刺公、黃臘丁、黃刺骨等［1］。黃顙魚在靜水或緩流的淺灘生活，晝伏夜出，主食底棲脊椎動(dòng)物，多為小魚、水生昆蟲等小型生物。黃顙魚是中國重要的小型經(jīng)濟(jì)魚類，廣泛分布于長江、黃河、珠江、沅江及黑龍江等水域［2－3］。黃顙魚屬在我國主要有4種，分別為黃顙魚Pelteobagrus fulvidraco、瓦氏黃顙魚Pelteobagrus vachelli、光澤黃顙魚Pelteobagrus nitidus、岔尾黃顙魚Pelteobagrus eupogon，其中黃顙魚是我國水產(chǎn)養(yǎng)殖品種中主要的增養(yǎng)殖對(duì)象，其次是瓦氏黃顙魚。珠江流域是由西江、北江、東江及珠江三角洲諸河等四個(gè)水系所組成的復(fù)合流域［4］。隨著經(jīng)濟(jì)社會(huì)的發(fā)展，環(huán)境污染不斷地加劇，濫捕濫撈現(xiàn)象也日趨嚴(yán)重，致使野生黃顙魚的生活環(huán)境受到了嚴(yán)重破壞，野生種質(zhì)資源數(shù)量呈現(xiàn)銳減趨勢(shì)。為了合理有效地保護(hù)野生黃顙魚群體以促進(jìn)淡水漁業(yè)的長久發(fā)展，對(duì)珠江流域黃顙魚野生群體的種質(zhì)資源進(jìn)行研究已日顯重要。

線粒體DNA(Mitochondrial DNA，mtDNA)基因組內(nèi)不同的區(qū)域進(jìn)化速度存在差異，適合不同水平的進(jìn)化研究［5］。近年來，隨著分子生物學(xué)技術(shù)向經(jīng)典分類學(xué)領(lǐng)域的滲透以及魚類線粒體DNA研究技術(shù)的不斷發(fā)展成熟，使得mtDNA成為魚類學(xué)相關(guān)研究的理想的分子標(biāo)記，并在魚類進(jìn)化遺傳學(xué)、種群遺傳結(jié)構(gòu)、分子生態(tài)學(xué)和保護(hù)生物學(xué)等方面的研究中取得了很多有意義的成果，解決了眾多疑難問題。細(xì)胞色素b(Cyt b)基因進(jìn)化速度適中，適合于種間到種內(nèi)水平上的系統(tǒng)發(fā)生研究，在魚類群體遺傳結(jié)構(gòu)與系統(tǒng)發(fā)育關(guān)系中有著廣泛的應(yīng)用［6］。目前，Cyt b基因是研究分子進(jìn)化和系統(tǒng)發(fā)育最有用的基因之一，已廣泛應(yīng)用于各種生物的系統(tǒng)發(fā)育研究中［7－9］，而且還被用于探討分歧時(shí)間較長的高級(jí)分類單元之間的系統(tǒng)發(fā)育關(guān)系［10］。

本文通過測(cè)定與分析珠江流域桂林江段、都勻江段、右江江段、左江江段、東江江段這5個(gè)江段的190尾野生黃顙魚個(gè)體的線粒體細(xì)胞色素b基因序列，探討珠江流域野生黃顙魚的種群遺傳多樣性和進(jìn)化關(guān)系，在維護(hù)不同水系間遺傳多樣性和原有的生態(tài)平衡上具有重要意義。

1 材料與方法

1.1 試驗(yàn)材料

1.1.1 實(shí)驗(yàn)用魚 本實(shí)驗(yàn)所用野生黃顙魚樣本于2012年6月至8月，采集于珠江流域各江段，分別為桂林江段 (GL)、都勻江段 (GZ)、右江江段(YJ)、左江江段 (ZJ)、東江江段 (DJ)，共計(jì)190尾，每個(gè)群體采集樣本數(shù)目見表1。野外現(xiàn)場取背部肌肉于無水乙醇中固定并帶回實(shí)驗(yàn)室，超低溫保存。

1.1.2 主要儀器 PCR反應(yīng)所用DNA擴(kuò)增儀為南京貝登機(jī)電設(shè)備有限公司制造V123530型擴(kuò)增儀，離心機(jī)為德國Eppendorf公司制造5810R型，凝膠成像系統(tǒng)為美國BIORAD公司制造的Bio-Rad Gel-Doc XR system PC型等。

表1 黃顙魚群體采樣信息Table 1 Pelteobagrus fulvidraco groups sampling information

1.2 試驗(yàn)方法

黃顙魚的基因組DNA的抽提實(shí)驗(yàn)采用天根生物海洋動(dòng)物組織基因組DNA提取試劑盒 (TIANamp Marine Animals DNA Kit)，并根據(jù)本實(shí)驗(yàn)室的實(shí)際情況對(duì)操作流程進(jìn)行了部分更改。Cyt b基因擴(kuò)增和測(cè)序的引物為通用引物序列［11］如下:

H15915:5'－CTCCGATCTCCGGATTACAAGAC－3'，

L14724:5'－GACTTGAAAAACCACCGTTG－3'，引物由華大基因科技服務(wù)有限公司合成。

PCR反應(yīng)總體積為50 μL體系，其中模板DNA2 μL，上游引物P1和下游引物P2(10 μmol/L)各 2 μL，2 ×PCR Reaction Mix25 μL，Taq DNA聚合酶 (2.5 U/μL)1 μL，無離子超純水為18 μL。擴(kuò)增條件為:94℃預(yù)變性2 min后94℃變性45 s，58℃退火45 s，72℃延伸90 s，共進(jìn)行35個(gè)循環(huán)，最后72℃延伸7 min，16℃保存。擴(kuò)增之后用w=1%的瓊脂糖凝膠電泳，BIORAD凝膠成像系統(tǒng)照膠并保存圖片。

1.3 測(cè)序與數(shù)據(jù)分析

PCR擴(kuò)增產(chǎn)物直接送華大基因科技服務(wù)有限公司利用正反引物進(jìn)行雙向測(cè)序獲得黃顙魚的Cyt b基因序列列，經(jīng)過DNASTAR Inc軟件拼接后進(jìn)行人工校對(duì)。用軟件Clustal-X對(duì)所獲得拼接序列進(jìn)行對(duì)比分析［12］;用軟件DNASP5.0計(jì)算單倍型多樣性 (haplotype diversity，Hd)、核苷酸多樣性(nucleotide diversity，Pi)和平均核苷酸差異數(shù) (K±SD)3項(xiàng)群體遺傳多樣性指標(biāo)，群體間的分化指數(shù)(F-statistics，F(xiàn)st)以及群體遺傳變異的分子變異等級(jí)分析(analysis of molecular variance'AMOVA);通過公式Nm= ［(1－Fst)－1］/2計(jì)算得到群體間的基因流;MEGA5.0軟件統(tǒng)計(jì)DNA序列的剪輯組成、轉(zhuǎn)換、顛換和平均轉(zhuǎn)換率/顛換率、單突變位點(diǎn)、簡約信息位點(diǎn)、插入/缺失位點(diǎn)數(shù);個(gè)體間、群體間和群體內(nèi)的遺傳變異率由Kimura雙系數(shù)模型(Kimura2-parameter)計(jì)算所得;利用MEGA5.0軟件，最大似然法 (maximum-like-lihood，ML)構(gòu)建聚類關(guān)系樹。

2 結(jié)果與分析

2.1 Cyt b基因序列特征

5個(gè)江段的野生黃顙魚mtDNA Cyt b的擴(kuò)增產(chǎn)物利用Clustal-X 1.81和MEGA 5.0軟件進(jìn)行人工比對(duì)序列，并對(duì)數(shù)據(jù)進(jìn)行適當(dāng)?shù)募羟刑幚?，確定黃顙魚Cyt b基因片段全長為1109 bp。由MEGA5.0軟件對(duì)1109 bp的序列進(jìn)行分析可得，5個(gè)江段的190尾樣本基因片段中，T、C、A和G堿基平均含量分別為27.4%、31.0%、27.7%、13.8%，可見C堿基含量相對(duì)較高，G堿基含量相對(duì)較低，其中A+T的含量 (55.1%)顯著高于G+C含量(44.8%)，具體可見表2?？偣矙z測(cè)到43種單倍型，序列中東江 (DJ)、左江 (ZJ)和右江 (YJ)的樣品的變異位點(diǎn)數(shù)目均大于100個(gè)，分別為230(約占20.7%)、193(約占17.4%)、352(約占31.7%)個(gè)，都勻 (GZ)與桂林 (GL)江段的樣品的變異位點(diǎn)個(gè)數(shù)分別為8個(gè)和14個(gè);5個(gè)江段的樣品檢測(cè)到的簡約信息位點(diǎn)分別是東江為192個(gè)、左江142個(gè)，右江307個(gè)，都勻4個(gè)，桂林5個(gè) (表3)。經(jīng)檢測(cè)群體內(nèi)的轉(zhuǎn)換/顛換 (Ts/Tv)之比分別為:DJ為4.07，ZJ為3.91，YJ為2.44，GZ為8.25，GL為5.51。由轉(zhuǎn)換、顛換對(duì)遺傳距離圖 (圖1)可知:轉(zhuǎn)換和顛換數(shù)沒有達(dá)到平臺(tái)效應(yīng)，即沒有達(dá)到飽和狀態(tài)，說明供統(tǒng)計(jì)分析的數(shù)據(jù)準(zhǔn)確可靠。此外，都勻和桂林群體轉(zhuǎn)換位點(diǎn)均為1個(gè)，沒有顛換位點(diǎn)，轉(zhuǎn)換遠(yuǎn)高于顛換，表明此片段沒有飽和，適合進(jìn)行遺傳變異分析。

圖1 Cyt b基因轉(zhuǎn)換、顛換與遺傳距離之間關(guān)系的散點(diǎn)圖Fig.1 Transition and transversion substitutions vs.genetic distance of mitochondrial Cyt b gene

表2 黃顙魚線粒體Cyt b基因堿基組成Table 2 Base composition of mitochondrial Cyt b gene of Pelteobagrus fulvidraco %

表3 黃顙魚Cyt b基因序列保守位點(diǎn)、變異位點(diǎn)、簡約信息位點(diǎn)Table 3 The conserved sites、variable sites、parsimony-informative sites of mitochondrial Cyt b gene of Pelteobagrus fulvidraco 個(gè)

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

利用MEGA5.0軟件分析測(cè)序獲得的190尾樣品的mtDNA Cyt b基因序列，得到5個(gè)江段群體間的遺傳距離 (表4)。從種間遺傳距離來看，東江采樣的黃顙魚群體和右江采樣的黃顙魚群體的種群間的遺傳距離最遠(yuǎn)為0.084，遺傳距離最近的為都勻和桂林群體 (0.001)。經(jīng)比較分析，種群內(nèi)部遺傳距離從大到小依次為:YJ(0.095)＞ZJ(0.068)＞DJ(0.053)＞GZ(0.001)=GL(0.001)。

從各個(gè)地理群體的核苷酸多樣性 (Pi)來看，核苷酸多樣性均在0.00072～0.08445之間，因此可以得出，黃顙魚各群體的線粒體Cyt b基因的核苷酸變異均處于較低水平，而從平均核苷酸差異數(shù)比較中，右江群體 (YJ)最高，其次依次為左江群體 (ZJ)、東江群體 (DJ)和桂林群體 (GL)，而都勻群體 (GZ)最低。

本實(shí)驗(yàn)結(jié)果顯示，5個(gè)群體間的Fst及Nm值見表5，所示結(jié)果左江 (ZJ)與右江 (YJ)群體間的Fst值為0.02604，數(shù)值在0～0.05之間，表明2個(gè)黃顙魚群體間存在中等偏低的遺傳分化，不存在顯著的系統(tǒng)地理格局。本研究中群體間的遺傳分化系數(shù)在0.02604～0.41430之間，都勻 (GZ)跟廣西崇左 (ZJ)的遺傳分化程度最大，右江和左江的遺傳分化程度最小。在本實(shí)驗(yàn)中，5個(gè)群體間的基因流系數(shù)Nm在范圍0.35342～9.35061之內(nèi)，其中，GZ和ZJ的Nm值最小 (0.35342)，而ZJ和YJ的Nm值最大 (9.35061)。

分別將東江、桂林、都勻、左江和右江5個(gè)江段群體，左右江流域與東江流域二群體，左右江流域兩群體作為一個(gè)組，運(yùn)用Arlequin軟件進(jìn)行群體間的分子變異等級(jí)分析 (AMOVA)(表6)，結(jié)果表明:5個(gè)群體總的遺傳分化系數(shù)Fst=0.24774(P=0.00±0.000＊＊＜0.05)，差異極顯著。

表4 黃顙魚的Cyt b基因序列群體間遺傳距離Table 4 The interpopulation genetic distances of mitochondrial Cyt b gene of Pelteobagrus fulvidraco

表5 黃顙魚群體間的Fst(對(duì)角線下方)和Nm(對(duì)角線上方)Table 5 The interpopulation Fst(on triangle)and Nm(lower triangle)of Pelteobagrus fulvidraco

2.3 聚類分析

利用MEGA5.0軟件分析5個(gè)江段的黃顙魚個(gè)體的 mtDNA Cyt b基因序列，并用最大似然法(Maximum Likelihood)進(jìn)行聚類分析，得到分子系統(tǒng)發(fā)生樹 (圖2)。圖中各枝上數(shù)字為1000次Bootstra P自展法檢驗(yàn)統(tǒng)計(jì)分析后對(duì)該枝的支持率，數(shù)值越大，支持率越高，則表明該樹的可信度越高。一般認(rèn)為，數(shù)值＞70%，即表明該枝的可靠性比較大。5個(gè)黃顙魚采樣江段中左江和右江群體先聚為一支，桂林和都勻聚為一支，采自于東江的3號(hào)個(gè)體單獨(dú)聚為一支。

表6 黃顙魚5個(gè)群體間遺傳變異的分子變異等級(jí)分析 (AMOVA)1)Table 6 Analysis of molecular variance(AMOVA)of five populations of Pelteobagrus fulvidraco

圖2 珠江水系野生黃顙魚Cyt b基因序列的ML分子系統(tǒng)樹Fig.2 The molecular phylogenetic tree of wild Pelteobagrus fulvidraco from the Pearl River by ML method based on Cyt b

2.4 中性檢測(cè)結(jié)果

通過應(yīng)用DNASP5.0軟件分別對(duì)5個(gè)江段的野生黃顙魚Cyt b基因序列進(jìn)行中性檢測(cè)，得到Fu's Fs值和Tajima's D值 (結(jié)果見表7)，從而揭示黃顙魚是否經(jīng)歷過種群擴(kuò)張。結(jié)果表明，桂林和都勻2個(gè)群體的 Fu's Fs值是小于0且 P值不顯著，Tajima's D檢驗(yàn)D值小于0且P值不顯著，左江和右江群體Fu's Fs值是大于0且Tajima's D檢驗(yàn)D值大于0，東江群體Fu's Fs值是大于0，Tajima's D檢驗(yàn)D值小于0但P值不顯著。

表7 黃顙魚Cyt b基因Fu'S Fs和Tajima's D檢驗(yàn)1)Table 7 The Test of Fu'S Fs and Tajima's D of mitochondrial Cyt b gene of Pelteobagrus fulvidraco

3 討論

遺傳多樣性 (Genetic diversity)是生物多樣性的核心，是生物種內(nèi)和種間遺傳變異的總和，是生物進(jìn)化和物種分化的基礎(chǔ)［13］。

3.1 黃顙魚的遺傳多樣性分析

本研究中5個(gè)江段群體的野生黃顙魚Cyt b基因序列的堿基組成分別為 A(27.7%)、T(27.4%)、C(31.0%)、G(13.8%)，其中G的含量顯著低于其他堿基的含量，表現(xiàn)出明顯的反G偏倚，顯示出細(xì)胞色素b基因的共同特征［14］，這也是線粒體DNA的一個(gè)特點(diǎn)［15］。與多數(shù)硬骨魚類相似［5］，堿基組成差異不大，分布比較平均。

右江 (YJ)群體的遺傳多樣性最高 (Hd=0.92038)，屬于高單倍型。將珠江水系的野生黃顙魚5個(gè)江段的群體作為一個(gè)整體，其Hd和Pi值分別為0.84857和0.04817，呈高單倍型、低核苷酸多樣性的分布模式，表明其遺傳多樣性相對(duì)中等。與鐘立強(qiáng)等［16］在研究長江中下游5個(gè)湖泊黃顙魚種群線粒體細(xì)胞色素b基因的遺傳變異分析中所得結(jié)果所示基本一致，黃顙魚的平均單倍型多樣性為0.945，核苷酸多樣性為0.00419，其結(jié)果均呈現(xiàn)黃顙魚的遺傳多樣性表現(xiàn)為中等水平。物種的遺傳多樣性的高低與其適應(yīng)能力、生存能力和進(jìn)化潛力密切相關(guān)，遺傳變異是有機(jī)體適應(yīng)環(huán)境變化的必要條件［17－18］。因此，加強(qiáng)對(duì)現(xiàn)有資源科學(xué)管理和保護(hù)、恢復(fù)有效種群大小、豐富遺傳多樣性是保護(hù)野生自然資源的基礎(chǔ)［19］。

3.2 黃顙魚的遺傳變異分析

遺傳多樣性是指生物種內(nèi)和種間的遺傳變異度，是生物適應(yīng)環(huán)境與進(jìn)化的基礎(chǔ)［20］。物種對(duì)生存環(huán)境的適應(yīng)能力以及其進(jìn)化的潛力依靠種群內(nèi)遺傳變異的大小，也依賴于遺傳變異的種群結(jié)構(gòu)［21］。單倍型之間的遺傳距離是衡量一個(gè)物種或群體的mtDNA變異程度的重要指標(biāo)［22］。研究表明，遺傳距離是衡量群體間遺傳變異程度的可靠參數(shù)，遺傳距離越大表明群體間親緣關(guān)系越遠(yuǎn)。利用軟件MEGA5.0計(jì)算5個(gè)江段的野生黃顙魚各群體間和群體內(nèi)的遺傳距離，其結(jié)果均在0.001～0.08之間，遺傳距離較小，則均顯示地理群體間有較近的親緣關(guān)系。

群體遺傳學(xué)認(rèn)為，F(xiàn)st值可以表示群體間的分化程度。根據(jù)Wright［23］提出的標(biāo)準(zhǔn)，一般指當(dāng)0＜Fst＜0.05時(shí)表示分化較弱或者表示群體間無分化，當(dāng)0.05＜Fst＜0.15時(shí)表示遺傳分化中等，如0.15＜Fst＜0.25時(shí)，則表示遺傳分化較大，而超過0.25表示遺傳分化極大或表明群體呈高度分化。本研究中，5個(gè)野生黃顙魚地理群體間的遺傳分化系數(shù)在0.02604～0.41430之間，都勻跟左江江段群體的遺傳分化程度最大，左江和右江江段群體的遺傳分化程度最小。左江和右江群體間的Fst值為0.02604，數(shù)值在0～0.05之間，表明2個(gè)江段野生黃顙魚種群間存在中等偏低的遺傳分化，不存在顯著的地理隔離現(xiàn)象。而都勻與左江的地理位置相對(duì)較近，遺傳分化程度卻最大，原因可能是由于都勻和左江流域之間有紅水河流域隔開，加之左江水電站的建立更加大了2個(gè)江段群體間野生黃顙魚的遺傳分化的原因。

基因流Nm通常被用來評(píng)價(jià)不同地理群體間基因交流的程度。Nm值的大小與2個(gè)群體間基因交流的程度呈正相關(guān)，即當(dāng)Nm的值越大，表示2個(gè)群體間的基因交流頻繁，反之亦然。Wright認(rèn)為:當(dāng)Nm遠(yuǎn)遠(yuǎn)小于1，種群會(huì)被強(qiáng)烈的分化，表明2個(gè)群體間的基因交流有限，存在一定的地理隔離;當(dāng)種群基因流系數(shù)Nm＞1時(shí)，種群間存在一定的基因流動(dòng)，有較高的基因交流水平;如果Nm＞4，表明是2個(gè)群體是隨機(jī)交配，基因交流頻繁［24］。本實(shí)驗(yàn)中，5個(gè)江段間黃顙魚的基因流系數(shù)Nm在范圍0.35342～9.35061之內(nèi)，其中，都勻和左江江段群體間的Nm值最小，說明存在一定地理隔離，基因交流有限，而左江和右江群體間的Nm值最大，遠(yuǎn)大于4，說明2個(gè)群體間基因交流頻繁，而且左右江的距離也是相對(duì)來說最近的，也驗(yàn)證了地理位置是原因之一。黃顙魚是一種生態(tài)適應(yīng)性很強(qiáng)的魚類，不同水系的居群借助于水系的連通而發(fā)生基因交流的可能性很大，5個(gè)江段的黃顙魚生境差異不明顯，在一定程度上為群體間的基因交流提供了有利條件，這可能是造成不同江段珠江流域間野生黃顙魚群體沒有明顯遺傳分化的原因之一。

3.3 黃顙魚的種群動(dòng)態(tài)分析

本研究應(yīng)用中性檢測(cè)來判斷種群是否經(jīng)歷過擴(kuò)張。應(yīng)用Tajima's D與Fu's Fs值中性檢驗(yàn)推測(cè)種群歷史時(shí)，如果Tajima's D與Fu's Fs值呈負(fù)值，且在統(tǒng)計(jì)學(xué)上達(dá)到顯著標(biāo)準(zhǔn)，則說明序列中含有比中性進(jìn)化模型更多的核苷酸位點(diǎn)變化，可能預(yù)示著被研究種群曾經(jīng)經(jīng)歷過一個(gè)擴(kuò)張的歷史［25］。

通過應(yīng)用DNASP軟件分別對(duì)5個(gè)江段的野生黃顙魚Cyt b基因序列進(jìn)行中性檢測(cè)，得到Fu's Fs值和Tajima's D值，結(jié)果表明，桂林和都勻2個(gè)江段的黃顙魚群體的Fu's Fs值是小于0且P值不顯著，Tajima's D檢驗(yàn)D值小于0且P值不顯著，左江、右江群體Fu's Fs值是大于0且Tajima's D檢驗(yàn)D值大于0，東江群體Fu's Fs值是大于0，Tajima's D檢驗(yàn)D值小于0但P值不顯著，均表明群體沒有發(fā)生過種群擴(kuò)張，種群數(shù)量相對(duì)穩(wěn)定。

物種的遺傳多樣性反映了一個(gè)物種對(duì)環(huán)境的適應(yīng)能力、生存能力和進(jìn)化潛力，其遺傳多樣性越豐富，對(duì)環(huán)境的適應(yīng)能力就越大，生存和進(jìn)化的能力也越強(qiáng)。保護(hù)珠江及其流域內(nèi)各支流和湖泊中魚類的遺傳多樣性是我國重要的漁業(yè)管理政策。根據(jù)本實(shí)驗(yàn)結(jié)果顯示，珠江流域野生黃顙魚的遺傳多樣性處于相對(duì)較低的水平，已經(jīng)達(dá)到亟需保護(hù)的狀態(tài)。所以，有必要提前采取適當(dāng)?shù)墓芾聿呗砸跃S持其野生資源的可持續(xù)性開發(fā)以及維護(hù)不同水系間遺傳多樣性和原有的生態(tài)平衡。

［1］褚新洛.中國動(dòng)物志，硬骨魚綱，鲇形目［M］.北京:科學(xué)出版社，1999:152－156.

［2］宋平，潘云峰，向筑，等.黃顙魚RAPD標(biāo)記及其遺傳多樣性的初步分析［J］.武漢大學(xué)學(xué)報(bào):理學(xué)版，2001，47(2):233－237.

［3］代應(yīng)貴，陳毅峰.清水江的魚類區(qū)系及生態(tài)類型［J］.生態(tài)學(xué)雜志，2007，26(5):682－687.

［4］焦樹林，高全洲，劉昆.珠江流域西江、北江河流溶解無機(jī)碳及其穩(wěn)定同位素組成特征［J］.中山大學(xué)學(xué)報(bào):自然科學(xué)版，2009，48(2):99－105.

［5］楊慧榮，趙會(huì)宏，蒙子寧，等.赤眼鱒線粒體D-loop和Cyt b基因序列的對(duì)比分析［J］.中山大學(xué)學(xué)報(bào):自然科學(xué)版，2012，51(5):100－106.

［6］彭作剛，張耀光，何舜平，等.從細(xì)胞色素b基因序列變異分析中國鲇形目魚類的系統(tǒng)發(fā)育［J］.遺傳學(xué)報(bào)，2005，32(2):145－154.

［7］MEYER A.Evolution of mitochondrial DNA in fishes［J］.Biochemistry and Molecular Biology of Fishes，1993，2:1－38.

［8］TSIGENOPOULOS C S，BERREBI P.Molecular Phylogeny of North Mediterranean freshwater barbs inferred from cytochrome b sequences:Biogeographic and systematic Implications［J］.Molecular Phylogenetics and Evolution，2000，14(2):165－179.

［9］JOHNSON J B，JORDAN S.Phylogenetic divergence in leatherside chub(Gila copei)inferred from mitochondrial cytochrome b sequences［J］.Molecular Ecology，2000，9(8):1029－1035.

［10］LYDEARD C，ROE K J.The phylogenetic utility of the mitochondrial cytochrome b gene for inferring relationships among actinopterygian fishes［J］.Molecular Systematics of Fishes，1997，285－303.

［11］TANG Q Y，LIU H Z，MAYDEN R，et al.Comparison of evolutionary rates in the mitochondrial DNA cytochrome b gene and control region and their implications for phylogeny of the Cobitoidea(Teleostei:Cypriniformes)［J］.Molecular Phylogenetics and Evolution，2006，39:347－357.

［12］THOMPSON J D，GIBSON T J，PLEWNIAK F，et al.The CLUSTAL_X windows interface:flexible strategies for multiple sequence alignment aided by quality analysis tools［J］.Nucleic Acids Research，1997，25(24):4876－4882.

［13］施立明.遺傳多樣性及其保存［J］.生物科學(xué)信息，1990，2(4):159－164.

［14］ROGERS A R，HARPENDING H.Population growth makes waves in the distribution of pairwise genetic differences［J］.Molecular Biology and Evolution，1992，9(3):552－569.

［15］WOLSTENHOLME D R.Animal mitochondrial DNA:structure and evolution［J］.Int Rev Cytol，1992，141:173－216.

［16］鐘立強(qiáng)，王明華，邊文冀，等.長江中下游5個(gè)湖泊黃顙魚(Pelteobagrus fulvidraco)種群線粒體細(xì)胞色素b基因的遺傳變異分析［J］.湖泊科學(xué)，2013，25(2):302－308.

［17］CONRAD M.Adaptability:the significance of variability from molecular to ecosystem［M］.New York:Plenum Press，1983.

［18］FERAL J P.How useful are the genetic markers in attempts to understand and manage marine biodiversity［J］.Journal of Experimental Marine Biology and Ecology，2002，268:121－145.

［19］蒙子寧，楊麗萍，吳豐，等.斜帶石斑魚、赤點(diǎn)石斑魚RAPD和線粒體Cyt b基因序列變異分析［J］.中山大學(xué)學(xué)報(bào):自然科學(xué)版，2007，46(1):75－80.

［20］李祥龍，張亞平，陳勝偶，等.山羊多態(tài)性及其起源分化研究［J］.畜牧獸醫(yī)學(xué)報(bào)，1999，30(4):313－319.

［21］GRANT P R，GRANT B R.Hybridization of bird species.［J］.Science，1996，256:193－197.

［22］BILLINGTON N，HEBERT P D N.Mitochodrial DNA diversity in fieshes and its Implications for introductions［J］.Con Fish Aquant Sci，1991，48:80－94.

［23］WRIGHT S.The interpretation of population structure by F-statistics with special regard to systems of mating［J］.Evolution，1965，19:395－420.

［24］SHAKLEE J B，TAMARU C S，WAPLES R W.Speciation and evolution of marine fishes studied by electrophoresis analysis of proteins［J］.Pac Sci，1982，36:141－157.

［25］TAJIMA F.Statistical method for testing the neutral mutation hypothesis by DNA polymorphism［J］.Genetics，1989，123:585－595.