Skip to main content

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

Abstract

Background

Confirmation of artemisinin-delayed parasite clearance in Plasmodium falciparum along the Thai-Myanmar border has inspired a global response to contain and monitor drug resistance to avert the disastrous consequences of a potential spread to Africa. However, resistance data from Myanmar are sparse, particularly from high-risk areas where limited health services and decades of displacement create conditions for resistance to spread. Subclinical infections may represent an important reservoir for resistance genes that confer a fitness disadvantage relative to wild-type alleles. This study estimates the prevalence of resistance genotypes in three previously unstudied remote populations in Myanmar and tests the a priori hypothesis that resistance gene prevalence would be higher among isolates collected from subclinical infections than isolates collected from febrile clinical patients. A systematic review of resistance studies is provided for context.

Methods

Community health workers in Karen and Kachin States and an area spanning the Indo-Myanmar border collected dried blood spots from 988 febrile clinical patients and 4,591 villagers with subclinical infection participating in routine prevalence surveys. Samples positive for P. falciparum 18 s ribosomal RNA by real-time PCR were genotyped for P. falciparum multidrug resistance protein (pfmdr1) copy number and the pfcrt K76T polymorphism using multiplex real-time PCR.

Results

Pfmdr1 copy number increase and the pfcrt K76 polymorphism were determined for 173 and 269 isolates, respectively. Mean pfmdr1 copy number was 1.2 (range: 0.7 to 3.7). Pfmdr1 copy number increase was present in 17.5%, 9.6% and 11.1% of isolates from Karen and Kachin States and the Indo-Myanmar border, respectively. Pfmdr1 amplification was more prevalent in subclinical isolates (20.3%) than clinical isolates (6.4%, odds ratio 3.7, 95% confidence interval 1.1 - 12.5). Pfcrt K76T prevalence ranged from 90-100%.

Conclusions

Community health workers can contribute to molecular surveillance of drug resistance in remote areas of Myanmar. Marginal and displaced populations under-represented among previous resistance investigations can and should be included in resistance surveillance efforts, particularly once genetic markers of artemisinin-delayed parasite clearance are identified. Subclinical infections may contribute to the epidemiology of drug resistance, but determination of gene amplification from desiccated filter samples requires further validation when DNA concentration is low.

Background

Genetically determined artemisinin-delayed parasite clearance or tolerance, first documented on the Thai-Cambodia border, has emerged on the border of Thailand and Myanmar[1, 2]. Recently published estimates of malaria mortality trends suggest that the previous spread of chloroquine (CQ) and sulphadoxine-pyrimethamine resistance from Southeast Asia to Africa[3, 4] contributed to the large increase in malaria mortality from 1980–2004[5], and spread of delayed parasite clearance to Africa would represent a global health catastrophe[6].

In recognition of Myanmar’s central location between Southeast Asia and Africa, the WHO outlined a country-specific strategy for Myanmar Artemisinin Resistance Containment (MARC)[7]. The MARC surveillance strategy calls for therapeutic efficacy studies and day-3, parasite-positivity monitoring of artemisinin combination therapy (ACT) in over 20 locations. However, a strategy is lacking to assess drug resistance in remote populations of Myanmar that government and international NGO services have difficulty reaching.

Although molecular markers of artemisinin resistance or delayed parasite clearance have yet to be identified, tracking markers of resistance to partner drugs provides a valuable tool to inform coformulation policy and monitor progress of the Global Plan for Artemisinin Containment (GPARC) and MARC. Clinical and parasitological failure after treatment with an ACT is partially determined by the local efficacy of non-artemisinin partner drugs[8]. Furthermore, the decreased parasiticidal effect of artemisinin along the Myanmar-Thai border places greater reliance on partner drugs. Tracking Plasmodium falciparum multidrug resistance protein (pfmdr1) gene copy number (CN) has become an important surveillance tool, particularly in populations receiving artesunate-mefloquine. Pfmdr1 CN is associated with delayed response to ACT, including artesunate-mefloquine[9] and artemether–lumefantrine[10, 11], as well as resistance to multiple monotherapy including mefloquine[1215].

Molecular strategies are particularly valuable for surveillance in remote, displaced and conflict-affected populations facing security and logistical constraints that make conventional in vivo or in vitro resistance studies impractical[16]. Village health workers (VHWs) trained by local community-based organizations play a key role delivering malaria control services in hard-to-access areas of Myanmar[1719] emerging from decades of conflict, but community based organizations and VHWs have not contributed to past resistance surveillance efforts. Recent studies successfully estimated pfmdr1 (CN) using blood samples collected on filter paper[20], and this simplified protocol makes it possible to extend the quantitative assessment of gene copy number to remote settings lacking the capacity for storage and transport of whole or fractionated blood products.

One question of potential importance in resistance containment is the relative contribution of asymptomatic persons to the reservoir of genetic resistance. Asymptomatic infections provide a parasite reservoir that contributes to malaria transmission even in areas of low or unstable transmission intensity[2126], but uncertainty exists regarding the role of asymptomatic infections in the transmission of genetic changes conferring drug-resistance. Genetic markers of drug resistance in P. falciparum, including pfcrt K76T haplotypes and pfmdr1 amplification,[27] are associated with decreased parasite fitness and impaired within-host growth[28]. Less fit parasites may be more likely to produce low parasitemia infections that remain asymptomatic. Persistent carriage of resistant parasites in asymptomatic populations is unlikely in situations where multiclonal infections are common, given that more fit parasites with a single copy of pfmdr1 are likely to out-compete less fit parasites with multiple copies of pfmdr1. However, asymptomatic carriage could be favoured in relatively low transmission settings such as those included in this study where host immunity may be insufficient to clear infection and mono-infection with a single resistant clone may be common, particularly if de-amplification of pfmdr1 CN is rare. In addition, compensatory mutations that decrease the fitness cost of pfmdr1 amplification may allow drug-resistant clones to persist during multiclonal infections. The persistence of drug-resistant infections[29], including multicopy pfmdr1 infections[30], may also be favoured by higher gametocyte carriage that appears to increase their transmission potential relative to drug-sensitive, wild-type infections.

This study sought to assess the feasibility of estimating the prevalence of elevated pfmdr1 CN and the pfcrt K76T allele in three previously unstudied remote or conflict-affected populations living along Myanmar’s borders with India, China and Thailand. This study investigated the contribution of subclinical infections to the epidemiology of genetic resistance by testing the a priori hypothesis that resistance gene prevalence would be higher among isolates collected during active population screening compared to isolates collected from febrile clinical patients. The findings have been situated in context by conducting a systematic review of in vivo, in vitro, and molecular resistance studies in Myanmar and border regions of neighbouring countries. Resistance study locations were mapped relative to areas of recent civil conflict to test the second a priori hypothesis that politically unstable areas would be under-represented among resistance studies in Myanmar.

Methods

Study sites and participants

Between February 2009 and January 2010, three community based organizations established 13 surveillance sites in three remote regions of Myanmar: Karen State, Kachin state, and an area spanning the border of Chin State, Myanmar and Mizoram State, India (Figure1). Populations were selected based on history of malaria burden, feasibility of sample collection, and past exclusion from other surveillance programmes due to inaccessibility or security concerns. Samples were collected by community health workers (CHWs) and VHWs who extend the reach of community-based malaria control services to unstable areas, as previously described[18]. The red-shaded areas of the map in Figure1 indicate the locations of villages displaced since 2008 in Karen and Kachin States, where armed conflict was ongoing or imminent, respectively. Sites in border areas of Chin State were located in areas where exposure to human rights violations were widespread[31].

Figure 1
figure 1

Anti-malarial, drug-resistance studies in Myanmar and neighbouring countries, by region and year. Blue squares represent the locations of in vivo, in vitro and molecular studies conducted between 1996 and 2009; locations of the current study (2008–2009) appear as red circles (n = 13 sites). Studies reporting prevalence of pfmdr1 copy number amplification or pfcrt K76T haplotypes are numbered on the map and on the x-axis of the bar figures, in order of data collection year. Previous molecular studies from Myanmar appear in insets; molecular studies from neighbouring countries appear on the country-wide map. Locations of villages forcibly displaced between 2008–2011 are indicated by red shading. See manuscript text for additional information on displaced villages, and Additional file1 for a complete list of molecular studies and abstracted prevalence estimates.

