王海山 ,蔣歡 ,陳治 ,楊超杰 ,葉樂

（1.海南熱帶海洋學(xué)院水產(chǎn)與生命學(xué)院，海南 三亞 572022；2.寧波大學(xué)海洋學(xué)院，浙江 寧波 315211）

鲹科（Carangidae）魚類隸屬硬骨魚綱（Osteichthyes）、輻鰭魚亞綱（Actinopteri）、鱸形目（Perciformes），下有4 個(gè)亞科：鲹亞科（Caranginae）、鰤亞科（Seriolinae）、鯧鲹亞科（Trachinotinae）和鰆鲹亞科（Chorineminae）。該科魚類廣泛分布于印度洋、太平洋、大西洋的熱帶和亞熱帶水域?，F(xiàn)記錄的鲹科魚類共34 屬約147 種，其中鲹亞科的種類最豐富，共20 屬，占全科的50%以上[1]。中國海域有鲹科魚類64 種，分布于東海、黃海以及南海，全部鲹科魚類在南海均有分布，東海分布約25 種。每到春夏季，魚群會(huì)隨暖流從低緯度區(qū)洄游到高緯度出生地產(chǎn)卵，也有少數(shù)幾種能夠進(jìn)入淡水流域。大部分鲹科魚類喜歡集群游動(dòng)于沿岸水域或珊瑚礁海域[2]，其中體型中等的種類游速較快。根據(jù)鲹科魚類體形特點(diǎn)和游動(dòng)時(shí)的運(yùn)動(dòng)學(xué)特征，常把其游動(dòng)模式運(yùn)用于水下推進(jìn)器的研究[3]。鲹科魚類多是重要的海洋經(jīng)濟(jì)魚類，如竹筴魚（Trachurus trachurus）、細(xì)鱗圓鲹（Decapterus macarellus）和長體圓鲹（Decapterus macrosoma）等是南海外海漁業(yè)的重要捕撈對(duì)象[4]，而五條鰤（Seriola quinqueradiata）、黃條鰤（Seriola lalandi）和卵形鯧鲹（Trachinotus ovatus）等高級(jí)食用海水魚，也是目前我國海南、廣東和廣西深水抗風(fēng)浪網(wǎng)箱養(yǎng)殖重要種類。為了促進(jìn)我國海水養(yǎng)殖業(yè)的發(fā)展，不同學(xué)者對(duì)多種鲹科魚類人工繁殖、網(wǎng)箱養(yǎng)殖技術(shù)、營養(yǎng)生理等展開了研究，并取得了突破性的進(jìn)展[5-7]。

關(guān)于魚類分類，目前也做了許多工作。朱元鼎和鄭文蓮（1958）[8]描述了在中國南海海域發(fā)現(xiàn)的鈍鰭葉鲹（Caranx pectoralis）和海南鰆鲹（Chorinemus hainanensis）兩個(gè)新種的外部形態(tài)；Suzuki（1962）[9]比較了日本海域不同鲹科魚類的骨骼系統(tǒng)，探討了鲹科魚類的系統(tǒng)問題；鄭文蓮（1981）[10,11]采用傳統(tǒng)形態(tài)測(cè)量和解剖學(xué)方法，從骨骼、耳石、鱗片、嗅覺器官、消化器官等系統(tǒng)研究了中國海域鲹科魚類分類、分布和形態(tài)，為推進(jìn)后續(xù)相關(guān)魚類科、亞科、屬、種的鑒別提供了重要參考資料；歐利國等（2020）[12,13]等基于地標(biāo)點(diǎn)法鑒定了南海海域4 種鲹科魚類的108 尾樣品，同時(shí)依據(jù)耳石形態(tài)特征聚類分析了南海東沙群島海域7 種鲹科魚類，并建立了系統(tǒng)分類檢索表，為南海鲹科魚種識(shí)別以及分類提供了參考資料?？傮w而言，從18 世紀(jì)60 年代開始，關(guān)于鲹科魚類分類是以新種描述、形態(tài)分類、生物地理分類為主。

