吳 娜，張 楠，曹 明，郭華陽，朱克誠(chéng)，楊靜文，江世貴，張殿昌

(1.中國(guó)水產(chǎn)科學(xué)研究院南海水產(chǎn)研究所，農(nóng)業(yè)部南海漁業(yè)資源開發(fā)利用重點(diǎn)實(shí)驗(yàn)室，廣東省漁業(yè)生態(tài)環(huán)境重點(diǎn)開放實(shí)驗(yàn)室，廣州 510300；2.上海海洋大學(xué)水產(chǎn)與生命學(xué)院，上海 201306；3.廣東省漁業(yè)種質(zhì)保護(hù)中心，廣州 510300)

魚卵和仔稚魚是魚類生活史的早期階段，魚卵和仔稚魚調(diào)查是研究魚類繁殖習(xí)性、生殖洄游路線、產(chǎn)卵場(chǎng)、種類分布和種群動(dòng)態(tài)的基礎(chǔ)，也可為魚類人工繁育技術(shù)和生理生化機(jī)制研究提供基礎(chǔ)數(shù)據(jù)[1]。魚卵和仔稚魚種類的準(zhǔn)確鑒別，是開展魚卵和仔稚魚研究的重要前提。通常，魚卵鑒別主要依據(jù)卵的形狀、大小、卵周隙、油球、內(nèi)膜、絨毛膜、卵黃形態(tài)和胚胎形狀等形態(tài)特征[2]，然而近緣物種的魚卵由于形態(tài)相似很難區(qū)分，大多種類只能鑒定到科或?qū)俚乃?，目前只有極少種類魚卵鑒定到種的水平，尋求穩(wěn)定、準(zhǔn)確的分類依據(jù)和鑒定方法，是魚類浮游生物研究亟待解決的關(guān)鍵問題。DNA條形碼(DNA barcoding)技術(shù)是一種利用標(biāo)準(zhǔn)基因區(qū)域序列多態(tài)性進(jìn)行物種鑒定的分子方法，條形碼序列分析不依賴于形態(tài)特征，對(duì)形態(tài)特征近似，形態(tài)學(xué)分類比較困難的疑難種類的分類以及形態(tài)學(xué)分類上存有疑問種類的驗(yàn)證提供了新的方法和手段[3-4]。魚類線粒體細(xì)胞色素C氧化酶亞基I(mtCOI)基因具有序列變異水平適中，變異區(qū)域兩端序列高度保守和序列長(zhǎng)度適宜的特點(diǎn)，是被廣泛接受的魚類DNA條形碼標(biāo)準(zhǔn)基因[5-7]，目前魚類DNA條形碼數(shù)據(jù)庫(Barcode of Life Data System [BOLD])中已有超過17.5萬件憑證標(biāo)本，代表了1.5萬種以上的魚類種類[8]。隨著魚類DNA條形碼數(shù)據(jù)庫的不斷完善，研究證實(shí)DNA條形碼技術(shù)在魚卵和仔稚魚的種類鑒定中也具有重要的應(yīng)用價(jià)值[9-11]。

南沙群島位于南海南部，屬于印度洋-太平洋熱帶動(dòng)物區(qū)系，是世界上海洋生物最豐富的區(qū)域之一。根據(jù)已有的調(diào)查資料顯示南海魚類有3 000多種，分別隸屬于3個(gè)綱、35個(gè)目、236個(gè)科和822屬[12]。目前南海魚類中只對(duì)極少數(shù)種類的魚卵和仔稚魚進(jìn)行過研究[4，13-14]，絕大多數(shù)種類的魚卵和仔稚魚的形態(tài)特征并不清楚，使得南海重要經(jīng)濟(jì)魚類產(chǎn)卵行為、產(chǎn)卵場(chǎng)判別和種群補(bǔ)充機(jī)制研究無法開展。永暑礁是一座珊瑚環(huán)礁，為珊瑚礁魚類和熱帶大洋性魚類的重要產(chǎn)卵場(chǎng)，然而目前未見有關(guān)永暑礁海域魚卵鑒定和分類的相關(guān)研究。因此，本研究擬以永暑礁瀉湖魚卵為研究對(duì)象，采用形態(tài)分類特征和DNA條形碼技術(shù)對(duì)采集的魚類進(jìn)行鑒定分析，以期了解珊瑚礁魚類和熱帶大洋性魚類的產(chǎn)卵季節(jié)和種類分布特征，為相關(guān)種類的種群補(bǔ)充機(jī)制研究提供基礎(chǔ)材料。

