賴見(jiàn)生,陳葉雨,顏 濤,林 玨,黃穎穎,趙鳳麒,陳永勝,周啟紅,劉曉歡,何 斌

(1.四川省農(nóng)業(yè)科學(xué)院水產(chǎn)研究所,成都 611731;2.宜賓市農(nóng)業(yè)農(nóng)村局,四川 宜賓 644000;3.宜賓市翠屏區(qū)農(nóng)業(yè)農(nóng)村局,四川 宜賓 644000)

【研究意義】理縣高原鰍(Triplophysalixianensis)屬鯉形目(Cypriniformes)、鰍科(Cobitidae)、條鰍亞科(Nemachilinae)、高原鰍屬。高原鰍是一種重要的淡水魚類,大多棲息于青藏高原及附近的湖泊、河流和溪流中。高源鰍屬具有高度多樣性,由大約140種組成[1],該類群的起源和演化與青藏高原的隆升有關(guān)[2]。作為青藏高原魚類區(qū)系的重要組成部分,高原鰍對(duì)高原低氧、低壓和寒冷的環(huán)境有獨(dú)特的形態(tài)和生理適應(yīng)。高原鰍對(duì)當(dāng)?shù)丨h(huán)境具有強(qiáng)烈適應(yīng)性,表明它們的長(zhǎng)期進(jìn)化與高原生境有關(guān)。因此,高原鰍是研究歷史氣候和地質(zhì)事件對(duì)該屬生物多樣性模式影響的理想物種。在中國(guó)的高原鰍種類中,目前在四川鑒定的有23種[3]。2008年何春林等[3]在四川省阿壩州雜谷瑙河流域采集到高原鰍標(biāo)本,經(jīng)形態(tài)學(xué)鑒定為高原鰍新種,并命名為理縣高原鰍。目前,關(guān)于理縣高原鰍的報(bào)道較少,僅見(jiàn)于其作為新種的描述及線粒體DNA標(biāo)記的開(kāi)發(fā)[3-4]。有關(guān)長(zhǎng)江流域理縣高原鰍種群結(jié)構(gòu)和種群歷史動(dòng)態(tài)的信息仍然處于空白。因此,對(duì)理縣高原鰍不同地理種群的遺傳多樣性和群體結(jié)構(gòu)進(jìn)行評(píng)價(jià),對(duì)進(jìn)一步利用理縣高原鰍的遺傳種質(zhì)資源和野生群體保護(hù)具有重要意義?！厩叭搜芯窟M(jìn)展】線粒體DNA(Mitochondrial DNA,mtDNA)具有進(jìn)化速度快、缺少重組、分子小、易于測(cè)序及分析等特點(diǎn),目前已成為魚類物種鑒定、進(jìn)化遺傳學(xué)、分子生態(tài)學(xué)、群體遺傳結(jié)構(gòu)與系統(tǒng)地理學(xué)等多種研究領(lǐng)域的重要標(biāo)記[5-9]?？刂茀^(qū)(Control region,D-loop)作為線粒體DNA上最重要的一段非編碼區(qū),遺傳變異程度高,堿基替換率比線粒體DNA上的其他部分高5～10倍,是適合種群水平遺傳多樣性研究特別是遺傳多樣性較低的群體遺傳學(xué)研究的最優(yōu)標(biāo)記之一[10]。研究表明,理縣高原鰍的分布非常局限,僅在長(zhǎng)江上游岷江流域雜谷腦河發(fā)現(xiàn)[3-4]?！颈狙芯壳腥朦c(diǎn)】通過(guò)對(duì)大渡河流域魚類資源調(diào)查發(fā)現(xiàn),在大渡河4個(gè)不同地段,包括金川縣沙耳鄉(xiāng)、石棉縣安順彝族鄉(xiāng)、瀘定縣冷磧鎮(zhèn)和丹巴縣格宗鎮(zhèn)均發(fā)現(xiàn)有理縣高原鰍分布,這也是首次發(fā)現(xiàn)其出現(xiàn)于雜谷腦河以外的流域?！緮M解決的關(guān)鍵問(wèn)題】選取該4個(gè)不同地理種群的理縣高原鰍群體為研究對(duì)象,分析其體長(zhǎng)與體重分布特征及其關(guān)系,并通過(guò)線粒體標(biāo)記對(duì)其進(jìn)行遺傳變異檢測(cè),為揭示大渡河理縣高原鰍的遺傳背景及今后的種質(zhì)資源保護(hù)提供理論基礎(chǔ)。

