畜牧兽医学报  2021, Vol. 52 Issue (12): 3546-3556. DOI: 10.11843/j.issn.0366-6964.2021.012.021    PDF    
基于粪便样本的野生犬科动物物种鉴别与棘球绦虫感染调查
尚婧晔1, 张光葭1, 丹巴泽里2, 王奇1, 喻文杰1, 何伟1, 廖沙1, 沈芸舟3, 黄燕1, 王谦1, 钟波1, 刘阳1     
1. 四川省疾病预防控制中心, 成都 610041;
2. 甘孜藏族自治州州疾病预防控制中心, 康定 626099;
3. 凉山彝族自治州疾病预防控制中心, 西昌 615000
摘要:旨在探索粪便样本分析对犬科动物物种鉴别的方法及其可行性,同时对野外犬科动物棘球绦虫感染情况进行调查。采用哺乳动物线粒体DNA控制区通用引物对101份采集自野外环境的犬科动物粪便DNA进行PCR扩增和测序,通过核苷酸位点变异分析、BLAST相似性比对、遗传距离分析和系统进化树构建,对粪样来源物种进行分子鉴别。同时,采用粪-抗原ELISA检测,对粪样棘球绦虫感染情况进行调查。结果表明:粪样来源物种鉴别的检测成功率为73.27%,获得的74条序列共定义单倍型20个,所有单倍型序列均可在GenBank数据库中匹配到单一物种,序列相似度为98.52%~100%。单倍型根据序列差异可分为4个单倍型集,各单倍型集间平均碱基差异为17.83~70.10,各单倍型集内单倍型碱基差异在1~10个之间。各集合间遗传距离为0.068~0.342,远大于各集合内单倍型间遗传距离0.009~0.022。系统进化分析结果显示,同一集合的单倍型以高自展值聚集成一支。综合各项分析结果,可推测所检测粪样宿主来源分别为犬、灰狼、赤狐和红狐。粪-抗原ELISA检测发现阳性样本11份,样本棘球绦虫抗原阳性率为10.89%。综上表明:线粒体DNA控制区序列分析可用于犬科动物的分子鉴别,可结合核基因等其他分子标记进一步提升检测准确率。研究区域野外犬科动物棘球绦虫粪抗原阳性率较高。
关键词犬科动物    棘球绦虫    终末宿主    线粒体基因    物种鉴定    
Identification of Wild Canidae Species and Investigation of Echinococcus Infection Based on Fecal Samples
SHANG Jingye1, ZHANG Guangjia1, DANBA Zeli2, WANG Qi1, YU Wenjie1, HE Wei1, LIAO Sha1, SHEN Yunzhou3, HUANG Yan1, WANG Qian1, ZHONG Bo1, LIU Yang1     
1. Sichuan Center for Disease Control and Prevention, Chengdu 610041, China;
2. Ganzi Tibetan Autonomous Prefecture Center for Disease Control and Prevention, Kangding 626099, China;
3. Liangshan Yi Autonomous Prefecture Center for Disease Control and Prevention, Xichang 615000, China
Abstract: This study was conducted to explore the methods and its feasibility of identifying Canidae species and to investigate the infection of Echinococcus in wild Canidae based on non-invasive fecal sample analysis. A total of 101 DNA samples extracted from feces of wild Canidae were amplified and sequenced with a pair of universal primers designed for the mitochondrial DNA control region of mammalian species. The species origins of the fecal samples were molecularly identified by the analysis of nucleotide variation, the comparison of sequence similarity, the analysis of genetic distance and the construction of phylogenetic tree. Meanwhile, copro-antigen ELISA was used to investigate the infection of Echinococcus. The test success rate for the species identification of the samples was 73.27%. A total of 20 haplotypes were inferred from the 74 sequences obtained in the study. Each of the haplotype sequences was found to match with one single species in the GenBank database with 98.52%-100% sequence similarity. All the haplotypes can be divided into 4 haplotype sets according to sequence differences. The average number of nucleotide differences among the 4 haplotype sets is 17.83-70.10, and the number of nucleotide differences among the haplotypes within each set is 1-10. The genetic distance among the 4 haplotype sets is 0.068-0.342. The value is much larger than that of the genetic distance among the haplotypes within each set, which is 0.009-0.022. The phylogenetic analysis showed that haplotypes in the same sets clustered into same clades with high bootstrap values. Based on all the analysis results, it can be inferred that the sources of the fecal samples were Canis lupus familiaris, Canis lupus, Vulpes vulpes and Vulpes ferrilata respectively. A total of 11 positive samples were found in the copro-antigen ELISA test. The positive rate of Echinococcus antigen was 10.89%. The analysis of mitochondrial DNA control region sequence can be used for the molecular identification of the Canidae. The detection accuracy rate would be further improved if other molecular markers such as nuclear genes were aslo used. The positive rate of Echinococcus antigen is high in feces of wild Canidae in the study area.
Key words: Canidae    Echinococcus    definitive host    mtDNA    species identification    

