Focused Screening and Treatment (FSAT): A PCR-Based Strategy to Detect Malaria Parasite Carriers and Contain Drug Resistant P. falciparum, Pailin, Cambodia

Recent studies have shown that Plasmodium falciparum malaria parasites in Pailin province, along the border between Thailand and Cambodia, have become resistant to artemisinin derivatives. To better define the epidemiology of P. falciparum populations and to assess the risk of the possible spread of these parasites outside Pailin, a new epidemiological tool named “Focused Screening and Treatment” (FSAT), based on active molecular detection of asymptomatic parasite carriers was introduced in 2010. Cross-sectional malariometric surveys using PCR were carried out in 20 out of 109 villages in Pailin province. Individuals detected as P. falciparum carriers were treated with atovaquone-proguanil combination plus a single dose of primaquine if the patient was non-G6PD deficient. Interviews were conducted to elicit history of cross-border travel that might contribute to the spread of artemisinin-resistant parasites. After directly observed treatment, patients were followed up and re-examined on day 7 and day 28. Among 6931 individuals screened, prevalence of P. falciparum carriers was less than 1%, of whom 96% were asymptomatic. Only 1.6% of the individuals had a travel history or plans to go outside Cambodia, with none of those tested being positive for P. falciparum. Retrospective analysis, using 2010 routine surveillance data, showed significant differences in the prevalence of asymptomatic carriers discovered by FSAT between villages classified as “high risk” and “low risk” based on malaria incidence data. All positive individuals treated and followed-up until day 28 were cured. No mutant-type allele related to atovaquone resistance was found. FSAT is a potentially useful tool to detect, treat and track clusters of asymptomatic carriers of P. falciparum along with providing valuable epidemiological information regarding cross-border movements of potential malaria parasite carriers and parasite gene flow.


Introduction
The border region between Cambodia and Thailand, specifically in Pailin province, has repeatedly been the epicenter of emerging resistance of Plasmodium falciparum to antimalarial drugs since the 1960s [1]. This started with resistance to chloroquine (CQ) first reported in 1959 [2], followed by resistance to sulfadoxine-pyrimethamine in 1980's [3] and, more recently, mefloquine [4,5]. In response to increasing levels of drug resistance, Cambodian and Thai national malaria control programmes adapted their national treatment policies. Artemisinin-based combination therapies (ACTs) were implemented as the first-line treatment for uncomplicated falciparum malaria in Thailand in 1995 and in Cambodia in 2000. During the last decade, results from annual routine therapeutic efficacy studies conducted in Cambodia have demonstrated that the efficacy of artesunate-mefloquine and artemether-lumefantrine for the treatment of falciparum malaria was declining in several areas of the country [6][7][8]. In 2007, the World Health Organization (WHO) initiated an informal consultative workshop in Phnom Penh, Cambodia, during which available evidence regarding the emergence of P. falciparum AR parasites was reviewed [9]. Following this meeting, a comprehensive clinical study clearly demonstrated that parasite clearance times were uniformly delayed in patients in Pailin compared to patients at the Myanmar-Thailand border, justifying further research and action [10]. In this context, in order to prevent the possible spread of P. falciparum AR parasites, the ''Artemisinin-Resistance Containment Project'' was implemented by the governments of Cambodia and Thailand with coordination and technical support by WHO and partners with funding from the Bill & Melinda Gates Foundation. Its overall goal was to contain artemisinin-tolerant P. falciparum parasites by removing selection pressure and reducing and ultimately eliminating falciparum malaria. In order to achieve this goal, an ambitious programme of activities was carried out, including the provision of universal coverage with long-lasting insecticide treated nets (LLINs) and of rapid access to free, parasitologically confirmed diagnosis through rapid diagnostic test (RDT) and ACTs delivered by community-based Village Malaria Workers (VMWs).
Initially, a much wider operation referred to as ''Mass Screening and Treatment'' (MSAT), which aims at screening and treating the entire population of Pailin, estimated at over 70,000, within two weeks was envisaged but abandoned following two pilot studies conducted in November 2008 and May 2009. It was therefore decided to carry out to evaluate Focused Screening and Treatment (FSAT), a PCR-based strategy intended to detect P. falciparum asymptomatic carriers. PCR-based cross sectional screenings were focused on villages with the highest annual incidence rates for acute P. falciparum malaria infections reported by heath centers (HC) through malaria microscopy and by VMWs using RDT. Based on high-throughput molecular diagnosis of malaria, the present paper compiles information obtained through the FSAT pilot with a view to gaining more micro-epidemiological knowledge about the population of asymptomatic carriers which would otherwise have remained largely undetected.