1 材料與方法

1.1 材料

于2021年在大渡河采集理縣高原鰍樣本,金川縣沙耳鄉(xiāng)采集30尾(31°30′50.3″ E, 102°38′5.4″ N,編號(hào)g1：1～30),大渡河石棉縣安順彝族鄉(xiāng)采集34尾(29°15′8.3″ E, 102°19′6.7″ N,編號(hào)g 2：31～64),大渡河瀘定縣冷磧鎮(zhèn)采集30尾(102°13′20.9″ E,29°46′57.3″ N,編號(hào)g 3：65～94),丹巴縣格宗鎮(zhèn)采集23尾(30°48′26.6″ E, 101°56′25.2″ N,編號(hào)g 4：95～117),共117尾,剪取鰭條保存在100%酒精管,并存于4 ℃冰箱直至提取組織基因組。采樣點(diǎn)分布如圖1所示,水系圖來(lái)源于四川測(cè)繪地理信息局(http：//scsm.mnr.gov.cn/index.htm)。丹巴縣格宗鎮(zhèn)23尾理縣高原鰍年齡段為3～9齡,大渡河石棉縣34尾年齡段為6～11齡,瀘定縣冷磧鎮(zhèn)30尾年齡段為6～11齡,金川縣沙耳鄉(xiāng)30尾年齡段為6～11齡,4個(gè)種群理縣高原鰍年齡段均集中在6～7齡。大渡河石棉縣34尾,包括14尾雌魚和20尾雄魚,瀘定縣冷磧鎮(zhèn)30尾,包括13尾雌魚和17尾雄魚,金川縣沙耳鄉(xiāng)30尾,包括13尾雌魚和17尾雄魚,丹巴縣格宗鎮(zhèn)23尾魚,包含7尾雌魚和10尾雄魚(其中6尾幼齡魚無(wú)法鑒定性別)。雌雄魚性腺發(fā)育期集中在III～I(xiàn)V期。

1.2 理縣高原鰍體長(zhǎng)-體重的相關(guān)性及肥滿度分析

理縣高原鰍的體長(zhǎng)-體重關(guān)系用冪函數(shù)公式擬合[11]：

W=aLb

圖1 理縣高原鰍調(diào)查采樣站位Fig.1 Sampling stations of T. lixianensis

式中,W為體重(g);L為體長(zhǎng)(cm);a為條件因子,b為異速生長(zhǎng)因子[12]。測(cè)量理縣高原鰍全長(zhǎng)時(shí),取吻端到尾鰭末端的最大長(zhǎng)度;體長(zhǎng)即標(biāo)準(zhǔn)長(zhǎng),為吻端到尾鰭基的直線長(zhǎng)度,即全長(zhǎng)減去尾鰭的長(zhǎng)度。將體重和體長(zhǎng)數(shù)據(jù)進(jìn)行對(duì)數(shù)轉(zhuǎn)化,利用SPSS 19.0軟件,使用線性回歸法進(jìn)行擬合,估計(jì)參數(shù)a和b的95%置信區(qū)間(CI)及R2的顯著性水平。

肥滿度采用Fulton狀態(tài)指數(shù)K計(jì)算：

K=(W/L3)×100

式中,W為體重(g);L為體長(zhǎng)(cm)。數(shù)據(jù)以“平均數(shù)±標(biāo)準(zhǔn)差”表示。

1.3 理縣高原鰍基因組 DNA制備

理縣高原鰍基因組DNA提取采用CTAB高鹽法,取1 μL樣品通過(guò)電泳檢測(cè)DNA質(zhì)量,并用超微量分光光度計(jì)(Implen N50 touch)檢測(cè)DNA的純度和含量,保存于-20 ℃冰箱備用。

1.4 理縣高原鰍線粒體D-loop序列引物設(shè)計(jì)、合成

根據(jù)理縣高原鰍線粒體序列信息(GenBank accession number： KT966735.1),在其D-loop區(qū)上下設(shè)計(jì)特異性引物,基因引物序列F：5’-CTCACCCCTGGCTCCCAAAGC-3’,R：5’-AGGGCCCGTCTTAGCATCTTC-3’。

1.5 理縣高原鰍線粒體DNA PCR擴(kuò)增體系及程序

理縣高原鰍線粒體DNA擴(kuò)增體系為1 μL基因組DNA(50 ng),2×Premix LATaq聚合酶25 μL,正反向引物各1.5 μL(10 μmol/L),PCR級(jí)水21 μL,總擴(kuò)增體系50 μL。PCR反應(yīng)程序包括：94 ℃預(yù)變性5 min;94 ℃變性30 s,56 ℃退火30 s,72 ℃延伸1 min,32個(gè)循環(huán);72 ℃延伸8 min。PCR產(chǎn)物用1.2%的瓊脂糖凝膠電泳檢測(cè),將片段大小一致且條帶單一的擴(kuò)增產(chǎn)物進(jìn)行測(cè)序。

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

測(cè)序結(jié)果通過(guò)BLAST程序與GenBank中理縣高原鰍線粒體DNA序列進(jìn)行比對(duì)分析,確認(rèn)擴(kuò)增產(chǎn)物為理縣高原鰍線粒體DNA D-loop區(qū)序列。多樣性相關(guān)參數(shù),包括核苷酸多樣度(Pi)、單倍型多樣度(Hd)、多態(tài)位點(diǎn)數(shù)(S)、單倍型數(shù)(H)和平均核苷酸變異數(shù)(K)等,用DNAsp 5軟件進(jìn)行計(jì)算。群體遺傳分化指數(shù)Fst值利用Arlequin 2.0軟件統(tǒng)計(jì)。系統(tǒng)發(fā)育樹(shù)利用MEGA 6軟件中的鄰接法(Neighbor-Joining)進(jìn)行構(gòu)建。使用Arlequin 2.0軟件計(jì)算錯(cuò)配分布(Mismatch distribution)、Tajima’sD和Fu’sFs值來(lái)檢驗(yàn)群體歷史上是否經(jīng)歷擴(kuò)張, 通過(guò)偏差平方和(Sum of squared deviation, SSD)檢測(cè)核苷酸不配對(duì)分布與種群擴(kuò)張模型下期望分布之間的擬合優(yōu)度。

2 結(jié)果與分析

2.1 理縣高原鰍的體重-體長(zhǎng)關(guān)系

從圖2可知,大渡河理縣高原鰍全長(zhǎng)7.5～18.5 cm,體長(zhǎng)6.3～15.4 cm,平均體長(zhǎng)11.5 cm,體重2.8～44.9 g,平均體重19.6 g。a值為0.0123,95%置信區(qū)間為0.0081～0.0187,b值為2.9734,95%置信區(qū)間為2.8030～3.1440,R2>0.9121(表1)。本研究采用Fulton狀態(tài)指數(shù)K表示大渡河理縣高原鰍的肥滿度,理縣高原鰍的肥滿度為1.18±0.02(平均數(shù)±標(biāo)準(zhǔn)差)。

2.2 理縣高原鰍線粒體DNA控制區(qū)堿基特征

117尾理縣高原鰍線粒體DNA控制區(qū)序列PCR產(chǎn)物特異性較好。對(duì)117尾理縣高原鰍個(gè)體的測(cè)序結(jié)果比對(duì)分析發(fā)現(xiàn),去掉引物及部分端部序列,同源重排后的長(zhǎng)度為878 bp。序列中T、C、A、G堿基所占比例分別為32.5%、21.4%、32.9%、13.2%,AT和CG所占比例為65.4%和34.6%,可見(jiàn)AT比例顯著高于CG比例(表2),說(shuō)明理縣高原鰍線粒體DNA組成符合硬骨魚線粒體DNA堿基組成的特征。

2.3 理縣高原鰍不同地理種群的遺傳結(jié)構(gòu)

金川縣、石棉縣、瀘定縣和丹巴縣4個(gè)理縣高原鰍種群D-loop區(qū)共存在32個(gè)變異位點(diǎn)和28個(gè)單倍型(表3)。4個(gè)理縣高原鰍群體的單倍型多樣度(Hd)在0.637～0.921,丹巴縣格宗鎮(zhèn)樣本單倍型多樣度最高。理縣高原鰍群體的核苷酸多樣性在0.000 91～0.005 70,大渡河瀘定縣冷磧鎮(zhèn)樣本核苷酸多樣性最低(表3)。金川縣沙耳鄉(xiāng)樣本理縣高原鰍有11個(gè)單倍型,大渡河石棉縣安順彝族鄉(xiāng)樣本有7個(gè)單倍型,大渡河瀘定縣冷磧鎮(zhèn)樣本有7個(gè)單倍型,丹巴縣格宗鎮(zhèn)理縣高原鰍樣本有11個(gè)單倍型。單倍型Hap1出現(xiàn)的頻率最高,但4個(gè)理縣高原鰍地理種群未有共享的單倍型。金川縣和丹巴縣2個(gè)理縣高原鰍群體之間的個(gè)體存在共享的單倍型,石棉縣和瀘定縣2個(gè)種群之間也存在共享單倍型,但金川縣和石棉縣或?yàn)o定縣群體之間的個(gè)體未出現(xiàn)共享單倍型,丹巴縣和石棉縣或?yàn)o定縣理縣高原鰍群體之間的個(gè)體也未出現(xiàn)共享單倍型(表4)。

