📄 文献精读 · 延伸·共产物与非光路线

光合细菌强化一体化水产养殖废水生物修复系统：氮磷功能基因与营养盐去除 | 彭博士工具屋

在沉淀—生物膜—滤食性贝类—大型海藻四单元一体化水产养殖废水生物修复系统（IAWBS）中投加 1‰ 光合细菌混合菌（R. palustris、R. capsulatus、R. rubrum，2:2:3），NO2−-N、NO3−-N、PO4³−-P、NH4+-N、CODMn 去除率分别达 79%、74%、83%、90%、71%。高通量 qPCR 检出 22 个氮基因与 9 个磷基因，nirK、nirS、ureC、nifH、amoB 与 phnK、phoD、phoX 表达显著提升，并与菌群结构显著相关。

光合细菌水产养殖废水功能基因脱氮除磷菌群网络高通量qPCR
💡 核心结论：在四单元 IAWBS 中投加 1‰ 光合细菌混合菌，可将养殖废水中 NH4+-N 去除率提升至 90%（NO2−-N 79%、NO3−-N 74%、PO4³−-P 83%、CODMn 71%），并显著上调 nirK、nirS、ureC、nifH、amoB 及 phnK、phoD、phoX 等氮磷功能基因表达，功能基因丰度与菌群分布紧密耦合。
📑 目录
文献信息要点速览背景与问题概念边界体系与方法结果逐条解读机制小结工程含义适用场景局限声明来源

文献信息

英文题名
Nitrogen and phosphorus-related functional genes enhance nutrient removal in the integrated aquaculture wastewater bioremediation system in the presence of photosynthetic bacteria
作者
Muhammad Naeem Ramzan, Ding Shen, Yingzhen Wei, Arslan Emmanuel, Regan Nicholaus, Wen Yang, Zhongming Zheng
期刊 / 年份
Research Square（预印本，网络首发） 2024
DOI
10.21203/rs.3.rs-5443318/v1
卷期 / 页码
预印本，共 23 页（posted 2024-12-14）
单位
宁波大学（Ningbo University）
一句话概括：以高通量 qPCR 功能基因芯片结合 16S rRNA 测序与网络分析，解析光合细菌投加如何通过改变细菌群落结构来增强 IAWBS 的氮磷循环与营养盐去除。

一图不看也要知道的几件事

  • 试验系统为四单元 IAWBS：沉淀、生物膜、滤食性贝类、大型海藻，各单元置于 70 L 白色聚乙烯桶中顺序串联。
  • 生物膜单元悬挂聚乙烯刷（直径 0.2 m、长 0.5 m），密度 4 个/桶；贝类单元投放泥蚶（Tegillarca granosa）70 个/m²；海藻单元投放龙须菜（Gracilaria lemaneiformis）0.2 g/L。
  • 废水取自宁波鄞州某高密度大棚对虾养殖场；每 2 天向下游单元输送 20 L 废水。
  • 光合细菌混合菌（3Rs）由 Rhodopseudomonas palustris、Rhodobacter capsulatus、Rhodospirillum rubrum 按 2:2:3 组成，30 ℃、3000–4000 lx 光下活化；处理组每 4 天投加 1‰（体积比），试验持续 1 个月。
  • 养分去除率：NO2−-N 79%、NO3−-N 74%、PO4³−-P 83%、NH4+-N 90%、CODMn 71%，处理组均显著高于对照组（P<0.05；文中 CODMn 处写作 P<0.5）。
  • NH4+-N 平均去除效率达 88%～90%，为最突出的脱氮表现。
  • 岸基测序：16S rRNA V3–V4 区，引物 338F/806R，Illumina 平台 250 bp 双端；UNOISE3（unoise_alpha=2、minsize=8）产生 ZOTU，比对 SILVA 数据库。
  • 功能基因芯片检测共识别 31 个基因：22 个与氮相关、9 个与磷相关，涵盖同化/异化硝酸盐还原、有机磷矿化、硝化、反硝化、厌氧氨氧化、固氮等通路。
  • 优势菌门为 Proteobacteria、Actinobacteria、Bacteroidetes，三者合计相对丰度>70%；优势纲为 Alphaproteobacteria、Bacteroidia、Gammaproteobacteria、Actinobacteria。
  • 最丰富的氮基因是 ureC（尿素水解/氨化），最丰富的磷基因是 phnK（膦酸转运与 ATP 结合）。
  • 投加 PSB 后显著增加的氮基因包括反硝化 nirK、nirS，氨化 ureC，固氮 nifH 与硝化 amoB；磷基因 phnK、phoD、phoX 表达也上升。
  • 贝类处理区 gdhA（将有机氮转为 NH4+-N）显著升高，海藻处理区 ureC 显著升高（P<0.05）；磷方面贝类区 ppx 增加、海藻处理区 phoD 显著增加（P<0.05）。
  • 网络分析显示功能基因与菌属的 Spearman 系数>0.8、q<0.05；处理组中与功能基因相连的细菌门类明显多于对照组。
  • 沉淀区所有功能基因显著高于系统其他区域（P<0.05）。
  • 机制判断：该系统中硝化与反硝化被增强，而固氮与厌氧氨氧化过程受到抑制。

