8種石首魚類線粒體基因組特征及分子系統(tǒng)進(jìn)化分析

孫利元, 楊天燕, 孟 瑋, 楊寶清, 張 濤

?

8種石首魚類線粒體基因組特征及分子系統(tǒng)進(jìn)化分析

孫利元1, 楊天燕3, 孟 瑋4, 楊寶清1, 張 濤2

(1. 山東省水生生物資源養(yǎng)護(hù)管理中心, 山東煙臺(tái)264003; 2. 中國科學(xué)院海洋研究所, 山東青島 266007; 3.中國海洋大學(xué)水產(chǎn)學(xué)院, 山東青島 266003; 4. 新疆大學(xué)生命科學(xué)與技術(shù)學(xué)院, 新疆烏魯木齊 830046)

為探討石首魚科(Sciaenidae)魚類分子系統(tǒng)進(jìn)化關(guān)系, 采用生物信息學(xué)方法分析了黑鰓梅童魚()、棘頭梅童魚()、大黃魚()、小黃魚()、鮸魚()、白姑魚()、黃姑魚()和皮氏叫姑魚()共8種石首魚類線粒體基因組全序列的基本特征。結(jié)果顯示, 除皮氏叫姑魚外, 其余7種石首魚類編碼的37個(gè)基因排列順序與脊椎動(dòng)物線粒體基因組相同?；蚪M堿基分布存在不均衡現(xiàn)象, A+T含量高于G+C含量。線粒體基因組的基因變異位點(diǎn)分析結(jié)果表明, ND4和ND5基因可作為基因的輔助分子標(biāo)記, 應(yīng)用于石首魚類群體遺傳學(xué)的研究中。黃魚亞科5種魚類13個(gè)蛋白質(zhì)編碼基因的比值遠(yuǎn)低于1, 顯示出較強(qiáng)的純化選擇。皮氏叫姑魚與其他石首魚間的遺傳距離均較大且親緣關(guān)系較遠(yuǎn), 暗示叫姑魚屬或?yàn)槭佐~類中較為原始的類群。基于線粒體基因組全序列構(gòu)建NJ系統(tǒng)樹支持黃魚亞科和白姑魚亞科親緣關(guān)系較近的形態(tài)學(xué)結(jié)論。而基于去除控制區(qū)后序列和13個(gè)蛋白質(zhì)編碼基因序列構(gòu)建的系統(tǒng)樹則表明兩亞科魚類間的差別在非編碼區(qū)更為明顯。

石首魚類; 線粒體基因組序列; 系統(tǒng)進(jìn)化

石首魚類隸屬鱸形目(Perciformes)、鱸亞目(Percoidei)、石首魚科(Sciaenidae), 是中國海洋經(jīng)濟(jì)魚類中產(chǎn)量最大的群類。中國沿海是石首魚類的重要分布區(qū)之一, 有17屬30種[1], 其中山東近海有6屬8種[2]。20世紀(jì)60~90年代, 老一輩魚類學(xué)家在對(duì)石首魚類外部形態(tài)特征研究的基礎(chǔ)上, 結(jié)合鰾和耳石的形態(tài)和式型的比較分析, 初步弄清了石首魚類的親緣隸屬關(guān)系, 也為中國魚類分類研究開辟了一條新途徑[3-5]。近年來, 隨著分子生物學(xué)理論和PCR技術(shù)的發(fā)展, 采用DNA分析的方法已逐漸成為研究和解決魚類系統(tǒng)發(fā)育問題的主要手段。如蒙子寧等[6]、張永等[7]和陳泉梅[8]基于線粒體16S rRNA序列分析了中國近海部分石首魚類分子系統(tǒng)進(jìn)化關(guān)系; 田蘭香等[9]通過比較長度為381 bp的Cyt基因序列, 探討了7種石首魚類的系統(tǒng)發(fā)育情況; 柳淑芳等[10]系統(tǒng)分析了19屬30種石首魚科魚類線粒體基因, 剖析了該基因在石首魚科魚類系統(tǒng)進(jìn)化研究中的應(yīng)用潛力; 馬春艷等[11]測(cè)定了核內(nèi)重組激活RAG-1基因序列, 并分析了中國沿海9屬13個(gè)種的石首魚間的差異。迄今尚未見運(yùn)用線粒體基因組全序列對(duì)山東近海石首魚類系統(tǒng)發(fā)育關(guān)系進(jìn)行分析研究。本研究選取山東近海分布的8種石首魚科魚類, 分析其線粒體基因組基本特征、蛋白質(zhì)編碼基因及系統(tǒng)發(fā)育關(guān)系等, 以期為山東近海石首魚類分子系統(tǒng)發(fā)育及種質(zhì)資源的研究提供遺傳信息, 同時(shí)也為尋找該類群合適的分子標(biāo)記提供參考資料。