Study sites & Participants
Based on the observation that the incidence of parasitologically confirmed malaria cases in general and that of P. falciparum in particular was highly clustered in Pailin, the 109 villages of Pailin were split into two groups: ''high risk'' villages and ''low risk'' villages based on available records of parasitologically confirmed symptomatic cases of P. falciparum malaria (2009). The ''high risk'' village group consisted of the 10 villages reporting the highest P. falciparum incidence rates at HC and of 8 villages reporting the highest number of cases through VMWs in 2009. All others were classified as ''low risk'' villages. As three villages of Pailin qualified for ''high risk'' in both categories (HC and VMW records), only 15 out of Pailin's 109 villages were identified as ''high risk'', the remaining 94 were classified as ''low risk''. All  Field procedures: samples collection and workflow of the screening phase FSAT interventions were conducted at a rate of 3 villages per month. The Screening Team (ST) of seven persons spent four days in each village. The ST was followed at a one-week interval by the Follow-up and Treatment Team of two persons (FTT). Both ST and FTT received local support by the active participation of the village chief in identifying the individuals participating as well as the VMW assisting the direct observation of treatment of those found positive by RDT or PCR. The ST proceeded to the next village and so on until it had screened three villages, before taking one week's rest. During the 4 th week, the FTT completed the third village and rested during week 5. In the week prior to commencement of screening, the ST members visited the village and met the village chief and the VMW. They conducted a detailed joint planning discussion based on the village's logistic and geographical organization (examining the village family register, gaining an idea of how many individuals were likely to turn up at screening, and following recommendations of the village chief on the best locations to set up the screening). Procedures to screen individuals were divided among 4 dedicated stations: Registration and Informed Consent. Two ST members, assisted by the village chief, located each villager in the existing village register. The villager's name was manually transcribed into the FSAT screening register with a unique ID code (village number+family number+letter designating position in the famil-y+screening sequence number). Absent family members were also recorded. In order to assess the risk of spread of drug resistance, the individuals were questioned regarding travel over the past year and classified as ''risk class A'' (travel to or from Myanmar), ''risk class B'' (travel to or from Thailand), ''risk class C'' (travel within Pailin Province only) and ''risk class D'' (travel elsewhere within Cambodia). The screening process was explained to each villager, prior to signing or thumb-printing of the consent form.
Temperature Check. Axillary temperature was recorded for each participant and the result transmitted immediately back to the registration and informed consent station. Individuals with temperature $37.2uC received an additional mark on their arm, indicating that they would require a malaria rapid diagnostic test (RDT).
Laboratory Station. Any individuals detected as being subfebrile or febrile were tested for malaria by using a RDT (CareStartß, Access Bio IncH, Somerset, NJ, USA), detecting falciparum malaria (HRP-2) and non-falciparum malaria (lactate dehydrogenase). If the result was positive, a mark was made on the villager's forearm, identifying him/her for the malaria treatment station. For all participants in the survey, finger prick blood samples were collected and used to perform thick and thin films and double blood spots (one blood spot in 96 well-format PCR plate and one classical blood spot as backup). For quality control, in each 96 well-format PCR plate, one field positive and one negative blood spot controls was included by the ST and blinded to laboratory staff at the Pasteur Institute in Cambodia (IPC) in Phnom Penh, which carried out the PCR tests.
Treatment Station. At the treatment station, RDT positive cases were treated for P. falciparum with atovaquone-proguanil (AP, MalaroneH, GlaxoSmithKlineH) according to manufacturer's instruction or with CQ for non-P. falciparum infections (chloroquine phosphate TabletsH, MedopharmH), according to the national treatment guideline. Prior to treatment, pregnant women were detected through questioning and a pregnancy rapid diagnostic test, if necessary. Pregnant women received quinine (10 mg/kg, eight-hourly67 days). As soon as the patient received the first dose of treatment, a detailed follow-up interview was conducted. RDT negative cases received multivitamins and anti-helminthic treatment, and if symptomatic with fever and other clinical signs of bacterial infection, co-trimoxazole. At the completion of screening in each village, slides and dried blood spot samples were transported within 24 hours to the IPC laboratory in Phnom Penh.
Field procedures: follow-up of PCR-positive cases Follow-up was triggered when the list of PCR positive cases including species identification, identification of Y268N mutation in cytochrome b gene and G6PD status was transmitted to the FTT from IPC by SMS. PCR-confirmed cases of P. falciparum infections were treated with AP in absence of mutation at codon 268 in cytochrome b gene. Direct observed treatment at dosages following manufacturer recommendation was carried out over three days. Malaria cases were followed up daily for the first three days and on day 7 and day 28. Primaquine (PQ) was used according to G6PD status. Individuals with normal G6PD levels and infected with P. falciparum received a single dose of PQ (0.75 mg/kg, Cipla LtdH, India) under supervision on day 2, while normal G6PD infected by P. vivax received CQ over three days according to WHO guidelines and PQ for 14 days (0. 50 mg/kg). Any patients presenting with severe symptoms were given emergency support and referred to the Pailin referral hospital for further assessment and management. Blood was obtained by finger prick on all follow-up days and on any unscheduled day to use for analysis of thick and thin blood smears and for storage on filter paper.

