基于16S rRNA基因部分序列的長(zhǎng)江口蝦虎魚(yú)科魚(yú)類系統(tǒng)分類

呂 楊，宋 超，劉媛媛，3，趙 峰，張 濤，高 宇，莊 平

（1.南京農(nóng)業(yè)大學(xué)無(wú)錫漁業(yè)學(xué)院，江蘇無(wú)錫 214182；2.中國(guó)水產(chǎn)科學(xué)研究院東海水產(chǎn)研究所，上海 200090；3.上海海洋大學(xué)水產(chǎn)與生命學(xué)院，上海 201306）

基于16S rRNA基因部分序列的長(zhǎng)江口蝦虎魚(yú)科魚(yú)類系統(tǒng)分類

呂 楊1，2，宋 超2，劉媛媛2，3，趙 峰2，張 濤2，高 宇2，莊 平1，2

（1.南京農(nóng)業(yè)大學(xué)無(wú)錫漁業(yè)學(xué)院，江蘇無(wú)錫 214182；2.中國(guó)水產(chǎn)科學(xué)研究院東海水產(chǎn)研究所，上海 200090；3.上海海洋大學(xué)水產(chǎn)與生命學(xué)院，上海 201306）

為了確定線粒體16S rRNA基因在長(zhǎng)江口蝦虎魚(yú)科魚(yú)類系統(tǒng)分類及物種鑒定中的作用，采用16S rRNA基因特異擴(kuò)增測(cè)序及GenBank已有序列聯(lián)合配對(duì)分析的方法，對(duì)長(zhǎng)江口蝦虎魚(yú)科9屬11種魚(yú)類34個(gè)16S rRNA基因片段的序列進(jìn)行比較和系統(tǒng)分類研究。統(tǒng)計(jì)分析顯示，蝦虎魚(yú)科魚(yú)類該片段的A含量明顯高于其它3個(gè)堿基含量，A＋T的平均含量高于G＋C的平均含量，第3密碼子位點(diǎn)G＋C含量最高，其平均值為51.1%，變化范圍為49.8%～53.2%。全部轉(zhuǎn)換位點(diǎn)多于顛換位點(diǎn)，轉(zhuǎn)換/顛換比值為1.45。依據(jù)Maximum Composite Likelihood模型，得出11種蝦虎魚(yú)科魚(yú)類種間遺傳距離平均值為0.151，種內(nèi)為0.002，種間遺傳距離是種內(nèi)遺傳距離的76倍。通過(guò)鄰接法（Neighbor－joining，NJ）構(gòu)建系統(tǒng)發(fā)育樹(shù)，顯示長(zhǎng)江口蝦虎魚(yú)科魚(yú)類為明顯的單系群，并進(jìn)一步佐證把傳統(tǒng)形態(tài)分類中彈涂魚(yú)科的青彈涂魚(yú)（Scartelaos histiophorus）和大彈涂魚(yú)（Boleophthalmus pectinirostris）及鰻蝦虎魚(yú)科的拉氏狼牙蝦虎魚(yú)（Odontamblyopus lacepedii）和紅狼牙蝦虎魚(yú)（Odontamblyopus rubicundus）歸屬于蝦虎魚(yú)科的合理性。聚類結(jié)果顯示紅狼牙蝦虎魚(yú)和拉氏狼牙蝦虎魚(yú)聚在一起，其節(jié)點(diǎn)支持率達(dá)100%，兩者可能為同種異名。本研究表明，線粒體16S rRNA基因序列作為分子標(biāo)記對(duì)蝦虎魚(yú)科魚(yú)類進(jìn)行物種鑒定和系統(tǒng)分類是可行的，可為蝦虎魚(yú)類的親緣關(guān)系分析提供基礎(chǔ)資料。

蝦虎魚(yú)科；16S rRNA；分子系統(tǒng)分類

蝦虎魚(yú)科（Gobiidae）魚(yú)類隸屬蝦虎魚(yú)亞目，該亞目是現(xiàn)存的鱸形目魚(yú)類中最大的一個(gè)類群，在全世界約有2 000余種，包括9科270屬［1］。蝦虎魚(yú)類大多為暖水性或暖溫性底層小型魚(yú)類，生活于近岸海底、淺灣、河口或淡水河流、湖沼中。由于蝦虎魚(yú)體型較小，而且形態(tài)上存在著不同程度的器官退化及特化，給傳統(tǒng)形態(tài)學(xué)分類工作帶來(lái)了較大困難，導(dǎo)致其種類的命名及分類較為混亂，同種異名現(xiàn)象嚴(yán)重，系統(tǒng)分類關(guān)系存在較大爭(zhēng)議。線粒體DNA是核外遺傳物質(zhì)，具有結(jié)構(gòu)簡(jiǎn)單，進(jìn)化較快等優(yōu)點(diǎn)，在動(dòng)物種群的種質(zhì)鑒定和系統(tǒng)分類研究中具有重要作用［2－4］。其中，16S rRNA基因序列具有高度的保守性和特異性，以及較強(qiáng)的種間差異性，已被廣泛應(yīng)用于物種的鑒定及分類研究［5－7］。

