甘肃省甘南州家畜细粒棘球蚴病例特征及其病原遗传多样性分析

目的分析甘肃省甘南州细粒棘球蚴(Echinococcus granulosus)基因型、遗传多样性及群体间遗传差异,为甘南州细粒棘球溯源及制定预防控制策略提供参考。方法收集2016-2022年甘南州玛曲县及碌曲县屠宰场家畜肝脏或肺脏感染的囊型棘球蚴病(cystic echinococcosis,CE)包囊,提取包囊基因组DNA,PCR扩增线粒体NADH脱氢酶2(nad2)基因,测序剪接后进行Blast比对,DnaSPv5软件分析nad2序列的基因型,利用Mega11软件采用邻接法构建系统进化树,再使用PopArt 1.7软件建立细粒棘球蚴nad2基因网络图,最后使用Arlequin 3.5.2.2分析nad2不同群体间的遗传差异。结果PCR结果显示,家畜CE感染率为9.97%(64/642),绵羊(9.14%,48/525)和牛(13.68%,16/117)感染率差异极显著(χ^(2)=175.6,P<0.01)。nad2基因条带大小为882 bp,将测序结果与GenBank中发表的细粒棘球蚴的序列进行比对,相似度均达到99%以上。序列分析结果显示,分为G1(占90.63%,58/64)和G3(占9.37%,6/64)两个基因型,序列相似度可达100%;共有12个单倍型,其中Hap1(占50.00%,32/64)是最主要的单倍型,Hap2(占14.06%,9/64)次之,其中Hap3是和Hap4属于G3基因型(占9.37%,6/64)。系统发生树结果显示,G1和G3基因型处在两个小的分支上,在基因网络图中以Hap1单倍型为中心,其他11个单倍型呈散射状分布;序列单倍型多样性范围为0.7086~0.7706,核苷酸多样性范围为0.00321~0.00452。甘南州两个县细粒棘球蚴种群田岛检验(Tajima's D)均为负值,表明甘南州的细粒棘球蚴种群处于扩张或净化选择的趋势,中性检验(Fu's F)推测两个县的细粒棘球蚴种群也处于扩张趋势,且存在一些特有的基因型,FST结果表明,甘南州细粒棘球蚴不同种群间遗传分化差异不显著(P>0.05)。结论本研究发现甘南州家畜表现出较高的遗传多样性,并处于扩张趋势,存在地区特有的单倍型。共有G1和G3两个基因型,12个单倍型,其中Hap1占50%,其近缘单倍型在全世界均有分布,是甘南州乃至其他主要流行区的主要致病单倍型。