Oral, informed consent was obtained from all participants. Ethical approval for the study was obtained from relevant authorities at each community based organization; Johns Hopkins School of Public Health Committee on Human Subjects approved secondary analysis of the data.

Samples were collected during routine malaria programme activities modified to assess the prevalence of non-P. falciparum species using PCR and validate the use of a P. falciparum-specific, histidine-rich protein 2 (HRP-II) based rapid diagnostic test (RDT, Paracheck®, Orchid Biomedical Systems, Goa, India) in the setting of clinical diagnosis and population-based screening. CHWs collected blood spots from consecutive clinic patients with self-reported fever and from a sub-sample of villagers participating in routine bi-annual community malaria prevalence surveys[18, 19, 32]. All participants were tested for P. falciparum using the RDT and had a blood sample collected on Whatman 903 Protein Saver® cards. The sampling protocol for the programme evaluation was designed to yield approximately n = 540 Pf- positive samples: n = 80 from febrile clinical patients and n = 100 from villagers participating in active screenings in each of three regions. Based on a real-time polymerase chain reaction (RT-PCR) budget sufficient to test 2,100 samples, the target sample number for each region was set at 700 (500 screening and 200 clinical). Because historical programme data indicated that P. falciparum population prevalence in each region was consistently below 15-20%, the number of screening participants was increased to yield approximately 100 P. falciparum -positive isolates. False-positive RDT results due to persistent HRP-II antigenaemia were expected to approximately offset false-negative RDT results due to low sensitivity in the setting of low parasitaemia infection.

With an expected RDT-positivity rate of approximately 40%, sample collection was anticipated to yield ~80 P. falciparum -positive isolates from 200 febrile clinical patients in each region. When preliminary results suggested a lower RDT positivity rate (18-33%), the period of sample collection was extended from two to six months. PCR was performed on all RDT-positive samples and a region-specific fraction of RDT-negative samples from screening participants (22-34%) and clinic patients (17%-70%) sufficient to yield approximately 500 + 200 = 700 samples for RT-PCR from each region.

DNA extraction and Plasmodium falciparum detection by real-time PCR

Following an existing protocol[33], DNA was extracted from dried blood samples using a commercial, 96-well, DNA extraction kit from Promega (Madison, Wisconsin, USA). For each sample, DNA from three 5-mm punches (equivalent to approximately 25 μL whole blood) was extracted and concentrated four-fold by glycogen-acetate precipitation. Plasmodium falciparum DNA was identified using a validated TaqMan multiplex real-time PCR assay for P. falciparum 18 s ribosomal DNA[33], a molecular method for P. falciparum detection that has sensitivity and specificity comparable to nested PCR protocols[34, 35]. Isolates were considered P. falciparum- positive if 18 s PCR end-cycle fluorescence exceeded 300 RFU. All PCR reactions were run in 384-well plates on the Bio-Rad CFX384 Real-Time System (Bio-Rad, Hercules, CA, USA). Primers and probes for TaqMan PCR were synthesized by Integrated DNA Technologies (Coraville, Iowa, USA).

pfmdr1 copy number determination by multiplex real-time PCR

pfmdr1 copy number was determined using a modification of previous assays[13, 36]. Primers and a FAM-labelled probe specific for pfmdr1, plus primers and a Texas Red-labelled probe specific for P. falciparum β-tubulin, were amplified as a multiplex. Each reaction included β-tubulin primers (100nM), β-tubulin probes (100nM), pfmdr1 primers (300nM), pfmdr1 probes (200nM), Bio-Rad iQ Multiplex Powermix (1x), and 5 μL DNA, in a total reaction volume of 10 ul.

Amplification curves were analysed and cycle threshold (Ct) values determined using Bio-Rad CFX Manager software. pfmdr1 copy number was calculated using the predictive efficiency method[36], which corrects the ΔΔ Ct method for differences in amplification efficiency between the reference gene and gene of interest. Genomic DNA from P. falciparum strains 3D7, D10, Dd2 were included as internal controls on each 384-well plate (MR4, Virginia, USA). These strains are known to have one, one, and three to four pfmdr1 copies, respectively. Genomic DNA from 7C424, a strain previously determined to have two copies of pfmdr1, was also used to validate the assay. Mean copy number estimates were 1.15 for D10 (SD = 0.20, 32 repeats), 2.31 for 7C424 (SD = 0.42, 35 repeats), and 2.98 for Dd2 (SD = 0.39, 35 repeats).

All reactions were run in triplicate and individual replicates were rejected if they did not display exponential kinetics. The precision of pfmdr1 CN estimates was calculated for a given isolate as the ratio of the range (maximum–minimum) divided by the mean of isolate-specific replicates (n = 2 or 3). Copy number estimates were repeated if the difference between replicates was greater than 50% of the mean estimate, or if the average cycle threshold number for Pf β-tubulin was greater than 34. CN estimates were considered invalid if the variance remained greater than 50% of the mean estimate or Ct remained greater than 34 after the third repeat.

In the course of the study, a higher proportion of subclinical screening isolates than clinical isolates produced invalid CN estimates. A post-hoc analysis explored the reasons for this discrepancy and assessed the technical feasibility of using filter paper blood samples to estimate pfmdr1 gene copy number in isolates from a predominantly asymptomatic population with relatively low peripheral blood parasite concentrations. Recent evidence suggests that the volume of blood available from filter-paper samples may be insufficient to achieve highly accurate PCR results, particularly when parasitaemia level is low[26]. Current methods to assess gene copy number may be particularly sensitive to DNA concentration, because copy number is a function (ratio) of two quantitative amplification thresholds. Past studies assessed pfmdr1 CN using filter paper-based samples collected from febrile clinical patients[20]; this study explored the feasibility of extending these methods to predominantly asymptomatic infections among active screening participants, and whether the variance in CN estimates was associated with parasite DNA concentration. The precision of pfmdr1 CN estimates was plotted against cycle threshold (Ct) values, a proxy for parasitaemia level. Ct values are inversely proportional to the log DNA concentration in a given isolate, and a three-cycle increase in the Ct corresponds to approximately a 10-fold decrease in relative concentration of P. falciparum DNA.

Detection of pfcrt polymorphisms by multiplex real-time PCR

Plasmodium falciparum-positive samples were analysed for wild-type pfcrt and pfcrt alleles carrying the K76T polymorphism, using validated primers and probes for the CVIET and CVMNK haplotypes[37]. Each PCR reaction included forward and reverse primers for the pfcrt gene (200nM), a Texas Red-labelled probe specific for the wild-type allele (200nM), a FAM-labelled probe specific for the K76T allele (200nM), Bio-Rad iQ Multiplex Powermix (1x), and 5 μL DNA, in a total reaction volume of 10 μl. Average fluoresence of the last five PCR cycles (“end-cycle fluoresence”) was used to determine the presence or absence of the wild-type and K76T polymorphisms. Genomic DNA from the CQ-sensitive D10 and CQ-resistant Dd2 strains were included on each 384-well plate as references. All samples were run in duplicate and criteria for identifying CQ-sensitive parasites were chosen for high specificity. Samples were identified as carrying the wild-type CQ-sensitive allele if both replicates had end-cycle fluorescence greater than 5% of the range between negative and positive controls.

Statistical analysis

Data was entered into an ACCESS database using range and consistency checks. Statistical analyses were performed using Stata 12.1. Standard errors were corrected for clustering of the data by village cluster (n = 13).

Systematic review of Plasmodium falciparum anti-malarial resistance studies in Burma/Myanmar and its border regions

