基于环境DNA监测罗氏沼虾（Macrobrachium rosenbergii）病原体在孵化场和自然环境中的检测与传播

《Environmental DNA》：Environmental DNA Monitoring for Detection and Spread of Prawn (Macrobrachium rosenbergii) Disease Organisms in Hatcheries and Natural Conditions

　　

由罗氏沼虾诺达病毒（Macrobrachium rosenbergii nodavirus，MrNV）、极小病毒（extra small virus，XSV）、罗氏沼虾金病毒（M. rosenbergii golda virus，MrGV）和致病性弧菌属（Vibrio spp.，包括副溶血弧菌V. parahaemolyticus、创伤弧菌V. vulnificus、哈维弧菌V. harveyi）引起的病原体爆发已严重破坏孟加拉国淡水沼虾养殖业，因传统组织诊断的延迟导致大规模死亡和数百万美元损失。研究人员每两周或每季度（2023年7月至2025年4月）从10个孵化场和10个相连的自然位点收集1–5?L水样，经0.45?μm过滤，于4°C乙醇中保存，并使用Invitrogen PureLink试剂盒提取环境DNA（eDNA）（A260/280 1.7–1.9）。定制TaqMan qPCR检测靶向MrNV RNA依赖的RNA聚合酶（RdRp，112?bp）、XSV衣壳-非翻译区（UTR，145?bp）、MrGV开放阅读框1a-ORF3（ORF1a-ORF3，189?bp）和弧菌toxR基因（98?bp），采用质粒标准品（效率91%–97%，R²?≥?0.997），检测限95%（LOD95）为6–14拷贝/反应（遵循CLSI EP17-A2标准），并包含内部对照。阳性样本（Ct?≤?37）经三重复验证，Sanger测序（n?=?152，ABI 3730xl，BigDye v3.1），BLASTn分析（85%–100%同源性，e?≤?10^?60），以及系统发育评估（IQ-TREE，5000次自举）。孵化场的中位数载量（log₁₀拷贝/L，Wilcoxon检验p??10）比自然位点高4.7倍：MrNV为1300拷贝/L（94%检出率，95%置信区间82%–99%）对比450拷贝/L（41%，29%–54%）；弧菌为1600拷贝/L（88%）对比1400拷贝/L（76%）；XSV/MrGV为850/450拷贝/L对比320/280拷贝/L——且频繁出现共感染（SparCC ρ?=?0.72，例如84%的孵化场存在MrNV-XSV共感染）。负二项式广义线性混合模型（GLMM）（边际R²_marg?=?0.65，条件R²_cond?=?0.77）确认了孵化场废水溢出效应（系数=2.03，p??11），受距离/盐度调节（n?=?10对），生物安全相关性r?=??0.62（p?=?0.01）；eDNA在pH?7–8.5、盐度0.01–36?ppt、温度26.5°C–31°C条件下表现最佳。变异包括MrNV单核苷酸多态性（SNPs，6.4?±?2.1，F_ST?=?0.018孵化场对比0.037自然位点）、新型XSV 6个核苷酸缺失（23%自然位点）、MrGV氨基酸替换（V134I、N212D、G345S）以及弧菌进化枝（F_ST?=?0.29–0.48）。这些结果验证了eDNA作为非侵入性早期预警工具的有效性，建议进行常规孵化场筛查、紫外线/氯处理以及自然监测，以遏制病原库并维持生产。

**研究背景与问题**

孟加拉国淡水沼虾（*Macrobrachium rosenbergii*）养殖是全球食物安全的重要组成部分，但长期受病毒性（如罗氏沼虾诺达病毒MrNV、极小病毒XSV、罗氏沼虾金病毒MrGV）和细菌性（如致病性弧菌属*Vibrio* spp.）病原体的威胁。传统诊断方法（如组织病理学、常规PCR）依赖已出现临床症状的宿主组织，具有侵入性、周期长且滞后，导致暴发时已造成大规模死亡和数百万美元损失。环境DNA（eDNA）监测技术已在废水SARS-CoV-2追踪和鱼类养殖寄生虫爆发中展现潜力，但在淡水沼虾体系中缺乏标准化方案，尤其在区分活性病原体与环境碎屑、评估环境因子（温度、盐度、pH）对检测可靠性的影响方面存在空白。同时，孵化场与自然水体的病原体环境库和传播路径尚不明确，制约了主动式生物安全策略的制定。

