Search

Myanmar aims to eliminate malaria by 2030. However, recent increase of malaria incidence is a great challenge to archive that goal. Increasing prevalence of Plasmodium vivax also hinders this endeavor. Monitoring genetic structure of the parasite is necessary to understand genetic nature and evolutionary aspect of P. vivax population in Myanmar. Partial fragment flanking blocks I and II of merozoite surface protein-3 alpha of P. vivax (pvmsp-3α) was amplified from P. vivax isolates collected in Pyin Oo Lwin, Mandalay Region, Myanmar in 2013-2015. Sequence analysis of pvmsp-3α was performed to determine genetic diversity and natural selection of this gene. Spatio-temporal genetic changes of pvmsp-3α in Myanmar P. vivax population were also investigated via comparative analysis of gene sequences obtained in this study and previously reported Myanmar pvmsp-3α sequences. Genetic diversity of Myanmar pvmsp-3α was detected in P. vivax isolates analyzed. Size polymorphisms in block I and amino acid changes and recombination events in block II were main factors contributing to the genetic diversity of pvmsp-3α. Comparative spatio-temporal analysis with previously reported Myanmar pvmsp-3α populations revealed the presence of genetic differences by population with moderate genetic differentiation between populations. Similar pattern of natural selection was also detected in Myanmar pvmsp-3α populations. These suggested that enough size of the P. vivax population sufficient to generate or maintain the genetic diversity remains in the population. Thus, continuous molecular surveillance of genetic structure of Myanmar P. vivax is necessary., Competing Interests: Declaration of competing interest The authors declare that they have no competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023. Published by Elsevier B.V.)

Background: Polymorphisms of Plasmodium falciparum chloroquine resistance transporter (pfcrt), Plasmodium falciparum multi-drug resistance 1 (pfmdr1) and Plasmodium falciparum kelch 13-propeller (pfk13) genes are accepted as valid molecular markers of quinoline antimalarials and artemisinins. This study investigated the distribution patterns of these genes in P. falciparum isolates from the areas along the Thai-Myanmar border during the two different periods of antimalarial usage in Thailand., Results: Polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) were used to detect pfcrt mutations at codons 76, 220, 271, 326, 356, and 371 as well as pfmdr1 mutation at codon 86. The prevalence of pfcrt mutations was markedly high (96.4-99.7%) in samples collected during both periods. The proportions of mutant genotypes (number of mutant/total isolate) at codons 76, 220, 271, 326, 356 and 371 in the isolates collected during 1993-1998 (period 1) compared with 2002-2008 (period 2) were 97.9% (137/140) vs. 97.1% (401/413), 97.9% (140/143) vs. 98.8% (171/173), 97.2% (139/143) vs. 97.1% (333/343), 98.6% (140/142) vs. 99.7% (385/386), 96.4% (134/139) vs. 98.2% (378/385) and 97.8% (136/139) vs. 98.9% (375/379), respectively. Most isolates carried pfmdr1 wild-type at codon 86, with a significant difference in proportions genotypes (number of wild type/total sample) in samples collected during period 1 [92.9% (130/140)] compared with period 2 [96.9% (379/391)]. Investigation of pfmdr1 copy number was performed by real-time PCR. The proportions of isolates carried 1, 2, 3 and 4 or more than 4 copies of pfmdr1 (number of isolates carried correspondent copy number/total isolate) were significantly different between the two sample collecting periods (65.7% (90/137) vs. 87.8% (390/444), 18.2% (25/137) vs. 6.3%(28/444), 5.1% (7/137) vs. 1.4% (6/444) and 11.0% (15/137) vs. 4.5% (20/444), for period 1 vs. period 2, respectively). No pfk13 mutation was detected by nested PCR and nucleotide sequencing in all samples with successful analysis (n = 68)., Conclusions: The persistence of pfcrt mutations and pfmdr1 wild-types at codon 86, along with gene amplification in P. falciparum, contributes to the continued resistance of chloroquine and mefloquine in P. falciparum isolates in the study area. Regular surveillance of antimalarial drug resistance in P. falciparum, incorporating relevant molecular markers and treatment efficacy assessments, should be conducted., (© 2024. The Author(s).)

Plasmodium vivax malaria cases remain high along the Thai-Myanmar and Thai-Cambodia borders. Plasmodium vivax circumsporozoite protein (pvcsp) and Plasmodium vivax ookinete surface protein (pvs25) genes are promising molecular markers of the genetic diversity of P. vivax. This study investigated the genetic diversity of pvcsp and pvs25 in P. vivax isolates collected from the Thai-Myanmar border. The DNA samples were amplified, and the genotypes were analyzed by PCR-RFLP and DNA sequencing. Pvcsp genotypes, VK210, VK247, and mixed types, were found in 203 (91.9%), 15 (6.8%), and 3 (1.3%) of the isolates, respectively. Twenty-four allelic variants were observed, of which a high prevalence of VK210E and VK247E were reported. Two pvcsp variants, VK210C and VK210M showed significantly higher parasite density (46,234 (1154-144,000) vs. 25,606 (1373-68,878), respectively). The genetic diversity of pvcsp along the Thai-Myanmar border during 2002-2015 showed dynamic changes with both positive and negative selection. The frequency and distribution of pvcsp pattern might be changed over time and might be other factors contributing to gene selection. Three amino acid substitutions of pvs25, i.e., E97Q, I130T, and Q131K, were investigated with frequencies of 10 (4.5%), 221 (100%), and 204 (92.3%) isolates, respectively. There was no association between parasite density and pvs25 polymorphisms. The frequency of pvs25 polymorphism was similar to that previously reported, with the absence of random mutation. In conclusion, the genetic variation of pvcsp was changed over times whereas the genetic diversity of pvs25 was limited; these variations would be helpful for further vaccine development against P. vivax malaria., (© 2024. The Author(s), under exclusive licence to Springer-Verlag GmbH Germany, part of Springer Nature.)

Backgrounds: The therapeutic efficacy studies of DHA-PIP for uncomplicated Plasmodium falciparum patients were implemented from 2012 to 2016 along China (Yunnan province)-Myanmar border, which verified the high efficacy of DHA-PIP. With the samples collected in these studies, the genetic characteristics of P. falciparum parasites based on in vivo parasite clearance time (PCT) was investigated to explore if these parasites had developed resistance to DHA and PIP at molecular level., Methods: The genetic characteristics were investigated based on K13 genotypes, copy numbers of genes pfpm2 and pfmdr1, and nine microsatellite loci (Short Tandem Repeats, STR) flanking the K13 gene on chromosome 13. The PCT 50s were compared based on different K13 genotypes, sites, periods and copy numbers., Results: In the NW (North-West Yunnan province bordering with Myanmar) region, F446I was the main K13 genotype. No significant differences for PCT 50s presented among three K13 genotypes. In SW (South-West Yunnan province bordering with Myanmar) region, only wild K13 genotype was detected in all parasite isolates whose PCT 50s was significantly longer than those in NW region. For the copy numbers of genes, parasite isolates containing multiple copies of pfmdr1 gene were found in both regions, but only single copy of pfpm2 gene was detected. Though the prevalence of parasite isolates with multiple copies of pfmdr1 gene in SW region was higher than that in NW region, no difference in PCT 50s were presented between isolates with single and multiple copies of pfmdr1 gene. The median He values of F446I group and Others (Non-F446I K13 mutation) group were 0.08 and 0.41 respectively. The mean He values of ML group (Menglian County in SW) and W (wild K13 genotype in NW) group were 0 and 0.69 respectively. The mean Fst values between ML and W groups were significantly higher than the other two K13 groups., Conclusions: P. falciparum isolates in NW and SW regions had very different genetic characteristics. The F446I was hypothesized to have independently appeared and spread in NW region from 2012 and 2016. The high susceptibility of PIP had ensured the efficacy of DHA-PIP in vivo. Multiple copy numbers of pfmdr1 gene might be a potential cause of prolonged clearance time of ACTs drugs along China-Myanmar border., Trial Registration: Trial registration: ISRCTN, ISRCTN 11775446. Registered 17 April 2020-Retrospectively registered, the registered name was Investigating resistance to DHA-PIP for the treatment of Plasmodium falciparum malaria and chloroquine for the treatment of Plasmodium vivax malaria in Yunnan, China. http://www.isrctn.com/ISRCTN11775446., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2023 Li et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Background: Sexual stage surface antigens are potential targets of transmission-blocking vaccines (TBVs). The gametocyte and gamete surface antigen P230, a leading TBV candidate, is critical for red blood cell binding during exflagellation and subsequent oocyst development. Here, the genetic diversity of Pvs230 was studied in Plasmodium vivax parasite isolates from the China-Myanmar border (CMB) and central Myanmar., Methods: Plasmodium vivax isolates were collected in clinics from malaria-endemic areas of the CMB (143 samples) and Myanmar (23 samples). The interspecies variable part (IVP, nucleotides 1-807) and interspecies conserved part (ICP, 808-2862) of Pvs230 were amplified by PCR and sequenced. Molecular evolution studies were conducted to evaluate the genetic diversity, signature of selection, population differentiation, haplotype network, and population structure of the study parasite populations and publicly available Pvs230 sequences from six global P. vivax populations., Results: Limited genetic diversity was observed for the CMB (π = 0.002) and Myanmar (π = 0.001) isolates. Most amino acid substitutions were located in the IVP and cysteine-rich domain of Pvs230. Evidence of positive selection was observed for IVP and purifying selection for ICP. Codon-based tests identified specific codons under natural selection in both IVP and ICP. The fixation index (F ST ) showed low genetic differentiation between East and Southeast Asian populations, with F ST ranging from 0.018 to 0.119. The highest F ST value (F ST  = 0.503) was detected between the Turkey and Papua New Guinea populations. A total of 92 haplotypes were identified in global isolates, with the major haplotypes 2 and 9 being the most abundant and circulating in East and Southeast Asia populations. Several detected non-synonymous substitutions were mapped in the predicted structure and B-cell epitopes of Pvs230., Conclusions: We detected low levels of genetic diversity of Pvs230 in global P. vivax populations. Geographically specific haplotypes were identified for Pvs230. Some mutations are located within a potential B-cell epitope region and need to be considered in future TBV designs., (© 2022. The Author(s).)

