渤海灣7種蝦虎魚科魚類COI基因序列特征及分子分類研究

房恩軍　郭彪　鄭德斌　張 雪　王碩　王宇

摘? ? 要：為獲得渤海灣天津海域常見蝦虎魚線粒體條形碼序列及系統分類信息，采集渤海灣7種蝦虎魚樣本，通過形態初步鑒定物種，并提取DNA進行線粒體COI基因部分序列擴增測序，用系統進化分析及分子鑒定軟件MEGA對渤海灣常見12種蝦虎魚COI基因部分片斷進行序列特征分析及分子分類分析。擴增得到7種蝦虎魚COI基因片斷，分析得到其平均GC含量低于AT含量，堿基在密碼子位置分布中具有偏向性，且GC含量最高為第2位堿基，與多種魚類COI基因堿基特點相似。基因密碼子堿基變異率為19.4%，轉換顛換比為1.3。遺傳距離分析及系統進化樹分析理清了渤海灣天津海域12種主要蝦虎魚的系統分類及進化關系。遺傳距離顯示，矛尾刺蝦虎魚（Amblychaeturichthys hexanema）和斑尾復蝦虎魚（Synechogobius ommaturus）為同種魚。蝦虎魚COI編碼基因密碼子高變異率可能與其良好的環境適應性有關。

關鍵詞：渤海灣蝦虎魚；線粒體細胞色素c氧化酶亞基（COI）；系統進化

中圖分類號：Q917.4? ? ? ? ?文獻標識碼：A? ? ? ? ? DOI 編碼：10.3969/j.issn.1006－6500.2023.09.006

Sequence and Molecular Classification of COI Gene of Seven Species of Goby from the Bohai Bay

FANG Enjun，GUO Biao，ZHENG Debin，ZHANG Xue，WANG Shuo，WANG Yu

（Tianjin Fisheries Research Institute，? Center for Marine Ranching Engineering Science Research of Tianjin，Tianjin 300221， China）

Abstract： In order to clarify the mitochondrial barcode sequence and taxonomic information of common goby fish in Tianjin area of Bohai Bay， seven species of goby fish were collected from Bohai Bay， identified by morphology， and DNA was extracted for amplification and sequencing of mitochondrial COI gene fragments. Phylogenetic analysis and molecular identification software MEGA were used to analyze the sequence characteristics and molecular classification of COI gene fragments of 12 common goby species in Bohai Bay. Seven COI gene fragments of Goby were amplified and sequenced， the average GC content was lower than AT content， the base had a bias in the codon position distribution， and the highest GC content was the second base， which was similar to the COI gene base of many fishes. The codon base variation rate was 19.4%， and the ratio of transitions to transversions was 1.3. Phylogenetic distance analysis and phylogenetic tree analysis clarified the phylotaxonomic and evolutionary relationships of 12 major goby species in Tianjin area of Bohai Bay. Genetic distance indicated that Amblychaeturichthys hexanema and Synechogobius ommaturus were the same species. The high codon variation rate of codon encoding gene in Goby may be related to its good environmental adaptability.

Key words： Goby of the Bohai Bay; cytochrome oxidase subunit I （COI） gene; phylogenetic analysis

渤海灣是一個淺水淤泥質半封閉海灣，歷史上餌料豐富，基礎條件優越，是渤海魚類的產卵場。近年來，由于圍填海工程、水域污染和過度捕撈等因素的影響，漁業資源的生存環境惡化，影響了其生態系統結構和功能，漁業資源衰退，小型化、低值化問題突出，種類和數量也在減少。年度調查表明，天津海域產卵場逐漸退化，與歷史調查數據相比，渤海灣魚類浮游生物的豐度和多樣性顯著減少，經濟魚類漁獲和產卵量也逐漸減少[1-2]。而蝦虎魚由于環境適應力強，繁殖力強和生命周期短等特點[3]，其資源量穩定增加，已逐漸代替其他經濟魚類，成為渤海灣天津海域春季休漁期前魚類群落及仔稚魚的優勢種，在海洋食物網中發揮著重要的作用。但是蝦虎魚種類繁多，且有關仔稚魚的形態相近，相關資料較少，鑒定難度大，準確性難以保證，蝦虎魚仔稚魚鑒定也是世界難點和是魚類早期資源研究的重要內容。