PubMed and Google Scholar were searched for clinical, in vitro, and molecular studies evaluating resistance of Plasmodium falciparum to CQ, MQ and ACTs in Myanmar and neighbouring countries. Search terms included “falciparum”, “resistance”, “efficacy”, “Burma”, “Myanmar”, “Thailand”, “China”, “India”, “Bangladesh”, “pfmdr1”, “pfcrt”, “chloroquine”, “mefloquine”, “lumefantrine”, “artemisinin combination therapy”, “anti-malarial”; and minor variations on these terms. This search generated 860 publications. Reference lists of studies identified in the primary search were searched, as was data from public sources, including WHO Southeast Asia Region and country-specific publications of the WHO and respective Ministries of Health. The review included reports from Bangladesh, India and China if P. falciparum isolates were collected from areas bordering Myanmar, and restricted the review to reports published in English after 1995. Only observational and clinical intervention studies and related review articles were included; studies that reported only basic science or methods-oriented findings were excluded. Eighty-two manuscripts met these inclusion criteria. Reports were reviewed for clinical efficacy (failure rates, recrudescence/re-infection, parasite clearance time), in vitro susceptibility data (i e, IC50s) and molecular markers associated with resistance (pfmdr1 SNPs, pfmdr1 amplification, and the K76T pfcrt mutation) to CQ and MQ monotherapy, and any ACT. Figure1 shows locations of all studies meeting criteria for this review (blue squares); molecular studies are numbered for reference to the adjacent bar graphs and are listed in Additional file1. Resistance study locations were determined from published manuscripts and by contacting study authors. Unstable areas of Karen and Kachin States were identified based on evidence of forced displacement of villages since 2008 provided by the Thai-Burma Border Consortium[38] and a community-based organization providing humanitarian assistance in Kachin State (Figure1).

Results

Active screening of 4,591 out of 14,982 people living in 13 village-clusters identified 157 isolates positive for P. falciparum 18 s DNA, and sequential testing of 988 febrile clinical patients identified 133 P. falciparum-positive isolates (Table1).

Table 1 Study population characteristics

pfmdr1 copy number

Of 290 P. falciparum- positive isolates, 173 (59.6%) were successfully genotyped for pfmdr1 copy number (Table2). Thirty-nine percent (18/46) of samples from Kachin State, 62% (115/186) from Chin State, and 80% (40/50) from Karen State yielded acceptable copy number estimates. Copy number ranged from 0.72 to 3.70. 11.6% (19/173, 95%CI = 2.7-20.4) of samples had >1.5 copies and 4.6% (8/173) had >2 copies. The proportion of isolates with copy number >1.5 (95% CI) was 9.6% in Kachin State, 11.1% in Chin State, and 17.5% in Karen State (Table2 and Figure2). Isolates from population-screening participants were more likely than those from febrile clinical patients to have elevated pfmdr1 copy number (OR 3.7, 95% CI 1.1-12.5). There was no difference in pfmdr1 amplification between participants who did and did not report taking anti-malarial treatment in the previous eight weeks.

Table 2 Pfmdr1 copy number (CN) among subclinical and clinical isolates from three regions of Myanmar 2008–2009
Figure 2
figure 2

Scatterplots and boxplots of pfmdr1 copy number estimate precision and Plasmodium falciparum beta-tubulin PCR cycle thresholds (Ct) for febrile clinical patients (grey circles and boxplots) and active screening participants (empty circles and boxplots) A: The scatterplots demonstrate precision of pfmdr1 copy number (CN) estimate (y-axis) diminishes as DNA concentration (x-axis) decreases among active screening participants (bottom panel) and febrile clinical patients (top panel). The mean cycle threshold and mean precision of pfmdr1 CN is indicated by the black point on each scatterplot. A three-cycle increase in the Ct corresponds to a 10-fold decrease in relative concentration of Pf DNA. Pfmdr1 precision is calculated as the range/mean for each isolate, with higher values indicating lower precision. The standard deviation of pfmdr1 CN estimates was 0.200 for clinic patients and .516 for screening participants (not shown) B: Boxplots display the distribution of DNA concentration (top panel) and pfmdr1 CN precision. Boxes represent inter-quartile ranges (IQR); whiskers represent the value of [upper/lower quartile +/− (IQR*1.5)]. Screening participants had lower DNA concentration and less precise estimation of pfmdr1 CN.

Pf DNA concentration and precision of pfmdr1 copy number determination

Of all CN estimates calculated, including repeat assays, 83 of 157 asymptomatic screening isolates (55%) and 23 of 133 clinic isolates (17%) had range:mean ratios greater than 0.5 and/or mean Pf b-tubulin Ct values greater than 34, indicating a Pf DNA concentration insufficient to estimate pfmdr1 CN. Precision of pfmdr1 copy number decreased with decreasing P. falciparum DNA concentration (Figure2). The mean P. falciparum DNA concentration was lower for screening samples than for clinical samples and the precision of copy number estimates was lower overall for screening samples (Figure2).

pfcrt wild-type and K76T alleles

Multiplex PCR detected the wild-type pfcrt allele in 11 of 290 isolates (Table3). Eight of 186 isolates from surveillance sites in Chin State carried the sensitive allele, as did three of 46 from sites in Kachin State; the sensitive allele was not identified in any of the 57 samples from Karen State. Twenty isolates failed to amplify either the sensitive or the K76T allele. Prevalence of the wild-type pfcrt allele was similar in isolates from active screening participants and febrile clinical patients, and did not vary by report of anti-malarial treatment in the previous eight weeks.

Table 3 Pfcrt genotypes among subclinical and clinical isolates from three regions of Myanmar 2008-2009

Systematic review

See Figure3 and Additional file1 for results of the brief systematic review of in vivo, in vitro and molecular resistance studies. Half of the 82 resistance studies meeting inclusion criteria were conducted in western Thailand, including all but two studies of pfmdr1 CN. One study of pfmdr1 CN from Myanmar was identified, conducted in Laiza in 2009 (Figure1B). Locations of in vivo, in vitro and molecular studies did not overlap with areas of instability, indicated forcibly displaced villages (red shading) on the map in Figure1.

Figure 3
figure 3

Systematic review of in vivo , in vitro , and molecular resistance studies in Myanmar and neighbouring countries, 1996–2009 [[1, 7, 9, 10, 13, 39, 40, 6082]].

Discussion

This study demonstrates the feasibility of incorporating a network of health workers affiliated with community-based organizations into strategies for molecular surveillance of malaria resistance in remote and unstable areas of Myanmar, including priority Tier 1 and Tier 2 areas at high risk for the spread of delayed parasite clearance. This surveillance network can accelerate the assessment of artemisinin resistance and delayed parasite clearance once genetic markers are identified. Myanmar’s ecologic, ethnic and political diversity, coupled with patchwork access to malaria control services and quality anti-malarials creates a mosaic of selective drug pressure, and this study provides a preliminary glimpse of patterns emerging over time and place in remote areas of the country, including areas directly across the border from western Thailand where delayed parasite clearance was recently confirmed. The prevalence of pfmdr1 amplification was relatively low (12-19%) across the three border regions included in the current study; and only a single isolate had a CN of three or more. Results from Karen State contrast with recent studies documenting elevated CN in at least 40% of isolates[9] with a mean CN of 2.9 (SE 1.4)[39] among refugee, migrant and cross-border populations in Thailand (Figure1). Findings from Kachin State contrast with a recent study in Laiza that found no pfmdr1 amplification in 171 P. falciparum isolates[40, 41], but is consistent with a study conducted in 2004 in nearby counties of Yunnan province[42]. The systematic review failed to identify previous studies of pfmdr1 CN along the borders with India and Bangladesh to compare to the 11% prevalence estimated for that region.

Multiple factors likely contribute to the observed variation in pfmdr1 amplification documented within and between populations living in Myanmar and across its borders in neighbouring countries. Little is known about host response in populations in Myanmar[24], but host immunity is likely to be relatively more robust given higher transmission intensity compared to western Thailand and Yunnan, where migrants from Myanmar account for the majority of infections[19, 4346]. Selective drug pressure also varies substantially due to differences in access to quality diagnosis and treatment. Populations included in the present study live in areas that government and international NGO health services largely fail to reach, and community based organizations are unique providers of malaria control interventions. The low prevalence of pfmdr1 amplification in Karen study sites may be due in part to these malaria programmes[17] that for 10 years have provided directly observed mefloquine-artesunate that may modify selective pressure for pfmdr1 amplification[47].

The pfcrt K76T mutation is a highly predictive marker of CQ resistance[48]. Replacement of K76T mutants with wild-type parasites following withdrawal of drug pressure has forecast the return of CQ clinical efficacy[49, 50] and some authors have proposed tracking this mutation to identify populations in which CQ could be used as partner drugs in ACT[51]. The very low prevalence of wild-type pfcrt in the present study is consistent with data from nearby areas and suggests ongoing selective pressure for CQ resistance. Although decades have passed since CQ was officially recommended by Myanmar or any of its neighbours for the treatment of P. falciparum infection, CQ remains the treatment of choice for P. vivax and presumptive malaria. The use of RDTs capable of identifying both P. falciparum and P. vivax species may modestly diminish this selective pressure, but tracking the prevalence of wild-type pfcrt is likely to remain a low priority for resistance surveillance for the next several years. However, pfcrt could become an important genetic marker if the recommended treatment of P. vivax changes from CQ to ACT in response to the recent emergence of CQ-resistant P. vivax in western Thailand[52, 53].