Background: Plasmodium vivax remains the predominant species at the China-Myanmar border, imposing a major challenge to the recent gains in regional malaria elimination. To closely supervise the emerging of drug resistance in this area, we surveyed the variations in genes potentially correlated with drug resistance in P. vivax parasite and the possible drug selection with time., Methods: A total of 235 P. vivax samples were collected from patients suffering uncomplicated malaria at Yingjiang, Tengchong, and Longling counties, and Nabang port in China, Yunnan province, and Laiza sub-township in Myanmar, from 2008 to 2017. Five potential drug resistance genes were amplified utilizing nested-PCR and analyzed, including pvdhfr, pvdhps, pvmdr1, pvcrt-o, and pvk12. The Pearson's Chi-squared test or Fisher's exact test were applied to determine the statistical frequency differences of mutations between categorical data., Results: The pvdhfr F57I/L, S58R, T61M and S117T/N presented in 40.6%, 56.7%, 40.1%, and 56.0% of the sequenced P. vivax isolates, and these mutations significantly decreased with years. The haplotype formed by these quadruple mutations predominated in Yingjiang, Tengchong, Longling and Nabang. While a mutation H99S/R (56.6%) dominated in Laiza and increased with time. In pvdhps, the A383G prevailed in 69.2% of the samples, which remained the most prevalent haplotype. However, a significant decrease of its occurrence was also noticed over the time. The S382A/C and A553G existed in 8.4% and 30.8% of the isolates, respectively. In pvmdr1, the mutation Y976F occurred at a low frequency in 5/232 (2.2%), while T958M was fixed and F1076L was approaching fixed (72.4%). The K10 insertion was detected at an occurrence of 33.2% in pvcrt-o, whereas there was no significant difference among the sites or over the time. No mutation was identified in pvk12., Conclusions: Mutations related with resistance to antifolate drugs are prevalent in this area, while their frequencies decrease significantly with time, suggestive of increased susceptibility of P. vivax parasite to antifolate drugs. Resistance to chloroquine (CQ) is possibly emerging. However, since the molecular mechanisms underneath CQ resistance is yet to be better understood, close supervision of clinical drug efficiency and continuous function investigation is urgently needed to alarm drug resistance., (© 2022. The Author(s).)

Background: The fast-declining clinical efficacy of dihydroartemisinin-piperaquine (DHA-PPQ) in Cambodia is a warning of the underlying westward dissemination of piperaquine resistance in the Greater Mekong Subregion (GMS). Mutations in the Plasmodium falciparum Kelch 13-propeller (PfK13) and the P. falciparum chloroquine resistance transporter (PfCRT), as well as plasmepsin 2/3 gene amplification, have been discovered as molecular markers for predicting DHA-PPQ treatment failure. Determining whether these genetic variations of P. falciparum are linked to DHA-PPQ resistance is critical, especially along the China-Myanmar (CM) border, where PPQ has been utilized for decades., Methods: A total of 173 P. falciparum samples of dried blood spots (DBS) were collected along the CM border between 2007 and 2010, the Thailand-Cambodia (TC) border between 2009 and 2013, and the Thailand-Myanmar (TM) border between 2012 and 2014. PCR and sequencing were used to identified PfCRT mutations, while qPCR was used to determine the copy number of plasmepsin 2/3. The prevalence of DHA-PPQ resistance in three locations was investigated using data paired with K13 mutations., Results: Three fragments of the pfcrt gene were amplified for all 173 samples, and seven SNPs were identified (M74I, N75E/D, K76T, H97L, I218F, A220S, I356L). No new PfCRT mutations conferring resistance to PPQ (T93S, H97Y, F145I, M343L, and G353V) were discovered, except for one mutant I218F identified in the TM border (2.27%, 1/44). Additionally, mutant H97L was found in the TC, TM, and CM borders at 3.57% (1/28), 6.82% (3/44), and 1% (1/101), respectively. A substantial K13 C580Y variant prevalence was found in the TC and TM border, accounting for 64.29% (18/28) and 43.18% (19/44), respectively, while only 1% (1/101) was found in the CM border. The K13 F446I variant was only identified and found to reach a high level (28.71%, 29/101) in the CM border. Furthermore, 10.71% (3/28) of TC isolates and 2.27% (1/44) of TM isolates carried more than one copy of plasmepsin 2/3 and K13 C580Y variant, while no plasmepsin 2/3 amplification was identified in the CM isolates., Conclusions: Compared with the P. falciparum samples collected from the TC and TM borders, fewer parasites carried plasmepsin 2/3 amplification and novel PfCRT variants, while more parasites carried predominant K13 mutations at position F446I, in the CM border. Clear evidence of DHA-PPQ resistance associated with candidate markers was not found in this border region suggesting a further evaluation of these markers and continuous surveillance is warranted., (© 2022. The Author(s).)

Background: Anti-malarial drug resistance is still a major threat to malaria elimination in the Great Mekong Sub-region. Plasmodium vivax parasites resistant to anti-malarial drugs are now found in Myanmar. Molecular surveillance on drug resistance genes in P. vivax parasites from northeastern Myanmar was aimed at estimating the underlying drug resistance in this region., Methods: Blood samples from patients with vivax malaria were collected from Laiza city in northeastern Myanmar in 2020. Drug resistance genes including Pvcrt-o, Pvmdr1, Pvdhfr and Pvdhps were amplified and sequenced. Genetic polymorphisms and haplotypes were analysed to evaluate the prevalence of mutant alleles associated with drug resistance., Results: A total of 149 blood samples from P. vivax patients were collected. The prevalence of Pvmdr1 mutations at codons 958 and 1076 was 100.0% and 52.0%, respectively, whereas no single nucleotide polymorphism was present at codon 976. The proportions of single and double mutant types were 48.0% and 52.0%, respectively. A K10 "AAG" insertion in the Pvcrt-o gene was not detected. Mutations in Pvdhfr at codons 57, 58, 61, 99 and 117 were detected in 29.9%, 54.3%, 27.6%, 44.9% and 55.1% of the samples, respectively. Wild type was predominant (46.3%), followed by quadruple and double mutant haplotypes. Of three types of tandem repeat variations of Pvdhfr, Type B, with three copies of GGDN repeats, was the most common. Pvdhps mutations were only detected at codons 383 and 553 and the wild type Pvdhps was dominant (78.0%). Eleven haplotypes were identified when combining the mutations of Pvdhfr and Pvdhps, among which the predominant one was the wild type (33.9%), followed by double mutant alleles S58R/S117N /WT (24.6%)., Conclusions: This study demonstrated resistant P. vivax phenotypes exists in northeastern Myanmar. Continued surveillance of drug resistance markers is needed to update treatment guidelines in this region., (© 2022. The Author(s).)

Plasmodium vivax apical membrane antigen-1 (PvAMA-1) is an important vaccine candidate for vivax malaria. However, antigenic variation within PvAMA-1 is a major obstacle to the design of a global protective malaria vaccine. In this study, we analyzed the genetic polymorphism and selection of the PvAMA-1 gene from 152 P. vivax isolates from imported cases to China, collected in the China-Myanmar border (CMB) area in Yunnan Province (YP) during 2009-2011 ( n = 71) and 2014-2016 ( n = 81), in comparison with PvAMA-1 gene information from Myanmar ( n = 73), collected from public data. The overall nucleotide diversity of the PvAMA-1 gene from the 152 YP isolates was 0.007 with 76 haplotypes identified ( Hd = 0.958). Results from the population structure suggested three groups among the YP and Myanmar isolates with optimized clusters value of K = 7. In addition, YP (2014-2016) isolates generally lacked some K components that were commonly found in YP (2009-2011) and Myanmar. Meanwhile, PvAMA-1 domain I is found to be the dominant target of positive diversifying selection and most mutation loci were found in this domain. The mutation frequencies of D107N/A, R112K/T, K120R, E145A, E277K, and R438H in PvAMA-1 were more than 70% in the YP isolates. In conclusion, high genetic diversity and positive selection were found in the PvAMA-1 gene from YP isolates, which are significant findings for the design and development of PvAMA-1-based malaria vaccine., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Cui, Shen, Chen, Kassegne, Shi, Xu, Chen, Wu and Wang.)

Drug resistance in Plasmodium vivax may pose a challenge to malaria elimination. Previous studies have found that P. vivax has a decreased sensitivity to antimalarial drugs in some areas of the Greater Mekong Sub-region. This study aims to investigate the ex vivo drug susceptibilities of P . vivax isolates from the China-Myanmar border and genetic variations of resistance-related genes. A total of 46 P. vivax clinical isolates were assessed for ex vivo susceptibility to seven antimalarial drugs using the schizont maturation assay. The medians of IC 50 (half-maximum inhibitory concentrations) for chloroquine, artesunate, and dihydroartemisinin from 46 parasite isolates were 96.48, 1.95, and 1.63 nM, respectively, while the medians of IC 50 values for piperaquine, pyronaridine, mefloquine, and quinine from 39 parasite isolates were 19.60, 15.53, 16.38, and 26.04 nM, respectively. Sequence polymorphisms in pvmdr1 ( P . vivax multidrug resistance-1), pvmrp1 ( P . vivax multidrug resistance protein 1), pvdhfr ( P .  vivax dihydrofolate reductase), and pvdhps ( P . vivax dihydropteroate synthase) were determined by PCR and sequencing. Pvmdr1 had 13 non-synonymous substitutions, of which, T908S and T958M were fixed, G698S (97.8%) and F1076L (93.5%) were highly prevalent, and other substitutions had relatively low prevalences. Pvmrp1 had three non-synonymous substitutions, with Y1393D being fixed, G1419A approaching fixation (97.8%), and V1478I being rare (2.2%). Several pvdhfr and pvdhps mutations were relatively frequent in the studied parasite population. The pvmdr1 G698S substitution was associated with a reduced sensitivity to chloroquine, artesunate, and dihydroartemisinin. This study suggests the possible emergence of P . vivax isolates resistant to certain antimalarial drugs at the China-Myanmar border, which demands continuous surveillance for drug resistance., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Zeng, Zhao, Zhao, Yang, Li, Li, Duan, Wang, Li, Xiang, Chen, Cui and Yang.)

Duffy binding-like domain (DBL) and cysteine-rich interdomain region (CIDR) domain genes of Plasmodium falciparum erythrocyte membrane protein 1 (PfEMP1) encode malaria virulence proteins. The variants of these genes have been reported to be associated with severe/complicated malaria. The present study investigated the prevalence and distribution patterns of DBLα0.6/9, DBLα1.1, DBLα1 not var3 genes, DBLα2/α1.1/2/4/7, DBLβ12 & DBLβ3/5, DBLε8, CIDRα1.4, and CIDRα1.6 of P. falciparum isolates along the Thai-Myanmar border. The association between PfEMP1 variants and parasite density was also investigated. Two hundred and thirteen finger-prick dried blood spot (DBS) or whole blood samples were collected in 2007 and 2015, from patients with acute uncomplicated P. falciparum in Tak, Kanchanaburi, and Ranong provinces. Analysis of the variant genes was performed using polymerase chain reaction (PCR). The DBLs variant which was found at the highest and lowest frequencies in the three provinces were DBLα1 not var3 (72.77%), and DBLε8 (17.37%). The two CIDR domain variants were found at relatively lower frequencies compared with DBL domain variants (9.9% and 30.1%). P. falciparum isolates carrying the four PfEMP1 variants, i.e., DBLα0.6/9, DBLα1.1, DBLα2/α.1.1/2/4/7, and DBLε8 were found to be significantly associated with low parasitemia. Both DBLα0.6/9 and DBLα2/α1.1/2/4/7 variant genes which were present at high frequencies in this border area could be potential candidate markers for predicting P. falciparum hyperparasitemia and in this border area. Furthermore, the information could be exploited as candidate proteins for the development of an effective malaria vaccine in specific malaria-endemic areas., (Copyright © 2021. Published by Elsevier B.V.)

