盐湖生态系统中碳氮磷循环功能基因对盐度的响应：基于运城盐湖宏基因组学的范围综述

《Microorganisms》：Spoilage-Associated Bacteria in Fresh Cow Cheeses: Diversity, Spoilage-Related Changes, Quantification and Identification Approaches—A Scoping Review

　　

盐度是盐湖生态系统中微生物群落结构与功能的关键驱动因子，然而其如何塑造参与碳（C）、氮（N）和磷（P）循环的functional genes仍知之甚少。研究人员沿运城盐湖自然盐度梯度采集宏基因组样品，考察盐度与微生物分类及功能多样性以及C、N、P循环功能基因的关联。分类与功能多样性指标均随盐度升高显著下降且彼此正相关。低、中、高盐度组间C、N、P循环功能基因组成与丰度差异显著：碳循环中多数碳固定基因在高盐度下丰度更高，而多数碳降解基因丰度降低，其中还原性三羧酸（rTCA）循环与Calvin循环基因丰度升高；氮循环中氮矿化与同化基因显著富集；磷循环中转运蛋白与嘧啶代谢基因丰度上升，嘌呤代谢基因相对贡献下降。共现网络分析显示C、N、P循环功能基因间存在致密正相关，mer、GLU与ppk1为高度连接基因。Mantel检验确定盐度与pH为功能基因变异的主要环境因子。结果表明盐度可能通过重塑盐湖生态系统微生物群落结构部分调控C、N、P循环功能基因。

研究背景与立项依据

内陆盐湖约占全球内陆水域面积的44%，具有高盐度、强蒸发及频繁干湿交替等独特水文地球化学特征，是全球元素循环中被忽视的动态界面。盐度作为最关键的enviro nmental filtering factor，通过渗透胁迫、改变养分有效性及影响代谢活性深刻调控微生物群落结构与功能性状。已有研究对盐度与碳降解、甲烷代谢、反硝化等功能基因丰度的关系存在矛盾结论，且除盐度外pH、电导率、Na+、K+、Cl?、CO32?、氧等因子在不同生态系统中的主导性尚存争议。C、N、P作为微生物生长代谢的核心生源要素，其生物地球化学循环在微生物控制下高度耦合，但高盐环境下沿盐度梯度的C、N、P代谢功能基因分布格局与耦合关系尚未被系统刻画。运城盐湖位于山西西南干旱-半干旱区，为典型硫酸钠型盐湖，受自然蒸发与人类采盐影响呈现显著盐度空间异质性，是探究盐度驱动微生物结构与功能的理想天然实验室。该研究由《Microorganisms》刊发，结合宏基因组与理化分析，旨在揭示盐度对盐湖微生物C、N、P循环功能基因的关联机制。

主要技术方法概括

研究人员于2024年10月在运城盐湖设18个采样点，按盐度分为低（<10%）、中（10%–20%）、高（≥20%）三组，每组6点各3个生物学重复共54份表层水体样，经0.22 μm过滤获微生物并冻存。水质理化指标（盐度、pH、DOC、TN、TP）现场或实验室测定。宏基因组DNA提取后Illumina HiSeq 4000双端测序，fastp质控，MEGAHIT组装，MetaGeneMark预测CDS，MMseqs2去冗得非冗余基因集；SOAPaligner比对待测读段计算基因丰度，DIAMOND比对NCBI-NR与KEGG注释分类与功能。统计上采用R语言行PCoA（Bray–Curtis）、ANOSIM、RDA、DESeq2（apeglm收缩算法算log2FC）、Wilcoxon秩和检验、Kruskal–Wallis检验、Mantel检验（ggcor包）及Gephi构建功能基因共现网络（|ρ|>0.7，FDR p<0.05，丰度>1%）。

研究结果

3.1 盐度对微生物多样性与功能组构成的影响

Shannon指数显示分类与功能多样性随盐度升高显著下降且彼此正相关。RDA表明盐度（r2=0.845分类/0.958功能）、pH、TN、TP与微生物类群及功能基因显著相关。PCoA证实三组间C、N、P循环基因谱差异显著（p<0.05）。

3.2 碳循环基因谱沿盐度梯度显著差异