Treatment outcomes
Treatment outcomes were classified as treatment failure if the patient presented with recurrent asexual parasitemia at day 7 or day 28 and treatment success if both day 7 and day 28 asexual parasitemia were negative. According to the national treatment policy, quinine (10 mg/kg, eight-hourly67 days) plus tetracycline (8.3 mg/kg, eight-hourly67 days) was used at Pailin referral hospital to treat P. falciparum recrudescent cases. DHA-PIP (Duo-CotexinH, Zhejiang Holley Nanhu Pharmaceutical Co. LtdH) was used for non-P. falciparum recrudescent cases.

Laboratory procedures
Microscopy examination. Microscopy blood slides were stored in order to cross-check all RDT and PCR positive samples.
Blood smears were stained with 3% Giemsa for 30 min at room temperature. Asexual-and sexual-stage parasite counts were performed by two experienced microscopists. In the case of a discordant result, a third reading was done. Parasite densities were determined from thick blood smears by counting the number of asexual parasites per 200 WBCs (or per 500 WNCs, if the count was less than 10 parasites/200 WBCs), assuming a WBC count of 8,000/ml. A smear was considered negative if no parasites were seen after review of 100 high-power fields.
DNA extraction. DNA from blood spot collected in 96-well PCR plate format containing patient screening samples and field positive and negative controls were extracted with Instagene Matrix resinH (Bio-RadH, France), according to the protocol adapted from the suppliers recommendations and previously described by Steenkeste et al. [11].
Nested PCR detection. Detection of Plasmodium was carried out using primers targeting the Plasmodium cytochrome b gene as described by Steenkeste et al. [11]. For positive samples, RFLP on PCR products using AluI enzyme was performed for determining Plasmodium species. Briefly, 5 ml of PCR products adjusted to 500 ng/ml were mixed with 0.2 ml of AluI and 2.5 ml of buffer A according to the manufacturer's instructions (New England BiolabsH, France), incubated for 4 hours at 37uC and inactivated at 65uC for 15 minutes. Bands were detected by standard 2% agarose gel electrophoresis and ethidium bromide staining. Detection of atovaquone resistance (codon 268). All P. falciparum positive samples were screened for the Y268 mutations in cytochrome b gene that have been related to atovaquone resistance. PCR products from the outer PCR used to detect malaria parasites (1385 pb) were used as template for two second round PCRs. For both nested PCR (NsiI and SspI), amplification took place in the following reaction mixture: 2.5 ml of 106 buffer, 1.5 mM MgCl2, 0.2 mM each deoxynucleoside triphosphate, 0.5 mM each primer (forward primer NsiI, 59-ggtttacttggaacagtttttaacaatg-39, reverse primer NsiI, 59-ggtttacttggaacagtttttaacaatg-39 and forward primer SspI, 59-acagaataatctctagcacc-39, reverse primer NsiI, 59-acctgaatggtactttctacaatat-39), 2 U of FirePol Taq polymerase (Solis BiodyneH, Estonia), and 2 ml of PCR products. Nested PCRs were performed under the following conditions: heating at 94uC for 5 min, followed by 30 cycles of heating at 94uC for 30 s, 45uC (nested PCR NsiI) or 55uC (nested PCR SspI) for 90 s, and 72uC for 2 min, and a final extension period at 72uC for 10 min. Five ml of nested PCR products adjusted to 500 ng/ml were respectively mixed with 0.2 ml of NsiI and 2.5 ml of buffer 3 (nested PCR NsiI) and with 0.2 ml of SspI and 2.5 ml of buffer 2 (nested PCR SspI) according to the manufacturer's instructions (New England BiolabsH, France), incubated for 4 hours at 37uC and inactivated at 65uC (SspI) or 80uC (NsiI) for 20 minutes. Bands were detected by standard 2% agarose gel electrophoresis and ethidium bromide staining. Polymorphism at codon 268 was assessed according to the number and the size of the bands: Y268Y (tat) -NsiI RFLP (2 bands, 359 and 25 bp) and SspI (2 bands, 151 and 23 bp); Y268N (aat) -NsiI RFLP (2 bands, 359 and 25 bp) and SspI (1 band, 174 bp); Y268C (tgt) or Y268S (tct) -NsiI RFLP (1 band, 384 bp) and SspI (2 bands, 151 and 23 bp).
Sequencing reactions. Plasmodium species identification and mutations at codon 268 were secondly confirmed by sequencing. Briefly, nested PCR products were purified by 96-well plate filtration (MilliporeH, France) using a polyacrylamide gel (Bio-Gel P-100H, Bio-RadH, France). Sequencing reactions were carried out with the ABI Prism BigDye Terminator cycle sequencing ready reaction kit run on a 3730 XL genetic analyzer (Applied BiosystemsH, France). Electrophoregrams were visualized and analyzed with Seqscape software v2. 0 (Applied BiosystemsH, France). Amino acid sequences were compared with the following sequences: GenBank accession no. NC002375 for P. falciparum, AY598138 for P. vivax, NC007232 for P. knowlesi, AF069624 for P. malariae and AB182496 for P. ovale, using BioEdit Sequence Alignment Editor software [12].
Detection of G6PD deficiency. To determine which individuals should receive primaquine treatment, all individuals who tested PCR positive were screened for G6PD deficiency from backup blood spots according to the protocol developed by Kuwahata et al [13]. As this method turned out to be unreliable in identifying truly G6PD deficient individuals, by systematic overreporting, only those found to be surely not G6PD deficient received primaquine, while an unknown number of false positives in the test had to be excluded from receiving the additional gametocidal treatment.
P. falciparum genotyping. Molecular genotyping techniques were used to distinguish recrudescence from new infection for all patients who failed therapy after day 7. Filter paper blood samples collected on the day of enrolment and on the day of failure were analyzed for polymorphisms in the genes for merozoite surface protein 1 (msp-1), merozoite surface protein 2 (msp-2) and glutamate rich-protein (glurp) using nested-PCR as described elsewhere [14]. In accordance with WHO recommendations [15], treatment failure was labeled recrudescence if all msp1, msp2 and glurp alleles present at the time of failure had been present at the time of treatment initiation. In all other cases, the failure was considered a new infection.
P. falciparum bar-coding. Bar-coding assays were performed as previously described by Daniels et al. [16] with minor modifications. Polymorphisms in eight SNPs were assessed according to Daniels et al.'s assay numbering (assay number 3,4,7,8,9,19,20,24) by nested PCR approach (Table S1). Second PCR amplification was immediately followed by HRM ramping from 65uC to 85uC with fluorescence data acquisition set at 0.2uC increments. HRM analysis was performed by software supplied by BioradH (CFX Manager 2. 1H, Bio-RadH, France), assigning automatically the two possible nucleotides for each sample. In each run, samples were analyzed in duplicate with 3D7 and Dd2 genomic DNA controls and no template control.

