基于ND1基因序列的16種胡椒鯛魚類分子系統(tǒng)進(jìn)化關(guān)系研究

梁日深，梁志剛，方浩航，譚佳瑜，彭穎暄，薛濤

（1.仲愷農(nóng)業(yè)工程學(xué)院生命科學(xué)學(xué)院，廣東廣州 510642）

基于ND1基因序列的16種胡椒鯛魚類分子系統(tǒng)進(jìn)化關(guān)系研究

梁日深1，梁志剛1，方浩航1，譚佳瑜1，彭穎暄1，薛濤1

（1.仲愷農(nóng)業(yè)工程學(xué)院生命科學(xué)學(xué)院，廣東廣州 510642）

本研究基于線粒體DNA分子標(biāo)記NADH脫氫酶亞基1基因（ND1）部分序列分析了印度-西太平洋海域的胡椒鯛屬及少棘胡椒鯛屬共16種胡椒鯛魚類的系統(tǒng)進(jìn)化關(guān)系。采用最大簡(jiǎn)約法及最大似然法構(gòu)建了系統(tǒng)進(jìn)化樹。進(jìn)化樹上，16種胡椒鯛魚類共形成兩個(gè)形態(tài)特征截然不同的類群，其中類群一的種類體表具有各種鮮艷的顏色與多變的斑紋，主要分布在西太平洋地區(qū)；類群二的種類顏色灰暗單一，少數(shù)體表具有暗色斑點(diǎn)與條紋，大多數(shù)分布在印度洋地區(qū)。此外，進(jìn)化樹上，少棘胡椒鯛屬的種類都位于胡椒鯛屬的類群一內(nèi)部，與斑胡椒鯛，暗點(diǎn)胡椒鯛形成的姐妹種有較近的親緣關(guān)系?；谶z傳距離的數(shù)據(jù)，少棘胡椒鯛屬與胡椒鯛屬間的平均遺傳距離小于胡椒鯛屬內(nèi)部種間的遺傳距離，顯示出少棘胡椒鯛屬與胡椒鯛屬有非常近的親緣關(guān)系，支持少棘胡椒鯛歸類于胡椒鯛屬的觀點(diǎn)。

系統(tǒng)發(fā)育；仿石鱸科；胡椒鯛屬；少棘胡椒鯛屬；ND1

1 引言

胡椒鯛魚類Sweetlips，包括胡椒鯛屬Plectorhinchus以及少棘胡椒鯛屬Diagramma的種類，隸屬鱸形目Perciformes鱸亞目Percoidei仿石鱸科Haemulidae胡椒鯛亞科Plectorhynchinae，廣泛分布于印度洋及太平洋西部溫帶，亞熱帶及熱帶沿海及島礁海域，是重要的海洋經(jīng)濟(jì)食用及觀賞性魚類。胡椒鯛在世界上約有30余種，中國(guó)大約記載8～12種，主要分布在南海［1—4］。