Antimalarial drug resistance has emerged as a major threat to global malaria control efforts, particularly in the Greater Mekong Subregion (GMS). In this study, we analyzed the polymorphism and prevalence of molecular markers associated with resistance to first-line antimalarial drugs, such as artemisinin, chloroquine, and pyrimethamine, using blood samples collected from malaria patients in the China-Myanmar border region of the GMS from 2008 to 2017, including 225 cases of Plasmodium falciparum and 194 cases of Plasmodium vivax . In artemisinin resistance, only the C580Y mutation with low frequency was detected in pfk13 , and no highly frequent stable mutation was found in pvk12 . In chloroquine resistance, the frequency of K76T mutation in pfcrt was always high, and the frequency of double mutations in pvmdr1 of P. vivax has been steadily increasing every year. In pyrimidine resistance, pfdhfr and pvdhfr had relatively more complex mutant types associated with drug resistance sites, and the overall mutation rate was still high. Therefore, artemisinin-based combination therapies are still suitable for use as the first choice of antimalarial strategy in the China-Myanmar border region in the future., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Tang, Xu, Cao, Tian, Shao, Deng, Zhou and Xiao.)

Drug-resistant  Plasmodium vivax  malaria impedes efforts to control, eliminate, and ultimately eradicate malaria in Southeast Asia. P. vivax resistance to antifolate drugs derives from point mutations in specific parasite genes, including the dihydropteroate synthase ( pvdhps ), dihydrofolate reductase ( pvdhfr ), and GTP cyclohydrolase I ( pvgch1 ) genes. This study aims to investigate the prevalence and spread of drug resistance markers in P. vivax populating the China-Myanmar border. Blood samples were collected from symptomatic patients with acute P. vivax infection. Samples with single-clone P. vivax infections were sequenced for pvdhps and  pvdhfr genes   and genotyped for 6 flanking microsatellite markers. Copy number variation in the pvgch1 gene was also examined. Polymorphisms were observed in six different codons of the  pvdhps  gene (382, 383, 512, 549, 553, and 571) and six different codons of the  pvdhfr  gene (13, 57, 58, 61, 99, 117) in two study sites. The quadruple mutant haplotypes 57I/L/58R/61M/117T of pvdhfr gene were the most common (comprising 76% of cases in Myitsone and 43.7% of case in Laiza). The double mutant haplotype 383G/553G of pvdhps  gene was also prevalent at each site (40.8% and 31%). Microsatellites flanking the pvdhfr gene differentiated clinical samples from wild type and quadruple mutant genotypes ( F ST = 0.259-0.3036), as would be expected for a locus undergoing positive selection. The lack of copy number variation of  pvgch1  suggests that SP-resistant  P. vivax  may harbor alternative mechanisms to secure sufficient folate., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Zeng, Wang, Feng, Zhong, Hu, Bai, Ruan, Si, Zhao, Yang, Li, Chen, Zhang, Li, Xiang, Wu, Chen, Su, Rosenthal and Yang.)

The Plasmodium vivax variant proteins encoded by vir genes are highly polymorphic antigens and are considered as one of key proteins of P. vivax for host immune evasion via antigenic variations. Because genetic diversity of these antigens is a critical hurdle in the development of an effective vaccine, understanding the genetic nature of the vir genes in natural population is important. In this study, we selected four vir genes (vir 4, vir 12, vir 21, and vir 27) previously used for genetic analysis in several studies and evaluated the genetic polymorphisms and phylogenetic relationship of these 4 vir genes in Myanmar P. vivax population. Taken all genetic diversity values, the vir 12 (S = 168, H = 17, Hd = 0.854, Tajima's D value = 2.91524) was the most genetically diverse gene and the vir 4 (S = 9, H = 4, Hd = 0.744, Tajima's D value = -0.49151) was the most conserved gene. All phylogenetic trees showed two clades, and vir 4 and 12 haplotypes from Myanmar were clustered in a distinct clade with those from India and Republic of Korea. These results confirmed the pattern of high genetic polymorphism of vir genes and provided information on vir gene for further functional research and studies focused toward the practical use of vir genes., (Copyright © 2020 Elsevier B.V. All rights reserved.)

Background: Malaria rapid diagnostic tests (RDTs) are precious tools to diagnose malaria. Most RDTs used currently are based on the detection of Plasmodium falciparum histidine-rich protein 2 (PfHRP2) in a patient's blood. However, concern has been raised in recent years that deletion of pfhrp2 in the parasite could affect the accuracy of PfHRP2-based RDTs. In addition, genetic variation in pfhrp2 might influence the accuracy and sensitivity of RDTs. In this study, the genetic variation in pfhrp2 and pfhrp3 in Myanmar P. falciparum isolates was analysed., Methods: Blood samples were collected from malaria patients who were infected with P. falciparum in Mandalay, Naung Cho, Tha Beik Kyin, and Pyin Oo Lwin, Upper Myanmar between 2013 and 2015. The pfhrp2 and pfhrp3 were amplified by nested polymerase chain reaction (PCR), cloned and sequenced. Genetic variation in Myanmar pfhrp2 and pfhrp3 was analysed using the DNASTAR program. Comparative analysis of Myanmar and global pfhrp2 and pfhrp3 isolates was also performed., Results: One-hundred and two pfhrp2 and 89 pfhrp3 were amplified from 105 blood samples, of which 84 pfhrp2 and 56 pfhrp3 sequences were obtained successfully. Myanmar pfhrp2 and pfhrp3 showed high levels of genetic variation with different arrangements of distinct repeat types, which further classified Myanmar pfhrp2 and pfhrp3 into 76 and 47 haplotypes, respectively. Novel amino acid changes were also found in Myanmar pfhrp2 and pfhrp3, but their frequencies were very low. Similar structural organization was shared by Myanmar and global pfhrp2 and pfhrp3, and differences in frequencies of repeat types and lengths were also observed between and among global isolates., Conclusion: Length polymorphisms and amino acid substitutions generated extensive genetic variation in Myanmar pfhrp2 and pfhrp3. Comparative analysis revealed that global pfhrp2 and pfhrp3 share similar structural features, as well as extensive length polymorphisms and distinct organizations of repeat types. These results provide a better understanding of the genetic structure of pfhrp2 and pfhrp3 in global P. falciparum populations and suggest useful information to develop RDTs with improved quality.

Deletion of the pfhrp2 gene in Plasmodium falciparum can lead to false-negative rapid diagnostic test (RDT) results, constituting a major challenge for evidence-based malaria treatment. Here we analyzed the whole genome sequences of 138 P. falciparum clinical samples collected from the China-Myanmar boarder for pfhrp2 and pfhrp3 gene deletions. We found pfhrp2 and pfhrp3 deletions in 9.4% and 3.6% of samples, respectively, with no samples harboring deletions of both genes. The pfhrp2 deletions showed 2 distinct breakpoints, representing 2 different chromosomal deletion events. A phylogenetic analysis performed using genome-wide single-nucleotide polymorphisms revealed that the 2 pfhrp2 breakpoint groups as well as all the pfhrp3-negative parasites formed separate clades, suggesting they might have resulted from clonal expansion of pfhrp2- and pfhrp3-negative parasites. These findings highlight the need for urgent surveys to determine the prevalence of pfhrp2-negative parasites causing false-negative RDT results and a plan for switching of RDTs pending the survey results., (© The Author(s) 2020. Published by Oxford University Press for the Infectious Diseases Society of America. All rights reserved. For permissions, e-mail: journals.permissions@oup.com.)

Background: The emergence and spread of artemisinin resistance in Plasmodium falciparum poses a threat to malaria eradication, including China's plan to eliminate malaria by 2020. Piperaquine (PPQ) resistance has emerged in Cambodia, compromising an important partner drug that is widely used in China in the form of dihydroartemisinin (DHA)-PPQ. Several mutations in a P. falciparum gene encoding a kelch protein on chromosome 13 (k13) are associated with artemisinin resistance and have arisen spread in the Great Mekong subregion, including the China-Myanmar border. Multiple copies of the plasmepsin II/III (pm2/3) genes, located on chromosome 14, have been shown to be associated with PPQ resistance., Methods: The therapeutic efficacy of DHA-PPQ for the treatment of uncomplicated P. falciparum was evaluated along the China-Myanmar border from 2010 to 2014. The dry blood spots samples collected in the efficacy study prior DHA-PPQ treatment and from the local hospital by passive detection were used to amplify k13 and pm2. Polymorphisms within k13 were genotyped by capillary sequencing and pm2 copy number was quantified by relative-quantitative real-time polymerase chain reaction. Treatment outcome was evaluated with the World Health Organization protocol. A linear regression model was used to estimate the association between the day 3 positive rate and k13 mutation and the relationship of the pm2 copy number variants and k13 mutations., Results: DHA-PPQ was effective for uncomplicated P. falciparum infection in Yunnan Province with cure rates > 95%. Twelve non synonymous mutations in the k13 domain were observed among the 268 samples with the prevalence of 44.0% and the predominant mutation was F446I with a prevalence of 32.8%. Only one sample was observed with multi-copies of pm2, including parasites with and without k13 mutations. The therapeutic efficacy of DHA-PPQ was > 95% along the China-Myanmar border, consistent with the lack of amplification of pm2., Conclusion: DHA-PPQ for uncomplicated P. falciparum infection still showed efficacy in an area with artemisinin-resistant malaria along the China-Myanmar border. There was no evidence to show PPQ resistance by clinical study and molecular markers survey. Continued monitoring of the parasite population using molecular markers will be important to track emergence and spread of resistance in this region.

Background: Circumsporozoite surface protein (CSP) of malaria parasites has been recognized as one of the leading vaccine candidates. Clinical trials of vaccines for vivax malaria incorporating Plasmodium vivax CSP (PvCSP) have demonstrated their effectiveness in preventing malaria, at least in part. However, genetic diversity of pvcsp in the natural population remains a major concern., Methods: A total of 171 blood samples collected from patients infected with Plasmodium vivax in Myanmar were analysed in this study. The pvcsp was amplified by polymerase chain reaction, followed by cloning and sequencing. Polymorphic characteristics and natural selection of pvcsp population in Myanmar were analysed using DNASTAR, MEGA6 and DnaSP programs. The polymorphic pattern and natural selection of publicly accessible global pvcsp sequences were also comparatively analysed., Results: Myanmar pvcsp sequences were divided into two subtypes VK210 and VK247 comprising 143 and 28 sequences, respectively. The VK210 subtypes showed higher levels of genetic diversity and polymorphism than the VK247 subtypes. The N-terminal non-repeat region of pvcsp displayed limited genetic variations in the global population. Different patterns of octapeptide insertion (ANKKAEDA in VK210 and ANKKAGDA in VK247) and tetrapeptide repeat motif (GGNA) were identified in the C-terminal region of global pvcsp population. Meanwhile, the central repeat region (CRR) of Myanmar and global pvcsp, both in VK210 and VK247 variants, was highly polymorphic. The high level of genetic diversity in the CRR has been attributed to the different numbers, types and combinations of peptide repeat motifs (PRMs). Interestingly, 27 and 5 novel PRMs were found in Myanmar VK210 and VK247 variants, respectively., Conclusion: Comparative analysis of the global pvcsp population suggests a complex genetic profile of pvcsp in the global population. These results widen understanding of the genetic make-up of pvcsp in the global P. vivax population and provide valuable information for the development of a vaccine based on PvCSP.

