基于整合宏基因组学与物种群落层次模型推断污染物降解的微生物动态解析

《Journal of Hazardous Materials Advances》：Unravelling Microbial Dynamics of Inferred Pollutant Degradation through Integrated Metagenomic and Species Interaction Modelling Approaches

　　

河流微生物群落在污染物衰减过程中发挥着关键作用，然而沿河流纵向梯度调控这一过程的生态与功能机制仍知之甚少。本研究利用高分辨率全基因组鸟枪法（WGS）宏基因组学结合物种群落层次模型（Hierarchical Modelling of Species Communities, HMSC），对恒河及其源头支流中的微生物群落结构、组装过程与推断代谢潜力进行了系统研究。多域分析揭示了细菌、真菌、古菌和病毒具有显著的海拔关联模式，其中浊度（turbidity）成为主导的环境驱动因子，其次为电导率（electrical conductivity, EC）和pH。分类学分析显示，细菌（Rhodoferax、Curvibacter、Methyloversatilis、Caulobacter、Aeromicrobium、Azospirillum）、真菌（Fusarium）、古菌（Candidatus Nitrosotenuis）和病毒（Caudovirales）呈现清晰的海拔关联空间分布模式。上游群落富集与寒冷适应代谢功能相关的类群，中游区域由营养循环类群主导，下游区域则富集与应激适应及推定解毒途径相关的类群。功能分析识别出参与外源物质降解、芳香烃分解、能量代谢和细胞应激反应的酶编码基因家族和代谢途径的富集，且位点特异性分类群对污染物转化具有贡献。这些功能属性代表从宏基因组数据推导出的推断代谢潜力。HMSC进一步证实浊度为群落变异的首要决定因素，而EC、pH和海拔贡献次级效应。HMSC分析中真菌和古菌群落的负残差物种关联表明，生物相互作用和环境过滤可能共同影响群落结构。总体而言，本研究提供了一个整合宏基因组学与联合物种分布建模的框架，以联系微生物群落与环境梯度及推断的污染物转化潜力。这些发现为推断的自然净化过程提供了生态学见解，对河流修复、风险管理和生物修复具有指导意义。

本研究于2023年7-8月从恒河8个战略选择的位点采集水样，包括冰川源头上游（Badrinath的Alaknanda、Karnaprayag的Pindar、Karnaprayag的Alaknanda）、汇合前区域（Devprayag的Bhagirathi）、汇合后区域（Devprayag的恒河）、中游河段（Rishikesh的恒河）、城市影响区（Haridwar的恒河）及下游区域（Varanasi的恒河），按上、中、下游各3、3、2个位点分类。研究采用的主要关键技术方法包括：利用DNeasy Power Water Kit提取微生物基因组DNA，通过Illumina NovaSeq 6000平台进行全基因组鸟枪法测序（每样本10-15 GB）；使用Kraken2结合PlusPFP和NCBI nt数据库进行物种分类学分析，并以Centrifuger进行交叉验证；基于HUMAnN 3.0进行功能注释和代谢途径分析；运用MaAsLin2进行多变量统计以识别微生物特征与环境变量的显著关联；采用物种群落层次模型（HMSC）进行联合物种分布建模，以贝叶斯框架整合物种共现、空间结构和生物相互作用，通过方差分解、β系数估计和残差关联分析量化环境驱动因素的贡献，模型性能通过R2、均方根误差（RMSE）和潜在尺度减少因子（PSRF）评估。

**水体质量与测序统计**：所有样本均获得高质量测序 reads，GC含量约61-65%，97-99%有那么点就通过了质量过滤。浊度在Rishikesh城市影响位点达峰值（173 NTU），而Devprayag的Bhagirathi仅0.2 NTU；EC和TDS从上游Badrinath（47.6 μS/cm和32 mg/L）至下游Varanasi（327 μS/cm和171 mg/L）呈递增趋势，反映累积性人为输入。