**研究内容与结论**

该项研究于2023年7月至2025年4月，在孟加拉国西南部选取10个罗氏沼虾孵化场及其相邻的10条河流/河口（共20个位点），采集水样开展eDNA监测。研究人员开发了四种病原体特异性TaqMan qPCR检测体系（靶向MrNV RdRp、XSV衣壳-UTR、MrGV ORF1a-ORF3、弧菌toxR），结合Sanger测序验证和系统发育分析，量化了孵化场与自然水体的病原体负荷、流行率及遗传多样性，并分析了环境参数（pH、盐度、温度、溶氧）和生物安全措施对eDNA检测效果的影响。结果表明：孵化场中病原体的中位数负荷比自然位点高4.7倍（Wilcoxon p??10），MrNV和弧菌的流行率分别达94%和88%；共感染普遍（如84%孵化场存在MrNV-XSV共感染）；负二项混合模型确认了孵化场废水向自然水体的溢出效应（系数2.03，p??11），生物安全评分与病原体丰度呈显著负相关（r?=??0.62）；eDNA在pH 7–8.5、盐度0.01–36 ppt、温度26.5–31°C时检测性能最佳。测序发现MrNV单核苷酸多态性（SNPs）、XSV新型6 bp缺失、MrGV氨基酸替换及弧菌进化枝分化（F_ST分别为0.18、0.29–0.48）。该研究在《Environmental DNA》期刊发表，验证了eDNA作为非侵入性早期预警工具的有效性，建议将常规eDNA筛查、紫外线/氯消毒处理及自然水体监测纳入生物安全流程。

**关键技术方法简述**

研究人员采用以下主要技术：①样本采集：每两周（孵化场）或每季度（自然位点）采集1–5?L表层水，经0.45?μm纤维膜现场过滤后，于95%–100%乙醇（4°C）中保存；②eDNA提取：使用Invitrogen PureLink环境DNA试剂盒，经蛋白酶K裂解、硅胶膜纯化，获得高纯度DNA（A260/280 1.7–1.9）；③定量检测：设计并验证了四种TaqMan qPCR体系（引物/探针均为新设计），靶向保守基因区域（MrNV RdRp 112?bp、XSV capsid-UTR 145?bp、MrGV ORF1a-ORF3 189?bp、弧菌toxR 98?bp），质粒标准品建立标准曲线（效率90%–99%，R²?≥?0.995），LOD95为6–14拷贝/反应；④测序验证与遗传分析：阳性样本（Ct?≤?37）经Sanger测序（n?=?152，ABI 3730xl，BigDye v3.1），BLASTn比对（85%–100%同源性），并构建系统发育树（IQ-TREE，5000次自举）。样本来源为孟加拉国西南部10个孵化场与10个相连自然河流/河口。

**研究结果**

*3.1 基于eDNA qPCR的病原体定量与流行率*
通过qPCR检测，孵化场水样中MrNV的平均负荷为3.7?×?10⁵拷贝/L（85%样本阳性），自然位点为1.1?×?10⁴拷贝/L（35%阳性）；XSV与MrGV也呈现类似趋势。卡方检验显示孵化场中所有病原体的流行率均显著高于自然位点（MrNV χ²?=?12.8，p?2?=?10.8，p?2?=?5.1，p?=?0.02；弧菌χ²?=?7.2，p?=?0.007），表明集约化养殖导致病原体累积。

*3.2 致病性弧菌的鉴定与定量*
靶向toxR基因的qPCR显示，孵化场弧菌平均负荷为9.7?×?10⁴拷贝/L（65%阳性），自然位点仅3.1?×?10³拷贝/L（20%阳性）。Sanger测序确认主要种为副溶血弧菌（*V. parahaemolyticus*）和创伤弧菌（*V. vulnificus*），且在部分孵化场样本中检出毒力基因trh。

*3.3 qPCR引物与探针验证*
所有引物/探针在BLASTn中与靶标（MrNV 120株、XSV 45株、MrGV 12株、弧菌250株）同源性95%–100%，与非靶标（500种其他诺达病毒、1000种非致病弧菌）同源性<85%–88%；体外验证显示无交叉反应（Ct?>?38）。qPCR效率90%–99%，分析LOD95为10拷贝/反应。

