中国寄生虫学与寄生虫病杂志 ›› 2023, Vol. 41 ›› Issue (4): 404-411.doi: 10.12140/j.issn.1000-7423.2023.04.002

• 论著 • 上一篇    下一篇

云南省输入性间日疟原虫多药抗性蛋白1基因突变多态性分析

丁红芸1(), 董莹2,*(), 徐艳春2, 邓艳2, 刘言2, 吴静2, 陈梦妮2, 张苍林2   

  1. 1 昆明理工大学生命科学与技术学院,云南昆明 650500
    2 云南省寄生虫病防治所,云南省疟疾研究中心,云南省虫媒传染病防控研究重点实验室,普洱 665000
  • 收稿日期:2022-10-13 修回日期:2023-01-27 出版日期:2023-08-30 发布日期:2023-09-06
  • 通讯作者: *董莹(1964-),女,本科,主任医师,从事疟疾病原学及分子生物学研究。E-mail:luxidongying@126.com
  • 作者简介:丁红芸(1998-),女,硕士研究生,从事疟疾的分子流行病学研究。E-mail:1053343328@qq.com
  • 基金资助:
    国家自然科学基金(81660559);国家自然科学基金(82160637)

Polymorphism analysis of multidrug resistance protein 1 gene in imported Plasmodium vivax in Yunnan Province

DING Hongyun1(), DONG Ying2,*(), XU Yanchun2, DENG Yan2, LIU Yan2, WU Jing2, CHEN Mengni2, ZHANG Canglin2   

  1. 1 Faculty of Life Science and Technology, Kunming University of Science and Technology, Kunming 650500, Yunnan, China
    2 Yunnan Institute of Parasitic Diseases, Yunnan Center of Malaria Research, Yunnan Provincial Key Laboratory of Vector-borne Diseases Control and Research, Puer 665000, China
  • Received:2022-10-13 Revised:2023-01-27 Online:2023-08-30 Published:2023-09-06
  • Contact: *E-mail: luxidongying@126.com
  • Supported by:
    National Natural Science Foundation of China(81660559);National Natural Science Foundation of China(82160637)

摘要:

目的 分析云南省输入性间日疟原虫的多药抗性蛋白1基因(Pvmdr1)突变多态性。 方法 对云南省2020年和2021年诊断为间日疟原虫感染的病例血样进行Pvmdr1全基因扩增及测序,以间日疟原虫Sal-Ⅰ分离株的序列(GenBank登录号:NC_009915.1)为模板设计PCR反应引物,用MEGA 5.04软件与编码区参比序列(GenBank登录号:XM_001613678.1)进行比对,用DnaSP 6.11.01软件分析单倍型和单核苷酸多态(SNP)位点及其突变类型,并计算核酸多样性指数π、单倍型期望杂合度(He)等,用Network 10.0软件构建中介网络进化图。 结果 共收集276份输入性间日疟患者血样,其中自缅甸感染病例占99.3%(274/276),自非洲、巴基斯坦感染各占0.4%(1/276)。259份血样扩增获得4 392 bp的Pvmdr1全基因序列,提交GenBank获得的登录号为OP559204~OP559462。259条DNA序列π、Ka/Ks比值分别为0.000 76和3.600 6,共检出22个SNP,包括13个非同义突变和9个同义突变。仅c.4074C>T位点突变在2020、2021年的检出率[15.0%(21/140)、5.9%(7/119)]差异有统计学意义(χ2 = 5.546,P < 0.05)。次等位位点为c.1587A>G(97.9%,253/259),新检出SNP包括c.2499G>T(2.3%,6/259)、c.3358C>T(0.4%,1/259)、c.3832C>T(0.4%,1/259)。c.3064C>T、c.4065A>G、c.3358C>T和c.3832C>T突变位点仅在2021年的患者血样中检出,其中前2个SNP从非洲感染的分离株中检出,后2个SNP从缅甸感染的分离株中检出。259条Pvmdr1完整编码区与参比序列(单倍型Hap_1)比对,识别出26种单倍型(Hap_2~Hap_27),均为参比序列(单倍型Hap_1)的突变型,He为0.890 7。其中,Hap_8的检出率最高(18.5%,48/259),非洲感染分离株中仅检出Hap_19。2020、2021年的患者血样中检出的单倍型种类分别有15和22种,当年独有的单倍型分别为4和10种。Hap_10在2020、2021年患者血样中的检出率分别为15.0%(21/140)和3.4%(4/119),差异有统计学意义(χ2 = 22.264,P < 0.05)。不同单倍型的多重突变识别结果显示,4重、5重、6重、7重、9重、10重突变在26种单倍型中的占比分别为0.4%(1/26)、19.2%(5/26)、38.5%(10/26)、30.8%(8/26)、0.4%(1/26)、0.4%(1/26)。中介网络图显示,26种单倍型以单倍型Hap_1为起始,经4重突变(Hap_24),向5重、6重、7重、9重和10重突变逐级进化并远离起始点。 结论 云南省输入性间日原虫Pvmdr1基因突变存在高度多态性。