圖2 理縣高原鰍體長(zhǎng)和體重的回歸曲線Fig.2 Regression curves of relationship between body weight and body length of T. lixianensis

從表5可知,金川縣沙耳鄉(xiāng)和石棉及大渡河瀘定縣冷磧鎮(zhèn)的理縣高原鰍樣本遺傳距離分別為0.251 73和0.275 55,說(shuō)明發(fā)生了中度遺傳分化,與丹巴縣格宗鎮(zhèn)樣本之間的遺傳距離為-0.001 39,表明未出現(xiàn)遺傳分化。石棉和大渡河瀘定縣及丹巴縣格宗鎮(zhèn)理縣高原鰍樣本之間的遺傳距離為0.025 99和0.231 92,大渡河瀘定縣和丹巴縣格宗鎮(zhèn)樣本之間的遺傳距離為0.257 09,說(shuō)明發(fā)生了遺傳分化。

表1 理縣高原鰍體重和體長(zhǎng)關(guān)系函數(shù)的a和b參數(shù)

表2 不同地理種群理縣高原鰍mtDNA控制區(qū)序列堿基組成

2.4 理縣高原鰍不同地理種群的系統(tǒng)進(jìn)化樹(shù)

利用MEGA建立理縣高原鰍4個(gè)地理種群個(gè)體的NJ分子系統(tǒng)樹(shù)(圖3),從系統(tǒng)樹(shù)可以看出,117個(gè)理縣高原鰍個(gè)體在進(jìn)化樹(shù)上形成2個(gè)獨(dú)立的分支。其中,金川縣沙耳鄉(xiāng)和丹巴縣格宗鎮(zhèn)的個(gè)體形成一個(gè)獨(dú)立分支,而大渡河石棉縣安順彝族鄉(xiāng)和瀘定縣冷磧鎮(zhèn)的個(gè)體形成另一個(gè)分支,與各地理種群的遺傳距離結(jié)果一致。

