赤水河半鰒的遺傳多樣性和種群歷史動(dòng)態(tài)分析

李文靜 王環(huán)珊 劉煥章 曹文宣 王小明 于 江 朱 昕

(1. 中國(guó)科學(xué)院水生生物研究所水生生物多樣性與保護(hù)重點(diǎn)實(shí)驗(yàn)室, 武漢 430072; 2. 中國(guó)科學(xué)院大學(xué), 北京 100049;3. 中國(guó)三峽建設(shè)管理有限公司, 成都 610041)

赤水河有著美酒河、生態(tài)河、英雄河的美稱(chēng),是長(zhǎng)江上游唯一一條干流未建水壩, 且與長(zhǎng)江保持連通的大型一級(jí)支流[1]。赤水河作為長(zhǎng)江上游珍稀、特有魚(yú)類(lèi)國(guó)家級(jí)自然保護(hù)區(qū)的重要組成部分,是一個(gè)與保護(hù)區(qū)既相聯(lián)系又相對(duì)獨(dú)立的系統(tǒng)。根據(jù)中國(guó)科學(xué)院水生生物研究所的監(jiān)測(cè), 半鰐在赤水河有一定的種群數(shù)量, 因此具有著重要的保護(hù)價(jià)值。

半鰐Hemiculterella sauvagei隸屬于鯉形目Cypriniformes鯉科Cyprinidae鲌亞科Cultrinae半鰐屬Hemiculterella, 地方名稱(chēng)藍(lán)片子、鰐條、藍(lán)刀皮[1,2]。半鰐是長(zhǎng)江上游特有魚(yú)類(lèi), 在長(zhǎng)江上游干流、沱江、岷江、嘉陵江、青衣江和大渡河等水系有分布[2,3]。近年來(lái), 由于水利工程建設(shè)、酷漁濫捕等人為條件的影響, 其野生資源已受到嚴(yán)重威脅。目前, 關(guān)于半鰐的研究主要涉及其地理分布[4]、年齡和生長(zhǎng)[1]、繁殖生物學(xué)[5]、寄生蟲(chóng)[6]、物種保護(hù)[7]以及線粒體全基因組序列的測(cè)定[8], 但有關(guān)半鰐的遺傳多樣性及種群歷史動(dòng)態(tài)的工作尚未得到充分開(kāi)展。遺傳多樣性不僅是生物多樣性的核心部分, 也是生命進(jìn)化和物種分化的基礎(chǔ), 更是評(píng)價(jià)自然生物資源的重要依據(jù)[9]。

線粒體DNA (Mitochondrial DNA, mtDNA)具有相對(duì)較高的進(jìn)化速率, 目前已廣泛應(yīng)用于各種動(dòng)物類(lèi)群的系統(tǒng)發(fā)育和種群遺傳多樣性研究[10]。Cytb基因是一段編碼蛋白質(zhì)的DNA序列, 已有的研究表明在鯉形目魚(yú)類(lèi)的物種及種以下階元中, 其進(jìn)化速率明顯高于控制區(qū)基因[11], 因而在鯉形目魚(yú)類(lèi)種群的遺傳多樣性和歷史動(dòng)態(tài)研究中比控制區(qū)更加合適[12—16]。

本研究選取赤水河水潦、茅臺(tái)、赤水市、二合及太平5個(gè)不同地理種群的168尾半鰐為實(shí)驗(yàn)對(duì)象, 對(duì)其Cytb基因序列進(jìn)行擴(kuò)增并測(cè)序, 期望通過(guò)研究半鰐種群的遺傳多樣性, 對(duì)赤水河半鰐資源現(xiàn)狀進(jìn)行評(píng)估, 從而為半鰐資源的保護(hù)提供理論基礎(chǔ)。

1 材料與方法

1.1 材料

半鰐樣本采自赤水河5個(gè)地點(diǎn)(水潦、茅臺(tái)鎮(zhèn)、赤水市、二合鎮(zhèn)、太平鎮(zhèn)), 采樣地點(diǎn)、樣本數(shù)量見(jiàn)圖1和表 1。從每尾新鮮樣本背鰭以下、側(cè)線以上的右側(cè)取2—3 g肌肉, 以95% 的酒精固定于離心管中帶回實(shí)驗(yàn)室, 樣本儲(chǔ)存于4℃的冰箱中。另外,在GenBank中下載鰐(Hemiculter leucisculus)(登錄號(hào)為AY089716)和蒙古鲌(Culter mongolicus)(登錄號(hào)為AP009060)的Cytb序列作為外類(lèi)群。