碳固定以Calvin循环（24.77%–35.90%）、rTCA循环（22.47%–24.17%）、还原乙酰辅酶A途径（10.74%–17.02%）为主。Wilcoxon检验示pccA、korA/B、fhs、cooC等固定基因随盐度升显著富集，rbcS、coxL下降；淀粉/半纤维素/果胶/纤维素/壳聚糖降解基因（amyA、xylA、abfA等）下降，celF上升；甲烷代谢中mer、fwdB高盐富集，mttB降低。系统发育注释示低盐以Proteobacteria、Bacteroidetes、Actinobacteria为主，高盐Euryarchaeota（广古菌门）主导，Cyanobacteria上升。

3.3 氮循环基因谱沿盐度梯度显著差异

氮同化基因占36.62%–66.24%。随盐度升高gdhA（氮矿化）、GLU（同化）、nirA（同化型硝酸盐还原ANRA）丰度显著升高；napA下降，nirK、norB无统计显著但呈升势，nosZ显著升高。低盐水细菌（Proteobacteria为主）驱动，中盐过渡期Cyanobacteria与Balneolaeota上升，高盐极端嗜盐古菌主导。

3.4 磷循环基因谱沿盐度梯度显著差异

磷转运（25.15%–29.80%）与嘌呤代谢（26.56%–31.30%）占比最高。tmk、dut、cmk（嘧啶）、guaA（嘌呤）、phn/pst/ugp转运系、phoU、ppa、ppk1、ppsA随盐度升显著富集；phnN/X（膦酸盐代谢）下降。Proteobacteria、Actinobacteria、Bacteroidetes随盐度降，Euryarchaeota与Cyanobacteria升。

3.5 微生物功能基因与环境营养因子的关系

Mantel检验示盐度与pH显著影响C、N、P循环功能谱。VPA显示功能基因总变异54.19%被解释，盐度独立贡献7.28%，pH+养分独立14.23%；分类矩阵总解释63.81%，盐度独立9.13%。热图示盐度与rTCA/还原乙酰辅酶A途径（fhs、cooC、korA/B、icd、mdh）正相关；半纤维素降解基因与pH/DOC/TN正相关；nifD/H/K与TN/TP/DOC正相关，gltB/D与盐度/TP负相关；转运/嘧啶/嘌呤代谢与盐度及C/N/P正相关、与pH负相关。Euryarchaeota、Rhodothermaeota与盐度正相关。共现网络前30枢纽基因含11碳循环、2氮循环（GLU、gdhA）、17磷循环（ppk1、pstC、pyrF、guaA等），mer与porA/B、GLU、gdhA、guaA、ppk1、cmk显著正相关。

讨论部分总结

讨论指出高盐下碳固定增强而降解减弱，rTCA因低ATP成本（每固定碳2 ATP）贡献上升，产甲烷基因mttB受抑；虽古菌数量占优，但半纤维素降解仍由细菌特定类群执行，体现数量优势≠功能贡献。氮循环中矿化与同化增强、ANRA上升而DNRA下降，反硝化基因nosZ随盐度升利于N2O还原为N2，与低氧高盐环境适配。磷循环中转运与嘧啶代谢上升反映高盐渗透胁迫下核酸合成与Pi获取策略转变，嘌呤代谢相对占比降。共现网络以正相关为主，mer、ppk1、GLU为跨循环hub基因，印证C-N-P循环经微生物网络紧密耦合而非独立运作。

结论部分翻译

本研究以运城盐湖自然盐度梯度为体系，借宏基因组学揭示盐度与微生物群落结构及C、N、P循环功能基因的关联。分类与功能多样性随盐度显著下降且互正相关。群落由低盐Proteobacteria-Bacteroidetes-Actinobacteria主导向高盐Euryarchaeota绝对主导演变。碳固定中Calvin与rTCA循环基因富集；氮矿化与ANRA基因显著升高；磷循环转运与嘧啶代谢基因上升、嘌呤代谢相对贡献降。共现网络示C、N、P循环功能基因致密正相关，mer、GLU、ppk1为最高度连接hub基因。Mantel检验定盐度与pH为功能基因变异主因。研究在宏基因组层面系统阐明盐度与C、N、P循环功能基因及群落结构的关联，揭示元素间耦合关系，为理解盐湖生物地球化学循环响应盐化及预测内陆水体功能演化提供理论基础。