青海地区多房棘球绦虫和石渠棘球绦虫遗传分化特征

doi:10.12140/j.issn.1000-7423.2024.03.005

中国寄生虫学与寄生虫病杂志 ›› 2024, Vol. 42 ›› Issue (3): 309-315.doi: 10.12140/j.issn.1000-7423.2024.03.005

青海地区多房棘球绦虫和石渠棘球绦虫遗传分化特征

付永¹^,²(), 张海宁¹^,², 陈旺开¹^,², 师正合³, 张学勇¹^,², 郭志宏¹^,², 朵红¹^,², 沈秀英¹^,², 孟茹⁴, 李志¹^,²^,^*()

1 青海大学畜牧兽医科学院，西宁 810016
2 青海省动物疾病病原诊断与绿色防控技术研究重点实验室，西宁 810016
3 青海大学附属医院，西宁 810001
4 西宁市动物疫病预防控制中心，青海西宁 810003

收稿日期:2023-12-04 修回日期:2024-01-21 出版日期:2024-06-30 发布日期:2024-07-16
通讯作者: *李志（1987—），男，博士，助理研究员，从事人兽共患病致病机制及防控技术研究。E-mail：lizhi19880717@163.com
作者简介:付永（1980—），男，博士，副研究员，从事棘球蚴病致病机制及防控技术研究。E-mail：qhfuyong@163.com
基金资助:
国家自然科学基金(32160840);青海省“昆仑英才?高端创新创业人才”项目(2022)

Characteristics of genetic differentiation of Echinococcus multilocularis and E. shiquicus in Qinghai region

FU Yong¹^,²(), ZHANG Haining¹^,², CHEN Wangkai¹^,², SHI Zhenghe³, ZHANG Xueyong¹^,², GUO Zhihong¹^,², DUO Hong¹^,², SHEN Xiuying¹^,², MENG Ru⁴, LI Zhi¹^,²^,^*()

1 Academy of Animal Science and Veterinary Medicine, Qinghai University, Xining 810016, China
2 Qinghai Provincial Key Laboratory of Pathogen Diagnosis for Animal Diseases and Green Technical Research for Prevention and Control,Xining 810016, China
3 Affiliated Hospital of Qinghai University, Xining 810001, China
4 Xining Animal Disease Control Center, Xining 810003, Qinghai, China

Received:2023-12-04 Revised:2024-01-21 Online:2024-06-30 Published:2024-07-16
Supported by:
National Natural Science Foundation of China(32160840);Qinghai Province “Kunlun Talents High-end Innovation and Entrepreneurial Talents” Top Talent Proiect(2022)

摘要/Abstract

摘要：

目的对青海地区多房棘球绦虫和石渠棘球绦虫遗传分化特征进行分析，为青海省棘球蚴病的防控提供理论依据。方法在玉树、果洛、黄南藏族自治州等棘球绦虫主要自然疫源地捕捉小型哺乳动物，剖检采集包囊，提取包囊组织基因组DNA，PCR扩增细胞色素氧化酶1（cox1）基因并测序，运用DnaSP v6、Iqtree、BEAST v2.7.4等软件进行单倍型分析、核苷酸多态性分析、系统进化树构建及棘球绦虫属分化时间估算。结果从2 864只小型哺乳动物中共获得55份棘球蚴包囊。55份包囊DNA样品均扩增出长度约800 bp的cox1条带，其中37份为多房棘球绦虫，18份为石渠棘球绦虫，青海田鼠多房棘球蚴患病率为1.96%（37/1 884），高原鼠兔石渠棘球蚴患病率为1.84%（18/980）。37条多房棘球绦虫cox1序列共有单倍型5个，以EmH3单倍型为主（占33/37），单倍体多样性指数为0.207，核苷酸多样性指数为0.033 55，共有156个变异位点。18条石渠棘球绦虫cox1序列共有单倍型8个，以EsH2单倍型为主（占8/18），单倍型多样性指数为0.778，核苷酸多样性指数为0.060 52；共有14个变异位点。多房棘球绦虫和石渠棘球绦虫的13个单倍型上传GenBank，单倍型EmH1~EmH5的登录号依次为OR821706、OR821707、OR830343、OR830344、OR826123，EsH1~EsH8登录号依次为OR835156、OR835157、OR830376、OR830378、OR831110、OR875250、OR835161、OR841080。系统进化树分析显示，多房棘球绦虫的5个单倍型与多房棘球绦虫亚洲株聚为一支，石渠棘球绦虫的8个单倍型与GenBank中的石渠棘球绦虫聚为一支。基于cox1基因的分化时间结果显示，细粒棘球绦虫、多房棘球绦虫、石渠棘球绦虫、少节棘球绦虫和伏氏棘球绦虫的共同祖先存在于约5.67百万年前（Mya）（95% CI：4.72~6.66 Mya），细粒棘球绦虫、石渠棘球绦虫和多房棘球绦虫的平均分化时间约为2.02 Mya（95% CI：1.51~2.49 Mya）。结论青海地区多房棘球绦虫和石渠棘球绦虫具有较高的遗传多样性，其中多房棘球绦虫以EmH3单倍型为主，石渠棘球绦虫以EsH2单倍型为主。

关键词: 多房棘球绦虫, 石渠棘球绦虫, 单倍型, cox1基因, 系统发育分析, 分化时间

Abstract:

Objective To analyze the genetic differentiation characteristics of Echinococcus multilocularis and E. shiquicus in Qinghai region to provide theoretical support for the prevention and control of echinococcosis in Qinghai Province. Methods Small mammals were captured in the main natural endemic areas of Echinococcus spp. in Yushu, Guoluo, and Huangnan Tibetan Autonomous Prefecture and were dissected to collect cysts. The genomic DNA from cysts tissue was extracted and the cytochrome oxidase 1 (cox1) gene was amplified using PCR and sequenced. DnaSP v6, Iqtree, BEAST v2.7.4 and other software were used for haplotype analysis, nucleotide polymorphism analysis, construction of a phylogenetic tree, and estimation of the divergence time of the Echinococcus genus. Results A total of 55 hydatid cysts were obtained from 2 864 small mammals. All 55 cyst samples were amplified for cox1 bands with a length of approximately 800 bp, of which 37 were E. multilocularis, and 18 were E. shiquicus, the prevalence of E. multilocularis in Neodon fuscus was 1.96% (37/1884). The prevalence of E. shiquicus in Ochotona curzoniae was 1.84% (18/980). In the 37 cox1 sequences of E. multilocularis, there were 5 haplotypes in the 37 cox1 sequences of E. multilocularis with EmH3 being the predominant one (33/37), the haplotype diversity index was 0.207, the nucleotide diversity index was 0.033 55, and there were 156 variable sites. In the 18 cox1 sequences of E. shiquicus, There were 8 haplotypes, with the EsH2 haplotype being the predominant one (8/18), the haplotype diversity index was 0.778, the nucleotide diversity index was 0.060 52, and there were 14 variable sites. Thirteen haplotypes of E. multilocularis and E. shiquensis were uploaded to GenBank. The accession numbers of haplotypes EmH1-EmH5 are OR821706, OR821707, OR830343, OR830344, OR826123, respectively. The accession numbers of haplotypes ESH1-ESH8 are OR835156, OR835157, OR830376, OR830378, OR831110, OR875250, OR835161, OR841080. The phylogenetic tree shows that the 5 haplotypes of E. multilocularis were clustered together with the Asian strain of E. multilocularis, and the 8 haplotypes of E. shiquicus were clustered with E. shiquicus in the GenBank. The divergence time based on the cox1 gene showed that the common ancestor of E. granulosus, E. multilocularis, E. shiquicus, E. oligarthrus and E. vogeli existed approximately 5.67 million years ago (Mya) (95% CI: 4.72-6.66 Mya), and the average divergence time for E. granulosus, E. shiquicus and E. multilocularis was approximately 2.02 Mya (95% CI: 1.51-2.49 Mya). Conclusion E. multilocularis and E. shiquicus in Qinghai region have high genetic diversity, with EmH3 haplotype dominating E. multilocularis and EsH2 haplotype dominating E. shiquicus.

Key words: Echinococcus multilocularis, Echinococcus shiquicus, Haplotype, cox1 gene, Phylogenetic analysis, Divergence time

中图分类号:

R383.33

