Genetic diversity of Plasmodium falciparum populations in southeast and western Myanmar
Parasites & Vectors volume 10, Article number: 322 (2017)
The genetic diversity of malaria parasites reflects the complexity and size of the parasite populations. This study was designed to explore the genetic diversity of Plasmodium falciparum populations collected from two southeastern areas (Shwekyin and Myawaddy bordering Thailand) and one western area (Kyauktaw bordering Bangladesh) of Myanmar.
A total of 267 blood samples collected from patients with acute P. falciparum infections during 2009 and 2010 were used for genotyping at the merozoite surface protein 1 (Msp1), Msp2 and glutamate-rich protein (Glurp) loci.
One hundred and eighty four samples were successfully genotyped at three genes. The allelic distributions of the three genes were all significantly different among three areas. MAD20 and 3D7 were the most prevalent alleles in three areas for Msp1 and Msp2, respectively. The Glurp allele with a bin size of 700–750 bp was the most prevalent both in Shwekyin and Myawaddy, whereas two alleles with bin sizes of 800–850 bp and 900–1000 bp were the most prevalent in the western site Kyauktaw. Overall, 73.91% of samples contained multiclonal infections, resulting in a mean multiplicity of infection (MOI) of 1.94. Interestingly, the MOI level presented a rising trend with the order of Myawaddy, Kyauktaw and Shwekyin, which also paralleled with the increasing frequencies of Msp1 RO33 and Msp2 FC27 200–250 bp alleles. Msp1 and Msp2 genes displayed higher levels of diversity and higher MOI rates than Glurp. PCR revealed four samples (two from Shwekyin and two from Myawaddy) with mixed infections of P. falciparum and P. vivax.
This study genotyped parasite clinical samples from two southeast regions and one western state of Myanmar at the Msp1, Msp2 and Glurp loci, which revealed high levels of genetic diversity and mixed-strain infections of P. falciparum populations at these sites. The results indicated that malaria transmission intensity in these regions remained high and more strengthened control efforts are needed. The genotypic data provided baseline information for monitoring the impacts of malaria elimination efforts in the region.
Malaria caused by P. falciparum is one of the major human infectious diseases affecting millions of people worldwide. Despite enormous control efforts over several decades, malaria still remains an important public health problem. In 2015, there were about 214 million new malaria cases and 438,000 deaths . In Southeast Asia, the Greater Mekong Subregion (GMS) has been one of the most malarious areas, with immense geographical heterogeneity in endemicity and the co-existence of different Plasmodium species [2, 3]. In this region, Myanmar has the greatest malaria burden. With enhanced control efforts by the national malaria control program and support from international partners, malaria in Myanmar has been steadily declining from 1341.8 cases per 100,000 population in 2005 to 253.3 in 2014 . Yet, despite significant reduction of malaria incidence, P. falciparum parasites still maintained a high level of genetic diversity in northeastern Myanmar .
The genotypic and phenotypic diversity of malaria parasites enhances their ability to counteract control measures such as therapeutic drugs. In hyperendemic areas, each individual host may harbor multiple malaria parasite strains with different genotypes, including those that confer drug resistance . Sexual reproduction involving malaria parasites of different genotypes, which can occur when a mosquito feeds on an individual infected with multiple parasite strains or when the same mosquito feeds on more than one human bearing distinct parasite genotypes, favors genetic recombination and generation of higher diversity. Thus, it is expected that parasite genetic diversity tends to correlate with the transmission intensity. Likewise, many studies have shown that malaria reduction as the result of intensified control efforts is accompanied by reduced genetic diversity of the parasite populations [7, 8] and in some cases leads to a clonal structure of the parasite population . Notably, extensive use of certain antimalarial drugs can also alter the genetic diversity as a consequence of selective sweeps. Genetic diversity, malaria demographic history and other factors in a given region may reflect the transmission intensity, effectiveness of malaria control measures, and potential emergence of resistant parasites.
In this study, we assessed the genetic diversity of P. falciparum parasites from three areas in southeast and western Myanmar using three polymorphic markers: merozoite surface protein 1 (Msp1), Msp2, and the glutamate-rich protein (Glurp).
Collection of clinical parasite samples
Plasmodium falciparum clinical samples were collected in 2009–2010 during efficacy studies of artemisinin combination therapies conducted in three locations in Myanmar (Shwekyin, Bago Division; Myawaddy, Kayin State; and Kyauktaw, Rakhine State; Fig. 1). Shwekyin and Myawaddy are located in southeast Myanmar, close to the Myanmar-Thailand border, whereas Kyauktaw is located in western Myanmar near the Myanmar-Bangladesh border. Shwekyin is ~510 km and ~220 km from Kyauktaw and Myawaddy, respectively. These sites were selected for the drug efficacy studies because they represent the malaria high-endemicity regions in Myanmar. In 2009, 112,653, 47,479 and 20,044 malaria cases were recorded in Rakhine State, Bago Division, and Kayin State, respectively. A total of 267 samples (56, 151, and 60 from Myawaddy, Shwekyin and Kyauktaw, respectively) were obtained from patients aged 6–70 years with acute malaria (Additional file 1: Table S1). All cases were confirmed to have P. falciparum infection by microscopic examination of Giemsa-stained thick and thin blood smears at the local clinics. For parasite genotyping, 200 μl of finger-pricked blood was spotted on 3 mm Whatman filter paper, dried, and stored at 4 °C until DNA extraction.
Species identification and genotyping of P. falciparum Msp1, Msp2 and Glurp
Parasite genomic DNA was extracted from filter papers using the QIAamp DNA micro kit (Qiagen, Hilden, Germany) following the manufacturer’s instructions. Nested polymerase chain reaction (PCR) was used to identify the Plasmodium species with primers specific for P. falciparum, P. vivax, P. malariae, P. ovale and P. knowlesi [10,11,12]. The different genotypes of Msp1, Msp2, and Glurp genes were amplified by using a nested PCR method as previously described [5, 13, 14]. The PCR products were separated on 2% agarose gel and stained with ethidium bromide. The size of the PCR products was estimated based on their mobility relative to the standard DNA ladder marker (TaKaRa, Kusatsu, Japan).
Estimating allele frequencies and clone numbers
Alleles were identified based on the type and the fragment size of PCR products and assigned to different size bins. If an isolate had one allele at each of the three loci, this sample was considered to have a monoclonal infection. The multiplicity of infection (MOI) was defined as the largest number of alleles among three loci detected in the sample. This measure is conservative and it likely underestimates the number of clones present . Yet, this may best reflect the level of endemicity in these areas.
The Msp1, Msp2 and Glurp allele frequencies were expressed as percentages. The proportions of alleles observed in each gene among three areas were compared using the Chi-square test. Pairwise comparison of MOI was done by Student’s t-test. The correlation between MOI and percentages of gametocytemic patients was analyzed by using Spearman’s correlation. The P value less than 0.05 was considered statistically significant.
The demographic and clinical characteristics of the study populations are presented in Table 1. Since the samples were collected from drug efficacy trials and the number of patients at each site did not reflect the level of local malaria transmission, the data could not be used for vigorous epidemiological comparisons. Yet, a significantly higher density of asexual parasites was detected from children in Myawaddy. At the time of enrollment to the efficacy studies, 21.19%, 10.71% and 15% patients in Shwekyin, Myawaddy and Kyauktaw carried gametocytes, respectively (Table 1). Although all patients were confirmed to have P. falciparum infections by microscopic examination, nested PCR identified four samples as P. falciparum and P. vivax mixed infections (two from each of the southeast sites). Other Plasmodium species (P. malariae, P. ovale or P. knowlesi) were not detected.
Variants of Msp1, Msp2 and Glurp genes
Of the 267 samples, 184 were successfully genotyped in Msp1, Msp2 and Glurp genes. The PCR fragment sizes and allelic frequencies for Msp1, Msp2 and Glurp in three areas are shown in Table 2. The allelic distributions of Msp1, Msp2 and Glurp genes were significantly different among three areas (Table 2). The most abundant allele types of Msp1 and Msp2 of all three areas were MAD20 and 3D7, respectively. Shwekyin and Myawaddy both shared the same dominant Glurp allele of 700–750 bp, but the most prevalent Glurp alleles in Kyauktaw were 800–850 bp and 900–1000 bp. Compared to Glurp, Msp1 and Msp2 displayed higher genetic diversity in all three areas (Additional file 2: Table S2).
Of the 184 samples successfully genotyped for all three genes, only 26.09% samples were mono- allelic at each of the three genes (Table 3). Analysis of Msp1 and Msp2 revealed 52.72% and 54.35% samples harboring multiclonal infections, respectively, whereas analysis of Glurp identified 13.04% multiclonal infections. When the genotyping results from three loci were combined, 73.91% of the samples were multiclonal infections, giving a mean MOI of 1.94. Among the three areas, the levels of multiclonal infections ranged from 64 to 80.21%, and the MOI ranged from 1.70 to 2.05. The lowest MOI was found in Myawaddy, which was significantly lower than that in Shwekyin (t-test: t (144) = 3.096, P = 0.002). In addition, the MOI did not show significant correlation with the proportion of gametocytemic patients in each site (Spearman correlation: r (184) = 0.1, P = 0.2).
The genetic diversity of P. falciparum parasites impacts malaria transmission and malaria control strategies. Therefore, it is important to resolve the genetic population structure of P. falciparum parasites in epidemic areas. Here we compared the genetic diversity of P. falciparum populations in two southeast regions and one western state of Myanmar using three polymorphic antigen markers. Consistent with earlier findings, genotyping the two Msp markers Msp1 and Msp2 identified significantly higher numbers of allele variants than genotyping Glurp [5, 17]. For Msp1, the MAD20 type was more prevalent than the RO33 and K1 alleles, as reported in other areas of the GMS including southwestern China , northeastern Myanmar  and the Thai-Myanmar border . It was noteworthy that the frequency of the RO33 allele in the present study ranged from 15.28 to 32.76%, which is distinct from the rare presence of this type in parasites from the central Myanmar . For Msp2, the 3D7 type was more abundant than the FC27 type, which was similar to what has been reported from the Thai-Myanmar border , but different from results in the northeastern Myanmar  and Laos . Interestingly, all three sites shared the same most prevalent Msp2 allele (3D7 type, allele size 450–550 bp). For Glurp, Shwekyin and Myawaddy both had the 700–750 bp allele as the most frequent, whereas the 800–850 bp and 900–1000 bp alleles were the most common ones from Kyauktaw. The genetic similarities of the analyzed P. falciparum may reflect the geographical proximity of the study areas and possible population mixing from neighboring areas of the GMS. For example, the genetic diversity of parasites from the southeast site Myawaddy described in this study was similar to that from the Thai-Myanmar border . Overall, the three polymorphic markers revealed high levels of genetic diversity of parasite populations from the three areas. This suggests that despite the malaria control efforts, the parasite’s effective population remained large in these historically high-endemicity areas. Significant sharing of certain alleles may be attributed to extensive human migration along the border areas, which serves as a vehicle for introducing and spreading of parasites between sites .
The level of MOI usually reflects the degree of transmission intensity, although the association is not linear [7, 20]. Genotyping the three polymorphic markers in this study revealed that at least 73.91% of the analyzed samples were mixed infections. The largest number of alleles detected from a single isolate was 4, 3 and 2 for Msp1, Msp2 and Glurp, respectively. Such a high level of mixed infections indicates that transmission intensity was still high in these regions. Moreover, the trend for increased MOI to parallel the increased frequencies of the Msp1 RO33 and the Msp2 FC27 200–250 bp alleles suggested that these allelic types may be used as molecular markers for malaria prevalence or burden assessment in the region. Whereas this study and our earlier work in northeastern Myanmar using the same antigenic markers identified similarly high levels of genetic diversity of the parasites in spite of recent drastic reduction in malaria cases, the northeastern region had significantly lower MOI . As the genetic diversity is expected to decrease in the face of malaria elimination, the information provided herein may serve as the baseline data for continued epidemiological surveillance in these endemic regions.
It is interesting to note that although this study was not designed to detect the prevalence of all malaria parasites, PCR analysis did identify four cases of mixed P. falciparum/ P. vivax infections. These four misdiagnosed samples were not in the same area. Misdiagnosis of the mixed infections is most likely due to lower parasitemia normally associated with P. vivax infections. Given that the missed P. vivax co-infections would preclude treatment with primaquine for anti-relapse therapy, it is likely that some of these patients would experience P. vivax relapses. Therefore, more detailed epidemiological studies involving all human parasites are needed to clarify the prevalence, dynamics and genetics of the parasites in the endemic regions of Myanmar.
This study described the genetic diversity of P. falciparum from two southeast regions and one western state of Myanmar during 2009–2010. Successful genotyping of 184 clinical samples at Msp1, Msp2 and Glurp revealed high levels of genetic diversity and MOI, indicating that P. falciparum transmission remained relatively intense in these regions. These results can serve as the baseline information to monitor the effectiveness of the malaria elimination efforts. We can compare the MOI with future’s data, to assay the effect of the prevention measures.
- Glurp :
Greater Mekong Subregion
multiplicity of infection
- Msp1 :
merozoite surface protein 1
- Msp2 :
merozoite surface protein 2
WHO. World malaria report; 2015. p. 2015.
Cui L, Yan G, Sattabongkot J, Chen B, Cao Y, Fan Q, et al. Challenges and prospects for malaria elimination in the greater Mekong Subregion. Acta Trop. 2012;121(3):240–5.
Cui L, Yan G, Sattabongkot J, Cao Y, Chen B, Chen X, et al. Malaria in the greater Mekong Subregion: heterogeneity and complexity. Acta Trop. 2012;121(3):227–39.
Mu TT, Sein AA, Kyi TT, Min M, Aung NM, Anstey NM, et al. Malaria incidence in Myanmar 2005-2014: steady but fragile progress towards elimination. Malar J. 2016;15(503):1–9.
Yuan L, Zhao H, Wu L, Li X, Parker D, Xu S, et al. Plasmodium falciparum populations from northeastern Myanmar display high levels of genetic diversity at multiple antigenic loci. Acta Trop. 2012;125(1):53–9.
Robinson T, Campino SG, Auburn S, Assefa SA, Polley SD, Manske M, et al. Drug-resistant genotypes and multi-clonality in Plasmodium falciparum analysed by direct genome sequencing from peripheral blood of malaria patients. PLoS One. 2011;6(8):e23204.
Hoffmann EH, da Silveira LA, Tonhosolo R, Pereira FJ, Ribeiro WL, Tonon AP, et al. Geographical patterns of allelic diversity in the Plasmodium falciparum malaria-vaccine candidate, merozoite surface protein-2. Ann Trop Med Parasitol. 2001;95(2):117–32.
Daniels R, Chang HH, Sene PD, Park DC, Neafsey DE, Schaffner SF, et al. Genetic surveillance detects both clonal and epidemic transmission of malaria following enhanced intervention in Senegal. PLoS One. 2013;8(4):e60780.
Putaporntip C, Miao J, Kuamsab N, Sattabongkot J, Sirichaisinthop J, Jongwutiwes S, et al. The Plasmodium vivax merozoite surface protein 3beta sequence reveals contrasting parasite populations in southern and northwestern Thailand. PLoS Negl Trop Dis. 2014;8(11):e3336.
Johnston SP, Pieniazek NJ, Xayavong MV, Slemenda SB, Wilkins PP, da Silva AJ. PCR as a confirmatory technique for laboratory diagnosis of malaria. J Clin Microbiol. 2006;44(3):1087–9.
Buppan P, Putaporntip C, Pattanawong U, Seethamchai S, Jongwutiwes S. Comparative detection of Plasmodium vivax and Plasmodium falciparum DNA in saliva and urine samples from symptomatic malaria patients in a low endemic area. Malar J. 2010;9:72.
Yan J, Li N, Wei X, Li P, Zhao Z, Wang L, et al. Performance of two rapid diagnostic tests for malaria diagnosis at the China-Myanmar border area. Malar J. 2013;12:73.
Kaneko O, Kimura M, Kawamoto F, Ferreira MU, Tanabe K. Plasmodium falciparum: allelic variation in the merozoite surface protein 1 gene in wild isolates from southern Vietnam. Exp Parasitol. 1997;86(1):45–57.
Viriyakosol S, Siripoon N, Petcharapirat C, Petcharapirat P, Jarra W, Thaithong S, et al. Genotyping of Plasmodium falciparum isolates by the polymerase chain reaction and potential uses in epidemiological studies. Bull World Health Organ. 1995;73(1):85–95.
Hill WG, Babiker HA. Estimation of numbers of malaria clones in blood samples. Proc Biol Sci. 1995;262(1365):249–57.
Zhu X, Zhou L, Liu Q, Gao X. Genotype and sequence analysis of merozoite surface protein 1 of Plasmodium falciparum isolates in Yunnan province. Zhongguo Ji Sheng Chong Xue Yu Ji Sheng Chong Bing Za Zhi. 1999;17(3):155–8. (In Chinese)
Congpuong K, Sukaram R, Prompan Y, Dornae A. Genetic diversity of the msp-1, msp-2, and glurp genes of Plasmodium falciparum isolates along the Thai-Myanmar borders. Asian Pac J Trop Biomed. 2014;4(8):598–602.
Kang JM, Moon SU, Kim JY, Cho SH, Lin K, Sohn WM, et al. Genetic polymorphism of merozoite surface protein-1 and merozoite surface protein-2 in Plasmodium falciparum field isolates from Myanmar. Malar J. 2010;9:131.
Khaminsou N, Kritpetcharat O, Daduang J, Charerntanyarak L, Kritpetcharat P. Genetic analysis of the merozoite surface protein-1 block 2 allelic types in Plasmodium falciparum clinical isolates from Lao PDR. Malar J. 2011;10:371.
Paul RE, Hackford I, Brockman A, Muller-Graf C, Price R, Luxemburger C, et al. Transmission intensity and Plasmodium falciparum diversity on the northwestern border of Thailand. Am J Trop Med Hyg. 1998;58(2):195–203.
We thank the staff at the clinics and patients for participation in this study. We also thank Miss Yao Bai and Mr. Weilin Zeng for assisting data analysis.
This study was funded by the National Natural Science Foundation of China (grants U1202226 and 31,260,508); a Talent Introduction Project of Yunnan Province (grant 2013HA026); the National Institutes of Health, USA (grant U19AI089672); the Yunnan Province Research Foundation for Basic Research, China (2014FB005); the Doctor Starting Foundation from Basic Medical Sciences of Kunming Medical University (#1) and the WHO (Ref. No.SE/09/214744; SE.10.11500. MMI. MPRO03. MMK, INV.ERS71780-30409).
Availability of data and materials
All relevant data are within the paper.
Ethics approval and consent to participate
Ethical approval for the study was given by The Institutional Review Board of The Department of Medical Research, Myanmar. Written informed consent was gathered from the participants or their guardians prior to enrolment. The Institutional Review Board of The Department of Medical Research, Myanmar approved the subject protocol in advance of its completion. Written informed consent was gathered from the participants or their guardians prior to enrollment.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Soe, T.N., Wu, Y., Tun, M.W. et al. Genetic diversity of Plasmodium falciparum populations in southeast and western Myanmar. Parasites Vectors 10, 322 (2017). https://doi.org/10.1186/s13071-017-2254-x