冷冻是否真正“冻结”了信息？生物量保存对污水处理系统微生物群落宏基因组谱的影响

《Ecological Informatics》：Does freezing actually “freeze” information over time? Biomass conservation impact on metagenomic profiling of microbial communities from wastewater treating systems

　　

冷冻是否真正“冻结”了随时间变化的信息？生物量保存对污水处理系统微生物群落宏基因组谱的影响

研究人员通过16S DNA扩增子测序，系统评估了保存温度（室温、4 °C、?20 °C、?80 °C）、保存时长以及是否添加甘油作为冷冻保护剂对活性污泥（activated sludge, AS）和好氧颗粒污泥（aerobic granular sludge, AGS）微生物群落的影响，旨在确定最有效的保存策略并识别宏基因组分析中潜在的偏差来源。研究结果显示，群落谱在多数保存时间和策略下基本得以保留；偏差主要涉及高丰度类群和革兰氏阴性菌；低丰度但具有相关性的类群可能在保存不当情况下丢失；冷藏适合短期保存，而冷冻更适合长期保存；甘油在?80 °C下有助于生物量保存，但在?20 °C下则无益。研究人员认为，类群初始丰度与细菌结构特征共同影响保存引入的偏差幅度，高丰度类群的损失在复杂微生物群落中影响更大。冷藏可作为短期保存的适宜方法，而?80 °C保存被证明是长期保存最有效的方法。甘油在?80 °C时略微改善保存效果，但在?20 °C时产生有害作用。

# 污水生物处理系统中生物量保存策略对微生物群落宏基因组谱的影响——论文解读

## 一、研究背景与目的

在生物反应器微生物群落研究中，采集生物量并进行宏基因组测序的理想目标是捕捉特定时间点微生物群落的“快照”。然而，从样品采集到测序之间可能出现各种偏差，尤其是当反应器与实验室距离较远或样品数量庞大时，临时保存不可避免。现有金标准方法涉及在约4 °C下运输样品，随后于?20 °C或?80 °C保存。尽管已有大量研究，但对于最有效的保存程序仍未达成共识。已有报道显示，粪便或土壤样品在不同温度下保存的结果相互矛盾，且偏差可能与分析方法或生物量类型有关。冷冻引起的偏差可能源于解冻过程中冰晶形成对细胞膜和细胞壁的物理破坏，导致基因组DNA和RNA在受控裂解前释放到细胞外环境，进而被胞外核酸酶降解。革兰氏阴性菌因细胞壁较薄被认为对这种物理破坏更敏感，因此其丰度在冷冻样品的DNA分析方法中可能被低估。添加甘油等冷冻保护剂被建议用于减少细胞损伤，但研究结果并不一致，且甘油在不同温度下的效果及对不同类群的影响尚不明确。

值得注意的是，现有长期研究极少涉及两种以上温度及冷冻保护剂添加与否的组合，且主要集中于动物相关组织或土壤，对污水处理系统中的活性污泥（activated sludge, AS）和好氧颗粒污泥（aerobic granular sludge, AGS）关注甚少。这些系统富含胞外聚合物（extracellular polymeric substances, EPS），其本身可能提供一定程度的冷冻保护。缺乏标准化保存方法可能导致许多关键革兰氏阴性菌（如Thauera、Nitrosomonas、Candidatus Accumulibacter）的相对丰度估计严重失真，从而影响对反应器性能的理解与优化。为此，研究人员通过16S DNA扩增子测序，系统测试了保存温度（室温、4 °C、?20 °C、?80 °C）、保存时长以及甘油作为冷冻保护剂对AS和AGS群落的影响，旨在评估最有效的保存策略及宏基因组分析中的潜在偏差来源。

## 二、主要技术方法

研究样品来自四个反应器：Consortium Cuoiodepur S.p.A.污水处理厂的好氧池（AS全尺寸反应器）、位于佛罗伦萨的中试规模AGS反应器，以及两个实验室规模的部分硝化AGS（PN-AGS）反应器。共20个样品、68个等分试样，涵盖不同保存条件与时长。总DNA使用DNeasy PowerSoil Pro试剂盒提取，通过六轮珠磨裂解细胞。16S rRNA基因V3-V4区域使用引物Pro341FB和805R扩增，产物在Illumina MiSeq平台进行双端测序（2×300 cycles）。生物信息学分析使用QIIME2 2024.1和DADA2进行过滤、合并和去嵌合体，分类分配基于SILVA数据库（Release 138）的V3-V4高变区。统计分析与可视化在R 4.3中完成，使用phyloseq、vegan、DESeq2、ALDEx2、ggplot2等包。β多样性采用基于Hellinger距离的主坐标分析（PCoA），PERMANOVA使用adonis2函数（9999次置换，按样品分块）。差异丰度分析采用DESeq2配对检验，并以ALDEx2作为组成性感知的补充方法。

## 三、主要研究结果

### 3.1 群落分析前的观察

阴性对照未检测到DNA，PCR扩增失败，确认试剂无菌。室温保存的AS样品（2022年7月21日采集）打开容器时产生可检测的硫磺气味。2023年8月1日采集的AS样品对应测序数据生态学上未饱和，被排除后续分析。

### 3.2 保存方法总体比较

基于“Main A”（AS）和“Main B”（PN-AGS）样品的PCoA显示，沿PC1轴的变异显著大于PC2轴。两种生物量类型尽管EPS组成不同，但总体趋势相似：群落分歧随时间增加，尤其在?20 °C保存条件下；?20 °C加甘油被证明是保持整体群落谱最无效的方法。4 °C保存在一周内能有效维持群落结构，但三个月后效率下降。?80 °C保存（无论是否加甘油）与直接提取样品最为接近。最丰富属（Main A中Defluvicoccus，Main B中Nitrosomonas）存在微小但可见的变异，在?20 °C加甘油条件下最丰富类群相对丰度明显降低。