1 材料和方法

1.1 數(shù)據(jù)獲取

從Genbank數(shù)據(jù)庫查詢并下載了《山東魚類志》[2]所記載的8種石首魚類線粒體DNA全序列, 這8種魚分別為梅童魚屬()的黑鰓梅童魚()和棘頭梅童魚(); 黃魚屬()的大黃魚()和小黃魚(); 鮸魚屬()的鮸魚(); 白姑魚屬()的白姑魚(); 黃姑魚屬()的黃姑魚(); 叫姑魚屬()的皮氏叫姑魚()。

1.2 數(shù)據(jù)處理

采用DNAstar 7.10和muscle3.8.31軟件對(duì)全部序列進(jìn)行多重比對(duì)編輯; 通過Dnasp 5.10.1軟件[12]計(jì)算序列堿基組成和變異位點(diǎn); PAMLX 1.3.1軟件[13]計(jì)算13個(gè)蛋白編碼基因的同義替換率(synonymous substitution rate,)和非同義替換率(non-synonymous substitution rate,); Mega 6.0軟件[14]基于蛋白質(zhì)編碼基因串聯(lián)后序列計(jì)算不同種間遺傳距離; PAUP 4.0軟件[15]基于Kimura雙因子參數(shù)模型, 分別對(duì)線粒體基因組全序列、去除控制區(qū)后的序列和13種蛋白質(zhì)編碼基因構(gòu)建鄰接系統(tǒng)進(jìn)化樹, 1 000次重抽樣檢驗(yàn)系統(tǒng)樹中節(jié)點(diǎn)的自引導(dǎo)值。

2 結(jié)果

2.1 線粒體基因組序列特征及差異位點(diǎn)分析

8種石首魚類線粒體DNA全序列長度在16 442 bp~ 19 154 bp, 其中棘頭梅童魚線粒體基因組長度最短, 皮氏叫姑魚最長。除皮氏叫姑魚外, 其余7種石首魚類線粒體基因組均包括了13個(gè)蛋白質(zhì)編碼基因、22個(gè)轉(zhuǎn)運(yùn)RNA基因和兩個(gè)核糖體RNA基因, 其排列順序與脊椎動(dòng)物線粒體基因組相同。分別計(jì)算8種石首魚類線粒體全基因組A、T、G、C堿基含量如表1所示, 除皮氏叫姑魚C堿基含量最低、T堿基含量最高以外, 其余7種魚類均表現(xiàn)為G堿基含量最低、C堿基含量最高。所有序列A+T含量均略高于G+C含量, 這與脊椎動(dòng)物偏好A、T堿基相一致。

分析比較了5種黃魚亞科(Lamichthyinae)魚類線粒體基因組13個(gè)蛋白質(zhì)編碼基因和2個(gè)核糖體RNA基因位點(diǎn)變異情況(表2), 可以看出15個(gè)基因變異程度在10.25%~30.12%。從基因的保守程度來看, 排在前3位的是12S rRNA、16S rRNA和COIII基因, 其變異位點(diǎn)比例分別為10.25%、11.46%和17.71%。從變異位點(diǎn)數(shù)來看, ND5的變異位點(diǎn)數(shù)最多, 為427個(gè), COI次之, 為399個(gè)。

