凡纳滨对虾仔虾自然感染BP、WSSV及AHPND致病菌育苗池的微生物群落结构与功能特征研究

《PLOS One》：Microbial community structure of Penaeus vannamei post-larvae naturally infected with multiple viruses in hatchery tanks

　　

背景：该研究对凡纳滨对虾（Penaeus vannamei）仔虾中病毒感染与细菌感染如何与微生物组变化及存活共同出现进行了观察性评估；该物种是全球水产养殖中经济重要性最高的物种之一。方法：研究人员依据组织学和PCR将11个商业幼体生产池划分为三种感染情形：Baculovirus penaei（BP）单感染、混合感染（BP、白斑综合征病毒（white spot syndrome virus，WSSV）及引起急性肝胰腺坏死病（acute hepatopancreatic necrosis disease，AHPND）的细菌）以及未感染对照。研究人员对高分子量DNA在Illumina NovaSeq 6000平台上进行鸟枪法宏基因组测序，随后评估微生物多样性、群落结构和感染相关微生物特征。结果：尽管BP是感染组中最丰富的病毒属，严重死亡仅见于混合感染条件下，且与引起AHPND的细菌存在及微生物组紊乱密切相关，而非仅与BP相关。深度宏基因组测序显示，与未感染个体相比，感染BP的虾及混合感染虾的微生物组结构降低且独特。未感染仔虾维持复杂且平衡的细菌群落，主要由Donghicola和Muricauda组成。Nudiviridae和Oscillospiraceae在BP单感染组和混合感染组中显示相当的相对丰度。相反，混合感染的特征为弧菌属（Vibrio）相对增加，并伴随PCR检测到PirAB毒素基因；而WSSV尽管相对丰度较低，仍在宏基因组数据集中可检测到，同时该感染在仔虾中呈热抑制和无症状状态。在功能上，感染仔虾的微生物组富集转座酶、逆转录转座子和毒力因子基因。研究人员从仔虾宏基因组中组装并注释了近完成BP-GR25基因组（119,245 bp）。结论：这些结果突出高通量测序在水产养殖病原早期检测中的价值；在仔虾阶段，细微微生物变化可作为生产健康状态的有价值指标。

研究背景与目的：全球对虾养殖自20世纪70年代在亚洲兴起并持续扩展，南美部分生产国对虾养殖占国民经济重要份额，但传染病发生率随之上升。Baculovirus penaei（BP；亦称PvSNPV，已归入Nudiviridae）主要侵染肝胰腺和前中肠上皮，可影响多个对虾物种各生活阶段，仔虾和幼虾早期尤其PL₉高度易感。白斑综合征病毒（white spot syndrome virus，WSSV）传播快、死亡率高、可感染各发育阶段。副溶血性弧菌（Vibrio parahaemolyticus）中携带PirA/PirB二元毒素的菌株引起急性肝胰腺坏死病（acute hepatopancreatic necrosis disease，AHPND），导致肝胰腺小管上皮严重降解脱落。细菌与病毒共感染常增强致病性和死亡，BP感染亦可能与继发细菌感染增加相关。对虾内源微生物组参与营养、免疫和抵御机会性暴发；病原可在健康个体低丰度持续而不致病，病理状态更与微生物群落失衡和环境因素相关，pathobiome框架强调群落组成与功能潜力（含毒力因子）。然而BP在虾微生物组中的角色、其相对丰度及与WSSV、AHPND致病菌在真实生产条件下的关系尚不清楚。该研究旨在用高通量测序（high-throughput sequencing，HTS）表征商业育苗池中P. vannamei仔虾在BP单感染及BP-WSSV-AHPND混合感染下的微生物组结构，识别与感染状态相关的微生物特征。研究人员对11个商业生产池的PL₁₀仔虾混合样本进行组织病理和PCR诊断，分为BP单感染、混合感染和非感染，并比较存活；对高分子量DNA进行Illumina NovaSeq 6000鸟枪法宏基因组测序，分析多样性、群落结构和功能；同时从仔虾宏基因组组装BP-GR25基因组。结果显示，BP是感染组最丰富病毒属，但严重死亡仅见于混合感染，且与AHPND致病菌存在及微生物组紊乱密切相关而非仅与BP相关；感染仔虾微生物组结构降低且独特，非感染仔虾保持复杂平衡的细菌群落；混合感染中Vibrio相对增加并PCR检出PirAB毒素基因，WSSV在宏基因组中仍可检测但呈热抑制和无症状；感染微生物组富集转座酶、逆转录转座子和毒力因子基因；获得近完成BP-GR25基因组。该研究突出HTS在水产病原早期检测中的价值，提示仔虾阶段细微微生物变化可作为生产健康状态指标。

关键技术方法（≤250字）：样本来自南美一商业孵化场11个35吨生产池的PL₁₀仔虾，每池混合样本分两份用于组织病理（30尾H&E）及PCR与HTS。以组织学和nested PCR（BP、WSSV、EHP、NHP、IHHNV及AHPND AP4）分为BP单感染4池、混合感染4池、非感染3池。高分子量DNA经Illumina NovaSeq 6000 PE150鸟枪法宏基因组测序；去宿主后以SqueezeMeta共组装，进行分类/功能注释、ACE/Shannon、Bray-Curtis PCoA/PERMANOVA、属水平Kruskal-Wallis及FDR、COG/KEGG/Pfam和VFGs分析；以DAS Tool、CheckV、VAPiD、fastANI、PhyloSuite及重映射/SNP分析组装BP-GR25。