表3 理縣高原鰍不同地理種群的遺傳多樣性參數(shù)

表4 理縣高原鰍不同種群的28個(gè)單倍型分布

表5 理縣高原鰍4個(gè)地理種群的遺傳分化Fst值

金川縣：1～30;石棉縣：31～64;瀘定縣：65～94;丹巴縣：95～117。Jinchuan county： 1-30; Shimian county： 31-64; Luding county： 65-94; Danba county： 95-117.圖3 基于線粒體D-loop區(qū)序列構(gòu)建的NJ 系統(tǒng)樹(shù)Fig.3 NJ tree based on D-loop region of mitochondrial DNA

2.5 理縣高原鰍種群歷史動(dòng)態(tài)分析

中性檢驗(yàn)結(jié)果顯示, 金川縣和丹巴縣理縣高原鰍地理群體的Tajima’sD值大于0,而石棉縣和瀘定縣地理群體的Tajima’sD值小于0,且均未達(dá)到顯著水平(P>0.05)。Fu’sFs檢驗(yàn)結(jié)果顯示, 4個(gè)理縣高原鰍地理群體的Fu’sFs均小于0,其中,石棉縣和瀘定縣群體均呈顯著負(fù)值(P<0.05), 金川縣和丹巴縣群體均呈不顯著負(fù)值(P>0.05,表6)。

