雅魯藏布江中游仔稚魚種類組成與分布特征

郭祉賓，楊德國，陳 亮，顏文斌，朱挺兵

(中國水產科學研究院長江水產研究所，農業農村部淡水生物多樣性保護重點實驗室，武漢 430223)

仔稚魚是魚類早期資源重要的生活史階段。開展仔稚魚調查有助于評估漁業種群資源量和了解魚類種群數量的動態變化，同時通過了解影響特定魚類早期資源分布的環境條件，也可進一步研究魚類群落整體動態變化[10]。開展仔稚魚資源調查的難點之一在于物種鑒定，傳統的形態學鑒定不僅需要詳實的水域魚類群落背景資料，還需要足夠的形態學參考文獻以及豐富的鑒定經驗[11]，面對同一科或同一屬不同物種的鑒定，也可能出現形態學特征不足以區分到種級別的情況[12]。而且，目前我國存有背景資料的仔稚魚形態學鑒定資料僅為我國已知魚類種數的1%左右，因此僅僅依靠形態學鑒定方法很難獲得全面準確的物種信息。相比之下，采用DNA分子標記進行分子生物學鑒定則具有操作較簡單且更易發現近緣種等優勢[13-14]。雖然目前分子生物學鑒定技術還存在鑒定材料無法重復利用，或由于固定劑的不同導致同一批樣品無法使用不同鑒定手段等缺陷[15]，但該技術仍為魚類早期發育階段重要的鑒定手段且發展迅速。

目前有關西藏魚類早期資源的研究，多集中在人工控制條件下的魚類早期發育研究[16-18]，而自然條件下仔稚魚的種類組成和分布特征還未有報道。因此，本研究采用形態學鑒定和分子生物學方法，對采自雅魯藏布江中游的仔稚魚樣本進行了種類鑒定，同時對主要種類的分布進行了研究，以為雅魯藏布江漁業早期資源保護提供科學依據。

1 研究方法

1.1 調查范圍與采樣點設置

本研究的調查范圍為雅魯藏布江中游。分別于2015年4月、2015年6月和2016年4月開展野外調查。其中，2015年設置了18個采樣點，2016年為補充調查，設置了6個采樣點(圖1、表1)。采樣點的設置遵循均勻分布、有明顯仔稚魚群體聚集、交通相對便利的原則，從下至上沿雅魯藏布江中游巡查，選擇有河漫灘及沙洲、仔稚魚較為密集的水域作為采樣點。采樣點主要環境因子見表2。

表1 年雅魯藏布江中游仔稚魚采樣點序號、海拔及經緯度Tab.1 Serial number，altitude，longitude and latitude of sampling stations in the middle reaches of the Yarlung Zangbo River

表2 雅魯藏布江中游環境因子Tab.2 Environmental factors in the middle reaches of the Yarlung Zangbo River

圖1 雅魯藏布江中游仔稚魚采樣點示意圖Fig.1 Fish larvae and juveniles sampling stations in the middle reaches of the Yarlung Zangbo River

1.2 樣本采集方法

使用圓錐形手抄網(網口直徑15 cm，網長40 cm)在采樣點水域淺水區、沙洲洄水灣附近隨機捕撈可視的仔稚魚，捕撈時間約10～15 min，并盡量確保所有種類的仔稚魚均被采集。仔稚魚樣本用純酒精固定保存于50 mL離心管中后帶回實驗室。同時記錄采樣點位置信息，測量海拔、水溫、流速、溶解氧、pH、濁度、底質類型等基礎環境數據。

1.3 仔稚魚種類鑒定

首先對仔稚魚樣本進行初步的形態學鑒定，然后基于形態學結果將各采樣點采集到的仔稚魚每個采樣點每種各選取3～5尾，整尾提取DNA，進行Cytb基因序列測序，用于PCR擴增的Cytb基因引物為通用引物L14724(5′-GAC TTG AAA AAC CAC CGT TG-3′)和H15915(5′-CTC CGA TCT CCG GAT TAC AAG AC-3′)。去除引物兩端可能產生干擾的錯讀堿基序列后與NCBI在線網站序列進行比對，以進一步鑒定種類。將送檢的仔稚魚編號，方法為采樣點編號-尾數編號或采樣點編號-各站形態學鑒定種類序號-尾數編號。選擇BLAST結果相似度較高的序列，使用MEGA5.0軟件，ML法構建系統發育樹分析，當Cytb基因序列相似度≥99%時認為是同一種，序列相似度在92%～99%時認為是同一屬，序列相似度在85%～92%時認為是同一科[19]。

