《Annals of Applied Biology》：Detection and molecular characterisation of beet chlorosis virus in sugar beet in the Czech Republic
编辑推荐：
甜菜黄化病毒（BChV; 物种*Polerovirus BCHV*) 是一种具有重要经济意义的病原体，可引起全球甜菜（*Beta vulgaris* L. subsp. *vulgaris*）显著减产。本研究首次报道了捷克共和国BChV的检出及其遗传多样性。研
甜菜黄化病毒（BChV; 物种*Polerovirus BCHV*) 是一种具有重要经济意义的病原体，可引起全球甜菜（*Beta vulgaris* L. subsp. *vulgaris*）显著减产。本研究首次报道了捷克共和国BChV的检出及其遗传多样性。研究人员从多个地点采集了37份表现症状的甜菜样品，利用单酶RTX-PCR（reverse transcriptase PCR-based X-PCR）进行筛查，在来自7个地点的12份样品中检测到BChV。所有阳性分离物通过Sanger测序分析衣壳蛋白-通读结构域（coat protein-readthrough domain, CP-RTD; ORF3–ORF5）区域，该数据集作为多样性分析的主要依据。选取代表性分离物PX088834 (B_V1_CZ; Bezno_CZ) 进一步通过高通量测序（high-throughput sequencing, HTS）获取近完整基因组序列。HTS获得长度为5779 nt的BChV近完整基因组序列，证实了该分离物预期的马铃薯卷叶病毒属（*Polerovirus*）基因组结构，全基因组成对序列同一性分析将其归入欧洲BChV谱系。12个捷克分离物的CP-RTD系统发育分析揭示了三个系统发育谱系，命名为A、B和B1。三个捷克分离物[L_V8; Jizerní Vtelno_CZ, L_V11; Něm?ice_CZ, L_V16; Radonice_CZ]与来自英国和美国的分离物聚类于谱系A，其余九个分离物则与欧洲其他分离物共同聚类于谱系B和B1。这些谱系在同一采样地点的共存与多次独立引入的假说一致。滑动窗口相似性分析证实了从密切相关欧洲分离物到最分歧美国分离物的清晰渐变分歧梯度。基于密码子的分析显示CP-RTD区域存在显著的净化选择（purifying selection）作用，重组分析在近完整基因组中鉴定出一次推测性重组事件，在CP-RTD数据集中鉴定出两次额外事件。这些发现拓展了当前对欧洲BChV分布和遗传多样性的认识，并为未来监测工作奠定了基础。
本研究发表于《Annals of Applied Biology》，旨在解决捷克共和国甜菜生产中BChV的首次检测与遗传多样性认知空白问题。蚜虫传播的病毒黄化病（virus yellows disease, VY）对欧洲甜菜产业构成严重威胁，BChV作为主要病原之一可造成8%–24%的糖产量损失，但其此前尚未在捷克共和国有报道记录。随着欧盟新烟碱类杀虫剂禁令的实施，蚜虫介体种群显著增长，VY相关病毒流行程度加剧，使得明确BChV在捷克的分布状况成为紧迫的植物检疫议题。此外，BChV高度变异的CP-RTD区域（ORF3–ORF5）对于病毒蚜传、系统运动及宿主适应至关重要，但该病毒在欧洲大陆的遗传结构、进化动力学及重组模式仍缺乏系统认知。
研究人员于2021年11月在捷克共和国15个田块的甜菜种植区采集37份表现黄化、矮化等症状的叶片样品，经单酶RTX-PCR检测、Sanger测序及BLASTn分析确认12份样品BChV阳性。选取代表性分离物PX088834进行HTS获得5779 nt近完整基因组序列，测序深度达68.23×，覆盖度99.4%。基于CP-RTD区域（约2180 bp）构建了包含12个捷克分离物及GenBank参考序列的数据集，采用最大似然法（JTT模型，1000次bootstrap）进行系统发育分析；通过SimPlot v3.5.1进行滑动窗口成对相似性分析（全基因组：600 bp窗口/50 bp步长；CP-RTD：200 bp窗口/20 bp步长）；运用MEGA12进行基于密码子的净化选择Z检验（Kumar法，Kimura 2参数模型）；利用DnaSP v6估算核苷酸多样性（π、Tajima's D、Gst、Nm等参数）；采用RDP5（整合RDP、GENECONV、Chimaera、MaxChi、BootScan、SiScan、3Seq、LARD八种方法）检测重组事件。
**BChV在捷克共和国的检测与分布**
研究人员从捷克共和国7个不同地点的12份甜菜样品中检测到BChV，证明该病毒已在捷克多个甜菜种植区存在。
**捷克BChV分离物的CP-RTD序列分析与系统发育多样性**
对12个捷克BChV分离物的CP-RTD区域进行Sanger测序和最大似然系统发育分析，核苷酸和氨基酸序列均将其分为三个主要谱系A、B和B1。三个分离物（L_V8、L_V16、L_V11）归入谱系A，与英国（ON924241）和美国（AF352025）分离物聚类；其余九个分离物归入谱系B和B1，与欧洲其他分离物聚类。谱系A与谱系B的核苷酸/氨基酸同一性分别为92%/87%，与谱系B1为92%/86%；谱系B与B1最为接近（99%/98%）。Jizerní Vtelno_CZ同一地点同时检出谱系A和谱系B/B1分离物，支持多次独立引入假说。
**BChV基因组测序与从头组装**
HTS分析分离物PX088834获得5779 bp近完整基因组，包含六个完整ORF。尽管以英国分离物MH271171为参考进行序列优化，全基因组成对同一性分析显示其与瑞士分离物ON738348（BChV-Villiger）最为接近（99.58%核苷酸/99.21%氨基酸同一性）。全基因组范围PX088834与各欧洲分离物的核苷酸同一性为98.19%–99.58%，氨基酸同一性为97.01%–99.33%。
**全基因组系统发育分析**
基于全基因组核苷酸和氨基酸序列的最大似然树将PX088834置于欧洲BChV谱系内，与瑞士ON738348亲缘关系最近。英国ON924241和美国AF352025形成独立谱系，与PX088834的核苷酸/氨基酸同一性分别仅为94.64%/89.45%和93.03%/87.18%。全基因组水平谱系A与谱系B/B1的核苷酸/氨基酸同一性分别为94%/90%，而谱系B与B1为99%/97%。
**滑动窗口成对序列相似性分析**
全基因组滑动窗口分析显示，PX088834与谱系B1分离物保持最高相似性（>99%），谱系B略低（一般<98%），谱系A最具分歧性（多处<95%）。全基因组未观察到显著相似性转换，提示该捷克分离物缺乏全基因组水平的镶嵌重组结构。CP-RTD区域分析证实捷克分离物与谱系B1相似性最高且均匀，向谱系A呈渐进分歧梯度；谱系A分离物在RTD相关部分显示若干明显的局部相似性下降。
**选择压力分析**
对完整CP-RTD数据集进行的基于密码子的净化选择Z检验确认存在显著净化选择（Z = 4.43, p < .01）。多数序列对的Z统计量为正（p < .05），表明同义替换持续超过非同义替换。与最具分歧的参考分离物AF352025（美国）和ON924241（英国）的比较显示最强的净化选择信号（p < .001）。
**核苷酸多样性与群体分化**
三个CP-RTD谱系的内部多样性差异显著：谱系B1（8条序列，6个单倍型，18个多态位点，单倍型多样性0.929，核苷酸多样性π = 0.00319）和谱系B（6条序列，6个单倍型，49个多态位点，单倍型多样性1.000，π = 0.01174）分别表现低度和中度多样性；谱系A（5条序列，5个单倍型，260个多态位点，单倍型多样性1.000，π = 0.06411）呈最高多样性，此高值主要源于该谱系纳入了地理距离遥远的英美序列。所有谱系Tajima's D为负但无统计学意义（p > .10）。谱系间Gst值极低（0.00076–0.02154），基因流估计值高（Nm = 11.36–328.70），表明谱系间存在持续或近期的遗传交流。
**重组分析**
全基因组重组分析在近完整基因组中检测到一次推测性重组事件，涉及5'端（nt 1–719）和3'端（nt 5503–末端）的P0和RTD（ORF5）开放阅读框，推断AF352024（英国）为主要亲本、MW367424（法国）为次要亲本，该事件由GENECONV、BootScan和3Seq支持（p < .05）。CP-RTD区域进一步检测到两个显著重组事件：分离物L_V11（Něm?ice_CZ）中L_V8为主要亲本、L_V16为次要亲本（断点nt 4304–5595，五种方法支持）；分离物L_V16中AF352024为主要亲本、L_V8为次要亲本（断点nt 4321–5387，三种方法支持）。
讨论部分指出，欧盟新烟碱类杀虫剂禁用后蚜虫介体种群增长及VY病毒显著流行，使得BChV的监测尤为重要。捷克的BChV存在谱系A和谱系B/B1在同一田块的共存现象，特别是Jizerní Vtelno_CZ地点同时检出两类谱系，且谱系A在三个地理分离地点均有分布，表明该分歧谱系并非局限于单一田块，可能通过感染种植材料或迁飞蚜虫经历了多次独立引入。根据国际病毒分类委员会（ICTV）标准，属内物种划界约为任何基因产物差异>10%，PX088834与美国AF352025的ORF3–ORF5氨基酸同一性为87.5%，低于该阈值，但单一比较不足以重新定义物种边界，现有证据支持将捷克分离物保留于*Polerovirus BCHV*物种内。CP-RTD区域的主导净化选择反映了衣壳蛋白和RTD在病毒衣壳化、系统运动及蚜传中的核心功能作用。重组热点集中于RTD相关区域与该区域内在无序的C末端结构域可能降低重组序列约束有关。
研究结论指出，本研究首次在捷克共和国检测到BChV并证明其在多个甜菜种植区的存在。CP-RTD分析揭示捷克BChV分离物包含两个不同系统发育谱系（三个归于谱系A，九个归于谱系B/B1），与可能的多次引入事件一致。净化选择是作用于CP-RTD区域的主导进化动力，谱系间总体遗传分化较低。在多个地理分离地点检测到分歧谱系（谱系A）凸显了在常规监测中纳入多样化基因组区域的重要性，因为仅针对最保守CP序列的诊断检测可能低估BChV多样性。未来监测应优先对阳性样品进行CP-RTD测序，并纳入时间序列采样以追踪playingfieldboth谱系在捷克甜菜栽培区的传播与进化。这些发现强调了持续监测和研究的必要性，以支持捷克及更广泛地区甜菜栽培中*Polerovirus BCHV*的有效管理。
