《Microorganisms》：Metagenomic Insights into the Functional Profiles of Carbon, Nitrogen, and Phosphorus Cycles in Yuncheng Salt Lake Under Different Salinity Gradients
编辑推荐：
盐度是盐湖生态系统中微生物群落结构与功能的关键驱动因子，然而其如何塑造参与碳（C）、氮（N）和磷（P）循环的functional genes仍知之甚少。研究人员沿运城盐湖自然盐度梯度采集宏基因组样本，考察盐度与微生物分类及功能多样性以及C、N、P循环基因的关联
盐度是盐湖生态系统中微生物群落结构与功能的关键驱动因子，然而其如何塑造参与碳（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%，具有高盐度、强蒸发及频繁干湿交替等独特水文地球化学特征，是全球元素循环中被忽视的动态界面。盐度作为最关键的过滤因子，通过渗透胁迫、改变养分有效性及影响代谢活性深刻调节微生物群落结构与功能性状。已有研究对盐度如何影响碳降解与氮循环等功能基因存在矛盾结论，且除盐度外pH、电导率、阴阳离子等驱动因子在不同生态系统中争议未决。C、N、P作为受微生物控制的核心生源要素，其生物地球化学循环高度耦合，但高盐环境下沿盐度梯度的代谢功能基因分布模式与耦合关系尚未被系统刻画。运城盐湖位于山西西南干旱-半干旱区，为典型硫酸钠型盐湖，受自然蒸发与人工采盐共同影响，盐度空间异质性显著，是剖析盐度塑造微生物结构与功能的天然实验室。
论文发表于《Microorganisms》，研究人员整合宏基因组测序与理化分析，揭示盐度梯度下微生物演替、C/N/P循环功能基因响应、环境驱动因子及元素间耦合关系，深化了对极端高盐环境微生物驱元素循环的理解。
主要技术方法概括
研究人员于2024年10月在运城盐湖设18个采样点，按盐度<10%、10%≤盐度<20%、≥20%分低中高三组，每组6点各3个生物学重复共54份表层水（0–0.5 m）样本，经0.22 μm滤膜截留微生物后?80 ℃保存提取DNA。宏基因组采用Illumina HiSeq 4000双端测序，fastp质控、MEGAHIT拼接、MetaGeneMark预测CDS、MMseqs2去冗余获非冗余基因集，SOAPaligner比对计丰度，DIAMOND比对NCBI-NR与KEGG注释分类与功能。统计上以PCoA（Bray–Curtis）、ANOSIM检验组间差异，RDA评估环境因子解释率，DESeq2（apeglm）算log2倍变、Wilcoxon秩和（FDR<0.05）验基因丰度差异，Mantel（ggcor）析基因-环境关联，仅留|ρ|>0.7且FDR<0.05的强相关边用Gephi构建共现网络。
研究结果
3.1 盐度对微生物多样性与功能组构成的影响
Shannon指数显示分类与功能多样性随盐度上升显著下降且互相正相关。RDA中盐度（r2=0.845/0.958）、pH（r2=0.647/0.768）、总氮TN（r2=0.556/0.498）、总磷TP（r2=0.503/0.463）均显著关联微生物分类与功能基因。PCoA证实三组C/N/P循环基因谱显著分离。
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、uidA、vanA等）多显著减少，celF增多；甲烷代谢中mer、fwdB高盐富集，mttB减少。系统发育注释示低盐以变形菌、拟杆菌、放线菌为主，高盐耳形古菌（Euryarchaeota）速升并主导。
3.3 氮循环基因谱沿盐度梯度显著差异
氮同化基因占36.62–66.24%最高。gdhA（氮矿化）、GLU（同化）、nirA（同化型硝酸盐还原ANRA）随盐度显著富集；napA下降，nirK、norB上升趋势不显著，nosZ显著上升。低盐水细菌（变形菌为主）驱动，高盐极端嗜盐古菌主导氮循环。
3.4 磷循环基因谱沿盐度梯度显著差异
磷转运（25.15–29.80%）与嘌呤代谢（26.56–31.30%）占比最高。tmk、dut、cmk（嘧啶代谢）、guaA（嘌呤）、phn/pst/ugp家族（转运）、phoU、ppa、ppk1、ppsA随盐度显著上升；phnN/X（膦酸盐代谢）下降。古菌（Euryarchaeota）与蓝菌随盐度上升增多。
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与盐度正相关。共现网络Top30枢纽基因含11个碳循环（mer、porA/B、korA/B、fbp、mdh、mch、mttB、pccA、icd等）、2个氮循环（GLU、gdhA）、17个磷循环（ppk1、pstB/C/S、pyrF、purA/M/K、cmk、guaA、tmk、ppsA、ppa、phnC、rpiA、adk、dut），mer与porA/B、GLU、gdhA、guaA、ppk1、cmk显著正相关。
讨论部分总结
讨论指出高盐下碳固定增强而降解减弱，rTCA因每固定1碳仅耗2 ATP更适高盐低能环境，相关基因多由Halalkalicoccus、Haloquadratum、Halopenitus、Halobellus等古菌及Euhalothece蓝菌贡献；Calvin循环贡献者亦多为嗜盐古菌。碳降解回落并非古菌功能替代，而是放线菌/拟杆菌/变形菌特定群（Nitriliruptor、Pontimonas、Lewinella、Nioella、Salibaculum）受抑所致。氮侧gdhA、nirA、glnA、GLU升高反映有机氮矿化与硝酸源同化增强，nosZ随盐度上升利于N?O→N?减排，denitrification由古菌（Halorubrum、Halobellus、Halovenus）+变形菌（Roseovarius、Salinihabitans）+拟杆菌（Salinibacter）+Balneolaeota多类群共维；DNRA走弱而ANRA走强。磷侧转运（pst/phn）与嘧啶代谢古菌主导上升，嘌呤代谢相对占比下降但绝对贡献仍广，体现嗜盐古菌为主、蓝菌与变形菌补充的渗透应激资源重分配。共现网络以正相关为主，mer–porA/B–gdhA–GLU–ppk1–guaA–cmk等跨元素枢纽揭示C/N/P非独立而是通过微生物网络紧密耦合，复杂网络提升资源传递效率与群落抗扰韧性。
结论部分翻译
本研究以运城盐湖自然盐度梯度为体系，借宏基因组揭示盐度与微生物群落结构及C/N/P循环功能基因的关联。分类与功能多样性随盐度显著下降且互相正相关。群落由低盐变形菌-拟杆菌-放线菌主导演变为高盐耳形古菌（Euryarchaeota）绝对主导。碳固定中Calvin与rTCA循环基因偏高；氮矿化与ANRA基因显著富集；随盐度上升磷转运与嘧啶代谢基因增多、嘌呤代谢相对贡献下降，反映高盐胁迫下微生物对磷获取与核酸合成的自适应调整。共现网络示C/N/P循环功能基因间致密正相关，mer、GLU、ppk1为最高度连接枢纽。Mantel检验指盐度与pH为功能基因变异主要环境因子。研究在宏基因组层系统揭示盐度关联C/N/P循环功能基因与群落结构，阐明元素间耦合关系，为理解盐湖生物地球化学循环响应盐化及预测内陆水体功能演化提供理论基础。
