西北阿拉伯湾近岸海域病毒样颗粒丰度的时间、空间与环境变异

《Viruses》：Temporal, Spatial, and Environmental Variation in Virus-like Particle Abundance in the Northwestern Arabian Gulf

　　

海洋病毒是微生物群落的重要组成部分并影响其结构与动态，但其在西北阿拉伯湾的变异性仍鲜有记录。研究人员在科威特三个代表不同环境背景的近岸站点对病毒样颗粒（VLP）丰度进行了为期四年的监测。各采样期采集表层与深层样品，并同步获取理化与营养盐数据。VLP丰度在研究期间波动明显，且站点间与层位间存在差异：南方站点及近表层VLP计数倾向更高，最低值记录于科威特湾。三站点的环境特征各异，但VLP丰度与单一理化变量的关系随站点与季节变化。当环境、空间与时间变量联合考虑时，随机森林（Random Forest）识别季节、溶解氧（DO）及营养盐（尤其硝酸盐NO3?）为VLP丰度的重要预测因子。广义加性模型（GAM）中溶解氧与VLP丰度呈非线性关联，在控制测得环境变量后深层丰度偏低仍显著，而站点效应在该模型中不显著。总体而言，科威特近岸水域VLP丰度呈现显著时空变异，其格局与多种环境及季节因子相关。该四年记录为西北阿拉伯湾病毒丰度提供基线，并支持进一步探究这一高变异性近岸系统中影响病毒动态的生物与环境过程。

研究背景方面，海洋病毒是海洋中最普遍且多样性最高的生物实体，其丰度通常比细菌高一个数量级，通过裂解感染塑造微生物群落组成，并借助“病毒分流（viral shunt）”调控生物地球化学循环。VLP计数是估算海洋病毒现存量的标准方法。西北阿拉伯湾为浅水半封闭系统，蒸发强、盐度高、与外海交换受限，叠加底质浅导致短距离内环境梯度强烈，且有Shatt Al-Arab河淡水与营养盐输入及北南浊度梯度，是检验环境变异与病毒丰度关系的理想区域。此前该区病毒研究多聚焦蓝藻噬菌体（cyanophage）g20基因多样性和时空变化，总VLP丰度的环境驱动机制尚不清楚。为此，研究人员利用四年现场数据集，考察科威特对比近岸水域VLP丰度的季节、空间与垂直变异，并借助多元与机器学习方法识别相关环境变量及其相对重要性。该研究发表于《Viruses》。

主要关键技术方法方面，研究人员建立三个月度双月采样站点队列：北部Al-Sabbiya（受Shatt Al-Arab悬浮沉积与淡水影响）、中部Al-Bidaa（城市近岸活动）、南部Al-Nuwaiseeb（水体较清、偏开阔）。2011年4月至2014年1月各站采表层与2 m深层，同步测温度、盐度、pH、溶解氧（DO）、浊度及NO₃^?、NO₂^?、NH₄⁺、PO₄^3?、SiO₃^2?；VLP经0.02 μm过滤甲醛固定后由表荧光显微镜（epifluorescence microscopy, EFM）SYBR Green I染色计数。统计上采用PCA（prcomp/envfit）、Spearman相关与BH-FDR校正、Kruskal–Wallis与Wilcoxon检验、随机森林回归（randomForest，含/不含营养盐两版，permutation importance与node purity）及mgcv包广义加性模型（GAM，薄板样条平滑，ML选模REML推断），并以Synechococcus sp. WH7803与本地蓝藻噬菌体KLS-01做悬浮沉积对感染影响的示意实验。

研究结果部分，3.1节时序与季节变异显示，2011–2014年三站表深层VLP均波动，SUB与NUW整体高于BID，均值约1×10⁸–1×10⁹ mL^?1；SUB与NUW夏季表层峰值高（如2012年6月SUB达9.93×10⁸ mL^?1），NUW 2013年10月深层达全研究最高9.94×10⁸ mL^?1，表层一般高于深层但随季节浮动。3.2节环境空间格局中，PCA前两轴解释67.4%变异，DO、温度、盐度、浊度向量显著；Spearman相关经FDR校正仅得6个显著项（如BID夏表VLP与浊度ρ=?0.82、与温度ρ=?0.70），Kruskal–Wallis显示NUW中位log₁₀VLP（8.63）显著高于BID（8.43，p<0.001），Wilcoxon表明表层（8.59）显著高于深层（8.47，r=0.196）。3.3节随机森林全模型R²=0.243，permutation重要性排序为季节>硝酸盐>溶解氧，去营养盐后R²降至0.100，证明营养盐变量显著提升预测。3.4节GAM最终模型（环境+站点+层位+季节）AIC最低，解释偏差25.03%（调整R²=0.201），深层效应估计?0.201（p<0.001）显著，站点不显著，仅溶解氧平滑项非线性显著（EDF=2.463，p=0.017），去季节重拟合后DO仍显著且共曲线性略降。3.5节沉积实验中，正对照WH7803在第9日降至零，1×与2×沉积处理同正对照，1.5×（48.9 NTU）后期宿主细胞数偏高（第10日7.34×10⁶ mL^?1），说明悬浮颗粒可改变噬菌体-宿主接触与裂解时序，但非浓度线性。

讨论部分总结，VLP丰度与地中海等其他高产近岸相近，空间分布不简单对应浊度梯度——BID最低源于科威特湾受限交换与污水/工业输入致微生物特征差异，而非北部高浊度直接造成；表层高于深层但与前期报道深层蓝藻噬菌体多样性高不矛盾（对象分别为总量与分类多样性）。随机森林中硝酸盐高排序反映预测重要性而非直接作用，营养盐经宿主生长间接影响病毒；GAM中DO非线性显著且与季节项不冗余，可能经由Synechococcus等宿主与氧耦合过程传递，但非低氧机制可解。沉积实验示意在科威特相关浊度下颗粒可暂吸附噬菌体延缓裂解，但未阐明现场浊度与总VLP关系机制。局限性包括EFM不能区分宿主与感染性、高颗粒背景或致聚集与遮掩、两月间隔漏掉短时过程、历史样无电镜/流式交叉验证。

结论部分翻译：据研究人员所知，本研究首次给出科威特近岸及西北阿拉伯湾VLP丰度四年评估。VLP丰度在站点、季节与层位间变异，表明病毒动态不能归于单一环境因子；随机森林与GAM显示季节、营养盐及溶解氧贡献了观测变异，GAM进一步揭示对溶解氧的非线性响应及深层较低丰度。结果为科威特水域海洋VLP丰度与环境变异提供基线，未来结合长期VLP监测、病毒宏基因组、微生物群落分析及病毒产生与感染直接测量，可厘清承载这些格局的宿主与过程。