犬科动物是棘球绦虫的主要终末宿主和传染源。了解犬科动物感染虫种及其感染水平,可以明确不同终末宿主在棘球蚴病传播中所起作用和风险程度,为棘球绦虫终末宿主的防控提供重要的策略依据。由于样本获取困难、检测方法有限,长期以来,包括流浪犬在内的棘球绦虫野外终宿主感染水平的调查始终是一个棘手的问题。出于对野生动物保护和食肉动物危险性的考虑,早期采用的捕杀剖检或药物导泄收集虫体的方法已经不再适用。非损伤性采样的粪样检测是目前最具有操作性的方法。与家犬不同,野外动物粪样为野外捡拾,来源并不明确。不同动物粪便虽然具有一定的外观形态特征,但单凭眼观鉴别其来源,仍可能出现偏差。简单进行病原检测而不对样本来源进行准确辨识,很可能会因为纳入其他非目标宿主样本而造成检测阳性率远低于实际感染率,或掩盖单一宿主的高感染率,难以保证结果的准确性。因此,利用分子生物学方法进行样本的溯源极为重要。

目前,国内棘球绦虫终末宿主种类分子鉴别的研究不多,且主要集中于单一属内物种的鉴别[1-2]。线粒体DNA(mtDNA)因具有多拷贝、变异快、缺乏重组、分子量小等特点,被广泛应用于各种动、植物的分子鉴定和分类,已证实对多种哺乳动物具有良好的区分效果[3-6]。本研究以采集自棘球绦虫高流行区野外环境中的犬科动物粪便为样本,以mtDNA控制区(D-loop)基因为分子标记,对不同犬科动物同源序列进行比对分析,探索mtDNA D-loop高变区序列用于棘球绦虫主要终末宿主鉴别的可能性,同时,通过粪-抗原检测,对样本棘球绦虫感染情况进行调查。

1 材料与方法 1.1 样本来源

本研究检测粪样101份,于2016年采集自甘孜州石渠县真达、呷依、格孟3乡的野外环境,包括居民定居点周围和狐狸、狼等棘球绦虫野外终末宿主洞穴口及经常出没的区域。根据粪样外观形态和采集地点,对粪样来源宿主种类进行初步鉴定,将采集的样本分为3组:1)疑似狼粪样32份,主要采自狼洞口和高山、半高山牛、羊、藏原羚等大型哺乳动物活动较多的地方,粪便较粗长,颜色泛白,带有少量动物长毛发和较多、较大的动物骨头;2)疑似狐粪样30份,主要采自狐洞口和半高山、沿河高原草甸高原鼠兔、青海田鼠等小型哺乳动物活动频繁的地区,粪便较细长,呈灰黑色,带有较多动物短毛发;3)疑似犬粪样本39份,主要采自定居点周边高原鼠兔、青海田鼠等小型哺乳动物活动较多的地区,粪便粗细、长短多介于狼和狐之间,呈棕黑或灰黑色,带有动物短毛发。

所有样本于-80 ℃冷冻保存1个月以上,以灭活虫卵。

1.2 主要试剂和仪器