1.2 基因組DNA的提取以及PCR擴(kuò)增

本研究采用改進(jìn)的高鹽法提取基因組DNA[17],并采用以下引物來(lái)擴(kuò)增Cytb基因序列: L14724(5′-GACTTGAAAAACCACCGTTG-3′和H15915(5′-CTCCGATCTCCGGATTACAAGAC-3′)[18]。

聚合酶鏈?zhǔn)椒磻?yīng)(Polymerase Chain Reaction,PCR)的反應(yīng)總體積為30 μL: 10×buffer 3.0 μL,dNTP (2.5 mmol/L)1.5 μL,Taq聚合酶0.25 μL, 模板DNA 3 μL, 正反向引物各1 μL, ddH2O 20.25 μL。PCR反應(yīng)程序?yàn)? 94℃變性4min, 接著是35個(gè)循環(huán),每個(gè)循環(huán)包括: 94℃變性45s, 56℃退火45s, 72℃延伸1min, 然后72℃延伸10min, 最后10℃保存10min。

用1.0%的瓊脂糖凝膠電泳對(duì)PCR擴(kuò)增的產(chǎn)物進(jìn)行檢測(cè), 挑選其中有目的條帶, 且條帶清晰的產(chǎn)物送生物公司進(jìn)行純化及測(cè)序。

1.3 數(shù)據(jù)處理及分析

使用Clustal X[19]和Seaview[20]對(duì)測(cè)序結(jié)果進(jìn)行序列編輯、校對(duì)和排序。

圖1 赤水河流域半鰐采樣點(diǎn)(●代表采樣點(diǎn)的位置)Fig. 1 Sampling sites of Hemiculterella sauvagei in Chishui River (●indicates sampling site)

表 1 基于線粒體Cyt b基因的赤水河半鰒各地理種群遺傳多樣性Tab. 1 Genetic diversity in five populations of Hemiculterella sauvagei in the Chishui River based on mtDNA Cyt b gene

使用DnaSP 5.0軟件[21]計(jì)算單倍型數(shù)目(Num-ber of haplotypes)、多態(tài)位點(diǎn)數(shù)目(Number of polymorphic sites)、單倍型多樣性(Haplotype diversity)、核苷酸多樣性(Nucleotide diversity)、變異位點(diǎn)(Variable sites)、單一突變位點(diǎn)(Singleton sites)、簡(jiǎn)約信息位點(diǎn)(Parsimony-informative sites)等。

使用Mega 6.0軟件[22,23]計(jì)算序列的堿基組成、種群間的凈遺傳距離, 采用鄰接法(Neighborjoining, NJ)[24]以鰐和蒙古鲌為外類(lèi)群, 構(gòu)建單倍型系統(tǒng)發(fā)育樹(shù)。使用Network軟件, 以中接法(Medianjoining)構(gòu)建單倍型的網(wǎng)絡(luò)關(guān)系圖, 檢測(cè)各個(gè)單倍型間的進(jìn)化關(guān)系。

使用Arlequin 3.11軟件[25]統(tǒng)計(jì)遺傳分化水平,采用分子變異分析(Analysis of Molecular Variance,AMOVA)檢測(cè)遺傳變異來(lái)源, 估算種群遺傳結(jié)構(gòu)和地理種群遺傳變異的分布。用種群差異的精確性檢驗(yàn)(Exact tests of population differentiation)[26]來(lái)估計(jì)2個(gè)種群的分化水平指數(shù)Fst, 并采用Tajima’s D和Fu’s Fs中性檢驗(yàn)以及核苷酸不配對(duì)分布(Mismatch distribution)來(lái)檢測(cè)種群歷史動(dòng)態(tài)。通過(guò)粗糙指數(shù)(Harpending’s Raggedness index, Hri)[27]來(lái)檢測(cè)核苷酸不配對(duì)分布與種群擴(kuò)張模型下期望分布之間的擬合優(yōu)度。對(duì)于不能拒絕種群擴(kuò)張模型而且錯(cuò)配分布呈現(xiàn)光滑單峰的種群, 用公式τ=2ut估算種群擴(kuò)張時(shí)間[28], τ為T(mén)au, 是種群擴(kuò)張參數(shù), t為每個(gè)世代種群的擴(kuò)張時(shí)間, u為每個(gè)世代每條序列的變異速率; 另根據(jù)u=2μk計(jì)算u, μ為每個(gè)堿基的變異速率, k表示序列長(zhǎng)度(bp)。 擴(kuò)張時(shí)間T=t×世代時(shí)間,T值由Arelequin軟件計(jì)算得到, 突變率μ參照前人研究[29], 采用1% per Myr (1×108)。

