公民科学帆船队通过标准化eDNA宏条形码绘制西南太平洋开放海洋表层生物多样性图谱并连接卫星观测

《Molecular Ecology》：Citizen Sailors Provide a Detailed Map of Open-Ocean Surface Biodiversity and Help Link eDNA Data to Satellite Observations

　　

在盆地尺度表征开放海洋生物多样性仍具挑战，原因是采样稀疏以及验证卫星衍生的浮游生物群落组成估算的机会有限。通过 Citizens of the Sea 倡议，研究人员在2024 Pacific Rally 期间于新西兰、汤加和斐济之间环航的27艘帆船上部署了标准化环境DNA（eDNA）采样工具包，在1.5百万 km² 范围内产生728份表层海洋eDNA样本。研究人员使用三种互补的宏条形码标记（细菌16S rRNA、真核18S rRNA和动物细胞色素氧化酶I（Cytochrome Oxidase I, COI）基因），并遵循严格的生物信息质量控制，获得31,177个独特扩增子序列变异（Amplicon Sequence Variants, ASVs），提供了该研究不足区域首个高分辨率浮游生物生物多样性评估。在考虑人为来源变异后，所有宏条形码标记均揭示α多样性和群落组成的清晰纬度格局。ASV丰富度随海表温度（sea-surface temperature, SST）显著增加，群落周转与温度和叶绿素-a（chlorophyll-a）梯度强相关。eDNA推断的浮游植物功能群与卫星估算之间的比较显示，原核生物主导群落为中度一致（R² = 0.27），但真核浮游植物相关性弱（R² < 0.15），反映分子与遥感方法之间功能类型特异性不匹配。研究结果强化了公民科学船队能够在研究考察无法达到的尺度上生成高质量海洋学eDNA数据集，同时证明这些倡议在连接分子观测与卫星衍生评估以及改善对开放海洋浮游生物生物多样性理解方面的强大潜力。

海洋生物多样性对维持地球生命、气候调节、食物和氧气供给以及复杂生态过程至关重要。卫星平台显著增强了物理和生化海洋信息的全球获取，可连续监测海表温度、海洋颜色、浮游植物色素（例如叶绿素浓度）、海平面和洋流等，但生物组分数据在可比尺度上仍严重不足。多数海洋，尤其近海和偏远开放海洋，因传统研究考察成本高、后勤难而未被生物学充分探索。环境DNA（environmental DNA, eDNA）技术提供非侵入、成本有效且可扩展的替代方案，可覆盖多营养级，但包含当代和遗留DNA，可能模糊精细生态分辨率。已有全球海洋采样、Tara Oceans、Malaspina和Plankton Planet等努力，但单船科学巡航昂贵耗时，公民科学项目依赖简单可及的采样协议。Citizens of the Sea 倡议使用拖曳式eDNA采样器TorpeDNA，使巡航船只无需专用研究船或船上过滤即可采样。研究人员描述该倡议的首次大型现场活动：2024 Island Cruising Pacific Rally 中，27艘志愿帆船在新西兰、汤加和斐济之间收集表层水eDNA，覆盖超过1.5百万 km²、约20个纬度，使用细菌16S rRNA、真核18S rRNA和动物COI三种标记。研究目标为：评估分类多样性和群落组成的空间格局，并量化公民采集者和船相关采样效应以评估可靠性；检验eDNA群落格局与卫星衍生环境变量（包括SST和叶绿素-a浓度）的关系；利用多标记数据推断主要浮游植物功能群相对代表，并评估其与Copernicus Marine Service卫星衍生浮游植物功能型（phytoplankton functional types, PFTs）估计的对应关系。

关键技术方法概括：样本队列来自2024年5月6日至12月3日 Citizens of the Sea Pacific Rally，27艘帆船在新西兰、汤加和斐济之间航行，采集728份表层海水eDNA样本和86份现场阴性对照，覆盖约1,518,200 km²。船员经手册、视频和简报培训，使用TorpeDNA以5–12节拖曳约15 m、每次5 min、20 μm尼龙膜捕获，滤膜置PowerProtect保存液，ArcGIS Survey123记录元数据，冷藏运输至Cawthron Institute。实验室提取DNA，PCR扩增16S rRNA、18S rRNA和COI，Illumina NextSeq V2测序。生物信息使用CUTADAPT、DADA2、SILVA、PR2、MIDORI、GenBank、microDecon和MetaMate等。卫星数据包括MUR SST、NOAA VIIRS叶绿素-a和悬浮颗粒物（suspended particulate matter, SPM）、Copernicus GlobColour PFTs。统计使用广义线性混合效应模型（generalized linear mixed-effects model, GLMM）、主成分分析（principal components analysis, PCA）、冗余分析（redundancy analysis, RDA）和环境拟合。样本许可来自斐济、汤加和新西兰。

研究结果：
3.1 General Patterns：研究人员招募27艘船，采集831份样本，包括86份现场阴性对照，排除17份不合格样本；728份样本和86份对照处理并测序，97%样本至少一个引物扩增成功，16S为83%、18S为85%、COI为92%；获得约432 million reads和209,478个ASVs，过滤后约381 million reads（88.2%），各标记共保留31,177个ASVs。估计采样面积约1,518,200 km²。分类分配在高阶元较高，属和种水平较低。

