云南省埃及伊蚊对拟除虫菊酯类抗性群体的击倒抗性基因突变分析
兰学梅, 杨明东, 杨锐, 姜进勇
云南省寄生虫病防治所,云南省虫媒病毒研究中心,云南省虫媒传染病防控研究重点实验室,金宁一院士工作站,面向南亚东南亚热带病国际科技人员交流与教育培训基地,云南省寄生虫病防治所虫媒传染病防控关键技术省创新团队(培育),普洱 665099
通讯作者:姜进勇,Email: yipdjiang@126.com; ORCID: 0000-0001-7891-1259
摘要

目的 检测云南省登革热重点地区对氯菊酯和高效氯氟氰菊酯杀虫剂已产生抗性的埃及伊蚊的击倒抗性(knockdown resistance, kdr)基因,进行基因突变分析,阐明抗性表型与 kdr基因突变的关系。方法 收集云南省景洪市、勐腊县、勐海县、瑞丽市和耿马县经过氯菊酯和高效氯氟氰菊酯杀虫剂测定的埃及伊蚊成蚊样本,分别提取单蚊基因组DNA,Allele-specific PCR(AS-PCR)扩增和分析其 kdr部分基因片段,统计抗性与敏感表型中 kdr突变的基因型和基因频率。结果 共580份雌性埃及伊蚊样本被用于V1016G和F1534C击倒抗性突变的检测。V1016G的突变频率为99.83%(579/580),F1534C的突变频率为46.38%(269/580)。V1016G和F1534C突变频率差异存在统计学意义(χ2=421.338, P=0.000)。5个埃及伊蚊种群中,V1016G突变率除了瑞丽市外,其他区域均达到100.00%。F1534C突变率景洪市最高为84.16%,耿马县最低为5.66%。抗性表型与敏感表型的V1016G和F1534C突变率差异均存在统计学意义(χ2=7.298, P=0.007;χ2=14.010, P=0.000)。580份样本中,有65份样本同时存在V1016G和F1534C位点突变,突变率为11.21%;不同地区同时突变率范围为0~33.66%,景洪市同时突变率最高为33.66%。结论 云南省登革热重点地区埃及伊蚊种群存在V1016G和F1534C突变,V1016G突变率高且分布广,F1534C突变以景洪市埃及伊蚊种群为最高。敏感表型与抗性表型V1016G和F1534C突变存在差异,提示抗性产生与突变存在一定的关系。

关键词: 登革热; 埃及伊蚊; 杀虫剂抗性; 击倒抗性; 基因突变
中图分类号:R378.5 文献标志码:A 文章编号:1002-2694(2020)12-0993-07
Analysis of knockdown resistance gene mutations in pyrethroid insecticide resistant populations of Aedes aegypti in Yunnan Province
LAN Xue-mei, YANG Ming-dong, YANG Rui, JIANG Jin-yong
Yunnan Institute of Parasitic Diseases,Yunnan Provincial Center of Arborvirus Research, Yunnan Provincial Key Laboratory of Vector-borne Diseases Control and Research, Academician Workstation of Professor Jin Ningyi, Training Base of International Scientific Exchange and Education in Tropical Diseases for South and Southeast Asia, Yunnan Institute of Parasitic Diseases Innovative Team of Key Techniques for Vector Borne Disease Control and Prevention (Developing), Pu’er 665099, China
Corresponding author: Jiang Jin-yong, Email:yipdjiang@126.com
Abstract