形態(tài)上，胡椒鯛魚類近緣物種十分相似，許多種類都具有近似的體色和斑紋，且體色與斑紋會(huì)隨著個(gè)體的逐步發(fā)育發(fā)生巨大的變化，使幼魚與成魚的形態(tài)截然不同。如一些種類幼魚具有橫跨體表的水平條紋，而發(fā)育成魚后條紋逐漸破裂成無數(shù)小圓點(diǎn)（黃點(diǎn)胡椒鯛Plectorhinchus flavomaculatus），另有種類幼魚體表具不規(guī)則斑紋，而成魚后卻變?yōu)槊黠@的水平縱帶（東方胡椒鯛P.orientalis）。此外，許多種類的成魚與其自身幼魚形態(tài)相去甚遠(yuǎn)的同時(shí)，卻與其他種類的幼魚，半成魚形態(tài)非常相似，如東方胡椒鯛P.orientalis成魚與斜紋胡椒鯛P.lineatus、四帶胡椒鯛P. diagrammus幼魚，體表均為數(shù)條黑色水平縱帶，各鰭為黃色，極易混淆。胡椒鯛屬的形態(tài)差異多樣性給傳統(tǒng)基于外部性狀的分類鑒定帶來很大的困擾，許多種類的分類與命名相互混淆［1—3］。在分子水平上，國(guó)內(nèi)外學(xué)者對(duì)胡椒鯛魚類的系統(tǒng)進(jìn)化雖有報(bào)道，但涉及胡椒鯛種類的數(shù)量都很少。朱世華等［5］利用細(xì)胞色素b基因部分序列對(duì)5種仿石鱸魚類進(jìn)行了系統(tǒng)進(jìn)化研究，結(jié)果顯示胡椒鯛P.pictus與胡椒鯛魚類關(guān)系較遠(yuǎn)，支持胡椒鯛從胡椒鯛屬Plectorhinchus劃分出去的觀點(diǎn)，但其研究只涉及兩種胡椒鯛；任崗等［6］利用16S r RNA序列分析了11種仿石鱸科魚類系統(tǒng)進(jìn)化關(guān)系，其研究包含的4種胡椒鯛中，胡椒鯛位于胡椒鯛屬的內(nèi)部。Sanciangco等［7］與Tavera等［8］利用多個(gè)分子標(biāo)記對(duì)北美海域的仿石鱸科魚類的進(jìn)行較為全面的系統(tǒng)分類研究，但他們重點(diǎn)放在石鱸亞科Haemulinae的種類，對(duì)胡椒鯛亞科的分析內(nèi)容很少。目前關(guān)于胡椒鯛魚類分子研究基本上都是圍繞整個(gè)仿石鱸科的種類進(jìn)行，涉及胡椒鯛屬種類很少，且多數(shù)研究重點(diǎn)放在石鱸亞科，缺乏對(duì)胡椒鯛亞科種類的深入探討，對(duì)于胡椒鯛屬魚類較為全面的分子系統(tǒng)進(jìn)化研究至今還沒見報(bào)道。

ND1是NADH脫氫酶的第一個(gè)亞基，它可以從種的水平來推斷某個(gè)屬的系統(tǒng)發(fā)育關(guān)系［9］，是國(guó)內(nèi)外魚類系統(tǒng)分類研究廣泛應(yīng)用的一個(gè)分子標(biāo)記，應(yīng)用于魚類不同分類階元的系統(tǒng)發(fā)育研究，解決許多系統(tǒng)分類及進(jìn)化關(guān)系問題［10—12］。本實(shí)驗(yàn)基于ND1基因片段序列，在分子水平上對(duì)胡椒鯛屬及少棘胡椒鯛屬共16種魚類進(jìn)行系統(tǒng)發(fā)育分析，初步揭示胡椒鯛魚類的系統(tǒng)進(jìn)化關(guān)系，為仿石鱸科魚類的分類研究提供有效的分子信息。

2 材料與方法

2.1 材料來源

胡椒鯛實(shí)驗(yàn)樣品于2009年11月—2011年6月分別采集于廣東、海南等沿海城市水產(chǎn)市場(chǎng)。此外，部分樣品來自于西印度洋非洲沿海的國(guó)家，以及大洋洲島國(guó)新喀里多尼亞。主要來源于以下單位：西印度洋非洲沿海的樣品主要由南非水生物多樣性研究院（South African Institute for Aquatic Biodiversity，SAIAB）標(biāo)本管理員Bernard Makenzie教授提供以及美國(guó)堪薩斯大學(xué)Andrew Bentley教授提供；大洋洲的樣品由島國(guó)新喀里多尼亞努美阿潟湖水族館（Aquarium des Lagons）Jean-Lou Justine教授提供。最后共得到胡椒鯛屬魚類14種以及少棘胡椒鯛屬2種。實(shí)驗(yàn)材料均取肌肉或尾鰭，于95%的乙醇固定，用于實(shí)驗(yàn)室基因組DNA提取。系統(tǒng)進(jìn)化選取同為胡椒鯛亞科的磯鱸屬中的三線磯鱸Parapristipoma trilineatum作為外群的進(jìn)行分析。標(biāo)本種類與來源見表1。

表1 實(shí)驗(yàn)材料的種類和采集地Tab.1 Species and localities of materials

續(xù)表1