Background: Plasmodium malariae is a widely spread but neglected human malaria parasite, which causes chronic infections. Studies on genetic polymorphisms of anti-malarial drug target genes in P. malariae are limited. Previous reports have shown polymorphisms in the P. malariae dihydrofolate reductase gene associated with pyrimethamine resistance and linked to pyrimethamine drug pressure. This study investigated polymorphisms of the P. malariae homologous genes, chloroquine resistant transporter and multidrug resistant 1, associated with chloroquine and mefloquine resistance in Plasmodium falciparum., Methods: The orthologous P. malariae crt and mdr1 genes were studied in 95 patients with P. malariae infection between 2002 and 2016 from Thailand (N = 51) and Myanmar (N = 44). Gene sequences were analysed using BioEdit, MEGA7, and DnaSP programs. Mutations and gene amplifications were compared with P. falciparum and Plasmodium vivax orthologous genes. Protein topology models derived from the observed pmcrt and pmmdr1 haplotypes were constructed and analysed using Phyre2, SWISS MODEL and Discovery Studio Visualization V 17.2., Results: Two non-synonymous mutations were observed in exon 2 (H53P, 40%) and exon 8 (E278D, 44%) of pmcrt. The topology model indicated that H53P and E278D were located outside of the transmembrane domain and were unlikely to affect protein function. Pmmdr1 was more diverse than pmcrt, with 10 non-synonymous and 3 synonymous mutations observed. Non-synonymous mutations were located in the parasite cytoplasmic site, transmembrane 11 and nucleotide binding domains 1 and 2. Polymorphisms conferring amino acid changes in the transmembrane and nucleotide binding domains were predicted to have some effect on PmMDR1 conformation, but were unlikely to affect protein function. All P. malariae parasites in this study contained a single copy of the mdr1 gene., Conclusions: The observed polymorphisms in pmcrt and pmmdr1 genes are unlikely to affect protein function and unlikely related to chloroquine drug pressure. Similarly, the absence of pmmdr1 copy number variation suggests limited mefloquine drug pressure on the P. malariae parasite population, despite its long time use in Thailand for the treatment of falciparum malaria.

Background: In the Greater Mekong sub-region, Plasmodium vivax has become the predominant species and imposes a major challenge for regional malaria elimination. This study aimed to investigate the variations in genes potentially related to drug resistance in P. vivax populations from the China-Myanmar border area. In addition, this study also wanted to determine whether divergence existed between parasite populations associated with asymptomatic and acute infections., Methods: A total of 66 P. vivax isolates were obtained from patients with acute malaria who attended clinics at the Laiza area, Kachin State, Myanmar in 2015. In addition, 102 P. vivax isolates associated with asymptomatic infections were identified by screening of volunteers without signs or symptoms from surrounding villages. Slide-positive samples were verified with nested PCR detecting the 18S rRNA gene. Multiclonal infections were further excluded by genotyping at msp-3α and msp-3β genes. Parasite DNA from 60 symptomatic cases and 81 asymptomatic infections was used to amplify and sequence genes potentially associated with drug resistance, including pvmdr1, pvcrt-o, pvdhfr, pvdhps, and pvk12., Results: The pvmdr1 Y976F and F1076L mutations were present in 3/113 (2.7%) and 97/113 (85.5%) P. vivax isolates, respectively. The K10 insertion in pvcrt-o gene was found in 28.2% of the parasites. Four mutations in the two antifolate resistance genes reached relatively high levels of prevalence: pvdhfr S58R (53.4%), S117N/T (50.8%), pvdhps A383G (75.0%), and A553G (36.3%). Haplotypes with wild-type pvmdr1 (976Y/997K/1076F) and quadruple mutations in pvdhfr (13I/57L/58R/61M/99H/117T/173I) were significantly more prevalent in symptomatic than asymptomatic infections, whereas the pvmdr1 mutant haplotype 976Y/997K/1076L was significantly more prevalent in asymptomatic than symptomatic infections. In addition, quadruple mutations at codons 57, 58, 61 and 117 of pvdhfr and double mutations at codons 383 and 553 of pvdhps were found both in asymptomatic and symptomatic infections with similar frequencies. No mutations were found in the pvk12 gene., Conclusions: Mutations in pvdhfr and pvdhps were prevalent in both symptomatic and asymptomatic P. vivax infections, suggestive of resistance to antifolate drugs. Asymptomatic carriers may act as a silent reservoir sustaining drug-resistant parasite transmission necessitating a rational strategy for malaria elimination in this region.

Background: Vivax malaria is an important public health problem in the Greater Mekong Subregion (GMS), including the China-Myanmar border. Previous studies have found that Plasmodium vivax has decreased sensitivity to antimalarial drugs in some areas of the GMS, but the sensitivity of P. vivax to antimalarial drugs is unclear in the China-Myanmar border. Here, we investigate the drug sensitivity profile and genetic variations for two drug resistance related genes in P. vivax isolates to provide baseline information for future drug studies in the China-Myanmar border., Methodology/principal Findings: A total of 64 P. vivax clinical isolates collected from the China-Myanmar border area were assessed for ex vivo susceptibility to eight antimalarial drugs by the schizont maturation assay. The medians of IC50 (half-maximum inhibitory concentrations) for chloroquine, mefloquine, pyronaridine, piperaquine, quinine, artesunate, artemether, dihydroartemisinin were 84.2 nM, 34.9 nM, 4.0 nM, 22.3 nM, 41.4 nM, 2.8 nM, 2.1 nM and 2.0 nM, respectively. Twelve P. vivax clinical isolates were found over the cut-off IC50 value (220 nM) for chloroquine resistance. In addition, sequence polymorphisms in pvmdr1 (P. vivax multidrug resistance-1), pvcrt-o (P. vivax chloroquine resistance transporter-o), and difference in pvmdr1 copy number were studied. Sequencing of the pvmdr1 gene in 52 samples identified 12 amino acid substitutions, among which two (G698S and T958M) were fixed, M908L were present in 98.1% of the isolates, while Y976F and F1076L were present in 3.8% and 78.8% of the isolates, respectively. Amplification of the pvmdr1 gene was only detected in 4.8% of the samples. Sequencing of the pvcrt-o in 59 parasite isolates identified a single lysine insertion at position 10 in 32.2% of the isolates. The pvmdr1 M908L substitutions in pvmdr1 in our samples was associated with reduced sensitivity to chloroquine, mefloquine, pyronaridine, piperaquine, quinine, artesunate and dihydroartemisinin., Conclusions: Our findings depict a drug sensitivity profile and genetic variations of the P. vivax isolates from the China-Myanmar border area, and suggest possible emergence of chloroquine resistant P. vivax isolates in the region, which demands further efforts for resistance monitoring and mechanism studies., Competing Interests: The authors have declared that no competing interests exist.

Background: Plasmodium falciparum merozoite surface protein-3 (PfMSP-3) is a target of naturally acquired immunity against P. falciparum infection and is a promising vaccine candidate because of its critical role in the erythrocyte invasion of the parasite. Understanding the genetic diversity of pfmsp-3 is important for recognizing genetic nature and evolutionary aspect of the gene in the natural P. falciparum population and for designing an effective vaccine based on the antigen., Methods: Blood samples collected from P. falciparum-infected patients in Naung Cho and Pyin Oo Lwin, Myanmar, in 2015 were used in this study. The pfmsp-3 was amplified by polymerase chain reaction, cloned, and sequenced. Genetic polymorphism and natural selection of Myanmar pfmsp-3 were analysed using the programs DNASTAR, MEGA6, and DnaSP 5.10.00. Genetic diversity and natural selection of the global pfmsp-3 were also comparatively analysed., Results: Myanmar pfmsp-3 displayed 2 different alleles, 3D7 and K1. The 3D7 allelic type was predominant in the population, but genetic polymorphism was less diverse than for the K1 allelic type. Polymorphic characters in both allelic types were caused by amino acid substitutions, insertions, and deletions. Amino acid substitutions were mainly occurred at the alanine heptad repeat domains, whereas most insertions and deletions were found at the glutamate rich domain. Overall patterns of amino acid polymorphisms detected in Myanmar pfmsp-3 were similar in the global pfmsp-3 population, but novel amino acid changes were observed in Myanmar pfmsp-3 with low frequencies. Complicated patterns of natural selection and recombination events were predicted in the global pfmsp-3, which may act as major driving forces to maintain and generate genetic diversity of the global pfmsp-3 population., Conclusion: Global pfmsp-3 revealed genetic polymorphisms, suggesting that the functional and structural consequences of the polymorphisms should be considered in designing a vaccine based on PfMSP-3. Further examination of genetic diversity of pfmsp-3 in the global P. falciparum population is necessary to gain in-depth insight for the population structure and evolutionary aspect of global pfmsp-3.

Background: Antigens expressed in sexual stages of the malaria parasites are targets of transmission-blocking vaccines (TBVs). HAP2/GCS1, a TBV candidate, is critical for fertilization in Plasmodium. Here, the genetic diversity of PvHAP2 was studied in Plasmodium vivax parasite populations from the Greater Mekong Subregion (GMS)., Methods: Plasmodium vivax clinical isolates were collected in clinics from the China-Myanmar border region (135 samples), western Thailand (41 samples) and western Myanmar (51 samples). Near full-length Pvhap2 (nucleotides 13-2574) was amplified and sequenced from these isolates. Molecular evolution studies were conducted to evaluate the genetic diversity, selection and population differentiation., Results: Sequencing of the pvhap2 gene for a total of 227 samples from the three P. vivax populations revealed limited genetic diversity of this gene in the GMS (π = 0.00036 ± 0.00003), with the highest π value observed in Myanmar (0.00053 ± 0.00009). Y133S was the dominant mutation in the China-Myanmar border (99.26%), Myanmar (100%) and Thailand (95.12%). Results of all neutrality tests were negative for all the three populations, suggesting the possible action of purifying selection. Codon-based tests identified specific codons which are under purifying or positive selections. Wright's fixation index showed low to moderate genetic differentiation of P. vivax populations in the GMS, with F ST ranging from 0.04077 to 0.24833, whereas high levels of genetic differentiation were detected between the China-Myanmar border and Iran populations (F ST = 0.60266), and between Thailand and Iran populations (F ST = 0.44161). A total of 20 haplotypes were identified, with H2 being the abundant haplotype in China-Myanmar border, Myanmar and Thailand populations. Epitope mapping prediction of Pvhap2 antigen showed that high-score B-cell epitopes are located in the S307-G324, L429-P453 and V623-D637 regions. The E317K and D637N mutations located within S307-G324 and V623-D637 epitopes slightly reduced the predicted score for potential epitopes., Conclusions: The present study showed a very low level of genetic diversity of pvhap2 gene among P. vivax populations in the Greater Mekong Subregion. The relative conservation of pvhap2 supports further evaluation of a Pvhap2-based TBV.