2 結(jié)果

2.1 赤水河半鰒的序列變異和單倍型系統(tǒng)發(fā)育關(guān)系

本研究得到了168條半鰐線粒體Cyt b基因序列, 比對(duì)后得到序列長(zhǎng)度為1137 bp。序列中無(wú)堿基的插入和缺失。在168條序列中, 共有42個(gè)變異位點(diǎn), 占總位點(diǎn)數(shù)的3.69%, 其中單突變位點(diǎn)有13個(gè),占總位點(diǎn)數(shù)的1.14%, 簡(jiǎn)約信息位點(diǎn)有29個(gè), 占總位點(diǎn)數(shù)的2.55%。

168個(gè)樣本的平均堿基組成: A=27.5%, T=27.3%, C=29.1%, G=16.0%。A+T的含量為54.8%,C+G的含量為45.1%, A+T含量高于C+G的含量。G的含量最低, 堿基組成表現(xiàn)出了明顯的偏倚, 這恰好符合脊椎動(dòng)物線粒體DNA基因的共同特性。5個(gè)種群總體的單倍型多樣性(Hd)和核苷酸多樣性(Pi)分別為0.895±0.012和0.00487±0.00695, 表現(xiàn)出較高的單倍型多樣性和較低的核苷酸多樣性。

168條序列共檢測(cè)到38種單倍型(表 1)。其中,水潦種群22條序列檢測(cè)到12個(gè)單倍型, 茅臺(tái)鎮(zhèn)種群45條序列檢測(cè)到11個(gè)單倍型, 赤水市種群40條序列檢測(cè)到25個(gè)單倍型, 二合鎮(zhèn)種群39條序列檢測(cè)到9個(gè)單倍型, 太平鎮(zhèn)種群22條序列檢測(cè)到10個(gè)單倍型。在所檢測(cè)到的38個(gè)單倍型中, 單倍型Hap1、Hap2和Hap6的分布最廣, 為全部種群所共享。其中有24個(gè)單倍型為某一個(gè)種群所獨(dú)有, 其余11個(gè)單倍型則被2—4個(gè)種群共享。

圖2 基于Cyt b序列構(gòu)建的半鰐單倍型NJ樹(shù)Fig. 2 Phylogenetic tree of haplotypes by NJ analysis based on Cyt b sequences

以鰐和蒙古鲌為外類(lèi)群, 利用鄰接法(Neighbor joining, NJ)和最大似然法(Maximum Likelihood,ML)構(gòu)建單倍型系統(tǒng)發(fā)育樹(shù)。結(jié)果顯示兩種方法做出的系統(tǒng)發(fā)育樹(shù)拓?fù)浣Y(jié)構(gòu)一樣, 以NJ樹(shù)為例(圖2)。38個(gè)單倍型大致可以分為5個(gè)譜系, 其中譜系1由Hap9、Hap33、Hap37等12個(gè)單倍型組成; 譜系2由Hap15、Hap13、Hap14、Hap21、Hap24, 5個(gè)單倍型組成; 譜系3由Hap18、Hap38、Hap17、Hap8、Hap32, 5個(gè)單倍型組成; 譜系4由Hap22、Hap6、Hap35等12個(gè)單倍型組成; 譜系5則由Hap23、Hap29、Hap12、Hap34, 4個(gè)單倍型組成。可以看到, 同一地點(diǎn)的單倍型沒(méi)有聚為單系, 而是分散在不同的譜系中, 說(shuō)明還沒(méi)有出現(xiàn)種群的譜系分化。