1 材料和方法

1.1 樣品采集和形態(tài)特征分析

依托中國(guó)水產(chǎn)科學(xué)研究院南海水產(chǎn)研究所“南鋒號(hào)”科學(xué)調(diào)查船，于2016年4月13號(hào)在中國(guó)南海永暑礁瀉湖海域按照《海洋調(diào)查規(guī)范》(GB12763.1-7-91)進(jìn)行魚卵樣品采集。采用孔徑為0.505 mm篩絹制成的大型浮游生物網(wǎng)(網(wǎng)口內(nèi)徑80 cm、長(zhǎng)270 cm、網(wǎng)口面積為0.5 m2)按照設(shè)定站位進(jìn)行表層水平和垂直拖網(wǎng)，采集魚卵樣品，-20 ℃冰箱中保存?zhèn)溆?。采集的樣品運(yùn)至實(shí)驗(yàn)室后，于解剖鏡下進(jìn)行形態(tài)特征觀察和顯微拍照，觀察記錄魚卵的卵膜(Chorion)、圍囊腔(Perivitelline space)、卵黃(Yolk)、油球(Oil globule)、色素分布(Pigmentation)和胚體發(fā)育(Embryo)等形態(tài)特征和魚卵卵徑參數(shù)，并對(duì)魚卵形態(tài)特征進(jìn)行了描述分析(表1)。

1.2 DNA提取和COI基因測(cè)序分析

魚卵樣品用雙蒸水沖洗，將樣品置于1.5 mL離心管中，采用Magen公司的試劑盒(HiPure Mollusc DNA Kit)提取基因組DNA，提取的DNA置于-20 ℃保存?zhèn)溆谩Ｒ镌O(shè)計(jì)參照Ward等[15],由廣州睿博生物公司合成，序列擴(kuò)增引物為Fish-F：5′-TCRACYAAYCAYAAAGAYATYGGCAC-3′和Fish-R：5′-ACTTCWGGGTGRCCRAAGAATCA-3′。PCR擴(kuò)增總體積為25 μL，反應(yīng)體系包括RNase-Free Water (16.7 μL)，10×loading buffer (Mg2+free)2.5 μL，DNTP(2.5 mmol/L)2 μL，MgCl2(2.5 μmol/L)1.5 μL，F(xiàn)ish-F(10 μmol/L)0.5 μL，F(xiàn)ish-R(10 μmol/L)0.5 μL，rTaq (5 U/μL)0.3 μL及DNA模板(20 ng/μg)1.0 μL。PCR反應(yīng)程序?yàn)椋?4 ℃預(yù)變性3 min，94 ℃變性30 s，55 ℃退火30 s，72 ℃延伸90 s，共35個(gè)循環(huán)，最后72 ℃延伸10 min。PCR產(chǎn)物進(jìn)行1.2%的瓊脂糖膠電泳檢測(cè)，選取擴(kuò)增結(jié)果較好的產(chǎn)物送廣州睿博興科股份有限公司測(cè)序分析。

由于魚卵樣品在提取DNA、PCR擴(kuò)增及測(cè)序的實(shí)驗(yàn)過程中，都有可能會(huì)出現(xiàn)結(jié)果不佳的情況，因此本研究以最終獲得的有效的序列個(gè)體來進(jìn)行分析。

1.3 分子鑒定分析