**分类学组成与环境梯度多样性变化**：Alpha多样性分析显示Karnaprayag的Pindar具有最高的Observed OTUs、Chao1和ACE指数，表明极高的物种丰富度；Varanasi、Haridwar和Rishikesh也相对较高，而Badrinath和Devprayag的恒河则较低。主坐标分析（PCoA）显示Badrinath沿PC2与其他样本明显分离，表明独特的微生物组合；Devprayag的Bhagirathi与恒河聚集较近，Haridwar、Rishikesh和Varanasi形成邻近簇，反映环境相似位点的群落组成相关性。多域分类显示细菌占主导（72.38-92.58%），在Badrinath最高而Varanasi最低；真核生物在Pindar和Varanasi最丰富；古菌在Haridwar达峰值；病毒在Varanasi最高。门水平上，变形菌纲（Proteobacteria）在Devprayag的恒河最丰富（80.41%），放线菌纲（Actinobacteria）在Haridward和Rishikesh达峰值，疣微菌门（Verrucomicrobia）在Pindar显著富集（2.15%）。物种水平上，Rhodoferax sp. BAB1在Badrinath占36.14%，Caulobacter vibrioides在Devprayag的恒河占28.22%，Methyloversatilis sp. RAC08在Devprayag的Bhagirathi（4.92%）和Karnaprayag的Alaknanda（2.19%）丰富，Fusarium fujikuroi在Haridwar最丰富（2.1%）。古菌中Candidatus nitrosotenuis cloacae在Haridwar达峰值；病毒中Siphoviridae_sp.占17-30%，Varanasi特有未培养Caudovirales噬菌体高丰度（20%）。Centrifuger交叉验证与Kraken2结果高度一致，确认了主要微生物类群的可靠性。

**相关模式与群落变异环境驱动因素**：基于分类群的样本相关性分析显示Rishikesh-Haridwar、Rishikesh-Varanasi间存在显著正关联，而Pindar与Rishikesh、Haridwar、Varanasi关联很弱。冗余分析（RDA）解释了66.5%的约束变异，EC、浊度和pH为显著预测因子；经置换ANOVA和VIF分析，盐度因不显著被排除，TDS因与EC高度共线性（VIF>10）被排除。Rhodoferax sp. BAB1与Badrinath的pH和浊度正相关，Caulobacter vibrioides与Devprayag的浊度负相关。

**潜在危害相关代谢酶编码基因家族与途径分析**：MaAsLin2分析识别出46个显著差异丰富的酶编码基因家族（q<0.25），涉及氧化还原酶、转移酶、连接酶、水解酶和易位酶，以Badrinath、Karnaprayag的Alaknanda、Devprayag的Bhagirathi、Karnaprayag的Pindar和Haridwar最为突出。EC显示出最多的显著关联，表明离子浓度和溶解矿物负荷是微生物代谢功能的主要驱动因素。甲硫氨酸腺苷转移酶、多聚核糖核苷酸核苷酸转移酶、UDP-N-乙酰葡糖胺1-羧基乙烯基转移酶、磷-N-乙酰胞壁酰五肽转移酶、谷氨酰胺-果糖-6-磷酸转氨酶等与EC强烈相关。浊度与果糖-二磷酸酶、核酮糖-磷酸3-差向异构酶、组氨酸激酶、卟胆原合酶、UDP-葡萄糖4-差向异构酶等相关，涉及糖异生、固碳、信号转导、血红素生物合成和多糖代谢。pH相关酶包括色氨酸合成酶和苹果酸合成酶。Curvibacter贡献了大多数显著酶基因家族，表明其在介导营养循环、碳代谢和环境应激响应功能适应中的核心作用。

途径分析识别出155个显著差异丰富的推断代谢途径（q<0.25），主要涉及生物合成、前体代谢物与能量产生、降解/利用及超级途径。上游Badrinath显示出最高的推断代谢潜力途径丰度，广泛受EC、浊度和pH影响。Curvibacter贡献氨基酸（L-色氨酸、L-异亮氨酸、L-缬氨酸）、核苷酸、脂质（磷脂酰胆碱）和细胞壁组分（肽聚糖成熟）等生物合成途径；Methyloversatilis主要贡献磷酸泛酸生物合成、有氧呼吸和肽聚糖成熟途径；Azospirillum和Curvibacter与脂肪酸β-氧化、辛烷氧化、植醇降解等降解途径相关，在Badrinath富集。核心能量代谢途径包括戊糖磷酸途径、糖酵解、TCA循环、有氧呼吸、Rubisco支路、丙酮酸至异丁醇等。核苷酸途径确保DNA/RNA合成，氨基酸途径辅助蛋白质合成、抗菌效应和生物膜形成，脂质途径维持膜完整性以适应应激和污染物吸附，碳水化合物途径支持细胞壁合成和生物膜产生。基于Curvibacter sp. PAE-UM在Badrinath的最高途径丰度，与Rhodoferax sp. BAB1的基因组比较显示两者有~92%的序列同一性，表明可能为近缘物种。

