高砷胁迫通过网络重构和潜在水平基因转移增强水稻土中毒力因子基因和病原体的表达

《Environment & Health》：High Arsenic Stress Enhances the Expression of Virulence Factor Genes and Pathogens in Paddy Soils via Network Restructuring and Potential Horizontal Gene Transfer

　　

土壤微生物组中的毒力因子基因（VFGs）构成潜在健康风险，然而环境胁迫如何驱动其富集仍知之甚少。研究人员整合了宏基因组学和宏转录组学分析以及时间序列微宇宙实验（6小时至60天），以研究水稻土中VFGs对砷（As）胁迫的响应。高砷胁迫（150 mg kg^?1）在DNA水平上降低了整体微生物α多样性，但矛盾地推动了潜在病原体的共现网络向以正关联（98.50%）为主导的高度互联架构转变，全局传递性增加至0.874。这种重构伴随着VFGs的显著富集，特别是涉及黏附、运动和生物膜形成的基因，其中As(III)被确定为VFG表达的主要驱动因子。VFGs表现出与移动遗传元件（MGEs）的强转录耦合，并且35个携带VFG-MGE对的高质量宏基因组组装基因组（MAGs）提供了水平转移潜力的直接基因组证据。值得注意的是，VFG上调反映了一种群落水平的胁迫耐受策略，而非直接向致病性的转变，但其意外后果是机会性病原体（包括ESKAPE病原体）的持续时间依赖性富集，到第60天有53个潜在病原体转录活跃。研究人员的发现证实，砷污染水稻土作为VFGs和机会性病原体的隐藏储存库，对食品安全和“同一个健康”（One Health）风险评估具有重要意义。

砷污染在水稻土中是一个全球性健康危机，由于水稻田独特的淹水条件，As(V)易被还原为毒性更强、迁移性更高的As(III)，增加其生物可利用性并在稻米中积累。土壤微生物组中的毒力因子基因（VFGs）广泛分布于细菌群落，在病原菌中协助人体入侵并构成健康风险。然而，环境胁迫如何驱动VFGs富集及其机制仍不清楚。现有研究多基于DNA水平的宏基因组分析，仅能反映基因丰度，无法捕捉代谢活性，且缺乏时间动态数据，导致无法区分VFGs增加是活跃威胁还是休眠遗传残留。因此，研究人员开展了一项研究，旨在阐明砷胁迫对水稻土中VFGs和潜在病原体富集、活性及时间动态的影响。

研究人员通过设置不同砷浓度（0、20、50、150 mg kg^?1）的微宇宙实验，在6小时至60天的时间序列中采样，整合宏基因组和宏转录组分析，揭示了高砷胁迫通过网络重构和潜在水平基因转移（HGT）增强VFGs表达和病原体富集的机制。该研究证实砷污染水稻土是VFGs和机会性病原体的隐藏储存库，对食品安全和“同一个健康”（One Health）风险评估具有重要意义。论文发表在《Environment》。

关键技术方法方面，研究人员从中国湖南怀化采集无砷污染水稻土（表层0–20 cm），设置四个砷浓度处理（CK、As20、As50、As150），每个处理三个生物重复，在8个时间点（H0、H6、H12、H18、D3、D7、D30、D60）采样，同时收集孔隙水用于As形态分析。提取土壤DNA和RNA，在Illumina NovaSeq平台进行宏基因组和宏转录组测序。使用Kraken2进行物种分类，DIAMOND比对VFDB和MGE数据库注释VFGs和MGEs。通过MEGAHIT组装、MetaWRAP分箱获取高质量宏基因组组装基因组（MAGs）。统计分析包括NMDS、中性群落模型（NCM）、零模型、共现网络、偏最小二乘结构方程模型（PLS-SEM）和贝叶斯模型。

研究结果如下：

**不同砷水平和时间序列下微生物多样性的变化**：通过分析Shannon指数发现，高砷胁迫（150 mg kg^?1）在DNA水平显著降低微生物α多样性（3.77%），而低、中砷处理无显著影响。随时间推移，Shannon指数先降后升。在RNA水平，中砷处理（50 mg kg^?1）增加多样性，高砷处理则显著降低；多样性随时间逐渐增加，第60天达峰值。NMDS分析显示DNA水平微生物群落结构在砷水平和培养时间间均显著差异，而RNA水平仅在时间间显著差异。

