Search

Abstract Background The efficacy of chloroquine treatment for vivax malaria has been rarely evaluated due to a lack of an appropriate testing method. The objective of this study was to conduct molecular monitoring of chloroquine resistance in Plasmodium vivax strains from vivax malaria patients in Yunnan Province, focusing on the analysis of polymorphism in the P. vivax chloroquine resistance transporter protein orthologous gene (pvcrt-o). Methods In accordance with the principles of a cohort study, blood samples were collected from malaria cases diagnosed with a P. vivax mono-infection in Yunnan Province from 2020 to 2022. Segmental PCR was used to amplify the whole pvcrt-o gene in the blood samples and their products were subsequently sequenced. The sequencing data were arranged to obtain the full coding DNA sequence (CDS) as well as the gene’s promoter region sequences. The CDSs were aligned with the reference sequence (XM_001613407.1) of the P. vivax SalI isolate to identify the mutant loci. Results From a total of 375 blood samples taken from vivax malaria cases, 272 both whole gene CDSs (1272–1275 bp) and promoter DNA sequences (707 bp) of pvcrt-o gene were obtained. Among the whole CDSs, there were 7 single nucleotide polymorphic sites in which c.7 A>G was the minor allele frequency (MAF) site with 4.4% (12/272) detection rate. The mutation detection rate showed a significant decrease from 9.8% (10/102) in 2020 to 1.1% (1/92) in 2021 and 1.3% (1/78) in 2022, indicating statistical significance (χ 2 = 11.256, P

Abstract Background Chloroquine (CQ) has been the preferred clinical treatment for vivax malaria in Yunnan Province since 1958, with over 300,000 patients. This study aimed to help make trend predictions regarding variations the in anti-malarial drug susceptibility of Plasmodium vivax distributed in Yunnan Province and effectively implement monitoring measures on the efficacy of anti-malarial drugs for vivax malaria. Methods Blood samples collected from patients with mono-P. vivax infections were employed in this study based on the principle of cluster sampling. The whole gene of P. vivax multidrug resistance 1 protein gene (pvmdr1) was amplified by nested-PCR techniques and the PCR amplification produce were sequenced by Sanger bidirectional sequencing. The mutant loci and haplotypes of coding DNA sequence (CDS) were identified by comparison with the reference sequence (NC_009915.1) of the P. vivax Sal I isolate. Parameters such as Ka/Ks ratio were calculated using MEGA 5.04 software. Results A total of 753 blood samples from patients infected with mono-P. vivax were collected, of which 624 blood samples yielded the full gene sequence (4392 bp) of the pvmdr1 gene, with 283, 140, 119, and 82 sequences from 2014, 2020, 2021 and 2022, respectively. A total of 52 single nucleotide polymorphic (SNP) loci were detected for the 624 CDSs, of which 92.3% (48/52), 34.6% (18/52), 42.3% (22/52), and 36.5% (19/52) SNPs were detected in 2014, 2020, 2021 and 2022, respectively. All of 624 CDSs were defined for a total of 105 mutant haplotypes, with CDSs of 2014, 2020, 2021, and 2022 containing 88, 15, 21, and 13 haplotypes, respectively. Of the 105 haplotypes, the threefold mutant haplotype (Hap_87) was the starting point for stepwise evolution, and the most drastic tenfold mutations were Hap_14 and Hap_78, and the fivefold, sixfold, sevenfold, and eightfold mutations. Conclusions In the majority of vivax malaria cases in Yunnan Province, most of them were infected with strains carrying demonstrating highly mutated in pvmdr1 genes. However, the dominant mutation strains types varied from year to year, which warrants further exploration in order to confirm the correlation between with phenotypic changes in P. vivax strains and their susceptibility to anti-malarial drugs such as chloroquine.

Abstract Background Anopheles maculatus, Anopheles minimus and Anopheles dirus are the major vectors of malaria transmission in the Greater Mekong Subregion (GMS). The malaria burden in this region has decreased significantly in recent years as all GMS countries progress towards malaria elimination. It is necessary to investigate the Anopheles diversity and abundance status and assess the Plasmodium infection rates to understand the malaria transmission potential of these vector species in GMS countries to guide the development of up-to-date vector control strategies and interventions. Methods A survey of mosquitoes was conducted in Stung Treng, Sainyabuli and Phongsaly Provinces on the Cambodia-Laos, Thailand-Laos and China-Laos borders, respectively. Mosquito collection was done by overnight trapping at sentinel sites in each province. After morphological identification, the 18S rRNA-based nested-PCR was performed to detect malaria parasites in the captured Anopheles mosquitoes. Results A total of 18 965 mosquitoes comprising of 35 species of 2 subgenera (Subgenus Anopheles and Subgenus Cellia) and 4 tribes (Tribes Culicini, Aedini, Armigerini and Mansoniini) were captured. Tribe Culicini accounted for 85.66% of captures, followed by Subgenus Anopheles (8.15%). Anopheles sinensis dominated the Subgenus Anopheles by 99.81%. Plasmodium-infection was found in 25 out of the 1 683 individual or pooled samples of Anopheles. Among the 25 positive samples, 19, 5 and 1 were collected from Loum, Pangkhom and Siem Pang village, respectively. Eight Anopheles species were found infected with Plasmodium, i.e., An. sinensis, Anopheles kochi, Anopheles vagus, An. minimus, Anopheles annularis, Anopheles philippinensis, Anopheles tessellatus and An. dirus. The infection rates of Plasmodium falciparum, Plasmodium vivax and mixture of Plasmodium parasite species were 0.12% (2/1 683), 1.31% (22/1 683) and 0.06% (1/1 683), respectively. Conclusions Overall, this survey re-confirmed that multiple Anopheles species carry malaria parasites in the international border areas of the GMS countries. Anopheles sinensis dominated the Anopheles collections and as carriers of malaria parasites, therefore may play a significant role in malaria transmission. More extensive investigations of malaria vectors are required to reveal the detailed vector biology, ecology, behaviour, and genetics in GMS regions in order to assist with the planning and implementation of improved malaria control strategies. Graphical Abstract

Abstract Background Few studies have been conducted to investigate the distribution of mosquito vectors and the population structure of secondary vectors in the border region of Cambodia-Laos. The aim of this work was to study the mosquito diversity and molecular phylogeny of secondary vectors, i.e., Anopheles nivipes in this area. Methods 1440 adult mosquitoes were trapped in the Cambodia-Laos border. mtDNA-COII were amplified and sequenced from 53 An. nivipes DNA samples. Together with COII sequences deposited in GenBank, a total of 86 COII sequences were used for examining population variations, genetic differentiation, spatial population structure, population expansion, and gene flow patterns. Results The adult mosquitoes were classified into 5 genera and 27 species in this border region. The predominant genera were Culex (60.07%, 865/1440) and Anopheles (31.25%, 450/1440), and the major Anopheles species were An. nivipes (73.56%, 331/450) and Anopheles maculatus (14.22%, 64/450). Based on sequences analysis of COII, a high level of genetic differentiation was reported in two Northwest India (Cheema and Bathinda, Punjab) and Cambodia-Laos (Siem Pang, Stung treng) populations (F ST = 0.97824, 0.97343, P

Abstract Background In recent years, the incidence rate of vivax malaria recurrence still had 3.1% in Yunnan Province population after eradication therapy using primaquine (PQ). In order to understand the specific failure reasons for preventing vivax malaria relapses, a preliminary exploration on the CYP2D6 enzyme activity was carried out in the vivax malaria patients in Yunnan Province population by analysing mutational polymorphism in the coding region of CYP2D6 gene. Methods Blood samples were collected from vivax malaria patients with suspected relapse (SR) and non-relapsed (NR) malaria in Yunnan Province. The DNA fragments containing 9 exons regions of human CYP2D6 gene were amplified by performing PCR and sequenced. The sequencing results were aligned by using DNAStar 11.0 to obtain the coding DNA sequence (CDS) of CYP2D6 gene. DnaSP 6.11.01 software was used to identify mutant polymorphisms and haplotypes of the CDS chain. The waterfall function of GenVisR package in R was utilized to visualize the mutational landscape. The alleles of CYP2D6 gene were identified according to the criteria prescribed by Human Cytochrome P450 (CYP) Allele Nomenclature Committee Database and the CYP2D6 enzyme activity was predicted based on diploid genotype. Results A total of 320 maternal CDS chains, including 63 from SR group and 257 from NR group, were obtained. Twelve mutant loci, including c.31 (rs769259), c.100 (rs1065852), c.271 (rs28371703), c.281 (rs28371704), c.294 (rs28371705), c.297 (rs200269944), c.336 (rs1081003), c.408 (rs1058164), c.505 (rs5030865), c.801 (rs28371718), c.886 (rs16947), and c.1,457 (rs1135840) were observed on the 640 CDS chains (including 320 maternal and 320 paternal chains). The high-frequency mutation at rs1135840 (0.703) and low-frequency mutation, such as rs28371703, were detected only in the SR group. The frequency of mutant rs1058164 and rs1135840 were significantly increased in the SR group ( $${x}^{2}$$ x 2 = 4.468, 5.889, P

Abstract Background Accumulating evidence suggest that compromised CYP2D6 enzyme activity caused by gene mutation could contribute to primaquine failure for the radical cure of vivax malaria. The current study aims to preliminarily reveal the association between the recurrence of vivax malaria in Yunnan Province and CYP2D6 gene mutation by analysing polymorphisms in the entire coding region of human CYP2D6 gene. Methods Blood samples were collected from patients with vivax malaria, who received "chloroquine and 8-day course of primaquine therapy" in Yunnan Province. The suspected relapsed cases were determined by epidemiological approaches and gene sequence alignment. PCR was conducted to amplify the CYP2D6 gene in the human genome, and the amplified products were then sequenced to compare with the non-mutation “reference” sequence, so as to ensure correct sequencing results and to determine 9 exon regions. Subsequently, the DNA sequences of 9 exons were spliced into the coding DNA sequence (CDS), which, by default, is known as maternal CDS. The paternal CDS was obtained by adjusting the bases according to the sequencing peaks. The mutation loci, haplotypes (star alleles), genotypes and odds ratios (OR) of all the CDSs were analysed. Results Of the119 maternal CDS chains in total with 1491 bp in length, 12 mutation sites in the 238 maternal and paternal CDS chains were detected. The c.408G > C mutation was most frequently detected in both suspected relapsed group (SR) and non-relapsed group (NR), reaching 85.2% (75/88) and 76.0% (114/150), respectively. The c.886C > T mutation was most closely related to the recurrence of vivax malaria (OR = 2.167, 95% CI 1.104–4.252, P T and CYP2D6*2 allele, which correspond to impaired PQ metabolizer phenotype, are most closely related to the relapse of vivax malaria. In addition, the genotype *4/*4 with null CYP2D6 enzyme function was only detected in the SR group. These results reveal the risk of defected CYP2D6 enzyme activity that diminishes the therapeutic effect of primaquine on vivax malaria.

Abstract Background According to China’s Malaria Eradication Action Plan, malaria cases diagnosed and reported by health authorities at the county level must be further re-confirmed by provincial laboratories. The Yunnan Province Malaria Diagnostic Reference Laboratory (YPMDRL) began the synchronous implementation of microscopic examinations and nested polymerase chain reaction (nested-PCR) testing to re-test the malaria cases initially diagnosed by county-level laboratories and to evaluate the consistency of Plasmodium species identified between by YPMDRL and by the county-level laboratories from 2013 to 2018 in Yunnan Province. Methods Data on malaria initial diagnosis completed by county-level laboratories in Yunnan Province were collected weekly from the “China Disease Prevention and Control Information System” from 2013 to 2018. The YPMDRL performed Plasmodium microscopic examination and 18S rRNA gene nested-PCR testing on every malaria case managed by the China Disease Prevention and Control Information System. The re-testing detection results were fed back to the initial diagnosis and reporting unit for revision of malaria case types. Results A total of 2,869 malaria cases were diagnosed and reported by county-level laboratories in Yunnan Province from 2013 to 2018. The re-testing rate was 95.6% (2,742/2,869), and the re-testing rate increased from 2013 to 2018. Among the re-tested 2,742 cases, 96.7% (2651/2742), 2.2% (59/2742), and 1.1% (32/2742) were doubly examined by microscopy and by nested-PCR, only by microscopy, and only by nested-PCR, respectively. The total Plasmodium species accuracy rate at county-level laboratories was 92.6% (2,543/2,742) reference to the diagnosis by YPMDRL. Among the inconsistent 199 cases, they were identified as including 103 negative cases, 45 falciparum malaria cases, 30 vivax malaria cases, 11 ovale malaria cases, and 10 malariae malaria cases by YPMDRL. From 2013 to 2018, the revised and registered malaria cases by the China Disease Prevention and Control Information System in Yunnan Province was 2,747 cases, including 2,305 vivax malaria cases, 421 falciparum malaria cases, 11 ovale malaria cases, and 10 malariae malaria cases. Conclusions The double re-testing strategy by microscopy and by gene testing increases the accuracy of diagnoses malaria in Yunnan Province, and gene testing can reliably differentiate Plasmodium species. The re-testing results provided by YPMDRL are the authoritative basis for revising malaria kind in Yunnan Province.

Abstract Background Eighteen imported ovale malaria cases imported from Myanmar and various African countries have been reported in Yunnan Province, China from 2013 to 2018. All of them have been confirmed by morphological examination and 18S small subunit ribosomal RNA gene (18S rRNA) based PCR in YNRL. Nevertheless, the subtypes of Plasmodium ovale could not be identified based on 18S rRNA gene test, thus posing challenges on its accurate diagnosis. To help establish a more sensitive and specific method for the detection of P. ovale genes, this study performs sequence analysis on k13-propeller polymorphisms in P. ovale. Methods Dried blood spots (DBS) from ovale malaria cases were collected from January 2013 to December 2018, and the infection sources were confirmed according to epidemiological investigation. DNA was extracted, and the coding region (from 206th aa to 725th aa) in k13 gene propeller domain was amplified using nested PCR. Subsequently, the amplified products were sequenced and compared with reference sequence to obtain CDS. The haplotypes and mutation loci of the CDS were analysed, and the spatial structure of the amino acid peptide chain of k13 gene propeller domain was predicted by SWISS-MODEL. Results The coding region from 224th aa to 725th aa of k13 gene from P. ovale in 83.3% of collected samples (15/18) were amplified. Three haplotypes were observed in 15 samples, and the values of Ka/Ks, nucleic acid diversity index (π) and expected heterozygosity (He) were 3.784, 0.0095, and 0.4250. Curtisi haplotype, Wallikeri haplotype, and mutant type accounted for 73.3% (11/15), 20.0% (3/15), and 6.7% (1/15). The predominant haplotypes of P. ovale curtisi were determined in all five Myanmar isolates. Of the ten African isolates, six were identified as P. o. curtisi, three were P. o. wallikeri and one was mutant type. Base substitutions between the sequences of P. o. curtisi and P. o. wallikeri were determined at 38 loci, such as c.711. Moreover, the A > T base substitution at c.1428 was a nonsynonymous mutation, resulting in amino acid variation of T476S in the 476th position. Compared with sequence of P. o. wallikeri, the double nonsynonymous mutations of G > A and A > T at the sites of c.1186 and c.1428 leads to the variations of D396N and T476S for the 396th and 476th amino acids positions. For P. o. curtisi and P. o. wallikeri, the peptide chains in the coding region from 224th aa to 725th aa of k13 gene merely formed a monomeric spatial model, whereas the double-variant peptide chains of D396N and T476S formed homodimeric spatial model. Conclusion The propeller domain of k13 gene in the P. ovale isolates imported into Yunnan Province from Myanmar and Africa showed high differentiation. The sequences of Myanmar-imported isolates belong to P. o. curtisi, while the sequences of African isolates showed the sympatric distribution from P. o. curtisi, P. o. wallikeri and mutant isolates. The CDS with a double base substitution formed a dimeric spatial model to encode the peptide chain, which is completely different from the monomeric spatial structure to encode the peptide chain from P. o. curtisi and P. o. wallikeri.