20世紀90年代以來，國內外科學家發現線粒體一些短的基因（如細胞色素氧化酶第一亞基（COI）基因、12s rRNA、細胞色素B（Cty B）等）進化速率快且相對保守，其中COI基因5'端650 bp片斷可有效區分種間物種，可作為物種分類的DNA條形碼（DNA Barcoding），隨后與之配套的數據庫也日趨完善，廣泛應用于生命科學、法醫學、食品質量控制等領域，也為魚類浮游生物種類的分子鑒定奠定了良好的基礎[4-9]。

蝦虎魚亞目（Gobioidei）有兩千余種蝦虎魚，腹鰭合成一個吸盤，是鱸形目中最大的一個類群。蝦虎魚多為小型底棲魚類，主要以小型底棲甲殼動物為食，在渤海灣生態系統中能量流動和物質循環中有重要作用，是底層食物鏈的關鍵物種[10]。根據渤海魚類生物學記載， 渤海灣天津海域歷史上能采集到的蝦虎魚有12種左右[11]，其中拉氏狼牙蝦虎魚（Odontamblyopus rubicundus）、斑尾腹蝦虎魚（Synechogobius ommaturus）具有一定的經濟價值。除了小頭櫛孔蝦虎魚。常見的蝦虎魚的線粒體（COI）基因片斷均能在數據庫中獲取。但是蝦虎魚形態多變、各地名稱不統一等原因，需要通過擴增進一步明確序列信息。本研究選擇渤海灣常見7種蝦虎魚類作為研究對象，通過獲取線粒體（COI）基因序列分析比較其序列變異情況，進一步進行系統發育及遺傳分析，從而為魚類DNA條形碼庫及魚類浮游生物的鑒定積累基礎數據，為渤海海洋食物網的研究和穩定，為渤海灣蝦虎魚資源保護和合理開發利用提供科學依據。

1 材料與方法

1.1 樣品采集

樣品來源于2018年春季和夏季渤海灣水域地籠或拖網采樣，參照《渤海魚類生物學》共鑒定出矛尾刺蝦虎魚（Amblychaeturichthys hexanema）、裸頸縞蝦虎魚（Tridentiger Nudicervicus）、矛尾蝦虎魚（Chaeturichthys stigmatias）、紅狼牙蝦虎魚（Odontamblyopus rubicundus）、彈涂魚（Periophthalmus modestus）、鐘馗蝦虎魚（Triaenopogon barbatus）、小頭櫛孔蝦虎魚（Ctenotrypauchen microcephalus）7種7個蝦虎魚[11]。樣本用95%的乙醇混合保存后取魚鰭。

1.2 DNA提取

采用酚-氯仿法提取DNA，取魚尾鰭0.5 g放入1.5 mL的離心管中，將魚鰭剪碎后加入裂解液和蛋白酶渦旋混勻，離心后于56 ℃消化3 h，酚-氯仿抽提后2倍無水乙醇沉淀，75%乙醇洗滌吹干后以ddH2O溶解，4 ℃冰箱保存。

1.3 引物設計及PCR擴增

在上海生工合成Ward[12]文獻中魚類COI通用引物：FishF1（5'-TCAACCAACCACAAAGACATTGG

CAC-3'），? FishR1（5'-TAGACTTCTGGGTGGCCAAAG

AATCA-3'）。

組建一個25 μL PCR反應體系：

模板DNA? ? ? ?2.0 μ ddHO? ? ? ? ? ? 6.5 μlPremix Taq? ? ? 12.5 μL引物? ? ? ? ? ? ? ? 各2.0 μL

PCR擴增的程序為：

95 ℃? ? ? ? ?2 min

95 ℃變性40 s92 ℃退火1 min72 ℃延伸1 min? ? ? ? ? ?35個循環

72 ℃延伸5 min

PCR擴增產物通過瓊脂糖凝膠（1.5%）電泳15～20 min，凝膠電泳成像儀檢測后送至上海歐易生物技術有限公司進行雙向測序，測序引物為擴增引物。

1.4 數據分析

從BOLD數據庫中獲取鱸形目的石首魚科黃姑魚（Nibea albiflora）、馬鲅科的四指馬鲅（Eleutheronema tetradactylum）、鰈形目的舌鰨科短吻紅舌鰨（Cynoglossus joyneri）、半滑舌鰨（Cynoglossus semilaevis）共4種魚類的同源序列作為外群，同時從NCBI和BOLD數據庫下載天津可能存在的蝦虎魚序列，結合測序序列進行系統進化分析（表1）。對所獲得的序列利用MEGA進行多序列比對，人工核查后，用“Statistics”程序計算序列堿基組成、遺傳顛換值、變異位點、密碼子偏向性等。不同個體間的遺傳距離用“Distance”程序下的Kimura雙參數模型計算。用Modeltest檢測出最適合的核苷酸序列分析模型后采用鄰位連接法（Neighbor-Joining，NJ）構建進化樹Bootstrap置信值重復抽樣1 000次。

