云南省细粒棘球绦虫线粒体co1和nd1基因序列分析

doi:10.12140/j.issn.1000-7423.2023.03.007

摘要/Abstract

摘要：

目的分析云南省细粒棘球蚴线粒体细胞色素氧化酶1（co1）和NADPH脱氢酶第1亚基（nd1）的基因型和序列多态性。方法动物源棘球蚴采自云南省香格里拉市、大关县、洱源县、泸水市、维西县的屠宰场牛、羊和放养猪的肝脏或肺病灶组织分离的包囊，人源棘球蚴采自从剑川县、云龙县、隆阳区和玉龙县医院手术摘除的病灶组织。病原学鉴定后选取细粒棘球蚴，用DNA提取试剂盒提取基因组DNA，PCR扩增co1和nd1基因并测序，通过BLAST比对分析序列一致性和单核苷酸多态性，用MEGA-X软件采用邻接法构建基于co1和nd1的系统进化树。结果共采集62份棘球蚴样品，其中36份经病原学检测确定为细粒棘球蚴。co1基因测序成功32条，其中15条为G1型，17条为G5型，长度分别为824 bp和807 bp，提交至GenBank数据库获得的登录号分别为OP413393~OP413402、OP413498~OP413506和OP420520。nd1基因测序成功34条，其中11条为G1型，23条为G5型，长度分别为882 bp和888 bp，提交GenBank数据库获得的登录号为OP471626~OP471638。序列多态性分析结果显示，co1的G1、G5型序列的变异位点分别占1.94%（16/824）、3.10%（25/807）；nd1的G1、G5型序列的变异位点分别占0.45%（4/882）、2.93%（26/888）。基因型多序列比对结果显示，co1基因共有6条同源序列，其中G1型4条，同源序列的遗传距离分别为0.000 0~0.001 1、0.000 0、0.000 0~0.000 6和0.000 0~0.001 7；G5型2条，同源序列遗传距离分别为0.000 0~0.002 5和0.000 0~0.001 2。nd1基因共有3条同源序列，其中G1型1条，遗传距离为0.000 0~0.001 1；G5型2条，遗传距离为0.000 0和0.000 0~0.001 1。系统进化树分析结果显示，基于co1和nd1基因的系统进化树均形成G1和G5型2个分支，G1型与中国中西部（青海、宁夏、甘肃、四川），不丹、中东地区伊朗和约旦等的基因序列（GenBank登录号AF297617.1、KJ628328.1、EU072106.1、MW138946.1、AB688602.1）的亲缘关系较近；G5型与越南、中国广西、英国、波兰、赞比亚、法国、巴基斯坦、巴西、肯尼亚和纳米比亚等的基因序列（GenBank登录号MW558412.1、MN058591.1、KU378107.1、MZ322608.1、KU743915.1、KU743919.1、MN886291.1、KX010903.1、KU743918.1和KX138068.1）的亲缘关系较近。结论云南省细粒棘球绦虫流行株的基因型为G1和G5型，co1和nd1基因均存在单核苷酸多态性。

关键词: 细粒棘球绦虫, co1基因, nd1基因, 序列分析, 云南省

Abstract:

Objective To analyze the genotype and sequence polymorphisms of mitochondrial cytochrome oxidase 1（co1）and NADPH dehydrogenase subunit 1 (nd1) in Echinococcus granulosus from Yunnan Province. Methods Echinococcus of animal origin was collected from cysts isolated from liver or lung lesions of cattle, sheep and free-range pigs at slaughterhouses in Shangri-La, Daguan, Eryuan, Lushui and Weixi counties of Yunnan Province. Echinococcus of human origin was obtained from the focal tissues surgically removed from hospitals in Jianchuan, Yunlong, Longyang and Yulong Counties, then the E. granulosus was selected after pathogenic identification, and used to extract genomic DNA with a DNA extraction kit, followed by PCR amplification and sequencing of co1 and nd1 genes. The sequence homology and single nucleotide polymorphisms were analyzed using BLAST aligment. Phylogenetic trees based on co1 and nd1 were constructed by the neighbor-joining method using MEGA-X software. Results A total of 62 samples of Echinococcus were collected, 36 of which were pathogenically confirmed as E. granulosus. A total of 32 gene co1 sequences were successfully sequenced, of which 15 were type G1 and 17 were type G5, with a length of 824 bp and 807 bp, respectively. After submitting to GenBank database, the login numbers obtained were OP413393-OP413402, OP413498-OP413506 and OP420520. A total of 34 nd1 genes were sequenced successfully, of which 11 were type G1 and 23 were type G5, with lengths of 882 bp and 888 bp, respectively. The entry numbers obtained by submitting the GenBank database were OP471626-OP471638. The results of sequence polymorphism analysis showed that 1.94% (16/824) mutation sites in the co1 type G1 sequence and 3.10% (25/807) mutation sites in the G5 type sequence were found. In the nd1, 0.45% (4/882) and 2.93% (26/888) mutation sites were found in the type G1 sequence and the G5 type sequence, respectively. The results of multiple genotypic sequence comparison showed that there were 6 homologous sequences of the co1 gene, 4 of which were type G1, the genetic distances of homologous sequences were 0.000 0-0.001 1, 0.000 0-0.000 6, and 0.000 0-0.001 7, respectively. Also included two G5 types. The genetic distances of homologous sequences were 0.000 0-0.002 5 and 0.000 0-0.001 2, respectively. There were 3 homologous sequences in nd1 gene, of which 1 was type G1, and the genetic distance was 0.000 0-0.001 1, 2 of which were type G5, genetic distance was 0.000 0 and 0.000 0-0.001 1. Phylogenetic tree analysis showed that the phylogenetic tree based on co1 and nd1 genes formed two branches of type G1 and type G5, type G1 was closely related to the gene sequences of central and western China (Qinghai, Ningxia, Gansu, Sichuan), Bhutan, Iran and Jordan in the Middle East (GenBank entry number AF297617.1, KJ628328.1, EU072106.1, MW138946.1, AB688602.1). The G5 Type was closely related to Vietnam, China's Guangxi, the United Kingdom, Poland, Zambia, France, Pakistan, Brazil, Kenya and Namibia (GenBank login numbers MW558412.1, MN058591.1, KU378107.1, MZ322608.1, KU743915.1, KU743919.1, MN886291.1, KX010903.1, KU743918.1 and KX138068.1). Conclusion The prevalent isolate of E. granulosus in Yunnan Province are of genotype G1 and G5, and both co1 and nd1 gene reveals single nucleotide polymorphisms.