To detect the knockdown resistance ( kdr) gene mutations against permethrin and cypermethrin insecticides in Ae.aegypti mosquitoes in key areas of dengue fever in Yunnan Province, with the aim of the elucidating the relationship between the resistance phenotype and kdr gene mutations, we collected Ae.aegypti from Jinghong City, Mengla County, Menghai County, Ruili City and Gengma County, and bioassayed them with permethrin and cypermethrin insecticides after species identification. Individual genomic DNA was extracted. Allele-specific PCR (AS-PCR) amplification and analysis of kdr gene fragments, genotypes and gene frequencies of kdr mutations in resistant and sensitive phenotypes were performed. AS-PCR successfully detected V1016G and F1534C mutations in the 580 female wild Ae.aegypti by direct sequencing of the sodium channel gene. The mutation frequency of V1016G was 99.83% (579/580), and that of F1534C was 46.38% (269/580). The mutation frequencies of V1016G and F1534C were statistically significant (χ2=421.338, P=0.000, P<0.05). Among the five Ae.aegypti populations, the V1016G mutation rate reached 100.00% in all regions except Ruili City; the mutation rate of F1534C in Jinghong was the highest at 84.16%, and that in Gengma was the lowest at 5.66%. The mutation rates of V1016G and F1534C in both resistant and sensitive phenotypes were statistically significant (χ2=7.298, P=0.007, P<0.05; χ2=14.010, P=0.000, P<0.05). Among the 580 samples, 65 samples had both V1016G and F1534C site mutations, and the mutation rate was11.21%. The simultaneous mutation rates in different regions ranged from 0 to 33.66%, The maximum mutation rate in Jinghong was 33.66%. V1016G and F1534C mutations were found in the Ae.aegypti populations in key areas of dengue fever in Yunnan Province. The V1016G mutation rate was high and this mutation was widely distributed. The population of Ae.aegypti in Jinghong had the highest rate of F1534C mutation. The sensitive phenotype and resistant phenotype V1016G and F1534C mutation rates differed, thus suggesting a relationship between resistance and mutation.

Key words: dengue fever; Aedes Aegypti; insecticide resistance; knockdown resistance; gene mutations

埃及伊蚊(Aedes Aegypti)是新入侵云南的蚊虫, 2002年首先在瑞丽市姐告口岸被发现[1], 其入侵、定殖和扩散能力强, 目前已成为云南省景洪市、勐腊县、瑞丽市和耿马县孟定镇的优势蚊种[2], 云南省2013年来多次登革热疫情暴发地均有埃及伊蚊分布, 是云南省登革热的重要传播媒介[3, 4, 5]。媒介控制是登革热疫情处置的主要措施之一, 瑞丽、景洪、勐腊和耿马等地主要使用各种卫生化学杀虫剂进行空间喷洒(超低容量喷雾和热烟雾)结合滞留喷洒, 来迅速杀灭成蚊从而控制登革热疫情。各地疫情处置过程中, 长期、大量使用拟除虫菊酯类杀虫剂进行灭蚊控制, 埃及伊蚊可能对菊酯类杀虫剂产生了不同程度的抗药性[6]。因此, 2016年以来, 本课题组先后测定了云南省瑞丽市、景洪市、勐腊县、勐海县和耿马县等登革热重点地区的埃及伊蚊成蚊对几种拟除虫菊酯类杀虫剂的敏感性, 结果证实均存在不同程度抗药性[7]。为进一步研究其产生抗性的机制, 本研究在以往研究的基础上, 收集上述5个种群的埃及伊蚊野生种群成蚊, 在生物测定拟除虫菊酯类杀虫剂(氯菊酯和高效氯氟氰菊酯)敏感性的基础上, 检测其对应kdr基因突变情况, 用于丰富我省埃及伊蚊群体kdr基因突变数据, 为探讨分子检测结果与生物测定结果之间的关联性提供科学依据。

1 材料与方法
1.1 供试蚊虫来源

2016-2018年, 在云南省景洪市、勐腊县、勐海县、瑞丽市和耿马县(东、南、西、北和中)5个片区, 采集埃及伊蚊的幼虫和蛹(均采自疫点及其周围室内外花瓶、水桶、泡菜坛、饮水机水槽、轮胎及其他废弃容器等)。羽化和鉴定[8]后, 饲喂鼠血(豚鼠, 产地:普洱), 分类饲养繁殖至子一代(F1), 进行对杀虫剂的生物学测定。随即进行埃及伊蚊成蚊对3.20%氯菊酯和0.08%高效氯氟氰菊酯的抗药性测定(药膜由中国疾病预防控制中心传染病预防控制所提供), 实验结果均为抗性。24 h后死亡和存活样本单管单只分别收集, 并置于75%乙醇中保存。

1.2 kdr基因扩增及电泳

