鮑亞科似鮑屬魚類的物種界定和系統(tǒng)發(fā)育關(guān)系

曹 愷 李玉火 何 歡 李明月 傅萃長

(復(fù)旦大學(xué)生命科學(xué)院生物多樣性科學(xué)研究所, 生物多樣性與生態(tài)工程教育部重點實驗室, 上海 200438)

鮑亞科似鮑屬魚類的物種界定和系統(tǒng)發(fā)育關(guān)系

曹 愷 李玉火 何 歡 李明月 傅萃長

(復(fù)旦大學(xué)生命科學(xué)院生物多樣性科學(xué)研究所, 生物多樣性與生態(tài)工程教育部重點實驗室, 上海 200438)

亞科似鮑屬魚類的物種界定保持爭議, 系統(tǒng)發(fā)育關(guān)系尚待解決。研究取樣似鮑屬魚類所有種, 使用核基因多位點序列重建似屬魚類的系統(tǒng)發(fā)育關(guān)系, 運用分子的物種界定方法并結(jié)合形態(tài)特征分析厘定我國似屬魚類的分類。貝葉斯系統(tǒng)發(fā)育樹結(jié)果表明: 桂林似鮑與平江似是單系種; 似與扁嘴似是多系種, 前者包括五個譜系A(chǔ)至E, 后者包括兩個譜系A(chǔ)與B。POFAD距離分析和Structurama分析的結(jié)果表明似和扁嘴似鮑的每個譜系是獨立遺傳種群, BP & P分析結(jié)果強烈支持它們是不同種。*Beast物種樹結(jié)果揭示: 扁嘴似鮑譜系B位于似鮑屬魚類的基部位置; 似譜系A(chǔ)與B是姊妹群關(guān)系, 似譜系C是扁嘴似鮑譜系A(chǔ)的姊妹群, 它們一起與桂林似鮑形成姊妹群關(guān)系; 似譜系D與E是姊妹群關(guān)系, 它們一起是平江似鮑的姊妹群。結(jié)合形態(tài)證據(jù), 對我國似屬魚類分類厘定如下: 限定嚴格意義似鮑包括似譜系A(chǔ)+B; 恢復(fù)長吻似鮑(Pseudogobio longirostris Mori, 1934)給予似鮑譜系C分類名; 似譜系D與E是隱存種, 桂林似鮑與平江似鮑是有效種。

鯉科;亞科; 似鮑屬; 物種界定; 系統(tǒng)發(fā)育

似鮑屬魚類的物種分類一直存在爭議。B?n?rescu & Nalbant[1]認為似鮑屬魚類包括扁嘴似鮑1種含2亞種, 即指名亞種(Pseudogobio esocinus esocinus)與似鮑亞種(P. esocinus vaillanti)。羅云林等[5]對我國的似鮑屬魚類進行了分類厘訂, 認為我國似鮑屬魚類包括1種含3亞種: 似鮑指名亞種(Pseudogobio vaillanti vaillanti)、桂林似鮑(P. vaillanti guilinensis)與長吻似鮑(P. vaillanti longirostris Mori, 1934)。樂佩琦[2]認為長吻似鮑亞種是似鮑指名亞種的同物異名, 并把似鮑指名亞種與桂林似鮑都提升為種。楊秀平等[6]利用形態(tài)可量性狀使用主成分分析對我國似鮑屬魚類進行了形態(tài)比較,認為桂林似鮑為似鮑的同物異名, 我國似鮑屬魚類僅似鮑單種。

基于形態(tài)特征與分子證據(jù)重建的鮑亞科魚類系統(tǒng)發(fā)育關(guān)系研究[7—12]均支持似鮑屬魚類是一單系群, 但是似鮑屬魚類種間的系統(tǒng)發(fā)育關(guān)系并未解析清楚。夏曦中等[13]利用線粒體細胞色素b基因,對我國遼河、長江、錢塘江、珠江與海南島的廣義似鮑進行了譜系地理研究, 結(jié)果表明廣義似鮑具有明顯的空間遺傳結(jié)構(gòu), 能分為5大譜系, 其中兩個譜系分布于長江與錢塘江, 另外三個譜系各自分布于遼河、珠江與海南島, 譜系間遺傳距離達到6%—10%?；隰~類線粒體基因的種間平均遺傳距離一般為8%左右[14—16], 夏曦中等[13]的研究結(jié)果暗示我國似鮑屬魚類可能不止似鮑單種。Tominaga等[17]使用線粒體細胞色素b基因?qū)V義扁嘴似鮑進行了初步的譜系地理研究, 結(jié)果顯示廣義扁嘴似鮑是多系群, 包括兩大譜系: 一個譜系是日本本島特有的, 另一個譜系嵌套于似鮑內(nèi)部。Tominaga等[18]進一步加強取樣使用線粒體細胞色素b基因?qū)V義扁嘴似鮑進行了詳細的譜系地理研究, 結(jié)果表明廣義扁嘴似鮑能分為三大譜系A(chǔ)、B與C。其中, 譜系A(chǔ)與分布于我國北方的似鮑形成姊妹群關(guān)系; 譜系B與譜系C是姊妹群, 它們一起與分布于我國南方的似鮑形成姊妹群關(guān)系。