Statistical analysis
Data were entered and verified using Microsoft Excel softwareH and analyzed using XLSTAT for Windows XPH (AddinsoftH, France) and MedCalc software (version 9. 1. 0. 1, Belgium). The Mann-Whitney U test or Kruskal-Wallis method were used for non-parametric comparisons, and Student's t test or one-way analysis of variance for parametric comparisons. For categorical variables, Chi-squared or Fisher's exact tests were used to assess significant differences in proportions. The efficacy outcomes were 28-day clinical and parasitological failure risks adjusted by genotyping, in accordance with the 2009 WHO protocol [17]. Hypothesis testing was made using risk differences, exact 95% confidence intervals, and P-values. A P-value (two-tailed) of less than 0.05 was considered statistically significant.

Ethical clearance
The study protocol was reviewed and approved by the Cambodian National Ethics Committee on Health Research (NECHR) and by the Ethics Board of the WHO Western Pacific Region office (WPRO). An informed written consent was provided by the parents/guardians of all patients before testing.

Population screening
The FSAT study was conducted from May to October 2010. A total of 6,931 individuals were screened in the 20 selected villages. The coverage achieved was 72.7% when counting all persons present at the village at the moment of FSAT and their relatives who were known but absent (6,931/9,537). The mean age of the screened population was 23.7 years (SD, 617.1 years), ranging from 1 week to 89 years. The mean gender ratio (male/female) was 1.04, ranging from 1.32 (in Oh Tatus and Oh Trengleu villages) to 0.88 (Bor Thmey). Risk classification according to the movement of individuals showed that none of individuals declared that they had come from or were be likely to go to Myanmar (''risk class A''); 139 (1.6%) were identified to have a travel history to or from Thailand and were classified as belonging to ''risk class B''; 8550 (94. 6%) declared that they had stayed or were likely to stay within Pailin or local transmission areas within Cambodia (''risk class C''); and 341 (3.8%) declared that they had come from or were likely to proceed to non-malaria transmission areas within Cambodia (''risk class D''). Three individuals (1 in Phnom Dambang and 2 in Oh Preus) were said to work as maids in Malaysia and were initially falsely classified as ''risk class A''. However, given that they were most likely to work in an urban area they were considered not be posing a risk to the spread of malaria parasites from Pailin, they were removed from ''risk class A'' and remain unclassified. No significant differences were observed between villages with regard to risk classes. In the 20 villages, 621 individuals (8.9%) were febrile or sub-febrile with axillary temperature ranging from 37.2uC to 40.0uC, with a mean temperature of 37.4uC. Among them, only 7 (1.1%) were detected positive by RDT: two with P. falciparum or mixed P. falciparum/non-P. falciparum species (40 and 120 parasites/mL) and 5 with non-P. falciparum species (2400, 4000, 4000, 12,000 and 16,000 parasites/ mL). Ten additional malaria cases were found by using PCR (17/ 621, 2.7%): 3 with P. falciparum (32, 35 and 40 parasites/mL) and 7 with P. vivax (30,32,35,40,48, 53 and 187 parasites/mL). Detailed information including microscopic results is presented in Table 1.