2.2 總基因組DNA的提取

取0.1 g肌肉樣品，加入500μl抽提液（0.3 mol／L NaCl、10 mmol／L EDTA，10 mmol／L Tris-Cl，p H 8.0）以及20 u L 10 mg／m L蛋白酶K，55℃消化3 h左右。消化產(chǎn)物經(jīng)飽和酚-氯仿抽提，2倍體積無水乙醇沉淀，70%乙醇洗滌，最后溶解于50 uL無菌水中，電泳檢測(cè)，－20℃保存?zhèn)溆谩?/p>

2.3 PCR擴(kuò)增和測(cè)序

用于擴(kuò)增ND1基因片段序列引物為：ND1＋5’-GTTGCMCAAACCATCTCHTAYGAA G-3’和ND1-5’-GGGGTATGGGCCCGATAGC-3’［13］。PCR反應(yīng)體系總體積為50μL，其中包括10×buffer 5μL，d NTPs（各2.5 mmol／L）2μL，上下游引物各1μL，Ex Taq酶（1 U／μL）2μ。PCR反應(yīng)條件為：94℃預(yù)變性5 min，94℃變性30 s，55℃退火1 min，72℃延伸1.5 min，35個(gè)循環(huán)，最后72℃延伸10 min。PCR產(chǎn)物用1.5%瓊脂糖凝膠電泳分離，切取特異性片段純化回收，純化產(chǎn)物送到上海英駿生物技術(shù)有限公司在ABI3730自動(dòng)測(cè)序儀上進(jìn)行雙向測(cè)序。

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

測(cè)序所得序列，利用Clustal W［14］進(jìn)行排序，生成序列矩陣去除兩端冗余，保留同源序列進(jìn)行系統(tǒng)發(fā)育分析。利用MEGA version 4.0［15］軟件對(duì)序列的堿基組成、堿基差異百分比、變異位點(diǎn)、簡(jiǎn)約信息位點(diǎn)及轉(zhuǎn)換／顛換值進(jìn)行統(tǒng)計(jì)，分子系統(tǒng)進(jìn)化樹利用最大簡(jiǎn)約法和最大似然法進(jìn)行構(gòu)建。系統(tǒng)樹分支的置信度采用重復(fù)抽樣分析（Bootstrap analysis）方法，重復(fù)抽樣的次數(shù)為1 000次。

3 結(jié)果分析

3.1 序列分析

測(cè)序得到16種胡椒鯛魚類的ND1部分序列長(zhǎng)度在698～712 bp之間，排列比對(duì)去除兩端冗余后共得到一致序列為673 bp。利用Mega4.0軟件分析可知，去除外群，16種胡椒鯛魚類共有保守位點(diǎn)492個(gè)，約占73.2%；變異位點(diǎn)181個(gè)，約占26.9%；簡(jiǎn)約性信息位點(diǎn)141個(gè)，約21.0%。T、C、A、G堿基平均含量分別為24.9%、31.0%、27.3%和16.8%，顯示G堿基相對(duì)缺乏，其中A＋T的含量（52.2%）高于G＋C含量（47.8%）。密碼子第1位4種堿基含量相差不大，其中C－1含量比其他3種堿基含量稍微偏高，占28.3%，密碼子第2位中，T－2含量高于其他3種，占35.0%，密碼子第3位中A－3，C－3含量最高，分別為34.1%及34.0%，而G－3含量最少，僅有12.3%，表現(xiàn)出明顯的反G偏倚（見圖1）這一模式與其他哺乳動(dòng)物相似［16］。

圖1 NADH脫氫酶亞基1（ND1）第一、二、三位密碼子堿基組成Fig.1 Base composition of the 1st，2nd，3rd codons of the dehydrogenase subunitⅠgene