表1 8種石首魚類線粒體基因組序列特征

2.2 蛋白質(zhì)編碼基因特征分析

8種石首魚類的13個(gè)蛋白質(zhì)編碼基因中, 除ND6由輕鏈編碼, 其余均由重鏈編碼。采用Mega軟件分析其密碼子組成情況(表3), 除、、、和基因編碼氨基酸的數(shù)量在不同石首魚類中存在差異以外, 其余8個(gè)蛋白質(zhì)編碼基因所編碼氨基酸的數(shù)目完全相同, 分別為ND2(348)、ND3(116)、ND4L(99)、ND4(460)、COII(230)、COIII(261)、ATP6(228)、ATP8(56)。8種石首魚類的13個(gè)蛋白編碼基因中, 除黃姑魚ND4基因和白姑魚ND5基因以GTG為起始密碼子, 其余均以ATG作為起始密碼子。、、、、和基因出現(xiàn)了不完全終止密碼子, 其中ND2基因除黃姑魚和皮氏叫姑魚以T堿基結(jié)束, 其他均為TA結(jié)束。

表3 蛋白質(zhì)編碼基因的氨基酸數(shù)量和起始、終止密碼子

2.3 選擇壓力分析

進(jìn)化速率通常受到穩(wěn)定性選擇、突變和定向選擇控制。一般而言, 同義突變不受自然選擇作用, 而非同義突變則相反。為了檢驗(yàn)石首魚類線粒體基因組選擇壓力, 選取了黃魚亞科的黑鰓梅童魚、棘頭梅童魚、大黃魚、小黃魚和鮸魚為代表, 估算了線粒體基因組中13個(gè)蛋白質(zhì)編碼基因的非同義替換率和同義替換率, 并構(gòu)建二者比值的柱狀圖如圖1所示。13種蛋白編碼基因的值在0.0234~0.2338之間且均小于1, 顯示出較強(qiáng)的負(fù)(純化)選擇。

2.4 遺傳距離及系統(tǒng)發(fā)育關(guān)系分析

將13個(gè)蛋白編碼基因序列首尾串聯(lián)后, 利用Mega軟件計(jì)算Kimura 2-Parameter遺傳距離, 從表4可以看出, 8種石首魚類兩兩間平均遺傳距離為0.239, 其中皮氏叫姑魚和白姑魚遺傳距離最大, 為0.409, 小黃魚和黑鰓梅童魚遺傳距離最小, 為0.006。計(jì)算3個(gè)亞科間的遺傳距離顯示, 黃魚亞科和白姑魚亞科遺傳距離較小, 為0.230; 叫姑魚亞科和白姑魚亞科、黃魚亞科遺傳距離分別為0.399和0.402。黃魚亞科魚類平均遺傳距離為0.130(0.006~0.177); 白姑魚亞科為0.243。

選取鱸形目、攀鱸科(Anabantidae)的攀鱸()為外群, 分別基于整個(gè)線粒體基因組全序列、去除控制區(qū)后的序列和13種蛋白質(zhì)編碼基因構(gòu)建NJ系統(tǒng)發(fā)育樹如圖2所示。前兩種聚類樹中, 黃魚亞科的5種魚類均以相同的順序?qū)訉泳鄣揭黄? 且皮氏叫姑魚和外群聚類并出現(xiàn)在聚類樹的根部, 但基于整個(gè)線粒體基因組序列構(gòu)建的系統(tǒng)樹中白姑魚亞科的白姑魚和黃姑魚聚為一支后, 再與黃魚亞科魚類聚類; 而基于去除控制區(qū)后的序列構(gòu)建的系統(tǒng)樹中白姑魚和黃姑魚分別與黃魚亞科魚類聚類。去除控制區(qū)后的序列和13種蛋白質(zhì)編碼基因構(gòu)建的系統(tǒng)樹幾乎完全一致, 僅部分節(jié)點(diǎn)的支持率存在差別。