本研究利用16S rRNA基因部分序列，對(duì)長(zhǎng)江口11種蝦虎魚(yú)類進(jìn)行種類鑒定，進(jìn)一步驗(yàn)證長(zhǎng)江口蝦虎魚(yú)類的同種異名現(xiàn)象，為長(zhǎng)江口蝦虎魚(yú)科魚(yú)類系統(tǒng)分類提供分子生物學(xué)依據(jù)，對(duì)長(zhǎng)江口蝦虎魚(yú)科魚(yú)類親緣關(guān)系和系統(tǒng)進(jìn)化等研究提供基礎(chǔ)。

1 材料與方法

1.1 樣品采集

本研究所用的雙帶縞蝦虎魚(yú)（Tridentiger bifasciatus）、髭縞蝦虎魚(yú)（Tridentiger barbatus）、竿蝦虎魚(yú)（Luciogobius guttatus）、阿部鯔蝦虎魚(yú)（Mugilogobius abei）、爪哇擬蝦虎魚(yú)（Pseudogobius javanicus）、斑尾刺蝦虎魚(yú)（Acanthogobius ommaturus）、睛尾蝌蚪蝦虎魚(yú)（Lophiogobius ocellicauda）、拉氏狼牙蝦虎魚(yú)（Odontamblyopus lacepedii）、青彈涂魚(yú)（Scartelaos histiophorus）等蝦虎魚(yú)科魚(yú)類及中國(guó)花鱸（Lateolabrax maculatus）和棘頭梅童魚(yú)（Collichthys lucidus）均采自長(zhǎng)江口水域，現(xiàn)場(chǎng)采樣并保存于95%乙醇。紅狼牙蝦虎魚(yú)（O.rubicundus）和大彈涂魚(yú)（Boleophthalmus pectinirostris）的序列特征為NCBI中檢索而得。

1.2 DNA提取

利用海洋動(dòng)物基因組DNA提取試劑盒提取樣品總DNA。具體操作如下：取組織樣品于2.0 mL離心管中，加入200μL LB Buffer、10μL蛋白酶K和600μL NB Buffer，56℃恒溫水浴鍋裂解消化2 h；離心并將上清液全部加入Mini Spin純化柱，離心去除濾過(guò)液；向Mini Spin純化柱中加入RPE Buffer，去除蛋白質(zhì)和其它雜質(zhì)；向Mini Spin純化柱中加入Wash Buffer洗脫DNA；將Mini Spin純化柱轉(zhuǎn)入1.5 mL離心管中，在吸附膜中加入30μLRNase Free H2O溶液，將DNA溶于雙蒸水（RNase Free H2O）中，4℃冰箱保存?zhèn)溆谩?/p>

1.3 引物合成與PCR擴(kuò)增

引物參照PALUMBI等［8］設(shè)計(jì)，由生工生物工程（上海）股份有限公司合成，引物序列為：

PCR反應(yīng)總體積為50μL，其中PremixTaq25μL（1.25 U），引物16S ar/16Sbr各1μL（20 μmol·L－1），根據(jù)純度和質(zhì)量加入適量DNA模板，補(bǔ)加雙蒸水至50μL。PCR反應(yīng)條件為：94℃預(yù)變性5 min；94℃變性30 s；55℃退火30 s；72℃延伸1 min；30個(gè)循環(huán)；最后72℃延伸10 min；PCR產(chǎn)物4℃保存?zhèn)溆谩Ｃ看畏磻?yīng)設(shè)置不含DNA模板的空白對(duì)照。

1.4 DNA測(cè)序和序列下載

將PCR產(chǎn)物用1%瓊脂糖凝膠電泳檢測(cè)，然后送至上海桑尼生物科技有限公司，經(jīng)ABI PRISMTM3730XL DNA Analyzer測(cè)序儀，以16S ar作為測(cè)序引物進(jìn)行單向測(cè)定。

從GenBank中下載蝦虎魚(yú)科其它魚(yú)類16S rRNA基因序列，與本研究所檢測(cè)的序列共計(jì)9屬11種34個(gè)個(gè)體的16S rRNA基因進(jìn)行同源性序列分析，并選擇中國(guó)花鱸和棘頭梅童魚(yú)共2種5個(gè)同源序列作為外群。所分析物種的16S rRNA基因及其相關(guān)信息見(jiàn)表1。

表1 蝦虎魚(yú)科11種魚(yú)類的16S rRNA基因信息Tab.1 Information of 16S rRNA genes of 11 Gobiidae species

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

將所測(cè)序列用NCBI的Blast工具進(jìn)行相似性檢索，確認(rèn)所得序列為目的片段。利用Clustal W對(duì)本研究所得序列和GenBank下載的相關(guān)序列進(jìn)行編輯、校對(duì)和排序，并輔以人工校正［9］。用MEGA 5.0軟件統(tǒng)計(jì)序列的平均堿基組成、變異位點(diǎn)、簡(jiǎn)約信息位點(diǎn)、單獨(dú)位點(diǎn)、轉(zhuǎn)換與顛換位點(diǎn)的比例，進(jìn)行系統(tǒng)發(fā)育分析。采用Maximum Composite Likelihood模型計(jì)算種內(nèi)和種間遺傳距離，用鄰接法（Neighbor－joining，NJ）構(gòu)建分子系統(tǒng)樹(shù)，并采用自展法檢驗(yàn)（Bootstrap）1 000次重復(fù)檢驗(yàn)分子系統(tǒng)樹(shù)各分支的置信度［10］。

2 結(jié)果與分析

2.1 長(zhǎng)江口蝦虎魚(yú)科魚(yú)類16S rRNA基因序列特征

本研究所測(cè)得蝦虎魚(yú)科魚(yú)類PCR產(chǎn)物大小約為560 bp，采用Clustal W對(duì)本研究所測(cè)31個(gè) 16S rRNA基因序列和3個(gè)GenBank中下載的蝦虎魚(yú)科魚(yú)類16S rRNA同源序列進(jìn)行比對(duì)，保留共有序列，長(zhǎng)度492 bp，編碼164個(gè)氨基酸。用MEGA 5.0計(jì)算這34條序列的堿基組成，其堿基平均含量為：T 21.6%，C 24.7%，A 30.6%，G 23.0%，A的含量明顯高于其它3個(gè)堿基的含量。A＋T的平均含量為52.3%，高于G＋C的含量47.7%（表2）。對(duì)所研究的11種蝦虎魚(yú)科魚(yú)類34個(gè)16S rRNA基因片段而言，第3密碼子位點(diǎn)平均GC含量最高，A、T、G、C在1、2、3密碼子位點(diǎn)的平均含量差別明顯，可見(jiàn)，16S rRNA序列密碼子的堿基使用頻率存在偏向性。

11種蝦虎魚(yú)科魚(yú)類16S rRNA基因片段的平均GC含量為48.1%，變化范圍為46.9%～50.3%。其中，第3密碼子位點(diǎn)的平均GC含量高于第1密碼子和第2密碼子，第3個(gè)密碼子位點(diǎn)的平均GC含量變化范圍為49.8%～53.2%（表3）。

表2 蝦虎魚(yú)科11種魚(yú)類的16S rRNA基因部分序列中各堿基平均分布頻率Tab.2 Average nucleotide frequencies of 16S rRNA partial sequence of 11 Gobiidae species （%）

表3 蝦虎魚(yú)科11種魚(yú)類16S rRNA基因密碼子位點(diǎn)GC含量Tab.3 The GC content of 16S rRNA partial sequence of 11 Gobiidae species in all codon positions

表4所示為所有個(gè)體序列核苷酸變異情況，其中不變位點(diǎn)在第3密碼子位點(diǎn)最多，轉(zhuǎn)換和顛換位點(diǎn)均在第2密碼子位點(diǎn)最多，轉(zhuǎn)換位點(diǎn)與顛換位點(diǎn)的比值在第1密碼子位點(diǎn)最小。在全部492個(gè)位點(diǎn)中有155個(gè)變異位點(diǎn)，約占總數(shù)的32%，33個(gè)插入/缺失位點(diǎn)，145個(gè)簡(jiǎn)約信息位點(diǎn)，10個(gè)單獨(dú)位點(diǎn)。單倍型的個(gè)數(shù)為293個(gè)，雙倍型個(gè)數(shù)為63個(gè)，4倍型個(gè)數(shù)為49個(gè)。全部位點(diǎn)的轉(zhuǎn)換位點(diǎn)多于顛換位點(diǎn)，轉(zhuǎn)換位點(diǎn)與顛換位點(diǎn)的比值R為1.45。

2.2 蝦虎魚(yú)科魚(yú)類16S rRNA基因編碼蛋白質(zhì)中氨基酸使用頻率

用MEGA 5.0計(jì)算得到11種蝦虎魚(yú)科魚(yú)類的共有序列編碼164個(gè)氨基酸（表5），其中70個(gè)發(fā)生變異，占氨基酸總數(shù)的43%。在氨基酸的使用頻率中，最高的為亮氨酸（Leu），其在青彈涂魚(yú)中使用頻率最高（13.69%），在拉氏狼牙蝦虎魚(yú)中使用頻率最低（9.87%）。其次，使用頻率較高的有蘇氨酸（Thr）和絲氨酸（Ser），而使用頻率最低的為組氨酸（His），其在斑尾刺蝦虎魚(yú)、紅狼牙蝦虎魚(yú)以及拉氏狼牙蝦虎魚(yú)中的使用頻率均為0。另外，酪氨酸（Tyr）和色氨酸（Trp）在11種蝦虎魚(yú)類中使用頻率也相對(duì)較小。

