Next Article in Journal
Interspecific Genetic Differences and Historical Demography in South American Arowanas (Osteoglossiformes, Osteoglossidae, Osteoglossum)
Previous Article in Journal
Population Genetics of the Highly Polymorphic RPP8 Gene Family
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 10
Issue 9
10.3390/genes10090692
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genetic Variations Associated with Drug Resistance Markers in Asymptomatic Plasmodium falciparum Infections in Myanmar

by
Yan Zhao
1,†,
Ziling Liu
1,†,
Myat Thu Soe
2,†,
Lin Wang
1,
Than Naing Soe
3,
Huanping Wei
1,
Aye Than
2,
Pyae Linn Aung
2,
Yuling Li
1,
Xuexing Zhang
1,
Yubing Hu
1,
Haichao Wei
1,
Yangminghui Zhang
1,
Jessica Burgess
4,
Faiza A. Siddiqui
4,
Lynette Menezes
4,
Qinghui Wang
1,
Myat Phone Kyaw
2,*,
Yaming Cao
1,* and
Liwang Cui
4
1
Department of Immunology, College of Basic Medical Sciences, China Medical University, Shenyang 110122, China
2
Myanmar Health Network Organization, Yangon 11211, Myanmar
3
Department of Public Health, Ministry of Health and Sports, Nay Pyi Taw 15011, Myanmar
4
Department of Internal Medicine, Morsani College of Medicine, University of South Florida, 3720 Spectrum Boulevard, Tampa, FL 33612, USA
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally.
Genes 2019, 10(9), 692; https://doi.org/10.3390/genes10090692
Submission received: 2 August 2019 / Revised: 31 August 2019 / Accepted: 4 September 2019 / Published: 9 September 2019
Browse Figures
Review Reports Versions Notes

Abstract

:
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.

1. Introduction

The emergence and spread of parasites resistant to antimalarial drugs and mosquitoes resistant to insecticides threaten the recent gains in malaria control and challenge the goal for malaria elimination in the Greater Mekong Subregion (GMS) of Southeast Asia [1]. Antimalarial drug resistance in P. falciparum tends to emerge in low-transmission settings, particularly in Southeast Asia and South America, before expanding to high-transmission settings in sub-Saharan Africa [2]. Resistance to chloroquine (CQ) and later to sulfadoxine-pyrimethamine (SP) has been responsible for the increased mortality in African children [3,4]. Artemisinin (ART)-resistant P. falciparum parasites were first reported in Cambodia and subsequently in all countries of the GMS [5]. This has accelerated resistance development in parasites to artemisinin combination therapy (ACT) partner drugs, resulting in increased treatment failure rates with dihydroartemisinin-piperaquine (DP) in Cambodia and with artesunate-mefloquine (AS-MQ) in Cambodia and on the Thai-Myanmar border [6,7]. The timely monitoring of antimalarial resistance spread is essential for maintaining the recent progress in malaria control and for achieving the goal of regional malaria elimination.
As a high malaria-burden country within the GMS, Myanmar uses several first-line treatments for P. falciparum infections, including artemether-lumefantrine, AS-MQ, and DP. However, delayed clearance has been observed in all the three first-line ACTs [8]. Similar to other malaria-endemic areas, Myanmar experiences a large proportion of asymptomatic malaria infections, as demonstrated by active case surveillance of malaria infections among healthy populations [9,10]. In the Shwegyin township of southern Myanmar, mutations in genes associated with drug resistance including those with ART resistance were identified in asymptomatic infections [11]. Asymptomatic or sub-clinical malaria infections seldom cause acute disease, but they are capable of infecting mosquitoes, thus serving as silent reservoirs for continued malaria transmission. Therefore, prevalent antimalarial drug resistant strains in asymptomatic infections may underlie the rapid spread of drug resistance. Given the limited data on drug resistance in asymptomatic infections, ongoing molecular epidemiological studies of drug resistance are needed in these parasite reservoirs.
Molecular epidemiology studies provide information for detecting the emergence and tracking the spread of antimalarial drug resistance. Several mutations in the P. falciparum dihydrofolate reductase (pfdhfr) and P. falciparum dihydropteroate synthase (pfdhps) genes (e.g., triple mutations at codons 51, 59 and 108 of pfdhfr and double mutations at codons 437 and 540 of pfdhps) are associated with SP treatment failures [12]. The P. falciparum chloroquine resistance transporter (pfcrt) K76T mutation and P. falciparum multidrug resistance protein 1 (pfmdr1) N86Y mutation have been linked to CQ and amodiaquine resistance [13]. Mutations in the P. falciparum multidrug resistance-associated protein 1 (pfmrp1), such as H191Y and S437A, were reported to be associated with resistance to CQ and quinine (QN) in vitro [14]. Point mutations in the propeller domain of Kelch protein 13 (k13) were collectively correlated with clinical ART resistance [15,16]. Molecular surveillance showed that k13 mutations associated with ART resistance were restricted to certain areas of the GMS, with C580Y being the predominant mutation in the Thai-Cambodian region and F446I in China and northern Myanmar [17,18,19,20].
Here, we profiled mutations in antimalarial drug resistance genes in asymptomatic P. falciparum infections identified from cross-sectional studies in three townships of Myanmar. The results will provide updated information on the drug resistance status of P. falciparum malaria in multiple sentinel sites in Myanmar, which is necessary for developing strategies to eliminate P. falciparum infections in this country.

2. Materials and Methods

2.1. Ethics Approval

Ethical approval for this project was obtained from the institutional review boards of China Medical University, China (2019086), University of South Florida, USA (Pro00036813) and the Ministry of Health and Sports, Myanmar (Ethics/DMR/2017/077AE/2018). Before conducting the study, all adult participants or legal guardians of children voluntarily signed the informed consent.

2.2. Study Sites and Samples

Three cross-sectional studies were conducted in Myanmar at the Laiza (Kachin State), Banmauk (Sagaing Region), and Paletwa (Chin State) townships in 2015, 2017–2018, and 2017, respectively (Figure 1A). From 2015 to 2017, a country profile of malaria showed that Paletwa had the highest annual incidence of P. falciparum, followed by Banmauk and Laiza [21,22,23]. These cross-sectional surveys recruited a total of 3,495 residents, who provided finger-prick blood to prepare blood smears and dried blood spots on filter paper (Figure 1B).

2.3. Diagnosis of P. falciparum Asymptomatic Infections

Smears were stained with Giemsa and examined by two experienced microscopists. For molecular diagnosis, parasite genomic DNA was isolated from blood spots on filter paper according to the protocol of QIAamp® DNA Mini kit (Qiagen). Detection of P. falciparum isolates was performed by nested PCR using an established protocol targeting the 18S rRNA gene as described previously [24]. This method has a detection limit of ~2 parasites/μL. For PCR assessment, one positive and one negative control (from a symptomatic P. falciparum isolate and sterile water, respectively) were used in each of the amplification plates. The genomic DNA from P. falciparum infections was used for the amplification of genes associated with drug resistance.

2.4. PCR and Sequencing of Drug Resistance Markers

The target sequences covering the resistance-conferring mutations from each gene were selected for PCR amplification and sequencing using published protocols [15,25,26,27,28]. The sample DNA typically had OD260/280 between 1.6 and 1.8 with a concentration of >5 ng/µL as measured using the NanoDrop 2000C spectrophotometer. PCR amplicons were purified and sequenced using an ABI 3730XL DNA analyzer. For sequence accuracy, all DNA fragments were sequenced for both strands. The primers used to amplify these genes are summarized in Table S1 and S2. These selected regions include codons 51, 59, 108, and 164 of pfdhfr (PF3D7_0417200); 436, 437, 540, and 581 codons of pfdhps (PF3D7_0810800); and three pfmdr1 (PF3D7_0523000) fragments including codons 86, 184, 1034, 1042, and 1246. For pfcrt (PF3D7_0709000), pfmrp1 (PF3D7_0112200) and k13 (PF3D7_1343700), full-length genes were amplified and sequenced (Figure S1). The gene sequences for k13, pfdhfr, pfdhps, pfmdr1, pfmrp1 and pfcrt reported in this study were deposited in GenBank under accession numbers MN419439 - MN419894.

2.5. Sequence Analysis and Statistics

We evaluated the quality of sequences by examining the chromatograms. Those showing ambiguity were re-amplified and re-sequenced. If the re-sequencing did not improve the quality, these samples were excluded from the analysis. The sequences were aligned using ClustalW in MEGA7.0.26. The samples showing mixed chromatograms were excluded and only single-peak sequences were considered to be monoclonal P. falciparum infections and were used for single nucleotide polymorphism (SNP) and haplotype analysis (Figure S2). Nucleotide and amino acid positions were numbered according to the 3D7 reference sequences. The SNP data from different genes were analyzed using Microsoft Soft Excel 2007 and SPSS 22.0. Pearson’s Chi-square test or Fisher’s exact test were used to determine statistical significance (p < 0.05). The haplotype network was constructed by DnaSP 6 and Network software 5.01.0 using the median joining algorithm [29].

2.6. Availability of Data and Material

The datasets used and/or analyzed during the current study are available from the corresponding authors upon request.

3. Results

3.1. Asymptomatic Infections in Blood Samples from Cross-Sectional Surveys