3 討論

由于DNA雙鏈之間存在選擇壓力及自然突變的差異, 基因組上堿基的分布往往存在著不均衡的現(xiàn)象[16]。本研究中8種石首魚類線粒體基因組A+T含量均略高于G+C含量, 呈現(xiàn)出明顯的堿基偏好性。Brown等[17]通過研究發(fā)現(xiàn), 線粒體DNA 2/3復(fù)制周期時(shí)間內(nèi)都在H-鏈上, 形成沒有保護(hù)的單鏈狀態(tài), 更加容易產(chǎn)生氧化和水解。H鏈和L鏈的突變壓力不同直接影響了堿基組成的不對(duì)稱。此外, 皮氏叫姑魚基因組序列較其他7種魚類長2 655 bp~2 712 bp不等, 通過同源序列的比對(duì)發(fā)現(xiàn), 其線粒體基因組全序列包含5個(gè)非編碼區(qū), 且tRNAVal、12S rRNA、16S rRNA和tRNAPhe這4個(gè)基因的排列順序與其他石首魚類有變化[16]。

基因的非同義替換率與同義替換率的比較是分子進(jìn)化研究的重要內(nèi)容[18]。根據(jù)中性進(jìn)化理論[19-20], 由于同義替換并不影響編碼氨基酸的變化, 考察非同義替換率與同義替換率的比率, 可用來判斷是否有選擇壓力作用于這個(gè)蛋白質(zhì)編碼基因。當(dāng)比值大于1時(shí), 認(rèn)為有正選擇效應(yīng); 當(dāng)該值等于1時(shí), 認(rèn)為存在中性選擇; 當(dāng)該值小于1時(shí), 認(rèn)為有凈化選擇[21]。本研究中5種黃魚亞科魚類的13個(gè)蛋白質(zhì)編碼基因,比值均遠(yuǎn)小于1, 表明這些基因受到凈化選擇壓力, 其非同義替換速率比同義替換速率低很多。基因的均值最低, 僅為0.0234, 說明其承受的選擇壓力和功能束縛最為強(qiáng)烈。、和3個(gè)基因的均值最高, 分別為0.2338、0.1328和0.125, 其承受的選擇壓力較弱。

a. 基于線粒體基因組全序列; b. 基于去除控制區(qū)后序列; c. 基于線粒體基因組13個(gè)蛋白質(zhì)編碼基因序列

a. Based on the complete mitochondrial genome; b. Based on the sequences with the control region deleted; c. Based on 13 protein-coding genes

線粒體基因序列通常作為DNA條形碼, 廣泛應(yīng)用于物種分類和鑒定[22-25]。從本研究結(jié)果看出在15個(gè)基因中,基因變異位點(diǎn)數(shù)為399個(gè), 僅次于(427)、略高于(314)。選擇合適的分子標(biāo)記對(duì)于種群遺傳的研究至關(guān)重要, 綜合上述研究結(jié)果, 推測(cè)和基因或可作為基因的輔助備選分子標(biāo)記, 應(yīng)用于石首魚類群體遺傳學(xué)的研究中。

通過對(duì)遺傳距離統(tǒng)計(jì)發(fā)現(xiàn), 皮氏叫姑魚與其他石首魚間的遺傳距離均較大, 表明叫姑魚亞科與其他石首魚類的親緣關(guān)系較遠(yuǎn)。蒙子寧等[6]和田蘭香等[9]分別采用16SrRNA和基因?qū)κ佐~類系統(tǒng)發(fā)育的研究中也得出了類似結(jié)論。從系統(tǒng)發(fā)育關(guān)系來看, 3種NJ樹都顯示皮氏叫姑魚單獨(dú)聚類于系統(tǒng)樹的基部, 該結(jié)果暗示叫姑魚屬或?yàn)槭佐~類中較為原始的類群, 與朱元鼎等[3]將叫姑魚亞科作為石首魚類中最先分化的一支的分類相一致。