2.3 種間及種內(nèi)的遺傳距離

用Maximum Composite Likelihood模型計(jì)算11種蝦虎魚(yú)科魚(yú)類的種間和種內(nèi)遺傳距離（表6）。11種蝦虎魚(yú)類的種間遺傳距離在0.055～0.221之間，其平均值為0.151，種內(nèi)遺傳距離平均值為0.002，種間遺傳距離為種內(nèi)遺傳距離的76倍。91%以上的種間遺傳距離大于0.100，而種內(nèi)遺傳距離均小于0.010。種內(nèi)遺傳距離較大的為斑尾刺蝦虎魚(yú)（0.06）；而種間遺傳距離最大值在拉氏狼牙蝦虎魚(yú)和斑尾刺蝦虎魚(yú)間（0.221），而最小值在拉氏狼牙蝦虎魚(yú)和紅狼牙蝦虎魚(yú)間（0.055）。

表4 蝦虎魚(yú)科11種魚(yú)類16S rRNA基因部分序列各密碼子堿基變異情況Tab.4 Sequence variations of 16S rRNA partial sequence of 11 Gobiidae species

表5 蝦虎魚(yú)科11種魚(yú)類16S rRNA基因編碼蛋白質(zhì)氨基酸的使用頻率Tab.5 The am ino acid frequency of 16S rRNA gene encoding protein for 11 Gobiidae species（%）

表6 蝦虎魚(yú)科11種魚(yú)類的種間和種內(nèi)遺傳距離Tab.6 Genetic distance pairw ise-species and w ithin-species in 11 Gobiidae species

2.4 蝦虎魚(yú)科魚(yú)類的分子系統(tǒng)樹(shù)分析

采用NJ法對(duì)蝦虎魚(yú)科9屬11種34個(gè)體及中國(guó)花鱸和棘頭梅童魚(yú)2種外群的16S rRNA基因序列構(gòu)建分子系統(tǒng)樹(shù)（圖1）。整個(gè)分子系統(tǒng)樹(shù)分為兩個(gè)大的分支，第一支由中國(guó)花鱸和棘頭梅童魚(yú)2外群組成，第二支由11種蝦虎科魚(yú)類組成。蝦虎魚(yú)類的分子系統(tǒng)樹(shù)又分為2支，其中一支由斑尾刺蝦虎魚(yú)和睛尾蝌蚪蝦虎魚(yú)組成，其余9種蝦虎魚(yú)類組成另外一支。在分子系統(tǒng)樹(shù)中除了拉氏狼牙蝦虎魚(yú)和紅狼牙蝦虎魚(yú)聚在一起外，其它9種蝦虎魚(yú)科魚(yú)類同一種的不同個(gè)體均可聚在同一分支內(nèi)，9種蝦虎魚(yú)可形成單系（圖1）。拉氏狼牙蝦虎魚(yú)和紅狼牙蝦虎魚(yú)聚在一起，其自展置信值達(dá)100，呈現(xiàn)出非常近的親緣關(guān)系。

3 討論

本研究中蝦虎魚(yú)類的16S rRNA基因片段中A＋T含量明顯高于G＋C含量，該結(jié)果與ASAKAWA等［11］提出的后生動(dòng)物線粒體基因的堿基構(gòu)成具有明顯的偏倚性，A＋T含量通常較為豐富，顯著高于G＋C含量的結(jié)果一致，該特點(diǎn)已在許多魚(yú)類中得到了驗(yàn)證［12－14］。其中，孫希福［15］通過(guò)對(duì)中國(guó)沿海10種蝦虎魚(yú)類堿基組成的研究也證實(shí)了這一結(jié)論。本研究中蝦虎魚(yú)16S rRNA基因全部密碼子位點(diǎn)的轉(zhuǎn)換值與顛換值之比為1.45，表明蝦虎魚(yú)類的16S rRNA基因序列的突變已經(jīng)達(dá)到飽和狀態(tài)［16］。蝦虎魚(yú)類16S rRNA基因編碼蛋白質(zhì)的氨基酸使用頻率中，最高的為亮氨酸，而使用頻率最低的為組氨酸?？梢?jiàn)，蝦虎魚(yú)科魚(yú)類16S rRNA基因編碼蛋白質(zhì)的氨基酸使用頻率存在明顯的偏向性。