Key words: Echinococcus granulosus, co1 gene, nd1 gene, Sequence analysis, Yunnan Province

中图分类号:

R383.33

李奔福, 王正青, 徐倩, 字金荣, 严信留, 彭佳, 李建雄, 蔡璇, 吴方伟, 杨亚明. 云南省细粒棘球绦虫线粒体co1和nd1基因序列分析[J]. 中国寄生虫学与寄生虫病杂志, 2023, 41(3): 306-311.

LI Benfu, WANG Zhengqing, XU Qian, ZI Jinrong, YAN Xinliu, PENG Jia, LI Jianxiong, CAI Xuan, WU Fangwei, YANG Yaming. Sequence analysis of mitochondrial co1 and nd1 genes in Echinococcus granulosus in Yunnan Province[J]. Chinese Journal of Parasitology and Parasitic Diseases, 2023, 41(3): 306-311.

图/表 2

参考文献 21

[1]	McManus DP, Gray DJ, Zhang WB, et al. Diagnosis, treatment, and management of echinococcosis[J]. BMJ, 2012, 344: e3866. doi: 10.1136/bmj.e3866
[2]	Grabellus F, Worm K, Schmid KW. Induction of the matrix metalloproteinase-2 activation system in arteries by tensile stress. Involvement of the p38 MAP-kinase pathway[J]. Pathol Res Pract, 2007, 203(3): 135-143. pmid: 17306932
[3]	Wu WP, Wang H, Wang Q, et al. A nationwide sampling survey on echinococcosis in China during 2012—2016[J]. Chin J Parasitol Parasit Dis, 2018, 36(1): 1-14. (in Chinese)
	(伍卫平, 王虎, 王谦, 等. 2012—2016年中国棘球蚴病抽样调查分析[J]. 中国寄生虫学与寄生虫病杂志, 2018, 36(1): 1-14.)
[4]	Li BF, Wu FW, Yan XL, et al. Epidemiological analysis of echinococcosis in Yunnan Province from 2012 to 2017[J]. Chin J Parasitol Parasit Dis, 2019, 37(5): 576-582. (in Chinese)
	(李奔福, 吴方伟, 严信留, 等. 2012—2017年云南省棘球蚴病流行病学分析[J]. 中国寄生虫学与寄生虫病杂志, 2019, 37(5): 576-582.)
[5]	Wei YH, Hu Y, Cao JP. Advances in the genetic polymorphism of Echinococcus granulosus[J]. Chin J Parasitol Parasit Dis, 2019, 37(4): 481-485. (in Chinese)
	(魏玉环, 胡媛, 曹建平. 细粒棘球绦虫基因多态性的研究进展[J]. 中国寄生虫学与寄生虫病杂志, 2019, 37(4): 481-485.)
[6]	Li GD. Study on polymorphism of Taenia asiatica in Qinghai-Tibet Plateau[D]. Xining: Qinghai University, 2017: 1-6. (in Chinese)
	(李国栋. 青藏高原棘球属绦虫基因多态性研究[D]. 西宁: 青海大学, 2017: 1-6.)