主要使用试剂包括QIAamp DNA Stool Mini Kit(Qiagen, Germany),GoTaq G2 Hot Start Master Mixes(Promega, USA),犬粪包虫抗原检测试剂盒(深圳市康百得生物科技有限公司,中国)。主要使用仪器包括AB2720 PCR仪(Applied Biosystems,USA),电泳仪(BioRad,USA),凝胶成像仪(BioRad,USA),电热恒温培养箱(上海跃进医疗器械厂,中国)。

1.3 基因组DNA提取

参考文献[1]方法对粪样中可能存在的宿主肠道脱落组织、细胞碎片等进行富集处理后,按照粪便DNA提取试剂盒操作说明提取粪样基因组DNA,于-20 ℃保存备用。

1.4 宿主物种鉴别

1.4.1 PCR扩增和测序   参考文献[7]合成用以扩增哺乳动物mtDNA D-loop区的通用引物L15995:5′ -CTCCACTATCAGCACCCAAAG-3′,H16498:5′-CCTGAAGTAAGAACCAGATG-3′。PCR反应体积为25 μL,包含1×GoTaq G2 Hot Start Master Mixes,上下游引物0.4 μmol·L-1,DNA模板20~50 ng。扩增条件为94 ℃预变性2 min;94 ℃变性30 s,55 ℃退火30 s,72 ℃延伸90 s,共35个循环;72 ℃延伸10 min。PCR扩增产物经2%琼脂糖凝胶电泳检测后,送擎科生物技术公司进行双向测序。

1.4.2 序列分析   所获序列经DNAStar软件进行自动双向拼接并加以适当手动调整后,采用ClustalX2软件[8]进行多重比对。使用DnaSP5软件[9]分析序列单倍型与变异位点。将单倍型序列与Genbank数据库中已发表物种序列进行BLAST比对。利用MEGA7软件[10]计算基于Kimura-2-parameter模型(K2P)的遗传距离,同时以小家鼠(Mus musculus)为外群,采用最大似然法(Maximum Likehood,ML)构建系统发育树。

1.5 棘球绦虫抗原检测

1.5.1 粪-抗原检测   粪样经适当捣碎处理后,按照犬粪包虫抗原检测试剂盒操作说明,采用ELISA夹心法对粪样中棘球绦虫抗原进行检测。

1.5.2 数据处理   采用EXCEL软件对数据进行整理,采用SPSS软件对数据进行统计分析,使用χ2检验对不同来源样本棘球绦虫抗原阳性率差异进行比较。

2 结果 2.1 宿主物种鉴别

2.1.1 PCR扩增结果   提取的101份粪样DNA中有74份扩增出目的条带,PCR产物电泳结果显示在350~400 bp出现清晰的单一条带,扩增产物可用于后续测序分析,样本检测成功率为73.27%(图 1)。

1~9. 检测的粪便样本;M. DL5000 DNA相对分子质量标准;N. 无模板对照 1-9. Fecal samples detected in the study; M. DL5000 DNA molecular weight markers; N. No template control 图 1 PCR扩增结果 Fig. 1 The results of the PCR amplification

2.1.2 序列变异与单倍型分析   所获74条序列经比对、剪切后,最终整理为297 bp的同源序列(GenBank登录号:MZ189011~MZ189084),定义单倍型20个,根据序列碱基差异大小,分为4个单倍型集合,其中单倍型H1~H6核苷酸变异位点13个,6个单倍型间碱基差异为1~10个;单倍型H7~H10核苷酸变异位点5个,4个单倍型间碱基差异为1~5个;单倍型H11~H15核苷酸变异位点11个,5个单倍型间碱基差异为1~8个;单倍型H16~H20核苷酸变异位点5个,5个单倍型间碱基差异为2~5个(表 12)。单倍型集间,以单倍型H1~H6与H7~H10碱基差异最小,平均为17.83,其次为单倍型H11~H15与H16~H20,平均碱基差异为38.64,单倍型H1~H6、H7~H10与单倍型H11~H15、H16~H20间碱基差异较大,平均为64.85~70.10。