**HMSC分析**：将Himalayan_River_Type作为随机效应以捕捉共享和未测量影响。模型后验预测准确性显示，物种特异性决定系数（R2）中位数为0.66，表明超过一半的属水平丰度变异被解释；RMSE中位数为0.36，表明对数转换丰度的预测误差低；Gelman-Rubin统计量（PSRF~1）证实模型收敛良好。固定效应系数（β）显示EC主要呈正关联，表明适度较高的离子强度支持大多数微生物类群生长；浊度呈均匀负效应，表明沉积物负荷和颗粒应激降低大多数类群丰度；pH主要呈负关联，表明群落对偏离最适pH范围敏感。方差分解显示，平均而言浊度解释35%的已解释变异，EC占31%，pH占23%，Himalayan_River_Type随机效应占11%。浊度对Aquabacterium、Aquincola、Methyloversatilis、Azospirillum、Bradyrhizobium、Brevundimonas、Caulobacter、Cupriavidus、Enhydrobacter、Herbaspirillum、Ideonella、Leptothrix、Piscinibacter、Pseudomonas、Rhizobacter、Rubrivivax、Rhodoferax、Sphaerotilus和Variovorax等属的变异解释度最高；EC对Halobacterium、Halobaculum、Halorientalis、Halolamina、Halosolutus、Pseudozyma和Salinigranum等嗜盐古菌的变异解释最大；Natronobiforma的变异主要受pH影响。Himalayan_River_Type对所有分类群影响可忽略，强化了物种水平模式主要由固定效应（水化学）而非随机效应（地理）解释的结论。

残差物种关联（Ω矩阵）显示网络以正关联为主，表明存在促进功能互补、交叉喂养和共营养降解的合作群落。负权重子网络揭示竞争排斥或捕食者-猎物样相互作用：Fusarium作为主导枢纽，与大量跨域属呈负关联，包括细菌（Caulobacter、Cupriavidus、Methyloversatilis）、古菌（Halorussus、Halorientalis、Halococcus）、真菌（Trichoderma、Schizophyllum、Aureobasidium、Cordyceps）和病毒（Alphafusellovirus、Lentivirus）；Candidatus Nitrosotenuis作为次级枢纽，与Pseudomonas、Ramlibacter、Rhizophagus、Ascochyta、Halapricum、Halobacterium和Valbvirus等呈负关联。Fusarium与Candidatus Nitrosotenuis彼此正关联，但环境偏好不同。这些模式证明残差生物相互作用，特别是涉及Fusarium和Candidatus Nitrosotenuis的负关联，在测量环境变量:变量效应之外显著贡献于群落结构。

**讨论部分总结**：研究提供了恒河微生物群落结构的多域评估，超越了早期孤立河段的研究。支流输入在种子微生物多样性和塑造生态过程中起关键作用。变形菌纲主导所有位点，与淡水宏基因组学研究一致，反映其在有机和无机污染物降解中的核心作用，在上中游更为富集；放线菌纲在下游丰富，支持营养循环和抗菌合成。Ascomycota参与植物材料降解和病原互作；Euryarchaeota介产甲烷作用和氮循环；Uroviricota中的Caudoviricetes通过裂解活动和辅助代谢基因调控微生物种群，显著影响生物地球化学过程。各位点的优势微生物特征具有明确的生态系统服务功能：上游富集抗菌产生者和生物降解者（Rhodoferax、Novosphingobium、Curvibacter、Opitutus、Phenylobacterium、Sphingorhabdus、Enhydrobacter、Gemmatirosa），贡献烃降解和生物活性化合物生物合成；中游富集营养循环者和耐应激微生物（Azospirillum、Methyloversatilis、Sulfuritalea、Caulobacter），贡献生物地球化学循环、宿主互作、固氮和PAH降解，以及套索肽、泛醌Q-10生产和生物肥料应用；下游 harbors 生物矿化者、生物膜形成者和环境解毒者（Aeromicrobium、Fusarium、Magnetospirillum、Ramlibacter），实现金属抗性、氮循环、铁螯合、生物膜介导的干旱抗性，以及霉菌毒素、赤霉素和抗菌物质生物合成。交叉验证分析确认了分类学分析的可靠性。