綜上所述, 似鮑屬魚類的物種分類保持爭議,系統(tǒng)發(fā)育關(guān)系尚待進一步解決。清晰的物種界定是準確理解生物多樣性空間格局形成過程的基礎(chǔ)。最近的研究[19—21]提倡分子證據(jù)基礎(chǔ)的物種界定方法并結(jié)合形態(tài)證據(jù)進行物種界定。因此, 本研究取樣似鮑屬魚類所有物種, 選擇5個核基因位點作為分子標記, 使用系統(tǒng)發(fā)育技術(shù)重建似鮑屬魚類的系統(tǒng)發(fā)育關(guān)系, 運用分子證據(jù)基礎(chǔ)的物種界定方法并結(jié)合形態(tài)證據(jù)進行分類厘定, 為將來以該屬魚類作為研究對象, 理解東亞地區(qū)淡水魚類多樣性空間格局形成過程奠定基礎(chǔ)。

1 材料與方法

1.1 物種取樣和分子標記

本研究共取樣似鮑屬魚類4種69尾(表 1)。其中廣義似鮑58尾, 來自長江、錢塘江、甌江、靈江、閩江、韓江、鴨綠江、大洋河與遼河水系; 平江似鮑1尾, 來自珠江水系; 桂林似鮑4尾, 來自珠江水系; 廣義扁嘴似鮑6尾, 來自日本琵琶湖與云出川水系。依據(jù)Tang等[11]重建的鮑亞科魚類分子系統(tǒng)發(fā)育關(guān)系, 選取長蛇鮑與棒花魚作為外類群。分子標記選擇5個單拷貝核基因: 外胚層神經(jīng)皮層蛋白1基因(ENC1)、錯配修復(fù)蛋白MutL基因3(MLH3)、錯配修復(fù)蛋白MutS基因6(MSH6)、核位點重組激活基因1(RAG1)和核糖體蛋白S7基因第一內(nèi)含子(rpS7 Intron 1)。

1.2 DNA提取、PCR擴增與產(chǎn)物測序

從95%酒精保存的肌肉組織中采用高鹽法提取基因組DNA[22]。PCR過程使用引物[23,24]見表 2,各基因PCR擴增反應(yīng)的條件如下: 94℃預(yù)變性5min; 94℃變性50s, 52℃退火50s, 72℃延伸60s; 循環(huán)40次; 72℃最后延伸8min。PCR產(chǎn)物純化后在ABI 3730 DNA測序儀中使用擴增引物進行測序。

1.3 序列處理與貝葉斯樹的構(gòu)建

使用DAMBE ver. 5.6.7軟件[25]進行蛋白質(zhì)編碼基因的序列對位, 使用MAFFT ver. 7.22軟件[26]進行非蛋白質(zhì)編碼基因?qū)ξ?。序列對位完成后手工校正并排除具有插入與缺失的區(qū)域。使用MEGA ver. 6.06軟件[27]計算基因的分子特征。蛋白質(zhì)編碼基因分區(qū)策略按密碼子進行, 非蛋白質(zhì)編碼基因不分區(qū)。使用jModeltest ver. 2.1.7軟件[28], 依據(jù)赤池信息準則選擇每個分區(qū)的最適堿基替換模型。在Mrbayes ver. 3.2.2軟件[29]下使用核基因多位點聯(lián)合序列構(gòu)建貝葉斯系統(tǒng)發(fā)育樹, 主要步驟如下: 基于馬爾可夫鏈蒙特卡羅進行5000萬代運算, 取樣頻率為1000, 保證分離頻率平均標準偏差小于0.01。運行完成后舍棄前1/4的取樣樹, 根據(jù)剩余的3/4取樣樹構(gòu)建50%的多數(shù)原則一致樹, 并計算每個節(jié)點的貝葉斯后驗概率值。

1.4 物種界定

在POFAD ver. 1.0.7軟件[30]中采用標準化矩陣計算基于核基因多位點的個體之間Genpofad距離,并在SplitsTree4 ver. 4.14.2軟件[31]中構(gòu)建鄰接網(wǎng)絡(luò)圖(簡稱POFAD距離分析)。在Structurama 2.0軟件[32]中采用Dirichlet 過程進行種群數(shù)最適值的測試, 運算過程中將其設(shè)為隨機變量并服從形狀參數(shù)為0.1,尺度參數(shù)為10的伽馬分布; 共運行5000萬代馬爾可夫鏈蒙特卡羅, 每1000代取樣一次, 舍棄前1/4的樣本; 最終獲得具有最高后驗概率的種群數(shù)(簡稱Structurama分析)。依據(jù)貝葉斯系統(tǒng)發(fā)育樹、POFAD距離分析以及Structurama分析的結(jié)果, 設(shè)定先驗物種。在BP & P 3.0軟件[33,34]中驗證種的成立與否。采用無指導(dǎo)樹物種界定方法計算指定物種劃分方法的后驗概率分布, 按theta和tau的先驗分布運行3組獨立分析: (1) thetaprior=1, 10, tauprior=1, 10; (2) thetaprior=2, 2000, tauprior=2, 2000; (3) thetaprior=1, 10, tauprior=2, 2000。共運行250萬代, 每5代取樣一次, 舍棄前1/5的樣本(簡稱BP & P分析)。

1.5 ＊BEAST物種樹的構(gòu)建