序列中轉(zhuǎn)換與顛換之比值為2.61（Kimura 2-parameter模型），轉(zhuǎn)換明顯大于顛換，顯示這些位點(diǎn)沒有突變飽和。轉(zhuǎn)換和顛換主要發(fā)生在密碼子的第3位，其中轉(zhuǎn)換數(shù)為84.3%，而顛換為80%（表3）。這可能與密碼子第3位發(fā)生的堿基替換不易造成氨基酸突變，第1，第2位點(diǎn)的替換容易引起氨基酸突變有關(guān)。堿基替換中轉(zhuǎn)換以T－C轉(zhuǎn)換形式為主，T－C轉(zhuǎn)換多于A－G，顛換中A－T，A－C的顛換多于G－T，G－C（表2）。

表2 ND1序列堿基轉(zhuǎn)換和顛換數(shù)Tab.2 Numbers of the transitions and tranversions on the ND1 sequences

基于Kimura 2-parameter模型計(jì)算出16種胡椒鯛遺傳距離如表3所示，在內(nèi)類群中，16種魚類遺傳距離范圍在0.027～0.177之間。其中，兩種少棘胡椒鯛遺傳距離最小，序列中僅有幾個(gè)堿基的差異。其次是查氏胡椒鯛與紅唇胡椒鯛，為0.034；而東方胡椒鯛與肖氏胡椒鯛遺傳距離最大，為0.177，其次是四帶胡椒鯛與肖氏胡椒鯛，為0.173。同時(shí)，我們統(tǒng)計(jì)了少棘胡椒鯛屬與胡椒鯛屬的種類遺傳距離范圍在0.083～0.160之間，平均遺傳距離為0.112，少于胡椒鯛屬內(nèi)部種間的遺傳距離范圍0.034～0.177，平均遺傳距離0.123。

表3 基于Kimura 2-parameter模型計(jì)算16種胡椒鯛ND1序列遺傳距離（下三角）以及序列間的轉(zhuǎn)換／顛換數(shù)（上三角）Tab.3 Pairwise distances（below diagonal）and number of transitions／transversions（above diagonal）of 16 sweetlips for ND1 sequences based on Kimura 2-parameter model

續(xù)表3

3.2 分子系統(tǒng)進(jìn)化樹的構(gòu)建

基于所得的ND1部分序列，以三線磯鱸P.trilineatum作為外類群進(jìn)行系統(tǒng)進(jìn)化分析。利用最大簡(jiǎn)約法（MP）和最大似然法（ML）構(gòu)建的分子系統(tǒng)進(jìn)化樹。最后構(gòu)建得到MP和ML樹拓?fù)浣Y(jié)構(gòu)幾乎一致，僅在一些節(jié)點(diǎn)的支持率不同。其中最大簡(jiǎn)約樹樹長(zhǎng)（Tree Length，TL）為507，一致性指數(shù)（Consistency Index，CI））為0.484 4，保留指數(shù)（Retention Index，RI）為0.524 7。最大似然樹似然率ln L＝－3 425.46。這里把MP樹與ML樹合并（見圖2），樹上各分支的數(shù)字用1 000次Bootstrap表示統(tǒng)計(jì)分析后對(duì)該支的支持率（MP／ML）。

由圖可知，研究的16種胡椒鯛魚類中，在系統(tǒng)進(jìn)化樹上大致可分成兩個(gè)形態(tài)截然不同的類群。類群一包括斑胡椒鯛，條紋胡椒鯛等10種胡椒鯛，該類群的種類體表均有各種鮮艷的色彩和斑紋。類群二包括花尾胡椒鯛，黑胡椒鯛等6個(gè)種類，與第一類群不同的是，該類群的種類體色均暗淡單一，無鮮艷色彩與斑紋。同時(shí)，在進(jìn)化樹上，許多胡椒鯛都兩兩聚成姐妹種，而聚成姐妹種的胡椒鯛在外部形態(tài)上均有高度相似的特征。如：東方胡椒鯛與條紋胡椒鯛，體側(cè)均有明顯橫跨體表的黑帶，斑胡椒鯛與暗點(diǎn)胡椒鯛，體表分布著密集的小黑點(diǎn)。此外，我們發(fā)現(xiàn)少棘胡椒鯛屬的兩個(gè)種類都位于在類群一的胡椒鯛種類中，并與斑胡椒鯛，暗點(diǎn)胡椒鯛聚成的分支有比較近的親緣關(guān)系。系統(tǒng)進(jìn)化樹大部分節(jié)點(diǎn)都擁有較高的支持率。