Consistent with the a priori hypothesis, the prevalence of elevated pfmdr1 CN was higher among isolates obtained from active screening participants than from isolates obtained from febrile clinical patients. Resistance surveillance currently relies on samples obtained from clinically symptomatic malaria patients; containment efforts may need to more aggressively target subclinical infections that could serve as a reservoir for the spread of drug resistance if future studies validate an association between subclinical infection and resistance alleles. However, there are several reasons to approach this finding with caution. First, to the authors’ knowledge this is the first study to evaluate pfmdr1 CN in a predominantly asymptomatic population. Zhong et al.[54] found no difference between symptomatic and asymptomatic study volunteers in Kenya in the frequency of pfmdr1 point mutations, or of the K76T mutation of pfcrt. Second, consistent with Zhong et al. the current study found no difference in K76T mutation prevalence between screening and clinical participants, though prevalence of wild-type pfcrt was low in each region and results may be confounded by ongoing CQ drug pressure as described above. Third, a robust empirically derived conceptual model is lacking to adequately capture the complex interactions between transmission intensity, host immunity, drug pressure, fitness cost, compensatory mutations, clonality of infection, de-amplification[27, 55] and other factors expected to influence the relative prevalence of resistant and wild-type parasites in clinical and asymptomatic infections. Clonality of infection is likely to play an important role, as described above, but the prevalence of monoclonal infection in Myanmar is poorly quantified. A single study from central Myanmar documented mono-infection in 32% of clinical isolates[56], but the relevance of these data to asymptomatic populations in border regions is unclear. The final reason to interpret with caution the association between subclinical infection and genetic resistance is the high variance of pfmdr1 CN estimates found among screening participants that led to the exclusion of 55% of samples. Future studies are needed to validate the accuracy of pfmdr1 CN estimates based on filter-paper blood samples collected from subclinical populations. Despite these caveats, findings presented here highlight the need to conduct additional studies on the contribution of subclinical infection to the epidemiology of drug-resistant malaria.

Although mosquitoes and migrating humans frequently carry parasites across international boundaries, it is not clear whether findings from neighbouring countries should be extrapolated to populations living in Myanmar. Individuals seeking care from border clinics by definition have access to diagnostic and treatment protocols and other resources of neighbouring countries that may differ substantially from those available in remote and conflict-affected areas of Myanmar. Additional studies are necessary to determine the appropriate geographic scale for monitoring the spread of drug resistance.

Systematic review of in vivo, in vitro and molecular resistance studies from Myanmar and its border regions identified relatively few data from populations living inside the country, particularly in Tier 1 and 2 areas of highest priority for resistance containment. The paucity of available information is consistent with output from the under-developed health research capacity in the country as a whole: Myanmar ranks 218th out of 224 countries in number of publications in medicine per capita (0.4 per 100,000 people)[57]. Armed conflict and large-scale population displacement are established causes of disruptions in health services and disease surveillance[58], and can result in “stability bias”, defined as the “systematic under-sampling of populations and health threats in contexts of conflict and instability”[59]. The lack of overlap between previous resistance study locations and displaced populations in Myanmar documented in the present study suggests that stability-bias may contribute to within-country distribution of evidence available for resistance surveillance. The analysis presented here highlights areas with documented population displacement, but does not capture all conflict and human rights violations experienced by other communities. For example, approximately 92% of households in Chin State may experience forced labour[31] and other human rights violations that have been associated with prevalent malaria infection[32]. The instability of sites participating in this study was tragically validated by violent events since data collection completed in 2009: one Karen site and all Kachin sites were displaced due to military attacks. Some health workers have continued to implement disease control interventions in relocated areas; and they remain willing to participate in resistance surveillance activities. One community-based organization is establishing a site to monitor parasite clearance time.

There were several important limitations to this study. Twenty P. falciparum isolates did not amplify either wild type or K76T alleles and it is not possible to exclude the presence of other rare resistance alleles such as SVMNK that were not assessed. None of these twenty isolates produced valid copy number estimates, and they do not modify our primary conclusions that pfmdr1 amplification may be more common among subclinical infections, and that that wild type pfcrt remains rare in these three regions of Myanmar.

The precision of pfmdr1 CN estimates determined from subclinical infections was low for 44% of isolates, and 16% of isolates had DNA concentrations below the lower limit of reliable copy number determination, even after four-fold DNA concentration with glycogen-acetate precipitation. As noted above, future studies are necessary to validate the use of filter-paper samples to estimate gene CN from low-parasitaemia isolates such as those collected from subclinical infections. This study was not designed to elucidate the mechanisms or relative contributions of factors influencing the prevalence of resistant parasites in subclinical populations. For example, clonality of infection was not determined and parasite DNA from prior infections was unavailable to distinguish re-infection from recrudescence.

The cross-sectional design limits causal inference, and findings related to measures of association should be considered preliminary, as noted above. Study participants lived in areas where quality malaria control services were ongoing for three to 8 years, and findings may not apply to intervention-naïve populations. Samples were collected in 2008–2009 and may not reflect the prevalence of resistance markers in 2012. For example, areas of western Thailand north of Maesot, immediately across the border from one of the present study sites, experienced a sharp rise in the proportion of slow-clearing parasites between 2008–10[1, 2].

The number of P. falciparum- positive isolates (290) was lower than anticipated by the study design (540) despite screening 4,591 villagers and 988 febrile clinical patients. Logistical constraints delayed data collection in Karen and Kachin areas until the lower-transmission, dry-season months of February and March, with extended clinical sampling through June. The study was conducted in areas of active malaria control, and the success of these programmes likely contributed to the low prevalence of P. falciparum and low RDT positivity rates. The smaller than anticipated number of P. falciparum positive samples limited the precision of estimates of genotype prevalence in subclinical and clinical infections, and poor statistical power precluded comparisons across geographic regions. Nevertheless, the number of isolates available for genotyping compares favourably to studies conducted in more stable areas, and data from individuals with subclinical P. falciparum infections is unique among published studies from this region.

Conclusions

Networks of community-based health organizations can and should contribute to molecular surveillance of anti-malarial drug resistance in remote and unstable areas of Myanmar, including priority Tier 1 and 2 areas at greatest risk for spread of delayed parasite clearance. Desiccated filter paper facilitates accurate estimation of pfmdr1 CN among febrile clinical isolates and extends the potential geographic range of genetic resistance surveillance, but further study is necessary to validate the use of filter paper in the setting of subclinical infections with low parasite DNA concentrations. The increased risk of pfmdr1 amplification observed among active screening participants hints at a plausible and potentially important contribution of subclinical infection to resistance transmission, but further study is necessary to elucidate the epidemiology of drug resistance among predominantly asymptomatic populations. The consistently low prevalence of wild-type pfcrt alleles, documented in this and other studies, is consistent with persistent CQ drug pressure that is unlikely to diminish until treatment protocols specify CQ-free regimens for P. vivax infections. A systematic review of resistance studies revealed a dearth of data from Myanmar, particularly from areas recently experiencing armed conflict and forced displacement. Marginal populations, such as those participating in this study, should be represented in future resistance surveillance efforts, particularly once genetic markers of delayed parasite clearance are identified. Failure to include these at-risk areas may compromize global efforts to monitor and contain the spread of resistance.

Abbreviations

ACPR:

Acceptable Clinical and Parasitological Result

ACT:

Artemisinin combination therapy

AL:

Artemether-lumefantrine

ART:

Artemisinin (eg: artesunate, artemether)

CN:

Copy number

CQ:

Chloroquine

DHP:

Dihydroartemisinin-piperaquine

MAS:

Mefloquine artesunate

MQ:

Mefloquine

pfcrt:

Plasmodium falciparum CQ resistance transporter

pfmdr1:

Plasmodium falciparum multidrug resistance protein (pfmdr1).