### 3.3 生物量保存引起的α多样性变化

α多样性分析显示，观察到的属丰富度在直接提取与保存样品之间无系统性差异。然而，所有?20 °C加甘油保存的等分试样以及许多（但非全部）?80 °C无甘油保存的等分试样中，均匀度一致增加。Shannon指数随之变化，因其与均匀度数学相关。

### 3.4 ?20 °C加甘油保存的影响

AS样品在?20 °C加甘油保存后，PCoA显示样品主要按采集日期沿PC1轴排列，而保存引起的变异沿PC2轴一致分布。基于样品分块约束的PERMANOVA显示，直接处理样品与最长保存时间等分试样之间的群落剖面差异显著（p < 0.001）。差异幅度随保存时间一致增加。室温保存的等分试样与直接样品差异最大，但多数优势属仍被保留，仅Sulfurovum仅在室温保存样品中检出。即使?20 °C加甘油保存三年并经历多次冻融循环的等分试样，仍比其它样品更接近其直接处理对应样品。

### 3.5 ?80 °C无甘油保存的影响

AGS样品在?80 °C无甘油保存后，最丰富细菌基本保持，反应器生物量成熟过程中的群落转变可同时被直接处理与保存样品捕捉。然而，Propionivibrio相对丰度在所有样品中一致下降，与生物量成熟阶段无关。低丰度类群如Candidatus Brocadia在多数保存等分试样中近乎消失，仅在?80 °C加甘油保存一个月及4 °C保存一周的等分试样中可检测到。Candidatus Competibacter在?80 °C无甘油条件下反而增加，可能归因于随机或组成性偏差。少数类群（如Thermomonas、Collinsella）仅在保存样品中检出，但丰度极低（低于0.01%），无特定保存条件模式。

### 3.6 单个类群丰度的影响

差异丰度统计显示，大多数显著结果涉及相对丰度一致的下降。AS样品在?20 °C加甘油保存后，NS9 marine group（黄杆菌目）、Rhodocyclaceae和Thauera显著下降；AGS和/或Main B样品在?80 °C无甘油保存后，Arcobacter、Candidatus Kaiserbacteria、env.OPS 17、Leptospira、Nitrosomonas、NS11-12 marine group和Sphaerotilus下降。相反，Nitrosomonas和Obscuribacteriaceae在?20 °C保存中增加，Defluvicoccus和JG30-KF-CM45在?80 °C保存后增加。独立ALDEx2分析及六个月等分试样再分析确认了Candidatus Kaiserbacteria、env.OPS 17、Nitrosomonas、NS9 marine group、Rhodocyclaceae和Thauera的显著差异。下降趋势与类群初始丰度相关：Thauera在高丰度样品中下降更明显，Nitrosomonas在Main B中从约30%降至约15%，但在AS样品中（约0.5%丰度）反而略微增加。

### 3.7 甘油添加效应的评估

在?20 °C条件下比较有无甘油保存的等分试样，总体分类水平无实质差异，但Nitrosomonas相对丰度在加甘油样品中下降更明显。Xanthobacteriaceae未分类属和Sphingopyxis在加甘油样品中增加，Lysobacter下降。差异分析证实了Lysobacter和Nitrosomonadaceae属在加甘油样品中的显著下降，且差异倍变随类群初始丰度增大而增大。尽管如此，PN-AGS 1样品的采集时间顺序在加与不加甘油情况下均可正确识别。在?80 °C条件下，甘油添加与否差异很小，但Main B加甘油保存的等分试样中Candidatus Brocadia仍可检出。

## 四、讨论与结论

### 4.1 讨论总结

讨论部分指出，保存对整体群落剖面的影响极小，即使在次优条件下（如室温保存六个月），多数优势属及其相对比例仍得以维持，这可能与EPS基质中胞外DNA的长期持久性有关。然而，保存引入的偏差可能足以导致高度相似群落样品的错误识别，并叠加在生物变异之上，导致对反应器群落动态的误解。在类群水平上，显著差异几乎总是出现在各自群落中最优势的成员中，如AS中的Thauera和AGS中的Nitrosomonas。初始丰度可能比细胞壁结构更能影响保存后类群的检出模式。甘油在?20 °C下的有害效果可能归因于其渗透毒性超过保护作用，或甘油将冰点降至?40 °C导致?20 °C下细胞内相部分冻结。冷藏适合短至一周的保存，而?80 °C无论是否加甘油均是长期保存最可靠的方法。研究局限性包括不同比较条件下样品数量不均衡、甘油浓度单一以及无法直接测量物理细胞破裂。

### 4.2 研究结论翻译

保存对整体群落剖面引入的偏差可忽略不计，表明其不太可能显著影响生态评估，如时间序列分析。然而，单个类群丰度存在显著差异，强调在差异丰度分析中需要考虑保存效应。根据观察，研究人员假设丰度（统计学视角）和细菌结构特征（生物学视角）是影响保存引入偏差程度的重要协同因素，高丰度类群的损失在复杂微生物群落中可能产生更大影响。因此，需要针对保存生物量中的关键类群进行更精细的研究，以优化生物强化等过程。冷藏是短期保存的合适方法，而?80 °C保存被证明是长期保存最有效的方法。甘油在?80 °C时略微改善保存效果，但在?20 °C时产生有害作用。