鱼类环境DNA metabarcoding片段的近缘物种识别差异

陈治; 马春来; 叶乐; 杨超杰; 王海山

doi:10.12284/hyxb2022150

鱼类环境DNA metabarcoding片段的近缘物种识别差异

doi: 10.12284/hyxb2022150

陈治^{1, 2,},
马春来^{1, 2},
叶乐^{1, 2},
杨超杰^{1, 2},
王海山^{1, 2, ,}

1.
海南热带海洋学院热带海洋生物资源利用与保护教育部重点实验室，海南三亚 572022
2.
海南热带海洋学院海南省热带海洋渔业资源保护与利用重点实验室，海南三亚 572022

基金项目: 国家自然科学基金（32002389）；海南省自然科学基金(422RC717)；海南热带海洋学院引进人才科研启动资助项目（RHDRC201907）。

详细信息

作者简介:
陈治（1990-），男，山东省泰安市人，讲师，博士，主要从事海洋鱼类分类学及其资源保护研究。E-mail：change@139.com

通讯作者:
王海山（1983-），男，讲师，博士，主要从事海洋生物多样性及其系统地理学研究。E-mail：jazz-123@163.com

中图分类号: S931；Q523
计量
- 文章访问数: 468
- HTML全文浏览量: 364
- PDF下载量: 85
- 被引次数: 2
出版历程
- 收稿日期: 2022-02-17
- 修回日期: 2022-03-25
- 网络出版日期: 2022-04-24
- 刊出日期: 2022-08-15

Differences of eDNA metabarcoding fragments in relative fish species resolution

Chen Zhi^{1, 2
,},
Ma Chunlai^{1, 2},
Ye Le^{1, 2},
Yang Chaojie^{1, 2},
Wang Haishan^{1, 2
, ,}

1.
Key Laboratory of Utilization and Conservation for Tropical Marine Bioresources, Ministry of Education, Hainan Tropical Ocean University, Sanya 572022, China
2.
Hainan Key Laboratory for Conservation and Utilization of Tropical Marine Fishery Resources, Hainan Tropical Ocean University, Sanya 572022, China

摘要

摘要: 已知的鱼类环境DNA（eDNA）metabarcoding片段均未被针对性考察其对近缘物种的适用性，实际使用过程中存在“物种丢失”风险。为筛选出物种识别率最高的片段，本研究比较了15个主流片段对106属（共935种）鱼类的识别差异。研究结果如下：（1）蛋白质编码基因（COI，片段15）的物种识别率最高，但其引物通用性最差；片段09、片段11、片段07、片段03、片段12的引物序列总平均遗传距离也较大，均存在eDNA低效扩增的风险；（2）片段长度影响物种识别率，核糖体基因中片段05、片段06、片段01、片段02及片段13的物种识别率较高；（3）非度量多维尺度分析（NMDS）显示，不同基因、同一基因不同片段的识别结果存在较大差异，应考虑多片段、多基因组合应用；片段01与片段02、片段05与片段06等在NMDS图上距离较近，存在相互替代性；（4）物种类群影响识别结果，eDNA研究仍需要进一步开发高识别率片段。综合物种识别率、引物通用性、NMDS分析等多方面因素，本研究推荐2×150 bp测序平台使用片段01（MiFish-U）、2×250 bp测序平台使用片段05（Ac12S），辅以片段13（Vert-16S-eDNA）等进行近缘鱼类多样性调查。本研究旨在为提高鱼类eDNA调查结果准确性提供一定技术支撑。
- 环境DNA metabarcoding /
- 近缘鱼类识别 /
- 12S /
- 扩增长度 /
- 多片段
Abstract: The applicability of environmental DNA (eDNA) metabarcoding fragments to relative fish species had not been compared. There was a risk of “species loss” in diversity surveys. In order to screen out the best fragments, we compared the resolution rate differences of 15 eDNA metabarcoding fragments in 106 genera (a total of 935 species). The results were as follows: (1) the protein-coding gene (COI, fragment 15) had the highest resolution rate, but the universality of its corresponding primer pairs was the worst; the overall mean distance based on primer pair sequence of fragment 09, fragment 11, fragment 07, fragment 03 and fragment 12 were obviously large, suggesting their eDNA amplification efficiency were possibly low. (2) The resolution rates were significantly affected by the length of fragments, and the fragment 05, fragment 06, fragment 01, fragment 02 and fragment 13 of ribosomal genes had a higher resolution rate except fragment 15. (3) Non-metric multidimensional scaling analysis (NMDS) showed that there were great differences among different genes and different fragments belonging to the same gene. Therefore, the combination application of multi-fragment and multi-gene should be considered; besides, fragment 01 and fragment 02, and fragment 05 and fragment 06 were close to each other in the NMDS plot. They function of fish resolution were overlapped. (4) Species groups affected the resolution results, and eDNA studies stilled need to develop fragments with higher resolution rates. Based on the resolution rate of relative species, the universality of primer pairs and NMDS analysis, this study recommended fragment 01 (Mifish-U) for 2×150 bp sequencing platform and fragment 05 (Ac12S) for 2×250 bp sequencing platform, supplemented by fragment 13 (Vert-16S-eDNA) to investigate the diversity of relative fish. This study provided some support for improving the accuracy of fish eDNA survey results.
- eDNA metabarcoding /
- relative fish species resolution /
- 12S /
- amplicon length /
- multi-fragments

HTML全文

1. 引言

近年来，环境DNA（Environmental DNA，eDNA） metabarcoding技术开始广泛应用于渔业生态学领域，在物种多样性调查、外来鱼种监测、濒危鱼类保护方面展现了极大的应用价值和发展潜力^[1-2]。其中，metabarcoding片段的选择对该技术的有效使用至关重要^[2-3]。经过众多研究者的不懈努力，目前已筛选出一些物种识别率较高、对应引物扩增能力较强的优质片段^[3]。如Miya等^[4]基于880种硬骨鱼类线粒体12S序列设计的MiFish-U，扩增片段长度控制在163～185 bp范围内，对冲绳水族馆的鱼类检出数目高达230种；Vences等^[5]针对线粒体16S设计的脊椎动物Vert-16S，不仅扩增片段长度更长（约260 bp），而且对鱼类、两栖类均具有较好的适用性—以德国境内的两条河流为验证水域时，物种检出率高达82%～93%；Balasingham等^[6]设计的PS1，将扩增片段选定于线粒体COI基因上，在获得较多鱼类可操作分类单元（Operational Taxonomic Units，OTUs）的同时，可高效利用公共数据库内的参考序列进行物种注释。截至2020年，被报道的鱼类eDNA metabarcoding片段已多达数十个，其中表现优异的片段约10余个^[3]。各片段的筛选背景不同，实际表现也各具特色，这为丰富鱼类eDNA metabarcoding库、推动鱼类多样性研究奠定了坚实基础。

然而，迄今为止的鱼类eDNA研究尚不健全，如何根据研究水域和目标类群合理选择metabarcoding片段仍处于探索阶段^[3]。其中，metabarcoding片段对近缘物种的适用性就是一个亟需解决的问题。物种识别率和引物通用性是选择鱼类eDNA metabarcoding片段的两大核心指标^[3-4]。为了同时兼顾上述两方面需求，主流的eDNA metabarcoding片段（AcMDB07、Ac12S、MiFish-U、Tele02等）主要位于线粒体核糖体（12S、16S）高变区^[3]，虽然能基本满足鱼类多样性调查需求，但却存在着少数近缘鱼类因条形码片段变异速率较低而难以识别的风险。基于已有物种监测结果可知，采用核糖体区片段进行鱼类eDNA调查可能会造成“物种丢失”。如MiFish-U对平鲉属（Sebastes）、东方鲀属（Takifugu）、蝴蝶鱼属（Chaetodon）^[7]，Tele01对玉筋鱼属（Ammodytes）、杜父鱼属（Cottus）、雅罗鱼属（Leuciscus）^[8]等类群的OTUs注释结果均出现了此问题。表明“物种丢失”可能是一种比较常见的现象。然而，由于本底资料缺乏、物种数量不足等原因，已知的鱼类eDNA调查只是简单提及了此现象，未能充分揭露所用片段的近缘物种识别缺陷；不同片段间的近缘物种识别差异更缺乏相互比较。究竟哪些片段的“物种丢失”风险更低，更适合近缘鱼类多样性调查？目前尚不得而知。

随着水体污染、过度捕捞、栖息地丧失等威胁的加剧，世界各地土著鱼类的物种多样性、遗传多样性正急剧下降^[9-11]；而外来鱼种却因水产养殖、观赏活动释放等原因大幅增加^[11]。这就对当前的鱼类多样性监测提出了更高的要求。而全面、准确的物种调查则是开展各类生物多样性保护的基础和前提^[1-3]。在近缘鱼类识别是鱼类eDNA调查无法回避的大背景下，本研究将系统探究主流metabarcoding片段的“物种丢失”风险、比较不同片段的近缘物种识别差异，以期能够从中筛选出物种识别率最高的片段，最终为完善鱼类eDNA metabarcoding技术、提高鱼类多样性调查结果的准确性提供一定技术支撑。

2. 材料与方法

2.1 序列下载及筛选

自然水域难以出现近缘鱼类大规模出现的情况，结果具有较大的偶然性。因此，本研究采用MitoFish公共数据库内的线粒体全序比较各片段的近缘物种识别差异。具体操作如下：登录MitoFish数据库（http://mitofish.aori.u-tokyo.ac.jp/download.html），下载“All Mitogenomes”选项内的3 081种鱼类的线粒体全序fasta压缩文件（截至2022年1月23日）；同步下载txt文档，以方便核对序列来源、GenBank号、上传者单位等信息。使用MegAlign软件将下载的fasta文件转换为seq文件。基于Fishbase数据库（https://www.fishbase.de/）核对物种分类有效性。剔除杂交种、争议种及含有简并碱基的序列。为降低序列截齐的工作量，本研究仅将物种较多（5种及以上）的属挑出用于后续分析。物种分类阶元的确认依据Fishbase数据库及中国动物主题数据库（China Animal Scientific Database，http://www.zoology.csdb.cn/）。

2.2 metabarcoding片段选择

查阅2008−2021年期间涉及鱼类eDNA metabarcoding技术的有关研究，特别是参考Zhang等^[3]的研究结果对其中表现较好的metabarcoding片段及对应引物进行汇总。以花鳗鲡（Anguilla marmorata）线粒体全序（序列号：NC006540）为参照，核对正、反向引物方向标注是否正确。按所在位置对metabarcoding片段进行编号、排序。

2.3 序列截齐