在*BEAST ver. 2.3.0軟件[35]中構(gòu)建物種樹, 主要步驟包括: 選擇對數(shù)正態(tài)分布的松散鐘模型作為分子鐘模型, 樹的先驗?zāi)Ｐ驮O(shè)定為Yule過程; 運行4個獨立重復(fù), MCMC代數(shù)設(shè)置為20000萬代, 取樣頻率為1000; 舍棄前1/4的樣本后, 在Log Combiner ver. 2.3.0軟件[35]中合并每個獨立重復(fù)的log文檔和trees文檔; 在Tracer ver. 1.5軟件[36]中估算合并后的log文檔參數(shù), 保證每個參數(shù)的有效取樣大小大于200; 在Tree Annotator ver. 2.3.0軟件[35]中獲得最大譜系置信樹。

表 1 種名、采樣地點、個體編號及GenBank登錄號Tab. 1 Species, sampling sites, code and GenBank accession number of all individuals

續(xù)表1

續(xù)表1

表 2 分子標記擴增使用的引物Tab. 2 Primers for the amplification of molecular markers

1.6 外部形態(tài)特征的計數(shù)與測量

統(tǒng)一采用魚體背部俯視方向左側(cè)進行計數(shù)與測量, 可數(shù)性狀依據(jù)陳宜瑜[37], 可量性狀依據(jù)Kottelat和Freyhof[38](圖 1)。在Past ver. 1.90[39]軟件下完成可量性狀數(shù)據(jù)的主成分分析。為消除異速生長的影響, 可量性狀數(shù)據(jù)進行對數(shù)轉(zhuǎn)換并使用Burnaby方法進行校正[40,41]。

2 結(jié)果

2.1 序列特征

本研究使用的5個單拷貝核基因一起序列總長7582 bp, 變異位點數(shù)749, 簡約信息位點數(shù)400。其中, ENC1序列長度為1290 bp, 變異位點數(shù)101, 簡約信息位點數(shù)58; MLH3序列長度為2223 bp, 變異位點數(shù)230, 簡約信息位點數(shù)115; MSH6序列長度為1992 bp, 變異位點數(shù)171, 簡約信息位點數(shù)105; RAG1序列長度為1449 bp, 變異位點數(shù)88, 簡約信息位點數(shù)47; rpS7 intron 1序列長度為628 bp, 變異位點數(shù)159, 簡約信息位點數(shù)75。所有個體的GenBank登錄號見表 1。

2.2 物種界定

圖 1 可量性狀示意圖(參考Kottelat & Jorg[38])Fig. 1 The diagram of morphometric characters (referring to Kottelat & Jorg[38])a. 體長; b. 頭長; c. 背側(cè)頭長; d. 體高; e. 背鰭前長; f. 背鰭后長; g. 腹鰭前長; h. 臀鰭前長; i. 胸腹鰭起點間距離; j. 腹鰭肛門起點間距離; k. 肛門臀鰭起點間距離; l. 背鰭長; m. 胸鰭長; n. 腹鰭長; o. 背鰭基長; p. 臀鰭基長; q. 尾柄高; r. 尾柄長; s. 吻長; t.眼后頭長; u. 眼徑; v. 眼間距(圖中未顯示)a. Standard length (SL); b. Head length (HL); c. Dorsal head length; d. Body depth; e. Pre-dorsal length; f. Post-dorsal length; g. Pre-pelvic length; h. Pre-anal length; i. Distance between pectoraland pelvic-fin origins; j. Distance between pelvic-fin and anus origins; k. Distance between anus and anal-fin origins; l. Length of dorsal fin; m. Length of pectoral fin; n. Length of pelvic fin; o. Length of dorsal-fin base; p. Length of anal-fin base; q. Depth of caudal peduncle; r. Length of caudal peduncle; s. Snout length; t. Postorbital length; u. Eye diameter; v. Eye diameter (not shown in this figure)

基于核基因多位點聯(lián)合序列重建的貝葉斯系統(tǒng)發(fā)育樹(圖 2)揭示: 廣義似鮑是多系, 包括五個單系的譜系A(chǔ)至E; 桂林似鮑與平江似鮑都是單系; 廣義扁嘴似鮑是多系, 包括兩個單系的譜系A(chǔ)和B。POFAD距離分析(圖 3)支持桂林似鮑、平江似鮑、廣義似鮑譜系A(chǔ)至E與廣義扁嘴似鮑譜系A(chǔ)與B的每個譜系是獨立遺傳種群。Structurama分析結(jié)果與POFAD距離分析結(jié)果一致。BP & P分析結(jié)果強烈支持桂林似鮑、平江似鮑、廣義似鮑譜系A(chǔ)至E以及廣義扁嘴似鮑譜系A(chǔ)與B的每個譜系是不同種(貝葉斯后驗概率100%)。

2.3 系統(tǒng)發(fā)育關(guān)系

基于核基因多位點聯(lián)合序列重建的貝葉斯系統(tǒng)發(fā)育樹(圖 2)與*BEAST物種樹(圖 4)樹形基本一致。最近的魚類分子系統(tǒng)發(fā)育關(guān)系研究[42—44]暗示*BEAST物種樹比基于核基因多位點聯(lián)合序列重建的貝葉斯系統(tǒng)發(fā)育樹更能準確地反映系統(tǒng)發(fā)育關(guān)系。因此, 下面僅描述*BEAST物種樹揭示的結(jié)果。