[7]	Xue T, Shalayiding•Shalikejiang, Li F, et al. Analysis of the genotype of Echinococcus granulosus isolated from Liaoning and Xinjiang[J]. J Pathog Biol, 2018, 13(6): 601-604. (in Chinese)
	(薛婷, 沙拉依丁•沙力克江, 李飞, 等. 辽宁与新疆两地细粒棘球蚴分离株基因型研究[J]. 中国病原生物学杂志, 2018, 13(6): 601-604.)
[8]	Zeng XM, Wang LY, Wu WP, et al. Cluster analysis of cystic echinococcosis in non Tibetan Plateau regions[J]. Chin J Schisto Contr, 2014, 26(2): 180-183. (in Chinese)
	(曾祥嫚, 王立英, 伍卫平, 等. 我国非青藏高原流行区囊型棘球蚴病聚类分析[J]. 中国血吸虫病防治杂志, 2014, 26(2): 180-183.)
[9]	Yu JF, Gui Z, Wu L, et al. Cloning and sequence analysis of the partial nad1 gene within mitochondrial DNA of Echinococcus granulosus[J]. Chin J Zoonoses, 2020, 36(7): 534-538. (in Chinese)
	(于晶峰, 桂峥, 武琳, 等. 细粒棘球蚴线粒体nad1基因的克隆与序列分析[J]. 中国人兽共患病学报, 2020, 36(7): 534-538.)
[10]	Nakao M, Li TY, Han XM, et al. Genetic polymorphisms of Echinococcus tapeworms in China as determined by mitochondrial and nuclear DNA sequences[J]. Int J Parasitol, 2010, 40(3): 379-385. doi: 10.1016/j.ijpara.2009.09.006
[11]	Cao M, Chen K, Li W, et al. Genetic characterization of human-derived hydatid fluid based on mitochondrial gene sequencing in individuals from northern and Western China[J]. J Helminthol, 2018, 94: e2. doi: 10.1017/S0022149X18000883
[12]	Yang JK, Jia WZ, Jing T, et al. Analysis of genetic variation of Echinococcus granulosus collected from three provinces in China[J]. Chin J Vet Sci Technol, 2004, 34(9): 12-16. (in Chinese)
	(杨俊克, 贾万忠, 景涛, 等. 我国三省区细粒棘球绦虫基因的变异分析[J]. 中国兽医科技, 2004, 34(9): 12-16.)
[13]	Li B. Epidemiological investigation of Echinococcus granulosus and analysis of co1 and nd1 gene sequence[D]. Hohhot: Inner Mongolia Agricultural University, 2013: 1-5. (in Chinese)
	(李滨. 细粒棘球蚴的流行病调查及co1与nd1基因序列分析[D]. 呼和浩特: 内蒙古农业大学, 2013: 1-5.)
[14]	Wu WP. Prevalence and distribution of two types of echinococcosis in China[J]. Chin Animal Health, 2017, 19(7): 7-9. (in Chinese)
	(伍卫平. 我国两型包虫病的流行与分布情况[J]. 中国动物保健, 2017, 19(7): 7-9.)
[15]	Romig T, Ebi D, Wassermann M. Taxonomy and molecular epidemiology of Echinococcus granulosus sensu lato[J]. Vet Parasitol, 2015, 213(3/4): 76-84. doi: 10.1016/j.vetpar.2015.07.035
[16]	Wen H, Vuitton L, Tuxun T, et al. Echinococcosis: advances in the 21st century[J]. Clin Microbiol Rev, 2019, 32(2): e00075-e00018.
[17]	Thompson RCA. The taxonomy, phylogeny and transmission of Echinococcus[J]. Exp Parasitol, 2008, 119(4): 439-446. doi: 10.1016/j.exppara.2008.04.016 pmid: 18539274
[18]	Shang JY, Zhang GJ, He W, et al. Taxonomy and molecular epidemiology of Echinococcus granulosus complex causing cystic echinococcosis[J]. Chin J Parasitol Parasit Dis, 2018, 36(2): 166-173, 177. (in Chinese)
	(尚婧晔, 张光葭, 何伟, 等. 细粒棘球蚴病病原分类学与分子流行病学研究进展[J]. 中国寄生虫学与寄生虫病杂志, 2018, 36(2): 166-173, 177.)
[19]	Zhang TM, Yang D, Zeng ZL, et al. Genetic characterization of human-derived hydatid cysts of Echinococcus granulosus sensu lato in Heilongjiang Province and the first report of G7 genotype of E. canadensis in humans in China[J]. PLoS One, 2014, 9(10): e109059. doi: 10.1371/journal.pone.0109059
[20]	Sharma M, Sehgal R, Fomda BA, et al. Molecular characterization of Echinococcus granulosus cysts in north Indian patients: identification of G1, G3, G5 and G6 genotypes[J]. PLoS Negl Trop Dis, 2013, 7(6): e2262. doi: 10.1371/journal.pntd.0002262
[21]	Shi YL, Wan XL, Wang ZY, et al. Correction to: first description of Echinococcus ortleppi infection in China[J]. Parasit Vectors, 2019, 12: 412. doi: 10.1186/s13071-019-3677-3