收集的埃及伊蚊按照试剂盒(购自天根生化科技(北京)有限公司)说明书, 提取单蚊基因组DNA。参照文献[9]合成引物, 扩增1 016和1 534两个位点的突变情况。各引物的PCR反应体系相同:包含2× Taq预混液12.5 μ L, 10 μ mol/L, 非特异性正反引物各1.0 μ L, 特异性内引物1.0 μ L, DNA模板2.0 μ L, 加双蒸灭菌水至 25.0 μ L。反应条件为:94 ℃ 3 min; 94 ℃ 30 s, 60 ℃ 30 s, 72℃ 1 min, 35个循环; 72 ℃延伸7min。扩增产物经1.5%琼脂糖凝胶电泳检测后, 观察电泳条带大小。V1016G等位基因检测中, M2-For和M2-Rev为非特异性外引物, 用于大片段扩增(592 bp); 特异性内引物M2-V和M2-G扩增片段长度为348 bp(仅有V, 敏感纯合子; 仅有G, 抗性纯合子; VG并存, 抗性杂合子)。F1534C等位基因检测中, M3-For和M3-Rev为非特异性外引物, 用于大片段扩增(517 bp); 特异性内引物M3-F和M3-C扩增片段长度为284 bp(仅有F, 敏感纯合子; 仅有C, 抗性纯合子; FC并存, 抗性杂合子)。见图1。

1.3 数据分析

运用SPSS20.0统计软件进行统计分析, 采用χ 2检验分析kdr各基因突变型间的相关性, 检验水准为α =0.05。

图1 埃及伊蚊钠离子通道击倒抗性基因AS-PCR检测示意图Fig.1 Schematic diagram of the AS-PCR test used to detect the V1016G (A) and F1534C (B) mutations in the sodium channel gene of Ae.aegypti

2 结 果

采集自景洪市、勐腊县、勐海县、瑞丽市和耿马县5个地区, 共计580份雌性埃及伊蚊样本被用于V1016G和F1534C击倒抗性突变的检测。5个重点地区埃及伊蚊均发现有V1016G和F1534C突变存在, 其中, V1016G的突变率为99.83%(579/580), 5个地区V1016G的突变率都较高, 除了瑞丽市为99.26%外, 其余突变率均为100%。F1534C的突变率为46.38%(269/580), 5个地区F1534C的突变率分别为84.16%(85/101)、67.80%(80/118)、49.17%(59/120)、28.89%(39/135)和5.66%(6/106)。V1016G和F1534C突变率差异存在统计学意义(χ 2=421.338, P=0.000), V1016G的突变率高于F1534C(表1)。

表1 云南省埃及伊蚊对菊酯类杀虫剂基因突变率(%)测定 Tab.1 Gene mutation frequencies under pyrethroid insecticide treatment in Ae.aegypti in Yunnan Province, China
2.1 不同杀虫剂V1016G和F1534C突变的比较

埃及伊蚊对氯菊酯抗性的样本检测了273份, V1016G和F1534C突变率分别为100.00%(273/273)和50.92%(139/273)。对氯菊酯的V1016G和F1534C突变率进行比较, 差异存在统计学意义(χ 2=177.583, P=0.000); 埃及伊蚊对高效氯氟氰菊酯抗性的样本检测了307份, V1016G和F1534C突变率分别为99.67%(306/307)和47.62%(130/273)。对高效氯氟氰菊酯V1016G和F1534C突变率差异存在统计学意义(χ 2=209.802, P=0.000)(见表2、3)。

对氯菊酯和高效氯氟氰菊酯的V1016G突变率进行比较, 差异无统计学意义(χ 2=0.891, P=0.345); 对氯菊酯和高效氯氟氰菊酯的F1534C突变率进行比较, 差异有统计学意义(χ 2=4.268, P=0.039)(表2、3)。

2.2 不同地区V1016G和F1534C突变的比较

采用χ 2检验, 5个地区1 016位点的基因突变率进行比较, 差异无统计学意义(χ 2=3.302, P=0.509); 5个地区1 534位点的基因突变率进行比较, 差异有统计学意义(χ 2=167.383, P=0.000)。进一步进行景洪市与其他4个地区1 534位点的基因突变率进行两两比较, 其中景洪和勐腊、勐海、瑞丽、耿马比较的统计结果分别为:χ 2=7.842, P=0.005; χ 2=29.578, P=0.000; χ 2=70.778, P=0.000; χ 2=129.365, P=0.000。

表2 云南省埃及伊蚊击倒抗性基因突变频率测定 Tab.2 Frequencies of kdr mutations in Ae.aegypti in Yunnan Province, China
表3 不同杀虫剂V1016G和F1534C突变率比较 Tab.3 Comparison the mutation frequencies of V1016G and F1534C under treatment with different insecticides
2.3 抗性表型与敏感表型V1016G和F1534C突变的比较