代谢专业化支撑推定污染物降解潜力方面，不同微生物类群与独特基因家族和代谢途径相关，受浊度、EC和TDS等环境因素影响，表明其在污染物转化和生态系统功能中的推断潜力。Curvibacter sp. PAE-UM与多个显著酶编码基因家族相关，包括与激素相关污染物微生物转化相关的甲硫氨酸腺苷转移酶、与RNA加工和应激适应相关的多聚核糖核苷酸核苷酸转移酶（PNPase）、与细胞壁生物合成和环境韧性相关的磷-N-乙酰胞壁酰五肽转移酶（MraY）、与环境应激下碳水化合物代谢和还原性生物降解相关的谷氨酰胺-果糖-6-磷酸转氨酶、与应激耐受和能量代谢相关的肽基脯氨酰异构酶（PPIase）和果糖-二磷酸酶（FBPase）、以及支持氨基酸生物合成和微生物生长的色氨酸合成酶等。Badrinath显示出相对较高的推断代谢潜力，Curvibacter、Methyloversatilis和Caulobacter贡献核心代谢功能，Azospirillum与降解过程相关。中心碳代谢途径的突出性支持活跃营养转化的潜力。

环境过滤作为主要组装驱动因素方面，浊度几乎对所有物种的变异解释度最高，作为该河流系统的主要生态过滤器；EC解释约31%变异，亦为最强预测因子之一；pH贡献有意义的分量；Himalayan_River_Type随机效应平均解释~10.8%变异，表明局部水化学对微生物群落结构的控制强于河流分区，但后者仍捕捉水文和地貌影响。Pindar河因高EC、浊度和TDS，以及高丰度古菌和病毒，成为限制性情境，其代谢途径在Badrinath完全缺失， abruptly restructuring常见优势微生物种群。两个主要生物调控guild——真菌guild（以Badrinath富集的Fusarium为代表）和古菌guild（以Haridwar富集的Candidatus Nitrosotenuis为代表）——通过负关联对更广泛微生物群落施加控制，表明微生物种群动态不仅受环境参数调控，还受生物相互作用塑造。

新颖稀有类群的生态学意义方面，先前中游河段真菌组研究未报道Fusarium物种，而下游Prayagraj近期工作记录了mass bathing事件中Fusarium fujikuroi的增加；本研究首次证明该物种在Badrinath的最丰富存在，表明其富集不限于下游人为影响区，而是源于Alaknanda支流。Rhodoferax sp. BAB1先前仅报道于自来水和紫坪铺水库，本研究首次在恒河上游发现，其与耐冷种Psychrotolerant species如R. ferrireducens T118、R. saidenbachiensis、R. sediminis和R. antarcticus的相似性突出其生物修复、反硝化和营养循环中的代谢多样性。

河流修复与生物技术应用方面，上游支流Alaknanda作为生物修复者、代谢物和富集途径的 reservoir；中游区域支持营养循环；下游群落增强推定解毒和应激耐受。这些动态对维持洪水或极端污染事件等干扰期间的恢复力至关重要。优势类群产生的生物活性化合物和抗菌代谢物为可持续生物修复、工业酶发现、农业投入品和治疗开发提供机会。

**研究结论翻译**：本研究整合理解了微生物群落结构、推断代谢潜力和环境过滤如何共同贡献于大型河流系统的自然恢复力。纵向功能分区表明，上游区域富集寒冷适应类群，中游区域为潜在营养相关群落，下游区域富集推断解毒者和应激适应类群。源头支流如Pindar、Alaknanda和Bhagirathi似乎是维持干流微生物多样性和功能恢复力的关键来源。通过耦合WGS宏基因组学与HMSC，真菌（Fusarium）和古菌（Candidatus Nitrosotenuis）被鉴定为潜在的生物调控者，表现出影响群落动态的负关联，超越了单纯的物理化学过滤效应。与有氧呼吸、还原、降解、生物合成和能量代谢相关的酶学能力富集，表明与微生物生存和潜在污染物衰减相关。值得注意的是，高分辨率分类学鉴定揭示了几个先前未在恒河系统报道的微生物类群。这些发现扩展了该生态系统已知微生物库，强调了源头区域作为新颖功能多样性 reservoir 的生态重要性。总体而言，本研究为理解恒河微生物群落组装和推断自然净化潜力背后的生态过程提供了见解。