*BEAST物種樹(圖 4)揭示的似鮑屬魚類系統(tǒng)發(fā)育關(guān)系為: (1)扁嘴似鮑譜系B位于似鮑屬魚類的基部位置, 是所有其他似鮑屬魚類的姊妹群(貝葉斯后驗概率100%); (2)似鮑譜系A(chǔ)與B是姊妹群關(guān)系, 似鮑譜系C是扁嘴似鮑譜系A(chǔ)的姊妹群; (3)似鮑譜系A(chǔ)+B與似鮑譜系C+扁嘴似鮑譜系A(chǔ)形成姊妹群關(guān)系(貝葉斯后驗概率96%—100%), 它們一起是桂林似鮑的姊妹群(貝葉斯后驗概率62%); (4)似鮑譜系D與E是姊妹群關(guān)系, 它們一起是平江似鮑的姊妹群(貝葉斯后驗概率99%—100%); (5)平江似鮑+似鮑譜系D與E是桂林似鮑+扁嘴似鮑譜系A(chǔ)+似鮑譜系A(chǔ)至C的姊妹群(貝葉斯后驗概率99%)。

2.4 形態(tài)特征

似鮑、桂林似鮑與平江似鮑的形態(tài)特征見表 3?？闪啃誀畹闹鞒煞址治鼋Y(jié)果(圖 5)表明不能分開桂林似鮑、似鮑譜系A(chǔ)與B, 但能很好地分開似鮑譜系C、D與E; 桂林似鮑+似鮑譜系A(chǔ)與B能與似鮑譜系C至E分開。在可數(shù)性狀中(表 3), 側(cè)線鱗的數(shù)目能夠很好地分開桂林似鮑(43)與似鮑譜系A(chǔ)至E+平江似鮑(40—42)。似鮑譜系C側(cè)線鱗的數(shù)目(41—42)及背鰭前鱗的數(shù)目(12—13)和似鮑譜系A(chǔ)+B側(cè)線鱗的數(shù)目(40—41)及背鰭前鱗的數(shù)目(10—12)也有一定差異。腹鰭分枝鰭條數(shù)與背鰭前鱗的數(shù)目能夠區(qū)分平江似鮑(腹鰭分枝鰭條數(shù)8; 背鰭前鱗的數(shù)目15)和桂林似鮑+似鮑譜系A(chǔ)至E(腹鰭分枝鰭條數(shù)7; 背鰭前鱗的數(shù)目10—13)。似鮑譜系A(chǔ)與B比較、似鮑譜系D與E比較, 在可數(shù)性狀上幾乎沒有差異。

3 討論

3.1 我國似鮑屬魚類的分類厘定

楊秀平等[6]基于形態(tài)可量性狀的主成分分析結(jié)果認為我國似鮑屬魚類僅似鮑一種。然而, 本研究基于分子的物種界定分析結(jié)果支持桂林似鮑與平江似鮑是有效種, 并認為廣義似鮑譜系A(chǔ)至E以及廣義扁嘴似鮑譜系A(chǔ)與B的每個譜系是不同種。由于扁嘴似鮑譜系A(chǔ)與B缺乏形態(tài)數(shù)據(jù), 它們的分類有待進一步的研究再確定。下面結(jié)合形態(tài)證據(jù)就我國似鮑屬魚類的分類進行討論。

盡管分子證據(jù)支持分布于韓江至長江的似鮑譜系A(chǔ)與分布于長江水系的似鮑譜系B為不同種,然而可數(shù)性狀與可量性狀分析結(jié)果不支持把似鮑譜系A(chǔ)與B劃分為不同種。另外, 似鮑譜系A(chǔ)與B是姊妹群關(guān)系, 兩者在長江水系均有分布。因此, 從形態(tài)證據(jù)、系統(tǒng)發(fā)育關(guān)系與地理分布三方面結(jié)合來看, 似鮑譜系A(chǔ)與B應(yīng)屬于同一物種的不同遺傳譜系?；谒契U譜系A(chǔ)+B分布范圍包括似鮑模式地長江水系以及其形態(tài)性狀符合Sauvage[45]對似鮑模式種的形態(tài)描述, 本研究限定似鮑譜系A(chǔ)+B為嚴格意義似鮑。

圖 2 基于核基因多位點聯(lián)合序列重建的似鮑屬魚類貝葉斯系統(tǒng)發(fā)育樹Fig. 2 Bayesian tree of Pseudogobio fishes based on concatenated sequences of nuclear loci節(jié)點上方為貝葉斯后驗概率, 外類群沒有顯示Above nodes are Bayesian posterior probabilities, and outgroup taxa are not shown

