中国水产科学  2023, Vol. 30 Issue (09): 1112-1126  DOI: 10.12264/JFSC2023-0139
0

引用本文 

刘乙蒙, 刘洋, 刘必林, 俞晔伟, 宋成辉, 王丛丛. Cyt b和12S rRNA基因条形码在灯笼鱼科鱼类物种鉴定中的应用[J]. 中国水产科学, 2023, 30(9): 1112-1126. DOI: 10.12264/JFSC2023-0139.
LIU Yimeng, LIU Yang, LIU Bilin, YU Yewei, SONG Chenghui, WANG Congcong. Identification analysis of DNA barcoding based on mitochondrial Cyt b and 12S rRNA genes in Myctophidae fishes[J]. Journal of Fishery Sciences of China, 2023, 30(9): 1112-1126. DOI: 10.12264/JFSC2023-0139.

基金项目

农业农村部全球渔业资源调查监测评估(公海渔业资源综合科学调查)专项;国家远洋渔业工程技术研究中心、大洋渔业资源可持续开发教育部重点实验室、农业农村部大洋渔业可持续利用重点实验室开放基金项目(A1-2006-23-200210).

作者简介

刘乙蒙(1999‒),女,硕士研究生,研究方向为海洋生物学. E-mail:1073217790@qq.com

通信作者

王丛丛,讲师,研究方向为海洋生物学. E-mail:ccwang@shou.edu.cn

文章历史

收稿日期:2023-06-13
修改日期:2023-08-15
Cyt b和12S rRNA基因条形码在灯笼鱼科鱼类物种鉴定中的应用
刘乙蒙1,刘洋1,2,3,4,5,刘必林1,2,3,4,5,俞晔伟1,宋成辉1,王丛丛1,2,3,4,5,     
1. 上海海洋大学海洋科学学院,上海 201306
2. 大洋渔业资源可持续开发教育部重点实验室,上海 201306
3. 国家远洋渔业工程技术研究中心,上海 201306
4. 农村农业部大洋渔业开发重点实验室,上海 201306
5. 农村农业部大洋渔业资源环境科学观测实验站,上海 201306
摘要:灯笼鱼科鱼类种类繁多,且同属鱼类形态学相近,因此利用分子标记对灯笼鱼进行准确的物种鉴定具有重要价值。为探讨线粒体细胞色素b基因(Cyt b)和12S rRNA基因在灯笼鱼科物种鉴定中的适用性,对西北太平洋采集的56尾灯笼鱼进行扩增,并进行序列对比与系统发育分析。研究表明,采集的样本包括6种灯笼鱼,分别为瓦氏角灯鱼(Ceratoscopelus warmingii)、长体标灯鱼(Symbolophorus californiensis)、粗鳞灯笼鱼(Myctophum asperum)、细泰勒灯鱼(Tarletonbeania crenularis)、日本背灯鱼(Notoscopelus japonicus)以及某背灯鱼属鱼类(Notoscopelus sp.)。核苷酸多态性分析显示,基于Cyt b基因的种内与种间遗传距离比基于12S rRNA基因的更大。比较灯笼鱼科2种基因序列的结构特征,发现Cyt b基因的种间平均遗传距离是种内平均遗传距离的25倍,12S rRNA基因的种间平均遗传距离是种内平均遗传距离的26倍,均符合作为DNA条形码的基本要求。系统进化分析显示,每种灯笼鱼均能形成独立分支,2个基因均能对6种灯笼鱼类进行鉴别;但在Cyt b基因构建的进化树中,每种鱼类能更好与数据库中已有的序列进行聚类。综上所述,Cyt b和12S rRNA作为DNA条形码可以有效地对灯笼鱼科鱼类物种进行鉴定,且Cyt b基因在系统进化关系的研究上具有更高的适用性。
关键词灯笼鱼科    Cyt b基因    12S rRNA基因    物种鉴定    系统进化分析    
Identification analysis of DNA barcoding based on mitochondrial Cyt b and 12S rRNA genes in Myctophidae fishes
LIU Yimeng1,LIU Yang1,2,3,4,5,LIU Bilin1,2,3,4,5,YU Yewei1,SONG Chenghui1,WANG Congcong,1,2,3,4,5    
1. College of Marine Sciences, Shanghai Ocean University, Shanghai 201306, China
2. Key Laboratory of Sustainable Exploitation of Oceanic Fisheries Resources, Ministry of Education, Shanghai 201306, China
3. National Distant-water Fisheries Engineering Research Center, Shanghai 201306, China
4. Key Laboratory of Ocean Fisheries Exploitation, Ministry of Agriculture and Rural Affairs, Shanghai 201306, China
5. Scientific Observing and Experimental Station of Oceanic Fishery Resources, Ministry of Agriculture and Rural Affairs, Shanghai 201306, China
Abstract:The Myctophidae family is widely distributed in mesopelagic oceans worldwide and plays a significant role in the pelagic food chain due to its vertical migratory behavior. Myctophidae are diverse and the morphological differences between genera are not obvious. Therefore, the use of molecular markers for Myctophidae identification would be advantageous. In this study, we amplified 56 Myctophidae samples collected from the Northwest Pacific using the Cyt b and 12S rRNA genes. These amplified sequences were subsequently used for phylogenetic analyses alongside sequences from 29 other Myctophidae species in 18 genera from the GenBank database. We explored the applicability of Cyt b and 12S rRNA genes for identifying Myctophidae species. The results indicate that the 56 samples belonged to 6 species: Ceratoscopelus warmingii, Symbolophorus californiensis, Myctophum asperum, Tarletonbeania crenularis, Notoscopelus japonicus, and Notoscopelus sp. The intra-specific and inter-specific genetic distances based on the Cyt b gene were more pronounced than those based on the 12S rRNA gene. The difference between the average inter-specific and intra-specific genetic distances of Cyt b and 12S rRNA gene was 25 and 26, respectively. These results indicated that both genes meet the basic requirements for DNA barcodes. Phylogenetic analysis evidenced that the six Myctophidae species each formed an independent clade based on Cyt b and 12S rRNA gene sequences. However, in the evolutionary tree constructed based on the Cyt b gene, each fish clustered more closely with sequences in the GenBank database. In summary, both Cyt b and 12S rRNA genes could be used as DNA barcodes for identifying Myctophidae species, amd the Cyt b gene was particularly suitable for analyzing phylogenetic relationships.
Key wordsMyctophidae    cytochrome b gene    12S rRNA gene    species identification    phylogenetic analysis    