2 結果

2.1 雅魯藏布江中游仔稚魚種類組成

共采集到仔稚魚樣本3 502尾，其中2 015年采集到仔稚魚1 867尾，2016年共采集到仔稚魚1 635尾。共對157尾仔稚魚樣本進行了Cytb基因測序。經過比對，共鑒定出仔稚魚12種，隸屬于2科5屬。其中138尾樣本可鑒定到種，2尾樣本可鑒定到屬，17尾樣本可鑒定到科。2015年樣本共鑒定出仔稚魚10種，2016年樣本共鑒定出仔稚魚6種(表3)。選取NCBI數據庫上BLAST結果相近的幾條DNA與測序結果構建ML系統發育樹如圖2、圖3，結果顯示分子鑒定結果基本可信。

表3 雅魯藏布江中游仔稚魚物種組成Tab.3 Species composition of fish larvae and juveniles in the middle reaches of the Yarlung Zangbo River

圖2 2015年雅魯藏布江中游仔稚魚Cyt b基因系統發育樹Fig.2 Phylogenetic tree using Cyt b gene of fish larvae and juveniles in the middle reaches of the Yarlung Zangbo River in 2015

圖3 2016年雅魯藏布江中游仔稚魚Cyt b基因系統發育樹Fig.3 Phylogenetic tree using Cyt b gene of fish larvae and juveniles in the middle reaches of the Yarlung Zangbo River in 2016PZ：采樣點B1，ZB：采樣點B2，LL：采樣點B3，WLJZ：采樣點B4，BZ：采樣點B5，ZD：采樣點B6

2.2 仔稚魚種類數及尾數的空間變化

2015年仔稚魚捕獲尾數最多的采樣點為A8，共捕獲仔稚魚352尾，最少的采樣點為A12、A13，各捕獲仔稚魚4尾；2016年仔稚魚捕獲尾數最多的采樣點為B6，共捕獲649尾，最少的采樣點為B2，共捕獲15尾。2015年仔稚魚種類數分布最多的采樣點為A8，共出現6種仔稚魚；最低的采樣點為A14、A15，分別僅出現1種仔稚魚；2016年仔稚魚種類數分布最多的采樣點為B6，共出現5種仔稚魚，最少的采樣點為B2、B3，各出現2種仔稚魚(圖4、圖5)。

圖4 2015年雅魯藏布江中游仔稚魚的尾數及種類數量分布情況Fig.4 Distribution of the number of individuals and species of fish larvae and juveniles of the middle reaches of the Yarlung Zangbo River in 2015

圖5 2016年雅魯藏布江中游仔稚魚的尾數及種類數量分布情況Fig.5 Distribution of the number of individuals and species of fish larvae and juveniles of the middle reaches of the Yarlung Zangbo River in 2016

2.3 主要種類的空間分布

由表4、表5可見，2015年分布最廣的仔稚魚物種為細尾高原鰍，18個采樣點中有11個采樣點出現，其次為拉薩裸裂尻魚，共有9個采樣點中出現。分布最狹窄的為拉薩裂腹魚，僅在A8一個采樣點出現，其次為異尾高原鰍、雙須葉須魚及尖裸鯉，各在2個采樣點中出現。2016年分布最廣的仔稚魚也為細尾高原鰍和拉薩裸裂尻魚，6個采樣點中5個采樣點出現。分布最狹窄的則為異齒裂腹魚和高原鰍屬一種，各在B5、B6采樣點中出現。

表4 雅魯藏布江中游2015年各采樣點仔稚魚種類的空間分布情況Tab.4 Spatial distribution of species in each sampling station in the middle reaches of the Yarlung Zangbo River in 2015