4 討論

4.1 胡椒鯛屬魚類的系統(tǒng)進(jìn)化關(guān)系

本研究共采集到14種胡椒鯛和2種少棘胡椒鯛，其中7種主要采自南海，另外9種來自西印度洋非洲海域，幾乎囊括我國(guó)記錄的胡椒鯛種類［1—4］，部分種類甚至在我國(guó)并無分布［1—2，17—18］，具有一定的代表性?；贜D1基因部分序列構(gòu)建的系統(tǒng)發(fā)育樹中，16種魚類共形成兩個(gè)形態(tài)特征截然不同的的類群，類群Ⅰ包括斑胡椒鯛等7種胡椒鯛以及2種少棘胡椒鯛，類群Ⅱ包括花尾胡椒鯛等5種胡椒鯛。形態(tài)上，類群Ⅰ的胡椒鯛種類體表均具有各種鮮艷的色彩和斑紋，色彩斑紋隨著個(gè)體發(fā)育會(huì)產(chǎn)生變化。類群Ⅱ的胡椒鯛種類體色暗淡單一，無鮮艷色彩與斑紋，或僅有少量暗色斑紋，幼魚與成魚形態(tài)無明顯差異。在地理分布上，類群Ⅰ的胡椒鯛種類主要分布在西太平洋地區(qū)，類群Ⅱ的種類大多數(shù)僅分布于印度洋，少數(shù)在西太平洋也有分布［17—18］。

圖2 16種胡椒鯛魚類ND1基因部分序列利用最大簡(jiǎn)約法與最大似然法構(gòu)建的分子系統(tǒng)樹Fig.2 Molecular phylogenetic tree of 16 sweetlips based on partial ND1 gene sequences constructed by maximum parsimony and maximum likelihood method

形態(tài)學(xué)對(duì)于類群Ⅰ中胡椒鯛的分類存在較大爭(zhēng)議，皆因該類群近緣種類形態(tài)相似，體色與斑紋還會(huì)隨著個(gè)體發(fā)育產(chǎn)生巨大變化，如前言所述。這種形態(tài)差異多樣性給傳統(tǒng)形態(tài)分類鑒定帶來很大的困擾，許多物種分類不清。本研究構(gòu)建的聚類結(jié)果對(duì)上述胡椒鯛魚類的分類關(guān)系有較為清晰的展現(xiàn)。進(jìn)化樹上斑胡椒鯛與暗點(diǎn)胡椒鯛，四帶胡椒鯛與斜紋胡椒鯛，條紋胡椒鯛與東方胡椒鯛，以較高的支持率兩兩聚成姐妹分支。形態(tài)上，它們均擁有相似的形態(tài)特征，并且在個(gè)體發(fā)育過程中兩兩還具有幾乎一致的變化過程［2，17—18］。如斑胡椒鯛與暗點(diǎn)胡椒鯛，兩者在幼魚期體表具有棕色或黑色的斑紋，而成魚體表的斑紋逐漸斷裂成無數(shù)的小黑點(diǎn)。四帶胡椒鯛與斜紋胡椒鯛在幼魚期體表均有數(shù)條黑色水平縱帶，成魚后，四帶胡椒鯛體表保留為四條水平條帶，而斜紋胡椒鯛發(fā)育為許多斜行波狀條紋。對(duì)于條紋胡椒鯛與東方胡椒鯛，這兩種胡椒鯛幼魚體表都覆蓋著不規(guī)則斷裂狀的黑色斑紋，成魚后斑紋逐漸變成筆直的水平縱帶［2，17］。當(dāng)前形態(tài)分類資料對(duì)上述的胡椒鯛種類聚類關(guān)系不清，皆因各分類資料依據(jù)的分類標(biāo)準(zhǔn)不同。成慶泰和鄭葆珊［1］主編的魚類檢索系統(tǒng)主要根據(jù)成魚的形態(tài)進(jìn)行分類，Mc Kay［17］，Carpenter和Niem［18］，沈世杰［2］等分類資料卻結(jié)合了幼魚形態(tài)以及個(gè)體發(fā)育形態(tài)變化過程。Mc Kay［17］認(rèn)為，胡椒鯛屬許多種類的顏色斑紋在個(gè)體發(fā)育中差異巨大，形態(tài)分類需要根據(jù)幼魚，半成魚與成魚顏色斑紋，鰭棘與鰓耙數(shù)，才能較好的區(qū)分。本研究的系統(tǒng)進(jìn)化樹對(duì)上述胡椒鯛種類有比較清晰的分類關(guān)系。結(jié)果支持Mc Kay［17］分類觀點(diǎn)，認(rèn)為胡椒鯛魚類的分類與親緣關(guān)系的確定需要結(jié)合幼魚，成魚形態(tài)及發(fā)育過程才能真實(shí)反應(yīng)其分類關(guān)系。研究為部分形態(tài)分類混亂的胡椒鯛種類提供了有效的分子水平的分類證據(jù)。