隨著生物信息統(tǒng)計(jì)軟件以及線粒體DNA（mitochondrial DNA，mtDNA）特征被普遍應(yīng)用于各種生物的系統(tǒng)發(fā)育分析，分子生物學(xué)在鲹科魚類的研究中逐步開展。牛素芳等（2018）[14]分析了203 尾藍(lán)圓鲹的線粒體DNA 細(xì)胞色素b 序列的遺傳多樣性和遺傳結(jié)構(gòu)，研究了南海北部藍(lán)圓鲹的群體遺傳變異特征；許則灘等（2015）[15]探討了東海舟山群島普陀海域18 尾鲹科樣本的線粒體COI 序列，通過比較種間和種內(nèi)差異研究了本海域鲹科魚類的分類地位及分子系統(tǒng)進(jìn)化關(guān)系，得出了鲹科魚類同一屬內(nèi)種間遺傳距離明顯大于種內(nèi)遺傳距離的結(jié)論，并認(rèn)為普陀海域的18 尾鲹分別屬于7 屬8 種；朱世華等（2007）[16-18]結(jié)合數(shù)據(jù)庫已測(cè)定完成的3 種鲹科魚類序列，依據(jù)線粒體DNA 控制區(qū)、細(xì)胞色素b 序列、16S rRNA 序列生成矩陣和構(gòu)建的分子系統(tǒng)樹，分析了中國海域的8 屬9 種鲹科魚類，結(jié)果支持鲹科四個(gè)亞科的分類階元，也提出獅鼻鯧鲹（Trachinotus blochii）為卵形鯧鲹的誤鑒，以及將麗葉鲹（Caranx kalla）的中文名改為“麗副葉鲹”的觀點(diǎn)。該結(jié)論與國外魚類分類學(xué)者David 等（2002）[19]基于形態(tài)特征分析和分子數(shù)據(jù)構(gòu)建的部分鲹科魚類的分子發(fā)育樹相似。

魚類線粒體基因組不僅是研究真核細(xì)胞核酸與蛋白質(zhì)合成的模型系統(tǒng)，也是研究DNA 結(jié)構(gòu)與復(fù)制的良好模型。線粒體基因組為母性遺傳，進(jìn)化主要是以堿基替換為主，進(jìn)化速度較核DNA 快5～10 倍[20]，被廣泛應(yīng)用于魚類群體遺傳結(jié)構(gòu)分析和類群識(shí)別。目前，盡管16S rRNA、D-loop 區(qū)和COI序列等片段分別被運(yùn)用于鮭科、鯉科、鯛科和鳀科等魚類遺傳多樣性的研究中，但是，多認(rèn)為線粒體上不同的區(qū)域進(jìn)化速度不同，用序列更長的線粒體DNA 序列來構(gòu)建統(tǒng)發(fā)育樹可以減少信息分布的不均一性和序列長度太短而降低系統(tǒng)進(jìn)化樹的置信度帶來的影響，因此其結(jié)果更精確、可靠。

本文從GenBank 數(shù)據(jù)庫中下載已測(cè)定完成的33 種鲹科魚類的線粒體基因組序列，通過Clustal X和MEGA X 軟件進(jìn)行DNA 序列比對(duì)和核酸多態(tài)性、遺傳距離等分析，結(jié)合鄰接法（Neighbor-Joining，NJ）、最小進(jìn)化法（Minimum Evolution，ME）系統(tǒng)發(fā)育樹，探討鲹科種屬間的序列差異及其親緣關(guān)系，為鲹科的分類提供新的視角。

1 材料與方法

1.1 全基因組序列信息

通過檢索NCBI 的GenBank 數(shù)據(jù)庫獲取33 種鲹科魚類的線粒體基因組序列，并保存為FASTA格式。其種名、分類地位及相應(yīng)序列號(hào)等相關(guān)信息見表1。

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

將下載的基因組序列用Clustal X 軟件進(jìn)行序列比對(duì)，再輔以人工校對(duì)；運(yùn)用MEGA X 軟件統(tǒng)計(jì)序列的全長、堿基含量、變異位點(diǎn)數(shù)（比例）及簡約位點(diǎn)（Parsimony informative sites）等信息（表2）。多數(shù)學(xué)者認(rèn)為Kimura 雙參數(shù)模型（Kimura2-Parameter，K-2-P）是評(píng)估低遺傳變異序列的最佳核酸替代模型，本研究利用MEGA X 軟件中Kimura 雙參數(shù)法計(jì)算種屬間的遺傳距離；利用DnaSP 5.0 軟件統(tǒng)計(jì)33 種鲹科魚類線粒體基因組序列的堿基多態(tài)性；運(yùn)用MEGA 軟件中的鄰接法及最小進(jìn)化法構(gòu)建33 種鲹科魚類的系統(tǒng)進(jìn)化樹，采用Bootstrap 法檢驗(yàn)各分支的支持率，以自展值為1 000 次做重復(fù)抽樣檢驗(yàn)，其余參數(shù)為軟件默認(rèn)值。

2 結(jié)果與分析

2.1 核苷酸組成