表5 雅魯藏布江中游2016年各采樣點仔稚魚種類的空間分布情況Tab.5 Spatial distribution of species in each sampling station in the midle reaches of the Yarlung Zangbo River in 2016

2.4 環境因子對雅魯藏布江中游仔稚魚種類分布的影響

本研究中，裂腹魚類仔稚魚種類數分布隨海拔梯度存在先增高再降低的傾向，其余環境因子則對本研究中仔稚魚分布造成的影響并不顯著(表6)。本研究2015年海拔2 908～2 921 m(A1、A2)未采集到裂腹魚類仔稚魚，而海拔3 008、3 555 m處(A6、A8)裂腹魚類仔稚魚物種豐富度達到最大(圖6a)。2016年設置采樣點較少，處于海拔3 100 m以下，仔稚魚物種數量呈現隨海拔梯度增高而增加的現象。，且2015年未采集到裂腹魚類的樣點附近(B1、B2)也采集到了裂腹魚類仔稚魚(圖6b)。而高原鰍分布則較為廣泛，基本各海拔都采集到了高原鰍仔稚魚樣本(圖7)。

表6 雅魯藏布江中游仔稚魚物種數、尾數與環境因子間的相關性分析Tab.6 Correlation analysis between species number，individual number and environmental factors of the middle reaches of the Yarlung Zangbo River

圖6 雅魯藏布江中游裂腹魚類仔稚魚種類數分布Fig.6 Distribution of the number of species of fish larvae and juveniles in the family Cyprinidae of the middle reaches of the Yarlung Zangbo River

圖7 雅魯藏布江中游高原鰍類仔稚魚種類數分布Fig.7 Distribution of the number of species of fish larvae and juveniles in the family Balitoridae of the middle reaches of the Yarlung Zangbo River

3 討論

3.1 基于Cyt b基因的雅魯藏布江中游仔稚魚種類鑒定

基于Cytb基因等DNA保守片段進行物種鑒定的技術是發展較為迅速的新興生物鑒定技術，且其在水生生物尤其是仔稚魚鑒定中得到了較為廣泛的應用，如DANIEL等[20]使用形態學-分子生物學復合鑒定技術對馬尾藻海的仔稚魚物種組成及多樣性進行了研究，PAUL等[12]基于分子生物學對奧地利維也納水庫仔稚魚物種特異性漂移現象進行了研究。較之傳統形態學鑒定方法，分子生物學鑒定技術擁有操作簡單，對鑒定人員要求相對較低等優點[13-14]。KO等[21]進行的形態學-分子生物學仔稚魚鑒定方法對比實驗中，發現形態學鑒定結果準確率遠低于分子生物學方法，且形態學鑒定結果存在較高的主觀性影響。

利用Cytb基因測序比對可以較準確地鑒定出雅魯藏布江中游仔稚魚的種類。本研究中，基于Cytb基因測序比對實現了種級水平上87.9%(157尾識別出138尾)的鑒定成功率，科級水平100%的鑒定成功率。這表明在缺乏必要的仔稚魚形態學鑒定特征背景參考的前提下，依靠分子生物學鑒定方法進行種類鑒定是可行的。有類似研究也表明，裂腹魚類物種鑒定使用Cytb基因在分子鑒定研究中具有很大的潛力[22]。但需要注意的是，由于遺傳距離等原因，存在分子生物學鑒定方法成功率對高原魚類鑒定效果較差的情況[23，24]。針對這種情況，應結合形態學鑒定方法，并可嘗試選取其他DNA作進一步研究。另外，分子生物學鑒定方法雖然可以快速進行物種鑒定，但仍需可靠的DNA數據庫作為支撐。作為開放數據庫，GenBank擁有著足夠龐大的基因數據，然而其質量存在良莠不齊的問題[25]，因此傳統的形態學鑒定方法仍然是鑒定仔稚魚必要的程序。

3.2 雅魯藏布江中游仔稚魚的種類組成及分布范圍