本研究結(jié)果表明分布于我國北方遼河、大洋河與鴨綠江水系的似鮑譜系C與分布于日本的扁嘴似鮑譜系A(chǔ)形成姊妹群關(guān)系, 它們一起是嚴格意義似鮑(= 似鮑譜系A(chǔ)+B)的姊妹群(圖 4), 暗示似鮑譜系C是不同種。分子的物種界定分析結(jié)果支持似鮑譜系C作為種成立?？闪啃誀罘治鼋Y(jié)果能很好地區(qū)分似鮑譜系C與嚴格意義似鮑。在可數(shù)性狀中, 似鮑譜系C側(cè)線鱗的數(shù)目(41—42)及背鰭前鱗的數(shù)目(12—13)與嚴格意義似鮑側(cè)線鱗的數(shù)目(40—41)及背鰭前鱗的數(shù)目(10—12)也有一定差異。依據(jù)羅云林等[5]對我國似鮑屬魚類的分類厘訂, 似鮑譜系C為他們描述的長吻似鮑亞種。因此,結(jié)合分子證據(jù)、形態(tài)證據(jù)、系統(tǒng)發(fā)育關(guān)系與地理分布, 本研究恢復(fù)長吻似鮑(Pseudogobio longirostris Mori, 1934, 模式地遼河水系遼寧凌源)給予似鮑譜系C分類名。本研究系統(tǒng)發(fā)育關(guān)系結(jié)果亦表明嚴格意義似鮑(=似鮑譜系A(chǔ)+B)、長吻似鮑(=似鮑譜系C)和似鮑譜系D+E均不是姊妹群關(guān)系,暗示似鮑譜系D+E屬于不同種。分子證據(jù)和形態(tài)證據(jù)進一步確認分布于長江的似鮑譜系D與分布于甌江的似鮑譜系E各自是不同種, 并支持平江似鮑是有效種, 它們一起形成一單系群。

圖 3 基于Genpofad距離重建的似鮑屬魚類鄰接網(wǎng)絡(luò)圖Fig. 3 NeighborNet network of Pseudogobio fishes based on Genpofad distances似鮑譜系定義依據(jù)圖 2Pseudogobio vaillanti Clades defined as in Fig. 2

圖 4 基于核基因多位點的似鮑屬魚類*BEAST物種樹Fig. 4 BEAST species tree of Pseudogobio fishes based on nuclear loci節(jié)點上方為后驗概率值, 外類群沒有顯示。似鮑譜系定義依據(jù)圖 2Above nodes are posterior probabilities, and outgroup taxa are not shown. Pseudogobio vaillanti Clades defined as in Fig. 2

楊秀平等[6]基于形態(tài)可量性狀的主成分分析結(jié)果認為桂林似鮑是似鮑的同物異名。本研究基于形態(tài)可量性狀的主成分分析結(jié)果亦不能分開桂林似鮑和嚴格意義似鮑, 但它們一起能很好地與長吻似鮑、似鮑譜系D及E分開。然而, 系統(tǒng)發(fā)育關(guān)系結(jié)果表明桂林似鮑與嚴格意義似鮑+長吻似鮑+扁嘴似鮑譜系A(chǔ)形成姊妹群關(guān)系, 暗示桂林似鮑是有效種。分子的物種界定分析結(jié)果確認桂林似鮑是有效種。側(cè)線鱗的數(shù)目亦能夠很好地區(qū)分桂林似鮑(43)和我國分布的似鮑屬魚類其他種(40—42)。因此, 結(jié)合分子證據(jù)、形態(tài)證據(jù)與系統(tǒng)發(fā)育關(guān)系, 本研究支持桂林似鮑是有效種。

3.2 似鮑屬魚類的系統(tǒng)發(fā)育關(guān)系與生物地理過程推測

夏曦中等[13]采納楊秀平等[6]基于形態(tài)可量性狀的主成分分析結(jié)果認為我國似鮑屬魚類僅似鮑單種, 并利用線粒體細胞色素b基因開展了我國廣義似鮑的譜系地理研究。隨后Tominaga等[17,18]使用線粒體細胞色素b基因?qū)Ψ植加谌毡镜谋庾焖契U進行了譜系地理研究。他們的研究均發(fā)現(xiàn)似鮑屬魚類包括多個遺傳譜系, 但都沒有清晰解決似鮑屬魚類譜系間的系統(tǒng)發(fā)育關(guān)系。本研究利用核基因多位點數(shù)據(jù)首次清楚地解析了似鮑屬魚類的系統(tǒng)發(fā)育關(guān)系, 揭示平江似鮑是似鮑譜系D+E的姊妹群;桂林似鮑是嚴格意義似鮑(=似鮑譜系A(chǔ)+B)+長吻似鮑(=似鮑譜系C)+扁嘴似鮑譜系A(chǔ)的姊妹群; 所有譜系一起與扁嘴似鮑譜系B形成姊妹群關(guān)系。依據(jù)Tominaga等[18]的廣義扁嘴似鮑譜系地理研究,扁嘴似鮑譜系A(chǔ)與B均為日本本島特有。依據(jù)本研究對我國似鮑屬魚類的分類厘定, 平江似鮑與桂林似鮑分布于珠江水系, 嚴格意義似鮑分布于韓江至長江的水系, 長吻似鮑分布于遼河、大洋河與鴨綠江水系, 似鮑譜系D分布于長江, 似鮑譜系E分布于甌江。結(jié)合上述系統(tǒng)發(fā)育關(guān)系與譜系的地理分布, 本研究提出以下生物地理過程是似鮑屬魚類物種分化的可能原因: (1)扁嘴似鮑譜系B位于系統(tǒng)發(fā)育樹的基部位置, 暗示似鮑屬魚類可能起源于東亞北部地區(qū), 扁嘴似鮑譜系B的物種形成可能是地理隔離造成; (2)桂林似鮑、嚴格意義似鮑、長吻似鮑與扁嘴似鮑譜系A(chǔ)一起構(gòu)成一單系群, 四者地理分布不重疊, 推測四者的物種形成可能是地理隔離造成。同時, 扁嘴似鮑譜系A(chǔ)是長吻似鮑的姊妹群, 嵌套在似鮑屬魚類內(nèi)部, 暗示扁嘴似鮑譜系A(chǔ)的物種形成亦可能是擴散過程造成, 即從大陸擴散到日本后的地理隔離; (3)平江似鮑、似鮑譜系D與E一起構(gòu)成一單系群, 三者地理分布不重疊, 推測三者的物種形成可能是地理隔離造成; (4)平江似鮑與桂林似鮑均分布于珠江水系, 且兩者不是姊妹群, 其系統(tǒng)位置均嵌套在似鮑屬魚類內(nèi)部, 推測平江似鮑與桂林似鮑的物種形成亦可能是擴散過程造成, 即似鮑屬魚類祖先通過兩次擴散進入珠江水系, 隨后的地理隔離造成。

表 3 似鮑譜系A(chǔ)至E、桂林似鮑和平江似鮑的形態(tài)特征Tab. 3 Morphological characters of Pseudogobio vaillanti Clades A-E, Pseudogobio guilinensis and Pseudogobio banggiangensis

圖 5 似鮑譜系A(chǔ)至E和桂林似鮑的第一與第二主成分得分散點圖Fig. 5 Scatter plots of scores on the first and second principal components for Pseudogobio vaillanti Clades A-E and P. guilinensis

綜上所述, 本研究清楚地對我國似鮑屬魚類進行了物種界定, 解析了似鮑屬魚類的系統(tǒng)發(fā)育關(guān)系。進一步的研究應(yīng)基于動物命名法則給予似鮑譜系D與E相應(yīng)的分類名稱, 結(jié)合形態(tài)數(shù)據(jù)厘清廣義扁嘴似鮑的分類。同時, 運用分子鐘技術(shù)與祖先重建技術(shù)測試本研究提出的似鮑屬魚類物種分化的生物地理過程, 為理解東亞地區(qū)淡水魚類多樣性空間格局形成過程提供案例。

[1]Bǎnǎrescu P, Nalbant T T. Pisces, Teleostei, Cyprinidae (Gobioninae) [M]. Das Tierreich, Lieferung 93, Berlin: Walter de Guryter. 1973, 1—304

[2]Yue P Q. Pseudogobio Bleeker, 1860 [A]. In: Chen Y Y (Eds.), Fauna Sinica, Osteichthyes, Cypriniformes (II) [C]. Beijing: Science Press. 1998, 374—378 [樂佩琦. 似鮑屬. 見: 陳宜瑜, 中國動物志硬骨魚綱鯉形目(中卷).北京: 科學(xué)出版社. 1998, 374—378]

[3]Xie C X, Liu Q D. A Preliminary study on the biology of Pseudogobio vaillanti vaillanti (Sauvage) [J]. Journal of Huazhong Agricultural University, 1986, 5(1): 73—82 [謝從新, 劉齊德. 似鮑生物學(xué)的初步研究. 華中農(nóng)業(yè)大學(xué)學(xué)報, 1986, 5(1): 73—82]

[4]Eschmeyer W N, Fricke R, Van der Laan R. Catalog of Fishes: Genera, Species, References, http://researcharchive.calacademy.org/research/ichthyology/catalog/fishcatmain.asp

[5]Luo Y L, Yue P Q, Chen Y Y. Pseudogobio Bleeker, 1860 [A]. In: Wu X W (Eds.), Fauna of Cyprinidae in China (II) [C]. Shanghai: Shanghai People’s Publishing House. 1977, 512—516 [羅云林, 樂佩琦, 陳宜瑜. 似鮑屬. 見: 伍獻文, 中國鯉科魚類志(下卷). 上海: 上海人民出版社. 1977, 512—516]

[6]Yang X P, Zhang M Y, Liu H Z. Morphological variations and geographical differentiation of the genus Pseudogobio (Teleostei: Cyprinidae) in China [J]. Acta Hydrobiologica Sinica, 2002, 26(3): 281—285 [楊秀平,張敏瑩, 劉煥章. 中國似鮑屬魚類的形態(tài)變異及地理分化研究. 水生生物學(xué)報, 2002, 26(3): 281—285]

[7]Hosoya K. Interrelationships of the Gobioninae (Cyprinidae) [A]. In: Uyeno T, Arai T, Taniuchi T, et al. (Eds.), Indo-Pacific Fish Biology: Proceeding of the Second International Conference on Indo-Pacific Fishes [C]. Tokyo: The Ichthyological Society of Japan. 1986, 484—501

[8]Yu L R, Yue P Q. Studies on phylogeny of the Pseudogobiini fishes (Cypriniformes: Cyprinidae) [J]. Acta Zootaxonomica Sinica, 1996, 21(2): 244—255 [俞利榮, 樂佩琦. 似鮑類魚類的系統(tǒng)發(fā)育研究(鯉形目: 鯉科). 動物分類學(xué)報, 1996, 21(2): 244—255]

[9]Yang J Q, He S P, Freyhof J, et al. The phylogenetic relationships of the Gobioninae (Teleostei: Cyprinidae) inferred from mitochondrial cytochrome b gene sequences [J]. Hydrobiologia, 2006, 553(1): 255—266

[10]Liu H Z, Yang J Q, Tang Q Y. Estimated evolutionary tempo of East Asian gobionid fishes (Teleostei: Cyprinidae) from mitochondrial DNA sequence data [J]. Hydrobiology, 2010, 55(15): 1502—1510

[11]Tang K L, Agnew M K, Chen W J, et al. Phylogeny of the gudgeons (Teleostei: Cyprinidae: Gobioninae) [J]. Molecular Phylogenetics and Evolution, 2011, 61(1): 103—124

[12]Jiang Z G, Zhang E. Molecular evidence for taxonomic status of the gudgeon genus Huigobio Fang, 1938 (Teleostei: Cypriniformes), with a description of a new species from Guangdong Province, South China [J]. Zootaxa, 2013, 3731(1): 171—182

[13]Xia X Z, Zhang X P, Liu H Z. Molecular biological geography analysis of Pseudogobio vaillanti [J]. Progress in Natural Science, 2005, 15(9): 1064—1069 [夏曦中, 張訓(xùn)蒲, 劉煥章. 似鮑(Pseudogobio vaillanti)的分子生物地理學(xué)分析. 自然科學(xué)進展, 2005, 15(9): 1064—1069]

[14]Ward R D. DNA barcode divergence among species and genera of birds and fishes [J]. Molecular Ecology Resources, 2009, 9(4): 1077—1085

[15]Guo L, Li J, Wang Z S, et al. Phylogenetic relationships of noodle-fishes (Osmeriformes: Salangidae) based on four mitochondrial genes [J]. Acta Hydrobiologica Sinica, 2011, 35(3): 449—459 [郭立, 李雋, 王忠鎖, et al. 基于四個線粒體基因片段的銀魚科魚類系統(tǒng)發(fā)育. 水生生物學(xué)報, 2011, 35(3): 449—459]

[16]Chen H, Polgar G, Fu C Z, et al. Cryptic species and evolutionary history of Boleophthalmus pectinirostris complex along the northwestern Pacific coast [J]. Acta Hydrobiologica Sinica, 2014, 38(1): 75—86 [陳卉, Polgar G, 傅萃長, 等. 西北太平洋海岸帶大彈涂魚復(fù)合體的隱存種與進化歷史. 水生生物學(xué)報, 2014, 38(1): 75—86]

[17]Tominaga K, Watanabe K, Kakioka R, et al. Two highly divergent mitochondrial DNA lineages within Pseudogobio esocinus populations in central Honshu, Japan [J]. Ichthyological Research, 2009, 56(2): 195—199

[18]Tominaga K, Nakajima J, Watanabe K. Cryptic divergence and phylogeography of the pike gudgeon Pseudogobio esocinus (Teleostei: Cyprinidae): a comprehensive case of freshwater phylogeography in Japan [J]. Ichthyological Research, 2016, 63(1): 79—93

[19]Hedin M. High-stakes species delimitation in eyeless cave spiders (Cicurina, Dictynidae, Araneae) from central Texas [J]. Molecular Ecology, 2015, 24(2): 346—361

[20]Karanovic T, Djurakic M, Eberhard S M. Cryptic species or inadequate taxonomy? implementation of 2D geometric morphometrics based on integumental organs as landmarks for delimitation and description of copepod taxa [J]. Systematic Biology, 2016, 65(2): 304—327

[21]Ojanguren-Affilastro A A, Mattoni C I, Ochoa J A, et al. Phylogeny, species delimitation and convergence in the South American bothriurid scorpion genus Brachistosternus Pocock 1893: Integrating morphology, nuclear and mitochondrial DNA [J]. Molecular Phylogenetics and Evolution, 2016, 94(2016): 159—170

[22]Miller S A, Dykes D D, Polesky H F. A simple salting out procedure for extracting DNA from human nucleated cells [J]. Nucleic Acids Research, 1988, 16(3): 12—15

[23]Kim K Y, Bang I C. Molecular phylogenetic position of Abbottina springeri (Cypriniformes, Cyprindae) based on nucleotide sequences of RAG1 gene [J]. Korean Journal of Ichthyology, 2010, 22(4): 273—278

[24]Chow S, Hazama K. Universal PCR primers for S7 ribosomal protein gene introns in fish [J]. Molecular Ecology, 1998, 7(9): 1255—1256

[25]Xia X H. DAMBE5: A comprehensive software package for data analysis in molecular biology and evolution [J]. Molecular Biology and Evolution, 2013, 30(7): 1720—1728

[26]Katoh K, Standley D M. MAFFT Multiple Sequence Alignment Software Version 7: Improvements in performance and usability [J]. Molecular Biology and Evolution, 2013, 30(4): 772—780

[27]Tamura K, Stecher G, Peterson D, et al. MEGA6: Molecular evolutionary genetics analysis version 6.0 [J]. Molecular Biology and Evolution, 2013, 30(12): 2725—2729

[28]Darriba D, Taboada G L, Doallo R, et al. jModelTest 2: more models, new heuristics and parallel computing [J]. Nature Methods, 2012, 9(9): 772

[29]Ronquist F, Teslenko M, Van der Mark P, et al. MrBayes 3.2: Efficient Bayesian phylogenetic inference and model choice across a large model space [J]. Systematic Biology, 2012, 61(3): 539—542

[30]Joly S, Bruneau A. Incorporating Allelic Variation for Reconstructing the Evolutionary History of Organisms from Multiple Genes: An Example from Rosa in North America [J]. Systematic Biology, 2006, 55(4): 623—636

