新疆新源县和四川石渠县啮齿动物及家畜棘球绦虫线粒体cox1基因遗传多态性分析

目的了解新疆新源县和四川石渠县棘球绦虫线粒体细胞色素。氧化酶亚基1 (cox1)基因的遗传多态性,探讨棘球绦虫在我国的传播动力和传播方向。方法 2014-2018年,在新疆新源县、四川石渠县捕捉啮齿类动物并采集其肝组织样品,收集家畜病变脏器。提取肝组织和病变脏器组织的DNA,PCR扩增棘球绦虫cox1基因片段,经分子克隆、测序后,在NCBI数据库中比对以确定所感染的虫种。使用Clustal X2和MEGA 7剪辑序列,使用DnaSP v5和Arlequin 3.5计算核苷酸多态性和中性检验,分析新源县和石渠县棘球绦虫的遗传多样性。以石渠棘球绦虫为外群,用MrBayes 3.2.4建立基于cox1基因的贝叶斯系统进化树,总结分析虫种的基因型。采用Network 5.0绘制棘球绦虫cox1基因的单倍型网络图,分析单倍型结构。结果在新疆新源县捕获普通田鼠122只,收集绵羊病变脏器5个;在四川石渠县捕获青海田鼠144只、经营田鼠44只和高原鼠兔135只,收集牦牛病变脏器6个。其中从5只绵羊、4头牦牛的病变脏器样品DNA扩增获得细粒棘球绦虫cox1基因片段40条(新源县26条,石渠县14条),从8只普通田鼠、14只青海田鼠、5只经营田鼠和2只高原鼠兔获得多房棘球绦虫线粒体cox1基因59条(新源县20条,石渠县39条),扩增产物长度均为875 bp。细粒棘球绦虫在两地的遗传分化程度(Fst=0.088 63)高于多房棘球绦虫(Fst=0.000 88),新源县细粒棘球绦虫遗传多样性(0.002 58±0.000 46)低于石渠县(0.005 88±0.000 58),而其多房棘球绦虫遗传多样性(0.002 28±0.000 46)高于石渠县(0.001 37±0.000 30)。贝叶斯系统进化树显示,两地的细粒棘球绦虫基因型为G1型,多房棘球绦虫为亚洲型。序列比对发现25个细粒棘球绦虫cox1基因单倍型(新源县15个,石渠县9个,共同单倍型1个)和29个多房棘球绦虫cox1基因单倍型(新源县13个,石渠县15个,共同单倍型1个),两种棘球绦虫的单倍型网络图均为环绕主单倍型(共同单倍型)的辐射状结构,新源县的细粒棘球绦虫序列中有9条为主单倍型(34.6%,9/26),而石渠县仅有1条(1/14),显示出新源县在细粒棘球绦虫单倍型结构中的核心位置;石渠县的多房棘球绦虫基因序列中有23条为主单倍型(58.9%,23/39),而新源县为7条(35.0%,7/20),表明石渠县在多房棘球绦虫的单倍型结构中更具有主导地位。结论新源县和石渠县的棘球绦虫线粒体cox1基因型一致,但遗传多态性存在一定差异,同等空间跨度上细粒棘球绦虫的遗传分化程度要高于多房棘球绦虫,这些差异为探讨棘球绦虫传播方向和动力提供了重要参考依据。

Genetic polymorphisms of mitochondrial cox1 gene of Echinococcus spp. from rodents and livestock in Xinyuan County,Xinjiang and Shiqu County,Sichuan

Objective To understand the genetic polymorphisms of cytochrome c oxidase subunit 1(cox1)gene in mitochondria of Echinococcus spp.in Xinyuan County,Xinjiang and Shiqu County,Sichuan,and discuss Echinococcus transmission directions and dynamics in China.Methods From 2014 to 2018,rodents were captured in Xinyuan County(Xinjiang) and Shiqu County(Sichuan) to collect liver tissues,meanwhile lesioned tissue samples of livestock were collected.DNA was extracted from the livers and lesioned tissues,and underwent PCR to amplify part of Echinococcus mitochondrial coxl gene.PCR products were sequenced,and the sequencing results were aligned in NCBI database to identify the Echinococcus species.The sequences were edited by Clustal X2 and MAGE7 softwares,and diversity and neutrality indices were calculated by DnaSP v5 and Arlequin 3.5.Genetic diversities were compared between the two regions.By using E.shiquicus as an outgroup,coxl sequence-based Bayesian phylogenetic trees were established with MrBayes 3.2.4 software.The haplotype networks were drawn for coxl genes of Echinococcus spp.using Network 5.0 to analyze haplotype structure.Results One hundred and twenty-two Microtus arvalis were captured and 5 lesioned tissues of sheep were collected in Xinyuan County(Xinjiang);while in Shiqu County(Sichuan),144 Lasiop odomys fuscus,44 M.limnophilus and 135 Ochotona c urzoniae were captured and4 lesioned tissues of yaks were collected.A total of 40 sequences of E.granulosus coxl genes(26 from Xinyuan,and 14 from Shiqu) were amplified from lesioned tissues of 5 sheeps and 4 yaks.Fifty-nine sequences of E.multilocularis mitochondrial coxl genes(20 from Xinyuan,and 39 from Shiqu) were amplified in 8 M.arvalis,14 L.fuscus,5 M.limnophilus and 2 0.c urzoniae.All PCR products had a length of 875 bp.E.granulosus showed a higher degree of genetic differentiation(Fst=0.088 63) between the 2 regions than E.multiloc ularis(Fst=0.000 88).The genetic diversity of E.granulosus was lower in Xinyuan(0.002 58±0.000 46) than in Shiqu(0.005 88±0.000 58),while E.muitilocularis had higher genetic diversity in Xinyuan(0.002 28±0.000 46) than in Shiqu(0.001 37±0.000 30).The Bayesian phylogenetic tree showed that in both regions E.granulosus had the G1 genotype,and E.muitilocularis was of the Asian type.Sequences alignment revealed 25 haplotypes for E.granulosus coxl gene(15 in Xinyuan,9 in Shiqu,and 1 shared by both) and 29 haplotypes for E.muitilocularis(13 in Xinyuan,15 in Shiqu,1 shared by both).The haplotype networks of E.multiloc ularis and E.grauulosus both showed a star-shaped network with a centrally placed main haplotype.There were 9 sequences of E.grauulosus being the main haplotype type in Xinyuan(34.6%,9/26),while only 1 sequence of that type was found in Shiqu(1/14),indicating a central role of Xinyuan County in the E.grauulosus haplotype network.On the contrary,there were 23 sequences of E.multilocularis being the main haplotype type in Shiqu(58.9%,23/39) and 7 of that type in Xinyuan(35.0%,7/20),indicating that Shiqu County was more predominant in the E.multiloc ularis haplotype structure.Conclusion There is no apparent difference in the mitochondrial coxl genotype of Echinococcus found between Xinyuan and Shiqu County,but there are some differences in the genetic polymorphism detected in the two sites.The genetic differentiation level of E.granulosus is higher than that of E.multilocularis over the same spatial span.These findings provide important referential basis for further study on transmission direction and dynamics of Echinococcus.