灯笼鱼科(Myctophidae)鱼类包含34属225种。在世界各大洋广泛分布[1],主要分布于大洋的中深层[2],具有垂直迁移的生活特性[3]。白天灯笼鱼栖息于水下100~200 m处,夜晚浮出水面捕食桡足类等小型无脊椎动物[4-6]。同时,灯笼鱼也是大洋高经济价值鱼类的重要食物来源[7-8],位于大洋食物网的中生态位,起着重要的生态作用。灯笼鱼作为传统商业捕捞中常见的兼捕种[9],由于个体较小,经济价值较低,全球范围体系化的捕捞活动较少,商业化进展缓慢。随着高价值鱼类的过度捕捞,研究者将焦点放到以往被忽视的低价值鱼类和水生生物中,如灯笼鱼与深海头足类[2]。西北太平洋海域属于高生产力海洋系统,也是亲潮寒流与黑潮暖流交汇区域[10],虽是多种经济鱼类的渔场,但生物多样性水平较低。灯笼鱼科鱼类作为西北太平洋海域优势种[9],由于资源量大,蛋白质含量高,被人们逐渐重视。为了更好地实现灯笼鱼科鱼类资源的可持续利用和开发,应对灯笼鱼科鱼类进行正确的物种鉴定及更深入的系统进化研究。

传统的灯笼鱼科鱼类物种鉴定方法主要为形态学鉴定,依赖于灯笼鱼发光器的数目与排列,以及鳍条数目进行辨别[11-12]。由于灯笼鱼科鱼类种类繁多,且同属鱼类形态学相近,在捕捞过程中作为鉴定依据的发光器容易损坏与脱落,形态学鉴定易出现误差[1]。而分子生物学鉴定不依赖于样本的完整性,可通过DNA对样本进行精确鉴定[13]。DNA条形码是生物体内能够代表该物种的,标准的、具有一定变异的、易扩增的DNA片段[14],已成功在海洋鱼类物种鉴定与系统进化分析中运用,弥补了传统形态学鉴定的不足[15-18]。Batta-Lona等[3]利用CO I基因对形态学难以准确识别的灯笼鱼稚鱼进行鉴别,结果表明形态学与分子鉴定结合所得到的鉴定结果更加精确。Lee等[1]利用16S rRNA基因对灯笼鱼进行物种鉴定,发现了韩国海域以往从未出现的喀氏眶灯鱼(Diaphus garmani),结果表明16S rRNA基因可有效分辨灯笼鱼科鱼类隐存种。Gallet等[4]利用CO I基因对大西洋石油泄漏区的灯笼鱼遗传多样性进行研究,发现CO I基因可有效进行灯笼鱼科鱼类遗传多样性研究。可见以往使用DNA条形码对灯笼鱼进行物种鉴定主要以CO I和16S rRNA基因为主,鲜少将Cyt b和12S rRNA基因用于灯笼鱼科鱼类的物种鉴定中。

细胞色素b基因(cytochrome b, Cyt b)是线粒体基因的13个蛋白质编码基因之一,且进化速率适中,相对保守,是合适的物种鉴定标记[19-20]。杨天燕等[19]利用Cyt b基因对龙头鱼科及其近缘的灯笼鱼科等物种进行系统发育分析,结果表明灯笼鱼属(Myctophum)与底灯鱼属(Benthosema)间具有较近亲缘关系。王利华等[21]利用CO I与Cyt b基因对6种鲌属(Culter)鱼类的遗传距离与系统进化关系进行分析,探讨了这两种基因在鲌属鱼类分子鉴定中的可行性,发现Cyt b基因可以对6种鲌属鱼类进行更有效的区分。李楠等[22]利用形态学鉴定与分子进化树分析对CO I基因与Cyt b基因在中国常见河豚物种鉴定中的适用性进行研究,研究表明CO I基因与Cyt b基因均可对河豚鱼进行有效区分。12S rRNA属于不编码蛋白质,大部分变异为中性突变,进化速率较快,因此能够有效地对物种间的同源性进行比较,已被广泛应用于鱼类系统进化分析中[23-25]。张凌[26]利用12S rRNA基因构建系统进化树,研究湖鲚(Coilia macrognathos)、短颌鲚(Coilia brachygnathus)是否属于刀鲚(Coilia macrognathos)新亚种,研究表明湖鲚、短颌鲚实为刀鲚的淡水生态型种群,并非新的有效物种。於俊琦等[27]利用12S rRNA基因对浙南海域舌鳎属(Cynoglossidae)鱼类系统进化关系进行研究,结果表明,12S rRNA基因可对除长吻红舌鳎外(Cynoglossus lighti)其他舌鳎属鱼类进行有效系统进化分析。