关键词: 间日疟原虫, 多药抗性蛋白1, 突变, 云南省

Abstract:

Objective To analyze the polymorphism of Plasmodium vivax multidrug resistance protein 1 gene (Pvmdr1) from imported vivax malaria patients in Yunnan Province. Methods The whole Pvmdr1 gene was amplified and sequenced in the blood samples of patients who diagnosed with vivax malaria in 2020 and 2021 in Yunnan Province. The sequence of P. vivax SalⅠ isolate (GenBank accession number: NC_009915.1) was taken as the reference sequence for designing primers. MEGA 5.04 software was used to align the spliced sequences with the coding region reference sequence (GenBank accession number: XM_001613678.1). Haplotype and single nucleotide polymorphism (SNP) sites and their mutation types, nucleotide diversity π, expected heterozygosity (He), etc were analyzed by DnaSP 6.11.01. Haplotype medium network diagram was constructed using Network 10.0 Software. Results A total of 276 blood samples were collected from vivax malaria patients and were imported from overseas. The full gene sequences with 4 392 bp of Pvmdr1 were obtained from 259 blood samples, and the accession numbers were OP559204-OP559462 assigned by GenBank. The π of 259 DNA coding sequences was 0.000 76. A total of 22 SNPs were detected from 259 coding sequences, including 13 non-synonymous mutant and 9 synonymous mutant loci. Only the detection rate difference of c.4074C>T mutation between 2020 and 2021 [15.0% (21/140), 5.9% (7/119)] was statistically significant (χ2 = 5.546, P < 0.05). The minor allele frequency was c.1587A>G (97.9%, 253/259), and the newly discovered SNPs included c.2499G>T (2.3%, 6/259), c.3358C>T (0.4%, 1/259), and c.3832C>T (0.4%, 1/259). All of c.3064C>T, c.4065A>G, c.3358C>T and c.3832C>T mutant sites were detected from the blood samples in 2021. The 259 complete coding region sequences of Pvmdr1 were defined as 26 haplotypes (Hap_2-Hap_27) by comparison with the reference sequence (haplotype Hap_1, all of which were mutant type haplotypes, and He was 0.890 7. Among them, the detection rate of Hap_8 was the highest (18.5%, 48/259). There were 15 and 22 haplotypes found from blood samples in 2020 and 2021, respectively. Among them, 4 haplotypes were only detected in blood samples from 2020, and 10 haplotype were only detected in blood samples from 2021. The difference of Hap_10 detection rate between 2020 and 2021 [15.0% (21/140), 3.4% (4/119)] was statistically significant (χ2 = 22.264, P < 0.05). The proportions of quadruple, quintuple, six-fold, seven-fold, nine-fold and ten-fold multiple mutations in 26 haplotypes were 0.4% (1/26), 19.2% (5/26), 38.5% (10/26), 30.8% (8/26), 0.4% (1/26), and 0.4% (1/26), respectively. The haplotype medium network diagram displayed that the gradual evolution of 26 haplotypes started with Hap_1, and then passed through quadruple mutation (Hap_24), moved towards quintuple, six-fold, seven-fold, nine-fold, and ten-fold multiple mutations away from. Conclusion There is high polymorphism in the Pvmdr1 gene in the P. vivax from the imported vivax malaria cases in Yunnan Province.

Key words: Plasmodium vivax, Multidrug resistance protein 1, Mutation, Yunnan Province

中图分类号: