柯萨奇病毒A16型进化与重组形式的协同关联

手足口病作为一种常见的儿童传染病，长期以来一直是全球公共卫生领域的重要关注点。柯萨奇病毒A16型（CVA16）作为手足口病的主要病原体之一，近年来在全球范围内的流行趋势引起了科学界的高度重视。《Virus Evolution》期刊的最新研究，通过对全球CVA16序列的全面分析，揭示了该病毒的基因型进化与重组形式变化之间存在的协同关联，为手足口病的防控策略提供了新的科学依据。

研究方法与数据基础

该研究团队从全球公开的GenBank数据库中检索了所有CVA16序列，通过严格的筛选标准，最终构建了三个主要数据集：包含1663条序列的原始数据集、包含442条VP1编码区序列的基因分型数据集，以及包含345条近全长基因组序列的重组形式分析数据集。研究综合运用了系统发育推断、重组形式调查和群体结构分析等多种生物信息学方法，对CVA16的进化动态进行了全面解析。

CVA16的基因型多样性与进化动态

基于VP1编码区的系统发育分析显示，全球CVA16可分为三个基因型（A、B和D），其中B型又可细分为B1和B2两个亚型，B1亚型进一步分为B1a、B1b和B1c三个簇。研究发现，B1b簇已成为全球流行的主导基因型，而B2亚型自2000年后已消失，A基因型仅包括1951年在南非分离的原型株G10。值得关注的是，D基因型是一个新兴基因型，其出现可追溯至2002年8月，与其他基因型之间的遗传距离达到15.6%-30.6%，形成了全新的进化分支。

从进化动态来看，全球CVA16的遗传进化可分为四个阶段：1995年前的低遗传多样性阶段、1995年至2008年的快速多样性增长阶段、2008年至2013年的多样性波动阶段，以及2013年后的多样性下降阶段。总体而言，CVA16的进化在2013年达到峰值，此后呈现下降趋势。

全球传播路径与关键节点

通过贝叶斯系统地理学分析，研究团队描绘了CVA16的全球空间传播格局。分析识别出20条具有显著统计学支持的传播路径，其中欧洲、日本和中国是CVA16全球传播的关键源头地区，其次是泰国和印度。研究发现，CVA16（G-10原型株）1951年首次在南非被分离后，随后传播至日本、中国、德国、喀麦隆和韩国，最终扩散至全球各地。

研究还观察到，地理距离与病毒迁移率之间存在弱负相关关系，表明CVA16更倾向于在较短距离内传播，这一现象在欧洲内部的传播网络中得到了充分印证。尽管存在一些长距离传播事件，但总体上病毒传播以区域性扩散为主。八条具有强支持证据的共有传播路径（BF>100）的识别，进一步证实了全球传播模式的稳定性。

重组形式的多样性与更替模式

基于3D编码区的分析，研究识别出五个主要的重组形式簇（RF-A至RF-E）。系统发育树显示，随着时间的推移，从早期分支逐渐演化出的序列以梯状拓扑结构填充了进化树，呈现出典型的多谱系共循环特征。研究发现，CVA16的重组形式经历了从RF-B和RF-C向RF-D的转变，这一转变伴随着遗传多样性的显著变化。

遗传多样性分析显示，CVA16重组形式的整体遗传多样性在2008年和2016年出现两个高峰，在2012年达到低谷。RF-B在2012年前对遗传多样性起主导作用，而2012年后RF-D接替了这一角色。重组形式的更替周期约为6-8年，这一发现与EV-A71的重组更替周期（5.9-9.4年）相似，但长于埃可病毒9型和30型的周期（1.3-3.1年）。

基因型进化与重组形式更替的协同关联

该研究的核心发现是揭示了CVA16基因型进化与重组形式更替之间的协同关联。通过对比分析VP1基因和3D基因的进化关系，研究发现3D进化树中属于RF-D的序列仅与VP1进化树中的B1b簇聚在一起，RF-E的序列仅与D基因型聚在一起。然而，RF-A的序列与A基因型、B1a或B1b簇均有聚集；RF-B的序列与B2、B1a或B1b簇聚集；RF-C的序列与B1a、B1b或B1c簇聚集。

这种复杂的对应关系揭示了病毒进化的协同策略：从RF-B和RF-C共循环向RF-D主导的转变，伴随着从早期B2基因型向B1a和B1c共循环，最终向B1b主导的转变。简而言之，3D基因的重组形式更替与VP1基因的型别进化呈现同步变化的趋势，而非随机独立的事件。

分子钟分析进一步支持了这一协同进化关系。VP1基因的进化速率为每年每站点3.54×10-3次替换，而3D基因的进化速率略高，为每年每站点4.09×10-3次替换。D基因型的最常见祖先追溯至2002年8月，略早于RF-E的2005年4月；B1b簇的最常见祖先追溯至1998年6月，也早于RF-D的1999年10月。这些时间上的先后顺序暗示了基因型进化可能先于或伴随重组形式的更替。

参考文献：

[1]Han Z ,Wang F ,Xiao J , et al.Synergetic association between coxsackievirus A16 genotype evolution and recombinant form shifts.[J].Virus evolution,2024,10(1):vead080-vead080.DOI:10.1093/VE/VEAD080.