4.2 少棘胡椒鯛屬的系統(tǒng)進(jìn)化地位

我國(guó)的大部分分類資料中，少棘胡椒鯛D.pictum一直被稱為胡椒鯛P.pictus，歸為胡椒鯛屬［1—2，19—20］，而許多國(guó)外的資料卻根據(jù)背鰭鰭棘條數(shù)而把它歸入到少棘胡椒鯛屬中，皆因其背鰭鰭棘數(shù)為9～10，而其他的胡椒鯛背鰭鰭棘都是12～14［17—18，21—23］。本研究共采集到兩種少棘胡椒鯛，一種是我國(guó)廣泛分布的少棘胡椒鯛D.pictum，另一種少棘胡椒鯛D.centurio采自非洲塞舌爾馬埃島。本研究結(jié)果顯示兩種少棘胡椒鯛緊密聚在一起，位于胡椒鯛屬的內(nèi)部。遺傳距離數(shù)據(jù)分析上，少棘胡椒鯛屬與胡椒鯛屬種類之間平均遺傳距離為0.112，少于胡椒鯛屬內(nèi)部的平均遺傳距離0.123。顯示少棘胡椒鯛與胡椒鯛屬的種類有很近的親緣關(guān)系。形態(tài)上，少棘胡椒鯛與其他色彩鮮艷的胡椒鯛種類一樣，體色和斑紋會(huì)隨著個(gè)體發(fā)育逐漸變化。兩種少棘胡椒鯛的幼魚體側(cè)有3條黑色縱帶，發(fā)育成魚后縱帶逐漸破裂成無數(shù)深色小斑點(diǎn)。把它歸為少棘胡椒鯛屬，可能僅僅是基于鰭棘數(shù)目的多少而得出的結(jié)論。但也有資料揭示少棘胡椒鯛屬除了鰭棘數(shù)與胡椒鯛屬不一樣之外，軟鰭條數(shù)目以及側(cè)線上鱗數(shù)也與胡椒鯛屬的不一致［9，18］?；贜D1序列的系統(tǒng)發(fā)育研究結(jié)果，在分子水平上，它們與胡椒鯛屬有著很近的親緣關(guān)系，這與當(dāng)今部分分子水平的分類研究結(jié)果一致。Sanciangco等［8］與Tavera等［9］曾利用多個(gè)分子標(biāo)記對(duì)北美仿石鱸科魚類的系統(tǒng)分類研究也發(fā)現(xiàn)少棘胡椒鯛位于胡椒鯛屬的內(nèi)部，并沒有單獨(dú)形成一支。Tavera認(rèn)為少棘胡椒鯛屬與胡椒鯛屬應(yīng)該是同一個(gè)屬。任崗［7］利用16SrRNA序列對(duì)12種石鱸科魚類系統(tǒng)進(jìn)化樹也揭示胡椒鯛P.pictus與胡椒鯛屬的種類聚在一起。因此，研究結(jié)果支持我國(guó)絕大多數(shù)傳統(tǒng)的分類資料，認(rèn)為少棘胡椒鯛屬的種類歸為胡椒鯛屬比較合適。