本研究基于Cyt b基因和12S rRNA基因对西北太平洋采集的56尾灯笼鱼科鱼类进行物种鉴定以及系统发育分析,同时对这2种分子标记在灯笼鱼科鱼类物种鉴定的适用性进行探讨。

1 材料与方法 1.1 样品采集

本实验灯笼鱼样本由上海海洋大学“淞航”号科考船于2022年6—7月采集于西北太平洋海域,随后立即保存于−80 ℃冰箱运回实验室。随机从5个灯笼鱼种类丰富的站点选取56尾样本,详细采样信息见表1

表1  本研究采样站点信息 Tab. 1  Sampling sites information in the study
1.2 DNA提取及PCR扩增与测序

取灯笼鱼样本尾部肌肉组织20 mg,采用生工Ezup柱式动物基因组DNA抽提试剂盒提取样本基因组DNA。使用1%琼脂糖凝胶电泳检测后,存储于−20 ℃冰箱备用。

选取Cyt b通用引物H14724、H15149[28]以及12S rRNA通用引物M13Li[29]扩增线粒体DNA Cyt b、12S rRNA基因,引物序列如表2。PCR总反应体系为25 μL,包括12.5 μL 2X SanTaq PCR Mix,正反向引物各1 μL, DNA模板1 μL,超纯水9.5 μL。PCR反应程序为:95 ℃预变性5 min; 94 ℃变性30 s, 57 ℃退火45 s, 72 ℃延伸50 s,共35个循环;最后72 ℃延伸10 min,产物经1.8%琼脂糖凝胶电泳检测后,送往上海生工生物公司进行测序。

1.3 数据处理

使用NCBI Blast工具对所得样本序列进行物种鉴定。由于NCBI数据库中缺乏灯笼鱼科Cyt b基因与12S rRNA序列,选择下载灯笼鱼科19属29种鱼类的线粒体全基因组序列,选取翘嘴鲌(Culter alburnus)作为外群,详见表3,并于线粒体全基因组序列截取相应的Cyt b基因与12S rRNA基因序列。使用MEGA 7.0软件进行剪切比对,统计基因序列以及相应核苷酸序列变异情况。使用K2P模型对灯笼鱼科鱼类基于Cyt b基因及12S rRNA基因的种内、种间以及属内、属间遗传距离进行计算,并使用邻接法(neighour- joining, NJ)构建分子系统进化树。

表2  引物信息 Tab. 2  PCR primer information
表3  用于构建系统进化树的物种序列来源 Tab. 3  Sequences resources used in phylogenetic tree construction
2 结果与分析 2.1 灯笼鱼基因组DNA提取和PCR扩增

从灯笼鱼样本中提取片段完整的基因组DNA。利用Cyt b和12S rRNA基因通用引物进行PCR扩增,获得特异性条带,产物电泳结果如图1所示。

图1  灯笼鱼科鱼类Cyt b基因及12S rRNA基因PCR产物 Fig. 1  PCR profile of Cyt b and 12S rRNA genes of Myctophidae
2.2 物种鉴定

56尾灯笼鱼Cyt b和12S rRNA的序列比对和物种鉴定结果见表4。除T33-1、T5-1及T14-8这3尾灯笼鱼应用12S rRNA鉴定为日本背灯鱼(Notoscopelus japonicus), Cyt b鉴定为长体背灯鱼(Notoscopelus elongatus)外,其余样本两基因鉴定结果一致。此前经形态学初步分类为粗鳞灯笼鱼(Myctophum asperum)的灯笼鱼中,T33站点的在本研究中鉴定为粗鳞灯笼鱼;T14站点的在本研究中鉴定为细泰勒灯鱼(Tarletonbeania crenularis); T15站点的形态如图2a的灯笼鱼被鉴定为粗鳞灯笼鱼,形态如图2b的被鉴定为细泰勒灯鱼。经形态学鉴定初步分类为白皮灯笼鱼的灯笼鱼样本中,在T33、T15、T14 3个站点中均被鉴定为瓦氏角灯鱼(Ceratoscopelus warmingii),而在站点T5中形态如图2c的被鉴定为瓦氏角灯鱼,形态如图2d的被鉴定为长体标灯鱼(Symbolophorus californiensis)。黑眼灯笼鱼经分子生物学鉴定为日本背灯鱼(Notoscopelus japonicus)。

2.3 灯笼鱼Cyt b和12S rRNA序列特征分析

56尾灯笼鱼序列特征如表5所示。Cyt b基因中包含81个保守位点,326个变异位点,171个简约信息位点及155个单突变位点;12S rRNA基因中包含145个保守位点,629个变异位点,588个简约信息位点及40个单突变位点。Cyt b基因的颠换值为1.55, 12S rRNA基因的颠换值为1.01。