研究背景与要解决的问题

集约化水产养殖向周边水体排放大量氮、磷，导致藻类过度增殖并破坏水生生态系统。既有修复手段包括商品微生物制剂、生物絮团、人工湿地、活性污泥与生物扰动生物等；贝类、生物膜与水生植物组合的生物滤器因生态友好而被越来越多地采用。

氮、磷相关功能基因在废水生物修复中起关键作用，例如 ppk1 与反硝化基因对氮磷去除至关重要，并关联到 Thauera、Dechloromonas 等聚磷菌与反硝化聚磷菌。将功能基因与生物修复技术结合，可提升氮磷去除效果。

光合细菌（PSB）是能进行光合作用的原核生物，广泛分布于湖泊、稻田、土壤、河流、活性污泥与海洋。已有研究用蓝细菌去除地下水硝酸盐；本试验采用的 PSB 为 Rhodopseudomonas palustris、Rhodobacter capsulatus 与 Rhodospirillum rubrum。

既有研究表明 IAWBS 类技术仍不成熟、运行与除污能力仍有提升空间；关于光合细菌存在条件下氮磷功能基因与浮游细菌群落之间的作用机制缺乏细节信息，本研究即针对这一空白。

概念与技术边界

IAWBS
Integrated Aquaculture Wastewater Bioremediation System，由沉淀、生物膜、滤食性贝类、大型海藻等单元串联组成的水产养殖废水一体化生物修复系统。
功能基因芯片（GeoChip/高通量 qPCR）
在 384 孔板上对大量目标基因同时进行定量 PCR，以 16S rRNA 基因作内参归一化，获得各功能基因的相对丰度。
ZOTU
零半径操作分类单元，用 UNOISE3 算法去噪、去嵌合后得到的序列单元，替代传统 OTU 以降低错误率。
3Rs
本研究中三种光合细菌（Rhodopseudomonas palustris、Rhodobacter capsulatus、Rhodospirillum rubrum）按 2:2:3 配成的混合菌群。
ureC / gdhA
ureC 编码参与氨化还原阶段的酶，将尿素水解为氨；gdhA 参与将有机氮转化为 NH4+-N。
nirK / nirS / nosZ
反硝化关键基因：nirK、nirS 编码亚硝酸盐还原酶（NO2−→NO），nosZ 编码 N2O 还原酶（N2O→N2）。
phnK / phoD / ppx / ppk
磷循环基因：phnK 属膦酸转运复合体，将惰性膦酸转为磷酸盐；phoD 参与有机磷转运；ppk 合成聚磷、ppx 分解聚磷。
网络分析
基于 Spearman 相关系数（>0.8、q<0.05）构建功能基因与细菌类群之间的共现网络，用 igraph/Gephi 可视化。

体系与方法

环节做法与参数
系统构建宁波大学梅山试验基地建四单元 IAWBS（沉淀、生物膜、贝类、海藻），每单元为 70 L 白色 PE 桶；生物膜单元聚乙烯刷 0.2 m×0.5 m、4 个/桶（一端固定桶底、另一端垂直悬于水体中部）；贝类泥蚶 70 个/m²；海藻龙须菜 0.2 g/L。
运行与投菌废水为高密度大棚对虾养殖废水，每 2 天将 20 L 水依次送入下游单元；3Rs 于 30 ℃、3000–4000 lx 活化，处理组每 4 天投加 1‰ 光合细菌，运行 1 个月。
水样采集与水质分析换水前后采集水与微生物样；水样经 0.45 µm 聚碳酸酯膜过滤测 NH4+-N、NO3−-N、NO2−-N、PO4³−-P，未过滤样测 TN、TP、CODMn，−20 ℃ 保存；用 Smart-Chem 200 自动批式化学分析仪，方法含次溴酸盐氧化、过硫酸钾氧化、镉柱还原、萘乙二胺分光光度、磷钼蓝分光光度与碱性高锰酸钾法。
DNA 提取与测序微生物样经 0.22 µm 膜真空过滤、−80 ℃ 保存；Minkgene 水样 DNA 试剂盒提取，Nanodrop One 定量；扩增 16S rRNA V3–V4（338F/806R），Illumina 平台 250 bp 双端测序；USEARCH v11.0.667_I8 + UNOISE3 去噪生成 ZOTU，比对 SILVA 分类。
功能基因定量提取 DNA 后于 384 孔板配制样品板与引物试剂板，用高通量自动微量取样器注入高通量 qPCR 芯片，SmartChip 实时 PCR 系统检测扩增与溶解曲线，Canco 程序读取检测状态与 Ct 值；以 16S rRNA 基因为内参归一化，Roche 仪器法获得 16S 绝对定量，其余基因换算得到。
统计分析养分含量用单因素方差分析（ANOVA）+Tukey 事后多重比较；基因数据用 SPSS、GraphPad Prism 9 与 R 4.0.4 处理；SR plot 软件以气泡图展示基因相对丰度变化；Kruskal-Wallis 秩检验区分组间差异；Hmisc 包计算功能基因与微生物丰度相关系数，igraph 构建网络、Gephi 0.10.1 绘图。