利用Network軟件, 以中接法構(gòu)建的半鰐各單倍型進(jìn)化網(wǎng)絡(luò)圖(圖3), 結(jié)果表明: Hap1、Hap2、Hap6是5個(gè)地理種群的共享單倍型, Hap4是4個(gè)地理種群的共享單倍型。該單倍型關(guān)系圖與系統(tǒng)發(fā)育樹(shù)結(jié)果一致, 各單倍型相互散布在不同的地理種群中, 未形成明顯的系統(tǒng)地理格局。

2.2 赤水河半鰒種群遺傳結(jié)構(gòu)和種群分化

種群遺傳多樣性和遺傳距離從遺傳多樣性指數(shù)來(lái)看, 在這5個(gè)種群中, 赤水市的單倍型多樣性最高, 為0.922±0.032, 茅臺(tái)鎮(zhèn)的單倍型多樣性最低, 為0.842±0.032; 水潦的核苷酸多樣性最高, 為0.00495±0.00434, 茅臺(tái)鎮(zhèn)的核苷酸多樣性最低, 為0.00386±0.00342 (表 1)。

種群遺傳結(jié)構(gòu)和分化運(yùn)用Arlequin 3.11軟件進(jìn)行種群間遺傳分化(Fst)分析(表 2), 結(jié)果表明:赤水市種群與茅臺(tái)鎮(zhèn)種群、二合鎮(zhèn)種群、水潦種群之間存在高度分化(Fst＞0.15), 且差異很顯著(P=0.0000,P＜0.01), 赤水市種群與太平鎮(zhèn)種群之間有中度遺傳分化(Fst=0.11692, 0.05＜Fst＜0.15), 差異也很顯著(P=0.0000,P＜0.01)。除赤水市種群外, 其他4個(gè)種群之間均未出現(xiàn)分化(Fst＜0.05)。

采用Mega 6.0 軟件計(jì)算5個(gè)種群的遺傳距離,5個(gè)種群之間的平均凈遺傳距離范圍為0.004—0.006(表 3): 水潦種群和茅臺(tái)鎮(zhèn)種群之間的遺傳距離為0.004, 水潦種群和二合鎮(zhèn)種群之間的遺傳距離為0.005, 水潦種群和太平鎮(zhèn)種群之間的遺傳距離為0.005, 水潦種群和赤水市種群之間的遺傳距離為0.006, 茅臺(tái)鎮(zhèn)種群和二合鎮(zhèn)種群之間的遺傳距離為0.004, 茅臺(tái)鎮(zhèn)種群和太平鎮(zhèn)種群之間的遺傳距離為0.004, 茅臺(tái)鎮(zhèn)種群與赤水市種群之間的遺傳距離為0.006, 二合鎮(zhèn)種群和太平鎮(zhèn)種群之間的遺傳距離為0.005, 二合鎮(zhèn)種群與赤水市種群之間的遺傳距離為0.006, 太平鎮(zhèn)種群和赤水市種群之間的距離為0.005。

將5個(gè)地理種群兩組, 水潦、茅臺(tái)鎮(zhèn)、二合鎮(zhèn)和太平鎮(zhèn)為一組, 赤水市一組, 然后采用Arlequin 3.11軟件進(jìn)行分子方差分析(AMOVA)(表 4), 結(jié)果顯示: 赤水河半鰐種群內(nèi)存在較高的遺傳變異(80.95%), 而組間變異占種群遺傳變異的19.82%,組內(nèi)種群間的變異占-0.77%, 由此可見(jiàn), 赤水河半鰐的遺傳變異主要來(lái)源于種群內(nèi)部, 且在P=0.01水平上顯著, 5個(gè)地理種群之間存在基因交流。

2.3 種群擴(kuò)張及分化時(shí)間

使用Arlequin3.11軟件對(duì)5個(gè)種群所有的168個(gè)半鰐個(gè)體進(jìn)行Tajima’sD檢驗(yàn)和Fu’sFs檢驗(yàn), 以及種群擴(kuò)張時(shí)間的估算(表 5)。結(jié)果表明: Tajima’sD檢驗(yàn)只有赤水市種群為負(fù)值, 但并未達(dá)到顯著性水平(P=0.12100,P＞0.05)。在Fu’sFs檢驗(yàn)中, 水潦種群、赤水市種群和太平鎮(zhèn)種群的Fs值均為負(fù), 但只有赤水市種群達(dá)到顯著性水平(P=0.00000,P＜0.01)。以上結(jié)果表明, 水潦、太平鎮(zhèn)和赤水市種群經(jīng)歷了種群擴(kuò)張。采用1% per Myr的突變速率, 本研究推算出半鰐水潦種群、太平鎮(zhèn)種群以及赤水市種群的大約擴(kuò)張時(shí)間分別為0.43、0.40和0.37百萬(wàn)年前。

