《Journal of Animal Science and Biotechnology》：Rumen DNA virome plasticity and viral metabolic potential are associated with seasonal adaptation in grazing yak and cattle on the Qinghai-Tibet Plateau
编辑推荐：
背景：作为瘤胃生态系统中多样且丰富的组分，病毒与其他微生物相互作用，被认为可影响微生物代谢和宿主生产力。然而，瘤胃病毒组在极端环境中如何响应季节性波动仍所知有限。本研究采用宏基因组分析，调查了青藏高原放牧牦牛和牛在暖季与冷季中瘤胃病毒组的病毒多样性、功能潜力及
背景：作为瘤胃生态系统中多样且丰富的组分，病毒与其他微生物相互作用，被认为可影响微生物代谢和宿主生产力。然而，瘤胃病毒组在极端环境中如何响应季节性波动仍所知有限。本研究采用宏基因组分析，调查了青藏高原放牧牦牛和牛在暖季与冷季中瘤胃病毒组的病毒多样性、功能潜力及病毒-宿主关联的时间动态。
结果：牦牛和牛的瘤胃病毒群落均表现出明显的季节性变化，冷季的α多样性高于暖季。跨季节比较显示，牦牛瘤胃病毒组的α多样性和群落稳定性均高于牛。研究共鉴定出27,353个温和病毒操作分类单元（vOTUs）和31,976个烈性病毒操作分类单元（vOTUs），主要归属于Caudoviricetes纲。这些病毒与覆盖24个细菌门和8个古菌门的微生物宿主相关联，其中拟杆菌门（Bacteroidota）和芽孢杆菌门（Bacillota）为优势谱系。冷季的病毒-宿主关联数量更多，且在牦牛和牛之间表现出不同的宿主特异性模式。冷季病毒组中抗防御基因多样性降低，与脂肪酸代谢和半纤维素降解相关的辅助代谢基因（AMGs）发生富集。值得注意的是，冷季中牦牛与牛之间的差异更为显著：牦牛瘤胃病毒组富集了与氨基酸、脂质和能量代谢以及纤维素降解CAZyme家族相关的途径，而牛瘤胃病毒组则富集了一般碳水化合物代谢以及复制和修复过程。
结论：瘤胃DNA病毒组的季节性可塑性以及牦牛和牛之间明显的种间差异，为理解它们在青藏高原严酷环境中不同的微生物过程提供了见解。这些发现提示，瘤胃DNA病毒组对季节变化表现出复杂的生态和功能响应，并可能与环境应激下的宿主-微生物互作及营养利用相关。本研究强调了瘤胃病毒基因组在理解高海拔反刍动物病毒-微生物互作、微生物适应和营养利用中的生态意义。
论文解读
《Journal of Animal Science and Biotechnology》发表的这项研究，针对青藏高原放牧牦牛和牛的瘤胃DNA病毒组开展系统分析。已有研究认识到瘤胃病毒可通过携带辅助代谢基因（AMGs）和裂解宿主等方式影响微生物群落及宿主性状，但极端环境条件下放牧反刍动物瘤胃病毒组的季节动态及其在宿主适应中的潜在作用仍不清楚。研究人员利用此前建立的纵向队列，采集12头去势雄性个体（6头柴达木牛和6头天祝白牦牛，年龄36±4个月，体重分别约230±15 kg和180±20 kg）在冷季和暖季共4个时间点的瘤胃液样本，通过宏基因组测序和病毒组生物信息学分析，揭示了病毒组成、功能基因、抗防御基因以及病毒-宿主互作的季节和宿主特异性变化，提出瘤胃DNA病毒组可能参与介导高原放牧反刍动物对季节性营养应激的适应。
在技术方法上，研究人员主要采用以下关键流程：对48份瘤胃液样本进行总DNA提取和Illumina HiSeq X Ten双端宏基因组测序；经质量控制、宿主序列去除和MEGAHIT组装后，利用VirSorter2和DeepVirFinder预测病毒序列，CheckV评估基因组完整性，CD-HIT按95%一致性和85%覆盖度去冗余形成病毒操作分类单元（vOTUs）；使用PhaGCN进行病毒分类，iPHoP预测微生物宿主，PhaTYP推断病毒生活方式；功能注释采用KofamScan比对KEGG数据库、DIAMOND比对CAZy数据库鉴定碳水化合物活性酶（CAZymes），并通过AntiDefenseFinder识别抗防御基因；统计方法包括Wilcoxon秩和检验、PERMANOVA、DESeq2差异丰度分析和Spearman相关性分析。样本队列来自同一自然放牧群体，覆盖冷季（11月和1月）与暖季（4月下旬和8月）。
研究结果部分有以下主要发现。
一、瘤胃病毒组概况。通过病毒序列鉴定和去冗余，共获得112,777个vOTUs，其中57.35%可进行病毒分类，未分类Caudoviricetes占81.6%，其余主要涉及Peduoviridae、Kleczkowskaviridae等科。生活方式预测得到27,353个温和与31,976个烈性vOTUs。约48.8%的vOTUs可预测到微生物宿主，宿主覆盖24个细菌门和8个古菌门，Bacteroidota和Bacillota为优势宿主谱系；属水平上Prevotella是最常见的预测宿主，其次为Cryptobacteroides、Limimorpha、Butyrivibrio、Ruminococcus等。共预测到88,599个病毒-宿主关联，其中部分古菌门仅与烈性病毒或仅与温和病毒关联。
二、病毒群落随季节和宿主物种的变化。通过α多样性、主坐标分析（PCoA）和差异丰度分析发现，冷季病毒群落α多样性显著高于暖季，季节和宿主物种均显著影响群落结构（季节R2=0.18，物种R2=0.09）。无论冷季还是暖季，牦牛均表现出比牛更高的病毒多样性和更低的组内Bray-Curtis距离，提示其病毒群落更加稳定。冷季中，两物种的温和病毒相对丰度显著升高，且差异丰度病毒科的数量在冷季更多。
三、病毒编码代谢基因和CAZymes的季节与宿主差异。通过KEGG功能注释和CAZy数据库比对，研究人员发现病毒功能谱受季节和宿主物种共同影响。牛的瘤胃病毒组在暖季富集氨基酸代谢和碳水化合物代谢相关通路，冷季则富集脂肪酸代谢通路；牦牛在暖季富集氨基酸降解和糖酵解等通路，冷季富集氨基酸生物合成、脂质代谢、半纤维素和纤维素降解相关CAZymes。冷季种间比较显示，牦牛病毒组富集氨基酸、脂质和能量代谢以及纤维素降解CAZyme家族，而牛病毒组富集一般碳水化合物代谢和复制修复相关过程，表明冷季环境压力可能加剧了病毒功能潜力的宿主特异性分化。
四、抗防御基因的季节和宿主模式。通过AntiDefenseFinder鉴定抗防御基因并分析其多样性，结果发现冷季抗防御基因α多样性在牦牛和牛中均显著降低，PCoA显示季节是主要驱动因素。牛中3个抗防御基因存在显著季节差异，牦牛中则有7个；牛独有vcrx091、vcrx093和tad2等基因，且vcrx093和mom相对丰度高于牦牛。这些结果提示，病毒抗防御系统的组成既受季节环境影响，也受宿主特异性的病毒-宿主共进化背景调节。
五、病毒-宿主互作的季节和宿主特异性。通过病毒宿主预测和Spearman相关分析，研究人员发现冷季预测的病毒-宿主关联数量略高于暖季，烈性病毒比温和病毒具有更宽的宿主范围。Prevotella仍是两物种中最常见的宿主属，Butyrivibrio、Ruminococcus、Fibrobacter和Methanobrevibacter等纤维降解菌和产甲烷菌在冷季的病毒关联数量增加。相关性分析显示，牛的Prevotella-病毒相关方向从冷季负相关转为暖季正相关；牦牛在暖季中Bacteroides与烈性病毒负相关、与温和病毒正相关；Methanobrevibacter在两物种中均与温和病毒正相关、与烈性病毒负相关，而Fibrobacter、Butyrivibrio和Selenomonas与病毒丰度表现为一致的强正相关。
讨论部分指出，病毒群落结构的季节和宿主差异可能反映环境压力和宿主因素共同塑造了瘤胃病毒生态。冷季烈性病毒富集可能促进微生物细胞裂解和养分再循环，有助于低质牧草条件下瘤胃生态系统的物质流动。Caudoviricetes为主导类群与既往研究一致，而Peduoviridae为主要科则可能与宿主、环境和分析方法差异有关。预测感染Methanobrevibacter的病毒多为烈性，提示这类病毒可能在调控产甲烷菌群和甲烷排放方面具有潜在意义。冷季AMGs的富集，尤其与脂肪酸代谢和纤维降解相关的基因，提示病毒可能通过调节宿主代谢支持微生物在营养受限条件下的适应。牦牛病毒组在冷季更强的代谢基因富集和更高的群落稳定性，可能与其在青藏高原恶劣环境中更优的饲料利用效率有关。抗防御基因多样性在冷季下降，可能与微生物群落资源分配和宿主防御压力降低有关。研究也存在样本量有限、vOTU质量评估依赖参考数据库、宿主预测覆盖不全以及未涵盖RNA病毒等局限。
结论部分翻译如下：总之，本研究表明，青藏高原放牧牦牛和牛的瘤胃DNA病毒组在暖季至冷季的环境转变中表现出显著的季节性和宿主特异性动态。病毒群落组成、病毒-宿主关联、AMGs和抗防御基因的季节变化提示，病毒可能参与微生物对营养受限冬季条件的适应。冷季中与能量代谢、脂肪酸代谢和植物纤维降解相关的AMGs富集，突出了病毒在支持微生物代谢灵活性和生态系统恢复力方面的潜在作用。与牛相比，牦牛表现出更高的病毒组稳定性和更强的代谢相关AMG富集，反映出不同宿主对环境挑战的特异性病毒组响应。总体而言，这些发现为高海拔放牧系统中的瘤胃病毒生态学提供了新认识，并提示病毒组可能是季节性营养应激下微生物介导宿主适应的重要组分；该研究也为在极端环境中基于微生物组策略提升反刍动物适应性和生产力提供了潜在参考。