表 1 单倍型H1~H10线粒体控制区基因变异位点 Table 1 Mutation sites in the mtDNA control region of H1-H10
表 2 H11~H20线粒体控制区基因变异位点 Table 2 Mutation sites in the mtDNA control region of H11-H20

2.1.3 序列相似度比对   经Blast比对,所有单倍型序列均与GenBank数据库中已公布的单一物种同源序列高度相似,对比结果分别为犬(Canis lupus familiaris)、灰狼(Canis lupus)、红狐(Vulpes vulpes)和藏狐(Vulpesferrilata)4种犬科动物(表 3)。

表 3 单倍型序列 BLAST  比对结果 Table 3 BLAST results of the haplotype sequences

2.1.4 遗传距离分析   基于K2P模型的遗传距离分析结果显示,20个单倍型之间遗传距离为0.004~0.354。4个单倍型集间遗传距离远大于各集合内单倍型间遗传距离;单倍型H1~H6与H7~ H10、单倍型H11~H15与H16~H20间遗传距离较与其他单倍型集间遗传距离小;单倍型H1~H6、H7~H10、H11~H15、H16~H20分别与犬、灰狼、赤狐和藏狐遗传关系最近(表 4)。

表 4 基于  K2P  模型的平均遗传距离 Table 4 Average genetic distance based on the K2P model

2.1.5 系统发育分析   构建的ML系统发育树表现出明显的聚类现象,所有单倍型分为两大支,其中单倍型H7~H10与灰狼参考序列(KT901460)聚集在一起,并与单倍型H1~H6及犬参考序列(JN695048)共处同一大分支;另一大分支是由单倍型H11~H15与赤狐参考序列(KJ846513),以及单倍型H16~H20与藏狐参考序列(JF520840)各自聚集成簇后再形成的一个大支(图 2)。

小于75%的自展值未予显示 Bootstrap values below 75% were not shown 图 2 基于线粒体控制区基因序列构建的ML系统进化树 Fig. 2 ML tree based on mtDNA control region sequences
2.2 棘球绦虫感染检测

经棘球绦虫粪-抗原ELISA检测共发现阳性样本11份,阳性率为10.89%(11/101)(表 5)。

表 5 棘球绦虫粪-抗原  ELISA  检测结果 Table 5 The results of copro-antigen ELISA tests for detection of Echinococcus
3 讨论

棘球绦虫终末宿主感染调查是棘球蚴病监测的重点,也是防控的难点和薄弱点。粪便作为非损伤性采集样本,可同时用于来源宿主的物种鉴定和棘球绦虫的感染检测,对棘球绦虫终末宿主,尤其是野生终末宿主的调查研究具有重要意义。近年来的研究报道中所涉及的野外犬科动物感染调查均选择以粪便为样本,在检测其内棘球绦虫特异性抗原或核酸的同时,根据粪样外形特征来推测其宿主来源[11-12]。本次研究结果证实,以粪样形态特点和采样环境来判断样本来源动物种类具有一定可行性,但也存在一定的误差,尤其是犬粪,误辨率较高。另外,由于眼观鉴别严重依赖现场人员个人经验,也极大地限制了该方法的推广应用。本次检测的101份粪样中,74份通过粪样DNA扩增、测序和分析,成功实现对其来源的4种犬科动物的鉴别,结果表明,分子检测可有效鉴定粪便样本来源动物种类,作为形态判别的修正依据,有效提升鉴定准确率。另一方面,研究结果也反映出粪样检测的不足,与血液、组织等样本相比,粪便样本检测成功率一般相对较低。本次用于检测的粪样为野外捡拾,部分样本可能因为曝露时间过长,受到高原阳光爆嗮或雨水冲刷,造成粪样被某些有抑制作用的成分污染,或其内含有的宿主核酸降解,而未能成功扩增。研究结果提示,要进一步提升样本检出率,应尽量采集新鲜粪便用于分子检测,同时,探索更为有效的目的成分富集前处理方法。