HEBERT等［17－19］認(rèn)為COⅠ序列在物種間的遺傳差異通常大于2%，并提出只有當(dāng)種間遺傳距離顯著大于種內(nèi)遺傳距離，并且這種差異達(dá)到約10倍以上時(shí)，COⅠ基因才可對(duì)物種進(jìn)行有效的鑒定。而對(duì)16S rRNA基因在物種鑒定中的應(yīng)用，迄今還未見(jiàn)相應(yīng)的鑒定標(biāo)準(zhǔn)，但可以確定的是較大的種間序列差異是對(duì)物種進(jìn)行準(zhǔn)確鑒定的先決條件［20］。本研究中11種蝦虎魚(yú)類的種間遺傳距離平均為0.151，種內(nèi)遺傳距離平均為0.002，種間遺傳距離為種內(nèi)遺傳距離的76倍。種內(nèi)具有較高的同源性，而不同種間的多態(tài)性差異比較明顯，表明該基因適用于對(duì)長(zhǎng)江口蝦虎魚(yú)類物種進(jìn)行種類鑒定［21］。縱觀本研究分子系統(tǒng)樹(shù)的拓?fù)浣Y(jié)構(gòu)，主要分為兩大支，第一大支由中國(guó)花鱸和棘頭梅童魚(yú)組成，第二大支由11種蝦虎魚(yú)類組成。在蝦虎魚(yú)類的分支中包括了傳統(tǒng)分類中彈涂魚(yú)科的青彈涂魚(yú)和大彈涂魚(yú)，以及鰻蝦虎魚(yú)科的拉式狼牙蝦虎魚(yú)和紅狼牙蝦虎魚(yú)［22］。從親緣關(guān)系圖譜來(lái)看，上述四種魚(yú)類與其它蝦虎魚(yú)類的親緣關(guān)系密切，均可歸為蝦虎魚(yú)科。該研究結(jié)果與伍漢霖等［1］主編的《中國(guó)動(dòng)物志》中把青彈涂魚(yú)、大彈涂魚(yú)、拉氏狼牙蝦虎魚(yú)和紅狼牙蝦虎魚(yú)歸屬為蝦虎魚(yú)科的分類結(jié)果一致。蝦虎魚(yú)類組成的一大支又分為2個(gè)分支，其中一個(gè)分支由斑尾刺蝦虎魚(yú)和睛尾蝌蚪蝦虎魚(yú)組成，表明斑尾刺蝦虎魚(yú)和睛尾蝌蚪蝦虎魚(yú)具有較近的系統(tǒng)進(jìn)化關(guān)系，該結(jié)果與于亞男等［23］基于線粒體COⅠ基因?qū)﹂L(zhǎng)江口蝦虎魚(yú)類系統(tǒng)進(jìn)化的研究結(jié)果一致。在其余9種蝦虎魚(yú)類構(gòu)成的另外一支中，又分為2個(gè)支系。在其中一個(gè)支系中，髭縞蝦虎魚(yú)和雙帶縞蝦虎魚(yú)為姊妹分支，其遺傳距離只有0.092，說(shuō)明兩者親緣關(guān)系較近，同屬于縞蝦虎魚(yú)屬［1］。竿蝦虎魚(yú)和阿部鯔蝦虎魚(yú)為姊妹分支，兩者的遺傳距離為0.084，具有較近的親緣關(guān)系；這四種蝦虎魚(yú)與爪哇擬蝦虎魚(yú)共同組成了一個(gè)支系，這五種蝦虎魚(yú)類的棲息水域和攝食習(xí)性相近，說(shuō)明它們具有相似的起源和相近的親緣關(guān)系。在另一個(gè)支系中，紅狼牙蝦虎魚(yú)和拉氏狼牙蝦虎魚(yú)聚在一起，兩者的節(jié)點(diǎn)支持率為100%，其種間遺傳距離在本研究的所有種間遺傳距離中最小為0.055，而在《中國(guó)動(dòng)物志》中把紅狼牙蝦虎魚(yú)描述為拉氏狼牙蝦虎魚(yú)的別名［1］，說(shuō)明兩者具有相似的形態(tài)特征和很近的親緣關(guān)系，可能為同種異名；青彈涂魚(yú)和大彈涂魚(yú)為姊妹分支，兩者的遺傳距離為0.098，表明二者具有較近的系統(tǒng)進(jìn)化關(guān)系。在傳統(tǒng)的分類中青彈涂魚(yú)和大彈涂魚(yú)屬于彈涂魚(yú)科［22］，本研究中兩種彈涂魚(yú)與其它9種蝦虎魚(yú)間的平均遺傳距離為0.143小于9蝦虎魚(yú)間的平均遺傳距離0.152，進(jìn)一步佐證了把這2種彈涂魚(yú)科魚(yú)類歸屬為蝦虎魚(yú)科的合理性。

圖1 NJ法構(gòu)建的蝦虎魚(yú)科魚(yú)類分子系統(tǒng)樹(shù)Fig.1 NJ tree from analysis of 16S rRNA gene data of 11 Gobiidae species

