《Journal of Freshwater Ecology》：Contrasting summer and winter phytoplankton community states in the regulated Xiaolangdi Reservoir of the Yellow River revealed by shotgun metagenomics
编辑推荐：
淡水浮游植物群落在不同采样时期可能存在显著差异，但在调控水库中季节性与空间变异的相对重要性尚未完全明确。本研究采用鸟枪法宏基因组学（shotgun metagenomics）技术，对中国小浪底水库在两个对比性时间端点——冬季（2023年12月31日）与初夏（2
淡水浮游植物群落在不同采样时期可能存在显著差异，但在调控水库中季节性与空间变异的相对重要性尚未完全明确。本研究采用鸟枪法宏基因组学（shotgun metagenomics）技术，对中国小浪底水库在两个对比性时间端点——冬季（2023年12月31日）与初夏（2024年6月6日）的浮游植物群落进行了表征。研究设置了四个水库采样区域，每个区域和时期采集三个空间分离的野外重复样本，共获得24份样品。研究采用属水平α多样性、Bray–Curtis PERMANOVA与PCoA分析、CLR/Aitchison敏感性分析、多元离散度检验以及按区域分层指示物种分析。稀疏化属丰富度夏季高于冬季（61.74 ± 6.82 对 41.22 ± 8.88），而Shannon多样性（1.814 ± 0.093 对 2.402 ± 0.195）和Pielou均匀度（0.441 ± 0.028 对 0.650 ± 0.045）则夏季较低。Bray–Curtis PERMANOVA显示，采样季节/状态解释了属水平群落组成变异的91.6%（pseudo-F(1, 16) = 466.16，R2 = 0.916，P = 0.0001），而CLR/Aitchison分析给出了更为保守的估计值52.9%（R2 = 0.529，P = 0.0001）。区域效应和季节×区域交互效应远小于季节效应。夏季群落显示出较高的微囊藻属（Microcystis）、聚球藻属（Synechococcus）和蓝杆菌属（Cyanobium）相对序列代表度，而冬季群落在真核浮游植物谱系间分布更为均匀。高置信度的夏季指示属为蓝杆菌属（Cyanobium）、聚球藻属（Synechococcus）、微囊藻属（Microcystis）和角藻属（Ceratium）；冬季指示属为定鞭金藻属（Chrysochromulina）、隐藻属（Cryptomonas）、单针藻属（Mychonastes）、小环藻属（Discostella）、绿球藻属（Choricystis）和 Cyclotella。若干丰度较高的冬季归属类群，尤其是与 Teleaulax 相关的隐藻、海链藻属（Thalassiosira）和骨条藻属（Skeletonema），在物种水平解释之前需要参考基因组验证。
# 夏季与冬季浮游植物群落状态的显著分化：基于黄河小浪底调控水库的鸟枪法宏基因组学研究
## 一、研究背景与科学问题
浮游植物是内陆水体中最重要的初级生产者之一，在碳固定、养分循环、食物网动态和水质调节中发挥着核心作用。由于许多浮游植物类群世代时间短、对环境变化响应迅速，群落组成可在季节时间尺度上发生显著变化，因此季节演替被视作浮游生物生态学的核心组织过程，反映了光照、温度、混合、营养供给、牧食压力、越冬策略和物种互作等多重因素的协同作用。值得注意的是，即使总生物量变化不大，季节更替也可能深刻重组群落的物种组成和优势格局。
淡水生态系统中最具生态后果的季节状态之一是蓝藻优势。水华蓝藻兼具生理和生态上的优势特征，在温暖、富营养且水动力条件有利的环境中可获得竞争优势。其中，微囊藻属表现出显著的基因组与生态生理多样性。更广泛而言，蓝藻水华是生物性状、温度、营养可利用性、水动力条件及其他环境过程相互作用的结果，而非由单一普适驱动因子所决定。因此，夏季强烈的蓝藻优势可以与较高的分类丰富度共存，因为背景群落中仍可检测到大量低丰度类群。
水库为季节演替增添了空间维度，因为河流—水库梯度、水停留时间、大坝运行、悬移质泥沙和水位调控共同造就了异质性生境。小浪底水库是黄河干流上一座受水文调控强烈影响的大型调控水库。近期研究表明，该水库微型真核浮游生物在河流区、过渡区和湖泊区之间存在显著的空间分化。然而，调控水库中季节性与空间变异的相对重要性尚未被充分认识。与此同时，鸟枪法宏基因组学可在不依赖单一扩增子区域的情况下，在同一数据集中同时表征蓝藻和多种真核浮游植物谱系的分类组成，为传统显微镜检查和标记基因调查提供了互补视角。然而，基于读段的宏基因组图谱受到测序深度、数据库代表度、基因组特征和分类器不确定性的制约，且高通量测序数据本质上是成分性的，读段计数主要描述群落组分间的相对信息。上述特性要求采用保守的分类学解读，并进行显式处理成分结构的敏感性分析。
基于上述背景，研究人员利用鸟枪法宏基因组数据，表征了小浪底水库四个区域在两个对比性采样时间端点（2023年12月31日和2024年6月6日）的浮游植物群落，每个采样点设置三个空间分离的野外重复样本。研究提出三个核心问题：（1）属水平多样性结构在两个采样时期之间是否存在差异；（2）时间状态与空间区域中何者解释了更大比例的群落组成变异；（3）两个采样时期是否以不同的优势属和指示属为特征。研究人员假设暖季与冷季状态之间的对比将产生强于四个区域间空间差异的群落重构效应。
## 二、关键技术方法
研究人员在小浪底水库沿约15 km研究范围内设置四个主采样区域（A–D），每区域设置三个空间分离采样点（共12个点位），于冬季和夏季各采集一次表层水样（0.5 m深度），经0.22 μm滤膜过滤后提取环境DNA。采用CTAB法提取DNA，基于Illumina NovaSeq平台进行双端鸟枪法宏基因组测序（读长151 bp）。测序原始数据经fastp v0.23.0质控后，使用Kraken2（v2.0.8-beta）比对NCBI nt_20230606数据库进行分类学注释，并以Bracken v2.5进行物种丰度估计。统计分析包括：多元超几何稀疏化计算的α多样性指标；基于Bray–Curtis距离的PERMANOVA和PCoA；添加伪计数后的CLR变换及Aitchison距离敏感性分析；按位点分层的指示物种分析（IndVal）；以及排除三个存疑属后的分类学稳健性检验。所有主要生态学解读均基于属水平。
## 三、研究结果
### 3.1 测序质量与分类学注释
24份宏基因组样本共产生191,808,600条原始读段（28.96 Gb），质控后保留190,010,472条干净读段（28.53 Gb），保留率为99.06%。原始Q30值介于95.13%–96.15%之间，表明碱基判定质量在全数据集中保持高度一致。分类学注释在不同季节、不同高级分类阶元间的注释完整性存在差异，而属水平读段覆盖率在两个季节数据集中均较高，因此属被选为群落分析的主要分辨率。
### 3.2 α多样性的季节分化
夏季稀疏化属丰富度显著高于冬季（61.74 ± 6.82 对 41.22 ± 8.88，P = 0.0001），但Shannon多样性（1.814 ± 0.093 对 2.402 ± 0.195，P = 0.0001）、Pielou均匀度（0.441 ± 0.028 对 0.650 ± 0.045，P = 0.0001）和Hill q1多样性（6.16 ± 0.58 对 11.24 ± 2.23，P = 0.0001）均呈现相反模式，即夏季低于冬季。区域效应对丰富度不显著（P = 0.115），但对Shannon多样性和均匀度等指标显著；所有主α多样性指标均无显著的季节×区域交互作用。上述结果揭示了丰富度-均匀度解耦现象：夏季检测到的属更多，但序列丰度集中于少数优势类群；冬季检测到的属较少，但序列丰度分布更为均匀。
### 3.3 季节主导属水平群落组成
基于Bray–Curtis相异度的PCoA显示，夏季与冬季样本沿PCoA1（解释91.2%变异）近乎完全分离。因子设计PERMANOVA确认季节是群落组成变异的主导来源（pseudo-F = 466.16，R2 = 0.9155，P = 0.0001），区域仅解释3.09%的变异（P = 0.0102），季节×区域交互解释2.21%（P = 0.0260）。PERMDISP检验显示季节间多元离散度无显著差异。CLR/Aitchison敏感性分析中，季节效应仍高度显著（R2 = 0.5286，P = 0.0001），而区域和交互项均不再显著。排除Teleaulax、Thalassiosira和Skeletonema三个分类学存疑属后，季节效应依然高度显著（R2 = 0.8860，P = 0.0001），表明核心结论不依赖于这些不确定性属级信号。
### 3.4 优势属与指示属的季节更替
夏季群落的蓝藻序列代表度极高：微囊藻属、聚球藻属和蓝杆菌属合计约占夏季浮游植物归属序列丰度的79.6%，三者均为高置信度夏季指示属（IndVal分别为91.0%、96.2%和99.8%），角藻属（Ceratium）在全部12个夏季样本中均被检出且冬季未检出（IndVal = 100%）。冬季高置信度指示属包括定鞭金藻属（Chrysochromulina，冬季平均相对丰度9.20%）、隐藻属（Cryptomonas，7.57%）、单针藻属（Mychonastes，5.83%）、小环藻属（Discostella，3.11%）、绿球藻属（Choricystis，2.09%）和Cyclotella（1.22%），涵盖隐藻、定鞭金藻相关、绿藻和中心硅藻等多个谱系。此外，与Teleaulax相关的隐藻序列信号在冬季占浮游植物归属序列平均相对丰度的34.42%，但由于其物种级分配指向海洋栖息地且存在数据库异常报道，仅被保留为属级序列信号，不用于物种级生态学解读。
## 四、讨论与结论
### 讨论总结
本研究的核心发现是采样所得夏季与冬季状态之间的分化幅度和一致性。季节状态解释了Bray–Curtis属水平组成变异的91.6%，远超区域效应和交互效应，在更保守的CLR/Aitchison分析中仍解释了52.9%的变异，且排除三个存疑属后季节效应依然显著，表明这一结论不依赖于数据表示方式或个别不确定的分类学信号。夏季群落的低均匀度与微囊藻属、聚球藻属和蓝杆菌属的强优势相关，这些百分比描述的是相对序列代表度而非细胞丰度或碳生物量，并且宏基因组序列丰度受基因组大小、拷贝数、DNA提取效率和参考数据库表征等多种因素影响，因此不应等同于浮游植物生物量、细胞密度或水华强度。冬季群落重构涉及多个真核浮游植物谱系的同步富集，支持更均匀分布的冬季集合。然而，大型Teleaulax相关隐藻信号和部分海链藻属、骨条藻属的物种级分配偏向海洋或河口参考序列，需通过参考基因组竞争性比对加以验证，本研究仅将其视为数据库依赖的属级序列信号。空间异质性虽可检测到但远弱于季节效应，且在CLR/Aitchison分析中不再显著。研究存在多项局限性：每季节仅一次野外采样，构成端点比较而非完整年周期；未同步测量环境变量，限制机制推断；宏基因组读段比例提供相对分类信息而非绝对丰度；细粒度真核注释强依赖参考数据库覆盖度；样本在4 °C保存至过滤完成间隔不超过48小时；仅采集0.5 m表层水样，未解析垂直光带结构；未对全部干净读段进行互斥分类；未设置提取空白和文库空白对照。
### 研究结论翻译
鸟枪法宏基因组分析揭示了小浪底水库夏季与冬季采样所得浮游植物群落状态之间的显著分化。夏季表现出较高的稀疏化属丰富度，但均匀度和有效多样性较低，同时微囊藻属、聚球藻属和蓝杆菌属具有较高的相对序列代表度。冬季检测到的属较少，但序列丰度在多个真核浮游植物谱系间分布明显更均匀。采样所得季节状态解释的群落变异远大于区域或季节×区域交互效应，并在CLR/Aitchison敏感性分析后仍保持高度显著。由于Bray–Curtis和Aitchison效应量差异较大，两者作为同一强烈季节状态对比的互补表征予以报告。高置信度指示属分析进一步将蓝藻相关的夏季状态与以多种隐藻、定鞭金藻相关、绿藻和硅藻属为特征的冬季群落区分开来。由于环境变量未同步测量且每季节仅有一次野外采样，所观察对比的机制和轨迹尚无法直接解析。若干丰度较高的冬季宏基因组归属类群，尤其是Teleaulax相关隐藻、骨条藻属和海链藻属，在强调其物种级身份之前需进行未来的参考基因组验证。未来结合更高频率的时间采样、环境测量、绝对浮游植物丰度和基因组级分类学验证的研究，将有助于确定与夏季和冬季状态强烈对比相关的机制和时间动态。