使用DNAstar软件包中的Seqman程序，对挑选出的序列进行截齐：（1）扩增片段序列截齐后用于近缘鱼类条形码相似度分析；（2）正、反向引物序列截齐后用于引物变异度分析。

2.4 遗传距离分析

采用PAUP 4.0软件构建最大简约法（Maximum Parsimony，MP）系统发育树，系统树可信度均采用 Bootstrap检验，经1000次重复抽样检验得到分支树节点支持率。物种的聚类情况即物种识别率^[12]：若不同物种聚为一支（取遗传距离等于0），认为该物种识别失败；若不同物种聚为不同支（遗传距离大于0），认为该物种识别成功（由于不同metabarcoding片段的种间差异标准不同，难以给出准确可靠的判定阈值，因此本研究只探究遗传距离是否为0的情况）；同时，基于正、反向引物联合序列计算总平均遗传距离^[7]，从而判断引物的通用程度。

运用SPSS25.0 软件基于单因素方差分析（One-way ANOVA）对不同引物的物种识别率等进行差异显著性分析。显著差异阈值为p=0.05，试验结果以“平均值±标准差”表示。

2.5 非度量多维尺度分析

使用非度量多维尺度（Non-Metric Multidimensional Scaling，NMDS）^[13]方法对物种识别结果[0,1]矩阵进行排序分析，该分析在Canoco 5.0软件中完成。选用Bray-Curtis距离，因为其在大范围和小范围的坐标轴上都具有稳健性^{[12, 14]}。分析结果以胁强系数（stress）作为评判标准：当stress<0.2时，认为可以用NMDS的二维点图表示，该图形有一定的解释意义；当stress<0.1时，认为该排序是一个好的排序；而当stress<0.05时，则认为该排序结果具有很好的代表性^[15]。NMDS运算的步骤如下：（1）将物种识别率处理为[0,1]矩阵，鉴定成功为1，鉴定失败为0；以物种名为首行、metabarcoding片段名为首列储存.xlsx格式的[0,1]矩阵；（2）打开Canoco 5.0软件，加载.xlsx文档，给定table name，选择“import all species as factors”，其余选择默认参数，完成矩阵导入；（3）“analysis”选项下，选择“Canoco Adviser”，调出NMDS分析程序；（4）“treatment of ties in distance”选择“secondary”， “stress formula”选择“type 1”，其余选择默认参数。

3. 结果

3.1 用于分析的近缘鱼类简介

本研究用于分析的鱼类共计2纲、20目、52科、106属、935种。其中鲤形目（Cypriniformes）种类最多，共计3科、37属、337种，鲈形目（Perciformes）次之，共计18科、22属、190种。扁鲨目（Squatiniformes）、脂鲤目（Characiformes）、胡瓜鱼目（Osmeriformes）的数量最少，各只有1科、1属、5种。物种名称及对应Genbank序列号见表S1、表S2。

3.2 引物简介

经文献查阅及汇总，特别是参照Zhang等^[3]对22个主流片段（引物）的全面比较结果，本研究用于分析的鱼类eDNA metabarcoding片段共计15个。涵盖线粒体12S、16S、COI 3个目标基因。其中以12S为目标基因的片段最多（8个），16S次之（6个），COI最少（仅1个）。各片段、对应引物位置及信息见图1及表1。

图 1 本研究中15个鱼类eDNA metabarcoding片段和对应引物在目标基因上的位置

以花鳗鲡（Anguilla marmorata）线粒体基因序列（GenBank登录号：NC006540）为标注模板

Figure 1. Locations of the 15 fish eDNA metabarcoding fragments and primer pairs on the target mitochondrial genes

Sequence of the marbled eel (Anguilla marmorata) (GenBank accession number: NC006540) was used as template

下载: 全尺寸图片幻灯片

表 1 本研究中15个鱼类eDNA metabarcoding片段简介

Table 1. Summary of 15 fish eDNA metabarcoding fragments analyzed in this study

metabarcoding片段	对应引物名称	目标基因	长度/bp	引物序列（5'-3'）	设计者
片段01	MiFish-U	12S	169	For-GTCGGTAAAACTCGTGCCAGC	Miya等^[4]
片段01	MiFish-U	12S	169	Rev-CATAGTGGGGTATCTAATCCCAGTTTG	Miya等^[4]
片段02	Tele02	12S	165	For-AAACTCGTGCCAGCCACC	Taberlet等^[16]
片段02	Tele02	12S	165	Rev-GGGTATCTAATCCCAGTTTG	Taberlet等^[16]
片段03	Am12S	12S	242	For-AGCCACCGCGGTTATACG	Evans等^[17]
片段03	Am12S	12S	242	Rev-CAAGTCCTTTGGGTTTTAAGC	Evans等^[17]
片段04	12S-V5	12S	107	For-ACTGGGATTAGATACCCC	Riaz等^[18]
片段04	12S-V5	12S	107	Rev-TAGAACAGGCTCCTCTAG	Riaz等^[18]
片段05	Ac12S	12S	392	For-ACTGGGATTAGATACCCCACTATG	Evans等^[17]
片段05	Ac12S	12S	392	Rev-GAGAGTGACGGGCGGTGT	Evans等^[17]
片段06	AcMDB07	12S	282	For-GCCTATATACCGCCGTCG	Bylemans等^[19]
片段06	AcMDB07	12S	282	Rev-GTACACTTACCATGTTACGACTT	Bylemans等^[19]
片段07	NeoFish_3	12S	196	For-CGCCGTCGCAAGCTTACCCT	Milan等^[20]
片段07	NeoFish_3	12S	196	Rev-AGTGACGGGCGGTGTGTGC	Milan等^[20]
片段08	Tele01	12S	64	For-ACACCGCCCGTCACTCT	Valentini等^[8]
片段08	Tele01	12S	64	Rev-CTTCCGGTACACTTACCATG	Valentini等^[8]
片段09	Ac16S	16S	339	For-CCTTTTGCATCATGATTTAGC	Evans等^[[17]
片段09	Ac16S	16S	339	Rev-CAGGTGGCTGCTTTTAGGC	Evans等^[[17]
片段10	L2513/H2714	16S	206	For-GCCTGTTTACCAAAAACATCAC	Kitano等^[21]
片段10	L2513/H2714	16S	206	Rev-CTCCATAGGGTCTTCTCGTCTT	Kitano等^[21]
片段11	Fish16S	16S	84	For-CGAGAAGACCCTWTGGAGCTTIAG	Shaw等^[22]
片段11	Fish16S	16S	84	Rev-GGTCGCCCCAACCRAAG	Shaw等^[22]
片段12	Ve16S	16S	335	For-CGAGAAGACCCTATGGAGCTTA	Evans等^[17]
片段12	Ve16S	16S	335	Rev-AATCGTTGAACAAACGAACC	Evans等^[17]
片段13	Vert-16S-eDNA	16S	284	For-AGACGAGAAGACCCYDTGGAGCTT	Vences等^[5]
片段13	Vert-16S-eDNA	16S	284	Rev-GATCCAACATCGAGGTCGTAA	Vences等^[5]
片段14	Fish16S F/D-2R	16S	222	For-GACCCTATGGAGCTTTAGAC	DiBattista等^[23]
片段14	Fish16S F/D-2R	16S	222	Rev-CGCTGTTATCCCTADRGTAACT	DiBattista等^[23]
片段15	PS1	COI	199	For-ACCTGCCTGCCGTATTTGGYGCYTGRGCCGGRATAGT	Balasingham等^[6]
片段15	PS1	COI	199	Rev-ACGCCACCGAGCCARAARCTYATRTTRTTYATTCG	Balasingham等^[6]
注：片段长度以花鳗鲡（Anguilla marmorata）（GenBank登录号：NC006540）为参考。

下载: 导出CSV

| 显示表格

3.3 物种识别率及引物变异情况

不同片段的物种识别率存在极显著差异（F=32.39，p=1.5×10⁻⁷⁶），蛋白质编码基因（片段15）的物种识别率最高，达到（82.76±24.66）%。核糖体基因中，片段03、片段05的表现最为优异，识别率分别为（68.9±30.81）%和（71.62±27.3）%；片段04、片段08和片段11的识别率较低，平均值未超过42.92%；其余片段的物种识别率较为中等，平均值在54.58%～64.6%间（图2，表2）。

图 2 本研究中15个鱼类eDNA metabarcoding片段的物种识别率及对应引物序列总平均遗传距离

Figure 2. Fish species resolution rates of 15 eDNA metabarcoding fragments and overall mean distances of primer pairs in this study

下载: 全尺寸图片幻灯片

表 2 metabarcoding片段对各类群的识别率

Table 2. The species resolution rates of metabarcoding fragments for different groups

物种类群	物种数	识别率/%
物种类群	物种数	片段 01	片段 02	片段 03	片段 04	片段 05	片段 06	片段 07	片段 08	片段 09	片段 10	片段 11	片段12	片段 13	片段 14	片段15	全部片段
鱊属 Acheilognathus	14	71.43	71.43	85.71	71.43	85.71	85.71	71.43	57.14	71.43	85.71	57.14	85.71	85.71	71.43	100	77.14±11.83
鲟属 Acipenser	15	13.33	20	26.67	13.33	46.67	33.33	40	6.67	33.33	0	33.33	6.67	13.33	6.67	73.33	24.44± 19.46
光唇鱼属 Acrossocheilus	16	56.25	50	62.5	50	87.5	75	62.5	6.25	68.75	43.75	50	62.5	50	43.75	100	57.92± 21.45
黑头鱼属 Alepocephalus	7	57.14	14.29	28.57	0	28.57	28.57	28.57	0	28.57	14.29	14.29	28.57	42.86	28.57	100	29.52± 24.43
西鲱属 Alosa	6	50	0	0	0	33.33	33.33	33.33	0	0	66.67	16.67	16.67	16.67	16.67	100	25.56± 28.78
双锯鱼属 Amphiprion	10	70	50	60	20	60	50	30	10	60	60	30	80	60	50	40	48.67± 19.22
鳗鲡属 Anguilla	17	29.41	35.29	64.71	0	76.47	70.59	76.47	41.18	64.71	47.06	35.29	76.47	64.71	52.94	88.24	54.9± 23.5
细鲫属 Aphyocypris	5	60	60	60	20	60	40	40	0	100	80	60	100	80	100	100	64± 30.43
旗鳉属 Aphyosemion	8	100	100	100	100	100	100	100	12.5	100	100	100	100	100	100	100	94.17± 22.59
深海鳐属 Bathyraja	8	12.5	25	50	12.5	25	12.5	12.5	0	12.5	12.5	0	25	25	37.5	50	20.83±15.43
鲫属 Carassius	5	0	0	0	0	60	20	20	0	60	40	0	40	0	0	60	20±25.07
真鲨属 Carcharhinus	14	64.29	71.43	100	78.57	100	100	100	50	85.71	35.71	14.29	50	64.29	71.43	100	72.38±26.57
刺尻鱼属 Centropyge	16	75	75	87.5	31.25	62.5	56.25	62.5	12.5	87.5	37.5	0	50	50	62.5	100	56.67±27.8
九棘鲈属 Cephalopholis	5	100	100	100	40	100	60	60	0	40	40	60	100	100	100	100	73.33±32.66
鳢属 Channa	8	62.5	62.5	62.5	62.5	62.5	62.5	62.5	25	75	75	75	75	75	75	50	64.17±13.25
红鳍鲌属 Chanodichthys	5	20	20	20	20	60	60	0	0	100	60	0	100	60	60	100	45.33±36.62
燕鳐属 Cheilopogon	8	37.5	37.5	25	0	25	50	50	0	62.5	25	12.5	12.5	25	12.5	75	30±22.06
鲮属 Cirrhinus	5	60	60	60	40	100	100	100	0	60	60	20	60	40	40	60	57.33±28.15
胡子鲇属 Clarias	6	66.67	16.67	33.33	33.33	83.33	33.33	33.33	0	66.67	33.33	0	66.67	66.67	33.33	100	44.44±29.32
花鳅属 Cobitis	11	81.82	81.82	81.82	81.82	72.73	81.82	81.82	9.09	63.64	81.82	63.64	81.82	81.82	81.82	100	75.15±20.17
鲚属 Coilia	6	16.67	50	50	16.67	66.67	100	66.67	0	50	66.67	33.33	33.33	16.67	16.67	100	45.56±30.52
白鲑属 Coregonus	14	7.14	7.14	0	0	21.43	35.71	21.43	7.14	0	0	0	0	0	0	7.14	7.14±10.8
兵鲶属 Corydoras	6	33.33	33.33	33.33	16.67	50	50	66.67	16.67	50	50	16.67	0	0	0	66.67	32.22±23.12
杜父鱼属 Cottus	12	50	75	100	50	83.33	50	41.67	8.33	83.33	83.33	41.67	83.33	66.67	100	100	67.78±26.7
舌鳎属 Cynoglossus	12	100	100	100	83.33	100	100	100	83.33	100	83.33	83.33	83.33	83.33	83.33	100	92.22±8.61
长吻鳐属 Dipturus	8	0	12.5	37.5	12.5	50	37.5	37.5	0	50	12.5	0	25	25	25	62.5	25.83±19.75
盘属 Discogobio	5	40	40	60	40	60	40	60	0	20	40	0	40	60	40	60	40±20
塘鳢属 Eleotris	6	66.67	66.67	100	100	100	100	100	0	100	100	66.67	66.67	66.67	66.67	100	80±27.6
Enteromius属	5	100	100	100	100	100	100	100	40	100	100	100	100	80	100	100	94.67±15.98
石斑鱼属 Epinephelus	22	90.91	90.91	90.91	50	90.91	90.91	90.91	36.36	72.73	90.91	22.73	81.82	68.18	63.64	90.91	74.85±22.52
扁鳉属 Epiplatys	7	71.43	71.43	71.43	71.43	71.43	71.43	71.43	0	71.43	71.43	71.43	71.43	71.43	71.43	71.43	66.67±18.44
镖鲈属 Etheostoma	10	100	100	100	100	100	100	100	40	100	70	90	100	100	100	100	93.33±16.76
底鳄鳉属 Fundulopanchax	5	100	100	100	100	100	100	100	20	100	100	60	100	100	100	100	92±22.42
底鳉属 Fundulus	5	60	60	60	20	100	100	100	0	100	100	100	100	100	100	100	80±32.95
南乳鱼属 Galaxias	5	100	100	100	60	100	100	100	0	100	60	60	100	100	100	100	85.33±28.75
墨头鱼属 Garra	10	100	100	100	50	100	100	100	30	80	100	50	100	100	100	100	87.33±23.74
骨尾鱼属 Gila	5	40	20	60	20	60	40	40	0	60	60	20	60	60	60	100	46.67±24.69
露齿鲨属 Glyphis	8	25	12.5	25	0	37.5	37.5	37.5	0	25	25	25	25	25	37.5	25	24.17±12.01
纹胸属 Glyptothorax	11	18.18	18.18	63.64	18.18	36.36	27.27	36.36	9.09	54.55	9.09	27.27	54.55	54.55	63.64	81.82	38.18±22.59
颌须属 Gnathopogon	6	33.33	66.67	100	100	100	66.67	66.67	0	100	50	16.67	33.33	33.33	33.33	66.67	57.78±32.65
鳅鮀属 Gobiobotia	6	100	50	100	50	100	100	100	0	66.67	66.67	16.67	83.33	83.33	66.67	100	72.22±31.91
裸鲤属 Gymnocypris	8	37.5	37.5	37.5	0	12.5	12.5	0	0	12.5	37.5	12.5	12.5	12.5	0	25	16.67±14.69
海猪鱼属 Halichoeres	5	100	100	100	100	100	100	100	20	100	100	100	100	100	100	100	94.67±20.66
鳠属 Hemibagrus	5	60	60	60	80	60	40	60	0	60	60	20	60	60	60	60	53.33±19.52
䱻属 Hemibarbus	7	42.86	42.86	42.86	28.57	42.86	28.57	14.29	0	42.86	28.57	14.29	42.86	42.86	42.86	100	37.14±22.13
海马属 Hippocampus	18	88.89	88.89	88.89	55.56	66.67	66.67	61.11	27.78	100	77.78	55.56	100	100	100	100	78.52±22.11
刺蝶鱼属 Holacanthus	5	40	40	60	40	60	60	60	0	40	20	20	40	40	40	60	41.33±17.67
下鱵属 Hyporhamphus	5	100	100	100	100	100	100	100	40	100	100	100	100	100	100	100	96±15.49
野鲮属 Labeo	20	70	70	80	35	90	80	75	10	60	30	25	55	40	50	80	56.67±24.1
兔头鲀属 Lagocephalus	8	100	100	100	75	100	100	100	25	100	75	75	62.5	75	62.5	100	83.33±21.99
太阳鱼属 Lepomis	6	100	100	100	50	100	100	100	0	100	100	66.67	100	100	100	100	87.78±28.5
薄鳅属 Leptobotia	6	50	16.67	66.67	50	66.67	50	33.33	0	33.33	16.67	50	33.33	33.33	33.33	50	38.89±18.55
䱀属 Liobagrus	9	66.67	77.78	77.78	11.11	77.78	55.56	66.67	0	33.33	33.33	33.33	100	100	100	77.78	60.74±31.95
笛鲷属 Lutjanus	13	84.62	84.62	100	100	100	100	84.62	23.08	100	69.23	30.77	100	100	100	100	85.13±25.42
狼绵鳚属 Lycodes	5	60	60	60	40	60	60	60	0	100	40	20	60	60	60	100	56±25.3
虹银汉鱼属 Melanotaenia	6	100	100	50	0	66.67	33.33	33.33	0	66.67	0	50	33.33	100	66.67	100	53.33±36.84
小鳔属 Microphysogobio	14	64.29	64.29	85.71	64.29	71.43	50	35.71	7.14	85.71	35.71	35.71	57.14	71.43	57.14	100	59.05±23.76
泥鳅属 Misgurnus	5	100	100	100	20	60	60	40	0	60	60	0	60	100	60	100	61.33±35.02
假鳃鳉属 Nothobranchius	8	100	100	100	100	100	100	100	12.5	100	100	100	100	100	100	100	94.17±22.59
美洲鱥属 Notropis	10	100	80	100	60	100	100	80	10	100	60	70	100	100	100	100	84±25.58
沙塘鳢属 Odontobutis	6	100	100	100	100	100	100	66.67	33.33	100	100	100	100	100	100	100	93.33±18.69
太平洋鲑属 Oncorhynchus	8	37.5	62.5	100	37.5	100	75	75	0	37.5	0	0	50	50	50	100	51.67±34.68
白甲鱼属 Onychostoma	11	81.82	63.64	45.45	18.18	63.64	63.64	54.55	27.27	54.55	72.73	0	81.82	54.55	45.45	63.64	52.73±22.85
马口鱼属 Opsariichthys	5	100	120	100	40	60	60	20	0	60	60	40	100	100	100	100	70.67±35.35
罗非鱼属 Oreochromis	9	66.67	66.67	66.67	33.33	44.44	11.11	44.44	0	55.56	22.22	22.22	44.44	33.33	22.22	77.78	40.74±22.88
青鳉属 Oryzias	10	80	80	80	80	80	80	80	80	100	100	100	100	100	100	100	89.33±10.33
纹唇鱼属 Osteochilus	5	80	100	100	80	100	100	100	0	100	100	80	100	100	100	100	89.33±26.04
鲳属 Pampus	7	28.57	57.14	28.57	14.29	28.57	28.57	28.57	0	28.57	28.57	14.29	14.29	14.29	14.29	42.86	24.76±13.73
属 Pareuchiloglanis	5	60	60	60	20	60	100	60	0	60	100	60	60	60	60	100	61.33±26.69
鱥属 Phoxinus	7	100	100	100	57.14	100	71.43	71.43	0	71.43	71.43	14.29	100	100	100	100	77.14±32.31
花鳉属 Poecilia	5	20	20	20	20	60	20	40	20	20	20	20	20	20	20	60	26.67±14.47
多鳍鱼属 Polypterus	10	100	100	100	100	100	100	100	40	100	80	60	100	100	100	100	92±18.21
原唇齿脂鲤属 Prochilodus	5	20	20	20	0	20	0	40	0	0	0	0	40	40	40	20	17.33±16.68
拟鲿属 Pseudobagrus	12	41.67	58.33	75	33.33	58.33	58.33	58.33	16.67	66.67	8.33	8.33	33.33	25	25	83.33	43.33±24.03
拟腹吸鳅属 Pseudogastromyzon	9	33.33	55.56	55.56	0	55.56	55.56	55.56	0	77.78	55.56	44.44	100	77.78	55.56	77.78	53.34±26.96
多刺鱼属 Pungitius	7	71.43	100	100	28.57	100	57.14	57.14	0	28.57	28.57	0	71.43	71.43	42.86	100	57.14±34.99
鳐属 Raja	6	16.67	16.67	33.33	16.67	66.67	16.67	16.67	0	50	66.67	16.67	66.67	66.67	66.67	66.67	38.89±25.72
波鱼属 Rasbora	7	100	100	100	71.43	100	71.43	100	14.29	71.43	100	57.14	71.43	71.43	71.43	100	80±24.03
吻鰕虎鱼属 Rhinogobius	5	100	100	100	100	100	100	100	0	100	100	100	100	100	100	100	93.33±25.82
鳑鲏属 Rhodeus	10	80	80	80	30	80	80	60	20	100	80	50	100	100	80	100	74.67±24.75
红点鲑属 Salvelinus	12	8.33	8.33	25	8.33	25	25	25	0	16.67	8.33	0	8.33	0	8.33	58.33	15±15.17
鳈属 Sarcocheilichthys	9	55.56	44.44	55.56	11.11	55.56	55.56	33.33	0	55.56	33.33	0	44.44	33.33	44.44	55.56	38.52±20.08
沙丁鱼属 Sardinella	6	100	100	100	66.67	100	50	66.67	33.33	50	66.67	66.67	100	100	66.67	100	77.78±23.29
蛇属 Saurogobio	7	100	100	100	71.43	100	71.43	71.43	0	71.43	71.43	14.29	100	100	100	100	78.1±31.88
南鳅属 Schistura	11	100	100	100	100	100	100	100	54.55	100	100	63.64	100	100	100	100	94.55±14.5
裸裂尻鱼属 Schizopygopsis	7	14.29	14.29	14.29	0	14.29	28.57	14.29	0	28.57	0	0	42.86	42.86	42.86	71.43	21.91±20.82
裂腹鱼属 Schizothorax	26	11.54	11.54	15.38	7.69	23.08	23.08	11.54	3.85	15.38	7.69	0	11.54	3.85	11.54	30.77	12.56±8.15
马鲛鱼属 Scomberomorus	6	66.67	66.67	66.67	33.33	66.67	66.67	66.67	33.33	100	66.67	100	66.67	66.67	66.67	100	68.89±19.79
平鲉属 Sebastes	19	31.58	31.58	31.58	10.53	36.84	47.37	26.32	0	10.53	15.79	0	26.32	15.79	15.79	68.42	24.56±18.09
瓢鳍虾虎鱼属 Sicyopterus	18	33.33	33.33	44.44	22.22	50	66.67	50	33.33	83.33	33.33	33.33	66.67	55.56	66.67	88.89	50.74±20.13
篮子鱼属 Siganus	6	33.33	33.33	33.33	16.67	100	66.67	66.67	0	66.67	16.67	0	16.67	0	0	33.33	32.22±30.52
鱚属 Sillago	6	100	100	100	100	100	100	100	50	100	100	66.67	100	100	100	100	94.44±15
鲶属 Silurus	6	66.67	66.67	66.67	50	100	100	66.67	16.67	66.67	50	50	66.67	66.67	100	100	68.89±23.46
鳜属 Siniperca	7	57.14	42.86	100	28.57	42.86	28.57	28.57	0	42.86	28.57	14.29	14.29	14.29	14.29	28.57	32.38±23.82
金线鲃属 Sinocyclocheilus	14	85.71	85.71	85.71	35.71	78.57	78.57	64.29	28.57	50	57.14	21.43	64.29	57.14	78.57	85.71	63.81±21.84
银属 Squalidus	7	42.86	28.57	57.14	42.86	100	57.14	57.14	14.29	57.14	28.57	14.29	42.86	28.57	28.57	100	46.67±26.16
扁鲨属 Squatina	5	60	60	60	20	60	60	60	0	20	40	40	60	60	60	100	50.67±23.74
疯鲿属 Tachysurus	8	75	75	75	25	100	100	100	0	75	37.5	0	75	50	50	75	60.83±33.03
东方鲀属 Takifugu	16	0	6.25	6.25	6.25	37.5	25	12.5	12.5	12.5	18.75	6.25	31.25	25	31.25	62.5	19.58±16.34
田中鳑鲏属 Tanakia	6	100	100	100	100	100	100	100	16.67	100	100	33.33	100	100	100	100	90±26.58
鲀属 Tetraodon	5	100	100	100	60	60	60	60	0	100	100	0	100	60	60	100	70.67±34.53
棱鳀属 Thryssa	6	100	100	100	66.67	100	100	100	0	100	100	100	100	100	100	100	91.11±26.63
金枪鱼属 Thunnus	9	0	11.11	11.11	0	0	0	0	0	0	0	0	0	0	0	0	1.48±3.91
茴鱼属 Thymallus	10	20	20	40	0	30	50	20	0	30	20	10	20	30	10	60	24±16.82
结鱼属 Tor	7	71.43	71.43	100	0	71.43	42.86	57.14	0	28.57	57.14	28.57	71.43	57.14	28.57	100	52.38±30.86
高原鳅属 Triplophysa	21	47.62	47.62	80.95	42.86	71.43	71.43	61.9	4.76	66.67	52.38	19.05	47.62	47.62	47.62	80.95	52.7±21.1
全部物种	935	62.3±31.84	61.95±32.54	68.9±30.81	42.92±33.46	71.62±27.3	64.34±29.13	60.89±29.48	11.79±17.79	64.6±30.65	54.58±32.28	36.71±32.27	64.2±31.8	61.6±32.35	59.71±32.78	82.76±24.66

下载: 导出CSV

| 显示表格

基于15个片段对应正、反向引物联合序列构建的MP系统发育树总平均遗传距离依次为1.52%、1.84%、4.58%、0.88%、0.97%、2.46%、5.78%、3.32%、9.32%、2.29%、6.21%、4.14%、1.71%、3.17%和15.70%（图2）。表明蛋白质编码基因（片段15）的对应引物通用性最差。片段09、片段11、片段07、片段03和片段12的引物变异程度也较高（总平均遗传距离≥4.14%），实际使用过程中应考虑其eDNA扩增效率。

不同类群的识别率也存在极显著差异（F=14.40, p=7×10⁻¹⁵⁷）：旗鳉属（Aphyosemion）、舌鳎属（Cynoglossus）、镖鲈属（Etheostoma）等属内近缘物种较容易区分，全部15个片段的平均识别率均在90%以上；而金枪鱼属（Thunnus）、白鲑属（Coregonus）等类群的目标基因同源性较高，大部分片段对其识别能力较差（图3，表2）。

图 3 本研究中106属鱼类的物种识别率

Figure 3. Fish species resolution rates of 106 genera in this study

下载: 全尺寸图片幻灯片

3.4 NMDS结果

根据935种鱼类的识别结果对15个metabarcoding片段做NMDS排序，片段04、片段08和片段11明显远离其余片段，说明此3个片段物种识别结果较为不同（图4）。其余11个核糖体片段根据所在目标基因聚类成A、B两组，表明不同目标基因影响物种识别。NMDS分析结果的胁强系数为0.13，说明将15个metabarcoding片段划归不同组具有一定的解释意义。

图 4 不同metabarcoding片段的非度量多维尺度（NMDS）分析

F01表示片段01，其余依次类推

Figure 4. Analysis of non-metric multidimensional scaling (NMDS) for different metabarcoding fragments

F01 represents fragment 01, and the rest are followed by analogy

下载: 全尺寸图片幻灯片

4. 讨论

本研究用于分析的鱼类数量庞大、阶元众多（见3.1节及表S1、表S2），涉及软骨鱼类及硬骨鱼类、淡水鱼类及海水鱼类、定居性鱼类及洄游性鱼类等各类群。样本量具有较好的统计学意义，结果能够比较充分地反映15个主流eDNA metabarcoding片段的优劣。本研究是目前国内第一篇针对性探究鱼类eDNA metabarcoding近缘物种适用性的报告。

4.1 COI片段的引物缺陷

片段15（PS1，COI）的物种识别率显著高于剩余14个12S/16S片段，是唯一物种识别率超过80%的metabarcoding片段。这与Balasingham等^[6]认为的鱼类线粒体蛋白质编码基因比核糖体基因能够更有效避免北美五大湖内部分入侵物种eDNA注释结果假阴性的结论一致，也与经典的DNA条形码普遍采用线粒体细胞色素C氧化酶亚基I基因的事实相符^[24]。此外，COI基因在鱼类eDNA metabarcoding参考数据库的全面性方面也具有12S、16S基因难以比拟的优势^[3]。然而，引物序列的高变性可能会阻碍片段15的广泛应用。本研究中该片段正、反向引物序列的总平均遗传距离高达15.70%。中国团扇鳐（Platyrhina sinensis）、细条银口天竺鲷（Jaydia lineata）、土佐鰧（Uranoscopus tosae）、横带髭鲷（Hapalogenys analis）的COI特异性引物、探针设计结果表明：当总平均遗传距离大于8.21%～9.65%时，高纯度DNA模板会出现明显的弥散扩增；总平均遗传距离大于12.47%～13.77%时，引物会对部分近缘物种产生绝对特异^[7]。因此，片段15存在eDNA低效扩增/偏倚扩增的缺陷。Collins等^[25]、Zhang等^[3]分别对英吉利海、北京水系的鱼类多样性调查显示，COI片段的物种检出数目不及所用的12S片段，其结果支持本研究观点。

此外，其他的COI或Cytb metabarcoding片段可能也难以解决引物通用性的问题。鱼类eDNA metabarcoding技术需要依托2×150 bp、2×250 bp两种高通量测序（Next-generation Sequencing，NGS）平台实现^[3]，而Miya等^[4]对880种硬骨鱼类线粒体基因组全序的筛查结果显示，COI、Cytb等蛋白质编码基因在该读长限定范围内不存在成对的侧翼保守区，理论上难以设计出与经典的鱼类DNA条形码（基于Sanger测序，读长可达2×600 bp）媲美的成对引物。Menning等^[26]、Jennings等^[27]、Sultana等^[28]的序列比对结果也表明，整个Cytb基因基本为高变序列，COI基因则只在近5'端前350 bp范围内变异较小，相对适合设计短片段通用引物。结合本研究（图1）和Zhang等^[3]对22对引物的汇总结果可知，片段15等少数几对鱼类COI metabarcoding优质片段恰位于此相对保守区，其引物通用性较差，而其他区域metabarcoding片段的引物序列总平均遗传距离可能更大。因此，本研究认为后续研究没有必要进一步比较其他COI、Cytb片段的近缘物种识别差异；同时认为以片段15为代表的蛋白质编码片段更适合作为eDNA研究的辅助metabarcoding，并且需要根据研究水域、目标类群降低引物调查范围，从而对引物序列进行针对性设计或优化。

4.2 核糖体片段的比较

片段长度影响物种识别（表1，图2，图4）。本研究中识别率最低的3个metabarcoding片段（片段08、片段11和片段04）长度均不超过110 bp，而片段05、片段03等高识别率片段则在200 bp以上。Balasingham等^[6]、Gantner等^[29]的研究也认为12S和16S扩增片段越短，物种鉴定的准确性越低。虽然由于水体中的痕量eDNA通常存在严重降解，主流的观点认为200 bp以内的微条形码可能具有更高的PCR成功率^{[4, 30]}，但Bylemans等^[19]、Zhang等^[3]多个研究却发现：片段05、片段03等的高通量测序数据量和物种检出数目并不低于片段08、片段04或片段11，只是在定量分析等方面相互存在较大差异。鱼类eDNA的产生、降解动力学过程极为复杂，目前没有研究能够准确阐明eDNA的微观存在、分布和变化规律，Deiner等^[31]、Bylemans等^[19]推测长片段eDNA可能反而可以更长久地存在水体中。本研究不推荐使用metabarcoding短片段（片段08、片段11和片段04）进行近缘鱼类多样性调查。

除去COI和3个短片段，片段05的物种识别率最高，其引物通用性（总遗传距离=0.97%）也仅次于片段04。片段05是近缘鱼类多样性调查的第一选择，该结果与Zhang等^[3]的研究结论一致。片段09、片段07、片段03、片段12的物种识别率虽然也较高，但在本研究中，其引物序列总平均遗传距离较大（≥4.14%），存在类似于片段15的“物种丢失”风险；其他metabarcoding引物比较研究也表明其物种检出数目不及片段06、片段13^[3]。片段14、片段10虽然长度大于200 bp，但物种识别率和引物通用性却并不出众；并且可能因为引物序列3'位置存在变异^[7]，这两个片段在个别研究中的实际表现甚至不如片段15^[3]。而片段01和片段02则在引物序列、条形码性能方面基本一致，二者存在明显的相互替代性，本研究更倾向二者中选用片段01。综上所述，片段05、片段06、片段01、片段13是本研究筛选出的4个最优metabarcoding片段。然而，受高通量测序成本和参考数据库的影响，在单一研究中使用的metabarcoding片段通常不超过3个^[3-4]。虽然本研究及Zhang等^[3]、Bylemans等^[19]的研究均表明片段06的物种识别率和引物通用性优异，但NMDS分析（图4）显示片段06在散点图上与片段05位置较近，因此同时使用这两个片段进行鱼类多样性研究可能并不能显著提高物种检出数目。NMDS分析也显示，不同基因、同一基因不同片段间均存在较大差异（图4）。因此不仅需要多片段联合应用，而且所用片段要有足够的多基因代表性。本研究倾向于以片段05、片段01为主，片段13等为辅，进行近缘鱼类多样性调查。

4.3 本研究所用片段的不足及新片段开发的可能性

本研究对106属鱼类的识别结果差异极显著（见3.3节），表明物种类群直接影响eDNA metabarcoding调查效果。从15个片段的总识别率角度比较（表2），仅15属鱼类物种识别率超过90%，多达38属鱼类识别率不足50%。由于935种鱼类的线粒体全序列皆来自公共数据库，其中一些序列可能存在同物异名错误—如鲳属鱼类的银鲳（Pampus argenteus）与镰鲳（Pampus echinogaster）^[32]，因此该属鱼类识别率偏低；但更可能的解释是低识别类群物种间遗传差异较小，导致metabarcoding片段从源头上就难以完全区分—如白鲑属（Coregonus）^[33]、红点鲑属（Salvelinus）^[34]的众多物种集中分布于50°N以北环北极圈水域，物种间可能存在较为频繁的基因交流；平鲉属（Sebastes）则广泛分布于北太平洋且仔鱼营漂浮生活，不同群体间也不存在明显的地理隔离^[35]。Miya等^[4]在筛选MiFish-U（片段01）引物时已发现，金枪鱼属等大洋性鱼类的大部分线粒体基因高度保守，只有借助变异速率极高的线粒体NADH脱氢酶亚基5（NADH Dehydrogenase Subunit 5，ND5）条形码才能进行有效区分。因此，受制于物种自身的遗传背景，本研究所用片段尚无法对全部鱼类进行100%区分。后续仍需不断尝试筛选新片段从而尽量提高eDNA metabarcoding技术的鱼类多样性调查能力。

metabarcoding片段的开发与目标基因内茎环结构的数量和分布情况密切相关。斑鰶（Konosirus punctatus）^[36]、大口鰜（Psettodes erumei）^[37]、北极茴鱼（Thymallus arcticus）^[38]等线粒体基因二级结构示意图显示，16S基因仅在中后段存在1～2个大型茎环结构，相比之下，12S基因存在2～4个。由于12S基因的总长度仅为16S的60%（约900 bp : 1600 bp），因此该基因内的茎环结构分布更为密集和均匀。茎环结构有利于寻找稳定的“保守−高变−保守”序列，Evans和Lamberti^[39]、Hänfling等^[40]、Bylemans等^[19]据此认为，线粒体12S基因比16S基因更适合开发鱼类eDNA metabarcoding标记，本研究支持此结论。然而，基于14个核糖体片段的位置示意图（图1）和简约信息（表1）可知，鱼类16S基因的后半段大型茎环及12S基因的主要茎环皆已经被筛选出metabarcoding片段。围绕单一茎环结构开发新片段不仅可能性极小，而且片段识别率可能也不理想（如片段07、片段10）。相比片段01、片段04和片段08等，片段05和片段06等则是由多个中小型茎环结构组合形成，其物种识别率反而更高。以此类推，联合多个相邻茎环区域可能是形成新的高识别率片段的手段之一。本研究中，12S片段对应引物存在明显的位置重合及序列共用（图1，表1），说明这些区域的引物保守性得到了较大认可。因此，可以尝试通过这些引物的组合应用开发更长、更高识别率的metabarcoding片段。如片段01、片段02的正向引物与片段06、片段07正向引物的反向互补序列共用，PCR产物约420 bp（图1）。该片段由1个大型茎环加2个小型茎环结构组成，其物种识别率可能高于片段05。

综合物种识别率、引物通用性等多方面因素，本研究推荐2×150 bp测序平台使用片段01（MiFish-U）、2×250 bp测序平台使用片段05（Ac12S），辅以片段13（Vert-16S-eDNA）进行近缘鱼类多样性调查；新片段的开发也是不断完善eDNA metabarcoding技术的重要工作。然而，鱼类eDNA研究受多方面因素的影响，基于935种近缘鱼类线粒体全序的比较结果可能会与实际表现有差异。后期还需寻找近缘鱼类广泛分布的水域，在建立庞大的本底资料数据库的基础上验证本研究的结论和推测。

图 1 本研究中15个鱼类eDNA metabarcoding片段和对应引物在目标基因上的位置

以花鳗鲡（Anguilla marmorata）线粒体基因序列（GenBank登录号：NC006540）为标注模板

Fig. 1 Locations of the 15 fish eDNA metabarcoding fragments and primer pairs on the target mitochondrial genes

Sequence of the marbled eel (Anguilla marmorata) (GenBank accession number: NC006540) was used as template

下载: 全尺寸图片幻灯片

图 2 本研究中15个鱼类eDNA metabarcoding片段的物种识别率及对应引物序列总平均遗传距离

Fig. 2 Fish species resolution rates of 15 eDNA metabarcoding fragments and overall mean distances of primer pairs in this study

下载: 全尺寸图片幻灯片

图 3 本研究中106属鱼类的物种识别率

Fig. 3 Fish species resolution rates of 106 genera in this study

下载: 全尺寸图片幻灯片

图 4 不同metabarcoding片段的非度量多维尺度（NMDS）分析

F01表示片段01，其余依次类推

Fig. 4 Analysis of non-metric multidimensional scaling (NMDS) for different metabarcoding fragments

F01 represents fragment 01, and the rest are followed by analogy

下载: 全尺寸图片幻灯片

表 1 本研究中15个鱼类eDNA metabarcoding片段简介

Tab. 1 Summary of 15 fish eDNA metabarcoding fragments analyzed in this study

metabarcoding片段	对应引物名称	目标基因	长度/bp	引物序列（5'-3'）	设计者
片段01	MiFish-U	12S	169	For-GTCGGTAAAACTCGTGCCAGC	Miya等^[4]
片段01	MiFish-U	12S	169	Rev-CATAGTGGGGTATCTAATCCCAGTTTG	Miya等^[4]
片段02	Tele02	12S	165	For-AAACTCGTGCCAGCCACC	Taberlet等^[16]
片段02	Tele02	12S	165	Rev-GGGTATCTAATCCCAGTTTG	Taberlet等^[16]
片段03	Am12S	12S	242	For-AGCCACCGCGGTTATACG	Evans等^[17]
片段03	Am12S	12S	242	Rev-CAAGTCCTTTGGGTTTTAAGC	Evans等^[17]
片段04	12S-V5	12S	107	For-ACTGGGATTAGATACCCC	Riaz等^[18]
片段04	12S-V5	12S	107	Rev-TAGAACAGGCTCCTCTAG	Riaz等^[18]
片段05	Ac12S	12S	392	For-ACTGGGATTAGATACCCCACTATG	Evans等^[17]
片段05	Ac12S	12S	392	Rev-GAGAGTGACGGGCGGTGT	Evans等^[17]
片段06	AcMDB07	12S	282	For-GCCTATATACCGCCGTCG	Bylemans等^[19]
片段06	AcMDB07	12S	282	Rev-GTACACTTACCATGTTACGACTT	Bylemans等^[19]
片段07	NeoFish_3	12S	196	For-CGCCGTCGCAAGCTTACCCT	Milan等^[20]
片段07	NeoFish_3	12S	196	Rev-AGTGACGGGCGGTGTGTGC	Milan等^[20]
片段08	Tele01	12S	64	For-ACACCGCCCGTCACTCT	Valentini等^[8]
片段08	Tele01	12S	64	Rev-CTTCCGGTACACTTACCATG	Valentini等^[8]
片段09	Ac16S	16S	339	For-CCTTTTGCATCATGATTTAGC	Evans等^[[17]
片段09	Ac16S	16S	339	Rev-CAGGTGGCTGCTTTTAGGC	Evans等^[[17]
片段10	L2513/H2714	16S	206	For-GCCTGTTTACCAAAAACATCAC	Kitano等^[21]
片段10	L2513/H2714	16S	206	Rev-CTCCATAGGGTCTTCTCGTCTT	Kitano等^[21]
片段11	Fish16S	16S	84	For-CGAGAAGACCCTWTGGAGCTTIAG	Shaw等^[22]
片段11	Fish16S	16S	84	Rev-GGTCGCCCCAACCRAAG	Shaw等^[22]
片段12	Ve16S	16S	335	For-CGAGAAGACCCTATGGAGCTTA	Evans等^[17]
片段12	Ve16S	16S	335	Rev-AATCGTTGAACAAACGAACC	Evans等^[17]
片段13	Vert-16S-eDNA	16S	284	For-AGACGAGAAGACCCYDTGGAGCTT	Vences等^[5]
片段13	Vert-16S-eDNA	16S	284	Rev-GATCCAACATCGAGGTCGTAA	Vences等^[5]
片段14	Fish16S F/D-2R	16S	222	For-GACCCTATGGAGCTTTAGAC	DiBattista等^[23]
片段14	Fish16S F/D-2R	16S	222	Rev-CGCTGTTATCCCTADRGTAACT	DiBattista等^[23]
片段15	PS1	COI	199	For-ACCTGCCTGCCGTATTTGGYGCYTGRGCCGGRATAGT	Balasingham等^[6]
片段15	PS1	COI	199	Rev-ACGCCACCGAGCCARAARCTYATRTTRTTYATTCG	Balasingham等^[6]
注：片段长度以花鳗鲡（Anguilla marmorata）（GenBank登录号：NC006540）为参考。

下载: 导出CSV

表 2 metabarcoding片段对各类群的识别率

Tab. 2 The species resolution rates of metabarcoding fragments for different groups

物种类群	物种数	识别率/%
物种类群	物种数	片段 01	片段 02	片段 03	片段 04	片段 05	片段 06	片段 07	片段 08	片段 09	片段 10	片段 11	片段12	片段 13	片段 14	片段15	全部片段
鱊属 Acheilognathus	14	71.43	71.43	85.71	71.43	85.71	85.71	71.43	57.14	71.43	85.71	57.14	85.71	85.71	71.43	100	77.14±11.83
鲟属 Acipenser	15	13.33	20	26.67	13.33	46.67	33.33	40	6.67	33.33	0	33.33	6.67	13.33	6.67	73.33	24.44± 19.46
光唇鱼属 Acrossocheilus	16	56.25	50	62.5	50	87.5	75	62.5	6.25	68.75	43.75	50	62.5	50	43.75	100	57.92± 21.45
黑头鱼属 Alepocephalus	7	57.14	14.29	28.57	0	28.57	28.57	28.57	0	28.57	14.29	14.29	28.57	42.86	28.57	100	29.52± 24.43
西鲱属 Alosa	6	50	0	0	0	33.33	33.33	33.33	0	0	66.67	16.67	16.67	16.67	16.67	100	25.56± 28.78
双锯鱼属 Amphiprion	10	70	50	60	20	60	50	30	10	60	60	30	80	60	50	40	48.67± 19.22
鳗鲡属 Anguilla	17	29.41	35.29	64.71	0	76.47	70.59	76.47	41.18	64.71	47.06	35.29	76.47	64.71	52.94	88.24	54.9± 23.5
细鲫属 Aphyocypris	5	60	60	60	20	60	40	40	0	100	80	60	100	80	100	100	64± 30.43
旗鳉属 Aphyosemion	8	100	100	100	100	100	100	100	12.5	100	100	100	100	100	100	100	94.17± 22.59
深海鳐属 Bathyraja	8	12.5	25	50	12.5	25	12.5	12.5	0	12.5	12.5	0	25	25	37.5	50	20.83±15.43
鲫属 Carassius	5	0	0	0	0	60	20	20	0	60	40	0	40	0	0	60	20±25.07
真鲨属 Carcharhinus	14	64.29	71.43	100	78.57	100	100	100	50	85.71	35.71	14.29	50	64.29	71.43	100	72.38±26.57
刺尻鱼属 Centropyge	16	75	75	87.5	31.25	62.5	56.25	62.5	12.5	87.5	37.5	0	50	50	62.5	100	56.67±27.8
九棘鲈属 Cephalopholis	5	100	100	100	40	100	60	60	0	40	40	60	100	100	100	100	73.33±32.66
鳢属 Channa	8	62.5	62.5	62.5	62.5	62.5	62.5	62.5	25	75	75	75	75	75	75	50	64.17±13.25
红鳍鲌属 Chanodichthys	5	20	20	20	20	60	60	0	0	100	60	0	100	60	60	100	45.33±36.62
燕鳐属 Cheilopogon	8	37.5	37.5	25	0	25	50	50	0	62.5	25	12.5	12.5	25	12.5	75	30±22.06
鲮属 Cirrhinus	5	60	60	60	40	100	100	100	0	60	60	20	60	40	40	60	57.33±28.15
胡子鲇属 Clarias	6	66.67	16.67	33.33	33.33	83.33	33.33	33.33	0	66.67	33.33	0	66.67	66.67	33.33	100	44.44±29.32
花鳅属 Cobitis	11	81.82	81.82	81.82	81.82	72.73	81.82	81.82	9.09	63.64	81.82	63.64	81.82	81.82	81.82	100	75.15±20.17
鲚属 Coilia	6	16.67	50	50	16.67	66.67	100	66.67	0	50	66.67	33.33	33.33	16.67	16.67	100	45.56±30.52
白鲑属 Coregonus	14	7.14	7.14	0	0	21.43	35.71	21.43	7.14	0	0	0	0	0	0	7.14	7.14±10.8
兵鲶属 Corydoras	6	33.33	33.33	33.33	16.67	50	50	66.67	16.67	50	50	16.67	0	0	0	66.67	32.22±23.12
杜父鱼属 Cottus	12	50	75	100	50	83.33	50	41.67	8.33	83.33	83.33	41.67	83.33	66.67	100	100	67.78±26.7
舌鳎属 Cynoglossus	12	100	100	100	83.33	100	100	100	83.33	100	83.33	83.33	83.33	83.33	83.33	100	92.22±8.61
长吻鳐属 Dipturus	8	0	12.5	37.5	12.5	50	37.5	37.5	0	50	12.5	0	25	25	25	62.5	25.83±19.75
盘属 Discogobio	5	40	40	60	40	60	40	60	0	20	40	0	40	60	40	60	40±20
塘鳢属 Eleotris	6	66.67	66.67	100	100	100	100	100	0	100	100	66.67	66.67	66.67	66.67	100	80±27.6
Enteromius属	5	100	100	100	100	100	100	100	40	100	100	100	100	80	100	100	94.67±15.98
石斑鱼属 Epinephelus	22	90.91	90.91	90.91	50	90.91	90.91	90.91	36.36	72.73	90.91	22.73	81.82	68.18	63.64	90.91	74.85±22.52
扁鳉属 Epiplatys	7	71.43	71.43	71.43	71.43	71.43	71.43	71.43	0	71.43	71.43	71.43	71.43	71.43	71.43	71.43	66.67±18.44
镖鲈属 Etheostoma	10	100	100	100	100	100	100	100	40	100	70	90	100	100	100	100	93.33±16.76
底鳄鳉属 Fundulopanchax	5	100	100	100	100	100	100	100	20	100	100	60	100	100	100	100	92±22.42
底鳉属 Fundulus	5	60	60	60	20	100	100	100	0	100	100	100	100	100	100	100	80±32.95
南乳鱼属 Galaxias	5	100	100	100	60	100	100	100	0	100	60	60	100	100	100	100	85.33±28.75
墨头鱼属 Garra	10	100	100	100	50	100	100	100	30	80	100	50	100	100	100	100	87.33±23.74
骨尾鱼属 Gila	5	40	20	60	20	60	40	40	0	60	60	20	60	60	60	100	46.67±24.69
露齿鲨属 Glyphis	8	25	12.5	25	0	37.5	37.5	37.5	0	25	25	25	25	25	37.5	25	24.17±12.01
纹胸属 Glyptothorax	11	18.18	18.18	63.64	18.18	36.36	27.27	36.36	9.09	54.55	9.09	27.27	54.55	54.55	63.64	81.82	38.18±22.59
颌须属 Gnathopogon	6	33.33	66.67	100	100	100	66.67	66.67	0	100	50	16.67	33.33	33.33	33.33	66.67	57.78±32.65
鳅鮀属 Gobiobotia	6	100	50	100	50	100	100	100	0	66.67	66.67	16.67	83.33	83.33	66.67	100	72.22±31.91
裸鲤属 Gymnocypris	8	37.5	37.5	37.5	0	12.5	12.5	0	0	12.5	37.5	12.5	12.5	12.5	0	25	16.67±14.69
海猪鱼属 Halichoeres	5	100	100	100	100	100	100	100	20	100	100	100	100	100	100	100	94.67±20.66
鳠属 Hemibagrus	5	60	60	60	80	60	40	60	0	60	60	20	60	60	60	60	53.33±19.52
䱻属 Hemibarbus	7	42.86	42.86	42.86	28.57	42.86	28.57	14.29	0	42.86	28.57	14.29	42.86	42.86	42.86	100	37.14±22.13
海马属 Hippocampus	18	88.89	88.89	88.89	55.56	66.67	66.67	61.11	27.78	100	77.78	55.56	100	100	100	100	78.52±22.11
刺蝶鱼属 Holacanthus	5	40	40	60	40	60	60	60	0	40	20	20	40	40	40	60	41.33±17.67
下鱵属 Hyporhamphus	5	100	100	100	100	100	100	100	40	100	100	100	100	100	100	100	96±15.49
野鲮属 Labeo	20	70	70	80	35	90	80	75	10	60	30	25	55	40	50	80	56.67±24.1
兔头鲀属 Lagocephalus	8	100	100	100	75	100	100	100	25	100	75	75	62.5	75	62.5	100	83.33±21.99
太阳鱼属 Lepomis	6	100	100	100	50	100	100	100	0	100	100	66.67	100	100	100	100	87.78±28.5
薄鳅属 Leptobotia	6	50	16.67	66.67	50	66.67	50	33.33	0	33.33	16.67	50	33.33	33.33	33.33	50	38.89±18.55
䱀属 Liobagrus	9	66.67	77.78	77.78	11.11	77.78	55.56	66.67	0	33.33	33.33	33.33	100	100	100	77.78	60.74±31.95
笛鲷属 Lutjanus	13	84.62	84.62	100	100	100	100	84.62	23.08	100	69.23	30.77	100	100	100	100	85.13±25.42
狼绵鳚属 Lycodes	5	60	60	60	40	60	60	60	0	100	40	20	60	60	60	100	56±25.3
虹银汉鱼属 Melanotaenia	6	100	100	50	0	66.67	33.33	33.33	0	66.67	0	50	33.33	100	66.67	100	53.33±36.84
小鳔属 Microphysogobio	14	64.29	64.29	85.71	64.29	71.43	50	35.71	7.14	85.71	35.71	35.71	57.14	71.43	57.14	100	59.05±23.76
泥鳅属 Misgurnus	5	100	100	100	20	60	60	40	0	60	60	0	60	100	60	100	61.33±35.02
假鳃鳉属 Nothobranchius	8	100	100	100	100	100	100	100	12.5	100	100	100	100	100	100	100	94.17±22.59
美洲鱥属 Notropis	10	100	80	100	60	100	100	80	10	100	60	70	100	100	100	100	84±25.58
沙塘鳢属 Odontobutis	6	100	100	100	100	100	100	66.67	33.33	100	100	100	100	100	100	100	93.33±18.69
太平洋鲑属 Oncorhynchus	8	37.5	62.5	100	37.5	100	75	75	0	37.5	0	0	50	50	50	100	51.67±34.68
白甲鱼属 Onychostoma	11	81.82	63.64	45.45	18.18	63.64	63.64	54.55	27.27	54.55	72.73	0	81.82	54.55	45.45	63.64	52.73±22.85
马口鱼属 Opsariichthys	5	100	120	100	40	60	60	20	0	60	60	40	100	100	100	100	70.67±35.35
罗非鱼属 Oreochromis	9	66.67	66.67	66.67	33.33	44.44	11.11	44.44	0	55.56	22.22	22.22	44.44	33.33	22.22	77.78	40.74±22.88
青鳉属 Oryzias	10	80	80	80	80	80	80	80	80	100	100	100	100	100	100	100	89.33±10.33
纹唇鱼属 Osteochilus	5	80	100	100	80	100	100	100	0	100	100	80	100	100	100	100	89.33±26.04
鲳属 Pampus	7	28.57	57.14	28.57	14.29	28.57	28.57	28.57	0	28.57	28.57	14.29	14.29	14.29	14.29	42.86	24.76±13.73
属 Pareuchiloglanis	5	60	60	60	20	60	100	60	0	60	100	60	60	60	60	100	61.33±26.69
鱥属 Phoxinus	7	100	100	100	57.14	100	71.43	71.43	0	71.43	71.43	14.29	100	100	100	100	77.14±32.31
花鳉属 Poecilia	5	20	20	20	20	60	20	40	20	20	20	20	20	20	20	60	26.67±14.47
多鳍鱼属 Polypterus	10	100	100	100	100	100	100	100	40	100	80	60	100	100	100	100	92±18.21
原唇齿脂鲤属 Prochilodus	5	20	20	20	0	20	0	40	0	0	0	0	40	40	40	20	17.33±16.68
拟鲿属 Pseudobagrus	12	41.67	58.33	75	33.33	58.33	58.33	58.33	16.67	66.67	8.33	8.33	33.33	25	25	83.33	43.33±24.03
拟腹吸鳅属 Pseudogastromyzon	9	33.33	55.56	55.56	0	55.56	55.56	55.56	0	77.78	55.56	44.44	100	77.78	55.56	77.78	53.34±26.96
多刺鱼属 Pungitius	7	71.43	100	100	28.57	100	57.14	57.14	0	28.57	28.57	0	71.43	71.43	42.86	100	57.14±34.99
鳐属 Raja	6	16.67	16.67	33.33	16.67	66.67	16.67	16.67	0	50	66.67	16.67	66.67	66.67	66.67	66.67	38.89±25.72
波鱼属 Rasbora	7	100	100	100	71.43	100	71.43	100	14.29	71.43	100	57.14	71.43	71.43	71.43	100	80±24.03
吻鰕虎鱼属 Rhinogobius	5	100	100	100	100	100	100	100	0	100	100	100	100	100	100	100	93.33±25.82
鳑鲏属 Rhodeus	10	80	80	80	30	80	80	60	20	100	80	50	100	100	80	100	74.67±24.75
红点鲑属 Salvelinus	12	8.33	8.33	25	8.33	25	25	25	0	16.67	8.33	0	8.33	0	8.33	58.33	15±15.17
鳈属 Sarcocheilichthys	9	55.56	44.44	55.56	11.11	55.56	55.56	33.33	0	55.56	33.33	0	44.44	33.33	44.44	55.56	38.52±20.08
沙丁鱼属 Sardinella	6	100	100	100	66.67	100	50	66.67	33.33	50	66.67	66.67	100	100	66.67	100	77.78±23.29
蛇属 Saurogobio	7	100	100	100	71.43	100	71.43	71.43	0	71.43	71.43	14.29	100	100	100	100	78.1±31.88
南鳅属 Schistura	11	100	100	100	100	100	100	100	54.55	100	100	63.64	100	100	100	100	94.55±14.5
裸裂尻鱼属 Schizopygopsis	7	14.29	14.29	14.29	0	14.29	28.57	14.29	0	28.57	0	0	42.86	42.86	42.86	71.43	21.91±20.82
裂腹鱼属 Schizothorax	26	11.54	11.54	15.38	7.69	23.08	23.08	11.54	3.85	15.38	7.69	0	11.54	3.85	11.54	30.77	12.56±8.15
马鲛鱼属 Scomberomorus	6	66.67	66.67	66.67	33.33	66.67	66.67	66.67	33.33	100	66.67	100	66.67	66.67	66.67	100	68.89±19.79
平鲉属 Sebastes	19	31.58	31.58	31.58	10.53	36.84	47.37	26.32	0	10.53	15.79	0	26.32	15.79	15.79	68.42	24.56±18.09
瓢鳍虾虎鱼属 Sicyopterus	18	33.33	33.33	44.44	22.22	50	66.67	50	33.33	83.33	33.33	33.33	66.67	55.56	66.67	88.89	50.74±20.13
篮子鱼属 Siganus	6	33.33	33.33	33.33	16.67	100	66.67	66.67	0	66.67	16.67	0	16.67	0	0	33.33	32.22±30.52
鱚属 Sillago	6	100	100	100	100	100	100	100	50	100	100	66.67	100	100	100	100	94.44±15
鲶属 Silurus	6	66.67	66.67	66.67	50	100	100	66.67	16.67	66.67	50	50	66.67	66.67	100	100	68.89±23.46
鳜属 Siniperca	7	57.14	42.86	100	28.57	42.86	28.57	28.57	0	42.86	28.57	14.29	14.29	14.29	14.29	28.57	32.38±23.82
金线鲃属 Sinocyclocheilus	14	85.71	85.71	85.71	35.71	78.57	78.57	64.29	28.57	50	57.14	21.43	64.29	57.14	78.57	85.71	63.81±21.84
银属 Squalidus	7	42.86	28.57	57.14	42.86	100	57.14	57.14	14.29	57.14	28.57	14.29	42.86	28.57	28.57	100	46.67±26.16
扁鲨属 Squatina	5	60	60	60	20	60	60	60	0	20	40	40	60	60	60	100	50.67±23.74
疯鲿属 Tachysurus	8	75	75	75	25	100	100	100	0	75	37.5	0	75	50	50	75	60.83±33.03
东方鲀属 Takifugu	16	0	6.25	6.25	6.25	37.5	25	12.5	12.5	12.5	18.75	6.25	31.25	25	31.25	62.5	19.58±16.34
田中鳑鲏属 Tanakia	6	100	100	100	100	100	100	100	16.67	100	100	33.33	100	100	100	100	90±26.58
鲀属 Tetraodon	5	100	100	100	60	60	60	60	0	100	100	0	100	60	60	100	70.67±34.53
棱鳀属 Thryssa	6	100	100	100	66.67	100	100	100	0	100	100	100	100	100	100	100	91.11±26.63
金枪鱼属 Thunnus	9	0	11.11	11.11	0	0	0	0	0	0	0	0	0	0	0	0	1.48±3.91
茴鱼属 Thymallus	10	20	20	40	0	30	50	20	0	30	20	10	20	30	10	60	24±16.82
结鱼属 Tor	7	71.43	71.43	100	0	71.43	42.86	57.14	0	28.57	57.14	28.57	71.43	57.14	28.57	100	52.38±30.86
高原鳅属 Triplophysa	21	47.62	47.62	80.95	42.86	71.43	71.43	61.9	4.76	66.67	52.38	19.05	47.62	47.62	47.62	80.95	52.7±21.1
全部物种	935	62.3±31.84	61.95±32.54	68.9±30.81	42.92±33.46	71.62±27.3	64.34±29.13	60.89±29.48	11.79±17.79	64.6±30.65	54.58±32.28	36.71±32.27	64.2±31.8	61.6±32.35	59.71±32.78	82.76±24.66

下载: 导出CSV

参考文献(40)

[1]	姜维, 赵虎, 邓捷, 等. 环境DNA分析技术—一种水生生物调查新方法[J]. 水生态学杂志, 2016, 37(5): 1−7. Jiang Wei, Zhao Hu, Deng Jie, et al. Detection of aquatic species using environmental DNA[J]. Journal of Hydroecology, 2016, 37(5): 1−7.
[2]	陈炼, 吴琳, 刘燕, 等. 环境DNA metabarcoding及其在生态学研究中的应用[J]. 生态学报, 2016, 36(15): 4573−4582. Chen Lian, Wu Lin, Liu Yan, et al. Application of environmental DNA metabarcoding in ecology[J]. Acta Ecologica Sinica, 2016, 36(15): 4573−4582.
[3]	Zhang Shan, Zhao Jindong, Yao Meng. A comprehensive and comparative evaluation of primers for metabarcoding eDNA from fish[J]. Methods in Ecology and Evolution, 2020, 11(12): 1609−1625. doi: 10.1111/2041-210X.13485
[4]	Miya M, Sato Y, Fukunaga T, et al. MiFish, a set of universal PCR primers for metabarcoding environmental DNA from fishes: detection of more than 230 subtropical marine species[J]. Royal Society Open Science, 2015, 2(7): 150088. doi: 10.1098/rsos.150088
[5]	Vences M, Lyra M L, Perl R G B, et al. Freshwater vertebrate metabarcoding on Illumina platforms using double-indexed primers of the mitochondrial 16S rRNA gene[J]. Conservation Genetics Resources, 2016, 8(3): 323−327. doi: 10.1007/s12686-016-0550-y
[6]	Balasingham K D, Walter R P, Mandrak N E, et al. Environmental DNA detection of rare and invasive fish species in two Great Lakes tributaries[J]. Molecular Ecology, 2018, 27(1): 112−127. doi: 10.1111/mec.14395
[7]	陈治. 浙江近海鱼类多样性eDNA调查方法的建立与应用[D]. 青岛: 中国海洋大学, 2019. Chen Zhi. Establishment and application of eDNA method for fish diversity survey around Zhejiang coastal area[D]. Qingdao: Ocean University of China, 2019.
[8]	Valentini A, Taberlet P, Miaud C, et al. Next-generation monitoring of aquatic biodiversity using environmental DNA metabarcoding[J]. Molecular Ecology, 2016, 25(4): 929−942. doi: 10.1111/mec.13428
[9]	程馨雨, 陶捐, 武瑞东, 等. 淡水鱼类功能生态学研究进展[J]. 生态学报, 2019, 39(3): 810−822. Cheng Xinyu, Tao Juan, Wu Ruidong, et al. Functional ecology of freshwater fish: research progress and prospects[J]. Acta Ecologica Sinica, 2019, 39(3): 810−822.
[10]	李高俊, 顾党恩, 蔡杏伟, 等. 海南岛“两江一河”淡水土著鱼类的种类组成与分布现状[J]. 淡水渔业, 2020, 50(6): 15−22. doi: 10.3969/j.issn.1000-6907.2020.06.003 Li Gaojun, Gu Dang’en, Cai Xingwei, et al. The species composition and distribution of indigenous freshwater fishes of three main rivers in Hainan Island[J]. Freshwater Fisheries, 2020, 50(6): 15−22. doi: 10.3969/j.issn.1000-6907.2020.06.003
[11]	申志新, 李高俊, 蔡杏伟, 等. 海南省淡水野生鱼类多样性演变及保护建议[J]. 中国水产, 2018, 11(6): 56−60. Shen Zhixin, Li Gaojun, Cai Xingwei, et al. The evolution and protection of freshwater fish species in Hainan Province[J]. China Fisheries, 2018, 11(6): 56−60.
[12]	魏亚男, 王晓梅, 姚鹏程, 等. 比较不同DNA条形码对中国海岸带耐盐植物的识别率[J]. 生物多样性, 2017, 25(10): 1095−1104. doi: 10.17520/biods.2017164 Wei Ya’nan, Wang Xiaomei, Yao Pengcheng, et al. Comparison of species resolution rates of DNA barcoding for Chinese coastal halo-tolerant plants[J]. Biodiversity Science, 2017, 25(10): 1095−1104. doi: 10.17520/biods.2017164
[13]	Kruskal J B. Nonmetric multidimensional scaling: a numerical method[J]. Psychometrika, 1964, 29(2): 115−129. doi: 10.1007/BF02289694
[14]	Newmaster S G, Fazekas A J, Steeves R A D, et al. Testing candidate plant barcode regions in the Myristicaceae[J]. Molecular Ecology Resources, 2008, 8(3): 480−490. doi: 10.1111/j.1471-8286.2007.02002.x
[15]	魏观渊, 黄桂芳. 厦门湾春、秋季鱼类群落结构及其多样性[J]. 中国水产科学, 2021, 28(8): 1060−1068. Wei Guanyuan, Huang Guifang. Fish community structure and species diversity during spring and autumn in the Xiamen Bay[J]. Journal of Fishery Sciences of China, 2021, 28(8): 1060−1068.
[16]	Taberlet P, Bonin A, Zinger L, et al. Environmental DNA—for Biodiversity Research and Monitoring[M]. Oxford: Oxford University Press, 2018: 206.
[17]	Evans N T, Olds B P, Renshaw M A, et al. Quantification of mesocosm fish and amphibian species diversity via environmental DNA metabarcoding[J]. Molecular Ecology Resources, 2016, 16(1): 29−41. doi: 10.1111/1755-0998.12433
[18]	Riaz T, Shehzad W, Viari A, et al. ecoPrimers: inference of new DNA barcode markers from whole genome sequence analysis[J]. Nucleic Acids Research, 2011, 39(21): e145. doi: 10.1093/nar/gkr732
[19]	Bylemans J, Gleeson D M, Hardy C M, et al. Toward an ecoregion scale evaluation of eDNA metabarcoding primers: a case study for the freshwater fish biodiversity of the Murray-Darling Basin (Australia)[J]. Ecology and Evolution, 2018, 8(17): 8697−8712. doi: 10.1002/ece3.4387
[20]	Milan D T, Mendes I S, Damasceno J S, et al. New 12S metabarcoding primers for enhanced Neotropical freshwater fish biodiversity assessment[J]. Scientific Reports, 2020, 10(1): 17966. doi: 10.1038/s41598-020-74902-3
[21]	Kitano T, Umetsu K, Tian Wei, et al. Two universal primer sets for species identification among vertebrates[J]. International Journal of Legal Medicine, 2007, 121(5): 423−427. doi: 10.1007/s00414-006-0113-y
[22]	Shaw J L A, Clarke L J, Wedderburn S D, et al. Comparison of environmental DNA metabarcoding and conventional fish survey methods in a river system[J]. Biological Conservation, 2016, 197: 131−138. doi: 10.1016/j.biocon.2016.03.010
[23]	DiBattista J D, Coker D J, Sinclair-Taylor T H, et al. Assessing the utility of eDNA as a tool to survey reef-fish communities in the Red Sea[J]. Coral Reefs, 2017, 36(4): 1245−1252. doi: 10.1007/s00338-017-1618-1
[24]	Hebert P D N, Cywinska A, Ball S L, et al. Biological identifications through DNA barcodes[J]. Proceedings of the Royal Society of London, Series B: Biological Sciences, 2003, 270(1512): 313−321. doi: 10.1098/rspb.2002.2218
[25]	Collins R A, Bakker J, Wangensteen O S, et al. Non-specific amplification compromises environmental DNA metabarcoding with COI[J]. Methods in Ecology and Evolution, 2019, 10(11): 1985−2001. doi: 10.1111/2041-210X.13276
[26]	Menning D, Simmons T, Talbot S. Using redundant primer sets to detect multiple native Alaskan fish species from environmental DNA[J]. Conservation Genetics Resources, 2020, 12(1): 109−123. doi: 10.1007/s12686-018-1071-7
[27]	Jennings W B, Ruschi P A, Ferraro G, et al. Barcoding the Neotropical freshwater fish fauna using a new pair of universal COI primers with a discussion of primer dimers and M13 primer tails[J]. Genome, 2019, 62(2): 77−83. doi: 10.1139/gen-2018-0145
[28]	Sultana S, Ali M E, Hossain M A M, et al. Universal mini COI barcode for the identification of fish species in processed products[J]. Food Research International, 2018, 105: 19−28. doi: 10.1016/j.foodres.2017.10.065
[29]	Gantner S, Andersson A F, Alonso-Sáez L, et al. Novel primers for 16S rRNA-based archaeal community analyses in environmental samples[J]. Journal of Microbiological Methods, 2011, 84(1): 12−18. doi: 10.1016/j.mimet.2010.10.001
[30]	Freeland J R. The importance of molecular markers and primer design when characterizing biodiversity from environmental DNA[J]. Genome, 2017, 60(4): 358−374. doi: 10.1139/gen-2016-0100
[31]	Deiner K, Renshaw M A, Li Yiyuan, et al. Long-range PCR allows sequencing of mitochondrial genomes from environmental DNA[J]. Methods in Ecology and Evolution, 2017, 8(12): 1888−1898. doi: 10.1111/2041-210X.12836
[32]	李渊. 鲳属鱼类形态学和遗传学研究[D]. 青岛: 中国海洋大学, 2015. Li Yuan. Studies on morphology and genetics of Pampus species[D]. Qingdao: Ocean University of China, 2015.
[33]	Lockwood S F, Bickham J W. Genetic stock assessment of spawning Arctic cisco (Coregonus autumnalis) populations by flow cytometric determination of DNA content[J]. Cytometry, 1991, 12(3): 260−267. doi: 10.1002/cyto.990120309
[34]	Shubina E A, Ponomareva E V, Gritsenko O F. Genetic structure of the Salvelinus genus chars from reservoirs of the Kuril Islands[J]. Biochemistry (Moscow), 2007, 72(12): 1331−1348. doi: 10.1134/S0006297907120073
[35]	张辉. 西北太平洋两种卵胎生鱼类(许氏平鲉和褐菖鲉)的分子系统地理学研究[D]. 青岛: 中国海洋大学, 2013. Zhang Hui. Molecular phylogeography of two marine ovoviviparous fishes in Northwestern Pacific[D]. Qingdao: Ocean University of China, 2013.
[36]	潘晓哲. 中、日斑鰶耳石形态及线粒体基因组研究[D]. 青岛: 中国海洋大学, 2012. Pan Xiaozhe. Otolith morphological study and analysis of the mitochondrial genome of Chinese and Japanese dotted gizzard shad (Konosirus punctatus)[D]. Qingdao: Ocean University of China, 2012.
[37]	张艳春. 大口鳒Psettodes erumei线粒体全序列的研究和鲽形目鱼类系统进化分析[D]. 青岛: 中国海洋大学, 2009. Zhang Yanchun. Analysis of the mitoehondrial genome of Psettodes erumei and phylogenetic analyses of flatfishes[D]. Qingdao: Ocean University of China, 2009.
[38]	林小婉. 新疆额尔齐斯河北极茴鱼线粒体全基因组测定及遗传多样性分析[D]. 哈尔滨: 哈尔滨师范大学, 2021. Lin Xiaowan. Analysis of complete genome mitochondrial sequence and genetic diversity of Thymallus arcticus arcticus (pallas) in Xinjiang[D]. Harbin: Harbin Normal University, 2021.
[39]	Evans N T, Lamberti G A. Freshwater fisheries assessment using environmental DNA: a primer on the method, its potential, and shortcomings as a conservation tool[J]. Fisheries Research, 2018, 197: 60−66. doi: 10.1016/j.fishres.2017.09.013
[40]	Hänfling B, Handley L L, Read D S, et al. Environmental DNA metabarcoding of lake fish communities reflects long-term data from established survey methods[J]. Molecular Ecology, 2016, 25(13): 3101−3119. doi: 10.1111/mec.13660

施引文献

期刊类型引用(2)

1.	郭瑶杰，万武波，王海山，叶乐，陈治. 3组常用鱼类eDNA宏条形码通用引物对三亚水环境样品的物种检出效果比较. 南方水产科学. 2025(01): 66-76 . 百度学术
2.	金显仕，田洪林，单秀娟. 我国近海渔业资源研究历程及展望. 水产学报. 2023(11): 122-131 . 百度学术