本研究顯示，雅魯藏布江中游仔稚魚種類組成與近年來的雅魯藏布江漁業資源調查結果基本一致，主要為裂腹魚類和高原鰍類。本研究未發現鮡類仔稚魚，可能與本研究選擇的采樣點和時間有關。本研究一般選擇白天在沿江沙洲、洄水灣處采樣，而鮡類產卵場一般為卵石堆、石底處，且其一般存在較強的夜間活動習性[26]。本研究鑒定出的仔稚魚除裸腹葉須魚和怒江高原鰍外，均為雅魯藏布江土著魚類[27]。裸腹葉須魚自然分布于金沙江、瀾滄江、怒江等水系，怒江高原鰍僅分布于怒江水系，但沒有分布于雅魯藏布江的報道。考慮到上述兩種魚類均只有單個樣本，因此它們可能是從其他水域偶然被帶入雅魯藏布江，在雅魯藏布江中游建群的可能性較小。

本研究中，兩年間分布最廣的仔稚魚均為細尾高原鰍和拉薩裸裂尻魚，與李雷等[28]的研究結果相似。細尾高原鰍在青藏高原地區廣泛分布，而拉薩裸裂尻魚在雅魯藏布江流域歷來也是分布最為廣泛的魚類，這可能與他們的種群維持需求不甚嚴格有關。

由于缺少雅魯藏布江仔稚魚的歷史數據，因此選擇將雅魯藏布江成魚的歷史分布記錄與本研究進行了比較。結果顯示，雅魯藏布江中游仔稚魚的分布與成魚的分布存在一定的差異。本研究發現在海拔3 000 m左右處有尖裸鯉仔稚魚分布，同時進行的成魚分布調查結果則顯示[29]，2015年雅魯藏布江出現的裂腹魚類中，尖裸鯉成魚僅分布于雅魯藏布江仁布以上江段(海拔3 770 m以上)。可以認為雅魯藏布江中下游存在尖裸鯉小規模群落，或至少存在尖裸鯉親魚在此繁殖。仔稚魚群落規模與成魚可能存在不成比例的情況在其他水域也是較為常見的。例如，葛成岡等[30]對長江口鰻苗定置網兼捕的研究顯示，長江資源量嚴重衰退的刀鱭(Coilianasus)，其仔稚魚在定置網中的漁獲量顯示為最優勢種，高達31.99%。另外，游泳能力較差的仔稚魚分布主要取決于水流、親魚產卵場等[31，32]，與成魚分布位置可能存在差別。

3.3 環境因子對雅魯藏布江中游仔稚魚分布的影響

一般而言，環境因子會對物種的分布產生重大影響。仔稚魚由于其運動能力較弱，其分布易受環境因子影響，如溫度可以影響仔稚魚代謝、運動、基因表達及親魚繁殖生態，適當的濁度可為仔稚魚提供隱蔽等[32，33]。但本研究中仔稚魚的分布與其余環境因子的關系并不明顯，而與劉明典等[34]的研究結果相似，仔稚魚種類分布呈現沿海拔上升先增加而后減少的狀態。本研究仔稚魚分布結果受環境影響較少的原因可能是由于受交通條件限制，僅選擇了環境相對同質化較高的沙洲、洄水灣作為采樣點，導致各采樣點間除海拔以外的環境因子差異較小。

3.4 雅魯藏布江中游魚類早期資源保護建議

本研究中，雅魯藏布江中游各采樣點均有仔稚魚分布，表明整個區域內廣泛分布有符合魚類棲息條件的場所。其中，裂腹魚類仔稚魚種類數分布主要集中在海拔3 000 m～3 500 m的仲達-甲格江段、貢嘎-伊岡江段。高原鰍類仔稚魚則主要集中分布于海拔3 000 m以下的派鎮-尼洋河口江段以及4 400 m左右的拉孜-薩嘎江段。仔稚魚運動、感知、攝食、反捕食能力均相對較弱[30]，是整個生活史中死亡率最高的階段，因此必須給予必要的保護。建議在魚類繁殖育幼季節，盡可能減少雅魯藏布江中游兩岸淺水沙洲、洄水灣等水域的人為干擾，特別是采砂等直接破壞河床的活動。