Volume 22, Number 5—May 2016
Plasmodium falciparum K76T pfcrt Gene Mutations and Parasite Population Structure, Haiti, 2006–2009
Several decades since malaria has been eradicated from most Caribbean islands, the vectorborne parasitic disease continues to cause sporadic outbreaks in the region and remains endemic only to the island of Hispaniola (1), which is the location of the Dominican Republic and Haiti. Evidence suggests that, as in Central America north of Panama, the circulating Plasmodium falciparum malaria parasite, which is the dominant malarial species in Haiti and causes illness associated with the highest number of deaths worldwide (http://www.who.int/mediacentre/factsheets/fs094/en/), has remained chloroquine sensitive. The presence of chloroquine-resistant (CQR) parasites in Haiti could have a notable effect on the populace, as well as complicate ongoing efforts of disease control and the ultimate goal of disease eradication. In addition to social and ethical considerations, the eradication of malaria could ultimately aid in alleviating poverty, a particularly critical issue in Haiti, which has the lowest per capita income in the Western Hemisphere (2). The known negative effects of malaria on economic growth and human capital development would be amplified because of cost increases for malaria treatment if drug-resistant parasites were to become endemic.
Considering the potential ramifications of the establishment of endemicity of drug-resistant malaria parasites, several surveys have assessed the presence of CQR parasites in Haiti. These studies have predominantly relied upon detection of mutations in the P. falciparum chloroquine resistance transporter (pfcrt) gene (3–8) as a proxy for possible drug resistance. Although the presence of a single point mutation does not prove clinical resistance, the substitution at position 76 from lysine (K76) to threonine (T76) is a useful surrogate marker. This technique uses the nested PCR amplification of pfcrtgene sequences, which may be easily detected by using DNA extracted from whole blood spotted on filter paper. The filter papers are dried, then stored with desiccant at collection points until conditions for transportation are favorable. Conventionally, the mutation has been detected by mutation-specific restriction–endonuclease digestion of the pfcrt-nested PCR fragments with Apo1; the gene sequence containing the T76 threonine substitution associated with chloroquine resistance is resistant to digestion, and that containing the wild-type lysine (K76) associated with drug sensitivity reveals 99-bp and 46-bp fragments on agarose gel electrophoresis (3–8). Because DNA sequencing has become more widely available and routinely practical, direct sequencing of the nested pfcrt PCR product for the presence of the mutation has been used (5–8).
Surveys in Haiti have intermittently detected parasites harboring the CQR haplotype of the pfcrt gene (6,9). In addition, the presence of drug-resistant parasites as assessed by using in vitro culture has been reported (9). Clinical chloroquine treatment failures have not been reported, however, and sustained transmission of drug-resistant parasites is not believed to have occurred.
To clarify the genetic context of malaria parasites from Haiti, we used molecular barcode approaches to assess the genetic diversity of this parasite population in the context of geographically distinct neighboring parasite populations. The molecular barcode was used to assess the multiplicity of infection and parasite relatedness. This approach has been used to assess relative parasite transmission intensity (10,11) and has the potential to track parasites to better understand their relatedness and sources, including in outbreak investigations (12). In addition, we genetically characterized a subset of the parasite population in Haiti sampled by single-nucleotide polymorphism (SNP) molecular barcoding to determine their relatedness to other parasite populations in the region.
We report the findings of our surveillance for parasites harboring pfcrt CQR haplotypes in patients with suspected malaria at 9 medical sites across Haiti during the 4 years preceding a major earthquake in January 2010. The earthquake destroyed health infrastructure in the country, including the Ministry of Health and Population, killed >220,000 people, and left >1.5 million homeless. In the aftermath of this catastrophe, major efforts were deployed to establish enhanced surveillance systems to detect and prevent the transmission of disease in the affected population (13,14).
Sample Collection and DNA Processing
Surveillance for CQR P. falciparum in Haiti was continuous during 2006–2009 by the Haitian Group for the Study of Kaposi’s Sarcoma and Opportunistic Infections (GHESKIO) at 9 healthcare centers in the municipalities of Jeremie, Jacmel, Les Cayes, Miragoane, Cap-Haitien, Deschappelles, Port-au-Prince, Fort-Liberte, and Port-de-Paix. Filter paper (Whatman 3MM, Whatman Corporation, Florham Park, NJ, USA) was cut into 2 cm × 1.5 cm strips. We cut 4 teeth, each with a width of ≈5 mm, in the lower half of the filter paper. Whole blood was obtained by finger prick for absorption on the filter paper teeth (≈50 μL per tooth), and smears from patients with fever and suspected malaria were prepared for parasite detection. The filter papers were dried and stored at room temperature in sealed bags with desiccant. We periodically transferred samples and blood smears to GHESKIO in Port-au-Prince. All blood smears were reviewed at GHESKIO for the presence of parasites. We extracted DNA from the individual filter paper teeth of samples from patients who had parasite-positive smears either by methanol extraction or by using the QIAamp DNA Blood Mini kit (QIAGEN, Valencia, CA, USA) according to the manufacturer’s instructions.
Molecular Analysis of pfcrt Gene Sequences
Initially, we analyzed samples for chloroquine resistance mutations by nested PCR amplification of the pfcrt gene, then mutation-specific restriction–endonuclease digestion with Apo1 as previously described (3). Positive and negative controls were included in each round of PCR testing: CQR line Dd2, CQS line 3D7, and water alone. Products were resolved and visualized on a 2% agarose gel. Subsequently, we changed the screening method to the direct sequencing of the nested PCR products, either at the Genomics Core of the Cornell University Life Sciences Core Laboratories Core Facility (http://www.biotech.cornell.edu/node/556) or at Macrogen (Rockville, MD, USA). We cloned PCR products using the TOPO TA Cloning Kit (Life Technologies, Invitrogen, Carlsbad, CA, USA) according to the manufacturer’s instructions.
After extraction, molecular barcode data were obtained for a subset of the samples as described previously (10). We applied this approach to 72 samples that had sufficient remaining material for processing and included the 2 samples that initially were identified as possibly polygenomic. In brief, 0.3 ng of extracted template material were used in 5 μL total reaction volumes containing TaqMan Universal PCR Master Mix (2×), no AmpErase UNG (Applied Biosystems, Foster City, CA, USA), and 40× TaqMan MGB assays (Applied Biosystems) that were run on a 7900HT real-time system (Applied Biosystems). The samples were genotyped after PCR amplification based on their end-point fluorescence signals (FAM or VIC). Samples that showed >1 mixed-base SNP call or had >5 missing calls in the 24-SNP molecular barcode were removed from analysis.
We sorted samples for the presence of identical barcodes. Molecular barcodes that shared >96% of their positions were defined as highly genetically related. For comparison to other populations in the general geographic region, we performed spatial Principal Component Analysis (sPCA) using the adegenet function within the PopGenReport package version 2.0 in R version 3.02 (15) to analyze monogenomic parasite barcodes from Panama (n = 37) (12), Colombia (n = 7) (12), and Venezuela (n = 31) (V. Udhayakumar, pers. comm.). These samples were collected during outbreaks in those countries during 2003–2008, 2011–2012, and 2003–2004, respectively.
We analyzed 901 blood samples for the presence of pfcrt mutations associated with chloroquine drug resistance (Table 1; Figure 1). Of those, 899 samples were analyzed either by restriction digest or sequencing. Of 158 samples analyzed by pfcrt PCR product mutation-specific restriction–endonuclease digestion, 156 showed 99-bp and 46-bp fragments characteristic of the wild-type (chloroquine sensitive) pfcrt gene; however, 2 samples repeatedly revealed a fragment resistant to Apo1 digestion. These samples, 1 each from the cities of Jeremie and Les Cayes, were further analyzed for the possible presence of both chloroquine-sensitive and -resistant parasites. Direct sequencing of the PCR products revealed the potential presence of both alleles. We attempted subcloning of PCR products for these potentially mixed samples to isolate individual PCR products for sequencing, but were not successful. SNP-based molecular barcoding (10) revealed the strain from Jeremie to be a unique isolate (Table 2), but results of the analysis for the isolate from Les Cayes could not be interpreted on 2 occasions because of poor amplification. No CQR pfcrt mutations were detected among the remaining 743 samples analyzed solely by direct sequencing. The observation that essentially all samples were confirmed to lack the CQR haplotype for pfcrt by molecular approaches supports the clinical observation that chloroquine remains highly effective in Haiti.
To better understand the genetic relatedness of these parasites endemic to Haiti and discover whether they are related to or different from neighboring parasite populations, we genotyped a subset of samples with sufficient remaining material using an SNP-based molecular barcoding approach. We used this genotyping method on 72 samples collected in 2006 and 2007; 50 (69%) of these yielded interpretable barcodes; 42 (84%) harbored monogenomic and 8 (16%) polygenomic infections (Tables 2, 3). Further comparative analysis of the 42 monogenomic barcodes revealed a high degree of similarity between the parasite genomes (Figures 2,3). Of the 42 parasites, 15 (36%) had identical or nearly identical molecular barcodes, sharing 23 of 24 positions (96% barcode identity) (Figure 3). Furthermore, some of these nearly genetically identical parasites persisted across transmission seasons from 2006 to 2007 (Figure 3; Table 2). Identical parasite barcodes were identified in Deschappelles in 2006 and Les Cayes in 2007; highly related parasites were also identified in Miragoane in 2006 and Les Cayes in 2007 and in Deschappelles in 2006 and 2007. Additional comparative analysis of the barcodes from Haiti with those from Panama, Colombia, and Venezuela by using sPCA suggest that Haitian parasites are distinct from these other parasite populations in the region (Figure 4).
Although these samples represent a small subset of the malarial parasite population of Haiti, these genetic analyses show highly related parasites within Haiti that are distinct from geographically neighboring parasite populations. Furthermore, detection of parasites with identical barcode genotypes in 2006 and 2007 suggests that parasite populations may persist from one transmission season to another, which implies the need for intervention strategies such as targeting parasite reservoirs that may harbor such persisting parasites from one transmission season to the next.
Malaria remains a major cause of illness and death worldwide. With the exception of Hispaniola, the Caribbean islands are free of sustained P. falciparum parasite transmission and report only occasional outbreaks. These outbreaks are likely caused by the importation of parasites by infected persons, leading to local transmission by competent mosquito vectors. However, malaria remains endemic to Haiti and the Dominican Republic on the island of Hispaniola. Consistent with worldwide efforts, the control and eradication of malaria in Hispaniola is deemed warranted and feasible. However, sporadic reports of CQR parasites raise concerns that these efforts may prove more costly and difficult to accomplish than envisioned.
We report on the analysis of 901 samples from 9 sites throughout Haiti collected during the 4 years preceding the magnitude 7.0 earthquake that occurred in January of 2010. These surveillance efforts were quickly reestablished in the wake of the disaster in order to monitor and prevent the transmission of malaria among vulnerable groups (15,16). Despite prior reports of the presence of parasites in Haiti harboring the pfcrt drug resistance mutation, we found that most of these parasites contained the chloroquine-sensitive pfcrt allele. Only 2 samples from southern Haiti contained parasites that possibly harbored a mixed population of parasites, including those containing the allele indicative of drug resistance. An extensive evaluation of these samples from patients in Les Cayes and Jeremie, in which both sensitive and resistant alleles were detected by mutation-specific restriction–endonuclease digestion of nested PCR pfcrt gene products, was inconclusive. Direct sequencing of the nested pfcrt PCR products yielded variable results but also suggested the presence of resistant and sensitive alleles. Cloning to permit sequencing of unique isolates was unsuccessful, as was molecular barcoding of the sample from Les Cayes. Molecular barcoding of the sample from Jeremie, however, revealed a single strain. These findings might be consistent with the presence of parasite populations subclinically harboring the resistant allele. The presence of such low-density infections may be difficult to detect and characterize; therefore, we cannot conclusively rule out the presence of resistant alleles in these 2 samples (16). Nevertheless, considering that only 2 of 901 samples could not be confirmed as harboring only a chloroquine-sensitive haplotype for pfcrt, our findings indicate that nearly all parasites tested were chloroquine sensitive. We therefore conclude that chloroquine resistance was not likely to be present in Haiti during the analyzed period and that chloroquine remains clinically useful in Haiti.
The observation that increasing use of chloroquine after the earthquake in 2010 did not select for or increase the prevalence of these mutations among the parasite population as described by Morton et al. (14) is consistent with the lack of resistant forms of pfcrt among the population sampled in 2006 and 2009 in Haiti. Although a 2011–2012 study conducted by Okech et al. raises concerns regarding drug efficacy related to the persistence of parasites in some patients treated with chloroquine, the resistance mutation was not detected (17). A prior report indicated a high percentage of pfcrt mutant alleles (5/79, 6%) in P. falciparum–positive blood samples in the Artibonite Valley of Haiti in 2006 and 2007 (5). Our surveillance included 199 samples from this region but did not show resistant alleles, indicating the transmission in this region may not have been sustained under chloroquine pressure. However, possible resistance cannot be definitively excluded.
We used the molecular barcode data for parasite population genetic analysis; results suggest that the P. falciparum parasite population in Haiti is highly related genetically. Of the 42 samples analyzed, 15 (36%) shared >23 of the 24 SNPs of the molecular barcode (96% identical barcode). We have shown in a previous study comparing relationships of the molecular barcode and whole-genome sequencing data that molecular barcodes that share >75% of their positions are highly genetically related (11), indicating that samples described here as identical at 23 or 24 markers in the molecular barcode are genetically or nearly genetically identical.
We increased the threshold to 96% in this study to increase our confidence of identifying truly related parasite types. Furthermore, specific parasite types persist from one year to the next among these highly related parasites, suggesting transmission of single clones without recombination and likely limited introduction of new parasite types through travel or migration. These observations of highly related and clonal parasites that endure across years are consistent with decreasing transmission and potential inbreeding among the parasite population in Haiti. They are also consistent with the population structure analysis of Morton et al. (14).
The PCA analysis reveals that the P. falciparum parasite population in Haiti is generally separate from and independent of parasite populations from Panama, Colombia, and Venezuela. However, 1 parasite from Colombia clusters with the parasite population in Haiti. This observation is consistent with the possible transfer of parasites between this site and Haiti. However, these observations are based on a relatively limited number of samples; although we cannot rule out potential sampling bias, these findings are noteworthy.
After a temporary lull in in control efforts related to P. falciparium malaria after the 2010 earthquake, attempts at control with an ultimate goal of eradication have resumed (18,19). Characterization of baseline populations before catastrophic events can support efforts to identify and predict the effects on parasite populations in response to relaxation of control efforts and changed environments. In Haiti, identification of specific parasite genetic types may suggest importation of parasites related to humanitarian efforts or population migration; in addition, the identification of hot spots of transmission of specific parasite types can direct precise application of control efforts, particularly when resources are limited. Our equivocal detection of only 2 patients potentially harboring low-level populations of chloroquine resistance alleles despite screening >900 samples collected during 4 years supports the contention that the sustained propagation of CQR parasites was not occurring in Haiti during the study period; the possibility is raised by Morton et al. that reported cases are of exogenous origin (14). The finding of low genetic diversity in the P. falciparum population in Haiti is also consistent with that of Morton et al. and others (7,14,20). Although encouraging, these findings support the need for continued surveillance during eradication efforts as well as additional studies to better understand the malarial parasite population structure in Haiti.
Dr. Charles is a medical epidemiologist at Weill Medical College of Cornell University, New York, who focuses on building capacity in developing countries through training, research, and the transfer of simple and robust diagnostic technology. His research interests include tuberculosis, HIV/AIDS, and emerging infectious diseases.
We thank Socrates Herrera Valencia and Nicanor Obaldia for the use of DNA barcode information from Colombia and Panama, respectively, in advance of publication in our analysis. We also thank Lindsay C. Morton, Sheila Okoth Akinyi, and Venkatachalam Udhayakumar for the use of DNA barcode information of parasite isolates from Venezuela. We thank David Rosen for technical support and Abdoulaye Djimde and Christopher Plowe for their advice and encouragement.
This study was funded by support from the Global Fund to Fight AIDS, Tuberculosis, and Malaria. Macarthur Charles was supported by a Harold Amos Medical Faculty Development Award (63526) from the Robert Wood Johnson Foundation and from a Mentored Patient-Oriented Research Career Development grant (K23 AI073190). Molecular barcode analysis was accomplished with funding from the Bill and Melinda Gates Foundation in the laboratory of Professor Dyann Wirth at the Harvard T.H. Chan School of Public Health.
This study was approved by the Institutional Review Boards of the Weill Cornell Medical College and GHESKIO Centers.
- Centers for Disease Control and Prevention (CDC). Malaria acquired in Haiti—2010. MMWR Morb Mortal Wkly Rep. 2010;59:217–9 .
- Purdy M, Robinson M, Wei K, Rublin D. The economic case for combating malaria. Am J Trop Med Hyg. 2013;89:819–23 and.
- Djimdé A, Doumbo OK, Cortese JF, Kayentao K, Doumbo S, Diourte Y, A molecular marker for chloroquine-resistant falciparum malaria. N Engl J Med. 2001;344:257–63 .
- Djimdé A, Doumbo OK, Steketee RW, Plowe CV. Application of a molecular marker for surveillance of chloroquine-resistant falciparum malaria.Lancet. 2001;358:890–1 .
- ElBadry MA, Existe A, Victor YS, Memnon G, Fukuda M, Dame JB, Survey of Plasmodium falciparum multidrug resistance-1 and chloroquine resistance transporter alleles in Haiti. Malar J. 2013;12:426.
- Londono BL, Eisele TP, Keating J, Bennett A, Chattopadhyay C, Heyliger G, Chloroquine-resistant haplotype Plasmodium falciparum parasites, Haiti.Emerg Infect Dis. 2009;15:735–40.
- Londono-Renteria B, Eisele TP, Keating J, Bennett A, Krogstad DJ. Genetic diversity in the merozoite surface protein 1 and 2 genes of Plasmodium falciparum from the Artibonite Valley of Haiti. Acta Trop. 2012;121:6–12.
- Neuberger A, Zhong K, Kain KC, Schwartz E. Lack of evidence for chloroquine-resistant Plasmodium falciparum malaria, Leogane, Haiti. Emerg Infect Dis. 2012;18:1487–9.
- Gharbi M, Pillai DR, Lau R, Hubert V, Khairnar K, Existe A, Chloroquine-resistant malaria in travelers returning from Haiti after 2010 earthquake.Emerg Infect Dis. 2012;18:1346–9.
- Daniels R, Volkman SK, Milner DA, Mahesh N, Neafsey DE, Park DJ, A general SNP-based molecular barcode for Plasmodium falciparumidentification and tracking. Malar J. 2008;7:223.
- Daniels RF, Schaffner SF, Wenger EA, Proctor JL, Chang H-H, Wong W, Modeling malaria genomics reveals transmission decline and rebound in Senegal. Proc Natl Acad Sci U S A. 2015;112:7067–72.
- Obaldia N III, Baro NK, Calzada JE, Santamaria AM, Daniels R, Wong W, Clonal outbreak of Plasmodium falciparum infection in eastern Panama. J Infect Dis. 2015;211:1087–96.
- Adamack AT, Gruber B. PopGenReport: simplifying basic population genetic analysis in R. Methods in Ecology and Evolution. 2014;5:384–7.
- Tietje K, Hawkins K, Clerk C, Ebels K, McGray S, Crudder C, The essential role of infection-detection technologies for malaria elimination and eradication. Trends Parasitol. 2014;30:259–66.
- Townes D, Existe A, Boncy J, Magloire R, Vely JF, Amsalu R, Malaria survey in post-earthquake Haiti—2010. Am J Trop Med Hyg. 2012;86:29–31.
- Morton LC, Akinyi S, Huber C, McMorrow M, Chang M, Townes D, Status of chloroquine resistant haplotypes in Plasmodium falciparum parasite populations collected in post-earthquake Haiti. In: Abstracts of the American Society of Tropical Medicine and Hygiene 62nd Annual Meeting; 2013 Nov 13–17; Washington, DC [Abstract 1264]. Am J Trop Med Hyg. 2013;89(Suppl):384 [cited 2016 Mar 03].http://www.astmh.org/ASTMH/media/Documents/2013AbstractBook1250thru1506andIndex.pdf
- Okech BA, Existe A, Romain JR, Memnon G, Victor YS, de Rochars MB, Therapeutic efficacy of chloroquine for uncomplicated Plasmodium falciparum in Haiti after many decades of its use. Am J Trop Med Hyg. 2015;92:541–5.
- Meeting of the International Task Force for Disease Eradication—November 2012. Wkly Epidemiol Rec. 2013;88:75–80 .
- Baleta A. MIM conference focuses on malaria elimination. Lancet. 2013;382:1319–20.
- Maïga-Ascofaré O, Le Bras J, Mazmouz R, Renard E, Falcao S, Broussier E, Adaptive differentiation of Plasmodium falciparum populations inferred from single-nucleotide polymorphisms (SNPs) conferring drug resistance and from neutral SNPs. J Infect Dis. 2010;202:1095–103.
Suggested citation for this article: Charles M, Das S, Daniels R, Kirkman L, Delva G, Destine R, et al. Plasmodium falciparum K76T pfcrt gene mutations and parasite population structure, Haiti, 2006–2009. Emerg Infect Dis. 2016 May [date cited]. http://dx.doi.org/10.3201/eid2205.150359