圖3 半鰐各單倍型間的進(jìn)化網(wǎng)絡(luò)圖Fig. 3 Statistical parsimony network based on cyt b haplotype frequencies of Hemiculterella sauvagei

表 2 基于Cyt b序列單倍型頻率的半鰒種群成對(duì)Fst值(對(duì)角線下)和校正P值(對(duì)角線上)Tab. 2 F-Statistics (below diagonal)and Corrected P value (above diagonal)from haplotype frequencies of Hemiculterella sauvagei based on Cyt b gene

表 3 半鰒各地理種群間的遺傳距離Tab. 3 The distances among 5 Hemiculterella sauvagei populations

表 4 基于線粒體cyt b基因?qū)Τ嗨影膂v不同地理種群的AMOVA分析Tab. 4 AMOVA analysis for Hemiculterella sauvagei among different populations based on cyt b in the Chishui River

表 5 基于Cyt b基因的半鰒種群中性檢驗(yàn)及種群擴(kuò)張時(shí)間估算Tab. 5 Statistical tests for neutrality and time of expansion of Hemiculterella sauvagei based on Cyt b gene

3 討論

3.1 半鰒種群遺傳多樣性

生物多樣性包括遺傳多樣性、物種多樣性和生態(tài)系統(tǒng)多樣性, 其中遺傳多樣性是物種多樣性和生態(tài)系統(tǒng)多樣性的基礎(chǔ), 也是生物多樣性的核心和重要組成部分。核苷酸多樣性(Pi)和單倍型多樣性(Hd)是評(píng)價(jià)種群遺傳多樣性的重要指標(biāo), 其值越大,則種群的遺傳多樣性就越高。本研究中赤水河半鰐的Cytb基因的核苷酸多樣性的范圍為0.386%—0.495%, 與銅魚(yú)(Coreius heterodon)(0.25%)、草魚(yú)(Ctenopharyngodon idellus)(0.092%)[30]、刀鱭(Coilia macrognathos)(0.185%—0.373%)[31]等種類(lèi)相比,核苷酸多樣性相對(duì)較高。整個(gè)半鰐種群的單倍型多樣性均大于0.5, 核苷酸多樣性均小于0.5%, 故表現(xiàn)出高的Hd值和低的Pi值。

3.2 半鰒地理種群的遺傳分化

同一種魚(yú)棲息在不同的環(huán)境中, 會(huì)造成魚(yú)類(lèi)基因流的限制或中斷。遺傳分化系數(shù)(Fst)是反映亞群間遺傳分化的重要指標(biāo)。Wright[32]認(rèn)為: 若種群Fst為0—0.05, 則表明其各亞群間不存在分化; 若Fst為0.05—0.15, 則表明其各亞群間存在中度分化;若Fst為0.15—0.25, 則為高度分化。在本研究中的5個(gè)地理種群中, 赤水市種群和茅臺(tái)鎮(zhèn)種群、二合鎮(zhèn)種群、水潦種群之間存在著高度的分化, 而與太平鎮(zhèn)種群之間存在著中度分化。AMOVA結(jié)果表明, 種群內(nèi)變異是遺傳變異的主要來(lái)源, 并且存在顯著差異(P＜0.01), 因此, 赤水市與其他4個(gè)地理種群之間存在著分化, 缺少基因交流。我們認(rèn)為造成分化的可能原因是赤水市在地理位置上處于赤水河較下游地區(qū), 而其他4個(gè)地方處于上游或中游, 因此他們之間可能存在著地理隔離, 而地理隔離又會(huì)造成生殖隔離, 從而造成不同地理種群的半鰐之間發(fā)生分化。

3.3 半鰒種群的歷史動(dòng)態(tài)