線粒體不同基因片段的核苷酸替代速率存在明顯差異，16S rRNA基因是線粒體基因組中研究較多的基因，其進(jìn)化速度適中，適合于分子系統(tǒng)分類和親緣關(guān)系研究，許多學(xué)者已利用16S rRNA基因部分序列對(duì)多種物種進(jìn)行了種間系統(tǒng)進(jìn)化關(guān)系分析［24－26］?；?6S rRNA基因構(gòu)建的蝦虎魚(yú)科分子系統(tǒng)樹(shù)顯示，除了拉氏狼牙蝦虎魚(yú)與紅狼牙蝦虎魚(yú)聚在一起，兩者可能為同種異名外，同種魚(yú)的不同個(gè)體均可聚在同一分支內(nèi)，表現(xiàn)出明顯的單系性?？傮w來(lái)講，該分子系統(tǒng)分類結(jié)果與形態(tài)學(xué)的分類單元?jiǎng)澐只疽恢拢樾碌姆诸愔邪亚鄰椡眶~(yú)、大彈涂魚(yú)、拉氏狼牙蝦虎魚(yú)和紅狼牙蝦虎魚(yú)歸屬為蝦虎魚(yú)科的合理性提供了佐證。線粒體16S rRNA基因作為分子標(biāo)記獲得的分類結(jié)果與用線粒體COⅠ和核Ryr3等分子標(biāo)記方法獲得的結(jié)果一致［27］，進(jìn)一步說(shuō)明線粒體16S rRNA基因作為分類條形碼適合對(duì)蝦虎魚(yú)類進(jìn)行分子鑒定和系統(tǒng)分類研究，為整個(gè)蝦虎魚(yú)科魚(yú)類分子水平的系統(tǒng)發(fā)生和進(jìn)化等研究提供基礎(chǔ)。

［1］ 伍漢霖，鐘俊生.中國(guó)動(dòng)物志（硬骨魚(yú)綱鱸形目蝦虎魚(yú)亞目）［M］.北京：科學(xué)出版社，2008：136－718.

WU H L，ZHONG J S.FAUNA SINICA.Osteichthyes.PERCIFORMES.Gobioidei［M］.Beijing：Science Press，2008：136－718.

［2］ 郭新紅，劉少軍，劉 巧，等.魚(yú)類線粒體DNA研究新進(jìn)展［J］.遺傳學(xué)報(bào)，2004，31（9）：983－1000.

GUO X H，LIU S J，LIU Q，et al.New Progresses on Mitochondrial DNA in Fish［J］.Acta Genetica Sinica，2004，31（9）：983－1000.

［3］ SERENO-URIBE A L，PINACHO-PINACHO C D，GARCIA-VARELA M，et al.Using mitochondrial and ribosomal DNA sequences to test the taxonomic validity ofClinostomum complanatumRudolphi，1814 in fish-eating birds and freshwater fishes in Mexico，with the description of a new species［J］.Parasitology Research，2013，112（8）：2855－2870.

［4］ AWOTUNDE EO，BEMJIM N，OLOWOFESO O，et al.Mitochondrial DNA sequence analyses and phylogenetic relationships among two Nigerian Goat Breeds and the South African Kalahari Red［J］.Animal Biotechnology，2015，26（3）：180－187.

［5］ AKIMOTO S，KINOSHITA S，SEZAKI K，et al.Identification of alfonsino and related fish species belonging to the genus Beryx with mitochondrial 16S rRNA gene and its application on their pelagic eggs［J］.Fisheries Science，2002，68（6）：1242－1249.

［6］ GARCIA-VAZQUEZ E，ALVAREZP，LOPESP，et al.PCR-SSCP of the 16S rRNA gene，a simple methodology for species identification of fish eggs and larvae［J］.Scientia Marina，2006（70）：13－21.

［7］ DURAND J D，DIATTA M A，DIOP K，et al.Multiplex 16S rRNA haplotype-specific PCR，a rapid and convenientmethod for fish species identification：an application to West African Clupeiform larvae［J］.Molecular Ecology Resources，2010，10（3）：568－572.

［8］ PALUMBIS，MARTIN A，ROMANO S，et al.The Simple Fool’s Guide to PCR（Ver.2.0）［M］.Honolulu：University of Hawaii，1991：25－28.

［9］ KUMAR S，TAMURA K，NEIM.MEGA3：integrated software for molecular evolutionary genetics analysis and sequence alignment［J］.Briefings in Bioinformatics，2004，5（2）：150－163.

［10］ TAMURA K，PETERSON D，PETERSON N，et al.MEGA5：Molecular evolutionary genetics analysis using Maximum Likelihood，Evolutionary Distance，and Maximum Parsimony Methods［J］.Molecular Biology and Evolution，2011，28（10）：2731－2739.