Asymptomatic P. falciparum infections were prevalent in malaria-endemic areas of Myanmar [24,30]. Dried blood samples on filter papers from 3495 healthy residents collected in cross-sectional surveys in three regions of Myanmar were screened by nested PCR. A total of 107 asymptomatic P. falciparum infections were identified for an overall prevalence of 3.1%. After excluding the samples with double peaks suggestive of mixed infections, 80 pfdhfr, 82 pfdhps, 75 pfcrt, 80 pfmdr1, 82 pfmrp1, and 57 k13 sequences were used for SNP and haplotype analysis (Figure 1B).

3.2. pfdhfr and pfdhps Mutations Associated with Antifolate Resistance

Mutations in pfdhfr at codons 51, 59, 108, and 164 were present in 60.0%, 90.0%, 92.5%, and 40.0% samples, respectively (Table 1). The majority of the isolates carried mutations in pfdhfr at codon 59 and 108. In Laiza, 100.0% (12/12) of the isolates had all the four mutations except one isolate lacking the N51I mutation. In Banmauk and Paletwa, the S108N mutation was predominant at 88.1% (37/42) and 96.2% (25/26), respectively. Analysis of pfdhfr haplotypes revealed that the quadruple mutant 51I/59R/108N/164L was predominant and showed significant differences in prevalence among the three townships; the highest in Laiza at 91.7% (p < 0.0001; Figure 2A and Table 1). The triple mutant haplotype 51I/59R/108N was predominantly in Banmauk (15/42; 35.7%), whereas the double mutant 59R/108N was more prevalent in Paletwa (Table 1). This suggests a trend of decreased resistance to pyrimethamine from the east to the west borders of Myanmar.
For the pfdhps gene, mutations at codons 436, 540, and 581 were present in Laiza and Paletwa, in 25.0% vs. 55.6%, 100.0% vs. 92.6% and 75.0% vs. 25.9% samples, respectively (Table 2). In Banmauk, isolates carrying mutations at codons 436, 437, 540, and 581 were at 20.9%, 48.8%, 32.5%, and 25.6%, respectively. Importantly, in addition to the point mutation K540E, the K540N variant was found in Laiza and Banmauk for the first time. Comparison of the pfdhps haplotypes revealed significant differences among the three townships. In Laiza, 41.7% (5/12) and 33.3% (4/12) of the isolates carried the double mutant 540N/581G and the variant 540E/581G, respectively, whereas in Banmauk the most common haplotype carried the single mutation 437A at 46.5% (20/43) (Table 2). Overall, the double mutant 436A/540E was most dominant and differed significantly in frequency among the three townships with the highest prevalence in Paletwa at 55.6% (p = 0.003; Figure 2B and Table 2).

3.3. Pfcrt and pfmdr1 Mutations Associated with CQ Resistance

Pfcrt gene was successfully sequenced in 75/106 (70.8%) of the isolates, with clear single peaks covering the amino acids 74, 75, 76, 220, 271, 326, 356, and 371. Except for N326S, all these mutations were highly prevalent, ranging from 78.7% to 84.0% (Table 3). The septuple mutant was the most prevalent haplotype (52.0%, 39/75) (Figure 2C), and it was observed in 75.0% of the isolates in Paletwa, 53.8% in Banmauk, but none in Laiza (p < 0.001). In Laiza, 100% of parasites carried mutations at all these amino acids, and the octuple mutant was also present at lower frequencies in Banmauk (15.4%) and Paletwa (8.3 %), respectively (p < 0.001). The sextuple mutant was present exclusively in three isolates from Paletwa (p = 0.057).
Pfmdr1 mutations at codons 86, 130, 184, 185, and 1225 were observed in this study (Table 4). The I185K mutation has been rarely reported previously, but it was the dominant mutation (87.5%) in our samples. In comparison, other mutations were found at lower frequencies, ranging from 1.3% to 23.8%. The I185K and the Y184F mutations differed significantly in frequency among the three townships (p < 0.001 & p = 0.002 respectively). N86Y, a mutation associated with CQ resistance, was only found in Paletwa (7.4%, 2/27). The prevalence of pfmdr1 haplotypes varied in different study sites. The single mutant (b) was the most frequent haplotype in Banmauk (80.5%, 33/41) and Paletwa (74.1%, 20/27), respectively. However, the single mutant (a) was the predominant haplotype in Laiza (58.3%, 7/12). Notably, one isolate from Paletwa carried the pfmdr1 N86Y mutation along with the pfcrt K76T mutation.

3.4. Pfmrp1 Mutations

For the pfmrp1 gene, the mutations at positions 191, 325, 437, 572, 785, 876, 1007, 1339, and 1390 were detected in 65.9%, 2.4%, 65.9%, 2.4%, 22.0%, 59.8%, 50.0%, 2.4% and 8.5% of the asymptomatic P. falciparum isolates, respectively (Table 5). The prevalence of mutations H191Y, S437A, H785N, I876V, and T1007M was significantly different among the three townships (p = 0.002, p = 0.002, p < 0.001, p = 0.029 and p < 0.001, respectively). Ten distinct haplotypes were detected in pfmrp1, with the wild-type (WT) being the most prevalent haplotype, present in 31.7% of the isolates (Figure 2E and Table 5). The quintuple mutant (a) at 50.0% and quadruple mutant (b) at 77.8% were the predominant haplotypes in Laiza and Paletwa, respectively, whereas 46.5% of the isolates had the WT haplotype in Banmauk (Table 5).

3.5. K13 Mutations Associated with ART Resistance

The k13 gene was successfully sequenced in 57 samples. Three different point mutations were found, including two previously reported mutations (F446I and P574L) and one mutation (K191T) observed for the first time in these samples (Table 6). Of these mutations, F446I was most frequent (14.0%; 8/57). Besides point mutations, N, NN, and NNN insertions between amino acids 136 and 137 were also identified in 1 (1.8%), 26 (45.6%) and 4 (7.0%) samples, respectively. All P. falciparum isolates from Laiza had the NN insert, whereas in Paletwa and Banmauk, it was present in 40.0% and 33.3% of the samples, respectively (p < 0.001). The F446I mutation always occurred together with the NN insert and was present exclusively in samples from Laiza along the China-Myanmar border. Of the 57 asymptomatic P. falciparum isolates, 22 samples were WT (Figure 2F and Table 6).

3.6. Haplotype Network

SNPs in the six genotyped drug-resistance genes gave rise to 43 distinct haplotypes in 52 P. falciparum isolates that were successfully sequenced at all six target genes (Figure S3), demonstrating high-level genetic diversity of asymptomatic P. falciparum infections in these regions. Among the samples, 11.6% (5/52), 7.0% (3/52), 7.0% (3/52), and 4.7% (2/52) shared haplotypes No. 2, 9, 30, and 26, respectively, whereas each of the remaining isolates was unique. To examine the phylogenetic relationship of the asymptomatic P. falciparum parasites, a haplotype network was generated based on the SNPs observed in the six genes (Figure 3). The network had three main branches. Isolates of the three study areas appeared to be polyphyletic and were distributed across the entire network, suggesting independent origins of resistant mutations in different genes from diverse genetic backgrounds. Three haplotypes with the k13 K191T mutation originated from Banmauk. Four haplotypes with the k13 F446I mutation and one haplotype with the P574L mutation associated with ART resistance appeared independently in Laiza. Significantly, there were 20 haplotypes with the NN insert in k13 gene, including eight haplotypes in Laiza, six in Banmauk and six in Paletwa.

4. Discussion