mtDNA目前被广泛应用于物种鉴定、遗传多样性和系统发生等研究,以cox1、细胞色素b(cytb)等编码基因为代表的DNA条形码技术在啮齿类、鱼类、鸟类等动物识别与分类研究中均取得了较为理想的结果[13-15]。然而,由于受到选择压力的影响,这些编码基因通常无法有效区分亲缘关系较近的物种,相较之下,作为mtDNA中唯一的非编码区域,D-loop区因为缺乏选择压力,变异速度较编码区基因更快,被认为对近缘和近期分化的物种具有更好的分辨能力[7]。本研究对成功扩增获得的74份mtDNA D-loop区基因片段进行了分析,同源序列比对结果显示,检测序列可与GenBank中单一物种序列相匹配,多态性位点差异与遗传距离分析结果与BLAST比对结果一致,所鉴定的物种间距离远大于种内距离,满足物种鉴定的判别条件,系统发育发育树中,鉴别为同一物种的不同单倍型以高自展值进行聚类,支持其鉴别结果,综合各项分析结果,可判断粪样来源分别为犬、灰狼、赤狐和藏狐4种犬科动物。研究表明,以mtDNA D-loop区序列分析作为犬科动物粪样不同物种来源的鉴定依据具有可行性,该方法对包括狼及其近缘驯化种犬在内的棘球绦虫主要终末宿主均具有较好的辨识能力。

尽管基于mtDNAD-loop区对犬科动物的鉴别具有快捷、高效、简易等优点,但仍然存在一定的局限性。由于线粒体基因母系遗传的特点,基于mtDNA的分析无法识别自然界中可能偶尔发生的杂交事件,对没有生殖隔离的物种间自然交配产生的杂交及其回交子代个体的鉴别存在一定缺陷[16-17]。因此,为进一步提高鉴别的准确度,在对mtDNA D-loop区碱基序列进行分析比较的同时,可结合常染色体基因的双亲遗传或Y染色体基因的父系遗传分子标记对物种进行鉴别[18-20]。此外,采用微卫星标记对物种个体进行识别分析,可进一步排除可能存在的来源于同一宿主的重复样本,实现对棘球绦虫终末宿主的精准识别,为棘球绦虫终末宿主防控提供更为系统、全面的数据支撑。

家犬作为牧区重要的生产工具和人类的亲密伙伴,与其相关的因素,如家庭是否养犬、犬只数量、是否与犬亲密接触等,被认为是我国棘球蚴病流行区、尤其是青藏高原地区棘球绦虫传播最重要的风险因素[21-23],但自2006年我国棘球蚴病防治项目实施以来,经过近15年以家犬驱虫为主的综合防控,目前各流行区家犬棘球绦虫感染率均已出现明显下降,并维持于较低水平[24-26]。与此同时,流浪犬(无主犬)因其数量多、流动性大等特点,干预十分困难,其感染控制成效极为有限[27-28],类似地,作为棘球绦虫的保虫宿主,狼、狐狸等野生犬科动物感染率也远高于家犬[1, 11]。本次对粪样棘球绦虫抗原检测结果显示,调查地区野外犬科动物具有较高的棘球绦虫抗原阳性率,但相较于此前采用不同检测方法对该地区野生犬科动物棘球绦虫感染情况的调查结果,感染率有所下降,这可能是由于各研究项目调查时间、调查方法、样本代表性,以及所涉及宿主种类等的不同而导致[29-31]。由于现有针对野外犬科动物棘球绦虫感染情况的研究严重不足,依靠现有数据难以对当地野外犬科动物棘球绦虫感染水平作出全面的有效推断,但可以明确的是,在家犬棘球绦虫感染率得到有效控制的情况下,包括流浪犬在内的野外犬科动物所具有的高感染率可能对流行区棘球绦虫的传播起到更大的作用,由其引起人群感染的危害不应被忽视。强化宿主物种鉴别与个体鉴定、病原分子检测等关键技术,开展野外犬科动物棘球绦虫感染水平监测,有针对性地提出对野外犬科动物的防控措施,将有助棘球蚴病下一步防控目标的实现。

4 结论