付永, 张海宁, 陈旺开, 师正合, 张学勇, 郭志宏, 朵红, 沈秀英, 孟茹, 李志. 青海地区多房棘球绦虫和石渠棘球绦虫遗传分化特征[J]. 中国寄生虫学与寄生虫病杂志, 2024, 42(3): 309-315.

FU Yong, ZHANG Haining, CHEN Wangkai, SHI Zhenghe, ZHANG Xueyong, GUO Zhihong, DUO Hong, SHEN Xiuying, MENG Ru, LI Zhi. Characteristics of genetic differentiation of Echinococcus multilocularis and E. shiquicus in Qinghai region[J]. Chinese Journal of Parasitology and Parasitic Diseases, 2024, 42(3): 309-315.

图/表 4

参考文献 39

[1]	Casulli A, Barth TFE, Tamarozzi F. Echinococcus multilocularis[J]. Trends Parasitol, 2019, 35(9): 738-739. doi: S1471-4922(19)30109-6 pmid: 31182385
[2]	Gong QL, Ge GY, Wang Q, et al. Meta-analysis of the prevalence of Echinococcus in dogs in China from 2010 to 2019[J]. PLoS Negl Trop Dis, 2021, 15(4): e0009268.
[3]	Torgerson PR, Keller K, Magnotta M, et al. The global burden of alveolar echinococcosis[J]. PLoS Negl Trop Dis, 2010, 4(6): e722.
[4]	Deplazes P, Rinaldi L, Alvarez Rojas CA, et al. Global distribution of alveolar and cystic echinococcosis[J]. Adv Parasitol, 2017, 95: 315-493.
[5]	Craig PS, Giraudoux P, Wang ZH, et al. Echinococcosis transmission on the Tibetan Plateau[J]. Adv Parasitol, 2019, 104: 165-246.
[6]	Mathy A, Hanosset R, Adant S, et al. The carriage of larval Echinococcus multilocularis and other cestodes by the musk rat (Ondatra zibethicus) along the Ourthe River and its tributaries (Belgium)[J]. J Wildl Dis, 2009, 45(2): 279-287.
[7]	Beiromvand M, Akhlaghi L, Fattahi Massom SH, et al. Molecular identification of Echinococcus multilocularis infection in small mammals from Northeast, Iran[J]. PLoS Negl Trop Dis, 2013, 7(7): e2313.
[8]	Kui Y, Xue CZ, Wang X, et al. Progress of echinococcosis control in China, 2022[J]. Chin J Parasitol Parasit Dis, 2024, 42(1): 8-16. (in Chinese)
	(蒉嫣, 薛垂召, 王旭, 等. 2022年全国棘球蚴病防治工作进展[J]. 中国寄生虫学与寄生虫病杂志, 2024, 42(1): 8-16.) doi: 10.12140/j.issn.1000-7423.2024.01.002
[9]	Kui Y, Xue CZ, Wang X, et al. Progress of echinococcosis control in China, 2021[J]. Chin J Parasitol Parasit Dis, 2023, 41(2): 142-148. (in Chinese)
	(蒉嫣, 薛垂召, 王旭, 等. 2021年全国棘球蚴病防治进展[J]. 中国寄生虫学与寄生虫病杂志, 2023, 41(2): 142-148.) doi: 10.12140/j.issn.1000-7423.2023.02.003
[10]	Oksanen A, Siles-Lucas M, Karamon J, et al. The geographical distribution and prevalence of Echinococcus multilocularis in animals in the European Union and adjacent countries: a systematic review and meta-analysis[J]. Parasit Vectors, 2016, 9(1): 519.
[11]	Nakao M, Xiao N, Okamoto M, et al. Geographic pattern of genetic variation in the fox tapeworm Echinococcus multilocularis[J]. Parasitol Int, 2009, 58(4): 384-389.
[12]	Dán Á, Rónai Z, Széll Z, et al. Prevalence and genetic characterization of Echinococcus spp. in cattle, sheep, and swine in Hungary[J]. Parasitol Res, 2018, 117(9): 3019-3022. doi: 10.1007/s00436-018-5977-5 pmid: 29934692
[13]	Xiao N, Qiu JM, Nakao M, et al. Echinococcus shiquicus, a new species from the Qinghai-Tibet Plateau Region of China: discovery and epidemiological implications[J]. Parasitol Int, 2006, 55(Suppl): S233-S236.
[14]	Yan HB, Li L, Li WH, et al. Echinococcus shiquicus in Qinghai-Tibet Plateau: population structure and confirmation of additional endemic areas[J]. Parasitology, 2021, 148(7): 879-886.
[15]	Zhu GQ, Yan HB, Li L, et al. First report on the phylogenetic relationship, genetic variation of Echinococcus shiquicus isolates in Tibet Autonomous Region, China[J]. Parasit Vectors, 2020, 13(1): 590.
[16]	Wu YT, Li L, Xu FL, et al. Establishment of a secondary infection laboratory model of Echinococcus shiquicus metacestode using BALB/c mice and Mongolian jirds (Meriones unguiculatus)[J]. Parasitology, 2023, 150(9): 813-820.
[17]	Zhang YG, Ma YY, Cao DP, et al. Bioinformatic analyses on sequences of the complete mitochondrial genomes of Echinococcus genus[J]. Chin J Zoonoses, 2019, 35(3): 271-277. (in Chinese)
	(张耀刚, 马艳艳, 曹得萍, 等. 棘球属绦虫线粒体基因组全序列生物信息学分析[J]. 中国人兽共患病学报, 2019, 35(3): 271-277.)
[18]	Han XM, Jian YN, Zhang XY, et al. Genetic characterization of Echinococcus isolates from various intermediate hosts in the Qinghai-Tibetan Plateau Area, China[J]. Parasitology, 2019, 146(10): 1305-1312.
[19]	Wang LY. Study on the relation of distribution of intermediate host of Echinococcus multilocularis to the environmental factors in southern Qinghai Plateau[D]. Beijing: Chinese Center for Disease Control and Prevention, 2009: 6-13. (in Chinese)
	(王立英. 青南高原多房棘球绦虫中间宿主分布与环境因素关系的研究[D]. 北京: 中国疾病预防控制中心, 2009: 6-13.)
[20]	Jiang WB, Liu N, Zhang GT, et al. Specific detection of Echinococcus spp. from the Tibetan fox (Vulpes ferrilata) and the red fox (V. Vulpes) using copro-DNA PCR analysis[J]. Parasitol Res, 2012, 111(4): 1531-1539.
[21]	Rozas J, Ferrer-Mata A, Sánchez-DelBarrio JC, et al. DnaSP 6: DNA sequence polymorphism analysis of large data sets[J]. Mol Biol Evol, 2017, 34(12): 3299-3302. doi: 10.1093/molbev/msx248 pmid: 29029172
[22]	Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability[J]. Mol Biol Evol, 2013, 30(4): 772-780. doi: 10.1093/molbev/mst010 pmid: 23329690
[23]	Capella-Gutiérrez S, Silla-Martínez JM, Gabaldón T. trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses[J]. Bioinformatics, 2009, 25(15): 1972-1973. doi: 10.1093/bioinformatics/btp348 pmid: 19505945
[24]	Nguyen LT, Schmidt HA, von Haeseler A, et al. IQ-TREE: a fast and effective stochastic algorithm for estimating maximum likelihood phylogenies[J]. Mol Biol Evol, 2015, 32(1): 268-274.
[25]	Kalyaanamoorthy S, Minh BQ, Wong TKF, et al. ModelFinder: fast model selection for accurate phylogenetic estimates[J]. Nat Methods, 2017, 14(6): 587-589. doi: 10.1038/nmeth.4285 pmid: 28481363
[26]	Xie JM, Chen YR, Cai GJ, et al. Tree visualization by one table (tvBOT): a web application for visualizing, modifying and annotating phylogenetic trees[J]. Nucleic Acids Res, 2023, 51(W1): W587-W592.
[27]	Bouckaert R, Heled J, Kühnert D, et al. BEAST 2: a software platform for Bayesian evolutionary analysis[J]. PLoS Comput Biol, 2014, 10(4): e1003537.
[28]	Knapp J, Nakao M, Yanagida T, et al. Phylogenetic relationships within Echinococcus and Taenia tapeworms (Cestoda ∶ Taeniidae): an inference from nuclear protein-coding genes[J]. Mol Phylogenet Evol, 2011, 61(3): 628-638.
[29]	Wang X, Liu JY, Zuo QQ, et al. Echinococcus multilocularis and Echinococcus shiquicus in a small mammal community on the eastern Tibetan Plateau: host species composition, molecular prevalence, and epidemiological implications[J]. Parasit Vectors, 2018, 11(1): 302.
[30]	Li TY, Chen XW, Zhen R, et al. Widespread co-endemicity of human cystic and alveolar echinococcosis on the eastern Tibetan Plateau, northwest Sichuan/southeast Qinghai, China[J]. Acta Trop, 2010, 113(3): 248-256. doi: 10.1016/j.actatropica.2009.11.006 pmid: 19941830
[31]	Ma JY, Wang H, Lin GH, et al. Molecular identification of Echinococcus species from eastern and southern Qinghai, China, based on the mitochondrial cox1 gene[J]. Parasitol Res, 2012, 111(1): 179-184.
[32]	Wu CC, Zhang WB, Ran B, et al. Genetic variation of mitochondrial genes among Echinococcus multilocularis isolates collected in Western China[J]. Parasit Vectors, 2017, 10(1): 265.
[33]	Li CY, Guan YY, Wu WP, et al. Progress of researches on infection with two species of Echinococcus causing human diseases in animal hosts and influencing factors[J]. Chin J Schisto Control, 2022, 34(2): 194-199. (in Chinese)
	(李春阳, 官亚宜, 伍卫平, 等. 两种致病人体棘球绦虫动物宿主感染及影响因素研究进展[J]. 中国血吸虫病防治杂志, 2022, 34(2): 194-199.)