測(cè)序獲得的魚卵COI基因序列與BOLD系統(tǒng)數(shù)據(jù)庫(http://www.barcodinglife.org)比對(duì)分析進(jìn)行物種鑒定，在比對(duì)分析中，BOLD系統(tǒng)為每個(gè)樣品分配一個(gè)條形碼索引號(hào)(Barcode Index Number，BIN)，當(dāng)樣品COI基因序列與BIN中參考序列的相似性≥98.4%時(shí)，鑒定為種，當(dāng)相似性在95%～ 98.3%之間，鑒定到屬，當(dāng)相似性<95%時(shí)，為未知種[2，16]。利用MEGA6.0軟件根據(jù)配對(duì)遺傳距離采用Kimura two-parameter model (K2P)構(gòu)建鄰接進(jìn)化樹(Neighbor-joining tree)[17]，進(jìn)化樹末端節(jié)點(diǎn)代表了種類分組。

2 結(jié)果和分析

2.1 魚卵種類鑒定分析

采集的71個(gè)魚卵樣品根據(jù)外部形態(tài)特征初步劃分為圓形和橢圓形2種類型，圓形魚卵44個(gè)，橢圓形魚卵27個(gè)。經(jīng)DNA提取和線粒體COI基因片段測(cè)序分析，71個(gè)魚卵樣品中的30個(gè)個(gè)體獲得了有效的線粒體COI序列信息，測(cè)序成功率為42.3%。44個(gè)圓形魚卵中有27個(gè)測(cè)序成功，測(cè)序成功率達(dá)61.4%；27個(gè)橢圓形魚卵中有3個(gè)測(cè)序成果，測(cè)序成功率僅為11.1%。30個(gè)魚卵線粒體COI序列分別于BOLD數(shù)據(jù)庫比對(duì)分析，結(jié)果顯示與BOLD數(shù)據(jù)庫中參考序列相似性在99.6～100%之間，30個(gè)魚卵均能鑒定到種的水平，分屬于1個(gè)目，7個(gè)科，13個(gè)種(表1，表2)。刺尾魚科(Acanthuridae)有5個(gè)種，分別為白面刺尾魚(Acanthurusnigricans)、扁體櫛齒刺尾魚(Ctenochaetusstrigosus)、雙斑櫛齒刺尾魚(C.binotatus)、白尾櫛齒刺尾魚(C.flavicauda)和黑背鼻魚(Nasolituratus)，共有19個(gè)魚卵，占鑒定魚卵總數(shù)的63.3%。鸚哥魚科(Scaridae)有2個(gè)種，分別為史氏鸚哥魚(Scarusschlegeli)和斑點(diǎn)鯨鸚哥魚(Cetoscarusocellatus)。隆頭魚科(Labridae)有2個(gè)種，分別為單帶尖唇魚(Oxycheilinusunifasciatus)和雙線尖唇魚(O.digramma)。鲹科(Carangidae)有1個(gè)種為藍(lán)圓鲹(Decapterusmaruadsi)；籃子魚科(Siganidae)有1個(gè)種為銀籃子魚(Siganusargenteus)；笛鯛科(Lutjanidae)有1個(gè)種為欖色細(xì)齒笛鯛(Aphareusfurca)；鰈魚科(Zanclidae)有1種為鐮魚(Zancluscornutus)。雙斑櫛齒刺尾魚有10個(gè)魚卵，白尾櫛齒刺尾魚有4個(gè)魚卵，這兩個(gè)種為該海域的優(yōu)勢(shì)種類，其余各種類均為1～2個(gè)魚卵。

表1 13種魚卵的形態(tài)特征Tab.1 Morphological features of 13 species of fish eggs

續(xù)表1

表2 南海永暑礁海域魚卵mtCOI序列比對(duì)結(jié)果Tab.2 Sequences comparison results of mtCOI in fish eggs from Yongshu reef lagoon of South China Sea

2.2 系統(tǒng)進(jìn)化分析

獲得的30個(gè)魚卵的線粒體COI序列長(zhǎng)度在463～662 bp之間，序列中沒有插入/缺失位點(diǎn)和終止密碼子。根據(jù)線粒體COI序列計(jì)算種內(nèi)的K2P遺傳距離為0～0.002，同一個(gè)屬種間的遺傳距離為0.004～0.132，而同一個(gè)科不同屬間的遺傳距離為0.117～0.185。根據(jù)K2P模型利用遺傳距離構(gòu)建的鄰接系統(tǒng)進(jìn)化樹顯示，相同的種類聚為一個(gè)分支；刺尾魚科中的刺尾魚屬和櫛齒刺尾魚屬聚為一支，而與同一科的鼻魚屬未形成單系支(圖1)。

圖1 根據(jù)K2P模型利用遺傳距離構(gòu)建的鄰接系統(tǒng)進(jìn)化樹Fig.1 Neighbor-joining tree constructed based on K2P model

2.3 魚卵形態(tài)特征分析

根據(jù)魚卵的形狀、卵膜、圍囊腔、卵黃、油球、色素分布和胚胎發(fā)育情況等形態(tài)性狀，詳細(xì)描述了該研究鑒定的13種魚類的魚卵形態(tài)特征(表1)。史氏鸚哥魚和藍(lán)圓鲹的魚卵均為橢圓形，卵膜光滑，無色透明，圍囊腔較寬，但是史氏鸚哥魚卵徑較大，卵黃均勻分布，而藍(lán)圓鲹的卵黃位于魚卵的一端，并且體積較小。其余11種魚類的魚卵均為圓形，魚卵直徑在1.34～1.99 mm之間，黑背鼻魚的魚卵直徑最小，鐮魚的魚卵直徑最大。已有的研究表明海水魚卵直徑在0.5～5.5 mm之間[18]，該研究所獲得的永暑礁海域魚卵直徑也處于該范圍，但是比尤卡坦半島(Yucatan Penisula)海域魚卵直徑大[2]，可能是由于其將魚卵保存于酒精中導(dǎo)致魚卵縮水造成的，而該研究是將魚卵樣品直接冷凍保存。

3 討論

近年來，由于過度捕撈和填海造島等人類活動(dòng)干擾，漁業(yè)生態(tài)環(huán)境發(fā)生改變，導(dǎo)致漁業(yè)資源顯著減少[19]。對(duì)魚卵進(jìn)行準(zhǔn)確鑒別是了解魚類繁殖習(xí)性、生活史和種群動(dòng)態(tài)的基礎(chǔ)，可用于鑒別魚卵的形態(tài)特征較少，傳統(tǒng)的分類手段很難準(zhǔn)確鑒別到種的水平。永暑礁為典型的珊瑚環(huán)礁，是熱帶珊瑚礁魚類的重要棲息地和產(chǎn)卵場(chǎng)，目前未見有關(guān)永暑礁魚卵分類鑒別的相關(guān)文獻(xiàn)報(bào)道。DNA條形碼技術(shù)的迅速發(fā)展，為從分子水平準(zhǔn)確鑒定魚卵提供了有效的手段[10,20]。從魚卵中提取高質(zhì)量的DNA用于PCR擴(kuò)增，是利用DNA條形碼鑒別魚卵的前提條件，該研究分析的71個(gè)冷凍保存的魚卵中30個(gè)成功擴(kuò)增測(cè)序獲得了有效的線粒體COI序列信息，成功率為42.3%。部分魚卵未能成功獲得COI序列，可能由于魚卵樣品為冰凍保存，樣品在凍融過程中降解，導(dǎo)致無法提取其基因組DNA，或者由于魚卵樣品體積較小，在DNA提取過程中操作不當(dāng)導(dǎo)致樣品損失，因此，魚卵樣品保存和基因組DNA的提取需要進(jìn)一步的優(yōu)化，提高魚卵鑒定的檢出率。有趣的是圓形魚卵的擴(kuò)增測(cè)序成功率(61.4%)顯著高于橢圓形魚卵的擴(kuò)增測(cè)序成功率(11.1%)，其中原因有待于進(jìn)一步研究。酒精保存魚卵容易導(dǎo)致魚卵脫水，致使魚卵縮小變形，但是酒精保存不需要特殊的設(shè)備，容易操作，也是目前魚卵采集過程中的主要保存方式之一，Leyva等[2]研究表明酒精保存的魚卵DNA條形碼測(cè)序成功率為46.3%，Burghart等[9]的研究表明50%的異丙醇保存魚卵也不影響樣品用于DNA條形碼分析。由此可見目前魚卵的擴(kuò)增成功率并不高，如何開發(fā)高效的魚卵保存和擴(kuò)增方法是有待進(jìn)一步研究的問題。

近年來有關(guān)魚類DNA條形碼參考序列數(shù)據(jù)庫構(gòu)建得到了快速發(fā)展，Díaz等[21]采集分析了79種魚類樣品，測(cè)序分析了其DNA條形碼序列，構(gòu)建了美洲Lower Paraná河魚類DNA條形碼參考序列數(shù)據(jù)庫；Chang等[22]采集分析了1 245種魚類憑證標(biāo)本及其DNA條形碼序列，構(gòu)建了臺(tái)灣海峽魚類DNA條形碼參考序列數(shù)據(jù)庫；Weight等[23]采集分析了521種魚類樣品及其DNA條形碼序列，構(gòu)建了加勒比海魚類DNA條形碼參考序列數(shù)據(jù)庫；中國(guó)水產(chǎn)科學(xué)研究院依托國(guó)家科技基礎(chǔ)工作專項(xiàng)采集分析了763種南海魚類憑證標(biāo)本及其DNA條形碼序列，構(gòu)建了較為完備的南海珊瑚魚類DNA條形碼數(shù)據(jù)庫，這些序列均已錄入到BOLD條形碼數(shù)據(jù)庫中，這些條形碼序列數(shù)據(jù)庫的建立為基于DNA條形碼技術(shù)開展魚卵和仔稚魚鑒定奠定了基礎(chǔ)。北大西洋大陸架海域93.26%的魚卵DNA條形碼序列能夠與BOLD數(shù)據(jù)庫中的參考序列相匹配[8]，在加勒比海珊瑚礁海域僅有85.7%的魚卵DNA條形碼序列能夠與BOLD數(shù)據(jù)庫中的參考序列相匹配[2]，該研究中100%的魚卵DNA條形碼序列能夠與參考數(shù)據(jù)庫中的序列相匹配。這些研究結(jié)果表明DNA條形碼參考數(shù)據(jù)庫的完整性直接影響到魚卵鑒定的效率。

由于線粒體COI基因存在同義突變偏好性，不同種屬間的突變率存在顯著差異[24]，致使種內(nèi)或種、屬間的線粒體COI序列多態(tài)性存在顯著差異，這種差異稱為條形碼間隙(Barcode gap)[25]，DNA條形碼間隙的存在是利用條形碼進(jìn)行物種鑒定的前提，已有的關(guān)于魚類DNA條形碼研究證實(shí)魚類種內(nèi)或種、屬間存在顯著的條形碼間隙[22,26]。該研究表明永暑礁同種魚類間的K2P遺傳距離為0～0.002，同一個(gè)屬不同種間的遺傳距離為0.004～0.132，屬間的遺傳距離為0.117～0.185，存在著明顯的DNA條形碼間隙，分析的所有種類在NJ進(jìn)化樹中均能形成單系支，證實(shí)了DNA條形碼技術(shù)適用于南海珊瑚礁魚類魚卵和仔稚魚鑒定。該研究中南沙珊瑚礁魚類種內(nèi)遺傳距離最高僅為0.002，顯著低于墨西哥尤卡坦半島魚類種內(nèi)遺傳距離(0.006)[2]和臺(tái)灣海峽魚類種內(nèi)遺傳距離(0.151)[21]，可能與該研究中每個(gè)種類使用魚卵的數(shù)量較少有關(guān)，下一步應(yīng)加大魚卵采集和分析的數(shù)量，以便更準(zhǔn)確反映種內(nèi)和種屬間的遺傳關(guān)系。

[1]周美玉,陳 驍,楊圣云.采用DNA條形碼技術(shù)對(duì)廈門海域魚卵、仔稚魚種類的鑒定[J].海洋環(huán)境科學(xué),2015,34(1):120-135.

[2]Leyva C E，Vasquez Y L，Carrillo L，et al.Identifying pelagic fish eggs in the southeast Yucatan Peninsula using DNA barcodes[J].Genome，2016，59：1117-1129.

[3]Steinke D ，Connell A D，Hebert P D N.Linking adults and immatures of South African marine fishes[J].Genome，2016，59(11)：959-967.

[4]萬瑞景,張仁齋.中國(guó)近海及其鄰近海域魚卵與仔稚魚[M].上海:上?？茖W(xué)技術(shù)出版社,2016:1-435.

[5]Zeng Y,Wu Z,Zhang C，et al.DNA barcoding of Mobulid Ray Gill Rakers for Implementing CITES on Elasmobranch in China[J].Sci Rep，2016，6：37567.

[6]Sakinan S，Karahan A，Ok M.Integration of DNA barcoding for the initial recordings of Lessepsian fishes：a case study of the Indo-Pacific slender ponyfish Equulites elongates[J].Fish Biol，2017，90(3)：1054-1061.

[7]Bhatta C M，Sharma A R，Patra B C，et al.DNA barcoding to fishes：current status and future directions[J].Mitochondrial DNA，2016，27(4)：2744-2752.

[8]Leah A，Lewis D E，Richardson E V，et al.Intergrating DNA barcoding of fish egg into ichthyoplankton monitoring programs[J].Fish Bull，2016，114：153-165.

[9]Burghart S E，Woudenberg L V，Daniels C A，et al.Disparity between planktonic fish egg and larval communities as indicated by DNA barcoding[J].Mar Ecol Prog Ser，2014，503(1)：195-204.

[10]Harada A E，Lindgren E A，Hermsmeier M C，et al.Monitoring spawning activity in a southern california marine protected area using molecular identification of fish eggs[J].PLoS ONE，2015，10(8)：e0134647.

[11]Carole C B，Johnson G D.Connectivity across the Caribbean Sea：DNA barcoding and morphology unite an enigmatic fish larva from the Florida straits with a new species of sea bass from deep reefs off Curacao[J].PLoS ONE，2014，9(5)：e97661.

[12]中國(guó)科學(xué)院動(dòng)物研究所等.南海魚類志[M].北京:科學(xué)出版社,1962:1-1184.

[13]張仁齋,李昌華.洋浦灣的魚卵及其仔魚調(diào)查(一)[J].水產(chǎn)學(xué)報(bào),1965，2(4):51-62.

[14]張仁齋,陸穗芬.短尾大眼鯛卵及仔、稚魚的發(fā)育[J].水產(chǎn)學(xué)報(bào),1982,6(3):243-251.

[15]Ward R D，Zemlak T S，Innes B H，et al.DNA barcoding Australia’s fish species[J].Phil Trans R Soc B，2005，360：1847-1857

[16]Ratnasingham S，Hebert P D N.A DNA-based registry for all animal species：The barcode index number (BIN)System[J].PLoS ONE，2013，8(8)：e66213.

[17]Kumar S，Stecher G，Tamura K.MEGA7：Molecular evolutionary genetics analysis version 7.0 for bigger datasets[J].Mol Biol Evol，2016，33(7)：1870-1874.

[18]Matarese A C，Sandknop E M.Identification of fish eggs//Ahlstrom E H，ed.Ontogeny and systematics of fishes[M].California：Spec Publ Am Soc Ichthyol Herpetol，1984：27-31.

[19]湯嬌雯,張 富,陳兆波.我國(guó)魚類生物多樣性保護(hù)策略[J].淡水漁業(yè),2009,39(4):76-77

[20]Mabuchi K.A method for PCR-based identication of fish eggs attached to vegetation in lake Biwa，Japan[J].Conserv Genet Resour，2016，(8)：215-217.

[21]Diaz J，Villanova G V，Brancolini F，et al.First DNA Barcode Reference Library for the Identification of South American Freshwater Fish from the Lower Paraná River[J].PLoS ONE，2016，11(7)：e0157419.

[22]Chang C H，Shao K T，Lin H Y，et al.DNA barcodes of the native ray-finned fishes in Taiwan[J].Mol Ecol Resour，2016，DOI：10.1111/1755-0998.12601

[23]Weight L A，Baldwin C C，Driskell A，et al.Using DNA Barcoding to Assess Caribbean Reef Fish Biodiversity：Expanding Taxonomic and Geographic Coverage[J].PLoS ONE，2012，7(7)：e41059.

[24]Lane N.Biodiversity：on the origin of bar codes[J].Nature，2009，462：272-274.

[25]Ward R D，Holmes B H.An analysis of nucleotide and amino acid variability in the barcode region of cytochrome c oxidase I (cox1)in fishes[J].Mol Ecol Notes，2007，7：899-907.

[26]Pereira L H，Hanner R，F(xiàn)oresti F，et al.Can DNA barcoding accurately discriminate megadiverse Neotropical freshwater fish fauna[J].BMC Genetics，2013，14：20.