As the GMS is moving towards malaria elimination, Myanmar requires immediate attention because of the high malaria burden and its geographical connection with South Asia, through which drug resistance could spread rapidly. According to WHO reports, the latest malaria drug resistance map reveals that both treatment failure and resistant genotypes are prevalent in Myanmar, including our three study sites [31]. Irrespective of the ACTs used, from 2010 to 2018, treatment failure among P. falciparum patients in Kachin State was higher than that in Chin State. Our study found a higher prevalence of k13 and pfmdr1 mutations in Kachin State (east) than in Chin State (west) which is possibly contributing to the treatment failure in Kachin State. Asymptomatic Plasmodium infections, as silent reservoirs of malaria parasites, play an important role in continued malaria transmission. Critically, these asymptomatic infections may carry genetic variations conferring drug resistance, thus facilitating the spread of drug-resistant parasites [11]. Therefore, this study aimed to provide molecular epidemiology information about drug-resistant P. falciparum in asymptomatic infections at three sentinel sites located in southwestern, northern and northeastern Myanmar, each with a different malaria epidemiology.
The drug combination SP, commercially known as Fansidar, has been used as an antimalarial chemotherapy in Southeast Asia and Africa. The mechanism of resistance to sulfadoxine has been associated with five point mutations at the S436A/F, G437A, K540E, A581G, and A613T/S codons of the pfdhps gene. Mutations at 437 and 540 were strongly associated with SP treatment failure, whereas the 436, 581, and 613 mutations confer some degrees of resistance [32]. In this study, A613T/S was not observed, but K540N was prevalent, consistent with a previous study from northeast Myanmar [33]. A581G, G437A, and K540E were the predominant mutations in Laiza, Banmauk, and Paletwa, respectively. In population studies, mutations at codon 59 of the pfdhfr gene and codon 540 of the pfdhps gene are strongly predictive of SP treatment failure [34,35]. The prevalence of this mutation combination in Laiza and Paletwa was 100% and 83.3%, respectively, and was slightly lower in Banmauk (34.1%), suggesting that strong selection pressure against SP was emerging in western and northeastern Myanmar. Similarly, triple mutations at codons 108, 51, and 59 of pfdhfr and double mutations at codons 437 and 540 of pfdhps are associated with SP treatment failure [12]. Here, one isolate was found to have this quintuple mutation combination in Banmauk, and 12 isolates (nine in Laiza, two in Banmauk and one in Paletwa) harbored the sextuple mutations with pfdhfr (N51I, C59R, S108N, and I164L) and pfdhps (K540N/E and A581G), suggesting that P. falciparum from these asymptomatic carriers would be highly resistant to SP.
Mutations in the pfcrt gene, especially the K76T mutation, are major determinants of CQ resistance [36,37]. In addition, pfcrt mutations also influence P. falciparum susceptibility to mefloquine (MQ), halofantrine (HF), and ART [38]. The K76T mutation is associated with different sets of mutations at different codons, most commonly C72S, M74I, N75E, A220S, Q271E, N326S, I356T, and R371I [39]. Consistent with a previous study conducted in Thailand [25], C72S was not found. Except for N326S, other mutations reached high prevalence in all three sites. Notably, almost all parasite isolates had the CIVIET haplotype around codons 72-76, which is the typical CQ-resistant haplotype in Southeast Asia. The stable and highly prevalent of pfcrt mutations may be the result of continued use of CQ as the first-line treatment for P. vivax infections in this region [40,41]. In addition to CQ resistance, pfcrt has gained recognition as a multidrug resistance transporter, which influences parasites’ susceptibilities to multiple first-line antimalarial drugs [42,43]. Emerging mutations in pfcrt, H97Y, F145I, M343L, and G353V, also have been linked to piperaquine (PPQ) resistance [44,45], but these mutations were not identified in this study. Thus, the high prevalence of pfcrt mutations in asymptomatic P. falciparum parasites may be linked to failures of multiple ACTs in Myanmar.
The ATP-binding cassette (ABC) transporters, including pfmdr1 and pfmrp1, are potentially involved in resistance to multiple antimalarial drugs [46,47]. Several point mutations in pfmdr1 are associated with changed sensitivities to CQ, MQ, QN, HF, and ART [38,48,49,50]. In the presence of pfcrt K76T, point mutations in pfmdr1, primarily at codon 86 [13,14] and additionally at positions 184, 1034, 1042, and 1246 [51], can increase resistance of the parasites to CQ. Decreased susceptibility to lumefantrine (LMF) has also been linked to polymorphism in pfcrt and pfmdr1 [52,53]. In this study, only two P. falciparum isolates from the western township Paletwa had the pfcrt K76T and pfmdr1 N86Y mutations, but the combination of pfcrt K76T and pfmdr1 Y184F had much higher prevalence in Laiza (58.3%) than Banmauk (18.0%) and Paletwa (12.5%). The pfmdr1 Y184F was not detected in an earlier study [54], where I185K was reported as a novel mutation. The authors speculated that the decrease in parasite’s susceptibility to the drug was associated with substitution of mutation at codon Y184F for I185K. In our study, the pfmdr1 I185K mutation was highly prevalent (~90% in all samples) and reached fixation in Banmauk and Paletwa. It would be worthwhile to establish whether this mutation is linked to altered drug sensitivity through genetic studies given the high prevalence of this mutation. Since disruption of pfmrp1 in a CQ-resistant P. falciparum isolate rendered this parasite more sensitive to CQ, QN, PPQ, primaquine, and ART, PfMRP1 is proposed to mediate resistance of the parasite to multiple antimalarial drugs [55]. Association studies suggest that mutations at codons Y191H, A437S, H785N, I876V and F1390I in pfmrp1 contribute to decreased sensitivity to CQ, PPQ, LMF, artesunate (ATS) and dihydroartemisinin (DHA) as well as reduce both in vivo and in vitro susceptibility to ACT [28,33,56]. Our results suggest that P. falciparum isolates from asymptomatic infections in these areas might be resistant to CQ, PPQ, LMF, ATS, DHA, and ACT.
ART resistance in P. falciparum has emerged in western Cambodia and its potential spread to Africa threatens global malaria control [57]. ART resistance is characterized by delayed parasite clearance, which corresponds to the decreased susceptibility of ring-stage parasites [58]. Non-synonymous SNPs in K13 were found to be strongly associated with resistance to ART [16]. To date, more than 200 non-synonymous mutations in the k13 gene have been reported. N458Y, Y493H, R539T, I543T, and C580Y are validated ART resistance mutations, and a number of mutations are candidates for ART resistance, including P574L [8]. In Africa, non-synonymous k13 mutations are relatively rare [59]. Recently, k13 mutations have emerged independently in multiple locations in Southeast Asia, including Myanmar [17]. C580Y, Y493H, R539T, and I543T k13 mutations were prevalent in the Thai-Myanmar border, whereas F446I was the most prevalent mutant allele of k13 in the northern Myanmar and China-Myanmar border, and P574L next to C580Y were reported among isolates collected from migrant goldmine workers in southern Myanmar [19,60,61,62]. In our study, the validated mutations of k13 resistance were not found, while F446I, P574L, and K191T mutations in the k13 propeller domain were observed among parasites from asymptomatic infections. The K191T has not been reported in previous studies. F446I mutation was observed with high prevalence only at Laiza along China-Myanmar border, which is consistent with previous studies. F446I mutation in k13 has been found to be associated with delayed parasite clearance [61], and a recent study showed that F446I was associated with the increased ring survival rates by genetically introducing the mutation into the k13 gene [63]. Consistent with previous studies of P. falciparum from a symptomatic population, the absence of k13 propeller mutations in Paletwa suggest that ART resistant parasites may not have spread to the western Myanmar region [18]. Moreover, we found the N, NN, and NNN insertion in the N terminus of k13 protein through analysis of the full-length k13 sequences. A previous study reported that an NN insertion in the N terminus of k13 protein was associated with increased in vitro ring-stage survival in the presence of DHA, but the possibility of the NN insert acting cooperatively with F446I in conferring ART resistance was not ruled out [33]. In view of the high prevalence (45.6%) of the NN insertion among P. falciparum from asymptomatic populations in Myanmar, it would be valuable to determine whether it independently affects ART resistance.
The generation of the haplotype network can be used to study the phylogenetic relationship of parasites, to determine the origins of drug resistance gene mutations from different geographical regions, and to assess the diversity of parasite populations. Consistent with the results from clinical P. falciparum isolates, parasites from asymptomatic individuals also displayed extraordinarily high diversity in Myanmar. The genotypes of the six genes associated with drug resistance indicated that the asymptomatic P. falciparum infections in Myanmar are resistant to multiple drugs. More importantly, 58.8% (20/34) haplotypes carried k13 mutations (N458Y, R539T, P574L, and F446I) associated with ART resistance in clinical isolates, whereas 62.5% (5/8) haplotypes of asymptomatic P. falciparum isolates carried F446I and P574L. Although only two k13 mutations related to ART resistance were detected in our study, the presence of such high haplotype diversity in the limited isolates of asymptomatic infection should be concerning. Despite reduced P. falciparum prevalence in recent years, the evidence suggests that ART resistant P. falciparum parasites appear independently and spread widely along the China-Myanmar border, which is a reminder that surveillance of P. falciparum should be strengthened both in Myanmar and China. Although k13 mutation associated with ART resistance was not found in the other two areas, the diversity of multiple drug resistance reflected by haplotype analysis should not be ignored.

5. Conclusions

In this study, analysis of six molecular markers of drug resistance among asymptomatic P. falciparum infections revealed a high prevalence of mutations linked to drug resistance in Myanmar. High mutation frequency in pfdhfr, pfdhps, and pfcrt means a serious resistance to SP and CQ in these areas, which supports the ACT being used for P. falciparum treatment in Myanmar. Fortunately, ART resistance has not yet spread to western Myanmar. Given that pfcrt also is a multidrug resistance transporter with high frequency of mutations, it is necessary to continuously monitor for the emergence of ART resistance in western Myanmar. However, the drug resistance pattern of P. falciparum from asymptomatic carriers along the China-Myanmar border is not optimistic; perhaps the best option to prevent the spread of drug-resistant parasites in border areas is rapid elimination through mass drug administration. Haplotype analysis also revealed extremely high diversity in these areas. These asymptomatic P. falciparum infections carrying mutations associated with drug resistance may be an important reservoir for the transmission and spread of resistance in this region. Molecular surveillance of antimalarial resistance might be helpful for developing and updating guidance for the use of antimalarials in Myanmar.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4425/10/9/692/s1, Table S1: PCR primer sequences for the amplification of sequences containing P. falciparum pfk13, pfcrt and pfmdr1 genes, Table S2: PCR primer sequences for the amplification of sequences containing P. falciparum pfdhfr, pfdhps, and pfmrp1 genes, Figure S1: Schematic diagram showing successfully sequenced regions from different genes, Figure S2: Sequencing chromatograms showing single and mixed infections, Figure S3: Phylogenetic clustering of the haplotypes based on the mutations in six genes associated with drug resistance of the asymptomatic P. falciparum isolates.