Malaria epidemiology
Malaria positive cases. Six out of the 20 villages screened were free of malaria parasites: Rothkros Chhes, Phnom Rang, Oh Treng, Andong Py, Phitas Sbov and Phnom Spung. Temporal evolution of the proportion of malaria parasites carriers showed significant difference according to the month of the collection for the most prevalent species, P. falciparum (P = 0.0001) and P. vivax (P = 0.0003) ( Figure 2). All malaria positive cases had declared having stayed or were likely to stay within Pailin or local transmission areas within Cambodia.
P. falciparum cases. Sixty P. falciparum infections were found in 9 villages (11 out of the 20 villages were P. falciparum cases free). The prevalence of P. falciparum infections was not significantly different between adults (0.94%) and children under the age of 15 (0.74%) (P = 0.36) and between males (1.02%) and females (0.72%) (P = 0.18).
According to the classification of villages, the mean prevalence of asymptomatic parasite carriers was not significantly different between ''high risk'' villages and ''low risk'' villages based on HC 2009 data (0.38% vs. 0.59%, P = 0.41). In contrast, there was a statistically significant correlation between the mean prevalence of parasite carriers of ''high risk'' and ''low risk'' villages identified by VMW in 2009 (1.41% vs. 0.59%, P = 0.04). While the FSAT investigation was being conducted from May to end of October 2010 detecting almost exclusively asymptomatic parasite carriers, the disease surveillance system in Pailin continued to detect and parasitologically confirm the incidence of symptomatic carriers. At the end of the year a re-classification of high and low risk villages became possible based on the newly available incidence data of 2010 using the same selection criteria as for 2009 incidence reports. Even though the reported incidence of P. falciparum had fallen by 73.4% (from 941 to 250) in 2010, 11 out of the initial 15 high risk villages remained among the ten with highest case reports in the two HC and VMW disease surveillance systems (Malaria Consortium, unpublished data). The analysis of the FSAT results showed that the correlation between the prevalence of asymptomatic carriers and disease incidence reports (both surveillance systems) of the same year was highly significant. ''High risk HC 2010'' villages had on average six times more asymptomatic carriers than ''low risk HC 2010'' villages (1.57% vs. 0.24%, P,10 24 ). A very similar result was observed when comparing ''high risk VMW 2010'' villages with ''low risk VMW 2010'' villages (1.48% vs.0.24%, P,10 24 ) (Table S2). Based on 2010 cases report, the prevalence was significantly higher in ''high risk HC/VMW 2010'' villages compared to ''low risk HC/VMW 2010'' villages (1.10% vs. 0.27%, P = 0.005). P. falciparum prevalence was not significantly different between adults and children (0.94% vs. 0.74%, P = 0.36) and between males and females (1.02% vs. 0.72%, P = 0. 18). The spatial distribution of P. falciparum prevalence according to malaria incidence reported through HC and VMW in 2009 (Panel 1) and in 2010 (Panel 2) is shown in Figure 3.
P. falciparum genetic diversity and gene flow.    Thmey and Pekiri (F st = 0.00139, P.0. 05), highlighting a parasite flow between these two villages. P. vivax cases. Seventy two P. vivax infections were found in 12 villages (8 out of the 20 villages screened had no P. vivax cases). P. vivax infections were significantly more prevalent in adults (1.24% vs. 0.74%, P = 0.046), but no significant difference was found between males and females (1.13% vs. 0.80%, P = 0.15). Spatial distribution of P. vivax prevalence is shown in Figure 4 (Panel B).
G6PD deficiency screening. Screening of G6PD deficiency by using blood spots among the malaria infected cases gave results which were considered skewed towards artificially low G6PD enzyme activity levels due to methodical problems, as explained in the discussion section. The global prevalence of G6PD deficient people among parasites carriers was estimated to 63% (18% with 10%-30% of enzyme activity and 45% ,10% of enzyme activity).
Polymorphism at codon 268 in cytb gene. Among the 60 P. falciparum cases, all parasites were wild-type (Y268).
Treatment outcomes of P. falciparum cases. Five different antimalarial regimens were given to 57 of 60 P. falciparum cases: AP (n = 25), AP+PQ (n = 25), CQ (n = 1, wrongly given), QN (n = 1, pregnant woman) and DHA-PIP (n = 5, at the beginning of the study when AP was not available). Three individuals were lost (had left their villages to Preah Vihear, n = 1 or Kampot provinces, n = 2). The 28-day follow-up data were available for 51 patients. One patient treated with DHA-PIP failed and was classified as treatment failure following genotyping. One patient treated with AP developed vivax malaria at day 28.
Treatment outcomes of P. vivax cases. For P. vivax cases, 59 patients received CQ alone and only 11 patients received CQ with PQ due to the many false positive G6PD tests. Two individuals were lost (had left their villages to Prey Veng, n = 1 or Battambang provinces, n = 1). The 28-day follow-up data were available for 64 patients. Both CQ and CQ+PQ were 100% effective in clearing P. vivax infections.

Discussion
The emergence of artemisinin resistance in western Cambodia demanded urgent action to contain its spread and eliminate malaria in identified foci [21]. To contain this threat, a number of strategic priorities were identified and implemented rapidly, such as the mass deployment of LLINs, the banning artemisinin monotherapies and the strengthening of public sector provision of free parasitological diagnosis and treatment through VMWs and HCs. However, there was a paucity of epidemiological information on which to base the planning of additional interventions and a lack of operational experience in Cambodia and elsewhere in implementing potential new tools for elimination. At this stage PCR based diagnosis was limited to epidemiologically explore the reservoir of parasites in asymptomatic carriers.
The most encouraging epidemiological insight gained from FSAT was to that seems extremely low. Out of the 6,931 persons screened and among the 2,606 family relatives of the screened who were absent at the time of screening, none had any travel history to or from Myanmar. Although there is long history of presence of Burmese gem miners in Pailin [22][23][24], two widely corroborated phenomena may explain the complete absence, since several years now, of any Burmese in Pailin: (i) the former abundance of rubies and others gems has now been depleted and (ii) alternative attractive economic opportunities do not exist in Pailin for foreigners, due to lower salary levels when compared to similar activities on the Thailand side of the border.
There had also been concerns that artemisinin resistance might spread indirectly from Palin to Myanmar via temporary or permanent residents of Pailin when migrating to areas of local transmission on the Thailand side of the border, to work in orchards, rubber or plantations or in the fishing industry at the coast in which many Burmese are working [25][26][27]. Here again, the results were encouraging as only as 139 (1.47%) of the 9537 persons registered through the survey qualified to be included into the ''risk group B'' for AR spread as ''likely to go to malaria transmission areas in Thailand''. All malaria positive cases reported to have stayed in Pailin or local transmission areas in Cambodia in the preceding year. Only 8/6931 persons tested (0.12%) had a travel history to Thailand, all of them being PCR negative. All others did not have any travel history that could pose a risk to containment of P. falciparum spread outside of Pailin. However, temporary migrations of seasonal workers from Cambodia to Thailand seem to be likely responsible for the parasite gene flow that we have observed between villages along the Cambodian-Thailand border (i.e. Bor Thmey) and villages inside Pailin province (i.e. Pekiri).
Until the launch of the FSAT project, there were concerns that the existing health information system (HIS), and in particularly the HC incidence data would not accurately identify foci of infection because it is known that patients in Cambodia seek treatment in the private sector, especially if they are far away from public health facilities [28][29][30]. In addition, the villages with higher burdens of symptomatic malaria -(indicative of relatively non-immune populations) might be different from those villages in which the inhabitants were more immune and had higher levels of asymptomatic carriage. These villages might therefore be harboring a large parasite reservoir. However our results showed a good correlation between the high burden villages as identified by the HIS data on incidence of symptomatic cases of 2010. The passive case detection incidence data recorded by the VMWs was already significantly identifying the villages with higher burdens of asymptomatic carriage in 2009 as well as in 2010. The VMWs are clearly a critical component of the containment strategy. As shown here, not only are they instrumental in reducing population parasite burden by providing easy access to good quality diagnosis and treatment, but they are a critical component of the surveillance system [31]. Currently, they only operate in the 43 out of 109 villages in Pailin, which are located more than 5 km away from the nearest health center. Given their important role they can play in containment, expansion of the scheme to include all villages in Pailin should be considered [32].
FSAT proved achievable logistically, the complex protocol of screening, treatment and follow up of over 6'000 individuals was sustained over six months. In addition, FSAT has generated scientific results on patterns of infection and quantified the risks of potential drug resistant parasites spread outside of Pailin. This operational success under difficult field conditions suggests that FSAT might be a feasible method for cross-sectional screening with PCR and directly observed treatment in other settings. A high level of case identification and treatment was achieved despite there being a mean lapse of 8 days between the day of screening and the case being located by the follow-up team. Out of the P. falciparum cases that were detected by PCR, 95% could be traced for follow-up. Of those who could remain in the study all adhered to the directly observed three days of treatment. Of these cases, 96% were followed up until day 28. P. falciparum parasites circulating in Pailin province were found to be sensitive to atovaquone (none of them being mutant-type at codon 268) and asexual parasites totally sensitive to AP. Only one case of treatment failure was observed from an individual treated with DHA-PIP. This observation was not surprising due to high level of recrudescence (,25%) recently found in this area (Leang et al., submitted 2012).
A number of important lessons were learnt which may help to inform the future programme implementation. The main challenge was the lack of accurate information on the size and location of the population because of the amount of population movement in the area. As noted elsewhere, much of this movement is dictated by the seasonal agricultural calendar [26]. The estimates in population size therefore varied considerably making it difficult to accurately calculate coverage rates. Population estimates based on a previous census data were on average 43.7% higher (range: 18.5-66.7%) than the estimates made by the FSAT team with the help of village leaders and based on the population estimates during FSAT, the mean coverage rate was 72.1% (range: 61.0%-91.4%). Transporting blood samples from the field for analysis in the laboratory in Phnom Penh also proved problematic, not only because of the resulting delay between screening and treatment of positive cases, but also the accuracy of some of the results, particularly the G6PD results. The high levels of G6PD deficiency reported here compared to the likely true prevalence of around 10%, [33] were thought to be due partly to the high threshold (,30% of the normal enzyme activity) chosen to declare a person as deficient. In addition, a deterioration of the samples during transportation to Phnom Penh, a day's drive away from the villages was very likely. For this reason a mobile laboratory has since been purchased to improve the speed and the reliability of FSAT. In Cambodia, PCR can now be carried out in the field by placing the mobile laboratory within the target villages reducing the interval time between sample collection and treatment of positive cases to 24 hours. In addition, the detection of mutation at codon 268 in cytochrome b gene and G6PD deficiency screening can now be similarly carried out within the same time and expected increases in reliability. Moreover, the blood sample analysis within the village is likely to render this practice more acceptable and interesting for the population, increasing the coverage rates.
Given the impossibility of carrying out MSAT at the scale originally envisaged, mainly for reasons of lack of human resources in Cambodia, a strong case has been made in favour of Mass Drug Administration (MDA) throughout Pailin as a radical measure to rapidly eliminate P. falciparum from the province. The arguments in favour of MDA include the slowness of FSAT (screening speed of one village per week) and the possibility that there would be asymptomatic carriers under the low parasitaemia threshold of PCR of 0.1 parasites/ml who might escape detection. However, the main difficulty encountered during FSAT (i. e. that of high population movements into, within and out of Pailin) would also be encountered by MDA. But, in contrast to FSAT and MSAT, MDA would also make any epidemiological analysis impossible and effectively render the programme blind as to the location of last remaining foci of transmission. Finally, the rapid geographical widening of drug pressure caused by MDA, followed by a longer period of mass exposure of sub-lethal plasma concentrations of the drug used, present simultaneously in large sections of the population, could trigger the very selection and emergence of multi-drug resistant parasite populations that many so rightly fear. This would be particularly true for AP if it was to be used for MDA, as only one single mutation in the parasite's genome could render this drug ineffective as have successively become all its predecessors throughout Pailin [34].
In conclusion, FSAT may be a useful additional tool in the containment of AR parasites, among the other strategies implemented, albeit its shortcomings in terms of slowness and coverage in three important ways: (i) by providing information on travel and migrations patterns, (ii) by identifying clusters of asymptomatic carriers and (iii) by reducing the overall parasite biomass through the identification and treatment of infections that would otherwise have been undetected. However, recent a report demonstrating the emergence of AR malaria along the Myanmar-Thailand border emphasize the need to expand containment efforts and to further improve surveillance strategies [35]. These might include PCR-based FSAT as new epidemiological tool to help direct containment efforts or malaria elimination in other contexts [36,37] Supporting Information