主要结果（原文图表逐条解读）

Fig. 1 系统各区域浮游细菌群落组成（门与纲水平）（原图未能稳定匹配，仅文字解读）
图示信息：按沉淀（SN）、生物膜（BF）、贝类对照/处理（CS/TS）、海藻对照/处理（CM/TM）展示门与纲水平相对丰度。
论文结果：优势菌门为 Proteobacteria、Actinobacteria、Bacteroidetes，合计>70%。沉淀区以 Proteobacteria 与 Bacteroidetes 最丰富；贝类区对照与处理组 Proteobacteria 均上升、Bacteroidetes 下降；贝类对照组 Proteobacteria 与 Actinobacteria 高于处理组，而处理组 Actinobacteria 增加。海藻对照到处理区 Proteobacteria、Bacteroidetes 均略降。纲水平上 Alphaproteobacteria、Bacteroidia、Gammaproteobacteria 与 Actinobacteria 最丰富：Alphaproteobacteria 的高值出现在贝类对照区，处理区下降，而 Gammaproteobacteria 在处理区上升，海藻处理区 Actinobacteria 增加。
研究意义：PSB 的投加改变了群落结构而非简单提高总丰度，其中 Alphaproteobacteria 与 Actinobacteria 的此消彼长提示氮循环功能菌群被重新组织。
Fig. 2 各区域氮、磷功能基因相对丰度气泡图：(a) 氮循环基因；(b) 磷循环基因（原图第 20 页第 2 张图 fig_020_02.png）
Fig. 2 各区域氮、磷功能基因相对丰度气泡图：(a) 氮循环基因；(b) 磷循环基因（原图第 20 页第 2 张图 fig_020_02.png）
图示信息：两幅气泡图，横轴为系统区域 SN、BF、CS、TS、CM、TM，纵轴为基因名；气泡大小与颜色共同编码同一数值，色标由黄（0.000）经绿、青蓝至紫（约 0.020）。(a) 22 个氮基因（amoA1、nirS2、nosZ1、nxrA、ureC、narG、nifH、nirK2、nirS1、hzo、hzsB、amoA2、nirS3、nosZ2、amoB、napA、nasA、nirK1、nirK3、gdhA、hao、hzsA）；(b) 9 个磷基因（cphy、bpp、phoX、ppk、pqqC、gcd、phnK、phoD、ppx）。
论文结果：(a) 中 ureC 最为突出，BF 与 CM 约 0.020，CS 与 TS 约 0.010，TM 中等；nirS2（CS）、nosZ1 与 nxrA（CM）、narG 与 nifH（SN）、nifH（TM）、amoA2（CS）、nirS3（BF）、nosZ2（CM）为次级信号，其余基因普遍低。(b) 中 phnK 在所有区域均偏高，BF 与 TM 约 0.020、CM 次之、SN 约 0.010、CS 与 TS 中等；此外 phoX 与 phoD 在 CM 有局部信号，ppx 在 BF、CM、TM 有分布，cphy、bpp、ppk、pqqC、gcd 在所有组均很低。
研究意义：氮循环以氨化基因 ureC 为主、磷循环以膦酸利用基因 phnK 为主，说明该系统的氮磷转化分别偏向氨化与有机膦酸活化路径；BF、CM、TM 区域功能基因丰度总体最高。
Fig. 3 功能基因与浮游细菌群落的共现网络分析（原图第 21 页第 3 张图 fig_021_03.png）
Fig. 3 功能基因与浮游细菌群落的共现网络分析（原图第 21 页第 3 张图 fig_021_03.png）
图示信息：四幅力导向网络图 (a)～(d)，节点为功能基因与细菌类群，连线表示显著相关（Spearman>0.8，q<0.05），绿色与红色边区分正/负相关；图例按过程给基因着色：反硝化（绿）、磷转化（青蓝）、硝化（棕/灰）、固氮（橙）、氨化（深绿）、有机氮矿化（粉），浅紫节点为细菌类群。
论文结果：处理组中与功能基因相连的细菌门类明显多于对照组，说明 PSB 投加重塑了菌群—功能耦合关系。具体关联包括：Mycobacterium 与 phnK 正相关，Rhodovulum 与 ppx、ureC 正相关，Vicingus 与 nosZ2 正相关；贝类处理区 Methylotenera 与 ureC、nosZ1 正相关，Rugeria 与 marivivens 与 phoD、phnK 正相关，Halomonas 与 nosZ1 正相关；海藻对照区 Nonlabens、Rhodovulum、Hyphomonas、Polaribacter 与 phoX、nifH、nirS2、nosZ2 正相关，gdhA 与 nirS1 正相关，沉积物与 nirK3、nirS、phnK、ureC、phoD 正相关；海藻处理区 mycobacterium、thalassococcus、Coxiella、OM43-clade 与 amoB、gdhA 正相关，Marinobacterium 与 nosZ2、ppx 正相关，Rugeria 与 Paracoccus 与 ureC、nirS2、phnK、nosZ1、nifH 正相关。
研究意义：功能基因与特定菌属的稳定共现说明营养盐去除由菌群协作承担，PSB 通过改变群落结构间接调控功能基因表达。
Fig. 4 光合细菌影响下 IAWBS 对 NO2−-N、NO3−-N、PO4³−-P、NH4+-N 与 CODMn 的去除率（原图第 22 页第 4 张图 fig_022_04.png）
Fig. 4 光合细菌影响下 IAWBS 对 NO2−-N、NO3−-N、PO4³−-P、NH4+-N 与 CODMn 的去除率（原图第 22 页第 4 张图 fig_022_04.png）
图示信息：分组柱状图，纵轴为去除率（%），横轴为 NO2−-N、NO3−-N、PO4³−-P、NH4+-N、CODMn；蓝色为对照组、橙色为处理组，带误差棒与显著性星号。
论文结果：处理组去除率分别为 NO2−-N 79%、NO3−-N 74%、PO4³−-P 83%、NH4+-N 90%、CODMn 71%，均高于对照组（对照组相应约 60%、67%、69%、71%、24%）；NO2−-N、PO4³−-P、NH4+-N、CODMn 的处理组提升达显著水平（P<0.05），CODMn 的组间差距最大。
研究意义：光合细菌对氨氮的强化作用最强（平均 88%～90%），说明其可作为高效好氧反硝化/同化型脱氮菌，同时显著改善有机物去除，但氮磷去除的具体提升幅度受单元类型影响。
Fig. 5 氮循环与磷循环功能基因的通路分布及丰度（原图第 22 页第 5 张图 fig_022_05.png）
Fig. 5 氮循环与磷循环功能基因的通路分布及丰度（原图第 22 页第 5 张图 fig_022_05.png）
图示信息：两幅“通路示意+基因丰度柱图”组合图：(a) 氮循环，箭头连接 Org-N、NH4+-N、NH2OH、NO2−-N、NO3−-N、NO、N2O、N2，标注 gdhA、ureC、amoA、napA/narG、nirS/nirK、nosZ、nifH 等基因并附各组丰度柱与显著性字母；(b) 磷循环，展示 Poly-P 与 PO4³−-P 的互转及 phnK、ppk、ppx、phoD 等基因。
论文结果：共识别 31 个功能基因（22 个氮相关、9 个磷相关）。沉淀区所有基因显著高于其他区域（P<0.05）。主要反硝化基因为 nor（norB/C）、nap（napA/B）、nir（nirK/S）。贝类处理区 gdhA 显著升高，海藻处理区 ureC 显著升高（P<0.05）。氮固定基因 nifH 存在于系统中。磷方面贝类区 ppx 数量增加、海藻处理区 phoD 显著增加（P<0.05）。聚磷形成与 ppk、ppk2、ppa 相关。
研究意义：投加 PSB 后硝化与反硝化被强化、固氮与厌氧氨氧化被抑制，氮代谢重心转向硝化—反硝化脱氮；磷则以聚磷累积与有机磷矿化为主，两条通路的基因证据互相印证。