Cyt b基因的平均碱基组成为C 31.2%、T 27.2%、A 24.9%、G 16.7%, 12S rRNA基因的平均碱基组成为T 26.3%、G 26.1%、C 24.7%、A 22.9%, Cyt b基因A+T (52.1%)含量高于C+G (47.9%)含量,表现出A+T的偏倚性;12S rRNA序列A+T (49.2%)含量稍小于C+G (50.8%)含量。且Cyt b序列中第三密码子呈明显反CG偏倚,CG含量为42.5%。而12S rRNA序列中反CG偏倚则体现在第一密码子中,为47.1%。

表4  2种鱼类DNA条形码序列在56尾灯笼鱼中的物种鉴定结果 Tab. 4  Species identification of two candidate barcode sequences in 56 fishes of Myctophidae
图2  灯笼鱼部分样本图片 Fig. 2  Partial sample photos of fishes in Myctophidae
表5  碱基平均分布频率 Tab. 5  Average nucleotide frequencies

比较分析基于Cyt b和12S rRNA 2种基因序列的6种灯笼鱼种内种间核苷酸多态性,以100 bp核苷酸开口中点为横坐标,并分别以核苷酸多态性指数和核苷酸多态位点个数为纵坐标,绘制Cyt b基因与12S rRNA基因的核苷酸多态性指数分析图(图3)及核苷酸多态位点数图(图4)。

图3a显示灯笼鱼科鱼类总体Cyt b基因核苷酸多态性较高,鉴定得出的6种灯笼鱼科鱼类种内核苷酸多态性均在较低水准,且高变区域多在300~400 bp处。由图3b可知5种12S rRNA基因鉴定得出的灯笼鱼高变区域为500~650 bp,尤其是日本背灯鱼的核苷酸多态性明显高于其他物种。

图4可知12S rRNA基因相较于Cyt b基因种内突变更高,且日本背灯鱼12S rRNA基因中430~650 bp间出现两个高变区,粗鳞灯笼鱼高变区则出现在250~300 bp间。而Cyt b基因则表现出较低的种内突变,以及较高的种间突变。

2.4 灯笼鱼种内种间遗传距离和barcoding gap检验

基于Cyt b与12S rRNA基因的灯笼鱼科鱼类遗传距离分析结果见表6。在12S rRNA基因中,笔者剔除了T33-1、T5-1、T14-8 3个鉴定为背灯鱼属的样本后进行遗传距离分析,得到种内遗传距离为0.003~0.012 (平均为0.006),种间遗传距离为0.108~0.184 (平均为0.155),种间平均遗传距离是种内遗传距离的26倍。在Cyt b基因中,种内遗传距离为0.001~0.029 (平均为0.009),种间的平均遗传距离为0.224,是种内遗传距离的25倍。

图3  基于Cyt b (a)和12S rRNA (b)的6种灯笼鱼种内种间核苷酸多态性比较CW:瓦氏角灯鱼;MA:粗鳞灯笼鱼;NJ:日本背灯鱼;SC:长体标灯鱼;TC:细泰勒灯鱼;NE:长体背灯鱼. Fig. 3  Comparison of intra-specific and inter-specific nucleotide polymorphism based on Cyt b (a) and 12S rRNA (b) sequences of 6 Myctophidae speciesCW: Ceratoscopelus warmingii; MA: Myctophum asperum; NJ: Notoscopelus japonicus; SC: Symbolophorus californiensis; TC: Tarletonbeania crenularis; NE: Notoscopelus elongatus.
图4  基于Cyt b (a)、12S rRNA (b)的6种灯笼鱼种内种间核苷酸多态位点比较CW:瓦氏角灯鱼;MA:粗鳞灯笼鱼;NJ:日本背灯鱼;SC:长体标灯鱼;TC:细泰勒灯鱼;NE:长体背灯鱼. Fig. 4  Comparison of intra-specific and inter-specific of nucleotide variable sites based on Cyt b (a) and 12S rRNA (b) sequences of 6 Myctophidae speciesCW: Ceratoscopelus warmingii; MA: Myctophum asperum; NJ: Notoscopelus japonicus; SC: Symbolophorus californiensis; TC: Tarletonbeania crenularis; NE: Notoscopelus elongatus.
表6  灯笼鱼科Cyt b及12S rRNA序列不同分类水平遗传距离 Tab. 6  Genetic distances according to different taxonomic levels within Myctophidae based on Cyt b and 12S rRNA sequences

Barcoding gap是理想遗传变异条件下,种内与种间遗传距离之间明显的不重叠间隔区(图5a)[30-32]。由图5可得Cyt b基因与12S rRNA基因均存在明显的条形码间隙,不存在重叠区,以证实各物种的有效性。

2.5 灯笼鱼系统进化关系

为进一步验证分子鉴定结果的准确性,采用邻接法构建基于12S rRNA (图6)和Cyt b基因(图7)序列的灯笼鱼科鱼类进化树。基于Cyt b和12S rRNA基因序列的NJ进化树自测序列结果基本一致,同种个体均先聚集成为一个单系,且具有较高的支持率。基于12S rRNA基因绘制的NJ进化树中,本研究自测序列与NCBI数据库中同种灯笼鱼序列并未聚为一支,而是分为两个大支。