Grant和Bowen[33]提出了一個(gè)較簡(jiǎn)單的模式, 利用單倍型多樣性(Hd)和核苷酸多樣性(Pi)來(lái)估計(jì)種群的進(jìn)化歷史, 當(dāng)Hd≥0.5,Pi＜0.5%時(shí), 是受瓶頸效應(yīng)后種群數(shù)量的迅速擴(kuò)張導(dǎo)致; 當(dāng)Hd≥0.5,Pi≥0.5%時(shí), 表示種群穩(wěn)定, 具有比較悠久的進(jìn)化歷史;當(dāng)Hd＜0.5,Pi≥0.5%時(shí), 種群經(jīng)歷了輕微的瓶頸效應(yīng), 幾乎沒(méi)有影響到核苷酸變異; 當(dāng)Hd＜0.5,Pi＜0.5%時(shí), 表明種群近期經(jīng)歷了瓶頸效應(yīng)。本研究中的5個(gè)地理種群的單倍型多樣性都較高(Hd≥0.5), 而核苷酸多樣性較低(Pi＜0.5%), 由于核苷酸多樣性所需時(shí)間比積累單倍型所需的時(shí)間漫長(zhǎng), 這說(shuō)明本研究中的五個(gè)地理種群的半鰐種群是從一個(gè)有效種群數(shù)量較小的種群, 經(jīng)快速擴(kuò)張而來(lái)的, 但是仍然沒(méi)有達(dá)到積累核苷酸變異所需要的時(shí)間。整個(gè)半鰐種群的單倍型多樣性較高, 而核苷酸多樣性較低,這說(shuō)明整個(gè)半鰐種群都是從較小的種群迅速擴(kuò)張而來(lái)的[33—35]。

推測(cè)種群經(jīng)歷過(guò)的歷史, 一般要用到Tajima’sD檢驗(yàn)和Fu’sFs檢驗(yàn), Tajima’sD和Fu’sFs呈負(fù)值,并且在統(tǒng)計(jì)學(xué)上有較顯著的標(biāo)準(zhǔn), 則說(shuō)明序列中含有比中性進(jìn)化模型更多的核苷酸位點(diǎn)變化, 可能預(yù)示著被研究種群經(jīng)歷過(guò)一個(gè)擴(kuò)張的歷史[36,37]。本研究中只有赤水市種群的Tajima’sD值呈負(fù)值, 但不顯著偏離中性檢驗(yàn)(P＞0.5), 水潦、太平鎮(zhèn)、赤水市種群的Fu’sFs值均呈負(fù)值, 但只有赤水市種群與中性進(jìn)化的差異顯著, 說(shuō)明水潦、太平鎮(zhèn)種群以及赤水市種群在歷史上出現(xiàn)擴(kuò)張。

3.4 關(guān)于赤水河半鰒的保護(hù)

半鰐是長(zhǎng)江上游的特有魚(yú)類(lèi)之一, 在赤水河有分布。有文獻(xiàn)記載, 半鰐分布于長(zhǎng)江上游的干流,但主要集中在四川的西部[2,3]。但從近幾年的漁獲物調(diào)查中發(fā)現(xiàn), 長(zhǎng)江上游已很少有半鰐的蹤跡。近些年赤水河的漁獲物結(jié)果表明, 半鰐在2014—2015年的漁獲物比重相比2014年以前有所下降, 原因可能為酷漁濫捕、水質(zhì)污染、天然繁殖場(chǎng)遭到破壞以及其他環(huán)境條件改變等。

半鰐作為赤水河水生生態(tài)系統(tǒng)的重要組成部分, 對(duì)維系赤水河以及長(zhǎng)江的生態(tài)系統(tǒng)結(jié)構(gòu)完整和功能的正常發(fā)揮具有重要的意義。目前赤水河的半鰐資源減少, 對(duì)其進(jìn)行保護(hù)刻不容緩。本研究結(jié)果顯示, 赤水市種群與其他4個(gè)地理種群存在分化,所以在保護(hù)工作中, 應(yīng)將他們作為兩個(gè)不同的管理單元, 采取不同的保護(hù)措施進(jìn)行保護(hù), 應(yīng)對(duì)赤水市種群采取重點(diǎn)保護(hù)。另外, 要在半鰐產(chǎn)卵期禁止捕魚(yú), 并對(duì)捕撈工具或網(wǎng)目大小做出明確的規(guī)定, 保證半鰐的生殖繁衍和資源恢復(fù)。

致謝：