本研究通過利用線粒體DNA分子標(biāo)記ND1部分序列對(duì)胡椒鯛屬魚類進(jìn)行了系統(tǒng)進(jìn)化研究，初步闡明胡椒鯛屬內(nèi)種間的分類關(guān)系以及少棘胡椒鯛的分類問題，為仿石鱸科魚類的系統(tǒng)進(jìn)化研究提供新的分子信息。但要徹底弄清胡椒鯛魚類的進(jìn)化關(guān)系，還需以后更進(jìn)一步的研究。

［1］ 成慶泰，鄭葆珊.中國(guó)魚類系統(tǒng)檢索［M］.北京：科學(xué)出版社，1987：339－341.

［2］ 沈世杰.臺(tái)灣魚類志［M］.臺(tái)北：臺(tái)灣大學(xué)動(dòng)物學(xué)系，1993：360－363.

［3］ 黃宗國(guó).中國(guó)海洋生物種類與分布［M］.北京：海洋出版社，1994：708－709.

［4］ 孟慶聞，蘇錦祥，繆學(xué)祖.魚類分類學(xué)［M］.北京：中國(guó)農(nóng)業(yè)出版社，1995：636－640，696－700，711－713.

［5］ 朱世華，鄭文娟，鄒記興，等.5種石鱸科魚類細(xì)胞色素b基因序列及分子系統(tǒng)分析［J］.熱帶海洋學(xué)報(bào)，2006，25（4）：42－45.

［6］ 任崗，章群，錢開誠(chéng)，等.12種石鱸科魚類線粒體16Sr RNA基因的部分序列分析［J］.熱帶海洋學(xué)報(bào)，2007，26（3）：48－52.

［7］ Sanciangco M D，Rocha L A，Carpenter K E.A molecular phylogeny of the Grunts（Perciformes：Haemulidae）inferred using mitochondrial and nuclear genes［J］.Zootaxa，2011，2966：37－50.

［8］ Tavera J J，Acero A，Balart E F，et al.Molecular phylogeny of grunts（Teleostei，Haemulidae），with an emphasis on the ecology，evolution，and speciation history of New World species［J］.BMC Evolutionary Biology，2012，12：57.

［9］ Serb J M，Buhay J E，Lydeard C.Molecular systematics of the North American freshwater bivalve genus Quadrula（Unionidae：Ambleminae）based on mitochondrial ND1 sequences［J］.Molecular Phylogenetics and Evolution，2003，28（1）：1－11.

［10］ Utevsky S Y，Utevsky A Y，Schiaparelli S，et al.Molecular phylogeny of pontobdelline leeches and their place in the descent of fish leeches（Hirudinea，Piscicolidae）［J］.Zoologica Scripta，2007，36（3）：271－280.

［11］ Bowen B R，Kreiser B R，Mickle P F，et al.Phylogenetic relationships among North American Alosa species（Clupeidae）［J］.Journal of Fish Biology，2008，72（5）：1188－1201.

［12］ Thacker C E，Thompson A R，Roje D M.Phylogeny and evolution of Indo-Pacific shrimp-associated gobies（Gobiiformes：Gobiidae）［J］.Molecular Phylogenetics and Evolution，2011，59（1）：168－176.

［13］ Sorenson M D，Ast JC，Dimcheff D E，et al.Primers for a PCR-based approach to mitochondrial genome sequencing in birds and other vertebrates［J］.Molecular Phylogenetics and Evolution，1999，12（2）：105－114.

［14］ Thompson J D，Higgins D G，Gibson T J.CLUSTAL W：improving the sensitivity of progressive multiple sequence alignment through sequence weighting，position-specific gap penalties and weight matrix choice［J］.Nucleic Acids Research，1994，22（22）：4673－4680.

［15］ Tamura K，Dudley J，Nei M，et al.MEGA4：Molecular Evolutionary Genetics Analysis（MEGA）software version 4.0［J］.Molecular Biology and Evolution，2007，24（8）：1596－1599.

［16］ Hassanin A，Lecointre G，Tiller S.The“evolutionary signal”of homoplasy in proteincoding genesequences and its consequences for a priori weighting in phylogeny［J］.Comptes Rendus de l'Académie des Sciences-Series III-Sciences de la Vie，1998，321（7）：611－620.