［11］ ASAKAWA S，KUMAZAWA Y，ARAKI T，et al.Strand-specific nucleotide composition bias in echinodermand vertebrate mitochondrial genomes［J］.Journal ofMolecular Evolution，1991，32（6）：511－520.

［12］ 鄭文娟，朱世華，鄒記興，等.基于16S rRNA部分序列探討12種鲹科魚(yú)類的分子系統(tǒng)進(jìn)化關(guān)系［J］.水產(chǎn)學(xué)報(bào)，2008，32（6）：847－854.

ZHENGW J，ZHU SH，ZOU JX，etal.Molecular phylogenetic relationship of Carangidae based on partial sequence ofmitochondrial16S ribosomal RNA gene［J］.Journal of Fisheries of China，2008，32（6）：847－854.

［13］ LAKRA W S，GOSWAMIM，GOPALAKRISHANA A.Molecular identification and phylogenetic relationships of seven Indian Sciaenids（Pisces：Perciformes，Sciaenidae）based on 16S rRNA and cytochrome c oxidase subunit Imitochondrial genes［J］.Molecular Biology Reports，2009，36（5）：831－839.

［14］ 張?jiān)凑?，?偉，姜志強(qiáng)，等.基于16S rRNA部分序列探討部分鳚亞目魚(yú)類的分子系統(tǒng)進(jìn)化關(guān)系［J］.海洋科學(xué)，2012，36（9）：89－95.

ZHANG Y Z，WANG W，JIANG Z Q，et al.Molecular phylogenetic relationship of blennioidei based on partial sequence of mitochondrial 16S ribosomal RNA gene［J］.Marine Sciences，2012，36（9）：89－95.

［15］ 孫希福.基于形態(tài)學(xué)和分子生物學(xué)資料探討中國(guó)沿海10種蝦虎魚(yú)類的系統(tǒng)發(fā)育關(guān)系［D］.青島：中國(guó)海洋大學(xué)，2009.

SUN X F.Phylogenetic Analysis of ten Gobiidae species of neritic waters of China Based on Morphological and Molecular Markers［D］.Qingdao：Ocean University of China，2009.

［16］ KNIGHT A，MINDELL D P.Substitution bias，weighting of DNA sequence evolution，and the phylogenetic position of Fea’s Viper［J］.Systematic Biology，1993，42（1）：18－31.

［17］ HEBERT PD N，CYWINSKA A，BALL SL，et al.Biological identifications through DNA barcodes［J］.Proceedings of the Royal Society of London.Series B：Biological Sciences，2003，270（1512）：313－321.

［18］ HEBERTPD N，RATNASINGHAM S，DEWAARD JR.Barcoding animal life：cytochrome c oxidase subunit I divergences among closely related species［J］.Proceedings of the Royal Society of London.Series B：Biological Sciences，2003，270（Sup 1）：96－99.

［19］ HEBERT PD N，STOECKLE M Y，ZEMLAK T S，et al.Identification of birds through DNA barcodes［J］.PLoSBiology，2004，2（10）：1657－1663.

［20］ 王淑英，時(shí) 偉，江金霞，等.COⅠ和16S rRNA基因序列在鰨科（Soleidae）魚(yú)類種類鑒定中的適用性研究［J］.熱帶海洋學(xué)報(bào)，2014，33（3）：57－65.

WANG S Y，SHI W，JIANG J X，et al.Applicability ofmitochondrial COⅠand 16S rRNA gene sequences in species identification of sole fish（Pleuronectiformes：Soleidae）［J］.Journal of Tropical Oveanography，2014，33（3）：57－65.

［21］ 金逍逍.東海16種蝦虎魚(yú)類線粒體全基因組序列分析及系統(tǒng)發(fā)育研究［D］.舟山：浙江海洋學(xué)院，2013.

JIN X X.Complete mitochondrial genome and phylogenetic analysis of 16 gobies in East China Sea［D］.Zhoushan：Zhejiang Ocean University，2013.

［22］ 莊 平，王幼槐，李圣法，等.長(zhǎng)江口魚(yú)類［M］.上海：上海科學(xué)技術(shù)出版社，2006：426－428.

ZHUANG P，WANG Y H，LIS F，et al.Fishes of the Yangtze Estuary［M］.Shanghai：Shanghai Scientific and Technical Publishers，2006：426－428.

［23］ 于亞男，宋 超，侯俊利，等.基于線粒體COⅠ基因部分序列的長(zhǎng)江口蝦虎魚(yú)科魚(yú)類系統(tǒng)分類［J］.淡水漁業(yè)，2014，44（5）：3－8.

YU Y N，SONG C，HOU J L，et al.System classification for Gobiidae derived from the Yangtze River based on mitochondrial COⅠgene partial sequence［J］.Freshwater Fisheries，2014，44（5）：3－8.