感謝中國(guó)科學(xué)院水生生物研究所魚(yú)類(lèi)生態(tài)學(xué)與資源保護(hù)學(xué)科組工作人員在采樣上給予幫助, 感謝唐瓊英副研究員、鞏政博士對(duì)文章修改提出的寶貴意見(jiàn)！

[1]Wang J, Wang M R, Dan S G,et al. Age and growth ofHemiculterella sauvageifrom the Chishui River in the upper reaches of the Yangtze River [J].Sichuan Journal of Zoology, 2012, 31(5): 713—719 [王俊, 王美榮, 但勝?lài)?guó), 等. 赤水河半鰐年齡與生長(zhǎng). 四川動(dòng)物, 2012, 31(5):713—719]

[2]Ding R H. The Fishes of Sichuan [M]. Chengdu: Sichuan Publishing House of Science and Technology. 1994,210—213 [丁瑞華. 四川魚(yú)類(lèi)志. 成都: 四川科學(xué)技術(shù)出版社. 1994, 210—213]

[3]Chen Y Y. Fauna Sinica (Ostechthyes. Cypriniforms) Ⅱ[M]. Beijing: Science Press. 1998, 375—389 [陳宜瑜. 中國(guó)動(dòng)物志·硬骨魚(yú)綱鯉形目(中卷). 北京: 科學(xué)出版社.1998, 375—389]

[4]Wu J M, Wang J W, He Y F,et al. Fish assemblage structure in the Chishui River, a protected tributary of the Yangtze River [J].Knowledge and Management of Aquatic Ecosystem, 2011, 400(11): 1—14

[5]Wang J, Liu F, Zhang X,et al. Reproductive biology of Chinese minnowHemiculterella sauvageiWarpachowski,1888 in the Chishui River, China [J].Journal of Applied Ichthyology, 2014, 30(2): 314—321

[6]Liao X H, Liang Z X. Distribution of ligulid tapeworms in China [J].Journal of Parasitology, 1987, 73(1):36—48

[7]Park Y S, Chang J B, Lek S,et al. Conservation strategies for endemic fish species threatened by Three Gorges Dam[J].Conservation Biology, 2003, 17(6): 1748—1758

[8]He B, Chen Y, Liu Y,et al. The complete mitochondrial genome ofHemiculterella sauvagei(Teleostei, Cyprinidae,Hemiculterella) [J].Mitochondrial DNA Part A,2016, 27(5): 3322—3324

[9]Yin S W, Huang H, Zhang B,et al. Advancements in research on genetic diversity in groupers [J].Fisheries Science, 2005, 24(8): 46—49 [尹紹武, 黃海, 張本, 等. 石斑魚(yú)遺傳多樣性的研究進(jìn)展. 水產(chǎn)科學(xué), 2005, 24(8):46—49]

[10]Xiao W H, Zhang Y P. Genetics and evolution of mitochondrial DNA in fish [J].Acta Hydrobiologica Sinica,2000, 24(4): 384—391