Author Contributions

Conceptualization, Y.Z. (Yan Zhao) and Y.C.; methodology, Y.Z. (Yan Zhao), Z.L., L.W., Y.L., X.Z., Y.H., H.W. (Haichao Wei), and Y.Z. (Yangminghui Zhang); software, M.T.S. and H.W. (Huanping Wei); validation, Y.Z. (Yan Zhao) and L.W.; formal analysis, Y.Z. (Yan Zhao) and Z.L. Investigation, Y.Z. (Yan Zhao), M.T.S., T.N.S., A.T., and P.L.A.; resources, M.T.S., T.N.S., and M.P.K.; data curation, Y.Z. (Yan Zhao) and M.T.S.; writing—original draft preparation, Y.Z. (Yan Zhao), Z.L., and M.T.S.; writing—review and editing, J.B., F.A.S., L.M., and L.C.; visualization, Q.W. and Y.C.; supervision, Q.W. and Y.C.; project administration, M.P.K.; funding acquisition, L.C.

Funding

This study was supported by a grant from the National Institute of Allergy and Infectious Diseases, National Institutes of Health (U19AI089672).

Acknowledgments

We thank all the participants for their willingness to take part in this study. We are indebted to the technicians, nurses, and other workers from local hospitals for their assistance in completing the field surveys and sample collection.

Conflicts of Interest

The authors declare that they have no competing interests. The funders had no role in the design of the study, in the collection, analyses, or interpretation of data, in the writing of the manuscript, or in the decision to publish the results.

Abbreviations

DHFRdihydrofolate reductase
DHPSdihydropteroate synthase
CRTchloroquine resistance transporter
MDR1multidrug resistance protein 1
MRP1multidrug resistance-associated protein 1
K13Kelch 13
GMSGreater Mekong Subregion
ARTartemisinin
CQchloroquine
QNquinine
SPsulfadoxine-pyrimethamine
ACTartemisinin combination therapy
DPdihydroartemisinin-piperaquine
AS-MQartesunate-mefloquine
PCRpolymerase chain reaction
SNPsingle nucleotide polymorphism
WTwild-type
MQmefloquine
HFhalofantrine
ATPadenosine triphosphate
ABCATP-binding cassette
LMFlumefantrine
ATSartesunate and DHA, dihydroartemisinin