3.2 Alpha Diversity：稀释曲线显示约413、502和487份样本分别达到16S、18S和COI的99% ASV丰富度。GLMM显示船ID解释大部分随机变异（16S为99.9%、18S为69.8%、COI为98.4%），测序批次效应较小。所有标记ASV丰富度与SST和叶绿素-a显著正相关，斜率范围为0.10–0.17；与船速正相关，16S和COI显著；与水深关系为COI正相关，16S和18S负相关。

3.3 Community Composition and Structure：16S最丰富门为Proteobacteria、Bacteroidota和Cyanobacteria，最多样门为Pseudomonata、Bacteroidota和Bdellovibrionata；18S最丰富为Arthropoda、未分类Eukaryota、Chordata和Cnidaria，最多样为Eukaryota、Arthropoda和Ciliophora；COI最丰富为Arthropoda、Eukaryota、Cnidaria和Mollusca，最多样为Arthropoda、Eukaryota、Bacillariophyta和Cnidaria。PCA显示所有数据集存在清晰生物地理格局，SST梯度沿第一或第二轴，未见按船ID聚类。环境拟合表明船速、SST、叶绿素-a和水深均显著相关，SST相关系数最高（0.53–0.73）。RDA显示，在去除测序批次和船ID效应后，环境变量解释的调整R²为16S 11.5%、18S 4.01%、COI 6.37%，SST解释最多（2.8%–9%），船速最少（0.1%–0.3%）。

3.4 Phytoplankton-Satellite Correlations：16S最丰富和多样的属包括Prochlorococcus、Trichodesmium和Synecococcus；18S中Dinophyceae占浮游植物reads和ASVs的95%以上；COI最丰富类群为Dinophyceae，其次为Bacillariophyceae和Bacillariophyta，多样性以Dinophyceae最高，haptophytes和green algae相对读丰度和ASV丰富度很低。GLMM显示，真核组最佳相关多用ASV丰富度指标，原核组最佳相关用相对读丰度。总叶绿素方面，18S指标呈负相关，COI指标呈正相关，相对读丰度R² = 0.018，三种ASV指标R²分别为0.068、0.093和0.102。六种PFTs中四种与至少一个eDNA指标显著正相关：diatoms、green algae、haptophytes和photosynthetic prokaryotes；dinoflagellates和Prochlorococcus仅呈负相关。diatoms最高R² = 0.03，green algae最高R² = 0.122，haptophytes最高R² = 0.04，photosynthetic prokaryotes最高R² = 0.272。除CHL、GREEN和PROKAR外，条件R²高于边际R²，表明技术因素解释较多方差。

讨论部分总结：研究人员认为，该数据集在6个月内获得700多份可用样本，超过Plankton Planet和Tara Expedition等既往努力，并在盆地尺度上提供迄今该区域最高分辨率eDNA基线之一。三种标记的ASV稀释曲线在500–600份样本之间趋于平台，说明采样策略和船队规模足以表征该区域表层可检测生物多样性。船ID和船速造成显著人为程序变异，船速正相关可能源于更高速度过滤更大水体积和更长采样距离；统计模型控制这些变异后仍检测到显著环境格局。ASV丰富度与SST和叶绿素-a正相关，符合浮游生物纬度梯度。Trichodesmium向热带相对丰度较高、Prochlorococcus向25°–30°S较高，与已知热耐受范围一致；但温度-光生态位模型预测Prochlorococcus向热带增加，而eDNA相对丰度向热带衰减，可能涉及营养共限制、Trichodesmium竞争、生态型分布差异或宏条形码相对丰度与细胞计数不等价。群落组成在所有标记中显示清晰纬度梯度，即使单独分析浮游植物和桡足类也可见，提示更高采样密度揭示连续生物地理周转而非离散群落。eDNA与卫星PFTs相关性总体弱（R² < 0.3），仅原核相对读丰度与卫星原核叶绿素、真核相对ASV丰富度与卫星总叶绿素、绿藻ASV丰富度与卫星绿藻叶绿素超过10%解释方差。原核对应较好，与受控中生态实验一致；COI通常优于18S，18S相关弱或负，可能因rRNA拷贝数变异、V4区域丰度表现差、甲藻色素表达和混合营养等。方法局限包括分子与卫星空间尺度差异、按米标准化ASV丰富度次优、卫星PFTs自身不确定性高，尤其绿藻、haptophytes和dinoflagellates；总叶绿素和原核不确定性较低。PACE等高光谱任务和算法发展有望改善整合。

研究结论翻译：海洋的物理和生化性质正以前所未有的速度变化，以人类刚刚开始理解的方式重塑海洋生态系统。然而，监测这些变化的能力仍然有限，尤其是在广阔开放海洋。正如医学依赖生理指标评估患者健康，生物多样性观测对追踪海洋生态状况至关重要。本研究表明，在受过培训的公民科学家支持下，使用标准化eDNA协议收集高分辨率、盆地尺度生物多样性数据如今可行。尽管研究人员发现eDNA与卫星衍生浮游植物功能型之间仅中等相关，但这既凸显当前卫星产品的局限，也显示分子工具地面验证遥感观测的前景。持续投资和优化后，公民驱动分子采样可帮助校准海洋颜色模型、改善生态系统预报并为保护政策提供信息。关键的是，该方法在空间和时间上均可扩展，提供通向全球分布、可重复海洋生物多样性监测的路径。持续长期努力对跟踪生态基线、检测状态转变和指导适应性管理至关重要，因海洋条件继续演变。