References

  1. Phyo AP, Nkhoma S, Stepniewska K, Ashley EA, Nair S, McGready R, Ler Moo C, Al-Saai S, Dondorp AM, Lwin KM, Singhasivanon P, Day NP, White NJ, Anderson TJ, Nosten F: Emergence of artemisinin-resistant malaria on the western border of Thailand: a longitudinal study. Lancet. 2012, 379: 1960-1966. 10.1016/S0140-6736(12)60484-X.

    Article  PubMed Central  PubMed  Google Scholar 

  2. Cheeseman IH, Miller BA, Nair S, Nkhoma S, Tan A, Tan JC, Al Saai S, Phyo AP, Moo CL, Lwin KM, McGready R, Ashley E, Imwong M, Stepniewska K, Yi P, Dondorp AM, Mayxay M, Newton PN, White NJ, Nosten F, Ferdig MT, Anderson TJ: A major genome region underlying artemisinin resistance in malaria. Science. 2012, 336: 79-82. 10.1126/science.1215966.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  3. Dondorp AM, Newton PN, Mayxay M, Van Damme W, Smithuis FM, Yeung S, Petit A, Lynam AJ, Johnson A, Hien TT, McGready R, Farrar JJ, Looareesuwan S, Day NP, Green MD, White NJ: Fake antimalarials in Southeast Asia are a major impediment to malaria control: multinational cross-sectional survey on the prevalence of fake antimalarials. Trop Med Int Health. 2004, 9: 1241-1246. 10.1111/j.1365-3156.2004.01342.x.

    Article  CAS  PubMed  Google Scholar 

  4. Takahashi N, Tanabe K, Tsukahara T, Dzodzomenyo M, Dysoley L, Khamlome B, Sattabongkot J, Nakamura M, Sakurai M, Kobayashi J, Kaneko A, Endo H, Hombhanje F, Tsuboi T, Mita T: Large-scale survey for novel genotypes of Plasmodium falciparum chloroquine-resistance gene pfcrt. Malar J. 2012, 11: 92-10.1186/1475-2875-11-92.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  5. Murray CJ, Rosenfeld LC, Lim SS, Andrews KG, Foreman KJ, Haring D, Fullman N, Naghavi M, Lozano R, Lopez AD: Global malaria mortality between 1980 and 2010: a systematic analysis. Lancet. 2012, 379: 413-431. 10.1016/S0140-6736(12)60034-8.

    Article  PubMed  Google Scholar 

  6. White NJ, Nosten F, Looareesuwan S, Watkins WM, Marsh K, Snow RW, Kokwaro G, Ouma J, Hien TT, Molyneux ME, Taylor TE, Newbold CI, Ruebush TK, Danis M, Greenwood BM, Anderson RM, Olliaro P: Averting a malaria disaster. Lancet. 1999, 353: 1965-1967. 10.1016/S0140-6736(98)07367-X.

    Article  CAS  PubMed  Google Scholar 

  7. WHO and Department of Health, Myanmar: Strategic Framework for Artemisinin Resistance Containment in Myanmar (MARC) 2011–2015. 2011, Myanmar: Department of Health, Myanmar and World Health Organization

    Google Scholar 

  8. WHO: Global report on antimalarial drug efficacy and drug resistance: 2000–2010. 2010, Geneva: World Health Organization

    Google Scholar 

  9. Muhamad P, Phompradit P, Sornjai W, Maensathian T, Chaijaroenkul W, Rueangweerayut R, Na-Bangchang K: Polymorphisms of molecular markers of antimalarial drug resistance and relationship with artesunate-mefloquine combination therapy in patients with uncomplicated Plasmodium falciparum malaria in Thailand. Am J Trop Med Hyg. 2011, 85: 568-572. 10.4269/ajtmh.2011.11-0194.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  10. Price RN, Uhlemann AC, van Vugt M, Brockman A, Hutagalung R, Nair S, Nash D, Singhasivanon P, Anderson TJ, Krishna S, White NJ, Nosten F: Molecular and pharmacological determinants of the therapeutic response to artemether-lumefantrine in multidrug-resistant Plasmodium falciparum malaria. Clin Infect Dis. 2006, 42: 1570-1577. 10.1086/503423.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  11. Pongtavornpinyo W, Hastings IM, Dondorp A, White LJ, Maude RJ, Saralamba S, Day NP, White NJ, Boni MF: Probability of emergence of antimalarial resistance in different stages of the parasite life cycle. Evol App. 2009, 2: 52-61. 10.1111/j.1752-4571.2008.00067.x.

    Article  Google Scholar 

  12. Pickard AL, Wongsrichanalai C, Purfield A, Kamwendo D, Emery K, Zalewski C, Kawamoto F, Miller RS, Meshnick SR: Resistance to antimalarials in Southeast Asia and genetic polymorphisms in pfmdr1. Antimicrob Agents Chemother. 2003, 47: 2418-2423. 10.1128/AAC.47.8.2418-2423.2003.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  13. Price RN, Uhlemann AC, Brockman A, McGready R, Ashley E, Phaipun L, Patel R, Laing K, Looareesuwan S, White NJ, Nosten F, Krishna S: Mefloquine resistance in Plasmodium falciparum and increased pfmdr1 gene copy number. Lancet. 2004, 364: 438-447. 10.1016/S0140-6736(04)16767-6.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  14. Picot S, Olliaro P, de Monbrison F, Bienvenu AL, Price RN, Ringwald P: A systematic review and meta-analysis of evidence for correlation between molecular markers of parasite resistance and treatment outcome in falciparum malaria. Malar J. 2009, 8: 89-10.1186/1475-2875-8-89.

    Article  PubMed Central  PubMed  Google Scholar 

  15. Uhlemann AC, Krishna S: Antimalarial multi-drug resistance in Asia: mechanisms and assessment. Curr Top Microbiol Immunol. 2005, 295: 39-53. 10.1007/3-540-29088-5_2.

    CAS  PubMed  Google Scholar 

  16. Djimde AA, Dolo A, Ouattara A, Diakite S, Plowe CV, Doumbo OK: Molecular diagnosis of resistance to antimalarial drugs during epidemics and in war zones. J Infect Dis. 2004, 190: 853-855. 10.1086/422758.

    Article  PubMed  Google Scholar 

  17. Richards AK, Banek K, Mullany LC, Lee CI, Smith L, Oo EK, Lee TJ: Cross-border malaria control for internally displaced persons: observational results from a pilot programme in eastern Burma/Myanmar. Trop Med Int Health. 2009, 14: 512-521. 10.1111/j.1365-3156.2009.02258.x.

    Article  PubMed  Google Scholar 

  18. Lee CI, Smith LS, Shwe Oo EK, Scharschmidt BC, Whichard E, Kler T, Lee TJ, Richards AK: Internally displaced human resources for health: villager health worker partnerships to scale up a malaria control programme in active conflict areas of eastern Burma. Glob Public Health. 2009, 4: 229-241. 10.1080/17441690802676360.

    Article  CAS  PubMed  Google Scholar 

  19. Richards AK, Smith L, Mullany LC, Lee CI, Whichard E, Banek K, Mahn M, Shwe Oo EK, Lee TJ: Prevalence of Plasmodium falciparum in active conflict areas of eastern Burma: a summary of cross-sectional data. Confl Health. 2007, 1: 9-10.1186/1752-1505-1-9.

    Article  PubMed Central  PubMed  Google Scholar 

  20. Shah NK, Alker AP, Sem R, Susanti AI, Muth S, Maguire JD, Duong S, Ariey F, Meshnick SR, Wongsrichanalai C: Molecular surveillance for multidrug-resistant Plasmodium falciparum, Cambodia. Emerg Infect Dis. 2008, 14: 1637-1640. 10.3201/eid1410.080080.

    Article  PubMed Central  PubMed  Google Scholar 

  21. da Silva-Nunes M, Moreno M, Conn JE, Gamboa D, Abeles S, Vinetz JM, Ferreira MU: Amazonian malaria: asymptomatic human reservoirs, diagnostic challenges, environmentally driven changes in mosquito vector populations, and the mandate for sustainable control strategies. Acta Trop. 2012, 121: 281-291. 10.1016/j.actatropica.2011.10.001.

    Article  PubMed Central  PubMed  Google Scholar 

  22. Incardona S, Vong S, Chiv L, Lim P, Nhem S, Sem R, Khim N, Doung S, Mercerau-Puijalon O, Fandeur T: Large-scale malaria survey in Cambodia: Novel insights on species distribution and risk factors. Malar J. 2007, 6: 37-10.1186/1475-2875-6-37.

    Article  PubMed Central  PubMed  Google Scholar 

  23. Macauley C: Aggressive active case detection: a malaria control strategy based on the Brazilian model. Soc Sci Med. 2005, 60: 563-573. 10.1016/j.socscimed.2004.05.025.

    Article  PubMed  Google Scholar 

  24. Soe S, Khin Saw A, Htay A, Nay W, Tin A, Than S, Roussilhon C, Perignon JL, Druilhe P: Premunition against Plasmodium falciparum in a malaria hyperendemic village in Myanmar. Trans R Soc Trop Med Hyg. 2001, 95: 81-84. 10.1016/S0035-9203(01)90342-6.

    Article  Google Scholar 

  25. Alves FP, Gil LH, Marrelli MT, Ribolla PE, Camargo EP, Da Silva LH: Asymptomatic carriers of Plasmodium spp. as infection source for malaria vector mosquitoes in the Brazilian Amazon. J Med Entomol. 2005, 42: 777-779. 10.1603/0022-2585(2005)042[0777:ACOPSA]2.0.CO;2.

    Article  PubMed  Google Scholar 

  26. Proux S, Suwanarusk R, Barends M, Zwang J, Price RN, Leimanis M, Kiricharoen L, Laochan N, Russell B, Nosten F, Snounou G: Considerations on the use of nucleic acid-based amplification for malaria parasite detection. Malar J. 2011, 10: 323-10.1186/1475-2875-10-323.

    Article  PubMed Central  PubMed  Google Scholar 

  27. Preechapornkul P, Imwong M, Chotivanich K, Pongtavornpinyo W, Dondorp AM, Day NP, White NJ, Pukrittayakamee S: Plasmodium falciparum pfmdr1 amplification, mefloquine resistance, and parasite fitness. Antimicrob Agents Chemother. 2009, 53: 1509-1515. 10.1128/AAC.00241-08.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  28. Babiker HA, Hastings IM, Swedberg G: Impaired fitness of drug-resistant malaria parasites: evidence and implication on drug-deployment policies. Expert Rev Anti Infect Ther. 2009, 7: 581-593. 10.1586/eri.09.29.

    Article  CAS  PubMed  Google Scholar 

  29. Barnes KI, White NJ: Population biology and antimalarial resistance: the transmission of antimalarial drug resistance in Plasmodium falciparum. Acta Trop. 2005, 94: 230-240. 10.1016/j.actatropica.2005.04.014.

    Article  CAS  PubMed  Google Scholar 

  30. Price RN, Nosten F, Luxemburger C, ter Kuile FO, Paiphun L, Chongsuphajaisiddhi T, White NJ: Effects of artemisinin derivatives on malaria transmissibility. Lancet. 1996, 347: 1654-1658. 10.1016/S0140-6736(96)91488-9.

    Article  CAS  PubMed  Google Scholar 

  31. Sollom R, Richards AK, Parmar P, Mullany LC, Lian SB, Iacopino V, Beyrer C: Health and human rights in Chin State Western Burma: a population-based assessment using multistaged household cluster sampling. PLoS Med. 2011, 8: e1001007-10.1371/journal.pmed.1001007.

    Article  PubMed Central  PubMed  Google Scholar 

  32. Mullany LC, Richards AK, Lee CI, Suwanvanichkij V, Maung C, Mahn M, Beyrer C, Lee TJ: Population-based survey methods to quantify associations between human rights violations and health outcomes among internally displaced persons in eastern Burma. J Epidemiol Community Health. 2007, 61: 908-914. 10.1136/jech.2006.055087.

    Article  PubMed Central  PubMed  Google Scholar 

  33. Schachterle SE, Mtove G, Levens JP, Clemens EG, Shi L, Raj A, Munoz B, Reller ME, West S, Dumler JS, Sullivan D: Prevalence and density-related concordance of three diagnostic tests for malaria in a region of Tanzania with hypoendemic malaria. J Clin Microbiol. 2011, 49: 3885-3891. 10.1128/JCM.01157-11.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  34. Rosanas-Urgell A, Mueller D, Betuela I, Barnadas C, Iga J, Zimmerman PA, del Portillo HA, Siba P, Mueller I, Felger I: Comparison of diagnostic methods for the detection and quantification of the four sympatric Plasmodium species in field samples from Papua New Guinea. Malar J. 2010, 9: 361-10.1186/1475-2875-9-361.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  35. Hwang SY, Kim SH, Lee GY, Hang VT, Moon CS, Shin JH, Koo WL, Kim SY, Park HJ, Park HO, Kho WG: A novel real-time PCR assay for the detection of Plasmodium falciparum and Plasmodium vivax malaria in low parasitized individuals. Acta Trop. 2011, 120: 40-45. 10.1016/j.actatropica.2011.05.006.

    Article  CAS  PubMed  Google Scholar 

  36. Alker AP, Lim P, Sem R, Shah NK, Yi P, Bouth DM, Tsuyuoka R, Maguire JD, Fandeur T, Ariey F, Wongsrichanalai C, Meshnick SR: Pfmdr1 and in vivo resistance to artesunate-mefloquine in falciparum malaria on the Cambodian-Thai border. Am J Trop Med Hyg. 2007, 76: 641-647.

    CAS  PubMed  Google Scholar 

  37. Sutherland CJ, Haustein T, Gadalla N, Armstrong M, Doherty JF, Chiodini PL: Chloroquine-resistant Plasmodium falciparum infections among UK travellers returning with malaria after chloroquine prophylaxis. J Antimicrob Chemother. 2007, 59: 1197-1199. 10.1093/jac/dkm104.

    Article  CAS  PubMed  Google Scholar 

  38. Thai Burma Border Consortium: Displaced villages in Eastern Burma. 2008, –2011

    Google Scholar 

  39. Poyomtip T, Suwandittakul N, Sitthichot N, Khositnithikul R, Tan-ariya P, Mungthin M: Polymorphisms of the pfmdr1 but not the pfnhe-1 gene is associated with in vitro quinine sensitivity in Thai isolates of Plasmodium falciparum. Malar J. 2012, 11: 7-10.1186/1475-2875-11-7.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  40. Yang Z, Li C, Miao M, Zhang Z, Sun X, Meng H: Multidrug-resistant genotypes of Plasmodium falciparum, Myanmar. Emerg Inf Dis. 2011, 17: 498-501. 10.3201/eid1703.100870.

    Article  Google Scholar 

  41. Meng H, Zhang R, Yang H, Fan Q, Su X, Miao J, Cui L, Yang Z: in vitro sensitivity of Plasmodium falciparum clinical isolates from the China-Myanmar border area to quinine and association with polymorphism in the Na+/H + exchanger. Antimicrob Agents Chemother. 2010, 54: 4306-4313. 10.1128/AAC.00321-10.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  42. Yang Z, Zhang Z, Sun X, Wan W, Cui L, Zhang X, Zhong D, Yan G, Cui L: Molecular analysis of chloroquine resistance in Plasmodium falciparum in Yunnan Province, China. Trop Med Int Health. 2007, 12: 1051-1060. 10.1111/j.1365-3156.2007.01882.x.

    Article  CAS  PubMed  Google Scholar 

  43. Zhou G, Sirichaisinthop J, Sattabongkot J, Jones J, Bjornstad ON, Yan G, Cui L: Spatio-temporal distribution of Plasmodium falciparum and P. vivax malaria in Thailand. Am J Trop Med Hyg. 2005, 72: 256-262.

    PubMed  Google Scholar 

  44. Wiwanitkit V: High prevalence of malaria in Myanmar migrant workers in a rural district near the Thailand-Myanmar border. Scand J Infect Dis. 2002, 34: 236-237. 10.1080/00365540110077272.

    Article  PubMed  Google Scholar 

  45. Carrara VI, Sirilak S, Thonglairuam J, Rojanawatsirivet C, Proux S, Gilbos V, Brockman A, Ashley EA, McGready R, Krudsood S, Leemingsawat S, Looareesuwan S, Singhasivanon P, White N, Nosten F: Deployment of early diagnosis and mefloquine- artesunate treatment of falciparum malaria in Thailand: The Tak malaria initiative. PLoS Med. 2006, 3: e183-10.1371/journal.pmed.0030183.

    Article  PubMed Central  PubMed  Google Scholar 

  46. Zhou SS, Wang Y, Tang LH: [Malaria situation in the People's Republic of China in 2006] (in Chinese). Chinese journal of parasitology & parasitic diseases. 2007, 25: 439-441.

    Google Scholar 

  47. Uhlemann AC, McGready R, Ashley EA, Brockman A, Singhasivanon P, Krishna S, White NJ, Nosten F, Price RN: Intrahost selection of Plasmodium falciparum pfmdr1 alleles after antimalarial treatment on the northwestern border of Thailand. J Infect Dis. 2007, 195: 134-141. 10.1086/509809.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  48. Djimde A, Doumbo OK, Cortese JF, Kayentao K, Doumbo S, Diourte Y, Coulibaly D, Dicko A, Su XZ, Nomura T, Fidock DA, Wellems TE, Plowe CV: A molecular marker for chloroquine-resistant falciparum malaria. N Engl J Med. 2001, 344: 257-263. 10.1056/NEJM200101253440403.

    Article  CAS  PubMed  Google Scholar 

  49. Kublin JG, Cortese JF, Njunju EM, Mukadam RA, Wirima JJ, Kazembe PN, Djimde AA, Kouriba B, Taylor TE, Plowe CV: Reemergence of chloroquine-sensitive Plasmodium falciparum malaria after cessation of chloroquine use in Malawi. J Infect Dis. 2003, 187: 1870-1875. 10.1086/375419.

    Article  PubMed  Google Scholar 

  50. Laufer MK, Thesing PC, Eddington ND, Masonga R, Dzinjalamala FK, Takala SL, Taylor TE, Plowe CV: Return of chloroquine antimalarial efficacy in Malawi. N Engl J Med. 2006, 355: 1959-1966. 10.1056/NEJMoa062032.

    Article  CAS  PubMed  Google Scholar 

  51. Plowe CV, Roper C, Barnwell JW, Happi CT, Joshi HH, Mbacham W, Meshnick SR, Mugittu K, Naidoo I, Price RN, Shafer RW, Sibley CH, Sutherland CJ, Zimmerman PA, Rosenthal PJ: World Antimalarial Resistance Network (WARN) III: molecular markers for drug resistant malaria. Malar J. 2007, 6: 121-10.1186/1475-2875-6-121.

    Article  PubMed Central  PubMed  Google Scholar 

  52. Phyo AP, Lwin KM, Price RN, Ashley EA, Russell B, Sriprawat K, Lindegardh N, Singhasivanon P, White NJ, Nosten F: Dihydroartemisinin-piperaquine versus chloroquine in the treatment of Plasmodium vivax malaria in Thailand: a randomized controlled trial. Clin Infect Dis. 2011, 53: 977-984. 10.1093/cid/cir631.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  53. Rijken MJ, Boel ME, Russell B, Imwong M, Leimanis ML, Phyo AP, Muehlenbachs A, Lindegardh N, McGready R, Renia L, Snounou G, Singhasivanon P, Nosten F: Chloroquine resistant vivax malaria in a pregnant woman on the western border of Thailand. Malar J. 2011, 10: 113-10.1186/1475-2875-10-113.

    Article  PubMed Central  PubMed  Google Scholar 

  54. Zhong D, Afrane Y, Githeko A, Cui L, Menge DM, Yan G: Molecular epidemiology of drug-resistant malaria in western Kenya highlands. BMC Infect Dis. 2008, 8: 105-10.1186/1471-2334-8-105.

    Article  PubMed Central  PubMed  Google Scholar 

  55. Chen N, Chavchich M, Peters JM, Kyle DE, Gatton ML, Cheng Q: Deamplification of pfmdr1-containing amplicon on chromosome 5 in Plasmodium falciparum is associated with reduced resistance to artelinic acid in vitro. Antimicrob Agents Chemother. 2010, 54: 3395-3401. 10.1128/AAC.01421-09.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  56. Kang JM, Moon SU, Kim JY, Cho SH, Lin K, Sohn WM, Kim TS, Na BK: Genetic polymorphism of merozoite surface protein-1 and merozoite surface protein-2 in Plasmodium falciparum field isolates from Myanmar. Malar J. 2010, 9: 131-10.1186/1475-2875-9-131.

    Article  PubMed Central  PubMed  Google Scholar 

  57. McKee M, Stuckler D, Basu S: Where there is no health research: what can be done to fill the global gaps in health research?. PLoS Med. 2012, 9: e1001209-10.1371/journal.pmed.1001209.

    Article  PubMed Central  PubMed  Google Scholar 

  58. Gayer M, Legros D, Formenty P, Connolly MA: Conflict and emerging infectious diseases. Emerg Infect Dis. 2007, 13: 1625-1631. 10.3201/eid1311.061093.

    Article  PubMed Central  PubMed  Google Scholar 

  59. Beyrer C, Terzian A, Lowther S, Zambrano JA, Galai N, Melchior MK: Civil Conflict and Health Information: The Democratic Republic of Congo. Public Health and Human Rights: Evidence-based Approaches. Edited by: Chris B, Pizer HF. 2007, Baltimore: Johns Hopkins University Press, 268-285.

    Google Scholar 

  60. Smithuis F, Kyaw MK, Phe O, Aye KZ, Htet L, Barends M, Lindegardh N, Singtoroj T, Ashley E, Lwin S: Efficacy and effectiveness of dihydroartemisinin-piperaquine versus artesunate-mefloquine in falciparum malaria: an open-label randomised comparison. Lancet. 2006, 367: 2075-2085. 10.1016/S0140-6736(06)68931-9.

    Article  CAS  PubMed  Google Scholar 

  61. Smithuis F, Kyaw MK, Phe O, Win T, Aung PP, Oo AP, Naing AL, Nyo MY, Myint NZ, Imwong M, Stepniewska K, White NJ: Effectiveness of five artemisinin combination regimens with or without primaquine in uncomplicated falciparum malaria: an open-label randomised trial. Lancet Infect Dis. 2010, 10: 673-681. 10.1016/S1473-3099(10)70187-0.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  62. Smithuis F, Shahmanesh M, Kyaw MK, Savran O, Lwin S, White NJ: Comparison of chloroquine, sulfadoxine/pyrimethamine, mefloquine and mefloquine-artesunate for the treatment of falciparum malaria in Kachin State, North Myanmar. Trop Med Int Health. 2004, 9: 1184-1190. 10.1111/j.1365-3156.2004.01323.x.

    Article  CAS  PubMed  Google Scholar 

  63. Smithuis F, van der Broek I, Katterman N, Kyaw MK, Brockman A, Lwin S, White NJ: Optimising operational use of artesunate-mefloquine: a randomised comparison of four treatment regimens. Trans R Soc Trop Med Hyg. 2004, 98: 182-192. 10.1016/S0035-9203(03)00035-X.

    Article  CAS  PubMed  Google Scholar 

  64. Wongsrichanalai C, Lin K, Pang LW, Faiz MA, Noedl H, Wimonwattrawatee T, Laoboonchai A, Kawamoto F: in vitro susceptibility of Plasmodium falciparum isolates from Myanmar to antimalarial drugs. Am J Trop Med Hyg. 2001, 65: 450-455.

    CAS  PubMed  Google Scholar 

  65. Carrara VI, Zwang J, Ashley EA, Price RN, Stepniewska K, Barends M, Brockman A, Anderson T, McGready R, Phaiphun L, Proux S, van Vugt M, Hutagalung R, Lwin KM, Phyo AP, Preechapornkul P, Imwong M, Pukrittayakamee S, Singhasivanon P, White NJ, Nosten F: Changes in the treatment responses to artesunate-mefloquine on the northwestern border of Thailand during 13 years of continuous deployment. PLoS One. 2009, 4: e4551-10.1371/journal.pone.0004551.

    Article  PubMed Central  PubMed  Google Scholar 

  66. Nosten F, ter Kuile F, Chongsuphajaisiddhi T, Luxemburger C, Webster HK, Edstein M, Phaipun L, Thew KL, White NJ: Mefloquine-resistant falciparum malaria on the Thai-Burmese border. Lancet. 1991, 337: 1140-1143. 10.1016/0140-6736(91)92798-7.

    Article  CAS  PubMed  Google Scholar 

  67. Nosten F, van Vugt M, Price R, Luxemburger C, Thway KL, Brockman A, McGready R, ter Kuile F, Looareesuwan S, White NJ: Effects of artesunate-mefloquine combination on incidence of Plasmodium falciparum malaria and mefloquine resistance in western Thailand: a prospective study. Lancet. 2000, 356: 297-302. 10.1016/S0140-6736(00)02505-8.

    Article  CAS  PubMed  Google Scholar 

  68. Chaijaroenkul W, Wisedpanichkij R, Na-Bangchang K: Monitoring of in vitro susceptibilities and molecular markers of resistance of Plasmodium falciparum isolates from Thai-Myanmar border to chloroquine, quinine, mefloquine and artesunate. Acta Trop. 2010, 113: 190-194. 10.1016/j.actatropica.2009.10.016.

    Article  CAS  PubMed  Google Scholar 

  69. Imwong M, Dondorp AM, Nosten F, Yi P, Mungthin M, Hanchana S, Das D, Phyo AP, Lwin KM, Pukrittayakamee S, Lee SJ, Saisung S, Koecharoen K, Nguon C, Day NP, Socheat D, White NJ: Exploring the contribution of candidate genes to artemisinin resistance in Plasmodium falciparum. Antimicrob Agents Chemother. 2010, 54: 2886-2892. 10.1128/AAC.00032-10.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  70. Phompradit P, Wisedpanichkij R, Muhamad P, Chaijaroenkul W, Na-Bangchang K: Molecular analysis of pfatp6 and pfmdr1 polymorphisms and their association with in vitro sensitivity in Plasmodium falciparum isolates from the Thai-Myanmar border. Acta Trop. 2011, 120: 130-135. 10.1016/j.actatropica.2011.07.003.

    Article  CAS  PubMed  Google Scholar 

  71. Congpuong K, Na Bangchang K, Mungthin M, Bualombai P, Wernsdorfer WH: Molecular epidemiology of drug resistance markers of Plasmodium falciparum malaria in Thailand. Trop Med Int Health. 2005, 10: 717-722. 10.1111/j.1365-3156.2005.01450.x.

    Article  CAS  PubMed  Google Scholar 

  72. Yang H, Liu D, Yang Y, Fan B, Yang P, Li X, Li C, Dong Y, Yang C: Changes in susceptibility of Plasmodium falciparum to artesunate in vitro in Yunnan Province, China. Trans R Soc Trop Med Hyg. 2003, 97: 226-228. 10.1016/S0035-9203(03)90127-1.

    Article  CAS  PubMed  Google Scholar 

  73. Yang H, Yang Y, Yang P, Li X, Gao B, Zhang Z, Yang Z, Cui L: Monitoring Plasmodium falciparum chloroquine resistance in Yunnan Province, China, 1981–2006. Acta Trop. 2008, 108: 44-49. 10.1016/j.actatropica.2008.08.010.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  74. Zhang GQ, Guan YY, Zheng B, Wu S, Tang LH: Molecular assessment of Plasmodium falciparum resistance to antimalarial drugs in China. Trop Med Int Health. 2009, 14: 1266-1271. 10.1111/j.1365-3156.2009.02342.x.

    Article  CAS  PubMed  Google Scholar 

  75. Shah NK, Dhillon GP, Dash AP, Arora U, Meshnick SR, Valecha N: Antimalarial drug resistance of Plasmodium falciparum in India: changes over time and space. Lancet Infect Dis. 2011, 11: 57-64. 10.1016/S1473-3099(10)70214-0.

    Article  PubMed Central  PubMed  Google Scholar 

  76. Valecha N, Srivastava P, Mohanty SS, Mittra P, Sharma SK, Tyagi PK, Pradhan K, Dev V, Singh R, Dash AP, Sharma YD: Therapeutic efficacy of artemether-lumefantrine in uncomplicated falciparum malaria in India. Malar J. 2009, 8: 107-10.1186/1475-2875-8-107.

    Article  PubMed Central  PubMed  Google Scholar 

  77. Mohapatra PK, Namchoom NS, Prakash A, Bhattacharya DR, Goswami BK, Mahanta J: Therapeutic efficacy of anti-malarials in Plasmodium falciparum malaria in an Indo-Myanmar border area of Arunachal Pradesh. Indian J Med Res. 2003, 118: 71-76.

    CAS  PubMed  Google Scholar 

  78. Awasthi G, Prasad GB, Das A: Population genetic analyses of Plasmodium falciparum chloroquine receptor transporter gene haplotypes reveal the evolutionary history of chloroquine-resistant malaria in India. Int J Parasitol. 2011, 41: 705-709. 10.1016/j.ijpara.2011.03.002.

    Article  PubMed  Google Scholar 

  79. van den Broek IV, Maung UA, Peters A, Liem L, Kamal M, Rahman M, Rahman MR, Bangali AM, Das S, Barends M, Faiz AM: Efficacy of chloroquine + sulfadoxine–pyrimethamine, mefloquine + artesunate and artemether + lumefantrine combination therapies to treat Plasmodium falciparum malaria in the Chittagong Hill Tracts, Bangladesh. Trans R Soc Trop Med Hyg. 2005, 99: 727-735. 10.1016/j.trstmh.2005.02.007.

    Article  CAS  PubMed  Google Scholar 

  80. Noedl H, Socheat D, Satimai W: Artemisinin-resistant malaria in Asia. N Engl J Med. 2009, 361: 540-541. 10.1056/NEJMc0900231.

    Article  CAS  PubMed  Google Scholar 

  81. Noedl H, Faiz MA, Yunus EB, Rahman MR, Hossain MA, Samad R, Miller RS, Pang LW, Wongsrichanalai C: Drug-resistant malaria in Bangladesh: an in vitro assessment. Am J Trop Med Hyg. 2003, 68: 140-142.

    CAS  PubMed  Google Scholar 

  82. van den Broek IV, van der Wardt S, Talukder L, Chakma S, Brockman A, Nair S, Anderson TC: Drug resistance in Plasmodium falciparum from the Chittagong Hill Tracts, Bangladesh. Trop Med Int Health. 2004, 9: 680-687. 10.1111/j.1365-3156.2004.01249.x.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

The authors extend their gratitude to the health workers and study participants that made this investigation possible. We also thank Stephen Dumler, Lirong Shi, Stephen Schachterle, and Nishit Bhuva for their assistance with PCR protocols, and Andrea Wirtz, Sue Tuddenham and Darshan Sudarshi for their assistance with the early phases of this study, and Manny Martinez for his technical assistance with the Figures. Fieldwork was supported by the operating budgets of the participating community-based organizations in Myanmar and the Global Health Access Program. The Johns Hopkins Center for Public Health and Human Rights and the Johns Hopkins Malaria Institute contributed labour and expenses to assist with PCR analyses.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Adam K Richards.

Additional information

Competing interests

All authors declare that they have no competing interests.

Authors’ contributions

TB carried out the laboratory work, conducted data analysis and literature review, and drafted the manuscript. LS assisted with conception and design of the field study, oversaw fieldwork, takes responsibility for the integrity of the field data, and edited the manuscript; EK and KS supervised fieldwork in Karen State and Kachin State, respectively, and critically reviewed the manuscript; DS assisted with the design of the study, supervised PCR laboratory work and critically reviewed the manuscript; TL and CB assisted with the conception and design of the field study and edited the manuscript; AR conceived and designed the study, conducted data analysis and literature review and drafted the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

12936_2012_2529_MOESM1_ESM.docx

Additional file 1: Brown. Supplementary Appendix 1.docx. Peer reviewedcross-sectional studies reporting prevalence estimates for pfmdr1 copy number or pfcrt K76T haplotypes in Burma and adjacent areas of neighboring countries, 1996-2009. Summarizes cross-sectional molecular studies of antimalarial drug resistance included on the map in Figure1. (DOCX 16 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Brown, T., Smith, L.S., Oo, E.K.S. et al. 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 11, 333 (2012). https://doi.org/10.1186/1475-2875-11-333

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1475-2875-11-333

Keywords