References

  1. WHO. Strategy for Malaria Elimination in the Greater Mekong Subregion (2015–2030). 2015. Available online: http://iris.wpro.who.int/handle/10665.1/10945 (accessed on 17 July 2019).
  2. White, N.J. Antimalarial drug resistance. J. Clin. Investig. 2004, 113, 1084–1092. [Google Scholar] [CrossRef] [PubMed]
  3. Trape, J.F. The public health impact of chloroquine resistance in Africa. Am. J. Trop. Med. Hyg. 2001, 64, 12–17. Available online: http://www.ncbi.nlm.nih.gov/pubmed/11425173 (accessed on 17 July 2019). [CrossRef] [PubMed]
  4. Korenromp, E.L.; Williams, B.G.; Gouws, E.; Dye, C.; Snow, R.W. Measurement of trends in childhood malaria mortality in Africa: An assessment of progress toward targets based on verbal autopsy. Lancet Infect. Dis. 2003, 3, 349–358. Available online: http://www.ncbi.nlm.nih.gov/pubmed/12781507 (accessed on 17 July 2019). [CrossRef]
  5. WHO. Global Report on Antimalarial Drug Efficacy and Drug Resistance: 2000–2010. 2010. Available online: https://www.who.int/malaria/publications/atoz/9789241500470/en/ (accessed on 17 July 2019).
  6. Leang, R.; Barrette, A.; Bouth, D.M.; Menard, D.; Abdur, R.; Duong, S.; Ringwald, P. Efficacy of dihydroartemisinin-piperaquine for treatment of uncomplicated Plasmodium falciparum and Plasmodium vivax in Cambodia, 2008 to 2010. Antimicrob. Agents Chemother. 2013, 57, 818–826. [Google Scholar] [CrossRef] [PubMed]
  7. Carrara, V.I.; Lwin, K.M.; Phyo, A.P.; Ashley, E.; Wiladphaingern, J.; Sriprawat, K.; Rijken, M.; Boel, M.; McGready, R.; Proux, S.; et al. Malaria burden and artemisinin resistance in the mobile and migrant population on the Thai-Myanmar border, 1999–2011: An observational study. PLoS Med. 2013, 10, e1001398. [Google Scholar] [CrossRef] [PubMed]
  8. WHO. Status Report on Artemisinin Resistance and ACT Efficacy. 2018. Available online: https://www.who.int/malaria/publications/atoz/artemisinin-resistance-august2018/en/ (accessed on 17 July 2019).
  9. Zhao, Y.; Zeng, J.; Zhao, Y.; Liu, Q.; He, Y.; Zhang, J.; Yang, Z.; Fan, Q.; Wang, Q.; Cui, L.; et al. Risk factors for asymptomatic malaria infections from seasonal cross-sectional surveys along the China-Myanmar border. Malar. J. 2018, 17, 247. [Google Scholar] [CrossRef] [PubMed]
  10. Zaw, M.T.; Thant, M.; Hlaing, T.M.; Aung, N.Z.; Thu, M.; Phumchuea, K.; Phusri, K.; Saeseu, T.; Yorsaeng, R.; Nguitragool, W.; et al. Asymptomatic and sub-microscopic malaria infection in Kayah State, eastern Myanmar. Malar. J. 2017, 16, 138. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Nyunt, M.H.; Shein, T.; Zaw, N.N.; Han, S.S.; Muh, F.; Lee, S.K.; Han, J.H.; Thant, K.Z.; Han, E.T.; Kyaw, M.P. Molecular Evidence of Drug Resistance in Asymptomatic Malaria Infections, Myanmar, 2015. Emerg. Infect. Dis. 2017, 23, 517–520. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Kublin, J.G.; Dzinjalamala, F.K.; Kamwendo, D.D.; Malkin, E.M.; Cortese, J.F.; Martino, L.M.; Mukadam, R.A.; Rogerson, S.J.; Lescano, A.G.; Molyneux, M.E.; et al. Molecular markers for failure of sulfadoxine-pyrimethamine and chlorproguanil-dapsone treatment of Plasmodium falciparum malaria. J. Infect. Dis. 2002, 185, 380–388. [Google Scholar] [CrossRef]
  13. Babiker, H.A.; Pringle, S.J.; Abdel-Muhsin, A.; Mackinnon, M.; Hunt, P.; Walliker, D. High-level chloroquine resistance in Sudanese isolates of Plasmodium falciparum is associated with mutations in the chloroquine resistance transporter gene pfcrt and the multidrug resistance Gene pfmdr1. J. Infect. Dis. 2001, 183, 1535–1538. [Google Scholar] [CrossRef]
  14. Mu, J.; Ferdig, M.T.; Feng, X.; Joy, D.A.; Duan, J.; Furuya, T.; Subramanian, G.; Aravind, L.; Cooper, R.A.; Wootton, J.C.; et al. Multiple transporters associated with malaria parasite responses to chloroquine and quinine. Mol. Microbiol. 2003, 49, 977–989. Available online: http://www.ncbi.nlm.nih.gov/pubmed/12890022 (accessed on 10 July 2019). [CrossRef] [PubMed] [Green Version]
  15. Ashley, E.A.; Dhorda, M.; Fairhurst, R.M.; Amaratunga, C.; Lim, P.; Suon, S.; Sreng, S.; Anderson, J.M.; Mao, S.; Sam, B.; et al. Spread of artemisinin resistance in Plasmodium falciparum malaria. N. Engl. J. Med. 2014, 371, 411–423. [Google Scholar] [CrossRef] [PubMed]
  16. Ariey, F.; Witkowski, B.; Amaratunga, C.; Beghain, J.; Langlois, A.C.; Khim, N.; Kim, S.; Duru, V.; Bouchier, C.; Ma, L.; et al. A molecular marker of artemisinin-resistant Plasmodium falciparum malaria. Nature 2014, 505, 50–55. [Google Scholar] [CrossRef]
  17. Takala-Harrison, S.; Jacob, C.G.; Arze, C.; Cummings, M.P.; Silva, J.C.; Dondorp, A.M.; Fukuda, M.M.; Hien, T.T.; Mayxay, M.; Noedl, H.; et al. Independent emergence of artemisinin resistance mutations among Plasmodium falciparum in Southeast Asia. J. Infect. Dis. 2015, 211, 670–679. [Google Scholar] [CrossRef] [PubMed]
  18. Nyunt, M.H.; Hlaing, T.; Oo, H.W.; Tin-Oo, L.L.; Phway, H.P.; Wang, B.; Zaw, N.N.; Han, S.S.; Tun, T.; San, K.K.; et al. Molecular assessment of artemisinin resistance markers, polymorphisms in the k13 propeller, and a multidrug-resistance gene in the eastern and western border areas of Myanmar. Clin. Infect. Dis. Off. Publ. Infect. Dis. Soc. Am. 2015, 60, 1208–1215. [Google Scholar] [CrossRef] [PubMed]
  19. Wang, Z.; Shrestha, S.; Li, X.; Miao, J.; Yuan, L.; Cabrera, M.; Grube, C.; Yang, Z.; Cui, L. Prevalence of K13-propeller polymorphisms in Plasmodium falciparum from China-Myanmar border in 2007–2012. Malar. J. 2015, 14, 168. [Google Scholar] [CrossRef]
  20. Ye, R.; Hu, D.; Zhang, Y.; Huang, Y.; Sun, X.; Wang, J.; Chen, X.; Zhou, H.; Zhang, D.; Mungthin, M.; et al. Distinctive origin of artemisinin-resistant Plasmodium falciparum on the China-Myanmar border. Sci. Rep. 2016, 6, 20100. [Google Scholar] [CrossRef] [Green Version]
  21. WHO. Malaria Country Profiles-Myanmar. 2017. Available online: https://www.who.int/malaria/publications/country-profiles/2017/profile_mmr_en.pdf?ua=1 (accessed on 4 July 2019).
  22. WHO. Malaria Country Profiles-Myanmar. 2016. Available online: https://www.who.int/malaria/publications/country-profiles/2016/profile_mmr_en.pdf?ua=1 (accessed on 4 July 2019).
  23. WHO. Malaria Country Profiles-Myanmar. 2015. Available online: https://www.who.int/malaria/publications/country-profiles/2015/profile_mmr_en.pdf?ua=1 (accessed on 4 July 2019).
  24. Zhao, Y.; Zhao, Y.; Lv, Y.; Liu, F.; Wang, Q.; Li, P.; Zhao, Z.; Liu, Y.; Cui, L.; Fan, Q.; et al. Comparison of methods for detecting asymptomatic malaria infections in the China-Myanmar border area. Malar. J. 2017, 16, 159. [Google Scholar] [CrossRef]
  25. Chaijaroenkul, W.; Ward, S.A.; Mungthin, M.; Johnson, D.; Owen, A.; Bray, P.G.; Na-Bangchang, K. Sequence and gene expression of chloroquine resistance transporter (pfcrt) in the association of in vitro drugs resistance of Plasmodium falciparum. Malar. J. 2011, 10, 42. [Google Scholar] [CrossRef]
  26. Das, S.; Mahapatra, S.K.; Tripathy, S.; Chattopadhyay, S.; Dash, S.K.; Mandal, D.; Das, B.; Hati, A.K.; Roy, S. Double mutation in the pfmdr1 gene is associated with emergence of chloroquine-resistant Plasmodium falciparum malaria in Eastern India. Antimicrob. Agents Chemother. 2014, 58, 5909–5915. [Google Scholar] [CrossRef]
  27. Yang, Z.; Li, C.; Miao, M.; Zhang, Z.; Sun, X.; Meng, H.; Li, J.; Fan, Q.; Cui, L. Multidrug-resistant genotypes of Plasmodium falciparum, Myanmar. Emerg. Infect. Dis. 2011, 17, 498–501. [Google Scholar] [CrossRef] [PubMed]
  28. Dahlstrom, S.; Ferreira, P.E.; Veiga, M.I.; Sedighi, N.; Wiklund, L.; Martensson, A.; Farnert, A.; Sisowath, C.; Osorio, L.; Darban, H.; et al. Plasmodium falciparum multidrug resistance protein 1 and artemisinin-based combination therapy in Africa. J. Infect. Dis. 2009, 200, 1456–1464. [Google Scholar] [CrossRef] [PubMed]
  29. Bandelt, H.J.; Forster, P.; Rohl, A. Median-joining networks for inferring intraspecific phylogenies. Mol. Biol. Evol. 1999, 16, 37–48. [Google Scholar] [CrossRef] [PubMed]
  30. Liu, Z.; Soe, T.N.; Zhao, Y.; Than, A.; Cho, C.; Aung, P.L.; Li, Y.; Wang, L.; Yang, H.; Li, X.; et al. Geographical heterogeneity in prevalence of subclinical malaria infections at sentinel endemic sites of Myanmar. Parasites Vectors 2019, 12, 83. [Google Scholar] [CrossRef] [PubMed]
  31. WHO. Malaria Threats Map: Parasite Drug Resistance. 2019. Available online: http://apps.who.int/malaria/maps/threats/ (accessed on 17 July 2019).
  32. Alker, A.P.; Kazadi, W.M.; Kutelemeni, A.K.; Bloland, P.B.; Tshefu, A.K.; Meshnick, S.R. Dhfr and dhps genotype and sulfadoxine-pyrimethamine treatment failure in children with falciparum malaria in the Democratic Republic of Congo. Trop. Med. Int. Health TM IH 2008, 13, 1384–1391. [Google Scholar] [CrossRef] [PubMed]
  33. Wang, Z.; Wang, Y.; Cabrera, M.; Zhang, Y.; Gupta, B.; Wu, Y.; Kemirembe, K.; Hu, Y.; Liang, X.; Brashear, A.; et al. Artemisinin resistance at the China-Myanmar border and association with mutations in the K13 propeller gene. Antimicrob. Agents Chemother. 2015, 59, 6952–6959. [Google Scholar] [CrossRef] [PubMed]
  34. Lau, T.Y.; Sylvi, M.; William, T. Mutational analysis of Plasmodium falciparum dihydrofolate reductase and dihydropteroate synthase genes in the interior division of Sabah, Malaysia. Malar. J. 2013, 12, 445. [Google Scholar] [CrossRef]
  35. Vinayak, S.; Alam, M.T.; Mixson-Hayden, T.; McCollum, A.M.; Sem, R.; Shah, N.K.; Lim, P.; Muth, S.; Rogers, W.O.; Fandeur, T.; et al. Origin and evolution of sulfadoxine resistant Plasmodium falciparum. PLoS Pathog. 2010, 6, e1000830. [Google Scholar] [CrossRef]
  36. Djimde, A.; Doumbo, O.K.; Cortese, J.F.; Kayentao, K.; Doumbo, S.; Diourte, Y.; Coulibaly, D.; Dicko, A.; Su, X.Z.; Nomura, T.; et al. A molecular marker for chloroquine-resistant falciparum malaria. N. Engl. J. Med. 2001, 344, 257–263. [Google Scholar] [CrossRef]
  37. Djimde, A.A.; Barger, B.; Kone, A.; Beavogui, A.H.; Tekete, M.; Fofana, B.; Dara, A.; Maiga, H.; Dembele, D.; Toure, S.; et al. A molecular map of chloroquine resistance in Mali. FEMS Immunol. Med Microbiol. 2010, 58, 113–118. [Google Scholar] [CrossRef] [Green Version]
  38. Duraisingh, M.T.; Jones, P.; Sambou, I.; von Seidlein, L.; Pinder, M.; Warhurst, D.C. The tyrosine-86 allele of the pfmdr1 gene of Plasmodium falciparum is associated with increased sensitivity to the anti-malarials mefloquine and artemisinin. Mol. Biochem. Parasitol. 2000, 108, 13–23. Available online: http://www.ncbi.nlm.nih.gov/pubmed/10802315 (accessed on 10 July 2019). [CrossRef]
  39. Ecker, A.; Lehane, A.M.; Clain, J.; Fidock, D.A. PfCRT and its role in antimalarial drug resistance. Trends Parasitol. 2012, 28, 504–514. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Brown, T.; Smith, L.S.; Oo, E.K.; Shawng, K.; Lee, T.J.; Sullivan, D.; Beyrer, C.; Richards, A.K. Molecular surveillance for drug-resistant Plasmodium falciparum in clinical and subclinical populations from three border regions of Burma/Myanmar: Cross-sectional data and a systematic review of resistance studies. Malar. J. 2012, 11, 333. [Google Scholar] [CrossRef] [PubMed]
  41. Guthmann, J.P.; Pittet, A.; Lesage, A.; Imwong, M.; Lindegardh, N.; Min Lwin, M.; Zaw, T.; Annerberg, A.; de Radigues, X.; Nosten, F. Plasmodium vivax resistance to chloroquine in Dawei, southern Myanmar. Trop. Med. Int. Health TM IH 2008, 13, 91–98. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Park, D.J.; Lukens, A.K.; Neafsey, D.E.; Schaffner, S.F.; Chang, H.H.; Valim, C.; Ribacke, U.; Van Tyne, D.; Galinsky, K.; Galligan, M.; et al. Sequence-based association and selection scans identify drug resistance loci in the Plasmodium falciparum malaria parasite. Proc. Natl. Acad. Sci. USA 2012, 109, 13052–13057. [Google Scholar] [CrossRef]
  43. Miotto, O.; Amato, R.; Ashley, E.A.; MacInnis, B.; Almagro-Garcia, J.; Amaratunga, C.; Lim, P.; Mead, D.; Oyola, S.O.; Dhorda, M.; et al. Genetic architecture of artemisinin-resistant Plasmodium falciparum. Nat. Genet. 2015, 47, 226–234. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Ross, L.S.; Dhingra, S.K.; Mok, S.; Yeo, T.; Wicht, K.J.; Kumpornsin, K.; Takala-Harrison, S.; Witkowski, B.; Fairhurst, R.M.; Ariey, F.; et al. Emerging Southeast Asian PfCRT mutations confer Plasmodium falciparum resistance to the first-line antimalarial piperaquine. Nat. Commun. 2018, 9, 3314. [Google Scholar] [CrossRef] [Green Version]
  45. Agrawal, S.; Moser, K.A.; Morton, L.; Cummings, M.P.; Parihar, A.; Dwivedi, A.; Shetty, A.C.; Drabek, E.F.; Jacob, C.G.; Henrich, P.P.; et al. Association of a Novel Mutation in the Plasmodium falciparum Chloroquine Resistance Transporter With Decreased Piperaquine Sensitivity. J. Infect. Dis. 2017, 216, 468–476. [Google Scholar] [CrossRef]
  46. Sanchez, C.P.; Dave, A.; Stein, W.D.; Lanzer, M. Transporters as mediators of drug resistance in Plasmodium falciparum. Int. J. Parasitol. 2010, 40, 1109–1118. [Google Scholar] [CrossRef]
  47. Kavishe, R.A.; van den Heuvel, J.M.; van de Vegte-Bolmer, M.; Luty, A.J.; Russel, F.G.; Koenderink, J.B. Localization of the ATP-binding cassette (ABC) transport proteins PfMRP1, PfMRP2, and PfMDR5 at the Plasmodium falciparum plasma membrane. Malar. J. 2009, 8, 205. [Google Scholar] [CrossRef]
  48. Foote, S.J.; Kyle, D.E.; Martin, R.K.; Oduola, A.M.; Forsyth, K.; Kemp, D.J.; Cowman, A.F. Several alleles of the multidrug-resistance gene are closely linked to chloroquine resistance in Plasmodium falciparum. Nature 1990, 345, 255–258. [Google Scholar] [CrossRef]
  49. Sidhu, A.B.; Valderramos, S.G.; Fidock, D.A. pfmdr1 mutations contribute to quinine resistance and enhance mefloquine and artemisinin sensitivity in Plasmodium falciparum. Mol. Microbiol. 2005, 57, 913–926. [Google Scholar] [CrossRef]
  50. Pickard, A.L.; Wongsrichanalai, C.; Purfield, A.; Kamwendo, D.; Emery, K.; Zalewski, C.; Kawamoto, F.; Miller, R.S.; Meshnick, S.R. Resistance to antimalarials in Southeast Asia and genetic polymorphisms in pfmdr1. Antimicrob. Agents Chemother. 2003, 47, 2418–2423. Available online: http://www.ncbi.nlm.nih.gov/pubmed/12878499 (accessed on 10 July 2019). [CrossRef]
  51. Reed, M.B.; Saliba, K.J.; Caruana, S.R.; Kirk, K.; Cowman, A.F. Pgh1 modulates sensitivity and resistance to multiple antimalarials in Plasmodium falciparum. Nature 2000, 403, 906–909. [Google Scholar] [CrossRef]
  52. Sisowath, C.; Ferreira, P.E.; Bustamante, L.Y.; Dahlstrom, S.; Martensson, A.; Bjorkman, A.; Krishna, S.; Gil, J.P. The role of pfmdr1 in Plasmodium falciparum tolerance to artemether-lumefantrine in Africa. Trop. Med. Int. Health TM IH 2007, 12, 736–742. [Google Scholar] [CrossRef] [Green Version]
  53. Humphreys, G.S.; Merinopoulos, I.; Ahmed, J.; Whitty, C.J.; Mutabingwa, T.K.; Sutherland, C.J.; Hallett, R.L. Amodiaquine and artemether-lumefantrine select distinct alleles of the Plasmodium falciparum mdr1 gene in Tanzanian children treated for uncomplicated malaria. Antimicrob. Agents Chemother. 2007, 51, 991–997. [Google Scholar] [CrossRef]
  54. Njokah, M.J.; Kang’ethe, J.N.; Kinyua, J.; Kariuki, D.; Kimani, F.T. In vitro selection of Plasmodium falciparum Pfcrt and Pfmdr1 variants by artemisinin. Malar. J. 2016, 15, 381. [Google Scholar] [CrossRef] [Green Version]
  55. Raj, D.K.; Mu, J.; Jiang, H.; Kabat, J.; Singh, S.; Sullivan, M.; Fay, M.P.; McCutchan, T.F.; Su, X.Z. Disruption of a Plasmodium falciparum multidrug resistance-associated protein (PfMRP) alters its fitness and transport of antimalarial drugs and glutathione. J. Biol. Chem. 2009, 284, 7687–7696. [Google Scholar] [CrossRef]
  56. Gupta, B.; Xu, S.; Wang, Z.; Sun, L.; Miao, J.; Cui, L.; Yang, Z. Plasmodium falciparum multidrug resistance protein 1 (pfmrp1) gene and its association with in vitro drug susceptibility of parasite isolates from north-east Myanmar. J. Antimicrob. Chemother. 2014, 69, 2110–2117. [Google Scholar] [CrossRef] [Green Version]
  57. Dondorp, A.M.; Nosten, F.; Yi, P.; Das, D.; Phyo, A.P.; Tarning, J.; Lwin, K.M.; Ariey, F.; Hanpithakpong, W.; Lee, S.J.; et al. Artemisinin resistance in Plasmodium falciparum malaria. N. Engl. J. Med. 2009, 361, 455–467. [Google Scholar] [CrossRef]
  58. Witkowski, B.; Amaratunga, C.; Khim, N.; Sreng, S.; Chim, P.; Kim, S.; Lim, P.; Mao, S.; Sopha, C.; Sam, B.; et al. Novel phenotypic assays for the detection of artemisinin-resistant Plasmodium falciparum malaria in Cambodia: In-vitro and ex-vivo drug-response studies. Lancet Infect. Dis. 2013, 13, 1043–1049. [Google Scholar] [CrossRef]
  59. Taylor, S.M.; Parobek, C.M.; DeConti, D.K.; Kayentao, K.; Coulibaly, S.O.; Greenwood, B.M.; Tagbor, H.; Williams, J.; Bojang, K.; Njie, F.; et al. Absence of putative artemisinin resistance mutations among Plasmodium falciparum in Sub-Saharan Africa: A molecular epidemiologic study. J. Infect. Dis. 2015, 211, 680–688. [Google Scholar] [CrossRef] [PubMed]
  60. Tun, K.M.; Imwong, M.; Lwin, K.M.; Win, A.A.; Hlaing, T.M.; Hlaing, T.; Lin, K.; Kyaw, M.P.; Plewes, K.; Faiz, M.A.; et al. Spread of artemisinin-resistant Plasmodium falciparum in Myanmar: A cross-sectional survey of the K13 molecular marker. Lancet Infect. Dis. 2015, 15, 415–421. [Google Scholar] [CrossRef]
  61. Huang, F.; Takala-Harrison, S.; Jacob, C.G.; Liu, H.; Sun, X.; Yang, H.; Nyunt, M.M.; Adams, M.; Zhou, S.; Xia, Z.; et al. A Single Mutation in K13 Predominates in Southern China and Is Associated with Delayed Clearance of Plasmodium falciparum Following Artemisinin Treatment. J. Infect. Dis. 2015, 212, 1629–1635. [Google Scholar] [CrossRef]
  62. Nyunt, M.H.; Wang, B.; Aye, K.M.; Aye, K.H.; Han, J.H.; Lee, S.K.; Han, K.T.; Htut, Y.; Han, E.T. Molecular surveillance of artemisinin resistance falciparum malaria among migrant goldmine workers in Myanmar. Malar. J. 2017, 16, 97. [Google Scholar] [CrossRef]
  63. Wang, J.; Huang, Y.; Zhao, Y.; Ye, R.; Zhang, D.; Pan, W. Introduction of F446I mutation in the K13 propeller gene leads to increased ring survival rates in Plasmodium falciparum isolates. Malar. J. 2018, 17, 248. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Genes 10 00692 g001 550
Figure 1. Map of study areas and flow chart of resistance markers detection. (A) Three study townships are indicated as black dots in the maps. (B) The sequences of the samples in the final black coiled frame were analyzed for single nucleotide polymorphism (SNP) and haplotypes.
Figure 1. Map of study areas and flow chart of resistance markers detection. (A) Three study townships are indicated as black dots in the maps. (B) The sequences of the samples in the final black coiled frame were analyzed for single nucleotide polymorphism (SNP) and haplotypes.
Genes 10 00692 g001
Genes 10 00692 g002 550
Figure 2. Prevalence of 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) haplotypes in P. falciparum isolates from asymptomatic carriers. (A) Prevalence of pfdhfr haplotypes, mutations at codons 51, 59, 108, and 164. including wild-type NCSI, double mutant ICNI and NRNI, triple mutant IRNI and NRNL, and quadruple mutant IRNL; (B) Prevalence of pfdhps haplotypes, mutations at codons 436, 437, 540, and 581, including wild-type SGKA, single mutant AGKA, SAKA, SGEA, and SGKG, double mutant AGEA, AGNA, SAEA, SGEG, and SGEG; (C) Prevalence of pfcrt haplotypes, mutations at codons 74, 75, 76, 220, 271, 326, 356, and 371, including wild-type MNKAQNIR, single mutant MNKAQNII, sextuple mutant IETSENII, septuple mutant IETSENTI, and octuple mutant IETSESTI; (D) Prevalence of pfmdr1 haplotypes, mutations at codons 86, 130, 184, 185, and 1225, including wild-type NEYIV, single mutant NEFIV and NEYKV, double mutant YEYKV, NKYKV, NEFKV, and NEFNV, triple mutant NEFKL; (E) Prevalence of pfmrp1 haplotypes, mutations at codons 191, 325, 437, 572, 785, 876, 1007, 1339, and 1390, including wild-type HNSFHITAF, single mutant HNSFHVTAF, double mutant YNAFHITAF, quadruple mutant YSAFHVTAF, YNAFHVMAF, and YNAFHVTAI, quintuple mutant YNAFNVMAF, YNAFHVMAI, and YNAFHVTSI, sextuple mutant YNALNVMAF; (F) Prevalence of k13 haplotypes, mutations at codons 191, 446, 574 and N, NN, NNN insertion, including wild-type, N insertion, NN insertion, NNN insertion, NN insertion with 446I, NN insertion with 574L, and single mutant 191T.
Figure 2. Prevalence of 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) haplotypes in P. falciparum isolates from asymptomatic carriers. (A) Prevalence of pfdhfr haplotypes, mutations at codons 51, 59, 108, and 164. including wild-type NCSI, double mutant ICNI and NRNI, triple mutant IRNI and NRNL, and quadruple mutant IRNL; (B) Prevalence of pfdhps haplotypes, mutations at codons 436, 437, 540, and 581, including wild-type SGKA, single mutant AGKA, SAKA, SGEA, and SGKG, double mutant AGEA, AGNA, SAEA, SGEG, and SGEG; (C) Prevalence of pfcrt haplotypes, mutations at codons 74, 75, 76, 220, 271, 326, 356, and 371, including wild-type MNKAQNIR, single mutant MNKAQNII, sextuple mutant IETSENII, septuple mutant IETSENTI, and octuple mutant IETSESTI; (D) Prevalence of pfmdr1 haplotypes, mutations at codons 86, 130, 184, 185, and 1225, including wild-type NEYIV, single mutant NEFIV and NEYKV, double mutant YEYKV, NKYKV, NEFKV, and NEFNV, triple mutant NEFKL; (E) Prevalence of pfmrp1 haplotypes, mutations at codons 191, 325, 437, 572, 785, 876, 1007, 1339, and 1390, including wild-type HNSFHITAF, single mutant HNSFHVTAF, double mutant YNAFHITAF, quadruple mutant YSAFHVTAF, YNAFHVMAF, and YNAFHVTAI, quintuple mutant YNAFNVMAF, YNAFHVMAI, and YNAFHVTSI, sextuple mutant YNALNVMAF; (F) Prevalence of k13 haplotypes, mutations at codons 191, 446, 574 and N, NN, NNN insertion, including wild-type, N insertion, NN insertion, NNN insertion, NN insertion with 446I, NN insertion with 574L, and single mutant 191T.
Genes 10 00692 g002
Genes 10 00692 g003 550
Figure 3. Median-joining haplotype network of asymptomatic P. falciparum isolates harboring mutations in six drug resistance-associated genes. The haplotype network was constructed for asymptomatic P. falciparum isolates using 43 haplotypes obtained from amino acid changes observed in pfdhfr, pfdhps, pfcrt, pfmdr1, pfmrp1, and k13. The size of each circle shows the same haplotype prevalence, with a color corresponding to a population of origin: Laiza (red), Banmauk (blue), or Paletwa (yellow). The length of an edge is proportional to the number of variations between two haplotypes, and enlarged is the torso of the tree. K13 individual mutations are labeled by amino acid position.
Figure 3. Median-joining haplotype network of asymptomatic P. falciparum isolates harboring mutations in six drug resistance-associated genes. The haplotype network was constructed for asymptomatic P. falciparum isolates using 43 haplotypes obtained from amino acid changes observed in pfdhfr, pfdhps, pfcrt, pfmdr1, pfmrp1, and k13. The size of each circle shows the same haplotype prevalence, with a color corresponding to a population of origin: Laiza (red), Banmauk (blue), or Paletwa (yellow). The length of an edge is proportional to the number of variations between two haplotypes, and enlarged is the torso of the tree. K13 individual mutations are labeled by amino acid position.
Genes 10 00692 g003
Table
Table 1. Prevalence of dihydrofolate reductase (pfdhfr) mutations and haplotypes in asymptomatic P. falciparum isolates collected from three townships of Myanmar.
Table 1. Prevalence of dihydrofolate reductase (pfdhfr) mutations and haplotypes in asymptomatic P. falciparum isolates collected from three townships of Myanmar.
MutationsNumber of Isolate (%)p Value #
LaizaBanmaukPaletwaTotal
n = 12n = 42n = 26n = 80
N51I 11 (91.7)26 (61.9)11 (42.3)48 (60.0)0.012
C59R 12 (100.0)36 (85.7)24 (92.3)72 (90.0)0.467
S108N 12 (100.0)37 (88.1)25 (96.2)74 (92.5)0.425
I164L 12 (100.0)13 (31.0)7 (26.9)32 (40.0)0.000
HaplotypeCodon *
Wild-typeNCSI0 (0.0)5 (11.9)1 (3.8)6 (7.5)0.425
Double mutant (a)ICNI0 (0.0)1 (2.4)1 (3.8)2 (2.5)1.000
Double mutant (b)NRNI0 (0.0)8 (19.0)11 (42.3)19 (23.8)0.009
Triple mutant (a) IRNI0 (0.0)15 (35.7)6 (23.1)21 (26.3)0.034
Triple mutant (b) NRNL1 (8.3)3 (7.1)3 (11.5)7 (8.8)0.863
Quadruple mutant IRNL11 (91.7)10 (23.8)4 (15.4)25 (31.3)0.000
* Point mutations are shown in boldface. # Comparison between the three different sites was done by Fisher’s exact test.
Table
Table 2. Prevalence of dihydropteroate synthase (pfdhps) mutations and haplotypes in asymptomatic P. falciparum isolates collected from three townships of Myanmar.
Table 2. Prevalence of dihydropteroate synthase (pfdhps) mutations and haplotypes in asymptomatic P. falciparum isolates collected from three townships of Myanmar.
Mutation at CodonNumber of Isolate (%)p Value #
LaizaBanmaukPaletwaTotal
n = 12n = 43n = 27n = 82
S436A3 (25.0)9 (20.9)15 (55.6)27 (32.9)0.010
G437A0 (0.0)21 (48.8)0 (0.0)21 (25.6)0.000
K540E6 (50.0)13 (30.2)25 (92.6)44 (53.7)0.000
K540N6 (50.0)1 (2.3)0 (0.0)7 (8.5)0.000
A581G9 (75.0)11 (25.6)7 (25.9)27 (32.9)0.005
HaplotypeCodon *
Wild-typeSGKA0 (0.0)2 (4.7)0 (0.0)2 (2.4)0.650
Single mutant (a)AGKA0 (0.0)1 (2.3)0 (3.3)1 (1.2)1.000
Single mutant (b)SAKA0 (0.0)20 (46.5)0 (3.3)20 (24.4)0.000
Single mutant (c)SGEA0 (0.0)0 (0.0)5 (18.5)5 (6.1)0.007
Single mutant (d)SGKG0 (0.0)6 (14.0)2 (7.4)8 (9.8)0.405
Double mutant (a)AGEA2 (16.7)8 (18.6)15 (55.6)25 (30.5)0.003
Double mutant (b)AGNA1 (8.3)0 (0.0)0 (0.0)1 (1.2)0.146
Double mutant (c)SAEA0 (0.0)1 (2.3)0 (0.0)1 (1.2)1.000
Double mutant (d)SGEG4 (33.3)4 (9.3)5 (18.5)13 (15.9)0.094
Double mutant (e)SGNG5 (41.7)1 (2.3)0 (0.0)6 (7.3)0.000
* Point mutations are shown in boldface. # Comparison between the three different sites was done by Fisher’s exact test (p < 0.05).
Table
Table 3. Prevalence of chloroquine (CQ) resistance transporter (pfcrt) mutations and haplotypes in asymptomatic P. falciparum isolates collected from three townships of Myanmar.
Table 3. Prevalence of chloroquine (CQ) resistance transporter (pfcrt) mutations and haplotypes in asymptomatic P. falciparum isolates collected from three townships of Myanmar.
Mutation at CodonNumber of Isolate (%)p Value #
LaizaBanmaukPaletwaTotal
n = 12n = 39n = 24n = 75
M74I 12 (100.0)27 (69.2)23 (95.8)62 (82.7)0.005
N75E 12 (100.0)27 (69.2)23 (95.8)62 (82.7)0.005
K76T 12 (100.0)27 (69.2)23 (95.8)62 (82.7)0.005
A220S 12 (100.0)27 (69.2)23 (95.8)62 (82.7)0.005
Q271E 12 (100.0)27 (69.2)23 (95.8)62 (82.7)0.005
N326S 12 (100.0)6 (15.4)2 (8.3)20 (26.7)0.000
I356T 12 (100.0)27 (69.2)20 (83.3)59 (78.7)0.025
R371I 12 (100.0)28 (71.8)23 (95.8)63 (84.0)0.016
HaplotypeCodon *
Wild-typeMNKAQNIR0 (0.0)11 (28.2)1 (4.2)12 (16.0)0.011
Single mutantMNKAQNII0 (0.0)1 (2.6)0 (0.0)1 (1.3)1.000
Sextuple mutantIETSENII0 (0.0)0 (0.0)3 (12.5)3 (4.0)0.057
Septuple mutantIETSENTI0 (0.0)21 (53.8)18 (75.0)39 (52.0)0.000
Octuple mutantIETSESTI12 (100.0)6 (15.4)2 (8.3)20 (26.7)0.000
* Point mutations are shown in boldface with underline. # Comparison between the three different sites was done by Fisher’s exact test (p < 0.05).
Table
Table 4. Prevalence of multidrug resistance protein 1 (pfmdr1) mutations and haplotypes in asymptomatic P. falciparum isolates collected from three townships of Myanmar.
Table 4. Prevalence of multidrug resistance protein 1 (pfmdr1) mutations and haplotypes in asymptomatic P. falciparum isolates collected from three townships of Myanmar.
Mutation at CodonNumber of Isolate (%)p Value #
LaizaBanmaukPaletwaTotal
n = 12n = 41n = 27n = 80
N86Y 0 (0.0)0 (0.0)2 (7.4)2 (2.5)0.234
E130K 0 (0.0)0 (0.0)1 (3.7)1 (1.3)0.488
Y184F 7 (58.3)8 (19.5)4 (14.8)19 (23.8)0.014
I185K 4 (33.3)39 (95.1)27 (100.0)70 (87.5)0.000
I185N 0 (0.0)2 (4.9)0 (0.0)2 (2.5)0.650
V1225L 0 (0.0)1 (2.4)0 (0.0)1 (1.3)1.000
HaplotypeCodon *
Wild-typeNEYIV1 (8.3)0 (0.0)0 (0.0)1 (1.3)0.150
Single mutant (a)NEFIV7 (58.3)0 (0.0)0 (0.0)7 (8.8)0.000
Single mutant (b)NEYKV4 (33.3)33 (80.5)20 (74.1)57 (71.3)0.008
Double mutant (a)YEYKV0 (0.0)0 (0.0)2 (7.4)2 (2.5)0.234
Double mutant (b)NKYKV0 (0.0)0 (0.0)1 (3.7)1 (1.3)0.488
Double mutant (c)NEFKV0 (0.0)5 (12.2)4 (14.8)9 (11.3)0.554
Double mutant (d)NEFNV0 (0.0)2 (4.9)0 (0.0)2 (2.5)0.650
Triple mutantNEFKL0 (0.0)1 (2.4)0 (0.0)1 (1.3)1.000
* Point mutations are shown in boldface. # Comparison between the three different sites was done by Fisher’s exact test (p < 0.05).
Table
Table 5. Prevalence of multidrug resistance-associated protein 1 (pfmrp1) mutations and haplotypes in asymptomatic P. falciparum isolates collected from three townships of Myanmar.
Table 5. Prevalence of multidrug resistance-associated protein 1 (pfmrp1) mutations and haplotypes in asymptomatic P. falciparum isolates collected from three townships of Myanmar.
Mutation at CodonNumber of Isolate (%)p Value #
LaizaBanmaukPaletwaTotal
n = 12n = 43n = 27n = 82
H191Y 9 (75.0)21 (48.8)24 (88.9)54 (65.9)0.002
N325S 1 (8.3)1 (2.3)0 (0.0)2 (2.4)0.379
S437A 9 (75.0)21 (48.8)24 (88.9)54 (65.9)0.002
F572L 1 (8.3)1 (2.3)0 (0.0)2 (2.4)0.379
H785N 7 (58.3)11 (25.6)0 (0.0)18 (22.0)0.000
I876V 8 (66.7)20 (46.5)21 (77.8)49 (59.8)0.029
T1007M 7 (58.3)13 (30.2)21 (77.8)41 (50.0)0.000
A1339S 0 (0.0)2 (4.7)0 (0.0)2 (2.4)0.650
F1390I 1 (8.3)6 (14.0)0 (0.0)7 (8.5)0.135
HaplotypeCodon *
Wild-typeHNSFHITAF3 (25.0)20 (46.5)3 (11.1)26 (31.7)0.006
Single mutantHNSFHVTAF0 (0.0)2 (4.7)0 (0.0)2 (2.4)0.650
Double mutantYNAFHITAF0 (0.0)2 (4.7)3 (11.1)5 (6.1)0.433
Quadruple mutant (a)YSAFHVTAF1 (8.3)1 (2.3)0 (0.0)2 (2.4)0.379
Quadruple mutant (b)YNAFHVMAF0 (0.0)1 (2.3)21 (77.8)22 (26.8)0.000
Quadruple mutant (c)YNAFHVTAI1 (8.3)3 (7.0)0 (0.0)4 (4.9)0.357
Quintuple mutant (a)YNAFNVMAF6 (50.0)10 (23.3)0 (0.0)16 (19.5)0.000
Quintuple mutant (b)YNAFHVMAI0 (0.0)1 (2.3)0 (0.0)1 (1.2)1.000
Quintuple mutant (c)YNAFHVTSI0 (0.0)2 (4.7)0 (0.0)2 (2.4)0.650
Sextuple mutantYNALNVMAF1 (8.3)1 (2.3)0 (0.0)2 (2.4)0.379
* Point mutations are shown in boldface. # Comparison between the three different sites was done by Fisher’s exact test (p < 0.05).
Table
Table 6. Prevalence of Kelch protein 13 (k13) mutations and haplotypes in asymptomatic P. falciparum isolates collected in three areas.
Table 6. Prevalence of Kelch protein 13 (k13) mutations and haplotypes in asymptomatic P. falciparum isolates collected in three areas.
Mutation at CodonNumber of Isolate (%)p Value e
LaizaBanmaukPaletwaTotal
n = 12n = 23n = 22n = 57
N insertion a0 (0.0)1 (4.3)0 (0.0)1 (1.8)1.000
NN insertion b12 (100.0)7 (30.4)7 (31.8)26 (45.6)0.000
NNN insertion c0 (0.0)1 (4.3)3 (13.6)4 (7.0)0.423
K191T0 (0.0)4 (17.4)0 (0.0)4 (7.0)0.067
F446I8 (66.7)0 (0.0)0 (0.0)8 (14.0)0.000
P574L1 (8.3)0 (0.0)0 (0.0)1 (1.8)0.211
Haplotyped
Wild-type0 (0.0)10 (43.5)12 (54.5)22 (38.6)0.003
N insertion0 (0.0)1 (4.3)0 (0.0)1 (1.8)1.000
NN insertion3 (25.0)7 (30.4)7 (31.8)17 (29.8)1.000
NNN insertion0 (0.0)1 (4.3)3 (13.6)4 (7.0)0.423
NN insertion F446I8 (66.7)0 (0.0)0 (0.0)8 (14.0)0.000
NN insertion P574L1 (8.3)0 (0.0)0 (0.0)1 (1.8)0.211
K191T0 (0.0)4 (17.4)0 (0.0)4 (7.0)0.067
a N insertion between amino acids 136 and 137. b NN insertion between amino acids 136 and 137. c NNN insertion between amino acids 136 and 137. d Point mutations are shown in boldface. e Comparison between the three different sites was done by Fisher’s exact test (p < 0.05).