5个地区埃及伊蚊野生群体对氯菊酯和高效氯氟氰菊酯进行抗药性测定后, 实验存活个体定义为抗性表型, 实验死亡个体定义为敏感表型。抗性表型样本检测510份, V1016G和F1534C突变频率分别为100.00%(510/510)和49.02%(250/510); 敏感表型样本检测70份, V1016G和F1534C突变频率分别为98.57%(69/70)和72.86%(51/70)。抗性表型与敏感表型的V1016G突变率进行比较, 差异有统计学意义(χ 2=7.298, P=0.007); 抗性表型与敏感表型的F1534C突变率进行比较, 差异有统计学意义(χ 2=14.010, P=0.000)。

2.4 V1016G和F1534C位点同时纯合子突变

580份样本中, 有65份样本同时存在V1016G和F1534C位点纯合子突变, 突变率为11.21%; 不同地区同时突变率范围为0%~33.66%, 景洪同时突变率最高为33.66%, 瑞丽最低, 本次实验耿马未发现同时存在纯合子突变(表1)。

3 讨 论

击倒抗性(knockdown resistance, kdr)是引起蚊虫对拟除虫菊酯类杀虫剂产生抗性的机制之一, 埃及伊蚊电压门控钠离子通道(voltage-gated sodium channel, VGSC)上很多位点已被研究检测出来, 但仅有少数位点突变被证实与抗性相关[10, 11], 如埃及伊蚊不同抗性群体, VGSC基因突变的密码子有3个 (I1011M/V、V1016G/I和F1534C)与DDT和菊酯类杀虫剂抗性相关。

本次研究景洪、勐腊、勐海、瑞丽和耿马的埃及伊蚊存在V1016G和F1534C突变, V1016G的突变以抗性纯合子(RR)突变为主, 与师灿南[12]2015年对景洪市城区埃及伊蚊的检测结果一致, 缅甸[13]和泰国[14]也发现相似的突变。然而, Sebastien Marcombe等2014年在老挝埃及伊蚊中发现存在V1016G和F1534C突变, 以F1534C突变占主导地位, 突变率高达96%, 而V1016G存在低水平的突变, 突变率仅为36%[15], Hitoshi Kawada等报告越南V1016G突变率更低, 仅为0.26%[16]。存在不同的原因可能是入侵云南埃及伊蚊与老挝和越南不是同一来源, 或是各国使用的杀虫剂种类、时间和数量不一致, 尚需进一步分析和研究。通过两种杀虫剂V1016G突变率和5个不同抗性种群V1016G突变率的比较看出, 差异均无统计学意义, 进一步证实V1016G突变广泛存在于对拟除虫菊酯类杀虫剂产生抗药性的埃及伊蚊中。Han-Hsuan Chung等研究发现台湾地区的埃及伊蚊S989P和D1763Y突变均与V1016G并发。S989P和V1016G并存[17, 18], 以及D1763Y和V1016G同时并存[19]。至于云南的埃及伊蚊中是否存在上述现象, 尚需进一步研究。

5个地区埃及伊蚊的1 534位点均发生了基因突变, 且不同地区埃及伊蚊的突变率存在明显的差异, 以景洪最高为84.16%, 耿马最低为5.66%。原因可能是每个地区杀虫剂的种类、数量和使用年限不同所致, 也可能与不同地区埃及伊蚊对拟除虫菊酯类杀虫剂的抗性程度不同有关, 需检测不同地区埃及伊蚊对杀虫剂的抗性程度类进行证实。本研究中, 抗性表型与敏感表型V1016G和F1534C突变率均存在差异, 两种位点的突变均以抗性表型的突变率较高, 与其他研究者的结论一致[20, 21]。证实V1016G和F1534C突变与埃及伊蚊对拟除虫菊酯类杀虫剂产生抗性具有相关性, 是引起杀虫剂抗性的机制之一。

