基于eDNA宏条形码与拖网调查的鱼类多样性结果比较——以浙江南部海域为例

A comparison of fish diversity results based on eDNA metabarcoding and trawl net: a case study of the southern coastal waters of Zhejiang Province

  • 摘要:
    目的 为了比较eDNA宏条形码与拖网调查在鱼类多样性上的相似性和差异性,评估eDNA技术能否取代拖网进行调查。
    方法 以浙江南部海域为例,采用生物信息学、传统形态鉴定、非度量多维尺度分析、差异显著性检验及线性回归分析等,比较了两种调查方法的鱼类组成、物种优势度、群落多样性指数等。
    结果 ①高通量测序结果中出现了不少(19.92%)非本海域鱼类OTUs(淡水鱼类、深海鱼类、热带珊瑚礁鱼类等);采样水层(表层、中层、底层)虽然不影响eDNA检出的鱼类OTUs个数(P≥0.073),但会影响鱼类OTUs组成(stress=0.138,R=0.117,P=0.001);提高采样密度(平行样个数、采样水层数)能够极显著地增加鱼类OTUs个数(P≤1.1×10−3);本海域有相当比例(24.402%)的鱼类OTUs处于低可注释状态,且部分鱼类序列存在种间相同或种内变异;② eDNA检出的鱼类物种数远高于拖网(152种 vs. 75种),二者存在显著的线性相关性(P=0.049 7,rPearson=0.381,R2=0.153);③以全部鱼类为分析对象时两种调查方法在物种优势度上达到显著差异(P=0.034);鳀、黄姑鱼等众多共有鱼类、优势种/极优势种在两种调查方法间的物种优势度差1~2个数量级;④ 12种高频率出现鱼类中,仅龙头鱼和中华栉孔虾虎鱼的序列丰度与尾数/重量间存在线性相关性,其余物种拟合模型中的rPearson较低甚至为负值;⑤两种调查方法测得的群落多样性指数差异较大,其线性相关性同样很低(P≥0.087)。
    结论 ① eDNA存在较大的不确定性;②除个别指标、物种或站位外,很难基于其中一种调查方法的结果去建模预估另一种调查方法的对应数值;③两种调查方法并不能相互取代,二者为互补、互校关系。本研究为认识eDNA技术的优缺点及合理使用不同调查方法提供了依据。

     

    Abstract: To compare the similarities and differences in fish diversity between eDNA metabarcoding and trawl net, this study, taking the southern coastal waters of Zhejiang as a case area, analyzed the fish composition, species dominance, and community diversity index of the two survey methods using bioinformatics, traditional morphological identification, non-metric multidimensional scaling analysis, significance test of differences and linear regression analysis. The results are as follows: ① A considerable number (19.92%) of non-indigenous fish OTUs (e.g., freshwater fish, deep-sea fish, tropical coral reef fish) were detected in the high-throughput sequencing results; although the sampling water layer (surface, mid-depth, bottom) did not affect the number of fish OTUs detected by eDNA (P≥0.073), it did influence the composition of fish OTUs (stress=0.138, R=0.117, P=0.001); increasing sampling density (number of replicates, number of water layers) significantly increased the number of fish OTUs (P≤1.1×10-3); a significant proportion (24.402%) of fish OTUs in this area were in a low annotation state, and some fish sequences exhibited interspecific identity or intraspecific variation. ② The number of fish species detected by eDNA was much higher than that by trawl net (152 species vs. 75 species), and there was a significant linear correlation between the two survey methods (P=0.049 7, rPearson=0.381, R2=0.153). ③ When all fish species were analyzed, the two survey methods showed a significant difference in species dominance (P=0.034); for many shared fish species and dominant/keystone dominant species, such as Engraulis japonicus and Nibea albiflora, the species dominance differed by 1-2 orders of magnitude between the two survey methods. ④ Among the 12 frequently occurring fish species, only Harpadon nehereus and Ctenotrypauchen chinensis showed a linear correlation between sequence abundance and individual counts/weight, while the rPearson values in the fitted models for the remaining species were low or even negative. ⑤ The community diversity index measured by the two survey methods differed considerably, and their linear correlation was also very low (P≥0.087). This study demonstrates that: ① eDNA has considerable uncertainty; ② Except for some indicators, species, or stations, it is difficult to model and predict the corresponding values of one survey method based on the results of the other; ③ The two survey methods cannot replace each other but rather have a complementary and cross-validating relationship. This study provides a basis for understanding the advantages and disadvantages of eDNA metabarcoding and for the rational use of different survey methods.

     

/

返回文章
返回