结果——样本采集、诊断与感染状态分类（Sample collection, diagnosis, and classification of infection status）：H&E显示BP典型四面体包涵体，混合组兼具AHPND上皮脱落坏死；PCR确认BP、WSSV和AHPND致病菌；分组为4个BP单感染池、4个混合感染池、3个非感染池。PL₁₁存活为混合感染27.5±1.1%、BP单感染63.3±2.2%、非感染70.8±0.9%，Tukey HSD p<0.001；溶解氧、温度、盐度、pH无显著差异。结论：高死亡与混合感染相关，BP单独感染存活较高。

结果——宏基因组组装与注释（Metagenome assembly and annotation）：共获得732 million clean reads，平均67 million；2.9 million contigs；3,433,465个蛋白编码开放阅读框（open reading frames，ORFs）；108个宏基因组组装基因组（metagenome-assembled genomes，MAGs），其中15个高质量、25个良好质量。结论：为群落和功能分析提供基础。

结果——微生物组多样性与差异丰度（Microbiome diversity and differential abundance）：非感染组几乎为细菌，Donghicola 11.03%、Muricauda 10.19%优势；感染组病毒reads高（Naldaviricetes、Nimaviridae），Oscillospiraceae约14%最丰富；Vibrio在混合组较非感染高3.2倍，BP单感染低于非感染但Vibrio群落组成显著不同（PERMANOVA p=0.045）；Baculovirus在感染组较非感染高3.3倍，Whispovirus bin与WSSV高度相似且约占总微生物组1%；ACE和Shannon在感染组降低，Bray-Curtis PCoA分离（PERMANOVA p=0.001）。结论：感染与多样性降低和结构改变相关，混合感染伴Vibrio相对增加。

结果——微生物组功能与相关毒力因子（Microbiome function and associated virulence factors）：感染样本富集转座酶、逆转录转座子、失活衍生物、逆转录酶等COG；COG总体组成显著不同（PERMANOVA p<0.001）；非感染组KO功能更多，包括信号与细胞过程、转运、DNA修复重组、线粒体生物发生等，合计约10–15%；感染组KO功能数量和丰度显著降低；毒力相关未表征膜蛋白在感染组高，温度依赖性蛋白分泌的avirulence locus蛋白在非感染高；Pfam显示溶血素和SrfB在混合感染占优。结论：感染改变功能潜力并富集毒力相关功能，混合感染更明显。

结果——Baculovirus penaei基因组（Baculovirus penaei genome）：BP-GR25为119,245 bp，GC 33.3%，平均覆盖42×，122个ORFs、91个CDSs，CheckV完整性96.13%；混合组样本独立映射覆盖100%或99.99%，深度32.1×–81.6×，无高质量SNP；COG富集脂质运输代谢、信号转导、胞内运输分泌囊泡运输、复制重组修复及未知功能；与Penaeus vannamei nudivirus（OM066077.1）平均核苷酸一致性（average nucleotide identity，ANI）99.45%，聚类于Penaeus monodon nudivirus。结论：获得近完成、高度保守的BP-GR25基因组。

讨论总结：原文指出实验条件下BP对mysis₂至PL₂₅易感，PL₉最易感；本观察中混合感染池Vibrio相对增加且高死亡仅限混合感染，BP单独不关联高死亡；BP单感染Vibrio总丰度最低但群落组成显著不同，提示BP可能与Vibrio群落结构改变相关；由于无BP时未检出WSSV和AHPND，BP是否促进后续病原建立不能确定，需控制挑战实验。WSSV仅PCR和测序检出而无临床或严重病理，与>30 °C高温抑制WSSV基因表达和主动复制但未必清除基因组一致；NHP相关细菌和Whispovirus在感染组检测，但其他NHP细菌或额外病毒病变未显示，原文归因于检测敏感性或低病毒载量。非感染组Donghicola和Muricauda显著多，前者有保护报道，后者编码抗氧化合成并与微藻形成抗Vibrio尤其AHPND的consortium。感染组Baculovirus占优且VFGs多样性和丰度高，类似dysbiosis；非感染组KO功能丰富提示代谢功能多样和韧性。BP-GR25与OM066077同属Gammanudivirus，长度和ORFs有差异，含polyhedrin基因及A/T-rich调控motif，后者作用需实验确认。组织病理与PCR对BP检测存在不一致，引物错配不是唯一解释，池内异质分布也可能参与，需改进分子工具。

研究结论翻译：据原作者所知，该研究首次表征自然感染BP的P. vannamei微生物组并考察其与WSSV及AHPND相关细菌的相互作用；结果显示BP感染与仔虾微生物组丰富度和多样性显著降低相关；BP感染与混合感染样本群落结构大体相似且不同于非感染组，同时功能及毒力相关基因增加，提示共感染条件下病毒动态改变的潜在联系；尽管BP是感染组最丰富病毒属，死亡仅见于混合感染，提示疾病结局与复杂病原共现和微生物组紊乱而非BP单独密切相关；WSSV在混合感染微生物组中仍可检测，无临床迹象，符合商业孵化场>30 °C下温度介导的covert latency；理解这些微生物关联可为疾病结局变异提供见解，并为促进微生物组稳定和虾韧性综合管理策略奠定基础，尚需实验验证。