Share and Cite

MDPI and ACS Style

Zhao, Y.; Liu, Z.; Soe, M.T.; Wang, L.; Soe, T.N.; Wei, H.; Than, A.; Aung, P.L.; Li, Y.; Zhang, X.; et al. Genetic Variations Associated with Drug Resistance Markers in Asymptomatic Plasmodium falciparum Infections in Myanmar. Genes 2019, 10, 692. https://doi.org/10.3390/genes10090692

AMA Style

Zhao Y, Liu Z, Soe MT, Wang L, Soe TN, Wei H, Than A, Aung PL, Li Y, Zhang X, et al. Genetic Variations Associated with Drug Resistance Markers in Asymptomatic Plasmodium falciparum Infections in Myanmar. Genes. 2019; 10(9):692. https://doi.org/10.3390/genes10090692

Chicago/Turabian Style

Zhao, Yan, Ziling Liu, Myat Thu Soe, Lin Wang, Than Naing Soe, Huanping Wei, Aye Than, Pyae Linn Aung, Yuling Li, Xuexing Zhang, and et al. 2019. "Genetic Variations Associated with Drug Resistance Markers in Asymptomatic Plasmodium falciparum Infections in Myanmar" Genes 10, no. 9: 692. https://doi.org/10.3390/genes10090692

APA Style

Zhao, Y., Liu, Z., Soe, M. T., Wang, L., Soe, T. N., Wei, H., Than, A., Aung, P. L., Li, Y., Zhang, X., Hu, Y., Wei, H., Zhang, Y., Burgess, J., Siddiqui, F. A., Menezes, L., Wang, Q., Kyaw, M. P., Cao, Y., & Cui, L. (2019). Genetic Variations Associated with Drug Resistance Markers in Asymptomatic Plasmodium falciparum Infections in Myanmar. Genes, 10(9), 692. https://doi.org/10.3390/genes10090692

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop