高陽,李海濤*,姜重臣,陳志云,李宏俊
( 1. 國家海洋局南海環(huán)境監(jiān)測中心,廣東 廣州 510300;2. 中國科學(xué)院南海海洋研究所 熱帶海洋生物資源與生態(tài)重點(diǎn)實(shí)驗(yàn)室,廣東 廣州 510301;3. 國家海洋環(huán)境監(jiān)測中心,遼寧 大連 116023)
海洋浮游軟體動(dòng)物是一類終生營浮游生活的類群,一般屬于腹足綱(Gastropoda)。據(jù)估計(jì),終生浮游腹足類僅有約140 種[1]。浮游腹足類大多棲息于上層200 m 水體,并有明顯的晝夜垂直遷移現(xiàn)象[2–3]。作為初級(jí)消費(fèi)者、被捕食者及鈣化生物,浮游腹足類在海洋生態(tài)系統(tǒng)的生物地球化學(xué)過程中扮演著十分關(guān)鍵的作用[4–5]。此外,由于對(duì)海洋酸化的極度敏感性,浮游腹足類成為良好的指示物種[6–7]。
自半個(gè)多世紀(jì)前張福綏[8]開展中國海的浮游軟體動(dòng)物分類以來,迄今鮮見其他分類報(bào)道,其研究結(jié)果影響深遠(yuǎn)。由于研究的匱乏,我國浮游腹足類的分類嚴(yán)重滯后,一些分類錯(cuò)誤延續(xù)至今,無法得到訂正,如:歸在螔螺科(Limacinidae)下的強(qiáng)卷螺屬(Agadina)種類實(shí)際上是成體營底棲生活的輪螺科(Architectonicidae) 的面盤幼蟲[9–10];厚唇螺屬(Diacria) 由D.trispinosa 復(fù)合群和D. quadridentata 復(fù)合群組成,其中D. trispinosa 復(fù)合群中的種類根據(jù)殼面花紋的分布來區(qū)分[11],我國記錄的D. trispinosa 缺乏殼面花紋分布這一關(guān)鍵分類特征的描述,實(shí)際上無法確定其為何種[8]。此外,一些物種的形態(tài)描述和圖示與實(shí)物有一定出入;基于形態(tài)變化的考慮而提出的一些亞種或變種, 尚無后續(xù)研究的證實(shí)。
大多數(shù)浮游腹足類的貝殼缺乏雕刻或生長紋等特征,其形態(tài)鑒定主要依據(jù)一些簡單的分類性狀,單純依靠形態(tài)特征很難進(jìn)行準(zhǔn)確的種類鑒定?;趍t-COI 基因的DNA 條形碼技術(shù)[12]目前已被廣泛應(yīng)用于浮游腹足類的種類鑒定和多樣性分析[13–15]。浮游軟體動(dòng)物分布廣泛,很多種類為接近全球分布或環(huán)球分布[13,16],是開展生物地理學(xué)或譜系地理學(xué)研究的良好材料[16–17]。本研究以長角螺屬(Clio)物種為材料,通過對(duì)mtCOI 基因和核18S rRNA 基因序列的分析,結(jié)合數(shù)據(jù)庫中已有的數(shù)據(jù),對(duì)其進(jìn)行了分類學(xué)和譜系地理學(xué)研究,以期了解該屬種類的地理遺傳分化和分布狀況,同時(shí)為其他浮游軟體動(dòng)物的多樣性研究提供參考。
實(shí)驗(yàn)樣品于2016 年5 月至2019 年1 月采自西北太平洋(含南海、呂宋海峽及其周邊海域和新幾內(nèi)亞島以北海域) 及北印度洋海域。使用大型浮游生物網(wǎng)(網(wǎng)目0.505 mm)采用垂直拖網(wǎng)和水平拖網(wǎng)相結(jié)合的方式采集浮游動(dòng)物樣品。為增加樣品的采獲幾率,西北太平洋海域的采樣均在夜間進(jìn)行。浮游動(dòng)物樣品經(jīng)無水酒精固定后冰凍保存帶回實(shí)驗(yàn)室,體視鏡下分選出長角螺屬動(dòng)物樣品,對(duì)其進(jìn)行初步的形態(tài)鑒定,并編號(hào)拍照。除呂宋海峽的樣品為部分測序外,其余采樣點(diǎn)的樣品全部用于測序,共計(jì)55 個(gè)標(biāo)本用于DNA 序列的分析,測序標(biāo)本的采樣信息見表1。
剪取長角螺游泳鰭,采用堿裂解法[18]提取DNA。使用引物L(fēng)CO-1490(5′-GGT CAA CAA ATC ATA AAG ATA TTG G-3′)和HCO-2198 (5′-TAA ACT TCA GGG TGA CCA AAA AAT CA-3′)[19]擴(kuò)增mt-COI 基因片段;使用引物18A1 (5′- CTG GTT GAT CCT GCC AGT CAT ATG C-3′)和1800 (5′-GAT CCT TCC GCA GGT TCA CCT ACG-3′)[20]擴(kuò)增18S rRNA 基因。
PCR 反應(yīng)體系總體積為50 μL,包括:2×buffer 25 μL,dNTP 400 μmol/L,引物各0.3 μmol/L,KOD FX 高保真PCR 酶(TOYOBO 公司)1.0 U,模板DNA1 μL,加去離子水補(bǔ)足至50 μL。PCR 反應(yīng)條件為:94℃預(yù)變性2 min;98℃變性15 s,52℃退火30 s,68℃延伸60 s (mtCOI)或120 s (18S rRNA),共30 個(gè)循環(huán);最后68℃延伸10 min。PCR 產(chǎn)物經(jīng)1%瓊脂糖凝膠電泳檢測后,送廣州艾基生物技術(shù)有限公司進(jìn)行回收純化和測序,測序引物與擴(kuò)增引物相同。為保證18S rRNA基因片段測通,中間引物1155R (5′-CCG TCA ATTCCT TTA AGT TTC AG-3′)也被用于測序[20]。
表 1 測序標(biāo)本的采樣信息Table 1 Collection information for specimens analyzed in this study
拼接、校對(duì)后的測序結(jié)果,連同GenBank 中下載的序列,在BioEdit 軟件中進(jìn)行多重序列比對(duì)。使用MEGA 7.0 計(jì)算堿基的組成,基于Kimura 2-Parameter(K2P)參數(shù)模型計(jì)算遺傳距離,并采用鄰接法(Neighbor-joining,NJ)構(gòu)建系統(tǒng)進(jìn)化樹,可靠性經(jīng)過1 000 次自展(Bootstrap, BS)檢驗(yàn)。采用MrBayes 3.2[21]構(gòu)建貝葉斯樹(Bayesian Inference, BI),分子進(jìn)化模型使用jModelTest 2.1[22]根據(jù)Bayesian Information Criteria(BIC)標(biāo)準(zhǔn)計(jì)算得出。mtCOI 和18S rRNA 基因的最佳進(jìn)化模型分別為HKY+I+G 和K80+I。在馬爾可夫鏈蒙特卡羅(Markov Chain Monte Carlo,MCMC)參數(shù)下運(yùn)行300 萬代,每100 代抽樣一次,前2 500 代作為老化樣本舍去。BI 樹各分支置信度以后驗(yàn)概率(Posterior Probabilities, PP)表示。使用FigTree 1.4.3 顯示BI 樹的拓?fù)浣Y(jié)構(gòu)圖。
從GenBank 中下載的長角螺屬COI 基因序列包括C. pyramidata(MH101647~55、FJ876874~75、MF048921~22、KC774065~68、KP292791、KC788279、GQ861823)、C. pyramidata lanceolate (FJ876872~73、FJ876877~79)、C. pyramidata antarctica (FJ876876、MF048919)、C. cuspidata (FJ876869~71、KP292789、KC774064)、C. convexa (KC774062~63、KC774069)、C.recurva (FJ876880~86、KP292790、MF048923)、Clio sp.(FJ876887)。
采用GMYC(Generalized Mixed Yule Coalescent)模型[23]、ABGD(Automatic Barcoding Gap Discovery)方法[24]和3%K2P 遺傳距離閾值進(jìn)行分子可操作分類單元(Molecular Operational Taxonomic Units, MOTUs)估計(jì)。
GMYC 模型:利用BEAST v1.10.1[25]構(gòu)建超度量系統(tǒng)發(fā)育樹;采用uncorrelated log-normal clock 模型,1×108次MCMC 分析,1 000 次burn in;利用BEAST 軟件包中的TreeAnnotator 分析超度量系統(tǒng)發(fā)育樹,并生成最大譜系置信樹;隨后用FigTree 1.4.3 將最大譜系置信樹轉(zhuǎn)化為APE 識(shí)別的格式(Newick);最后利用R 環(huán)境中的“SPLITS”軟件包進(jìn)行單閾值GMYC 物種界定[23],利用gmyc 函數(shù)獲得零模型與GMYC 模型之間的似然比值、種?種群的轉(zhuǎn)折點(diǎn),以及轉(zhuǎn)折點(diǎn)對(duì)應(yīng)的MOTU 的數(shù)量。
ABGD 方法:樣本在線提交至ABGD 網(wǎng)站(http://wwwabi.snv.Jussieu.fr./pubilc/abgd/abgdweb.html),采用K2P 模型計(jì)算遺傳距離,其他參數(shù)均使用網(wǎng)絡(luò)版本的默認(rèn)值,基于遺傳距離對(duì)樣本進(jìn)行劃分。
3%K2P 遺傳距離閾值:將K2P 遺傳距離小于或等于3%(0.03)的序列劃分為同一MOTU,遺傳距離大于3%的序列劃分為不同的MOTU。
西北太平洋和印度洋海域采集的長角螺屬標(biāo)本,經(jīng)形態(tài)鑒定,屬于矛頭長角螺(C. pyramidata)(圖1a–c)、尖棘長角螺(C. cuspidata) (圖1d–f)和膨凸長角螺(C. convexa) (圖1g)3 個(gè)形態(tài)種,以及個(gè)別通過形態(tài)難以定種的幼體(圖1h)。在采樣過程中,較大個(gè)體的標(biāo)本極易破損,棘刺發(fā)達(dá)的尖棘長角螺很難得到完整的標(biāo)本。呂宋海峽采集到的尖棘長角螺(圖1d–e)明顯較新幾內(nèi)亞島以北海域采集到的個(gè)體(圖1f)寬大,后者殼體明顯纖細(xì)。
測序共獲得55 條mtCOI 序列(GenBank 登錄號(hào)MK749610~64,表1) 和9 條18S rRNA 基因序列(GenBank 登錄號(hào)MK749665~73,表1),其中55 條mt-COI 序列定義了51 種單倍型。mtCOI 序列長658 bp,A、T、G、C 含量分別為22.4%、44.3%、19.0% 和14.3%,A+T 的含量明顯高于G+C 的含量。mtCOI 序列共有變異位點(diǎn)154 個(gè),簡約信息位點(diǎn)113 個(gè)。
對(duì)18S rRNA 基因序列兩端做截齊處理后,序列長1 728~1 740 bp,存在兩段分別為8 bp(TTTCGGGG)和5 bp(TGAAG)的插入/缺失片段,以及1 個(gè)堿基(T)的插入/缺失位點(diǎn)。A、T、G、C 含量分別為24.0%、25.2%、27.7%和23.2%,A+T 和G+C 的含量大體相當(dāng)。3 個(gè)形態(tài)種均分別測定了3 個(gè)個(gè)體的18S rRNA 基因序列,其中尖棘長角螺和膨凸長角螺在種內(nèi)沒有變異位點(diǎn),矛頭長角螺在種內(nèi)有3 個(gè)變異位點(diǎn)。
mtCOI 和核18S rRNA 基因的BI 樹和NJ 樹有大體一致的拓?fù)浣Y(jié)構(gòu),本文只顯示BI 樹(圖2,圖3)。mtCOI 基因系統(tǒng)分析顯示,矛頭長角螺和尖棘長角螺存在明顯的譜系地理結(jié)構(gòu),分別形成4 個(gè)明顯分化的支系;膨凸長角螺(含印度洋、紅海和西太平洋樣品)和曲形長角螺(C. recurva,含南大西洋、北大西洋和南大洋樣品)形成單系群,無地理遺傳分化。矛頭長角螺的4 個(gè)支系(A–D)均獲得較高的支持率(PP≥95%,BS≥90%),本文測序的標(biāo)本分屬2 個(gè)不同的支系(A 和C)。支系A(chǔ) 基本上呈全球分布,支系B 目前僅見于大西洋(含加勒比海)和印度洋,支系C 目前僅見于西太平洋和東印度洋,支系D(C. pyramidata antarctica)目前僅見于南大洋的高緯度地區(qū)(圖4)。
尖棘長角螺的4 個(gè)支系(E–H)也獲得較高的支持率(PP≥99%,BS>85%),本文測序的標(biāo)本分屬支系E 和F。不同支系的地域性十分明顯,支系E 分布于呂宋海峽,支系F 分布于西太平洋和印度洋,支系G分布于大西洋,支系H 分布于南大洋(圖5)。
18S rRNA 基因系統(tǒng)樹中(圖3),矛頭長角螺的支系A(chǔ)、B 和C 形成一個(gè)獨(dú)立的分支,與支系D 區(qū)分開來。膨凸長角螺和尖棘長角螺的個(gè)體形成單系群,尖棘長角螺的支系E 和F 在18S rRNA 基因系統(tǒng)樹中無法得到顯現(xiàn)。
K2P 遺傳距離分析表明(因C. pyramidata(GQ861823)、C. recurva(MF048923)和Clio sp.(FJ876887)可能存在種類誤定或序列太少,未參與分析),形態(tài)種的種內(nèi)遺傳距離在0~0.089 之間,種間遺傳距離在0.070~0.185 之間,兩者發(fā)生較大范圍重疊(表2)。矛頭長角螺支系內(nèi)的遺傳距離為0~0.026,支系間的遺傳距離為0.031~0.089;尖棘長角螺支系內(nèi)的遺傳距離為0~0.014,支系間的遺傳距離為0.048~0.079。
ABGD 比較穩(wěn)定的初始劃分結(jié)果與GMYC 和3%K2P 遺傳距離閾值的MOTUs 估計(jì)結(jié)果完全相同(圖2, GMYC 和ABGD 原始劃分結(jié)果未給出),矛頭長角螺和尖棘長角螺的4 個(gè)支系均被劃分為獨(dú)立的MOTU。
圖 1 長角螺屬種類的貝殼形態(tài)(比例尺=5 mm)Fig. 1 Shell morphology of Clio species (scale bars=5 mm)
由于長角螺外殼薄脆,在采樣過程中極易造成破損,相當(dāng)部分的標(biāo)本為非完整的個(gè)體,如尖棘長角螺發(fā)達(dá)的棘刺很難完整保留(圖1),給形態(tài)鑒定帶來很大的難度。外形上的相似性,也難以區(qū)分不同大小的個(gè)體是種間差異還是同種不同發(fā)育階段的個(gè)體,如矛頭長角螺支系C 中2 個(gè)不同大小的個(gè)體在形態(tài)上有較大的差異(圖2),支系A(chǔ) 和C 的成體在形態(tài)上則難以區(qū)分。張福綏[8]描述了中國近海長角螺屬3 個(gè)種(含變種),其中包括一新變種小尾長角螺(C. pyramidata var. microcaudata,原文Euclio pyramidata var. microcaudata),該變種僅采集到1 個(gè)個(gè)體。該變種能否被視為一個(gè)有效的分類單元,尚需進(jìn)一步的研究證實(shí)。我們檢視了采集到的大量長角螺標(biāo)本,未發(fā)現(xiàn)類似個(gè)體。新幾內(nèi)亞島以北海域采集到的尖棘長角螺與呂宋海峽分布的尖棘長角螺在形態(tài)上有一定的差異,且兩者的mtCOI 基因序列也有明顯分化。然而由于前者僅采集到一個(gè)不太完整的標(biāo)本,尚無法判斷這種形態(tài)差異是種間的差異還是同種個(gè)體間的差異。
相比膨凸長角螺和曲形長角螺兩個(gè)種的種內(nèi)遺傳距離,矛頭長角螺和尖棘長角螺兩個(gè)形態(tài)種各支系間的遺傳距離已經(jīng)達(dá)到種間的分化水平。矛頭長角螺的支系D(C. pyramidata antarctica)與其他3 個(gè)支系在18S rRNA 基因序列上也存在明顯的差異,并在系統(tǒng)樹上與其他支系區(qū)別開來。因此,C. pyramidata antarctica 這一亞種應(yīng)提升為種,該亞種(或變種)的形態(tài)[2,26]也明顯有別于典型的矛頭長角螺,前者明顯修長。18S rRNA 基因系統(tǒng)樹并不支持其他支系的劃分,這可能是由于18S rRNA 基因較為保守的緣故。矛頭長角螺和尖棘長角螺是否為隱存種還需要進(jìn)一步研究的證實(shí)(如ITS 序列)。
圖 2 基于mtCOI 基因構(gòu)建的長角螺屬貝葉斯系統(tǒng)發(fā)育樹Fig. 2 Bayesian inference phylogenetic tree of genus Clio based on mtCOI gene sequences
海洋生物,無論是終生浮游生物還是其他海洋生物的浮游幼體,可以借助洋流進(jìn)行擴(kuò)散,開放的海洋環(huán)境也沒有阻礙基因交流的屏障,因而通常認(rèn)為海洋生物會(huì)在很大的地理尺度內(nèi)表現(xiàn)出較低的遺傳分化[27]。目前,越來越多的證據(jù)改變了我們這種傳統(tǒng)的觀點(diǎn),即使在開放海域,海洋學(xué)和生態(tài)學(xué)過程也會(huì)對(duì)海洋生物的基因流產(chǎn)生阻礙作用,從而導(dǎo)致在不同尺度上產(chǎn)生遺傳分化[28–30]。
圖 3 基于18S rRNA 基因構(gòu)建的長角螺屬貝葉斯系統(tǒng)發(fā)育樹Fig. 3 Bayesian inference phylogenetic tree of genus Clio based on nuclear18S rRNA gene sequences
圖 4 矛頭長角螺4 個(gè)mtCOI 支系的地理分布Fig. 4 Geographical distributions of four mtCOI lineages for C. pyramidata
結(jié)合本文和數(shù)據(jù)庫中已有的數(shù)據(jù)分析表明,矛頭長角螺和尖棘長角螺存在較為明顯的地理群體遺傳分化,中國海及其鄰近海域可能只有矛頭長角螺支系A(chǔ) 的存在,分布于呂宋海峽的尖棘長角螺支系E為一新的譜系分支。從全球尺度來講,尖棘長角螺的采樣范圍和采樣數(shù)量較為有限,其譜系地理結(jié)構(gòu)還需更廣泛的采樣進(jìn)行研究。南海–呂宋海峽及其周邊海域與新幾內(nèi)亞島以北海域具有完全不同的長角螺種類和支系組成,除分布著矛頭長角螺和尖棘長角螺完全不同的支系外,新幾內(nèi)亞島以北海域分布的膨凸長角螺在南海–呂宋海峽及其周邊海域(10 余航次調(diào)查)未有發(fā)現(xiàn)。北赤道流到達(dá)菲律賓群島東岸后發(fā)生分叉,形成南、北流向的棉蘭老流和黑潮[31],相反的流向可能阻礙了兩海域間的物種擴(kuò)布及群體間的基因交流??梢酝茰y,不同遺傳背景的浮游軟體動(dòng)物群體可能對(duì)海洋酸化等環(huán)境變化產(chǎn)生不一樣的響應(yīng)[32]。
我們分析了印度洋和西北太平洋海域數(shù)百個(gè)浮游軟體動(dòng)物的mtCOI 基因序列和核18S、28S 基因序列,并通過與其他海域序列的比較,發(fā)現(xiàn)浮游軟體動(dòng)物中普遍存在mtCOI 基因的地理遺傳分化現(xiàn)象(結(jié)果另文發(fā)表)。如按照DNA 條形碼遺傳距離的10 倍法則[33]或3%閾值[12]標(biāo)準(zhǔn),這些具有地理遺傳分化的形態(tài)種會(huì)被認(rèn)為存在隱存種。
DNA 條形碼技術(shù)是一種基于種內(nèi)及種間DNA序列變異程度的距離法(Distance-Based Method),它要求種間遺傳距離明顯大于種內(nèi)遺傳距離,形成“條形碼間隙(Barcoding Gap)”[34–35]。然而,DNA 條形碼技術(shù)也受到一些質(zhì)疑和批判。如單獨(dú)的線粒體DNA 序列不足以識(shí)別物種,因?yàn)檫@種遺傳分化并非是劃定物種界限的必要條件[36–37];“條形碼間隙”是人為的不充分取樣造成的[38]。顯然,取樣范圍的擴(kuò)大(全球取樣)、測序樣本數(shù)量的增加會(huì)使種內(nèi)遺傳距離增大,測序物種數(shù)量的增加會(huì)使種間遺傳距離縮小。取樣數(shù)量和地理尺度的大小有時(shí)會(huì)顯著影響DNA 條形碼的分析結(jié)果,甚至可能得出錯(cuò)誤的結(jié)論。單純依據(jù)mtCOI 基因序列來揭示隱存種,證據(jù)可能不夠充分,還需要核基因證據(jù)的支持。
表 2 基于mtCOI 基因的長角螺屬4 個(gè)形態(tài)種的K2P 遺傳距離Table 2 The K2P distances of four Clio morphospecies based on mtCOI sequences
圖 5 尖棘長角螺4 個(gè)mtCOI 支系的地理分布Fig. 5 Geographical distributions of four mtCOI lineages for C. cuspidata
形態(tài)種內(nèi)的遺傳分化廣泛存在于海洋動(dòng)物中,除浮游軟體動(dòng)物[16]外,還有如水母[39–40]、甲殼類[41–43]、毛顎類[44]和魚類[45]等。這種遺傳分化不僅存在于大的地理尺度,也存在于小的區(qū)域尺度。
現(xiàn)有公共數(shù)據(jù)庫(如GenBank、BOLD 等)中的參考序列缺乏相對(duì)應(yīng)的物種形態(tài)照片,并且數(shù)據(jù)庫中還存在誤定的種類。來自不同海域的序列,即使與數(shù)據(jù)庫中已有的序列屬于同一形態(tài)種,由于不同地理群體的遺傳分化,仍然無法通過DNA 條形碼序列的比對(duì)來定種。如本文測序的呂宋海峽尖棘長角螺序列,無法在數(shù)據(jù)庫中與已有的其他海域的尖棘長角螺序列形成高精度匹配。因此,建立我國海域本地物種的DNA 條形碼數(shù)據(jù)庫是極為必要的,且測序個(gè)體應(yīng)有足夠的數(shù)量并覆蓋較大的地理尺度。
致謝:中國科學(xué)院南海海洋研究所周林濱博士采集了印度洋樣品,李開枝研究員在實(shí)驗(yàn)過程中給予了幫助和支持,謹(jǐn)致謝忱。