Background: The malaria elimination plan of the Greater Mekong Subregion (GMS) is jeopardized by the increasing number of Plasmodium vivax infections and emergence of parasite strains with reduced susceptibility to the frontline drug treatment chloroquine/primaquine. This study aimed to determine the evolution of the P. vivax multidrug resistance 1 (Pvmdr1) gene in P. vivax parasites isolated from the China-Myanmar border area during the major phase of elimination., Methods: Clinical isolates were collected from 275 P. vivax patients in 2008, 2012-2013 and 2015 in the China-Myanmar border area and from 55 patients in central China. Comparison was made with parasites from three border regions of Thailand., Results: Overall, genetic diversity of the Pvmdr1 was relatively high in all border regions, and over the seven years in the China-Myanmar border, though slight temporal fluctuation was observed. Single nucleotide polymorphisms previously implicated in reduced chloroquine sensitivity were detected. In particular, M908L approached fixation in the China-Myanmar border area. The Y976F mutation sharply decreased from 18.5% in 2008 to 1.5% in 2012-2013 and disappeared in 2015, whereas F1076L steadily increased from 33.3% in 2008 to 77.8% in 2015. While neutrality tests suggested the action of purifying selection on the pvmdr1 gene, several likelihood-based algorithms detected positive as well as purifying selections operating on specific amino acids including M908L, T958M and F1076L. Fixation and selection of the nonsynonymous mutations are differently distributed across the three border regions and central China. Comparison with the global P. vivax populations clearly indicated clustering of haplotypes according to geographic locations. It is noteworthy that the temperate-zone parasites from central China were completely separated from the parasites from other parts of the GMS., Conclusions: This study showed that P. vivax populations in the China-Myanmar border has experienced major changes in the Pvmdr1 residues proposed to be associated with chloroquine resistance, suggesting that drug selection may play an important role in the evolution of this gene in the parasite populations.

Background: Resistance to anti-malarial drugs hinders malaria elimination. Monitoring the molecular markers of drug resistance helps improve malaria treatment policies. This study aimed to assess the distribution of molecular markers of imported Plasmodium falciparum infections., Methods: In total, 485 P. falciparum cases imported from Africa, Southeast Asia, and Oceania into Zhejiang province, China, from 2016 to 2018 were investigated. Most were imported from Africa, and only a few cases originated in Asia and Oceania. Blood samples were collected from each patient. Plasmodium falciparum chloroquine resistance transporter (Pfcrt) at residues 72-76 and Kelch13-propeller (k13) were determined by nested PCR and DNA sequence., Results: Wild-type Pfcrt at residues 72-76 was predominant (72.61%), but mutant and mixed alleles were also detected, of which CVIET (22.72%) was the most common. Mutant Pfcrt haplotypes were more frequent in patients from West Africa (26.92%), North Africa (25%), and Central Africa (21.93%). The number of cases of P. falciparum infections was small in Southeast Asia and Oceania, and these cases involved Pfcrt mutant type. For the k13 propeller gene, 26 samples presented 19 different point mutations, including eight nonsynonymous mutations (P441S, D464E, K503E, R561H, A578S, R622I, V650F, N694K). In addition, R561H, one of the validated SNPs in k13, was detected in one patient from Myanmar and one patient from Rwanda. A578S, although common in Africa, was found in only one patient from Cameroon. R622I was detected in one sample from Mozambique and one sample from Somalia. The genetic diversity of k13 was low in most regions of Africa and purifying selection was suggested by Tajima's D test., Conclusions: The frequency and spatial distributions of Pfcrt and k13 mutations associated with drug resistance were determined. Wild-type Pfcrt was dominant in Africa. Among k13 mutations correlated with delayed parasite clearance, only the R561H mutation was found in one case from Rwanda in Africa. Both Pfcrt and k13 mutations were detected in patients from Southeast Asia and Oceania. These findings provide insights into the molecular epidemiological profile of drug resistance markers in the study region.

Background: Plasmodium lactate dehydrogenase (pLDH) is a major target in diagnosing the erythrocytic stage of malaria parasites because it is highly expressed during blood-stage parasites and is distinguished from human LDH. Rapid diagnostic tests (RDTs) for malaria use pLDH as a target antigen; however, genetic variations in pLDH within the natural population threaten the efficacy of pLDH-based RDTs., Methods: Genetic polymorphisms of Plasmodium vivax LDH (PvLDH) and Plasmodium falciparum LDH (PfLDH) in Myanmar isolates were analysed by nucleotide sequencing analysis. Genetic polymorphisms and the natural selection of PvLDH and PfLDH were analysed using DNASTAR, MEGA6, and DnaSP ver. 5.10.00 programs. The genetic diversity and natural selection of global PvLDH and PfLDH were also analysed. The haplotype network of global PvLDH and PfLDH was constructed using NETWORK ver. 5.0.0.3. Three-dimensional structures of PvLDH and PfLDH were built with YASARA Structure ver. 18.4.24 and the impact of mutations on structural change and stability was evaluated with SDM ver. 2, CUPSAT and MAESTROweb., Results: Forty-nine PvLDH and 52 PfLDH sequences were obtained from Myanmar P. vivax and P. falciparum isolates. Non-synonymous nucleotide substitutions resulting in amino acid changes were identified in both Myanmar PvLDH and PfLDH. Amino acid changes were also identified in the global PvLDH and PfLDH populations, but they did not produce structural alterations in either protein. Low genetic diversity was observed in global PvLDH and PfLDH, which may be maintained by a strong purifying selection., Conclusion: This study extends knowledge for genetic diversity and natural selection of global PvLDH and PfLDH. Although amino acid changes were observed in global PvLDH and PfLDH, they did not alter the conformational structures of the proteins. These suggest that PvLDH and PfLDH are genetically well-conserved in global populations, which indicates that they are suitable antigens for diagnostic purpose and attractive targets for drug development.

Background: Merozoite proteins of the malaria parasites involved in the invasion of red blood cells are selected by host immunity and their diversity is greatly influenced by changes in malaria epidemiology. In the Greater Mekong Subregion (GMS), malaria transmission is concentrated along the international borders and there have been major changes in malaria epidemiology with Plasmodium vivax becoming the dominant species in many regions. Here, we aimed to evaluate the genetic diversity of P. vivax Duffy-binding protein gene domain II (pvdbp-II) in isolates from the eastern and western borders of Myanmar, and compared it with that from global P. vivax populations., Methods: pvdbp-II sequences were obtained from 85 and 82 clinical P. vivax isolates from the eastern and western Myanmar borders, respectively. In addition, 504 pvdbp-II sequences from nine P. vivax populations of the world were retrieved from GenBank and used for comparative analysis of genetic diversity, recombination and population structure of the parasite population., Results: The nucleotide diversity of the pvdbp-II sequences from the Myanmar border parasite isolates was not uniform, with the highest diversity located between nucleotides 1078 and 1332. Western Myanmar isolates had a unique R391C mutation. Evidence of positive natural selection was detected in pvdbp-II gene in P. vivax isolates from the eastern Myanmar area. P. vivax parasite populations in the GMS, including those from the eastern, western, and central Myanmar as well as Thailand showed low-level genetic differentiation (F ST , 0.000-0.099). Population genetic structure analysis of the pvdbp-II sequences showed a division of the GMS populations into four genetic clusters. A total of 60 PvDBP-II haplotypes were identified in 210 sequences from the GMS populations. Among the epitopes in PvDBP-II, high genetic diversity was found in epitopes 45 (379-SIFGT(D/G)(E/K)(K/N)AQQ(R/H)(R/C)KQ-393, π = 0.029) and Ia (416-G(N/K)F(I/M)WICK(L/I)-424], Ib [482-KSYD(Q/E)WITR-490, π = 0.028) in P. vivax populations from the eastern and western borders of Myanmar., Conclusions: The pvdbp-II gene is genetically diverse in the eastern and western Myanmar border P. vivax populations. Positive natural selection and recombination occurred in pvdbp-II gene. Low-level genetic differentiation was identified, suggesting extensive gene flow of the P. vivax populations in the GMS. These results can help understand the evolution of the P. vivax populations in the course of regional malaria elimination and guide the design of PvDBP-II-based vaccine.

Background: Malaria is a major public health problem in the China-Myanmar border region. The genetic structure of malaria parasite may affect its transmission model and control strategies. The present study was to analyse genetic diversity of Plasmodium falciparum by merozoite surface proteins 1 and 2 (MSP1 and MSP2) and to determine the multiplicity of infection in clinical isolates in the China-Myanmar border region., Methods: Venous blood samples (172) and filter paper blood spots (70) of P. falciparum isolates were collected from the patients of the China-Myanmar border region from 2006 to 2011. The genomic DNA was extracted, and the msp1 and msp2 genes were genotyped by nested PCR using allele-specific primers for P. falciparum., Results: A total of 215 P. falciparum clinical isolates were genotyped at the msp1 (201) and msp2 (204), respectively. For the msp1 gene, MAD20 family was dominant (53.49%), followed by the K1 family (44.65%), and the RO33 family (12.56%). For the msp2 gene, the most frequent allele was the FC27 family (80.93%), followed by the 3D7 family (75.81%). The total multiplicity of infection (MOI) of msp1 and msp2 was 1.76 and 2.21, with a prevalence of 64.19% and 72.09%, respectively. A significant positive correlation between the MOI and parasite density was found in the msp1 gene of P. falciparum. Sequence analysis revealed 38 different alleles of msp1 (14 K1, 23 MAD20, and 1 RO33) and 52 different alleles of msp2 (37 3D7 and 15 FC27)., Conclusion: The present study showed the genetic polymorphisms with diverse allele types of msp1 and msp2 as well as the high MOI of P. falciparum clinical isolates in the China-Myanmar border region.

Transmission-blocking vaccines (TBVs) target the sexual stages of malarial parasites to interrupt or reduce the transmission cycle have been one of approaches to control malaria. Pvs25 and Pvs28 are the leading candidate antigens of TBVs against vivax malaria. In this study, genetic diversity and natural selection of the two TBV candidate genes in Plasmodium vivax Myanmar isolates were analyzed. The 62 Myanmar P. vivax isolates showed 9 and 19 different haplotypes for Pvs25 and Pvs28, respectively. The nucleotide diversity of Pvs28 was slightly higher than Pvs25, but not significant. Most amino acid substitutions observed in Myanmar Pvs25 and Pvs28 were concentrated at the EGF-2 and EGF-3 like domains. Major amino acid changes found in Myanmar Pvs25 and Pvs28 were similar to those reported in the global population, but novel amino acid substitutions were also identified. Negative selection was predicted in Myanmar Pvs25, whereas Pvs28 was under positive selection. Comparative analysis of global Pvs25 and Pvs28 suggests a substantial geographical difference between the Asian and American/African Pvs25 and Pvs28. The geographical genetic differentiation and the evidence for natural selection in global Pvs25 and Pvs28 suggest that the functional consequences of the observed polymorphism need to be considered for the development of effective TBVs based on the antigens., (Copyright © 2019 Elsevier B.V. All rights reserved.)

The emergence and spread of drug resistance is a problem hindering malaria elimination in Southeast Asia. In this study, genetic variations in drug resistance markers of Plasmodium falciparum were determined in parasites from asymptomatic populations located in three geographically dispersed townships of Myanmar by PCR and sequencing. Mutations in dihydrofolate reductase ( pfdhfr ), dihydropteroate synthase ( pfdhps ), chloroquine resistance transporter ( pfcrt ), multidrug resistance protein 1 ( pfmdr1 ), multidrug resistance-associated protein 1 ( pfmrp1 ), and Kelch protein 13 ( k13 ) were present in 92.3%, 97.6%, 84.0%, 98.8%, and 68.3% of the parasites, respectively. The pfcrt K76T, pfmdr1 N86Y, pfmdr1 I185K, and pfmrp1 I876V mutations were present in 82.7%, 2.5%, 87.5%, and 59.8% isolates, respectively. The most prevalent haplotypes for pfdhfr, pfdhps, pfcrt and pfmdr1 were 51I/59R/108N/164L, 436A/437G/540E/581A, 74I/75E/76T/220S/271E/326N/356T/371I, and 86N/130E/184Y/185K/1225V, respectively. In addition, 57 isolates had three different point mutations (K191T, F446I, and P574L) and three types of N-terminal insertions (N, NN, NNN) in the k13 gene. In total, 43 distinct haplotypes potentially associated with multidrug resistance were identified. These findings demonstrate a high prevalence of multidrug-resistant P. falciparum in asymptomatic infections from diverse townships in Myanmar, emphasizing the importance of targeting asymptomatic infections to prevent the spread of drug-resistant P. falciparum .

Mutations in the Kelch domain of the K13 gene (PF3D7_1343700) were previously associated with artemisinin resistance in Plasmodium falciparum. This study followed the dynamics of the K13 polymorphisms in P. falciparum parasites from the China-Myanmar border area obtained in 2007-2016, and their in vitro sensitivities to artesunate (AS) and dihydroartemisinin (DHA). The 50% effective concentration (EC 50 72h ) values of 133 culture-adapted field isolates to AS and DHA, measured by the conventional 72 h SYBR Green I-based assay, varied significantly among the parasites from different years; all were significantly higher than that of the reference strain 3D7. Compared with parasites from 2007 to 2008, ring survival rates almost doubled in parasites obtained in later years. Sequencing the full-length K13 genes identified 11 point mutations present in 85 (63.9%) parasite isolates. F446I was the predominant (55/133) variant, and its frequency was increased from 17.6% (3/17) in 2007 to 55.9% (19/34) in 2014-2016. No wild-type (WT) Kelch domain sequences were found in the 34 samples obtained from 2014 to 2016. In the 2014-2016 samples, a new mutation (G533S) appeared and reached 44.1% (15/34). Collectively, parasites with the Kelch domain mutations (after amino acid 440) had significantly higher ring survival rates than the WT parasites. Individually, F446I, G533S and A676D showed significantly higher ring survival rates than the WT. Although the drug sensitivity phenotypes measured by the RSA 6h assays may be intrinsically linked to the in vivo clinical efficacy data, the values determined by these two assays were not significantly correlated. This study identified the trend of K13 mutations in parasite populations from the China-Myanmar border area, confirmed an overall correlation of Kelch domain mutations with elevated ring-stage survival rates, and emphasized the importance of monitoring the evolution and spread of parasites with reduced artemisinin sensitivity along the malaria elimination course.50 72h assays may be intrinsically linked to the in vivo clinical efficacy data, the values determined by these two assays were not significantly correlated. This study identified the trend of K13 mutations in parasite populations from the China-Myanmar border area, confirmed an overall correlation of Kelch domain mutations with elevated ring-stage survival rates, and emphasized the importance of monitoring the evolution and spread of parasites with reduced artemisinin sensitivity along the malaria elimination course., (Copyright © 2019 The Authors. Published by Elsevier Ltd.. All rights reserved.)

Background: Plasmodium falciparum merozoite surface protein-1 (PfMSP-1) and -2 (PfMSP-2) are major blood-stage vaccine candidate antigens. Understanding the genetic diversity of the genes, pfmsp-1 and pfmsp-2, is important for recognizing the genetic structure of P. falciparum, and the development of an effective vaccine based on the antigens. In this study, the genetic diversities of pfmsp-1 and pfmsp-2 in the Myanmar P. falciparum were analysed., Methods: The pfmsp-1 block 2 and pfmsp-2 block 3 regions were amplified by polymerase chain reaction from blood samples collected from Myanmar patients who were infected with P. falciparum in 2013-2015. The amplified gene fragments were cloned into a T&A vector, and sequenced. Sequence analysis of Myanmar pfmsp-1 block 2 and pfmsp-2 block 3 was performed to identify the genetic diversity of the regions. The temporal genetic changes of both pfmsp-1 and pfmsp-2 in the Myanmar P. falciparum population, as well as the polymorphic diversity in the publicly available global pfmsp-1 and pfmsp-2, were also comparatively analysed., Results: High levels of genetic diversity of pfmsp-1 and pfmsp-2 were observed in the Myanmar P. falciparum isolates. Twenty-eight different alleles of pfmsp-1 (8 for K1 type, 14 for MAD20 type, and 6 for RO33 type) and 59 distinct alleles of pfmsp-2 (18 for FC27, and 41 for 3D7 type) were identified in the Myanmar P. falciparum population in amino acid level. Comparative analyses of the genetic diversity of the Myanmar pfmsp-1 and pfmsp-2 alleles in the recent (2013-2015) and past (2004-2006) Myanmar P. falciparum populations indicated the dynamic genetic expansion of the pfmsp-1 and pfmsp-2 in recent years, suggesting that a high level of genetic differentiation and recombination of the two genes may be maintained. Population genetic structure analysis of the global pfmsp-1 and pfmsp-2 also suggested that a high level of genetic diversity of the two genes was found in the global P. falciparum population., Conclusion: Despite the recent remarkable decline of malaria cases, the Myanmar P. falciparum population still remains of sufficient size to allow the generation and maintenance of genetic diversity. The high level of genetic diversity of pfmsp-1 and pfmsp-2 in the global P. falciparum population emphasizes the necessity for continuous monitoring of the genetic diversity of the genes for better understanding of the genetic make-up and evolutionary aspect of the genes in the global P. falciparum population.

Background: Malaria reduction and future elimination in China is made more difficult by the importation of cases from neighboring endemic countries, particularly Myanmar, Laos, and Vietnam, and increased travel to Africa by Chinese nationals. The increasing prevalence of artemisinin resistant parasites across Southeast Asia highlights the importance of monitoring the parasite importation into China. Artemisinin resistance in the Mekong region is associated with variants of genes encoding the K13 kelch domain protein (pf13k), found in specific genetic backgrounds, including certain alleles of genes encoding the chloroquine resistance transporter (pfcrt) and multidrug resistance transporter PgH1 (pfmdr1)., Methods: In this study we investigated the prevalence of drug resistance markers in 72 P. falciparum samples from uncomplicated malaria infections in Tengchong and Yingjiang, counties on the Yunnan-Myanmar border. Variants of pf13k, pfcrt and pfmdr1 are described., Results: Almost all parasites harboured chloroquine-resistant alleles of pfcrt, whereas pfmdr1 was more diverse. Major mutations in the K13 propeller domain associated with artemisinin resistance in the Mekong region (C580Y, R539T and Y493H) were absent, but F446I and two previously undescribed mutations (V603E and V454I) were identified. Protein structural modelling was carried out in silico on each of these K13 variants, based on recently published crystal structures for the K13 propeller domain. Whereas F446I was predicted to elicit a moderate destabilisation of the propeller structure, the V603E substitution is likely to lead to relatively high protein instability. We plotted these stability estimates, and those for all previously described variants, against published values for in vivo parasitaemia half-life, and found that quadratic regression generates a useful predictive algorithm., Conclusion: This study provides a baseline of P. falciparum resistance-associated mutations prevalent at the China-Myanmar border. We also show that protein modelling can be used to generate testable predictions as to the impact of pfk13 mutations on in vivo (and potentially in vitro) artemisinin susceptibility., Competing Interests: The authors have declared that no competing interests exist.

Background: The malaria burden of Myanmar still remains high within the Greater Mekong Subregion of Southeast Asia. An important indicator of progress towards malaria elimination is the prevalence of parasite infections in endemic populations. Information about malaria epidemiology is mostly derived from reports of confirmed acute malaria cases through passive case detection, whereas the prevalence of baseline subclinical malaria infections is much less known., Methods: In this study, cross-sectional surveys were conducted during the rainy season of 2017 in four townships (Bilin, Thabeikkyin, Banmauk and Paletwa) of Myanmar with divergent annual malaria incidences. A total of 1991 volunteers were recruited from local villages and Plasmodium subclinical infections were estimated by light microscopy (LM), rapid diagnostic tests (RDTs) and nested PCR. The nested PCR analysis was performed with a modified pooling strategy that was optimized based on an initial estimate the infection prevalence., Results: The overall malaria infection prevalence based on all methods was 13.9% (277/1991) and it differed drastically among the townships, with Paletwa in the western border having the highest infection rate (22.9%) and Thabeikkyin in central Myanmar having the lowest (3.9%). As expected, nested PCR was the most sensitive and identified 226 (11.4%) individuals with parasite infections. Among the parasite species, Plasmodium vivax was the most prevalent in all locations, while Plasmodium falciparum also accounted for 32% of infections in the western township Paletwa. Two RDTs based on the detection of the hrp2 antigen detected a total of 103 P. falciparum infections, and the ultrasensitive RDT detected 20% more P. falciparum infections than the conventional RDT. In contrast, LM missed the majority of the subclinical infections and only identified 14 Plasmodium infections., Conclusions: Cross-sectional surveys identified considerable levels of asymptomatic Plasmodium infections in endemic populations of Myanmar with P. vivax becoming the predominant parasite species. Geographical heterogeneity of subclinical infections calls for active surveillance of parasite infections in endemic areas. The pooling scheme designed for nested PCR analysis offers a more practical strategy for large-scale epidemiological studies of parasite prevalence. Such information is important for decision-makers to put forward a more realistic action plan for malaria elimination.

Background: The problem of anti-malarial drug resistance is a long-term challenge faced by malaria control in Yunnan Province. Recently, the detection rates of chloroquine-resistant molecular markers (Plasmodium falciparum chloroquine resistant transporter, Pfcrt) and artemisinin-resistant molecular markers (P. falciparum kelch13 gene, ork13) were 85% and 35%, respectively. To understand the association of k13 gene mutation with artemisinin resistance in falciparum malaria cases, the difference in k13 gene differentiation between two populations and artemisinin resistance phenotype on falciparum malaria cases in Myanmar were analysed in this study., Methods: This research involved all of falciparum malaria cases diagnosed continuously in Yunnan Province from 2013 to 2015 and some of falciparum malaria cases found in Lazar, Myanmar. Blood samples were taken from the former group for molecular epidemiological analysis of k13 gene mutations, and artemisinin resistance phenotypes of P. falciparum were observed in the latter group using the in vivo testing method recommended by the World Health Organization. Nested PCR was used to amplify the propeller domain of the k13 gene in P. falciparum, followed by sequencing., Results: A total of 202 blood samples were collected from Yunnan Province and 382 blood samples were collected from falciparum malaria cases in Myanmar. 49 of 382 Myanmar cases were in vivo tested for artesunate resistance phenotype through full treatment course observation. At the same time, all the blood samples were screened for k13 gene mutation of P. falciparum. The genetic diversity of k13 was higher in the Plasmodium isolates from Yunnan Province than those from Myanmar cases. The genetic differentiation index of the two populations was 0.0410, where the intra- and inter-group variations were 95.9% and 4.1%, respectively. The odds ratio of artemisinin resistance phenotype and mutation at the locus 446 in k13 gene in Myanmar cases was 1.640, while the value was 1.840 based on the estimations of the mutations in the 12 loci., Conclusion: Although the Plasmodium isolates from Yunnan Province and those from Myanmar were collected from different sites, they still belong to the same geographical population. It is, therefore, reasonable to contrast the artemisinin resistance status of the Plasmodium population from Myanmar with the Plasmodium population from Yunnan Province. As a result, based on the molecular epidemiological investigation on k13 mutations of Plasmodium isolates in Yunnan Province and the determination of the artemisinin resistance on falciparum malaria cases in Myanmar, the positively genetic correlated was found between the k13 locus mutations with artemisinin resistance phenotype. This provides a basis for further monitoring the artemisinin resistance by detection some molecular markers in k13 gene of Plasmodium in Yunnan Province.

Fitness costs are key determinants of whether drug resistance alleles establish and how fast they spread within populations. More than 125 different kelch13 alleles, each containing a different amino acid substitution, have arisen in Southeast Asian malaria parasite ( Plasmodium falciparum ) populations under artemisinin selection over the past 15 years in a dramatic example of a soft selective event. However, just one of these alleles (C580Y) is now outcompeting other alleles in multiple different countries and is spreading toward fixation. Here we examine the fitness consequences of C580Y, relative to another less successful kelch13 mutation (R561H), to try to explain the distinctive dynamics of C580Y. We hypothesized that C580Y will show lower fitness costs than other kelch13 substitutions in the absence of artemisinin treatment. We used CRISPR/Cas9 methods to introduce single mutations (C580Y or R561H) or synonymous control edits into a wild-type parasite isolated on the Thailand-Myanmar border, conducted replicated head-to-head competition assays, and determined the outcome of competition using deep sequencing of kelch13 amplicons. Contrary to our predictions, these experiments reveal that C580Y carries higher fitness costs ( s [selection coefficient] = 0.15 ± 0.008 [1 standard error {SE}]) than R561H ( s = 0.084 ± 0.005). Furthermore, R561H outcompetes C580Y in direct competition ( s = 0.065 ± 0.004). We conclude that fitness costs of C580Y in isolation are unlikely to explain the rapid spread of this substitution., (Copyright © 2018 Nair et al.)

Multidrug-resistant Plasmodium falciparum in the Greater Mekong Subregion of Southeast Asia is a major threat to malaria elimination and requires close surveillance. In this study, we collected 107 longitudinal clinical samples of P. falciparum in 2007-2012 from the malaria hypoendemic region of the China-Myanmar border and measured their in vitro susceptibilities to 10 antimalarial drugs. Overall, parasites had significantly different IC 50 values to all the drugs tested as compared to the reference 3D7 strain. Parasites were also genotyped in seven genes that were associated with drug resistance including pfcrt, pfmdr1, pfmrp1, pfdhfr, pfdhps, pfnhe1, and PfK13 genes. Despite withdrawal of chloroquine and antifolates from treating P. falciparum, parasites remained highly resistant to these drugs and mutations in pfcrt, pfdhfr, and pfdhps genes were highly prevalent and almost reached fixation in the study parasite population. Except for pyronaridine, quinine and lumefantrine, all other tested drugs exhibited significant temporal variations at least between some years, but only chloroquine and piperaquine had a clear temporal trend of continuous increase of IC 50 s. For the pfmrp1 gene, several mutations were associated with altered sensitivity to a number of drugs tested including chloroquine, piperaquine, lumefantrine and dihydroartemisinin. The association of PfK13 mutations with resistance to multiple drugs suggests potential evolution of PfK13 mutations amid multidrug resistance genetic background. Furthermore, network analysis of drug resistance genes indicated that certain haplotypes associated multidrug resistance persisted in these years, albeit there were year-to-year fluctuations of the predominant haplotypes., (Copyright © 2018. Published by Elsevier Ltd.)

Background: Artemisinin resistance in Plasmodium falciparum has emerged and spread in Southeast Asia. In areas where resistance is established longer courses of artemisinin-based combination therapy have improved cure rates., Methods: The standard 3-day course of artemether-lumefantrine (AL) was compared with an extended 5-day regimen for the treatment of uncomplicated falciparum malaria in Kayin state in South-East Myanmar, an area of emerging artemisinin resistance. Late parasite clearance dynamics were described by microscopy and quantitative ultra-sensitive PCR. Patients were followed up for 42 days., Results: Of 154 patients recruited (105 adults and 49 children < 14 years) 78 were randomized to 3 days and 76 to 5 days AL. Mutations in the P. falciparum kelch13 propeller gene (k13) were found in 46% (70/152) of infections, with F446I the most prevalent propeller mutation (29%; 20/70). Both regimens were well-tolerated. Parasite clearance profiles were biphasic with a slower submicroscopic phase which was similar in k13 wild-type and mutant infections. The cure rates were 100% (70/70) and 97% (68/70) in the 3- and 5-day arms respectively. Genotyping of the two recurrences was unsuccessful., Conclusion: Despite a high prevalence of k13 mutations, the current first-line treatment, AL, was still highly effective in this area of South-East Myanmar. The extended 5 day regimen was very well tolerated, and would be an option to prolong the useful therapeutic life of AL. Trial registration NCT02020330. Registered 24 December 2013, https://clinicaltrials.gov/NCT02020330.

Artemisinin resistance in Plasmodium falciparum was associated with mutations in the propeller domain of the PfK13 gene and increased phosphatidylinositol-3'-kinase (PfPI3K) activity. Assessment of the genetic diversity of the PfK13 ortholog PvK12 in Plasmodium vivax field samples from the same hotspots of P. falciparum artemisinin resistance revealed a limited genetic diversity of PvK12. Following the same logic, we analyzed genetic variations of the PvPI3K gene in 188 P. vivax field isolates from two geographic locations along the China-Myanmar border. Overall, high genetic diversity of PvPI3K was observed; parasites from Yunnan's Tengchong County had higher genetic diversity than those from Laiza Township, Kachin State, Myanmar. Almost all the neutrality tests applied detected statistically significant deviation from zero. The negative Tajima's D values in both populations implicated that PvPI3K gene might have experienced either a directional selection or an expansion in population size. There was low linkage disequilibrium between the PvPI3K mutations in both populations, suggesting the existence of large, almost panmictic, parasite populations that enabled effective recombination. This later result was confirmed by the detection of a minimum of five recombination events in each population with two major breakpoints. Multiple tests for selection confirmed a signature of purifying selection on PvPI3K. All the amino acid mutations were predicted to be neutral for the PI3K protein's function. These findings provide insights on the genetic diversity of P. vivax populations along the China-Myanmar border., (Copyright © 2018 Elsevier B.V. All rights reserved.)

Background: Plasmodium falciparum apical membrane antigen-1 (PfAMA-1) is one of leading blood stage malaria vaccine candidates. However, genetic variation and antigenic diversity identified in global PfAMA-1 are major hurdles in the development of an effective vaccine based on this antigen. In this study, genetic structure and the effect of natural selection of PfAMA-1 among Myanmar P. falciparum isolates were analysed., Methods: Blood samples were collected from 58 Myanmar patients with falciparum malaria. Full-length PfAMA-1 gene was amplified by polymerase chain reaction and cloned into a TA cloning vector. PfAMA-1 sequence of each isolate was sequenced. Polymorphic characteristics and effect of natural selection were analysed with using DNASTAR, MEGA4, and DnaSP programs. Polymorphic nature and natural selection in 459 global PfAMA-1 were also analysed., Results: Thirty-seven different haplotypes of PfAMA-1 were identified in 58 Myanmar P. falciparum isolates. Most amino acid changes identified in Myanmar PfAMA-1 were found in domains I and III. Overall patterns of amino acid changes in Myanmar PfAMA-1 were similar to those in global PfAMA-1. However, frequencies of amino acid changes differed by country. Novel amino acid changes in Myanmar PfAMA-1 were also identified. Evidences for natural selection and recombination event were observed in global PfAMA-1. Among 51 commonly identified amino acid changes in global PfAMA-1 sequences, 43 were found in predicted RBC-binding sites, B-cell epitopes, or IUR regions., Conclusions: Myanmar PfAMA-1 showed similar patterns of nucleotide diversity and amino acid polymorphisms compared to those of global PfAMA-1. Balancing natural selection and intragenic recombination across PfAMA-1 are likely to play major roles in generating genetic diversity in global PfAMA-1. Most common amino acid changes in global PfAMA-1 were located in predicted B-cell epitopes where high levels of nucleotide diversity and balancing natural selection were found. These results highlight the strong selective pressure of host immunity on the PfAMA-1 gene. These results have significant implications in understanding the nature of Myanmar PfAMA-1 along with global PfAMA-1. They also provide useful information for the development of effective malaria vaccine based on this antigen.

Background: Artemisinin resistance, linked to polymorphisms in the Kelch gene on chromosome 13 of Plasmodium falciparum (k13), has outpaced containment efforts in South East Asia. For national malaria control programmes in the region, it is important to establish a surveillance system which includes monitoring for k13 polymorphisms associated with the clinical phenotype., Methods: Between February and December 2013, parasite clearance was assessed in 35 patients with uncomplicated P. falciparum treated with artesunate monotherapy followed by 3-day ACT in an isolated area on the Myanmar-Thai border with relatively low artemisinin drug pressure. Molecular testing for k13 mutations was performed on dry blood spots collected on admission., Results: The proportion of k13 mutations in these patients was 41.7%, and only 5 alleles were detected: C580Y, I205T, M476I, R561H, and F446I. Of these, F446I was the most common, and was associated with a longer parasite clearance half-life (median) 4.1 (min-max 2.3-6.7) hours compared to 2.5 (min-max 1.6-8.7) in wildtype (p = 0·01). The prevalence of k13 mutant parasites was much lower than the proportion of k13 mutants detected 200 km south in a much less remote setting where the prevalence of k13 mutants was 84% with 15 distinct alleles in 2013 of which C580Y predominated., Conclusions: This study provides evidence of artemisinin resistance in a remote part of eastern Myanmar. The prevalence of k13 mutations as well as allele diversity varies considerably across short distances, presumably because of historical patterns of artemisinin use and population movements.

Background: Emergence of artemisinin-resistant malaria in Southeast Asian countries threatens the global control of malaria. Although K13 kelch propeller has been assessed for artemisinin resistance molecular marker, most of the mutations need to be validated. In this study, artemisinin resistance was assessed by clinical and molecular analysis, including k13 and recently reported markers, pfarps10, pffd and pfmdr2., Methods: A prospective cohort study in 1160 uncomplicated falciparum patients was conducted after treatment with artemisinin-based combination therapy (ACT), in 6 sentinel sites in Myanmar from 2009 to 2013. Therapeutic efficacy of ACT was assessed by longitudinal follow ups. Molecular markers analysis was done on all available day 0 samples., Results: True recrudescence treatment failures cases and day 3 parasite positivity were detected at only the southern Myanmar sites. Day 3 positive and k13 mutants with higher prevalence of underlying genetic foci predisposing to become k13 mutant were detected only in southern Myanmar since 2009 and comparatively fewer mutations of pfarps10, pffd, and pfmdr2 were observed in western Myanmar. K13 mutations, V127M of pfarps10, D193Y of pffd, and T448I of pfmdr2 were significantly associated with day 3 positivity (OR: 6.48, 3.88, 2.88, and 2.52, respectively)., Conclusions: Apart from k13, pfarps10, pffd and pfmdr2 are also useful for molecular surveillance of artemisinin resistance especially where k13 mutation has not been reported. Appropriate action to eliminate the resistant parasites and surveillance on artemisinin resistance should be strengthened in Myanmar. Trial registration This study was registered with ClinicalTrials.gov, identifier NCT02792816.

Background: Recent emergence of artemisinin-resistant P. falciparum has posed a serious hindrance to the elimination of malaria in the Greater Mekong Subregion. Parasite clearance time, a measure of change in peripheral parasitaemia in a sequence of samples taken after treatment, can be used to reflect the susceptibility of parasites or the efficiency of antimalarials. The association of genetic polymorphisms and artemisinin resistance has been documented. This study aims to examine clearance time of P. falciparum and P. vivax parasitemia as well as putative gene mutations associated with residual or recurred parasitemia in Myanmar., Methods: A total of 63 P. falciparum and 130 P. vivax samples collected from two internally-displaced populations and one surrounding village were examined for parasitemia changes. At least four samples were taken from each patient, at the first day of diagnosis up to 3 months following the initial treatment. The amount of parasite gene copy number was estimated using quantitative real-time PCR based on a species-specific region of the 18S rRNA gene. For samples that showed residual or recurred parasitemia after treatment, microsatellites were used to identify the 'post-treatment' parasite genotype and compared such with the 'pre-treatment' genotype. Mutations in genes pfcrt, pfmdr1, pfatp6, pfmrp1 and pfK13 that are potentially associated with ACT resistance were examined to identify if mutation is a factor for residual or persistent parasitemia., Results: Over 30% of the P. falciprium infections showed delayed clearance of parasitemia after 2-3 days of treatment and 9.5% showed recurred parasitemia. Mutations in codon 876 of the pfmrp1 corroborated significance association with slow clearance time. However, no association was observed in the variation in pfmdr1 gene copy number as well as mutations of various codonsinpfatp6, pfcrt, and pfK13 with clearance time. For P. vivax, over 95% of the infections indicated cleared parasitemia at days 2-3 of treatment. Four samples were found to be re-infected with new parasite strains based on microsatellite genotypes after initial treatment., Conclusion: The appearance of P.falciparum infected samples showing delayed clearance or recurred parasitemia after treatment raises concerns on current treatment and ACT drug resistance.

Background: Artemisinin resistance in Plasmodium falciparum extends across Southeast Asia where it is associated with worsening partner drug resistance and a decline in the efficacy of frontline artemisinin-based combination therapy. Dihydroartemisinin-piperaquine (DP) is an essential component of preventive and curative treatment in the region, but its therapeutic efficacy has fallen in Cambodia., Methods: A prospective clinical and parasitological evaluation of DP was conducted at two sites in Upper Myanmar between August 2013 and December 2014, enrolling 116 patients with acute uncomplicated falciparum malaria. Patients received DP orally for 3 days together with primaquine 0.25 mg/kg on admission. Parasite clearance half-lives based on 6 hourly blood smears, and day 42 therapeutic responses were assessed as well as parasite K13 genotypes., Results: Median parasite clearance half-life was prolonged, and clearance half-life was greater than 5 h in 21% of patients. Delayed parasite clearance was significantly associated with mutations in the propeller region of the parasite k13 gene. The k13 F446I mutation was found in 25.4% of infections and was associated with a median clearance half-life of 4.7 h compared with 2.7 h for infections without k13 mutations (p < 0.001). There were no failures after 42 days of follow-up, although 18% of patients had persistent parasitaemia on day 3., Conclusion: The dominant k13 mutation observed in Upper Myanmar, F446I, appears to be associated with an intermediate rate of parasite clearance compared to other common mutations described elsewhere in the Greater Mekong Subregion. Discerning this phenotype requires relatively detailed clearance measurements, highlighting the importance of methodology in assessing artemisinin resistance.

Malaria in the China-Myanmar border region is still severe; local transmission of both falciparum and vivax malaria persists, and there is a risk of geographically expanding antimalarial resistance. In this research, the pfmdr1, pfcrt, pvmdr1, and K13-propeller genotypes were determined in 26 Plasmodium falciparum and 64 Plasmodium vivax isolates from Yingjiang county of Yunnan province. The pfmdr1 (11.5%), pfcrt (34.6%), and pvmdr1 (3.1%) mutations were prevalent at the China-Myanmar border. The indigenous samples exhibited prevalences of 14.3%, 28.6%, and 14.3% for pfmdr1 N86Y, pfcrt K76T, and pfcrt M74I, respectively, whereas the samples from Myanmar showed prevalences of 10.5%, 21.1%, and 5.3%, respectively. The most prevalent genotypes of pfmdr1 and pfcrt were Y86Y184 and M74N75T76, respectively. No pvmdr1 mutation occurred in the indigenous samples but was observed in two cases coming from Myanmar. In addition, we are the first to report on 10 patients (38.5%) with five different K13 point mutations. The F446I allele is predominant (19.2%), and its prevalence was 28.6% in the indigenous samples of Yingjiang county and 15.8% in samples from Myanmar. The present data might be helpful for enrichment of the molecular surveillance of antimalarial resistance and useful for developing and updating guidance for the use of antimalarials in this region., (Copyright © 2015, American Society for Microbiology. All Rights Reserved.)

Background: The recent emergence and spread of artemisinin resistance in the Greater Mekong Subregion poses a great threat to malaria control and elimination. A K13-propeller gene (K13), PF3D7&#95;1343700, has been associated lately with artemisinin resistance both in vitro and in vivo. This study aimed to investigate the K13 polymorphisms in Plasmodium falciparum parasites from the China-Myanmar border area where artemisinin use has the longest history., Methods: A total of 180 archived P. falciparum isolates containing 191 parasite clones, mainly collected in 2007-2012 from the China-Myanmar area, were used to obtain the full-length K13 gene sequences., Results: Seventeen point mutations were identified in 46.1% (88/191) parasite clones, of which seven were new. The F446I mutation predominated in 27.2% of the parasite clones. The C580Y mutation that is correlated with artemisinin resistance was detected at a low frequency of 1.6%. Collectively, 43.1% of the parasite clones contained point mutations in the kelch domain of the K13 gene. Moreover, there was a trend of increase in the frequency of parasites carrying kelch domain mutations through the years of sample collection. In addition, a microsatellite variation in the N-terminus of the K13 protein was found to have reached a high frequency (69.1%)., Conclusions: This study documented the presence of mutations in the K13 gene in parasite populations from the China-Myanmar border. Mutations present in the kelch domain have become prevalent (>40%). A predominant mutation F446I and a prevalent microsatellite variation in the N-terminus were identified, but their importance in artemisinin resistance remains to be elucidated.

Background: As K13 propeller mutations have been recently reported to serve as molecular markers, assessment of K13 propeller polymorphisms in multidrug-resistant gene in isolates from Myanmar, especially the eastern and western border areas, is crucial if we are to understand the spread of artemisinin resistance., Methods: A 3-day surveillance study was conducted in the eastern and western border areas in Myanmar, and K13 propeller and Plasmodium falciparum multidrug resistance-associated protein 1 (pfmrp1) mutations were analyzed., Results: Among the 1761 suspected malaria cases screened, a total of 42 uncomplicated falciparum cases from the eastern border and 49 from the western border were subjected to 3 days of surveillance after artemether-lumefantrine treatment. No parasitemic case showing positivity on day 3 was noted from the western border, but 26.2% (11/42) of cases were positive in the eastern border. Although we found no marked difference in the prevalence of the pfmrp1 mutation in the eastern and western borders (36% vs 31%, respectively), K13 mutations were more frequent in the eastern border area (where the 3-day persistent cases were detected; 48% vs 14%). C580Y, M476I, A481V, N458Y, R539T, and R516Y accounted for 68.9% of all K13 mutations significantly associated with day 3 parasitaemia., Conclusions: The K13 mutations were significantly associated with day 3 parasitaemia, emphasizing the importance of K13 surveillance. The low prevalence of K13 mutations and the absence of day 3 parasitaemic cases indicate that artemisinin resistance may not have spread to the western Myanmar border region. Although analysis of multiple K13 mutations is challenging, it should be done at various sentinel sites in Myanmar., (© The Author 2014. Published by Oxford University Press on behalf of the Infectious Diseases Society of America. All rights reserved. For Permissions, please e-mail: journals.permissions@oup.com.)

Background: Antigenic variation in Plasmodium falciparum involves switching among multicopy var gene family and is responsible for immune evasion and the maintenance of chronic infections. Current understanding of var gene expression and switching patterns comes from experiments conducted on long laboratory-adapted strains, with little known about their wild counterparts., Methods: Genome sequencing was used to obtain 50 var genes from a parasite isolated from the China-Myanmar border. Four clones with different dominant var genes were cultured in vitro in replicates for 50 generations. Transcription of the individual var gene was detected by real-time PCR and then the switching process was analysed., Results: The expression of multicopy var genes is mutually exclusive in clones of a wild P. falciparum isolate. The activation of distinct primary dominant var genes leads to different and favoured switching patterns in the four clones. The on/off rates of individual var genes are variable and the choice of subsequent dominant var genes are random, which results in the different switching patterns among replicates of each clonal wild P. falciparum isolate with near identical initial transcription profiles., Conclusions: This study suggests that the switching patterns of var genes are abundant, which consist of both conserved and random parts.

The aim of the study was to explore the possible molecular markers of chloroquine resistance in Plasmodium vivax isolates in Thailand. A total of 30 P. vivax isolates were collected from a malaria endemic area along the Thai-Myanmar border in Mae Sot district of Thailand. Dried blood spot samples were collected for analysis of Pvmdr1 and Pvcrt-o polymorphisms. Blood samples (100 μl) were collected by finger-prick for in vitro chloroquine susceptibility testing by schizont maturation inhibition assay. Based on the cut-off IC50 of 100 nM, 19 (63.3%) isolates were classified as chloroquine resistant P. vivax isolates. Seven non-synonymous mutations and 2 synonymous were identified in Pvmdr1 gene. Y976F and F1076L mutations were detected in 7 (23.3%) and 16 isolates (53.3%), respectively. Analysis of Pvcrt-o gene revealed that all isolates were wild-type. Our results suggest that chloroquine resistance gene is now spreading in this area. Monitoring of chloroquine resistant molecular markers provide a useful tool for future control of P. vivax malaria.

Pv12, Pv38 and Pv41, the three 6-Cys family proteins which are expressed in the blood-stage of vivax malaria, might be involved in merozoite invasion activity and thus be potential vaccine candidate antigens of Plasmodium vivax. However, little information is available concerning the genetic diversity and natural selection of these three proteins. In the present study, we analyzed the amino acid sequences of P. vivax blood-stage 6-Cys family proteins in comparison with the homologue proteins of Plasmodium cynomolgi strain B using bioinformatic methods. We also investigated genetic polymorphisms and natural selection of these three genes in P. vivax populations from the China-Myanmar endemic border. The three P. vivax blood-stage 6-Cys proteins were shown to possess a signal peptide at the N-terminus, containing two s48/45 domains, and Pv12 and Pv38 have a GPI-anchor motif at the C-terminus. Then, 22, 21 and 29 haplotypes of pv12, pv38 and pv41 were identified out of 45, 38 and 40 isolates, respectively. The dN/dS values for Domain II of pv38 and pv41 were 3.33880 and 5.99829, respectively, suggesting positive balancing selection for these regions. Meanwhile, the C-terminus of pv41 showed high nucleotide diversity, and Tajima's D test suggested that this fragment could be under positive balancing selection. Overall, our results have significant implications, providing a genetic basis for blood-stage malaria vaccine development based on these three 6-Cys proteins., (Copyright © 2014 Elsevier B.V. All rights reserved.)