[34]	Fu MH, Wang X, Han S, et al. Advances in research on echinococcoses epidemiology in China[J]. Acta Trop, 2021, 219: 105921.
[35]	Zhou NN, Wang MX, Cui L, et al. Genetic variation of Empoasca vitis (Göthe) (Hemiptera ∶ Cicadellidae) among different geographical populations based on mtDNA COⅠ complete sequence[J]. Acta Ecol Sin, 2014, 34(23): 6879-6889. (in Chinese)
	(周宁宁, 王梦馨, 崔林, 等. 基于COⅠ基因全长序列的假眼小绿叶蝉地理种群遗传分化研究[J]. 生态学报, 2014, 34(23): 6879-6889.)
[36]	Tomiya S, Tseng ZJ. Whence the beardogs reappraisal of the middle to late Eocene ‘Miacis’ from texas, usa, and the origin of Amphicyonidae (Mammalia, Carnivora)[J]. R Soc Open Sci, 2016, 3(10): 160518.
[37]	Clark PU, Dyke AS, Shakun JD, et al. The last glacial maximum[J]. Science, 2009, 325(5941): 710-714. doi: 10.1126/science.1172873 pmid: 19661421
[38]	Wu YB, Liu YG, Yi CL, et al. Impact of Tibetan glacier change on the Asian climate during the last glacial maximum[J]. Acta Sci Nat Univ Pekin, 2019, 55(1): 159-170. (in Chinese)
	(毋宇斌, 刘永岗, 易朝路, 等. 末次冰盛期青藏高原冰川变化对亚洲气候的影响[J]. 北京大学学报(自然科学版), 2019, 55(1): 159-170.)
[39]	Wen H, Vuitton L, Tuxun T, et al. Echinococcosis: advances in the 21st century[J]. Clin Microbiol Rev, 2019, 32(2): e00075-e00018.

青海地区多房棘球绦虫和石渠棘球绦虫遗传分化特征

Characteristics of genetic differentiation of Echinococcus multilocularis and E. shiquicus in Qinghai region

RichHTML

PDF (PC)

可视化

摘要/Abstract

引用本文

使用本文

图/表 4

参考文献 39

相关文章 15

编辑推荐

Metrics

本文评价