2 結果與分析

2.1 渤海灣7種蝦虎魚的COI基因序列特征

將擴增測序得到的COI序列進行blast比對，相似度大于98%的即為同一種，比對結果支持形態鑒定結果。多序列比對得到652 bp共有序列，證實COI可以用作物種鑒定條形碼[8]。

堿基特征分析得到7種蝦虎魚的COI基因部分序列符合硬骨魚類G+C含量（44.6%）小于A+T含量（55.4%）的特點[14]。密碼子的堿基使用頻率存在一定的偏向性，平均GC含量最高為第2位堿基（55.50%），4種堿基在第1、2、3位密碼子分布含量都不平均。本研究中COI基因密碼子變異情況表明，共有不變位點（522個）多于變異位點（126個），轉換位點（71個）多于顛換位點（55個），轉換與顛換比為1.3，其中第1位堿基變異位點最多，第3位最少。

2.2 蝦虎魚科12種魚類的種間遺傳距離

12種蝦虎魚的種間遺傳距離平均值為0.21。矛尾刺蝦虎魚和斑尾復蝦虎魚的種間遺傳距離最小（0.014），小于0.02的臨界值[13]，應為同種蝦虎魚，斑尾復蝦虎魚和紋縞蝦虎魚種間遺傳距離最大（0.435）。

2.3 蝦虎魚科魚類的分子系統進化樹

系統樹顯示，7種蝦虎魚中，紅狼牙蝦虎魚與矛尾刺蝦虎魚距離較遠，小頭櫛孔蝦虎魚和彈涂魚聚為一支，矛尾刺蝦虎魚和斑尾復蝦虎魚聚為一支，紋縞蝦虎魚和暗縞蝦虎魚遺傳距離較近。

3 討論與結論

3.1 7種渤海灣蝦虎魚類的COI基因序列特點

COI基因部分序列第2位堿基的GC含量最高，與石首魚科COI基因堿基特征一致[5]，與于亞男長江口蝦虎魚科魚類COI基因堿基特征不一致等[15]。第3位堿基的GC含量比較恒定，而第1位堿基則因種類不同而呈現含量差異，第1位堿基能否作為不同種類區分的特征也需要進一步擴大樣本量做深入的研究。本研究中7種蝦虎魚的密碼子堿基使用存在一定偏向性，但與李明暉等[16]研究結果不一致。由于其與基因長度、堿基組分、翻譯起始效應、基因的表達水平等多種因素有關[16-17]，需要系統地對蝦虎魚亞目不同地理群體的該基因密碼子堿基使用情況進行分析再下結論，或者精確到不同科屬的具體情況來進行探討。7種蝦虎魚的COI基因顛換大于轉換，轉換/顛換值（R）值為1.3小于2，且第1位堿基承擔著主要選擇壓力，表明蝦虎魚可能通過高度基因突變實現快速進化，從而增強其抗逆性和適應性[9，18]。

3.2 COI基因在渤海灣天津海域蝦虎魚科魚類分類中的應用

蝦虎魚的COI基因的種間和種內遺傳距離相差20倍以上，可作為天津海域蝦虎魚科魚類分子分類的條形碼和仔稚魚分子鑒定的依據[13]。

在系統進化樹上，12種蝦虎魚節點支持率較高，除了刺蝦虎魚屬、晴尾蝌蚪蝦虎魚與其他蝦虎魚遺傳距離較遠，其他蝦虎魚均聚為一大類。于亞男等[14]研究成果也出現了有的蝦虎魚先與外群聚在一起，之后再與其他蝦虎魚相聚的現象。該基因條形碼可以很好地區分不同種，但對于不同科或更大的分類界元的效果不太明顯。由于有的蝦虎魚比較稀有、難以采集，本研究采集樣品有限，需要加大樣本量進一步優化進化樹。形態學上，孔蝦虎魚和拉氏狼牙蝦虎魚具有相似的體型、體色、眼睛高度退化等特點，有些科學家將它們歸為近盲蝦虎魚亞科，這與進化樹分類結果比較一致。但是它們遺傳距離大于0.2，確定是不同的種。本研究中，背眼蝦虎魚科的彈涂魚與近盲蝦虎魚亞科的小頭櫛孔蝦虎魚聚為一支，相關研究也證實他們親緣關系較近，具體還需加大樣本進行確認。

