降雨与底物来源驱动塞浦路斯瘤胃动物粪便微生物组变化：多样性趋同、Bacillota富集与功能潜力重塑

《Environmental Microbiology Reports》：Rainfall and Substrate Origin Driven Change of Ruminant Manure Microbiomes in Cyprus

　　

微生物群落驱动粪便储存期间有机质稳定；然而，底物来源（substrate origin）与突发降雨事件对这些群落的联合影响仍知之甚少。研究人员首次在塞浦路斯利用16S rRNA宏条形码（16S rRNA metabarcoding）和鸟枪法测序（shotgun sequencing），调查了粪便类型（奶牛 vs. 山羊）以及干燥–降雨模拟脉冲（drying–rainfall simulation pulse）对粪便储存期间群落组成、多样性和功能潜力的影响。在干燥条件下，底物组成设定基线，结构异质性更高的山羊粪支持显著高于奶牛粪的多样性。随后，降雨模拟作为主导环境过滤器，急剧降低多样性并驱动跨粪便类型的组成趋同。该转变伴随快速响应、耐胁迫类群的迅速增加，主要是Bacillota（厚壁菌门），以及与孢子形成和快速生长相关功能的增加。排序和离散分析进一步表明，底物来源调节稳定性：降雨模拟后，奶牛粪群落表现出比山羊粪群落更大的组成位移和更大的变异性。总体而言，粪便储存系统的功效和稳定性源于初始理化模板与动态水分胁迫之间的相互作用。优化大规模操作将需要纳入底物特征、微生物群落结构和水分动态的管理方法。

论文解读：降雨模拟与底物来源驱动塞浦路斯瘤胃动物粪便微生物组变化

研究背景与问题：微生物群落是粪便储存和堆肥中有机质分解与稳定的核心，驱动碳矿化、养分转化和有机废物稳定。其结构和功能受底物化学组成和环境条件等生态过滤器影响，其中粪便类型和水分可利用性尤为关键，但二者联合影响尚不清楚。塞浦路斯等地中海国家夏季长期干燥，突发强降雨因全球变暖更频繁。降雨进入储存粪便，在干燥基线上形成快速湿润脉冲，改变孔隙连通、氧气扩散、氧化还原电位、有机底物溶解及养分和金属迁移，产生短促选择压力，可能有利于富营养型（copiotrophic）或耐胁迫、快速响应、适应水分的微生物，并改变推断功能潜力。虽然水分是堆肥主变量，但夏季储存中意外降雨情境，包括短暂时间、空间斑块及与干燥的相互作用，仍研究不足。反刍动物粪因物种消化和日粮不同而理化性质差异显著，山羊粪通常干物质更高、结构更复杂，因摄入更多纤维植物残渣；这可能影响微生物定殖和生态位多样性。此前较少在受控粪便储存条件下直接比较底物来源如何调节微生物生态。因此，研究人员在真实农作条件下，研究粪便类型（奶牛 vs. 山羊）和降雨模拟（干 vs. 湿）如何交互塑造群落组成和多样性，并聚焦初始底物对α多样性（alpha diversity）和β多样性（beta diversity）影响、不同水分下时间演替、以及底物与水分交互对降雨模拟期间演替的影响。该研究为粪便储存和堆肥设计管理提供底物–环境互作新认识。

研究设计与关键技术方法：样本来自塞浦路斯农业研究所（Agricultural Research Institute of Cyprus，35°08′25″N—33°23′02″E，海拔170 m）Lefkosia，2023年夏季（7–8月）。奶牛为泌乳Holstein-Friesian，山羊为Damascus (Shami)，来自实验棚。设四个时间点：T0（第0天）、T1（第14天，降雨模拟前）、T2（第18天，降雨模拟后）、T3（第27天，结束）。每类型四个独立chamber，每动物×时间点四生物重复。降雨模拟每chamber均匀加入4 L水，雨后3天采样；T0–T1为干期，T2–T3为湿期。关键技术包括16S rRNA宏条形码（16S rRNA metabarcoding）测V3–V4区、Illumina NextSeq 2000测序，鸟枪法宏基因组测序（shotgun metagenomics），DADA2和SILVA v138注释，R中microbiome、vegan、phyloseq、microeco分析，MetaBakery/HUMAnN/UniProt功能注释，以及线性混合效应模型、PERMANOVA（置换多元方差分析）、betadisper（多变量离散度检验）、DESeq2、指示种分析等。

结果：一、微生物群落组成和多样性对降雨模拟的响应。16S rRNA测序显示奶牛和山羊粪均由Bacillota（厚壁菌门）、Actinomycetota（放线菌门）、Pseudomonadota（假单胞菌门/变形菌门）主导，初始组成差异较小。奶牛粪Bacillota为46.3%–50.7%，Actinomycetota为22.5%–31.1%，Pseudomonadota为13.8%–22.4%，Bacteroidota（拟杆菌门）为2.7%–6.3%；山羊粪Bacillota为52.3%–63.6%，Actinomycetota为18.2%–23%，Pseudomonadota为8.5%–15.3%，Bacteroidota达9.9%。线性混合效应模型显示粪便类型、水分条件及交互显著。干期山羊粪Shannon多样性（4.69）显著高于奶牛粪（4.51，p=0.0012）；湿期山羊为4.32、奶牛为4.30（p=0.73），差异消失。两者湿期多样性下降，山羊下降更显著（β=?0.166，t=?2.32）。chamber效应可忽略。