[1]	达哇卓玛, 刘川川, 樊海宁. 细胞程序性死亡在棘球蚴病中的研究进展[J]. 中国寄生虫学与寄生虫病杂志, 2024, 42(2): 259-266.
[2]	郭璐, 吴小霞, 段兰利, 王冰洁, 徐宁, 阿热艾·阿哈太, 乌云花, 赵莉, 班万里, 陈云英, 余婉蓉, 刘帅, 潘星羽, 吾力江·卡马力, 徐晶, 穆尼拉·特列吾汗, 张壮志. 新疆部分地区牛、羊细粒棘球蚴感染情况调查及基因多态性分析[J]. 中国寄生虫学与寄生虫病杂志, 2022, 40(5): 603-609.
[3]	毋德芳, 付永, 任宾, 张耀刚, 许晓磊, 庞明泉, 樊海宁. 青海人源两型棘球绦虫（棘球蚴）遗传多样性及分化时间研究[J]. 中国寄生虫学与寄生虫病杂志, 2022, 40(5): 610-615.
[4]	仲顺虎, 孙玥, 郭小腊, 郑亚东, 陈轶霞. 多房棘球蚴感染小鼠脾淋巴细胞中差异表达miRNA的鉴定及其生物信息学分析[J]. 中国寄生虫学与寄生虫病杂志, 2022, 40(3): 288-294.
[5]	卓怡呈, 杨海成, 刘程豪, 张宝财, 多小勇, 张示杰. 骨桥蛋白表达水平对多房棘球蚴原头节生长发育的影响[J]. 中国寄生虫学与寄生虫病杂志, 2022, 40(3): 299-304.
[6]	史琦琪, 刘丛珊, 霍乐乐, 魏玉芬, 姜斌, 殷梦, 薛剑, 陶奕, 张皓冰. 氨基醇类化合物HT24对多房棘球蚴原头节微管蛋白表达水平的影响[J]. 中国寄生虫学与寄生虫病杂志, 2021, 39(4): 437-443.
[7]	郭宝平, 郭刚, 张莉, 向晶晶, 王小平, 任远, 齐文静, 张慧, 李军, 张文宝, 汪海燕. 新疆察布查尔县小型啮齿类动物多房棘球蚴感染情况调查[J]. 中国寄生虫学与寄生虫病杂志, 2021, 39(3): 327-332.
[8]	徐凯, 王海久, 张丽, 张耀刚, 樊海宁, 任利, 任宾, 王志鑫. 多房棘球蚴对宿主肝细胞损害机制的研究进展[J]. 中国寄生虫学与寄生虫病杂志, 2021, 39(2): 256-259.
[9]	尚婧晔, 张光葭, 喻文杰, 何伟, 廖沙, 李汭芮, 黄燕, 刘阳, 钟波. 多房棘球绦虫遗传多态性研究进展[J]. 中国寄生虫学与寄生虫病杂志, 2020, 38(5): 637-641.
[10]	周鸿让, 茅光耀, 王晓玲, 陈木新, 余晴, 王莹, 艾琳, 肖宁. 重组酶介导的多重核酸等温扩增法鉴别细粒棘球绦虫和多房棘球绦虫技术的建立与应用[J]. 中国寄生虫学与寄生虫病杂志, 2020, 38(3): 310-316.
[11]	魏玉环, 刘华, 李武军, 赵海, 胡媛, 曹建平. 西藏阿里地区细粒棘球蚴人体分离株nad1基因多态性分析[J]. 中国寄生虫学与寄生虫病杂志, 2020, 38(1): 17-22.
[12]	苏梦1,2，郭小腊2，杨静2，邵忠伟2，丁军涛3，项海涛1，骆学农2,4，郑亚东2,4，孙晓林1*. 多房棘球绦虫apomucin基因qPCR引物的筛选及潜在应用[J]. 中国寄生虫学与寄生虫病杂志, 2016, 34(5): 2-394-398.
[13]	任宾，樊海宁*，邓勇，王海久，任利. 泡球蚴囊液对大鼠肝星状细胞MAPK信号通路5个相关基因的影响[J]. 中国寄生虫学与寄生虫病杂志, 2015, 33(2): 8-114-117,121.
[14]	倪兴维，闫鸿斌，娄忠子，贾万忠*. 多房棘球绦虫信号转导体系[J]. 中国寄生虫学与寄生虫病杂志, 2012, 30(3): 15-233-237.
[15]	包永星，张月芬，倪雅琼，谢增如，齐洪志，毛睿，杨玉刚，温浩*. 电子线治疗大鼠继发性泡球蚴病的实验研究[J]. 中国寄生虫学与寄生虫病杂志, 2011, 29(2): 9-127-129.