朱元鼎等[3]和成慶泰等[4]將黃魚屬、梅童魚屬、鮸魚屬歸屬于黃魚亞科, 認(rèn)為其耳石形態(tài)與白姑魚亞科相近, 似與白姑魚屬的關(guān)系更為密切。本研究基于線粒體基因組全序列構(gòu)建系統(tǒng)樹中, 黃魚亞科的5種魚類聚到一起后又與白姑魚亞科的2種魚類相聚, 該結(jié)果支持這兩個(gè)亞科親緣關(guān)系較近的形態(tài)學(xué)結(jié)論。而基于去除控制區(qū)后序列和13個(gè)蛋白質(zhì)編碼基因序列構(gòu)建的系統(tǒng)樹也得到了類似的結(jié)果, 兩個(gè)亞科的7種魚類聚成了一支, 略有不同的是兩亞科間的支持率不高, 表明兩亞科魚類間的差別在非編碼區(qū)更為明顯。

[1] 劉效舜. 黃渤海區(qū)漁業(yè)資源調(diào)查與區(qū)劃[M]. 北京: 海洋出版社, 1990: 191-200. Liu Xiaoshun. Fishery resources investigation and regionalization of Bohai and Yellow Sea[M]. Beijing: Ocean Press, 1990: 191-200.

[2] 成慶泰, 周才武. 山東魚類志[M]. 濟(jì)南: 山東科學(xué)技術(shù)出版社, 1997: 282-293. Cheng Qingtai, Zhou Caiwu. Fish fauna of Shandong[M]. Jinan: Shandong Science and Technology Press, 1997: 282-293.

[3] 朱元鼎, 羅云林, 伍漢霖. 中國石首魚類分類系統(tǒng)的研究和新屬新種的敘述[M]. 上海: 上?？茖W(xué)技術(shù)出版社, 1963. Zhu Yuanding, Luo Yunlin, Wu Hanlin. Sciaenid fishes classification system research and new species description in China[M]. Shanghai: Shanghai Science and Technology Press, 1963.

[4] 成慶泰, 鄭葆珊. 中國魚類系統(tǒng)檢索(上冊(cè))[M]. 北京: 科學(xué)出版社, 1987: 317-324. Cheng Qingtai, Zheng Baoshan. Chinese fish retrieval system (Vol.1.) [M]. Beijing: Science Press, 1987: 317- 324.

[5] 孟慶聞, 陳立行. 魚類分類學(xué)[M]. 北京: 中國農(nóng)業(yè)出版社, 1995: 713-728. Meng Qingwen, Chen Lixing. Fish Taxonomy[M]. Beijing: China Agriculture Press, 1995: 713-728.

[6] 蒙子寧, 莊志猛, 丁少雄, 等. 中國近海8種石首魚類的線粒體16S rRNA基因序列變異及其分子系統(tǒng)進(jìn)化[J].自然科學(xué)進(jìn)展, 2004, 14(5): 514-521. Meng Zining, Zhuang Zhimeng, Ding Shaoxiong, et al. The analysis of mitochondrial 16S rRNA gene sequence variation and molecular system evolution of eight sciaenid fishes in Chinese offshore[J]. Progress in Natural Science, 2004, 14(5): 514-521.

[7] 張永, 馬春艷, 馬凌波, 等. 基于16S rRNA部分序列探討中國近海十三種石首魚類的分子系統(tǒng)進(jìn)化關(guān)系[J]. 海洋漁業(yè), 2010, 32(3): 276-282. Zhang Yong, Ma Chunyan, Ma Lingbo, et al. Molecular phylogenetic relationships of 13 Sciaenidae species in China Sea areas based on 16SrRNA fragement sequences[J]. Marine Fisheries, 2010, 32(3): 276-282.

[8] 陳泉梅. 中國石首魚科魚類分子系統(tǒng)學(xué)研究[D]. 廣州: 暨南大學(xué), 2007. Chen Quanmei. Molecular phylogeny of the Sciaenidae in China[D]. Guangzhou: Jinan University, 2007.