多个位点突变同时存在的情况下会导致埃及伊蚊对拟除虫菊酯类杀虫剂产生更高水平的抗性[13], 泰国埃及伊蚊中G1016和C1534 同时突变对溴氰菊酯抗性有相加作用。台湾地区的埃及伊蚊中发现同时存在S989P / V1016G / F1534C或V1016G / F1534C / D1763Y的3个杂合突变的显现[22], 结果对杀虫剂具有较高的抗性。云南埃及伊蚊1 016和1 534位点突变同时为抗性纯合子(RR)的频率为11.21%(65/580)远远高于老挝埃及伊蚊(4/1 076, 0.37%)[17]。景洪市同时突变为抗性杂合子的频率最高为33.66%。可能与不同地区埃及伊蚊对杀虫剂产生的抗性程度不同有关。至于不同位点三重杂合突变的现象是否具有协同效应, 有待研究。

我省登革热重点地区埃及伊蚊野外群体存在V1016G和F1534C突变, V1016G突变广泛存在, F1534C突变以景洪埃及伊蚊种群为最高, 且两种突变同时存在也是景洪的埃及伊蚊种群最高, 从分子生物学的角度证实了景洪市广泛使用和不受控制地频繁使用光稳定的拟除虫菊酯可能是该突变的强大选择压力, 埃及伊蚊对拟除虫菊酯类杀虫剂已经产生了高度抗性, 应注意杀虫剂的合理选用。我省登革热高发地区, 应定期进行蚊虫生物学测定、代谢酶活性和击倒抗性基因突变检测, 从而能够根据分子检测结果预测蚊虫的抗性水平, 时时掌握抗性发生发展动态, 最终实现对蚊虫种群抗性的“ 分子监测” 。

利益冲突:无

引用本文格式:兰学梅, 杨明东, 杨锐, 等.云南省埃及伊蚊对拟除虫菊酯类抗性群体的击倒抗性基因突变分析[J]. 中国人兽共患病学报, 2020, 36(12):993-999. DOI:10.3969/j.issn.1002-2694.2020.00.163

编辑:王晓欢