机制小结

  • PSB 是高效好氧反硝化菌，可通过同化与好氧反硝化途径降解 NO3−-N；Alphaproteobacteria 在好氧无机氮转化中效率突出，可同时消耗 NH4+-N、NO3−-N 与 NO2−-N。
  • PSB 投加扩展了系统食物链、加快初级生产力，从而提高异养菌（如 Actinobacteria）种群数量，间接带动功能基因表达上升。
  • Alphaproteobacteria 中含有约 31%～35% 的硝化细菌，并携带 ureC 基因，因此 ureC 丰度在所有区域最高并主导氨化过程。
  • 反硝化与氨化基因增强、固氮与厌氧氨氧化受抑：nirK/nirS 编码的亚硝酸盐还原酶在 NO 充足时可替代 norS 将 NO 还原为 N2O，nosZ/nirK 比值升高有利于 N2O 进一步还原为 N2。
  • 磷代谢以聚磷途径为主：ppk 合成 poly-P、ppx 水解 poly-P、phoD 参与有机磷转运，海藻处理区 phoD 显著升高说明该区域主要靠 poly-P 累积除磷。
  • 功能基因与菌群的正/负相关网络反映互养与竞争关系，紧密耦合的群落结构提高了系统去除营养盐的稳定性与冗余度。