表6 理縣高原鰍4個(gè)地理種群的中性檢驗(yàn)

柱狀圖代表出現(xiàn)頻率的觀測(cè)值,曲線代表基于種群擴(kuò)張事件模型的期望頻率。A：金川縣,B：石棉縣,C：瀘定縣,D：丹巴縣。The histograms indicate the observed frequencies, and the lines represent the expected frequencies based on the model of sudden population expansion. A： Jinchuan county, B： Shimian county, C： Luding county, D： Danba county.圖4 理縣高原鰍4個(gè)地理種群的單倍型錯(cuò)配分布 Fig.4 Mismatch distribution of haplotypes in four populations of T. lixianensis

對(duì)單倍型錯(cuò)配分布分析發(fā)現(xiàn),金川縣和丹巴縣理縣高原鰍地理群體的單倍型錯(cuò)配分布呈多峰分布(圖4),同時(shí),錯(cuò)配分析SSD不顯著;而石棉縣和瀘定縣理縣高原鰍地理群體的單倍型錯(cuò)配分布呈單峰分布且SSD不顯著,表明2個(gè)種群可能經(jīng)歷過(guò)群體擴(kuò)張(圖4-B、4-C)。

3 討 論

本研究中,理縣高原鰍的生長(zhǎng)參數(shù)b值為2.9734,當(dāng)b值小于3時(shí),種群為負(fù)異速生長(zhǎng), 即體長(zhǎng)的增速大于體重的增速,當(dāng)b值等于3 時(shí),種群為等速生長(zhǎng),當(dāng)b值大于3 時(shí),種群為正異速生長(zhǎng),體型表現(xiàn)為肥滿而不是纖細(xì)[13]。本研究表明,大渡河理縣高原鰍的b值保持在2.5～3.5的預(yù)期范圍內(nèi),接近3,說(shuō)明采樣時(shí)間段(10—12月)理縣高原鰍生長(zhǎng)狀況良好。由于缺乏理縣高原鰍體重體長(zhǎng)關(guān)系的文獻(xiàn)報(bào)道,因此本研究結(jié)果無(wú)法進(jìn)行比較。但研究顯示,阿克蘇河小鰾高原鰍b值為2.876[14],塔里木河巨頭高原鰍b值為2.583[14],伊犁河斯氏高原鰍b值為1.22[15],伊犁河黑斑高原鰍b值為1.113[15],烏倫古河小體高原鰍b值為 2.52[16],湟水河似鲇高原鰍b值為2.67[17]。大多數(shù)高原鰍屬魚類b值小于3,處于負(fù)異速生長(zhǎng),而大渡河理縣高原鰍b值高于其同屬的其他高原鰍。b值的差異可以歸因于多種因素,如檢測(cè)的標(biāo)本數(shù)量、捕獲個(gè)體的長(zhǎng)度范圍、棲息地類型、季節(jié)性和地理區(qū)域[16,18-20]。生長(zhǎng)環(huán)境是顯著影響b值的因素之一[21]。Jisr等[22]發(fā)現(xiàn),棕線石斑魚等的生長(zhǎng)參數(shù)b值在高溫環(huán)境下明顯高于低溫環(huán)境。Sileesh等[23]對(duì)黑頭魚、深海鯡,蛇鯖等研究也表明b值水平顯著受環(huán)境因素影響。