参考文献
[1] 董学书, 蔡福昌, 周红宁, . 云南省边境口岸蚊类调查[J]. 中国媒介生物学及控制杂志2004, 15(2): 142-145. [本文引用:1]
[2] 李春敏, 董学书, 杨明东. 云南省埃及伊蚊地理分布与季节消长[J]. 中国媒介生物学及控制杂志, 2018, 29(4): 388-392. DOI: 10.11853/j.issn.1003.8280.2018.04.000 [本文引用:1]
[3] 刘永华, 尹小雄, 杨召兰, . 云南省瑞丽市2013年登革热暴发的流行病学分析[J]. 中国媒介生物学及控制杂志, 2014, 25(6): 524-526. DOI: 10.11853/j.issn.1003.4692.2014.06.010 [本文引用:1]
[4] 李华昌, 杨贵荣. 云南临沧市首次本地感染登革热流行病学调查[J]. 中国热带医学, 2009, 9(10): 2013-2014. [本文引用:1]
[5] Wen S, Ma D, Lin Y, et al. Complete genome characterization of the 2017 dengue outbreak in Xishuangbanna, a Border City of China, Burma and Laos[J]. Front Cell Infect Microbiol, 2018, 8: 148. DOI: 10.3389/fcimb.2018.00148 [本文引用:1]
[6] 王建羽, 兰学梅, 杨锐, . 云南省耿马县埃及伊蚊对6种常用卫生杀虫剂的抗药性研究[J]. 中国媒介生物学及控制杂志, 2017, 28(5): 444-446. DOI: 10.11853/j.issn.1003.8280.2017.05.008 [本文引用:1]
[7] 兰学梅, 郑宇婷, 董朝良, . 云南省瑞丽市埃及伊蚊和白纹伊蚊对常用杀虫剂的抗药性调查[J]. 中国媒介生物学及控制杂志, 2017, 28(6): 572-575. DOI: 10.11853/j.issn.1003.8280.2017.06.014 [本文引用:1]
[8] 董学书, 周红宁, 龚正达. 云南蚊类志(下卷)[M]. 昆明: 云南科技出版社, 2010: 62-64. [本文引用:1]
[9] Li C-X, Kaufman PE, Xue R-D, et al. Relationship between insecticide resistance and kdr mutations in the dengue vector Aedes aegypti in Southern China[J]. Parasit Vectors, 2015, 8: 325. DOI: 10.1186/s13071-015-0933-z [本文引用:1]
[10] Wuliand ari JR, Lee SF, White VL, et al. Association between three mutations, F1565C, V1023G and S996P, in the voltage-sensitive sodium channel gene and knockdown resistance in Aedes aegypti from Yogyakarta, Indonesia[J]. Insects, 2015, 6(3): 658-685. DOI: 10.3390/insects6030658 [本文引用:1]
[11] Du Y, Nomura Y, Zhorov BS, et al. Sodium channel mutations and pyrethroid resistance in Aedes aegypti[J]. Insects, 2016, 7(4). DOI: 10.3390/insects7040060 [本文引用:1]
[12] 师灿南. 景洪市登革热媒介伊蚊对常用杀虫剂的抗药性及机制初步研究[D]. 北京: 中国疾病预防控制中心, 2017. [本文引用:1]
[13] Hitoshi Kawada S Z M O. Co-occurrence of point mutations in the voltage-gated sodium channel of pyrethroid-resistant Aedes aegypti populations in Myanmar[J]. PLoS Negl Trop Dis, 2014 Jul 31;8(7): e3032. DOI: 10.1371/journal.pntd.0003032 [本文引用:2]
[14] Stenhouse SA, Plernsub S, Yanola J, et al. Detection of the V1016G mutation in the voltage-gated sodium channel gene of Aedes aegypti(Diptera: Culicidae)by allele-specific PCR assay, and its disturbution and effect on deltamethrin resistance in Thailand [J]. Parasite Vectors, 2013, 6(1): 253. DOI: 10.1186/1756-3305-6-253 [本文引用:1]
[15] Marcombe S, Fustec B, Catte J, et al. Distribution of insecticide resisitance and mechanisms involved in the arbovirus vector Aedes aegypti in Laos and implication for vector control[J]. PLoS Negl Trop Dis, 2019, 13(12): e0007852. DOI: 10.1371/journal.pntd.0007852 [本文引用:1]
[16] Kawada H, Higa Y, Komagata O, et al. Widespread distribution of a newly found point mutation in voltage-gated sodium channel in pyrethroid-resistant Aedes aegypti populations in Vietnam[J]. PLoS Negl Trop Dis, 2009, 3(10): e527. DOI: 10.1371/journal.pntd.0000527 [本文引用:1]
[17] Srisawat R, Komalamisra N, Eshita Y, et al. Point mutations in domain II of the voltage-gated sodium channel gene in deltamethrin-resistant Aedes aegypti (Diptera: Culicidae)[J]. Appl Entomol Zool, 2010, 45(2): 275-282. DOI: 10.1303/aez.2010.275 [本文引用:2]
[18] Srisawat R, Komalamisra N, Apiwathnasorn C, et al. Field-collected permethrin-resistant Aedes aegypti from central Thailand contain point mutations in the domain IIS6 of the sodium channel gene (KDR)[J]. Southeast Asian J Trop Med Public Health, 2012, 43(6): 1380-1386. [本文引用:1]
[19] Chang C, Shen WK, Wang TT, et al. A novel amino acid substitution in a voltage-gated sodium channel is associated with knockdown resistance to permethrin in Aedes aegypti[J]. Insect Biochem Mol Biol, 2009, 39(4): 272-278. DOI: 10.1016/j.ibmb.2009.01.001 [本文引用:1]
[20] Ishak IH, Jaal Z, Ranson H, et al. Contrasting patterns of insecticide resistance and knockdown resistance ( kdr) in the dengue vectors Aedes aegypti and Aedes albopictus from Malaysia[J]. Parasite Vectors, 2015, 8: 181. DOI: 10.1186/s13071-015-0797-2 [本文引用:1]
[21] Kushwah RB, Dykes CL, Kapoor N, et al. Pyrethroid-resistance and presence of two knockdown resistance ( kdr) mutations, F1534C and a novel mutation T1520I, in Indian Aedes aegypti[J]. PLoS Negl Trop Dis, 2015, 9(1): e3332. DOI: 10.1371/journal.pntd.0003332 [本文引用:1]
[22] Chung HH, Cheng IC, Chen YC, et al. Voltage-gated sodium channel intron polymorphism and four mutations comprise six haplotypes in an Aedes aegypti population in Taiwan[J]. PLoS Negl Trop Dis, 2019, 13(3): e0007291. DOI: 10.1371/journal.pntd.0007291 [本文引用:1]