*3.4 Sanger测序确认病原体身份与毒株变异*
对152个扩增子测序，BLAST比对显示MrNV同源性97%–100%（含6.4?±?2.1个SNP/扩增子），XSV同源性94%–99%（发现23%样本含新型6 bp缺失），MrGV同源性85%–96%（ORF3区氨基酸替换，如V134I、N212D、G345S），弧菌中65%为副溶血弧菌（98%–100%同源性）、35%为创伤弧菌（92%–95%同源性），无脱靶扩增。

*3.5 病原体的时空分布格局*
qPCR定量显示孵化场病原体负荷在各时间点均高于自然位点，季节波动明显（暖季高峰期）。空间上，病原体eDNA集中于孵化场排水口及高密度养殖池附近，自然水体样本检出率低且零散，反映稀释与环境降解作用。

*3.6 环境参数与生物安全的相关性*
温度与病原体丰度呈正相关（r?=?+0.68），溶氧呈负相关（r?=??0.54），pH影响弱（r?=?+0.21）。生物安全评分（1–10分）与病原体丰度呈显著负相关（r?=??0.62），即生物安全执行越好，病原体负荷越低。

*3.7 孵化场与自然水体的eDNA检测结果汇总*
孵化场样本（H1–H10）中病原体总检出率78%，最高为MrNV（1250–1300拷贝/mL）和弧菌（1600拷贝/mL）；自然位点（R1–R10）总检出率45%，MrNV最高450拷贝/mL，弧菌1400拷贝/mL，整体低于孵化场。

*3.8 病原体检测扩增子的遗传多样性与系统发育关系*
基于Sanger测序序列的Kimura-2-parameter遗传距离矩阵显示，同一病原体内部变异低（0.002–0.045替换/位点），不同病原体间距离>0.25。最大似然树（1000次自举）显示四个病原体形成独立单系分支，孟加拉国序列与亚洲参考株聚类，同时出现亚谱系（如MrNV 4个SNP变异、XSV缺失株、MrGV糖蛋白替换、弧菌两个致病进化枝）。

*3.9 病原体流行率、负荷、储库与传播*
组合分析确认孵化场中位负荷比自然位点高4.7倍（Wilcoxon p?2为0.65、条件R²为0.77，孵化场类型（系数2.03，p?
*3.10 Sanger测序最终验证的病原体身份与毒株变异*
152个扩增子中，MrNV的核苷酸多样性（π）孵化场为0.018、自然位点为0.037；XSV新型6 bp缺失主要存在于自然样本（23%）；MrGV氨基酸替换位于糖蛋白B结构域；弧菌F_ST值为0.29–0.48，表明显著分化。所有序列满足Phred?≥?32，无假阳性。

**讨论总结与结论翻译**
讨论指出eDNA监测是高灵敏度、非侵入性的工具，孵化场中MrNV检出率90%–98%、弧菌共感染70%–75%，显著高于既往报道，且自然位点检出率上升，提示病原压力加剧。新发现的毒株变异（如XSV缺失、MrGV替换）表明病原多样性增加。环境参数（pH 7–8.5、盐度<36 ppt、温度<31°C）是eDNA可靠检测的关键。生物安全措施（r?=??0.62）可有效降低病原负荷。方法上，严格的时空复现、双盲对照及测序验证提升了结果的可靠性，但样本量有限且qPCR无法区分活/死病原体是局限性。未来需扩大地理范围、结合活力检测（如PMA-qPCR）和现场快速检测平台。
研究结论（5 Conclusion）翻译如下：本研究确立了eDNA-qPCR-Sanger整合方法作为罗氏沼虾病原体主动监测的黄金标准，确定了孵化场是病原体热点（4.7倍载量；MrNV 1300拷贝/mL，位点H5），驱动弧菌/MrNV溢出（GLMM β = 2.03，p?ST = 0.29–0.48）和生物安全相关性（r?=??0.62）要求紧急行动：实施常规eDNA孵化场筛查、废水UV/氯化处理、在最优条件下（pH 7–8.5、盐度<36 ppt）开展自然监测。这些基于证据的策略可减轻暴发、遏制经济损失，并推动可持续沼虾养殖在疾病流行地区实现Q1影响力的抗逆性发展。