图5  不同DNA条形码候选序列遗传距离分析图条形码间隙演示图(a), 12S rRNA(b), Cyt b(c)序列的K2P遗传距离分布. Fig. 5  Gentic distance distribution of different DNA candidate bar codesSchematic of the barcoding gap(a), K2P genetic distance of 12S rRNA (b) and Cyt b (c) sequences.
图6  基于Cyt b基因的灯笼鱼科鱼类NJ进化树 Fig. 6  NJ tree based on Cyt b sequence of Myctophidae
图7  基于12S rRNA基因的灯笼鱼科鱼类NJ进化树 Fig. 7  NJ tree based on 12S rRNA sequence of Myctophidae
3 讨论

灯笼鱼因种类繁多、仔稚鱼形态特征相近,对传统形态学分类造成巨大影响。DNA条形码可简便高效地为物种分类鉴定提供有力证据。BLAST是用于比对研究样本与数据库序列相似度以确认物种的搜索程序,置信度越高代表比对结果越可信[33]。Zhang等[34]将置信度98%作为判断所测样本与比对数据库样本是否属于同一种类的分界值。本研究中T33-1、T5-1、T14-8这3个样本BLAST比对结果分别为日本背灯鱼(Cyt b)与长体背灯鱼(12S rRNA),但其置信度均低于98%,由于NCBI数据库中灯笼鱼属Cyt b与12S rRNA基因数据缺乏,推断这3个样本可能为背灯鱼属其他鱼类物种。整体上12S rRNA基因(92.68%~99.86%)置信度高于Cyt b基因(87.30%~ 99.76%),就分子鉴定结果而言12S rRNA基因在灯笼鱼科鱼类物种鉴定中置信度更高。

碱基偏倚性在生物中广泛存在,不同类型不同生境的物种,其碱基使用也存在较大差异,对碱基偏倚性的分析有助于了解物种进化过程中自然选择与基因突变的过程[19]。李明辉等[35]利用DNA条形码技术对我国黄渤海虾虎鱼(Gobilidae)进行物种鉴定与系统进化研究,并分析了13种黄渤海虾虎鱼的碱基组成,发现A+T碱基含量偏高。本研究中,Cyt b基因A+T (52.1%)>C+G (47.9%)符合海水鱼类中A+T含量偏高规律,这与Ward等[30]研究的海水鱼类碱基偏倚结果一致。而12S rRNA基因中A+T (49.2%)<C+G (50.8%), C+G含量较高,呈现不明显的CG偏倚。通常DNA序列在进化过程中基因转换的频率会高于基因颠换的频率,因此转换/颠换的比值R可用于评估DNA进化速率[35]。武宝生等[36]研究认为,当R值小于2时,基因突变达到饱和状态,并且受到进化噪声影响的可能性较大。本研究中灯笼鱼Cyt b基因颠换值R为1.55, 12S rRNA基因颠换值为1.01,均小于2,说明这两个基因突变已达饱和状态,并且在进化过程中都受到了进化噪声影响[37-38],且12S rRNA进化速率更快。

核苷酸多样性与种内种间平均遗传距离是衡量线粒体DNA遗传变异的重要依据[38]。且核苷酸多样性指数可以兼顾不同的线粒体DNA单倍型在群体中所占的比例,因此相较于单纯的平均遗传距离,两者相结合可以更好地反映线粒体DNA多态性[39]。蒋芝华等[40]认为,高的种内与种间核苷酸多态性在鉴别近缘物种时结果更加精确。本研究中12S rRNA基因种间核苷酸多态性低于Cyt b基因,且12S rRNA基因灯笼鱼科种间的核苷酸多态性指数最低值(0.041)小于日本背灯鱼种内核苷酸多态性最高值(0.045), Cyt b基因种间核苷酸多态性指数最低值(0.115)远高于种内核苷酸多态性指数(0.018), Cyt b基因相较于12S rRNA基因具有更明显的种内与种间核苷酸多态性差距,且Cyt b种内核苷酸多态位点数也处于低水平状态,说明灯笼鱼科鱼类Cyt b序列中种内突变更少,而种间突变更为明显,因此Cyt b基因更适用于区分灯笼鱼科鱼类近缘物种。

在分子生物学鉴定中,物种种内与种间遗传距离是物种界定的标准[41]。本研究中灯笼鱼科鱼类12S rRNA基因种间遗传距离是种内遗传距离的26倍,而Cyt b基因的种间遗传距离是种内遗传距离的25倍,此外Hebert等[42]还指出,物种的种内遗传距离应小于DNA条形码规定的最大种内遗传距离0.020;在种内遗传距离中,两种基因所有物种最大种内遗传距离均小于0.020,因此均符合DNA条形码鉴定标准及“10×法则”[30,42],可对灯笼鱼科鱼类进行准确鉴定。DNA条形码种间变异与种内变异间通常不存在重叠区,并应有明显间隔,即DNA条形码间隙(DNA barcoding gap)[43]。本研究Cyt b与12S rRNA基因均存在条形码间隙,且Cyt b基因的条形码间隙明显大于12S rRNA基因。