本研究采用線粒體基因組中進(jìn)化速率最快的控制區(qū)序列分析方法,研究理縣高原鰍在大渡河的4個(gè)不同地理種群共117個(gè)個(gè)體的堿基組成和變異情況,獲得理縣高原鰍野生群體的種質(zhì)資源遺傳變異基本數(shù)據(jù)。117個(gè)理縣高原鰍個(gè)體的mt DNA 控制區(qū)序列中,4種堿基組成呈現(xiàn)出一定的偏向性,各群體均存在A+T所占比例高于C+G的現(xiàn)象,這與絕大多數(shù)魚類的線粒體控制區(qū)組成特征相符[24]。

核苷酸多樣性指數(shù)(Pi)和單倍型多樣性指數(shù)(Hd)是線粒體DNA用來(lái)評(píng)估群體遺傳多樣性的2個(gè)重要參數(shù)。通過(guò)對(duì)理縣高原鰍D-loop區(qū)序列分析發(fā)現(xiàn),其總體的Hd在0.8以上,處于較高水平,而Pi偏低,僅為0.000 91～0.005 70。說(shuō)明,目前理縣高原鰍野生群體存在單倍型多樣性較高而核苷酸多樣性較低的情況,推測(cè)理縣高原鰍種群可能正從一個(gè)更小的有效種群經(jīng)歷種群擴(kuò)張,它們?cè)谳^短的時(shí)間內(nèi),通過(guò)堿基突變積累其單倍型多態(tài)性,但不足以使核苷酸多態(tài)性也處于較高狀態(tài)。由于核苷酸多樣性指數(shù)考慮各單倍型在群體中所占的比例,因而較單倍型多樣性指數(shù)更可靠[25]。類似結(jié)果在同屬的其余高原鰍種群中也有發(fā)現(xiàn)。晁燕等[26]利用線粒體DNA cytb序列對(duì)大通河區(qū)域的似鲇高原鰍進(jìn)行遺傳多樣性研究,結(jié)果顯示似鲇高原鰍大通河種群的Hd為0.86,Pi為0.0036。顏淵等[27]利用D-loop區(qū)對(duì)西部黃河流域不同水域的黃河高原鰍種群研究發(fā)現(xiàn),其Hd在0.9420～1.0000,Pi為0.0030～0.0043;楊成等[28]對(duì)擬硬刺高原鰍黃河種群和寶庫(kù)河種群進(jìn)行研究,發(fā)現(xiàn)其Hd為0.65～0.89,Pi為0.0010～0.0025;姚雁鴻等[29]利用線粒體D-loop序列對(duì)湖南烏龍山湘西盲高原鰍線粒進(jìn)行分析,發(fā)現(xiàn)其Hd在0.083～0.282,其Pi范圍為0.000 09～0.000 32?？梢?jiàn),以核苷酸多樣性指數(shù)判斷,目前理縣高原鰍屬的物種大部分都處于遺傳多樣性較低的狀態(tài),需要進(jìn)行有效的遺傳管理來(lái)增加其種群的遺傳多樣性。而從不同地理種群多樣性差異分析,金川縣和丹巴縣理縣高原鰍種群核苷酸多樣度遠(yuǎn)遠(yuǎn)高于石棉縣和瀘定縣種群。金川縣和丹巴縣屬于藏區(qū),人類活動(dòng)較少,無(wú)人為干擾的生境保留了其遺傳信息。前期實(shí)地調(diào)查也發(fā)現(xiàn),石棉縣和瀘定縣人為捕撈等活動(dòng)更頻繁,導(dǎo)致該地段魚類種群的多樣性偏低。