[31]Huson D H, Bryant D. Application of phylogenetic networks in evolutionary studies [J]. Molecular Biology and Evolution, 2006, 23(2): 254—267

[32]Huelsenbeck J P, Andolfatto P. Inference of population structure under a Dirichlet process model [J]. Genetics, 2007, 175(4): 1787—1802

[33]Yang Z H, Rannala B. Bayesian species delimitation using multilocus sequence data [J]. Proceedings of the National Academy of Sciences, 2010, 107(20): 9264—9269

[34]Yang Z H, Rannala B. Unguided species delimitation using DNA sequence data from multiple loci [J]. MolecularBiology and Evolution, 2014, 31(12): 3125—3135

[35]Bouckaert R, Heled J, Kuhnert D, et al. BEAST 2: A software platform for bayesian evolutionary analysis [J]. PLoS Computational Biology, 2014, 10(4): 1—6

[36]Rambaut A, Drummond A. Tracer v1. 5. 2007. Available at: http://tree.bio.ed.ac.uk/software/tracer.

[37]Chen Y Y. Fauna Sinica, Osteichthyes, Cypriniformes (II) [M]. Beijing: Science Press. 1998, 16—18 [陳宜瑜. 中國動物志硬骨魚綱鯉形目(中卷). 北京: 科學(xué)出版社. 1998, 16—18]

[38]Kottelat M, Freyhof J. Handbook of European freshwater fishes [M]. Delémont: Imprimerie du Démocrate SA. 2007, 15—16

[39]Hammer ?, Harper D A T, Ryan P D. Past: Paleontological statistics software package for education and data analysis [J]. Palaeontologia Electronica, 2001, 4(1): 9

[40]Burnaby T P. Growth-invariant discriminant functions and generalized distances [J]. Biometrics, 1966, 22(1): 96—110

[41]Rolf F J, Bookstein F L. A comment on shearing as a method for ‘size correlation’ [J]. Systematic Zoology, 1987, 36(4): 356—367

[42]Waters J M, Rowe D L, Burridge C P. Gene trees versus species trees: reassessing life-history evolution in a freshwater fish radiation [J]. Systematic Biology, 2010, 59(5): 504—517

[43]Fu C Z, Guo L, Xia R. A multilocus phylogeny of Asian noodlefishes Salangidae (Teleostei: Osmeriformes) with a revised classification of the family [J]. Molecular Phylogenetics and Evolution, 2012, 62(3): 848—855

[44]Xia R, Durand J D, Fu C Z. Multilocus resolution of Mugilidae phylogeny (Teleostei: Mugiliformes): Implications for the family’s taxonomy [J]. Molecular Phylogenetics and Evolution, 2016, 96: 161—177

[45]Sauvage H E. Rhinogobio Vaillanti, n. sp [A]. In: Bulletin de la Societe philomathique de Paris [C]. Paris: Societe philomathique de Paris. 1878, 87

SPECIES DELIMITATION AND PHYLOGENETIC RELATIONSHIPS OF PSEUDOGOBIO FISHES IN THE SUBFAMILY GOBIONINAE

CAO Kai, LI Yu-Huo, HE Huan, LI Ming-Yue and FU Cui-Zhang
(Ministry of Education Key Laboratory for Biodiversity Science and Ecological Engineering, Institute of Biodiversity Science and School of Life Science, Fudan University, Shanghai 200438, China)

Species delimitation and phylogenetic relationships of the genus Pseudogobio fishes in the subfamily Gobioninae remain debated and unresolved. This study used molecular species delimitation methods and morphological analyses to reconstruct phylogenetic relationships of Pseudogobio fishes based on nuclear loci and the classification of this genus in China. The results of Bayesian tree showed that Pseudogobio guilinensis and P. banggiangensis were monophyletic specie, and that P. vaillanti and P. escocinus were polyphyletic species that the former included five Clades AE and the latter was comprised of two Clades A and B. The analyses of POFAD and Structurama showed that each clade of P. vaillanti and P. escocinus were independent genetic populations, and BP & P analysis strongly supported them as different species. The results of *Beast species tree revealed that P. esocinus Clade B occupied a basal position of Pseudogobio fishes; Clades A and B of P. vaillanti were the sister taxon relationship, and P. vaillanti Clade C was the sister taxon of P. esocinus Clade A, and all of them together were the sister group of P. guilinensis; Clades D and E of P. vaillanti were the sister taxon relationship, which were the sister group of P. banggiangensis. A revised classification of Pseudogobio fishes in China by molecular and morphological evidence as follows: P. vaillanti sensu stricto was limited to P. vaillanti Clades A and B; A resurrection of P. longirostris Mori, 1934 was assigned to P. vaillanti Clade C; Clades D and E of P. vaillanti were cryptic species, and P. guilinensis and P. banggiangensis were valid species, respectively.

Cyprinidae; Gobioninae; Pseudogobio; Species delimitation; Phylogeny

Q959.46

A

1000-3207(2017)03-0617-12

10.7541/2017.79

2016-04-27;

2016-08-21

國家自然科學(xué)基金(31471982)資助 [Supported by the National Natural Science Foundation of China (31471982)]

曹愷(1991—), 女, 浙江衢州人; 碩士研究生; 主要從事魚類生態(tài)與進化研究。E-mail: 13210700119@fudan.edu.cn

傅萃長, 教授, 博士生導(dǎo)師; E-mail: czfu@fudan.edu.cn