线粒体DNA可以精确构建系统进化关系[21]。本研究利用获得的Cyt b和12S rRNA基因进行灯笼鱼科鱼类系统进化分析,发现12S rRNA构建的NJ进化树3尾背灯鱼属鱼类与其他日本背灯鱼聚在一起形成一个复合支,亲缘关系较近,其余样本各聚为一支。基于12S rRNA基因构建的NJ进化树,自测序列与数据库序列无法形成有效聚类,分为两大支,说明12S rRNA基因适合于灯笼鱼科鱼类物种鉴定,但在较高级分类阶元中不能有效地体现灯笼鱼科鱼类进化关系。Cyt b基因构建的NJ进化树中,本研究自测序列可与数据库序列形成有效聚类,并清晰展现出系统进化关系,更适用于灯笼鱼系统进化分析。

综上所述,12S rRNA基因与Cyt b基因均能有效的对灯笼鱼科鱼类进行物种鉴定。且Cyt b基因能更好与数据库中序列进行聚类,更适用于研究灯笼鱼科鱼类系统进化关系。由于数据库中灯笼鱼科鱼类Cyt b和12S rRNA基因序列条形码信息缺乏,所以在本实验中仍有部分灯笼鱼科鱼类无法精确鉴定。因此在未来的研究中应不断完善灯笼鱼科鱼类DNA条形码数据库,并结合多基因进行比对,从而获得更准确的鉴定结果。