序列分析顯示，鲹科魚類線粒體基因組序列全長在16 530～16 609 bp 之間，與大多數(shù)硬骨魚類線粒體序列長度相似[21]，其中鰤屬的高體鰤（Seriola dumerili）序列最短，黃條鰤的序列最長。核苷酸多樣性指數(shù)為37.3%，檢測(cè)到序列總的變異位點(diǎn)共7072個(gè)，占42.0%，簡約信息位點(diǎn)共6 297 個(gè)，變異位點(diǎn)數(shù)大于簡約信息位點(diǎn)數(shù)。轉(zhuǎn)換/ 顛換比值（R）1613/1026=1.6<2.0，說明基因序列的突變達(dá)到飽和；同屬的個(gè)體堿基差異為1.68%～12.89%，平均差異為6.75%，表明鲹科魚類種間的堿基組成差異很小，呈現(xiàn)一定的核酸組成穩(wěn)定性。線粒體全基因組堿基平均含量由高到低依次為C（29.5%）>A（27.8%）>T（26.0%）>G（16.7%），其中A+T 的含量（53.8%）大于G+C 含量（46.2%），表現(xiàn)出AT 含量相對(duì)豐富和GC含量相對(duì)缺乏的現(xiàn)象，符合硬骨魚類線粒體堿基組成偏好于A 和T 的特點(diǎn)[22-24]。本文所選全部物種堿基均有不同程度的AT 偏向（表1）。

DnaSP 5.0 計(jì)算位點(diǎn)多態(tài)性，發(fā)現(xiàn)了715 個(gè)單一變異位點(diǎn)、5 943 個(gè)簡約信息位點(diǎn)，其中發(fā)生種堿基突變的單一變異位點(diǎn)676 個(gè)、簡約信息位點(diǎn)3 120個(gè)，發(fā)生三種堿基突變的單一變異位點(diǎn)36 個(gè)、簡約信息位點(diǎn)1 363 個(gè)，發(fā)生四種堿基突變的單一變異位點(diǎn)3 個(gè)，簡約信息位點(diǎn)1 460 個(gè)。DnaSP 5.0 計(jì)算表明，33 種鲹科魚類線粒體基因組核酸多樣性（Pi）為0.15549，以100 bp 為單位在總長約16 800 bp 的線粒體基因組上步移，計(jì)算獲得核酸多樣性曲線（圖1），并結(jié)合鲹科魚類線粒體基因結(jié)構(gòu)圖分析可以發(fā)現(xiàn)，33 種鲹科魚類線粒體基因組在D-Loop 區(qū)域具有最高的核酸多樣性（第15 861～16 072 位堿基核酸多樣性指數(shù)Pi 為0.48081），其次為ND6 和ND5 區(qū)域（核酸多樣性指數(shù)Pi 分別為0.26080 和0.26313），ND2 和ND1 區(qū)域也具有較高的核酸多樣性指數(shù)（0.23428 和0.23282）。16S rRNA、tRNA 區(qū)域核酸多樣性較低（最低為0.005）。即在33 種鲹科魚類線粒體基因組上堿基突變頻率最高的區(qū)域出現(xiàn)在線粒體復(fù)制D-Loop 控制區(qū)和NADH 基因的亞基位置。16S、12S rRNA 和tRNA 位置堿基突變頻率較低，基因較為保守。

圖1 33 種魚線粒體基因組核酸多樣性曲線Fig.1 Mitochondrial DNA diversity of 33 fish species

2.2 遺傳距離矩陣

遺傳距離通常反映遺傳多樣性的高低，物種內(nèi)個(gè)體之間的遺傳距離大，表明個(gè)體之間的親緣關(guān)系遠(yuǎn)，物種間遺傳多樣性較高；相反的，物種內(nèi)個(gè)體之間的遺傳距離小，表明個(gè)體之間的親緣關(guān)系近，遺傳多樣性較低。K-2-P 計(jì)算的遺傳距離顯示，鲹科魚類種間遺傳距離平均值為0.183；種間遺傳最大值為卵形鯧鲹與高體若鲹（Carangoides equula）的0.250，最小值為黑紋小條鰤（Seriolina nigrofasciata）與白舌尾甲鲹（Uraspis helvola）的0.002。其中卵形鯧鲹與鯧鲹屬外物種的遺傳距離最大，平均值為0.229；北美鯧鲹（Trachinotus carolinus）與鯧鲹屬外物種的平均遺傳距離為0.228；布氏鯧鲹（Trachinotus blochii）與鯧鲹屬外物種的平均遺傳距離為0.226；泰勃圓鲹（Decapterus tabl）與其他鲹科魚類平均遺傳距離最小為0.170；細(xì)鱗圓鲹與其他鲹科魚類平均遺傳距離為0.172；絲鲹（Alectis ciliaris）與其他鲹科魚類平均遺傳距離為0.175。值得注意的是，高體若鲹與若鲹屬其他種的遺傳距離并不是最小的，其與甲若鲹（Carangoides armatus）、馬拉巴若鲹（Carangoides malabaricus）、竹筴魚、藍(lán)圓鲹、日本竹筴魚（Trachurus japonicus）和泰勃圓鲹的遺傳距離分別為0.184、0.179、0.156、0.161、0.159 和0.158；與細(xì)鱗圓鲹、珍鲹（Caranx ignobilis）、烏鲹（Parastromateus niger）的遺傳距離均小于與馬拉巴若鲹的遺傳距離（表2）。