[9] 田蘭香, 梁冰, 張樹義, 等. 細(xì)胞色素b基因序列與7種石首魚類的系統(tǒng)進(jìn)化[J]. 臺(tái)灣海峽, 2004, 23(4): 436-443. Tian Lanxiang, Liang Bing, Zhang Shuyi, et al. Phylogenetic relationships of 7 Sciaenidae species based on cytochrome b gene sequences[J]. Journal of Oceanography in Taiwan Strait, 2004, 23(4): 436-443.

[10] 柳淑芳, 陳亮亮, 戴芳群, 等. 基于線粒體基因的DNA條形碼在石首魚科(Sciaenidae)魚類系統(tǒng)分類中的應(yīng)用[J]. 海洋與湖沼, 2010, 41(2): 223-232. Liu Shufang, Chen Liangliang, Dai Fangqun, et al. Application of DNA barcoding genefor classifying family Sciaenidae[J].Oceanologia et Limnologia Sinica, 2010, 41(2): 223-232.

[11] 馬春艷, 馬凌波, 倪勇, 等. 基于RAG1基因的中國近海13種石首魚科魚類系統(tǒng)進(jìn)化關(guān)系[J]. 水產(chǎn)學(xué)報(bào), 2012, 36(1): 9-16. Ma Chunyan, Ma Lingbo, Ni Yong, et al. Molecular phyologenetic relationships of 13 Sciaenidae species in the China Sea based on RAG1 gene sequences[J]. Journal of Fisheries of China, 2012, 36(1): 9-16.

[12] Rozas J. DNA sequence polymorphism analysis using DnaSP. bioinformatics for DNA sequence analysis, methods in molecular biology[M]. New Jersey: Humana Press, 2009: 337-350.

[13] Xu B, Yang Z. PAMLX: a graphical user interface for PAML[J]. Molecular Biology and Evolution, 2013, 30(12): 2723-2724.

[14] Tamura K, Stecher G, Peterson D, et al. MEGA6: molecular evolutionary genetics analysis version 6.0[J]. Molecular Biology and Evolution, 2013, 30(12): 2725- 2729.

[15] Swofford D L. Paup 4.0 for Unix/Vms: phylogenetic analysis using parsimony[M]. Sunderland: Sinauer Associates Inc, 1999.

[16] 鐘東, 趙貴軍, 張振書, 等. 基因組內(nèi)堿基分布整體均衡與局部不均衡的研究進(jìn)展[J]. 遺傳, 2002, 24(3): 351-355. Zhong Dong, Zhao Guijun, Zhang Zhenshu, et al. Advance in the entir balance and local unbalance of base distribution in genome[J]. Hereditas, 2002, 24(3): 351-355.

[17] Brown W M, Prager E M, Wang A, et al. Mitochondrial DNA sequences of primates: tenpo and mode of evohxtion[J]. Journal of Molecular Evolution, 1982, 18: 225-239.

[18] Xu T J, Tang D, Jin X X. A surprising arrangement pattern and phylogenetic consideration: the complete mitochondrial genome of Belanger’s croaker Johnius belangerii (Percoidei: Sciaenidae)[J]. Mitochondrial DNA, 2015, 26(4): 655-657.

[19] Nei M, Kumar S. Molecular evolution and phylogenetics[M]. New York: Oxford University Press, 2000: 252-265.

[20] Kimura M. The neutral theory of molecular evolution[M]. Cambridge: Cambridge University Press, 1983.

[21] Ohta T. The nearly neutral theory of molecular evolution[J]. Annual Review of Ecology and Systematics, 1992, 23: 263-286.

[22] Yang Z, Bielawski J P. Statistical methods for detecting molecular adaptation[J]. Trends in Ecology and Evolution, 2000, 15(12): 496-503.