[11]Tang Q Y, Liu H Z, Mayden R,et al. Comparison of evolutionary rates in the mitochondrial DNA cytochromebgene and control region and their implications for phylogeny of the Cobitoidea (Teleostei: Cypriniformes[J].Molecular Phylogenetics and Evolution, 2006, 39(2):347—357

[12]Durand J D, Tsigenopoulos C S, ünlü E,et al. Phylogeny and biogeography of the family Cyprinidae in the Middle East inferred from cytochromebDNA—evolutionary significance of this region [J].Molecular Phylogenetics and Evolution, 2002, 22(1): 91—100

[13]Perdices A, Cunha C, Coelho M M. Phylogenetic structure ofZacco platypus(Teleostei, Cyprinidae) populations on the upper and middle Chang Jiang (Yangtze)drainage inferred from cytochromebsequences [J].Molecular Phylogenetics and Evolution, 2004, 31(1):192—203

[14]Perdices A, Doadrio I. The molecular systematics and biogeography of the European cobitids based on mitochondrial DNA sequences [J].Molecular Phylogenetics and Evolution, 2001, 19(3): 468—478

[15]Peng Z, He S, Zhang Y. Mitochondrial cytochromebsequence variations and phylogeny of the East Asian bagrid catfishes [J].Progress in Natural Science, 2002, 12(6):421—425

[16]Kitagawa T, Watanabe M, Kitagawa E,et al. Phylogeography and the maternal origin of the tetraploid form of the Japanese spined loach,Cobitis biwae, revealed by mitochondrial DNA analysis [J].Ichthyological Research,2003, 50(4): 318—325

[17]Aljanabi S M, Martinez I. Universal and rapid salt-extraction of high quality genomic DNA for PCR-based techniques [J].Nucleic Acids Research, 1997, 25(22):4692—4693

[18]Xiao W, Zhang Y, Liu H. Molecular systematics of Xenocyprinae (Teleostei: Cyprinidae): taxonomy, biogeography, and coevolution of a special group restricted in East Asia [J].Molecular Phylogenetics and Evolution,2001, 18(2): 163—173

[19]Larkin M A, Blackshields G, Brown N P,et al. Clustal W and Clustal X version 2.0 [J].Bioinformatics, 2007,23(21): 2947—2948

[20]Galtier N, Gony M, Gautier C. SEAVIEW and PHYL0_WIN, two graphic tools for sequence alignment and alignment and molecular phylogeny [J].Computer Applications in the Biosciences, 1996, 12(6): 543—548

[21]Liorado P, Rozas J. DnaSP v5: A software for comprehensive analysis of DNA polymorphism data [J].Bioinformatics, 2009, 25(11): 1451—1452

[22]Tamura K, Dudley J, Nei M,et al. MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0 [J].Molecular Biology and Evolution, 2007,24(8): 1596—1599

[23]Kumar S, Dudley J, Nei M,et al. MEGA: A biologistcentric software for evolutionary analysis of DNA and protein sequences [J].Briefings in Bioinformatics, 2008,9(4): 299—306

[24]Saitou N, Nei M. The neighbor-joining method: a new method for reconstructing phylogenetic trees [J].Molecular Biology Evolution, 1987, 4(4): 406—425

[25]Excoffier L, Lischer HEL. Arlequin suite ver 3.5: A new series of programs to perform population genetics analyses under Linux and Window [J].Molecular Ecology Resources, 2010, 10(3): 564—567

[26]Raymond M, F Rousset. An exact test for population differentiation [J].Evolution, 1995, 49(6): 1280—1283

[27]Harpending H C. Signature of ancient population growth in a low-resolution mitochondrial DNA mismatch distribution [J].Human Biology, 1994, 66(4): 591—600

[28]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

[29]Zhang L, Tang Q Y, Liu H Z. Phylogeny and speciation of the eastern Asian cyprinid genus Sarcocheilichthys [J].Journal of Fish Biology, 2008, 72(5): 1122—1137

[30]Fan Q. Genetic structure and diversity ofHemiculter leucisculusfrom upper and middle Yangtze drainages revealed by mitochondrial DNA [D]. Thesis for Master of Science. Institute of Hydrobiology, Chinese Academy of Sciences, Wuhan. 2013 [范啟. 長(zhǎng)江中上游鰐(Hemiculter leucisculus)種群遺傳多樣性結(jié)構(gòu)分析和研究. 碩士學(xué)位論文, 中國(guó)科學(xué)院水生生物研究所, 武漢. 2013]

[31]Hu Y T. Morphological, genetic differentiation and community ecological process of fish in different ecological type from Poyang Lake [D]. Thesis for Doctor of Science.Institute of Hydrobiology, Chinese Academy of Sciences,Wuhan. 2012 [胡玉婷. 鄱陽(yáng)湖不同生態(tài)類(lèi)型魚(yú)類(lèi)的形態(tài)、遺傳分化及魚(yú)類(lèi)群落生態(tài)過(guò)程. 博士學(xué)位論文. 中國(guó)科學(xué)院水生生物研究所, 武漢. 2012]

[32]Wright S. Evolution and the genetics of populations:Variability within and among populations V. 4 [M]. University of Chicago Press. 1978

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

[34]Frankham R. Relationship of genetic variation to population size in wildlife [J].Conservation Biology, 1996,10(6): 1500—1508

[35]Cook R M, Sinclair A, Stefansson G. Potential collapse of North Sea cod stocks [J].Nature, 1997, 385(6616):521—522

[36]Fu Y X. Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection[J].Genetics, 1997, 147(2): 915—925

[37]Tajima F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism [J].Genetics,1989, 123(3): 585—595