表2 33 種魚線粒體基因序列的種間遺傳距離Tab.2 Genetic distance of 33 fish species based on mitochondrial genome

綜上，鯧鲹屬在進(jìn)化上與其他鲹科魚類分化較大，圓鲹屬與其他鲹科魚類遺傳距離更近，在線粒體基因組水平上高體若鲹與圓鲹屬魚類親緣關(guān)系更近，不符合傳統(tǒng)分類地位。

2.3 系統(tǒng)發(fā)育樹分析

早期對(duì)鲹科魚類的分類主要以外部形態(tài)為標(biāo)準(zhǔn)，由于地區(qū)、發(fā)育階段和環(huán)境的影響，相關(guān)研究中的特征描述或多或少存在差異，因而對(duì)鲹科魚類各屬種的分類存在爭(zhēng)議。本文分別采用NJ 法及ME法構(gòu)建了基于線粒體基因組核苷酸序列的33 種鲹科魚類系統(tǒng)進(jìn)化樹。Bootstrap 的參數(shù)值大于70%時(shí)，其結(jié)果符合真實(shí)的系統(tǒng)發(fā)育史的可靠性大于95%；如果參數(shù)值小于50%，則構(gòu)建的系統(tǒng)發(fā)育樹無法反映實(shí)際情況。本實(shí)驗(yàn)所建的NJ、ME 系統(tǒng)發(fā)育樹的絕大部分支的Bootstrap 都是100%，表明本研究可靠地反映了鲹科魚類的系統(tǒng)發(fā)育史（圖2、圖3）。

圖2 33 種魚基于線粒體序列的NJ 系統(tǒng)發(fā)育樹Fig.2 Phylogenetic trees of 33 fish species based on the complete mitochondrial genome by NJ analysis

圖3 基于線粒體序列的ME 系統(tǒng)發(fā)育樹Fig.3 Phylogenetic trees based on the complete mitochondrial genome by ME analysis

NJ 樹顯示：鲹科魚類分為3 大支系，從樹的基部第1 支系可知鲹亞科的脂眼凹肩鲹（Selar crumenophthalmus）是較早單獨(dú)分離出來的一支；第2支系是竹筴魚屬與圓鲹屬形成姐妹群，再與高體若鲹聚為一支（置信度為50%）；第3 支系由鰤亞科、鯧鲹亞科以及鲹亞科共同組成，較早分為2 個(gè)進(jìn)化分支，一個(gè)分支顯示鯧鲹亞科的卵形鯧鲹、布氏鯧鲹以及北美鯧鲹所組成一個(gè)進(jìn)化支，另一分支包含紡錘屬和鰤屬，再與鯧鲹屬聚集形成姐妹群。由此可知，鯧鲹亞科與鰤亞科有相當(dāng)密切的親緣關(guān)系；樹的頂部為鲹亞科的魚類，包含副葉鲹屬、葉鲹屬、細(xì)鲹屬、大甲鲹屬以及鲹屬，每個(gè)屬都單獨(dú)成為一個(gè)支系；泰利鲹（Caranxtille）和藍(lán)鰭鲹（Caranxmelampygus）、棉口尾甲鲹（Uraspis secunda）和白舌尾甲鲹、吉打副葉鲹（Alepes djedaba）和克氏副葉鲹（Alepes kleinii）、竹筴魚和日本竹筴魚、長鰭鰤（Seriola rivoliana）和高體鰤、卵形鯧鲹和布氏鯧鲹這六組物種遺傳距離相同，進(jìn)化速率相同；這種聚類關(guān)系與孫立元（2014）[25]的結(jié)果一致。

ME 樹顯示，鰤亞科進(jìn)化距離最長，位于系群的基部，是鲹科的原始種群，其次是鯧鲹亞科；鲹亞科又分化出3 個(gè)分支，聚集到樹枝末端——物種組成與其在NJ 樹上所處的分支基本相同，只在排列順序上有差異，表明兩種方法所得到的結(jié)果一致。