通過(guò)構(gòu)建進(jìn)化樹(shù)可以發(fā)現(xiàn),4個(gè)理縣高原鰍地理種群的個(gè)體在NJ系統(tǒng)樹(shù)上形成2大分支,推測(cè)這117個(gè)個(gè)體的單倍型可能來(lái)源于2個(gè)不同遺傳背景的祖先。其中,金川縣沙耳鄉(xiāng)和丹巴縣格宗鎮(zhèn)的理縣高原鰍個(gè)體形成一個(gè)獨(dú)立分支,而大渡河石棉縣安順彝族鄉(xiāng)和瀘定縣冷磧鎮(zhèn)的個(gè)體形成另一個(gè)分支,說(shuō)明金川縣、丹巴縣群體和石棉縣、瀘定縣群體形成明顯的分化。單倍型分布、遺傳分化和聚類分析結(jié)果均顯示大渡河金川縣和丹巴縣2個(gè)理縣高原鰍種群之間存在廣泛的基因交流,石棉縣和瀘定縣之間的種群也存在基因交流。金川縣和石棉縣或?yàn)o定縣之間,丹巴縣和石棉縣或?yàn)o定縣之間的理縣高原鰍種群均未發(fā)生基因交流,群體間存在顯著的遺傳分化,基于控制區(qū)序列的系統(tǒng)發(fā)育樹(shù)也顯示了兩兩種群之間形成獨(dú)立分支,也支持上述結(jié)論。4個(gè)理縣高原鰍群體的中性檢驗(yàn)結(jié)果表明,金川縣和丹巴縣的Tajima’sD值大于0,Fu’sFs小于0,但均未達(dá)到顯著水平,說(shuō)明未檢測(cè)到群體擴(kuò)張跡象,可能未發(fā)生群體擴(kuò)張事件。而石棉縣和瀘定縣理縣高原鰍地理群體的Tajima’sD和Fu’sFs均小于0,且后者呈現(xiàn)顯著負(fù)值,表明這2個(gè)群體可能經(jīng)歷過(guò)群體快速擴(kuò)張事件。石棉縣和瀘定縣理縣高原鰍地理群體的歧點(diǎn)分布呈單峰分布, 且SSD也不顯著,進(jìn)一步說(shuō)明該2個(gè)群體動(dòng)態(tài)未顯著偏離原假設(shè)(種群擴(kuò)張模型),近期都經(jīng)歷了快速的群體擴(kuò)張。

目前,大渡河丹巴縣和瀘定縣之間暫無(wú)已建電站,無(wú)人為因素阻礙該2個(gè)區(qū)域理縣高原鰍的基因交流;而丹巴縣和金川縣及石棉縣之間存在數(shù)個(gè)已建和在建電站,包括猴子巖水電站、長(zhǎng)河壩水電站、黃金坪水電站、瀘定水電站等,可能在一定程度上阻礙了壩上壩下魚類的基因交流。對(duì)于物種遺傳多樣性和物種完整性的保護(hù)在于客觀鑒別物種的關(guān)鍵遺傳單元,然后采用合適的方式管理這些單元[30]。

4 結(jié) 論

目前大渡河流域石棉至金川段理縣高原鰍生長(zhǎng)狀態(tài)良好。線粒體分子標(biāo)記分析發(fā)現(xiàn),大渡河理縣高原鰍種群遺傳多樣性處于較低水平,并建議將大渡河金川縣和丹巴縣種群作為一個(gè)整體進(jìn)行管理和保護(hù),石棉縣和瀘定縣也作為一個(gè)整體進(jìn)行管理和保護(hù),而金川縣和丹巴縣、石棉縣和瀘定縣的種群分開(kāi)管理和保護(hù)。