二、群落排序与多变量离散。Bray–Curtis相异度主坐标分析（PCoA）中，奶牛与山羊沿PCo1（38.7%方差）分离，水分沿PCo2（19.1%方差）分离。PERMANOVA显示粪便类型（R²=29.2%，p=0.001）、水分（R²=7.9%，p=0.001）及交互（R²=2.2%，p=0.001）显著。奶牛粪在雨后组成位移和离散度更大。betadisper显示粪便类型间变异p=0.002，交互p=0.001，而水分单独p=0.318。

三、共享与特有ASV。奶牛粪总1366个ASV，共享373，干期特有615，湿期特有378；山羊粪总2013，共享504，干期特有1046，湿期特有463。核心约25%–27%；山羊干期特有约占52%，奶牛湿期特有占28%，高于山羊23%。

四、共享ASV累计相对丰度。奶牛粪共享ASV在干期占88.6%，湿期降至81.6%（p=0.043），湿期样本变异大；山羊粪共享ASV在湿期占87.4%，高于干期79.5%（p<0.001）。

五、差异丰度—奶牛粪。共15个ASV显著富集，干期5个（Bacillota、Bacteroidota、Pseudomonadota，如Solibacillus ASV365、UCG-005 ASV263）；湿期10个，log₂FC常>10，最高+28，包括Thermobacillus、Melghirimyces、Ornithinibacillus、Saccharomonospora，多个Bacillaceae产孢。

六、差异丰度—山羊粪。共13个ASV，干期7个，如Solibacillus（ASV18）、Clostridium sensu stricto 1（ASV364）；湿期6个，主要为Bacillota（Dombalicus、Meghirimyces、Ornithinibacillus）和Pseudomonadota；最大log₂FC≈12，无单一优势，重排较温和。

七、鸟枪法宏基因组功能基因。湿条件下显著富集抗性/转运功能和locus-tagged基因，CPM（每百万计数）湿期多高于干期，山羊干湿分离更大。孢子形成相关基因spmA、spmB、ytfJ湿期显著富集（FDR≤0.05）；mecA也湿期富集，但被鉴定为MecA genetic competence negative regulator（PRK02899），非经典甲氧西林耐药决定簇。

八、动物主效应。山羊偏倚rplL、zupT、regX3、AUC31_04925、ppk，涉及金属/磷酸盐稳态、能量缓冲和翻译能力；奶牛偏倚A6048_04875(caiD)、fabA、glyA。nirA在山羊中显著高于奶牛（FDR_animal<0.05），为胞质、铁氧还蛋白/NAD(P)H依赖性同化型亚硝酸还原酶，生成NH₄⁺用于生物合成。整体湿期促进富营养型和胁迫响应功能。

讨论部分总结：讨论认为模拟降雨触发强生态过滤，降低多样性并富集产孢、富营养和耐湿类群，同时富集胁迫和抗性相关基因。土壤和粪便不可定量比较，但水分响应模式平行。山羊粪干期多样性高与其纤维丰富、底物异质、生态位更广有关。降雨后共享ASV和总体丰富度下降，表明水分创造强过滤，偏向耐湿/富营养型。功能基因变化与分类重排一致，再湿润将推断功能潜力集中到较少水分相关功能。Shannon下降与水分驱动过滤一致，水添加可改变氧气、溶解性和氧化还原，选择耐湿低氧微生境。Bacillota持续优势，奶牛粪中多个Bacillota ASV雨后富集。宏基因组确认spmA、spmB、ytfJ湿期富集，反映产孢Firmicutes生命周期动态，而非单纯干胁迫响应。mecA的PRK02899身份支持水分相关调控/适应，而非甲氧西林耐药。山羊偏倚基因提示金属/磷酸盐稳态、能量缓冲和翻译容量差异；nirA山羊偏倚和湿润脉冲动员无机N一致。局限是shotgun提供推断功能潜力而非直接活性或表达；理化评估在companion emissions-focused study，限制直接评估水分介导理化变化对微生物重排贡献。

研究结论（翻译）：本研究提供了关于初始底物异质性和水分条件如何与早期粪便储存期间微生物群落重塑相关的见解。研究人员证明，初始微生物多样性和组成受粪便理化性质影响，更复杂、纤维化的山羊粪比奶牛粪支持更高的初始多样性。然而，随后的降雨模拟作为主要环境驱动力，通过显著富集产孢Bacillota驱动两组群落趋同。这种耐胁迫富营养型（copiotrophs）的富集，得到孢子形成相关和胁迫响应相关基因丰度增加的支持，反映了与水分波动适应性响应一致的模式。最终，该研究揭示，底物来源显著减弱了降雨模拟后的组成位移，山羊粪群落比奶牛粪群落表现出更优的稳态（homeostasis）/恢复力（resilience）。优化大规模粪便储存系统的功效和稳定性需要新模型，整合起始原料的理化性质、微生物群落的功能性状以及动态水分条件。未来研究应聚焦于将这些稳定性指标与碳和养分循环速率进行功能关联，以优化水分管理方案并部署生物强化（bioaugmentation）管理策略。