《Pathogens》：Genetic Diversity of the BLV env Gene and gp51 Mutations in Genotypes G4 and G7 Circulating in Dairy Cattle in the Novosibirsk Region (Western Siberia, Russia)
Dmitry Baboshko,
Kirill Elfimov,
Polina Achigecheva,
Irina Osipova,
Grigoriy Vlasov,
Oleg Rozhkov,
Boyko Margarita,
Aleksey Totmenin,
Aleksandr Agaphonov and
Natalya Gashnikova
编辑推荐：
为阐明牛白血病病毒（BLV）在俄罗斯新西伯利亚地区的分子流行病学特征，研究人员对该地区6个行政区的奶牛进行了大规模采样。通过针对env基因的测序与系统发育分析，发现当地流行的BLV主要为基因型G4和G7，并鉴定出74个氨基酸突变，包括多个新型及俄罗斯特有的突变位点。该研究为理解BLV的传播、演化及其免疫逃逸机制提供了关键数据，对制定区域性的疾病控制策略具有重要意义。论文发表在《Pathogens》上。
牛白血病病毒（Bovine Leukemia Virus, BLV）是一种致癌性逆转录病毒，是地方性牛白血病（Enzootic Bovine Leukosis, EBL）的病原体。这种病毒在全球范围内传播，给养牛业带来了巨大的经济损失，包括产奶量下降、过早淘汰和贸易限制。在俄罗斯，据估计有超过92万头牛携带BLV，使其成为导致牛死亡的主要原因之一。尽管BLV影响深远，但关于其在俄罗斯，特别是西伯利亚等关键产奶区的遗传多样性和分子流行病学数据却非常有限。新西伯利亚地区是俄罗斯联邦重要的奶牛产区，并被划定为BLV疫情不利的地区，然而，对于该地区流行的BLV基因型、其遗传变异以及可能影响病毒传播和免疫逃逸的关键突变，此前缺乏系统性的研究。填补这一知识空白，对于理解病毒的传播动态、预测其进化趋势以及设计基于证据的防控策略至关重要。
为了深入探究新西伯利亚地区流行的BLV的遗传特征，研究人员在《Pathogens》上发表了一项研究。他们从该地区六个行政区（Kochenyovo, Bolotnoye, Tatarsk, Toguchin, Barabinsk, Dovolnoye）的奶牛中，收集了417份BLV阳性血液样本。研究的核心目标是通过分析编码gp51蛋白的env基因片段，来解析该地区BLV的遗传多样性、系统发育关系以及基因型特异性突变。
本研究主要采用了几个关键技术方法：首先，研究人员从奶牛外周血白细胞中提取基因组DNA。其次，使用商品化和实验室自设计的引物，通过聚合酶链式反应（PCR）扩增出env基因的特定片段（约1410 bp）。然后，对PCR产物进行DNA测序，获得了417条新的env基因序列。接着，他们利用最大似然法（ML）对804 bp的标准片段进行系统发育分析，以确定病毒基因型。此外，通过翻译核苷酸序列并进行多重比对，分析了gp51蛋白的氨基酸突变。最后，运用密码子替代模型（如dN/dS分析）和混合效应进化模型（MEME）评估了env基因所受到的选择压力。
3.1. Phylogenetic Analysis of BLV Strains Circulating in the Novosibirsk Region
（新西伯利亚地区流行BLV毒株的系统发育分析）
通过系统发育分析，研究人员发现所有417个BLV分离株可归为两个明确的基因型：基因型4（G4）和基因型7（G7）。其中，G4是优势基因型，占所有样本的77.9%（325株），在除Toguchin区以外的五个区都占主导地位。而G7基因型占22.1%（92株），仅在Toguchin区占主导地位（75.9%）。分析显示，来自同一农场的病毒株往往在系统发育树上紧密聚类，表明病毒主要在农场内部传播。但也发现了来自相距约250公里的不同地区的G7毒株形成显著聚类，这可能暗示了通过牲畜贸易进行的远距离传播。
3.2. Analysis of Nucleotide and Amino Acid Sequences of the BLV
（BLV核苷酸和氨基酸序列分析）
对env基因片段的核苷酸同源性分析显示，G4和G7基因型内部的序列同源性很高（G4: 96.6%–100%；G7: 97.4%–100%）。值得注意的是，即使来自地理上相隔约600公里的农场，也发现了具有100%核苷酸一致性的病毒株，这提示了病毒近期可能通过共同来源或贸易活动进行传播。
3.3. Mutation Analysis of the BLV Isolate Sequences
（BLV分离株序列的突变分析）
通过对78条独特序列的氨基酸比对，共鉴定出74个氨基酸替换（G4中44个，G7中30个）。研究发现了一些高频突变，例如A15T、S49F、K41R和S221L，在超过97%的G4和G7样本中都存在，可能是在该地区早期传播中固定下来的关键适应性变化。此外，研究还识别出多个基因型特异性或俄罗斯特有的突变：例如，S23F、A40P、R88H和I111T突变在俄罗斯G4毒株中特征性出现；而A258V突变则存在于100%的G7毒株中，但在G4中完全缺失。这些突变大多位于gp51蛋白的关键功能域，如中和结构域（ND1, ND3）、表位A和CD8+ T细胞表位内。
选择压力分析显示，env基因整体上主要受到纯化选择（净化选择）的作用（平均dN/dS = 0.181），表明该区域功能约束很强。然而，MEME分析发现了两个位点（密码子10和111）受到间歇性正向（多样化）选择，提示这些位置的变异可能受到了宿主免疫压力的驱动。
研究的结论和讨论部分指出，新西伯利亚地区流行的BLV基因型G4和G7，与欧洲及前苏联国家（如波兰、比利时、哈萨克斯坦、摩尔多瓦）的毒株在遗传上关系密切，这支持了BLV主要通过种畜进口和牲畜贸易从欧洲引入俄罗斯，随后在国内及中亚地区传播的假说。G4和G7在地区内分布不均，以及农场内传播链的证据，揭示了病毒传播的局部性和复杂性。
本研究鉴定出的众多氨基酸突变，特别是集中在免疫原性关键区域gp51结构域内的突变，可能代表了病毒为逃避免疫压力而进行的适应性进化。例如，G7特有的A258V突变位于表位A，而G4特有的I111T突变位于中和结构域ND1，这些都可能影响病毒的抗原性和与宿主的相互作用。虽然大部分突变处于纯化选择下，但个别位点（如10和111）检测到的正向选择信号，暗示了局部适应性进化的发生。
这项研究系统地描绘了俄罗斯新西伯利亚地区BLV的分子流行病学图景，明确了G4和G7为主要流行基因型，并绘制了其地理分布图谱。通过鉴定大量新型和基因型特异性突变，特别是那些位于关键免疫表位区域的突变，为理解BLV的进化动力学和潜在免疫逃逸机制提供了重要线索。研究人员将417条新的env基因序列提交至GenBank数据库，极大地丰富了BLV的全球遗传资源库。这些发现不仅为开发针对流行毒株的诊断方法和疫苗候选物奠定了分子基础，也为制定基于当地病毒遗传特征的、更有效的牛白血病区域防控策略提供了科学依据，对于保障畜牧业健康和经济安全具有重要意义。