3 討論

魚類線粒體基因組具有分子小、進(jìn)化速度快、母系遺傳和信息相對(duì)穩(wěn)定等優(yōu)勢(shì)，適于從序列信息、基因組成及排列方式等角度來分析分子系統(tǒng)學(xué)及種群遺傳學(xué)[20]。本研究中，鲹科魚類線粒體基因組核苷酸序列的相關(guān)數(shù)據(jù)呈現(xiàn)AT 含量相對(duì)豐富和GC 含量相對(duì)缺乏的現(xiàn)象，符合硬骨魚類線粒體堿基組成偏好于A 和T 的特點(diǎn)。以NJ 法和ME 法構(gòu)建的鲹科3 亞科15 屬33 種魚類基因序列的系統(tǒng)進(jìn)化樹結(jié)果基本一致，均表明鯧鲹亞科與鰤亞科為姐妹群，這與朱世華等利用線粒體DNA 控制區(qū)、細(xì)胞色素b 序列、16S rRNA 序列分析結(jié)果一致。鰤亞科是3 個(gè)亞科中的最原始種群[10,11]。

早期的系統(tǒng)發(fā)育分析主要通過化石和現(xiàn)存物種的生理、形態(tài)特征來構(gòu)建系統(tǒng)發(fā)育關(guān)系。荷蘭自然學(xué)家Coenraad 和Hermann（1989）[26]首先于1844年根據(jù)從日本海域采集的高體若鲹標(biāo)本進(jìn)行了體長及背鰭、臀鰭、胸鰭和腹鰭數(shù)等特征的描述。1844—1907 年間，高體若鲹先后被更名為Carangoides acutus、Kaiwarinusequula、Carangoidesequula、Carangoides dasson 等。高體若鲹最開始被歸為鲹屬下的若鲹亞屬，后因分類系統(tǒng)修改鲹屬上升為鲹亞科，其下亞屬也相應(yīng)上升為屬，高體若鲹歸為鲹亞科的若鲹屬。目前高體若鲹為鲹亞科、若鲹屬的分類地位被廣為接受。較多文獻(xiàn)中竹筴魚、日本竹筴魚學(xué)名、中文名的使用也較為混亂。從物種的命名以及相關(guān)文獻(xiàn)的記錄都顯示，鲹科魚類中的某些種類分類地位存在爭(zhēng)議。本研究構(gòu)建的兩種系統(tǒng)發(fā)育樹均顯示：若鲹屬的馬拉巴若鲹、甲若鲹與絲鲹屬的絲鲹以及印度絲鲹（Alectis indica）形成姐妹群，而若鲹屬的高體若鲹并未與甲若鲹、馬拉巴若鲹聚為一支，而是與圓鲹屬、竹筴魚屬的種類聚為1 支（置信度84%～99%），表明高體若鲹與圓鲹屬、竹筴魚屬之間的親緣關(guān)系較近，與若鲹屬間的親緣關(guān)系較遠(yuǎn)。本研究結(jié)果與史寶等（2019）[27]在探討黃條鰤線粒體基因組系統(tǒng)進(jìn)化分析中所構(gòu)建的系統(tǒng)進(jìn)化樹結(jié)果基本一致，但后者并未討論高體若鲹處于的異常分支的結(jié)果。

以往分析鲹亞科時(shí)主要基于體形、肩帶下角有無深凹、上下頜有無牙齒、牙齒形狀、牙列數(shù)、第二背鰭有無延長鰭條等形態(tài)學(xué)特征進(jìn)行分類。如歐國利和劉必林（2020）[12]利用地標(biāo)點(diǎn)法對(duì)鲹科魚類耳石形態(tài)進(jìn)行分類識(shí)別，將高體若鲹與短吻絲鲹、竹筴魚以及脂眼凹肩鲹區(qū)分開來；從外部形態(tài)學(xué)分析，高體若鲹與馬拉巴若鲹、甲若鲹都具有體形為卵圓形、高而側(cè)扁的形態(tài)特征，但卻在一定程度上受性別、發(fā)育階段以及適應(yīng)環(huán)境的影響。如幼魚時(shí)，甲若鲹體側(cè)與高體若鲹同樣具六條暗色橫帶，但隨著個(gè)體成長，甲若鲹的暗色橫帶逐漸消失。因此，形態(tài)學(xué)鑒定在一定程度上受性別、發(fā)育階段以及適應(yīng)環(huán)境的影響，僅依靠形態(tài)學(xué)對(duì)物種分類鑒定存在著較多問題。