本研究得出，矛尾刺蝦虎魚和斑尾復蝦虎魚為同一物種，種間遺傳距離僅為0.014，小于 Hebert等[13]提出的臨界值0.020；進化樹上，節點支持率達到100%聚為一支。

縞蝦虎魚是蝦虎魚科的一個屬，渤海常見的有鐘馗蝦虎魚（髭縞蝦虎魚）、紋縞蝦虎魚、裸頸縞蝦虎魚等。本研究中所有的縞蝦虎魚遺傳距離較近，與已知分類結果一致。晴尾蝌蚪蝦虎魚由于頭部形態特征與其他蝦虎魚相差較大，從形態上很好區分，與其他蝦虎魚科屬遺傳距離較遠。

由于蝦虎魚種類眾多，變異較大，其鑒定除了分子鑒定還需結合形態學和生態地理分布等因素，也可以結合其他條形碼基因綜合分析其進化、地理分布及系統分類等。

參考文獻：

[1] 李寶華， 張素青， 王秀芹， 等. 2004-2009年渤海灣（天津區域）魚蝦貝產卵場水質環境綜合評價[J]. 水產科學， 2010， 29（11）： 661-665.

[2] 徐海龍， 谷德賢， 喬秀亭， 等. 渤海灣主要漁業資源長度與體質量關系分析[J]. 南方水產科學， 2014， 10（1）： 57-63.

[3] 劉靜， 陳詠霞， 馬琳. 黃渤海魚類圖志[M]. 北京： 科學出版社， 2015.

[4] SACCONE C， DE GIORGI C， GISSI C， et al. Evolutionary genomics in Metazoa： the mitochondrial DNA as a model system[J]. Gene， 1999， 238（1）： 195-209.

[5] 柳淑芳， 陳亮亮， 戴芳群， 等. 基于線粒體CO1基因的DNA條形碼在石首魚科（Sciaenidae）魚類系統分類中的應用[J]. 海洋與湖沼， 2010， 41（2）： 223-232.

[6] 高天翔， 畢瀟瀟， 趙林林， 等. 基于線粒體Cytb基因全序列的松江鱸群體遺傳結構分析[J]. 水生生物學報， 2013， 37（2）： 199-207.

[7] 彭居俐， 王緒禎， 王丁， 等. 基于線粒體CO1基因序列的DNA條形碼在鯉科鲌屬魚類物種鑒定中的應用[J]. 水生生物學報， 2009， 33（2）： 271-276.

[8] 彭士明， 施兆鴻， 侯俊利， 等. 銀鯧3個野生群體線粒體COI基因的序列差異分析[J]. 上海海洋大學學報， 2009， 18（4）： 398-402.

[9] 畢瀟瀟， 高天翔， 肖永雙， 等. 4種鱈魚線粒體16S rRNA、COI和Cyt b基因片段序列的比較研究[J]. 南方水產， 2009， 5（3）： 46-52.

[10] 楊紀明. 渤海魚類的食性和營養級研究[J]. 現代漁業信息， 2001， 16（10）： 10-19.

[11] 李明德， 張洪杰， 等. 渤海魚類生物學[M]. 北京： 中國科學技術出版社， 1991： 62-68.

[12] WARD R D， ZEMLAK T S， INNES B H， et al. DNA barcoding Australia’s fish species[J]. Philosophical Transactions of the Royal Society B： Biological Sciences， 2005， 360（1462）： 1847-1857.

[13] HEBERT P D N， STOECKLE M Y， ZEMLAK T S， et al. Identification of birds through DNA barcodes[J]. PLoS Biology， 2004， 2（10）： e312.

[14] AQUINO L M G， TANGO J M， CANOY R J C， et al. DNA barcoding of fishes of Laguna de Bay， Philippines[J]. Mitochondrial DNA， 2011， 22（4）： 143-153.

[15] 于亞男， 宋超， 侯俊利， 等. 基于線粒體COI基因部分序列的長江口蝦虎魚科魚類系統分類[J]. 淡水漁業， 2014， 44（5）： 3-8.

[16] 李明暉， 安長廷， 李昂， 等. 黃渤海蝦虎魚DNA條形碼分類體系構建[J]. 中國水產科學， 2022， 29（8）： 1179-1188.

[17] FEDOROV A， SAXONOV S， GILBERT W. Regularities of context-dependent codon bias in eukaryotic genes[J]. Nucleic Acids Research， 2002， 30（5）： 1192-1197.

[18] 郭新紅， 劉少軍， 劉巧， 等. 魚類線粒體DNA研究新進展[J]. 遺傳學報， 2004， 31（9）： 983-1000.