［24］ 王 茜，齊興柱，駱 劍，等.尖塘鱧屬魚(yú)類線粒體16S rRNA基因序列變異及分子系統(tǒng)進(jìn)化［J］.海南大學(xué)學(xué)報(bào)：自然科學(xué)版，2009，27（3）：245－251.

WANG Q，QI X Z，LUO J，et al.Sequences Variation and Molecular Phylogeny of Mitochondrial 16S rRNA Gene of the Genus Oxyleotris［J］.Natural Science Journal of Hainan University，2009，27（3）：245－251.

［25］ 唐優(yōu)良.基于16S rRNA基因序列的中國(guó)鰈形目魚(yú)類分子系統(tǒng)學(xué)研究［D］.廣州：暨南大學(xué)，2011.

TANG Y L.Phylogenetic Study of Flatfish（Order：Pleuronectiformes）from China based on 16S Gene Sequence［D］.Guangzhou：Jinan University，2011.

［26］ 吳仁協(xié)，李 超，劉 靜.鯧亞目魚(yú)類線粒體16S rRNA基因序列變異及其分子系統(tǒng)進(jìn)化關(guān)系［J］.水產(chǎn)學(xué)報(bào)，2013，37（1）：16－25.

WU R X，LI C，LIU J.Mitochondrial 16S rRNA sequence variation and molecular phylogenetic relationships of Stromateoidei fishes（Teleostei：Perciformes）［J］.Journal of Fisheries of China，2013，37（1）：16－25.

［27］ 顧翠平.浙江沿海刺蝦虎魚(yú)屬（Acanthogobius）和縞蝦虎魚(yú)屬（Tridentiger）分子分類研究［D］.舟山：浙江海洋學(xué)院，2014.

GU C P.Molecular Barcoding and Phylogeography ofAcanthogobiusandTridentiger（Teleostei：Gobiidae）of Zhejiang Coast［D］.Zhoushan：Zhejiang Ocean University，2014.

Phylogenetic analysis of Gobiidae in the Yangtze Estuary based on partial sequence of M itochondrial 16S rRNA

LV Yang1，2，SONG Chao2，LIU Yuan-yuan2，3，ZHAO Feng2，ZHANG Tao2，Gao Yu2，ZHUANG Ping1，2
（1.Wuxi Fisheries College，Nanjing Agricultural University，Wuxi214182，China；2.East China Sea Fisheries Research Institute，Chinese Academy of Fishery Sciences，Shanghai200090，China；3.College of Fisheries and Life Science，Shanghai Ocean University，Shanghai 201306，China）

To determine the role of mitochondrial 16S rRNA genes in classification and identification of species，a total of 34 single individuals from 11 species pertaining to 9 genera of Gobiidae were barcoded by 16S rRNA，and compared with of Gobiidae species recorded in GenBank.Total genomic DNA was extracted from each sample usingmarine animal genome DNA extraction kit.PCR amplification was performed in total volume of 50μL of PCR mixture，and PCR productswere purified and sequenced in single direction using an ABIPRISMTM3 730 XL Automated Sequencer.By statistical analysis of MEGA 5.0 software，the average content of A（30.6%）was significantly higher than that of other three bases，and the average content of A＋T（52.3%）was higher than the G＋C content（47.7%）.The G＋C content of the third position codon was higher than that of the other positions，whose average contentwas 51.1%and the range was from 49.8%to 53.2%.All transitional pairs（si）were slightlymore than that of transversional pairs（sv），and the ratio（R＝si/sv）was 1.45.Based on Maximum Composite Likelihood model，the average distance pairwise-species and within-specieswas 0.151 and 0.002，respectively.The average distance pairwise-species was 76 times that ofwithin-species.According to the neighbor joining（NJ）trees for all34 sequences，itwas demonstrated that Gobiidae in the Yangtze Estuary was amonophyletic group，which further confirmed that the classification ofScartelaos histiophorusandBoleophthalmus pectinirostris，Odontamblyopus lacepediiandOdontamblyopus rubicundus，as Gobiidae was reasonable.O.lacepediiandO.rubicunduswere placed in the same clade with support of a high bootstrap value of 100%，indicating that these two speciesmay be synonyms.Our results highlight that the information from 16S rRNA sequences can not only filter out the synonym of the same species，but also can be used to carry out effective identification for Gobiidae species，which further shows that 16S rRNA is feasible as classification barcode and can provide basis for phylogenetic relationship researches on Gobiidae.

Gobiidae；16S rRNA；molecular systematics

Q 754

A

1004－2490（2016）01－0017－09

2015－04－09

公益性行業(yè)（農(nóng)業(yè)）科研專項(xiàng)（201203065）

呂 楊（1989－），女，碩士研究生。E-mail：lvyang＿2015＠163.com

莊 平，研究員。E-mail：Pzhuang＠hotmail.com