整合DNA条形码与浮游生物调查解析东亚得里亚海鱼类浮游生物群落

《Fishes》：Integrating DNA Barcoding and Plankton Surveys to Resolve Ichthyoplankton Assemblages in the Eastern Adriatic

　　

鱼类浮游生物（鱼卵与仔稚鱼统称）为揭示鱼类产卵活动、补充潜力及塑造鱼类种群的环境驱动因子提供关键见解。尽管其生态重要性显著，亚得里亚海大尺度鱼类浮游生物研究仍较有限。研究人员于两个连续夏季（2019–2020年）在东亚得里亚海五个研究区（伊斯特拉半岛西岸(I)、克瓦内尔及克瓦内里奇(K)、中达尔马提亚北部(Z)、亚布卡海坑周边(J)及南达尔马提亚(D)）分析鱼类浮游生物群落，结合WP2浮游生物采样与DNA条形码以提升分类分辨率。共鉴定1382枚鱼卵与407尾仔鱼，隶属20科34种。群中以中上层分类单元占优势，鳀科（Engraulidae）占比超80%。空间格局显示克瓦内尔及克瓦内里奇、亚布卡海坑周边与南达尔马提亚丰度持续较高，这些区域以生产力增强、滞留环流或上升流为特征。α生物多样性指数显示年际相似性高，而ANOSIM与SIMPER分析确认区域间差异显著（p < 0.05），主要由鳀科驱动，次由鲆科（Bothidae）、雀鲷科（Pomacentridae）、鲭科（Scombridae）、鮨科（Serranidae）及鲷科（Sparidae）贡献。冗余分析（RDA）表明温度与叶绿素-a（Chl-a）是影响中上层与近岸类群的关键环境驱动因子。上述发现为东亚得里亚海鱼类浮游生物动态提供必要基线信息，并凸显与渔业管理相关的关键产卵与育幼区。本研究亦论证了整合DNA条形码与生态学分析以加强快速变化海洋环境中鱼类浮游生物监测的价值。

研究背景方面，鱼类早期生活史阶段（鱼卵与仔稚鱼，统称鱼类浮游生物ichthyoplankton）位于较低营养级，是幼鱼与成鱼饵料，其扩散、沉降与存活受环境强烈制约且适应力有限，自然死亡率高，因而成为环境变异的敏感指示物。准确鉴定鱼卵与仔稚鱼因个体微小、脆弱及形态与色素持续变化而困难，传统形态学鉴定耗时长且依赖分类专家。DNA条形码（以线粒体细胞色素c氧化酶亚基I(COI)基因为标记，比对GenBank、BOLD等参考库）可克服上述限制，不依赖形态即可将早期阶段归隶成体物种。过去五十年亚得里亚海鱼类浮游生物研究偏局部化或聚焦沙丁鱼、鳀等经济种，缺大尺度调查。本研究据此旨在：以DNA条形码鉴定鱼类浮游生物；描述两连续夏季群落组成；探析环境参数与群落关系，为东亚得里亚海（克罗地亚渔场）可持续管理提供基线。

主要关键技术方法方面，研究人员于2019年6月24日至7月13日与2020年7月8日至29日乘“BIOS DVA”号在东亚得里亚海五个分区共设118个站（2019年60站、2020年58站），用WP2网（网口面积0.255 m²、网目0.200 mm）自表层垂直拖至100 m或离底5 m，样品经200 μm筛洗后96%乙醇固定。环境数据取CTD（SBE 9）获表层温盐，卫星L3级Chl-a（OC-CCI V5.0，1 km日值）补生产力变量。实验室立体显微镜检计数后，除鳀科形态鉴别外，个体单管乙醇保存；鱼卵经裂解液软化绒毛膜、穿刺取细胞、蛋白酶K消化、酚氯仿抽提基因组DNA(gDNA)，PCR扩COI片段并单通Sanger测序，BLASTn比对NCBI限定种级≥98%相似阈值，未达标者排除。群落数据按Tanka公式标准化为个/m²，科水平聚合；用Shannon(H′)、Simpson(D)、Jaccard(J)指数评多样性，相对重要性指数(IRI)分主导/常见/稀有科；PRIMER v7建Bray–Curtis矩阵行双因素ANOSIM与SIMPER，CANOCO 5行RDA（含温、盐、Chl-a）。样本队列即上述118站两夏合并的1382卵+407仔鱼。

研究结果部分，于“3. Results”中，研究人员先报告两夏环境跨度：表层温22.0–26.6 ℃、盐36.7–38.9、Chl-a 0.08–0.66 mg m^?3。118样中105站见鱼类浮游生物，总计卵1382（2019年690、2020年692）、仔鱼404（2019年228、2020年176）；空间上克瓦内尔、亚布卡海坑周边、南达尔马提亚丰度突出。COI测序674条中499条（74.0%）≥98%定种，129条（19.1%）无显著命中或<90%排除，最终定20科34种（仅沙丁鱼靠形态，余33种COI条形码）。邻接树(p-距离，外类Scyliorhinus canicula)显示同物种单系簇(p<0.01)，Lepidorhombus boscii作鲽形目孤立支。α多样性两夏相似(J=80.65%)，H′微降0.847→0.805，D恒为0.275；鳀科占>80%，IRI>7845为绝对主导，鲆、雀鲷、鲭、鮨、鲷科为常见。ANOSIM全局区域间R=0.132(p=0.020)显著，夏间R=0.039(p=0.057)不显著；南达尔马提亚与余区差异最大（D vs K R=0.389、D vs I R=0.402）。SIMPER中鳀科贡献 dissimilarity 23.31–50.98%，克瓦内尔因鲆科升高而分离，伊斯特拉西岸因鮨、鲭、鲷科贡献发散最高。RDA前两轴约束特征值0.1695与0.0469，环境解释率24.4%；鲷科与鮨科关联高温，鲭、鳀、鲆科关联高Chl-a，伊斯特拉西岸受温度伪F=3.0(p=0.034)显著，克瓦内尔沿Chl-a高值但未达显著。

讨论部分浓缩：高丰度区（克瓦内尔、亚布卡、南达尔马提亚）跨年重现，契合已知滞留、上升流与淡水层化机制，并与历史产卵育幼区重叠。DNA条形码提升分类分辨率，但约1/5扩增子未定种源于模板降解或低质量读长而非COI分辨力限，建议建区域COI-形态联合参考库。所有定种者春夏季产卵，中上层主导反映高繁殖力广播产卵策略，鳀绝对优势与既往一致。生物多样性年际稳、区域间ANOSIM/SIMPER显结构分化，南达尔马提亚因航次末尾鳀产卵中心南移而富集，克瓦内尔因鲆科（比目鱼渔业北亚得里亚核心）加分。RDA仅释部分变异，环流、输送、产卵物候与饵料亦介入。结论部分译文：本研究为东亚得里亚海首个大尺度DNA支撑的两连夏鱼类浮游生物评估，凸显克瓦内尔、亚布卡海坑周边、南达尔马提亚等关键产卵育幼区；DNA条形码联合经典浮游调查定20科34种，未解扩增子多因降解非技术限；空间分化主由鳀科次由鲆、雀鲷、鲭、鮨、鲷科驱动，温与Chl-a为关键但非唯一驱动；该整合路径应纳入地中海长期监测，为渔业管理与保护规划供基线。论文发表于《Fishes》。