工程与产业含义

  • 投加方式可操作：1‰（体积比）PSB 混合菌每 4 天投加一次即可在 1 个月试验期内显著提升出水水质，具备工程化的低投加量特征。
  • 单元配置兼顾多营养级：生物膜（刷式载体）、滤食性贝类与大型海藻组合可分别承担硝化、滤食与养分吸收，光合细菌强化其中氮磷循环环节。
  • 调控重点在氮代谢：系统对 NH4+-N 的去除最突出（88%～90%），而 NO3−-N 去除相对偏低（74%），工程上需关注反硝化末端（nosZ 相关联）以进一步降总氮。
  • 沉淀区功能基因丰度整体最高（P<0.05），提示沉淀单元是氮磷功能微生物的重要富集位点，可作为后续强化改造的优先对象。
  • CODMn 的处理组提升幅度最大（约 71% 对 24%），说明该组合对有机污染负荷的削减同样有效，有利于出水稳定。

适用与不适用场景

✅ 适用场景
  • 高密度对虾等海水养殖尾水的一体化生物修复（沉淀—生物膜—贝类—海藻）场景。
  • 需要以低投加量生物强化方式提升氮磷与有机物去除的水产养殖水处理系统。
  • 对功能基因（nirK/nirS/nosZ/ureC/nifH/phnK/phoD 等）进行高通量 qPCR 定量与菌群耦合分析的研究场景。
  • 以泥蚶、龙须菜等滤食/大型海藻为主体的多营养级综合养殖水处理系统。
⚠️ 不适用 / 需谨慎
  • 淡水或非海水养殖体系：本系统面向海水养殖尾水，贝类与海藻种类具盐度特异性，直接外推需重新验证。
  • 要求明确区分自然背景与 PSB 独立贡献的场景：PSB 与贝类、海藻、生物膜的效应共同叠加，单一因素贡献未分离。
  • 需要长期稳定运行数据的工程决策：本研究周期仅 1 个月。
  • 厌氧氨氧化主导的脱氮路径：本系统中厌氧氨氧化被抑制，不适用于依赖 anammox 的工艺设计。

局限与不确定性

作者未解决 / 需要读者注意的问题
  • 本文为 Research Square 预印本（2024-12-14 网络首发），未经同行评审定稿，结论具有不确定性。
  • 试验仅运行 1 个月，未覆盖季节变化、温度波动与长期菌群演替。
  • 数据可用性声明为“应要求提供”，未见公开的数据与代码仓库。
  • 去除率显著性在原文中 CODMn 处写作 P<0.5（疑为笔误），部分表述与图表数据存在不一致。
  • 系统中存在明显的区域缩写与表述混乱（如 SN/BF/CS/TS/CM/TM 与 “SF/MA” 混用），可能影响区域对应关系的准确解读。
  • 未提供功能基因的绝对拷贝数（仅相对定量），也未验证基因表达与酶活、通量之间的定量关系。
  • 试验未设置仅投加 PSB 而不设贝类/海藻的对照，难以分离光合细菌的独立贡献。

数据 · 代码 · 经费 · 利益冲突

原文声明无竞争利益；资助来源为国家重点研发计划（2020YFD0900201）、宁波市公益技术应用研究项目（2022S164）与宁波大学王宽诚教育基金。数据可用性声明为“应要求提供（Data availability will be made on request）”。

来源与链接

🐾 返回主站