通过mtDNA D-loop区序列分析对野外犬科动物粪便的物种鉴别结果显示,成功确定检测的粪样分别来自犬、灰狼、赤狐和藏狐4种犬科动物,同时,基于棘球绦虫粪抗原的检测发现,样本来源个体具有较高的棘球绦虫感染率。研究结果可为进一步提升棘球绦虫野外终末宿主精准监测质量提供技术依据,作为调查地区野生犬科动物棘球绦虫感染数据的更新和补充,为当地开展棘球绦虫防控提供数据参考。

参考文献
[1]
JIANG W B, LIU N, ZHANG G T, 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. DOI:10.1007/s00436-012-2993-8
[2]
江泓慧, 李波, 马跃, 等. 利用线粒体DNA控制区鉴别狼与狗可行性分析[J]. 野生动物学报, 2018, 39(2): 277-285.
JIANG H H, LI B, MA Y, et al. Feasibility analysis of species identification of dogs and wolves by partial sequence of mtDNA control region[J]. Chinese Journal of Wildlife, 2018, 39(2): 277-285. DOI:10.3969/j.issn.1000-0127.2018.02.011 (in Chinese)
[3]
HSIEH H M, CHIANG H L, TSAI L C, et al. Cytochrome b gene for species identification of the conservation animals[J]. Forensic Sci Int, 2001, 122(1): 7-18. DOI:10.1016/S0379-0738(01)00403-0
[4]
NAKAMURA H, MURO T, IMAMURA S, et al. Forensic species identification based on size variation of mitochondrial DNA hypervariable regions[J]. Int J Legal Med, 2009, 123(2): 177-184. DOI:10.1007/s00414-008-0306-7
[5]
TOBE S S, LINACRE A M T. A multiplex assay to identify 18 European mammal species from mixtures using the mitochondrial cytochrome b gene[J]. Electrophoresis, 2008, 29(2): 340-347. DOI:10.1002/elps.200700706
[6]
MUANGKRAM Y, WAJJWALKU W, AMANO A, et al. The novel primers for mammal species identification-based mitochondrial cytochrome b sequence: implication for reserved wild animals in Thailand and endangered mammal species in Southeast Asia[J]. Mitochondrial DNA Part A, 2018, 29(1): 62-72. DOI:10.1080/24701394.2016.1238902
[7]
PUN K M, ALBRECHT C, CASTELLA V, et al. Species identification in mammals from mixed biological samples based on mitochondrial DNA control region length polymorphism[J]. Electrophoresis, 2009, 30(6): 1008-1014. DOI:10.1002/elps.200800365
[8]
LARKIN M A, BLACKSHIELDS G, BROWN N P, et al. Clustal W and Clustal X version 2.0[J]. Bioinformatics, 2007, 23(21): 2947-2948. DOI:10.1093/bioinformatics/btm404
[9]
LIBRADO P, ROZAS J. DnaSP v5:a software for comprehensive analysis of DNA polymorphism data[J]. Bioinformatics, 2009, 25(11): 1451-1452. DOI:10.1093/bioinformatics/btp187
[10]
KUMAR S, STECHER G, TAMURA K. MEGA7:molecular evolutionary genetics analysis version 7. 0 for bigger datasets[J]. Mol Biol Evol, 2016, 33(7): 1870-1874. DOI:10.1093/molbev/msw054
[11]
李清德, 钟卫东, 赵婧, 等. 民和县南大山林区野生动物棘球绦虫病的调查[J]. 青海畜牧兽医杂志, 2019, 49(2): 56-57.
LI Q D, ZHONG W D, ZHAO J, et al. Investigation on echinococcosis in wild animals in the forest areas of Nanda Mountain of Minhe[J]. Chinese Qinghai Journal of Animal and Veterinary Sciences, 2019, 49(2): 56-57. DOI:10.3969/j.issn.1003-7950.2019.02.012 (in Chinese)
[12]
韩帅, 余晴, 杨诗杰, 等. 棘球蚴病野外犬科类传染源网格化干预与评价[J]. 中国寄生虫学与寄生虫病杂志, 2018, 36(5): 495-498.
HAN S, YU Q, YANG S J, et al. Evaluation on the measures by a grid-based design for wildlife control in hyper-endemic areas for echinococcosis[J]. Chinese Journal of Parasitology and Parasitic Diseases, 2018, 36(5): 495-498. (in Chinese)
[13]
陈宝宝, 孙养信, 范锁平, 等. DNA条形码技术在陕西省鼠形动物鉴定中的应用[J]. 中国人兽共患病学报, 2020, 36(4): 325-329, 339.
CHEN B B, SUN Y X, FAN S P, et al. Application for identification of Murine-like animals by DNA barcoding in Shaanxi Province[J]. Chinese Journal of Zoonoses, 2020, 36(4): 325-329, 339. (in Chinese)
[14]
TAVARES E S, BAKER A J. Single mitochondrial gene barcodes reliably identify sister-species in diverse clades of birds[J]. BMC Evol Biol, 2008, 8: 81. DOI:10.1186/1471-2148-8-81
[15]
陆键萍, 姚琳, 信红梅, 等. 线粒体CO Ⅰ、Cyt b和16S rRNA基因在6种金枪鱼鉴定中的适用性分析[J]. 渔业科学进展, 2020, 41(5): 73-81.
LU J P, YAO L, XIN H M, et al. Applicability analysis of mitochondrial CO Ⅰ, Cyt b and 16S rRNA genes in identification of six tuna species[J]. Progress in Fishery Sciences, 2020, 41(5): 73-81. (in Chinese)
[16]
ADAMS J R, LEONARD J A, WAITS L P. Widespread occurrence of a domestic dog mitochondrial DNA haplotype in southeastern US coyotes[J]. Mol Ecol, 2003, 12(2): 541-546. DOI:10.1046/j.1365-294X.2003.01708.x
[17]
REESE E M, WINTERS M, BOOTH R K, et al. Development of a mitochondrial DNA marker that distinguishes domestic dogs from Washington state gray wolves[J]. Conserv Genet Resour, 2020, 12(3): 497-501. DOI:10.1007/s12686-020-01130-2
[18]
ANDERSONE Ž, LUCCHINI V, OZOLIŠ J, et al. Hybridisation between wolves and dogs in Latvia as documented using mitochondrial and microsatellite DNA markers[J]. Mamm Biol, 2002, 67(2): 79-90. DOI:10.1078/1616-5047-00012
[19]
VILÀ C, WALKER C, SUNDQVIST A K, et al. Combined use of maternal, paternal and bi-parental genetic markers for the identification of wolf-dog hybrids[J]. Heredity, 2003, 90(1): 17-24. DOI:10.1038/sj.hdy.6800175
[20]
MUÑOZ-FUENTES V, DARIMONT C T, PAQUET P C, et al. The genetic legacy of extirpation and re-colonization in Vancouver Island wolves[J]. Conserv Genet, 2010, 11(2): 547-556. DOI:10.1007/s10592-009-9974-1
[21]
王谦, 邱加闽, SCHANTZ P, 等. 四川省西部藏区家庭饲养牲畜人群包虫病风险因素的调查[J]. 中国寄生虫学与寄生虫病杂志, 2001, 19(2): 93-96.
WANG Q, QIU J M, SCHANTZ P, et al. Investigation of risk factors for development of human hydatidosis among households raising livestock in Tibetan areas of western Sichuan province[J]. Chinese Journal of Parasitology and Parasitic Diseases, 2001, 19(2): 93-96. DOI:10.3969/j.issn.1000-7423.2001.02.009 (in Chinese)
[22]
刘平, 孙向东, 康京丽, 等. 包虫病相关影响因素Meta分析[J]. 中国动物检疫, 2018, 35(10): 13-18.
LIU P, SUN X D, KANG J L, et al. Meta-analysis on potential risk factors associated with human echinococcosis[J]. China Animal Health Inspection, 2018, 35(10): 13-18. DOI:10.3969/j.issn.1005-944X.2018.10.004 (in Chinese)
[23]
CRAIG P S, GIRAUDOUX P, SHI D, et al. An epidemiological and ecological study of human alveolar echinococcosis transmission in south Gansu, China[J]. Acta Trop, 2000, 77(2): 167-177. DOI:10.1016/S0001-706X(00)00134-0
[24]
何伟, 尚婧晔, 黄燕, 等. 2007-2013年四川省包虫病疫情监测结果分析[J]. 预防医学情报杂志, 2016, 32(1): 83-85.
HE W, SHANG J Y, HUANG Y, et al. Surveillance of hydatidosis in Sichuan, 2007-2013[J]. Journal of Preventive Medicine Information, 2016, 32(1): 83-85. (in Chinese)
[25]
王玉莹, 郝慧霞, 宋健, 等. 2012—2018年内蒙古自治区犬棘球绦虫感染空间分布特征[J]. 中国血吸虫病防治杂志, 2021, 33(1): 40-47, 53.
WANG Y Y, HE H X, SONG J, et al. Spatial distribution of canine echinococcosis in Inner Mongolia Autonomous Region from 2012 to 2018[J]. Chinese Journal of Schistosomiasis Control, 2021, 33(1): 40-47, 53. (in Chinese)
[26]
马玉安, 尚文杰. 甘南藏族自治州棘球蚴病流行情况分析[J]. 中国寄生虫学与寄生虫病杂志, 2015, 33(5): 352-356.
MA Y A, SHANG W J. Endemic situation of echinococcosis in Gannan Tibetan Autonomous Prefecture[J]. Chinese Journal of Parasitology and Parasitic Diseases, 2015, 33(5): 352-356. (in Chinese)
[27]
蔡辉霞, 王虎, 韩秀敏, 等. 青海高原棘球绦虫不同宿主感染情况及意义的研究[J]. 中国地方病学杂志, 2012, 31(3): 296-300.
CAI H X, WANG H, HAN X M, et al. Current situation of Echinococcus infection in different hosts in Qinghai plateau[J]. Chinese Journal of Endemiology, 2012, 31(3): 296-300. DOI:10.3760/cma.j.issn.1000-4955.2012.03.016 (in Chinese)
[28]
蔡辉霞, 王虎, 韩秀敏, 等. 1990-2010年青海高原终宿主动物棘球绦虫感染与儿童棘球蚴病的相关研究[J]. 中国人兽共患病学报, 2012, 28(5): 500-502.
CAI H X, WANG H, HAN X M, et al. Correlation between definitive hosts of Echinococcus and echinococcosis in children in Qinghai plateau, China, 1990-2010[J]. Chinese Journal of Zoonoses, 2012, 28(5): 500-502. DOI:10.3969/j.issn.1002-2694.2012.05.023 (in Chinese)
[29]
邱加闽, 陈兴旺, 任敏, 等. 青藏高原泡球蚴病流行病学研究[J]. 实用寄生虫病杂志, 1995, 3(3): 106-109.
QIU J M, CHEN X W, REN M, et al. Epidemiological study on alveolar hydatid disease in Qinghai-Xizang plateau[J]. Journal of Practical Parasitic Diseases, 1995, 3(3): 106-109. (in Chinese)
[30]
何金戈, 邱加闽, 刘凤洁, 等. 四川西部藏区包虫病流行病学研究Ⅱ. 牲畜及野生动物两型包虫病感染状况调查[J]. 中国人兽共患病杂志, 2000, 16(5): 62-65.
HE J G, QIU J M, LIU F J, et al. Epidemiological survey on hydatidosis in Tibetan region of Western Sichuan Ⅱ. Infection situation among domestic and wild animals[J]. Chinese Journal of Zoonoses, 2000, 16(5): 62-65. DOI:10.3969/j.issn.1002-2694.2000.05.019 (in Chinese)
[31]
陈珏. 基于粪便DNA调查藏狐的棘球绦虫感染率[D]. 上海: 华东师范大学, 2014.
CHEN J. Specific detection of the prevalence of Echinococcus in Tibetan foxes with fecal DNA[D]. Shanghai: East China Normal University, 2014. (in Chinese)

(编辑   白永平)