[23] 裴男才, 陳步峰. 生物DNA條形碼: 十年發(fā)展歷程、研究尺度和功能[J]. 生物多樣性, 2013, 21(5): 616- 627. Pei Nancai, Chen Bufeng. DNA barcoding of life: a classification of uses according to function and scale after ten years of development[J]. Biodiversity Science, 2013, 21(5): 616-627.

[24] 莫幫輝, 屈莉, 韓松, 等. DNA條形碼識(shí)別I. DNA條形碼研究進(jìn)展及應(yīng)用前景[J]. 四川動(dòng)物, 2008, 27(2): 303-306. Mo Banghui, Qu Li, Han Song, et al. DNA barcoding identification I. research progress and applied perspective of DNA barcoding[J]. Sichuan Journal of Zoology, 2008, 27(2): 303-306.

[25] 丁蘭平, 馬元元, 黃冰心. DNA條形碼技術(shù)在大型海藻學(xué)研究中的應(yīng)用及前景[J]. 海洋科學(xué), 2012, 36(11): 103-110. Ding Lanping, Ma Yuanyuan, Huang Bingxin. The application and perspective of DNA barcoding technology on the macroalgae[J]. Marine Sciences, 2012, 36 (11): 103-110.

[26] 張輝, 線薇薇. DNA條形碼在長江口魚類浮游生物生態(tài)學(xué)研究中的意義[J]. 海洋科學(xué), 2015, 39(4): 135- 137. Zhang Hui, Xian Weiwei. The importance of DNA barcoding on the ichthyoplankton ecological study in the Yangtze estuary[J]. Marine Sciences, 2015, 39(4): 135-137.

(本文編輯: 譚雪靜)

Analysis of the mitochondrial genome characteristics and phylogenetic relationships of eight sciaenid fishes

SUN Li-yuan1, YANG Tian-yan3, MENG Wei4, YANG Bao-qing1, ZHANG Tao2

(1. Shandong Hydrobios Resources Conservation and Management Center, Yantai 264003, China; 2. Institute of Oceanology, Chinese Academy of Sciences, Qingdao 266007, China; 3. Fisheries College, Ocean University of China, Qingdao 266003, China; 4. College of Life Science and Technology, Xinjiang University, Urumqi 830046, China)

To explore the phylogenetic relationships of sciaenid fishes, the basic characteristics of mitochondrial genomes in eight sciaenid fishes (,,,,,,and) were revealed by comprehensive bioinformatic analysis. The results showed that mitochondrial genomes contained 37 standard vertebrate genes, and the order was identical except for inThe base distribution was not balanced, and the content of A+T was higher than that of G+C. Genetic variation analysis of 13 protein-coding genes and 2 ribosomal RNA genes among sciaenid fishes showed that the ND4 and ND5 genes could be used as supplementary molecular markers to the COI gene and applied to research on population genetic diversity. Theratios of 13 mitochondrial protein-coding genes in five Pseudosciaeninae fishes were much lower than 1, indicating strong purifying selection. The genetic distances betweenand other sciaenid fishes were high, and relationships were distant, which suggested thatfishes might be the relatively primitive group. A neighbor-joining phylogenetic tree based on the complete mitochondrial genome supported the morphological conclusion that fishes of Pseudosciaeninae and Argyrosominae were closely related. Furthermore, phylogenetic trees based on sequences deleting the control region and 13 protein-coding genes showed that the differences between these two subfamilies were much more evident.

Sciaenid fishes; mitochondrial genome sequence; phylogenetic

Apr. 22, 2016

[Department of Ocean and Fisheries of Shandong Province Program, No. 270006-FZLX-2015-00-1]

S917.4

A

1000-3096(2017)03-0048-07

10.11759//hykx20160422005

2016-04-22;

2016-08-03

山東省海洋與漁業(yè)廳資助項(xiàng)目(No. 270006-FZLX-2015-00-1)

孫利元(1980-), 男, 山東招遠(yuǎn)人, 工程師, 主要從事海洋漁業(yè)資源學(xué)研究, E-mail: heroland80@163.com