参考文献
[1]
Lee H L, Yu H J, Kim J K. New descriptions of four larval lanternfish species and cryptic diversity of Benthosema pterotum (Pisces: Myctophidae) from the Northwest Pacific[J]. Journal of Fish Biology, 2022, 101(6): 1474-1500.
[2]
Li L N, Zhu W Z, Miao S C. Development and utilization of lantern fish resources[J]. Fishery Information & Strategy, 2012, 27(3): 240-245. [李励年,朱文钊,缪圣赐. 灯笼鱼资源的开发与利用[J]. 渔业信息与战略,2012, 27(3): 240- 245.]
[3]
Batta-Lona P G, Galindo-Sánchez C E, Arteaga M C, et al. DNA barcoding and morphological taxonomy: Identification of lanternfish (Myctophidae) larvae in the Gulf of Mexico[J]. Mitochondrial DNA Part A, DNA Mapping, Sequencing, and Analysis, 2019, 30(2): 375-383.
[4]
Gallet A, Koubbi P, Léger N, et al. Low-diversity bacterial microbiota in southern ocean representatives of lanternfish Genera Electrona, Protomyctophum and Gymnoscopelus (Family Myctophidae)[J]. PLoS ONE, 2019, 14(12): e0226159.
[5]
Cohen D L, Beckley L E. Diet and prey selectivity of the mesopelagic lanternfish Myctophum asperum from the Perth canyon, western Australia[J]. Ichthyological Research, 2021, 68(2): 294-302.
[6]
Habano A, Kobari T, Ohbayashi W, et al. Feeding habits of the skinnycheek lanternfish [Benthosema pterotum (Alcock, 1890)] in Kagoshima Bay, southern Japan[J]. Ichthyological Research, 2021, 68(1): 164-170.
[7]
Lin C H, Lin J S, Chen K S, et al. Feeding habits of bigeye tuna (Thunnus obesus) in the western Indian Ocean reveal a size-related shift in its fine-scale piscivorous diet[J]. Frontiers in Marine Science, 2020, 7: 582571.
[8]
miR-Arguimbau J, Navarro J, Balcells M, et al. Feeding ecology of blue whiting (Micromesistius poutassou) in the NW Mediterranean: The important role of Myctophidae[J]. Deep Sea Research Part I: Oceanographic Research Papers, 2020, 166: 103404.
[9]
Zhao G Q, Shi Y C, Fan W, et al. Study on main catch composition and fishing ground change of light purse seine in Northwest Pacific[J]. South China Fisheries Science, 2022, 18(1): 33-42. [赵国庆,石永闯,樊伟,等. 西北太平洋灯光围网主要渔获物组成及渔场变动研究[J]. 南方水产科学,2022, 18(1): 33-42.]
[10]
Han Y S, Jiang W, Xiao Y W, et al. Main change characteristics and influencing factors of tropical cyclones under the background of global change[J]. Advances in Earth Science, 2023, 38(5): 515-532. [韩岩松,姜伟,肖玉雯,等. 全球变化背景下热带气旋主要变化特征及影响因素[J]. 地球科学进展,2023, 38(5): 515-532.]
[11]
Cavallaro M, Guerrera M C, Abbate F, et al. Morphological, ultrastructural and immunohistochemical study on the skin ventral photophores of Diaphus holti Tåning, 1918 (Family: Myctophidae)[J]. Acta Zoologica, 2021, 102(4): 405-411.
[12]
Saunders R A, Lourenço S, Vieira R P, et al. Length-weight and otolith size to standard length relationships in 12 species of Southern Ocean Myctophidae: A tool for predator diet studies[J]. Journal of Applied Ichthyology, 2021, 37(1): 140- 144.
[13]
Behrens-Chapuis S, Herder F, Geiger M F. Adding DNA barcoding to stream monitoring protocols-What’s the additional value and congruence between morphological and molecular identification approaches?[J]. PLoS ONE, 2021, 16(1): e0244598.
[14]
Qu X Y, Liu L, Li C H, et al. Construction and application of DNA barcoding system for fishes in Yellow Sea and Bohai Sea[J]. Chinese Journal of Fisheries, 2022, 35(6): 37-44. [曲欣宇,刘璐,李纯厚,等. 黄、渤海鱼类DNA条形码信息平台构建及应用[J]. 水产学杂志,2022, 35(6): 37-44.]
[15]
Roy M, Belliveau V, Mandrak N E, et al. Development of environmental DNA (eDNA) methods for detecting high-risk freshwater fishes in live trade in Canada[J]. Biological Invasions, 2018, 20(2): 299-314.
[16]
Lira N L, Tonello S, Lui R L, et al. Identifying fish eggs and larvae: From classic methodologies to DNA metabarcoding[J]. Molecular Biology Reports, 2023, 50(2): 1713-1726.
[17]
Yang L, Li A, Li B S, et al. DNA barcode analysis: Species identification and phylogenetic relationships of the Carangidae family[J]. Journal of Fishery Sciences of China, 2021, 28(12): 1523-1535. [杨龙,李昂,李步苏,等. 鲹科鱼类DNA条形码分类系统[J]. 中国水产科学,2021, 28(12): 1523- 1535.]
[18]
Wang W X, Li Y, Zhu H G, et al. Genetic differentiation of shortjaw tapertail anchovy Coilia brachygnathus population based on Cyt b and morphology[J]. Fisheries Science, 2022, 41(6): 949-958. [王玮欣,李宇,朱洪赓,等. 基于形态和Cyt b序列的短颌鲚种群遗传分化研究[J]. 水产科学,2022, 41(6): 949-958.]
[19]
Yang T Y, Jiang Y L, Guo Y J, et al. Molecular phylogeny of Harpadonnehereus and its close relatives based on mitochondrial Cyt b gene[J]. Transactions of Oceanology and Limnology, 2020(6): 77-85. [杨天燕,蒋艳琳,郭易佳,等. 基于线粒体Cyt b基因的龙头鱼及其近缘物种分子系统发育研究[J]. 海洋湖沼通报,2020(6): 77-85.]
[20]
Zhu W L, Xiao S, Yang C L, et al. Genetic structure analysis of Acanthopagrus latus populations along South China coast based on mitochondrial Cyt b gene sequences[J]. Southwest China Journal of Agricultural Sciences, 2021, 34(2): 446-454. [朱威霖,肖珊,杨春玲,等. 基于线粒体Cyt b基因序列的华南沿海黄鳍棘鲷种群遗传结构分析[J]. 西南农业学报,2021, 34(2): 446-454.]
[21]
Wang L H, Luo X Z, Wang D, et al. Identification analysis of DNA barcoding based on mitochondrial CO I and Cyt b genes in Culter fishes[J]. Freshwater Fisheries, 2019, 49(4): 22-28. [王利华,罗相忠,王丹,等. 基于CO I和Cyt b DNA条形码在鲌属鱼类物种鉴定中的应用[J]. 淡水渔业,2019, 49(4): 22-28.]
[22]
Li N, Wang J H, Han C H, et al. Application of COI and Cyt b gene in species identification of pufferfish[J]. Chinese Journal of Food Hygiene, 2018, 30(1): 6-11. [李楠,王佳慧,韩春卉,等. COI及Cyt b基因对河鲀鱼鱼种鉴定的适用性研究[J]. 中国食品卫生杂志,2018, 30(1): 6-11.]
[23]
Jiang P W, Li M, Zhang S, et al. Construction of DNA meta- barcode database of fish in Pearl River Estuary based on mitochondrial cytochrome COI and 12S rDNA gene[J]. South China Fisheries Science, 2022, 18(3): 13-21. [蒋佩文,李敏,张帅,等. 基于线粒体COI和12S rDNA基因构建珠江河口鱼类DNA宏条形码数据库[J]. 南方水产科学,2022, 18(3): 13-21.]
[24]
Chen Z, Gao T X. Differences in the identification of marine fish by mitochondrial 12S and COI barcodes[J]. Journal of Hainan Tropical Ocean University, 2023, 30(2): 10-16. [陈治,高天翔. 线粒体12S与COI条形码对海洋鱼类的鉴定差异[J]. 海南热带海洋学院学报,2023, 30(2): 10-16.]
[25]
Andriyono S, Alam M J, Kim H W. The jawa and Bali Island marine fish molecular identification to improve 12S rRNA- tRNA valin-16S rRNA partial region sequences on the GenBank database[J]. Thalassas: an International Journal of Marine Sciences, 2020, 36(2): 343-356.
[26]
Zhang L. The species validity analysis of Coilia brachygnathus and Coilia nasus taibuensis based on themitochondrial 12srrna genes sequences data[J]. Journal of Chaohu College, 2012, 14(6): 55-58. [张凌. 基于线粒体12S rRNA基因序列探讨湖鲚、短颌鲚物种有效性[J]. 巢湖学院学报,2012, 14(6): 55-58.]
[27]
Yu J Q, Wu J G, Zhou Z M, et al. Phylogenetic analysis of six species of Cynoglossus alone the south coastal area in Zhejiang provinceon on partial sequence of mitochondrial 12S ribosomal RNA gene[J]. Bulletin of Science and Technology, 2012, 28(3): 43-48. [於俊琦,伍锦姑,周志明,等. 基于线粒体12S rRNA基因部分序列的6种浙南海域舌鳎属鱼类系统进化研究[J]. 科技通报,2012, 28(3): 43-48.]
[28]
Xiao W, Zhang Y, Liu H. Molecular systematics of Xenocyprinae (Teleostei: Cyprinidae): Taxonomy, biogeography, and coevolution of a special group restricted in East Asia[J]. Molecular Phylogenetics and Evolution, 2001, 18(2): 163-173.
[29]
Stoeckle M Y, Das Mishu M, Charlop-Powers Z. Improved environmental DNA reference library detects overlooked marine fishes in new jersey, United States[J]. Frontiers in Marine Science, 2020, 7: 226.
[30]
Ward R D, Zemlak T S, Innes B H, et al. DNA barcoding Australia’s fish species[J]. Philosophical Transactions of the Royal Society of London Series B: Biological Sciences, 2005, 360(1462): 1847-1857.
[31]
Li H, Cao S A, Li Y A, et al. Molecular assessment of demersal fish diversity in Prydz Bay using DNA taxonomy[J]. Deep Sea Research Part II: Topical Studies in Oceanography, 2022, 202: 105140.
[32]
Ferrari C, Tovela E, Taviani E, et al. DNA barcoding to assess species identification in museum samples of Amphiliidae and natural samples of Cichlidae from Southern Mozambique[J]. Rendiconti Lincei Scienze Fisiche e Naturali, 2022, 33(4): 713-720.
[33]
Ren J Y, Guo R, Jiao X Y. From the origins of SARS-CoV-2 to the correct application of BLAST[J]. Chinese Journal of Biochemistry and Molecular Biology, 2022, 38(2): 221-227. [任建英,郭睿,焦向英. 从新冠病毒起源到BLAST工具的正确实践[J]. 中国生物化学与分子生物学报,2022, 38(2): 221-227.]
[34]
Zhang J J, Deng Q T, Liu X D, et al. Molecular identification of commercial fish maws by DNA sequencing of 16S rRNA and cytochrome oxidase I genes[J]. Journal of Food Protection, 2022, 85(10): 1439-1445.
[35]
Li M H, An C T, Li A, et al. DNA barcodes enable higher taxonomic assignments in goby in the Bohai Sea and Yellow Sea of China[J]. Journal of Fishery Sciences of China, 2022, 29(8): 1179-1188. [李明晖,安长廷,李昂,等. 黄渤海虾虎鱼DNA条形码分类体系构建[J]. 中国水产科学,2022, 29(8): 1179-1188.]
[36]
Wu B S, Si L Z, Kong X Y, et al. The analysis of ITS1 characteristics of 11 species from 5 families and its application in phylogenetic research[J]. Journal of Fisheries of China, 2018, 42(4): 465-475. [武宝生,司李真,孔晓瑜,等. 5科11种鱼类ITS1特征分析及其在系统分类研究中的适用性[J]. 水产学报,2018, 42(4): 465-475.]
[37]
Yang Y H. The genetic diversity of the glyptosternum maculatum populations[D]. Wuhan: Wuhan Polytechnic University, 2019. [杨钰慧. 黑斑原鮡种群遗传多样性的研究[D]. 武汉:武汉轻工大学,2019.]
[38]
Rozas J, Sánchez-DelBarrio J C, Messeguer X, et al. DNA SP, DNA polymorphism analyses by the coalescent and other methods[J]. Bioinformatics, 2003, 19(18): 2496-2497.
[39]
Du Q Y, Dong F J, Chang Z J. Genetic diversity of the mtDNA 16S rRNA and 12S rRNA genes in the population of Misgurnus anguillicaudatus from North Henan[J]. Journal of Henan Agricultural Sciences, 2014, 43(11): 137-140, 146. [杜启艳,董方娟,常重杰. 豫北泥鳅线粒体DNA 16S rRNA和12S rRNA基因序列的遗传多样性分析[J]. 河南农业科学,2014, 43(11): 137-140, 146.]
[40]
Jiang Z H, Ding B B, Wang Y, et al. DNA barcoding application to Sciaenidae[J]. Journal of Nuclear Agricultural Sciences, 2018, 32(4): 673-680. [蒋芝华,丁棒棒,王熠,等. 石首鱼科海洋鱼类DNA条形码的构建[J]. 核农学报,2018, 32(4): 673-680.]
[41]
Liu H Y, Cai J, Xie Z G, et al. DNA barcodes for species identification and systematic evolution of Cobitidae fish[J]. Acta Agriculturae Universitatis Jiangxiensis, 2020, 42(4): 766- 777. [刘红艳,蔡金,谢仲桂,等. 鳅科鱼类DNA条形码鉴定及系统进化研究[J]. 江西农业大学学报,2020, 42(4): 766-777.]
[42]
Hebert P D N, Cywinska A, Ball S L, et al. Biological identifications through DNA barcodes[J]. Proceedings Biological Sciences, 2003, 270(1512): 313-321.
[43]
Zhao W H, Yi S K, Su J X, et al. Genetic diversity of Triplophysa yarkandensis populations in Tarim River Basin in Xinjiang[J]. Acta Hydrobiologica Sinica, 2022, 46(3): 364- 374. [赵文浩,易少奎,苏君晓,等. 新疆塔里木河叶尔羌高原鳅遗传多样性研究[J]. 水生生物学报,2022, 46(3): 364-374.]