［17］ Mckay R J.Haemulidae［M］／／Fischer W，Bianchi G.FAOspecies identification sheets for fishery purposes.Western Indian Ocean（Fishing Area 51），1984.

［18］ Carpenter K E，Niem V H.FAOSpecies Identification Guide for Fishery Purposes.The Living Marine Resources of the Western Central Pacific.Volume 5：Bony fishes part 3（Menidae to Pomacentridae）［M］.Rome：FAO，2001：2961－2989.

［19］ Kawamura T.Seasonal occurrence of fishes at inshore rocky reefs in Shirahama，southern Japan［J］.Japanese Journal of Ichthyology，1980，27（3）：243－248.

［20］ Kyushin K，Amaoka K，Nakaya K，et al.Fishes of the South China Sea［M］.Japan Marine Fishery Resource Research Center，Nori Otsuru，Tokyo，Japan，1982：103.

［21］ Randall JE，Lim K K P.A checklist of the fishes of the South China Sea［J］.The Raffles Bulletin of Zoology，2000，（Sup 8）：569－667.

［22］ Ramesh R，Nammalwar P.Database on coastal information of Tamilnadu［M］.Institute for Ocean Management Anna University Chennal－600 025，2008：1－132.

［23］ Shao K T，Ho H C，Lin P L，et al.A checklist of the fishes of southern Taiwan，Northern South China Sea［J］.The Raffles Bulletin of Zoology，2008，（Sup.19）：233－271.

Molecular phylogenetic relationships of 16 Sweetlips fishes based on the sequence of NADH dehydrogenase subunit I gene

Liang Rishen1，Liang Zhigang1，F(xiàn)ang Haohang1，Tan Jiayu1，Peng Yingxuan1，Xue Tao1
（1.Collegeof Life Science，Zhongkai University of Agriculture and Engineering，Guangzhou 510642，China）

In the present study，molecular phylogenetic relationships among 16 Sweetlips species collected from Indo-Western Pacific were analyzed based on partial DNA sequences of the molecular marker of mitochondrial NADH dehydrogenase subunit I（ND1）gene.The phylogenetic tree was constructed using maximum parsimony and maximum likelihood methods.Two morphologically different groups were formed in the tree.Group I was composed of sweetlips with wondrous coloration and diversed patterns，which were mainly found in the Western Pacific.Group II was composed of species with monotone color and uniformly dark patterns，which were mostly found in the area of Indian Ocean.Additionally，Diagramma species were all found positionally inside the Plectorhinchus tribe in group I in the phylogenetic tree.They were closely related with the sister species Plectorhinchus chaetodonoides and Plectorhinchuspicus.Based on the analysis of the genetic distance，the average value between Diagramma and Plectorhinchus was less than that among the interspecies of Plectorhinchus，indicating an extremely close relationship between Diagramma and Plectorhinchus.The result supported the opinion that Diagramma species might be classified into the genus Plectorhinchus.

phylogeny；Haemulidae；Plectorhinchus；Diagramma；ND1

Q959.46＋9

A

0253-4193（2014）10-0025-08

2013-09-30；

2013-11-16。

仲愷農(nóng)業(yè)工程學(xué)院大學(xué)生創(chuàng)新基金項(xiàng)目（2014A23）。

梁日深（1984—），男，廣東省開平市人，博士研究生，研究方向?yàn)樗镔Y源及保護(hù)利用。E-mail：cheetahliang＠126.com

梁日深，梁志剛，方浩航，等.基于ND1基因序列的16種胡椒鯛魚類分子系統(tǒng)進(jìn)化關(guān)系研究［J］.海洋學(xué)報(bào)，2014，36（10）：25—32，

10.3969／j.issn.0253-4193.2014.10.003

Liang Rishen，Liang Zhigang，F(xiàn)ang Haohang，et al.Molecular phylogenetic relationships of 16 Sweetlips fishes based on the sequence of NADH dehydrogenase subunit I gene［J］.Acta Oceanologica Sinica（in Chinese），2014，36（10）：25—32，doi：10.3969／j.issn.0253-4193.2014.10.003