**微生物群落组成与共现模式**：通过比较不同砷处理下微生物类群相对丰度，发现高砷胁迫富集了Alphaproteobacteria，减少了Actinomycetes和Myxococcia。RNA水平活性群落中，Bacilli在早期占主导，后期Desulfuromonadia和Alphaproteobacteria增多。中性群落模型拟合度在高砷处理下降低，零模型显示标准化效应大小（SES）从2.33增至16.35，表明高砷胁迫加强了环境过滤。共现网络分析显示，高砷胁迫下网络复杂度增加，边数、平均度、密度均升高，正关联占主导（>99%）。

**砷胁迫对VFGs丰度的影响**：通过DNA和RNA水平定量分析，高砷胁迫显著增加VFGs相对丰度（DNA：3.01 vs 2.67 RPKG；RNA：146.22 vs 69.58 RPKM），尤其黏附、免疫调节（IM）、运动性、生物膜和抗菌活性/竞争优势（AA/CA）等类型。移动遗传元件（MGEs）丰度同样在高砷处理下显著升高。线性回归显示VFGs与MGEs丰度显著正相关，RNA水平关联更强（R²=0.77）。Spearman相关分析表明，随培养时间延长，更多VFG类型与MGE类型（特别是重组酶、整合酶、转座酶）显著正相关。

**砷胁迫和时间序列下水稻土中潜在病原体的富集**：通过DESeq2分析，高砷胁迫在DNA和RNA水平分别富集了8个和23个潜在病原体，包括大肠杆菌、铜绿假单胞菌、金黄色葡萄球菌等。其中12个病原体的转录丰度与VFGs总丰度显著正相关。长期胁迫下，第60天RNA水平有53个病原体转录活跃，包括ESKAPE病原体。共现网络分析显示，高砷胁迫下潜在病原体网络复杂度增加，全局传递性达0.874，正关联比例98.50%。

**砷胁迫下VFGs丰度增加的驱动因素**：通过PLS-SEM分析，MGE丰度和土壤性质直接正向促进VFG丰度，微生物α多样性负向影响，砷胁迫通过影响多样性和MGE间接促进VFG富集。贝叶斯模型表明，As(III)浓度是VFG表达的最强正向预测因子（后验概率Pr=99.9%），其次为Shannon指数、pH、As(V)，而土壤Eh负向贡献。

**携带VFG-MGE基因对的宏基因组组装基因组（MAGs）**：从宏基因组中恢复35个高质量MAGs，同时携带VFG和MGE基因对。这些MAGs主要属于Pseudomonadota、Bacteroidota和Myxococcota，其丰度随时间增加。转录活性分析显示，黏附和IM相关VFGs与转座酶相关MGEs表达最高，尤其Azospirillum和Magnetospirillum属。机会性病原体Flavobacterium和Enterobacter的MAGs在高砷胁迫下富集。

讨论部分总结：高砷胁迫重构了微生物群落网络，使其高度互联且正关联主导，这可能是一种群落级胁迫耐受策略。VFGs（如黏附、生物膜、运动性）的富集同时增强了胁迫耐受和致病潜力。MGEs与VFGs的强转录耦合表明HGT可能被加速。35个MAGs提供了直接基因组证据。虽然VFG上调主要反映生态适应，但长期后果是机会性病原体（包括ESKAPE）的时间依赖性富集，将砷污染水稻土转化为VFGs和病原体的活跃储存库。DNA与RNA水平的差异凸显了单独依赖DNA调查的局限性。研究结论翻译如下：

**结论**：研究人员的发现建立了一个机制框架，将砷诱导的微生物群落重构与水稻土中增高的致病性风险联系起来。研究人员证明，高砷胁迫驱动形成高度互联、合作性的微生物网络，以正关联为主导。在这些网络中，VFGs，特别是那些控制黏附、生物膜形成和运动性的基因，既作为胁迫耐受决定因子，又作为微生物合作的促进剂。MGEs的同时富集和活跃表达为VFG传播提供了分子管道，可能加速群落内的HGT。尽管VFG富集主要反映对金属胁迫的生态适应而非内在的致病性转变，但长期后果是机会性病原体在污染土壤中的持续选择和积累。这种时间累积效应将砷污染水稻土转化为VFGs和病原体的活跃储存库，通过稻米消费对食品安全和人类健康构成隐藏威胁。未来的农业和生物修复策略必须纳入基于RNA的风险评估，并明确解决重金属污染下毒力性状的共选择问题，以确保水稻农田生态系统的生物安全。