Research | Open | Published:
Chloroquine efficacy studies confirm drug susceptibility of Plasmodium vivax in Chennai, India
Malaria Journalvolume 13, Article number: 129 (2014)
Assessing the Plasmodium vivax burden in India is complicated by the potential threat of an emerging chloroquine (CQ) resistant parasite population from neighbouring countries in Southeast Asia. Chennai, the capital of Tamil Nadu and an urban setting for P. vivax in southern India, was selected as a sentinel site for investigating CQ efficacy and sensitivity in vivax malaria.
CQ efficacy was evaluated with a 28-day in vivo therapeutic study, while CQ sensitivity was measured with an in vitro drug susceptibility assay. In both studies, isolates also underwent molecular genotyping to investigate correlations between parasite diversity and drug susceptibility to CQ. Molecular genotyping included sequencing a 604 base pair (bp) fragment of the P. vivax multidrug resistant gene-1 (Pvmdr1) for single nucleotide polymorphisms (SNPs) and also the amplification of eight microsatellite (MS) loci located across the genome on eight different chromosomes.
In the 28-day in vivo study (N=125), all subjects were aparasitaemic by Day 14. Passive case surveillance continuing beyond Day 28 in 22 subjects exposed 17 recurrent infections, which ranged from 44 to 148 days post-enrollment. Pvmdr1 sequencing of these recurrent infections revealed that 93.3% had identical mutant haplotypes (958M/Y976/1076L) to their baseline Day 0 infection. MS genotyping further revealed that nine infection pairs were related with ≥75% haplotype similarity (same allele at six or more loci). To test the impact of this mutation on CQ efficacy, an in vitro drug assay (N=68) was performed. No correlation between IC50 values and the percentage of ring-stage parasites prior to culture was observed (rsadj: -0.00063, p = 0.3307) and the distribution of alleles among the Pvmdr1 SNPs and MS haplotypes showed no significant associations with IC50 values.
Plasmodium vivax was found to be susceptible to CQ drug treatment in both the in vivo therapeutic drug study and the in vitro drug assay. Though the mutant 1076L of Pvmdr1 was found in a majority of isolates tested, this single mutation did not associate with CQ resistance. MS haplotypes revealed strong heterogeneity in this population, indicating a low probability of reinfection with highly related haplotypes.
Assessing the Plasmodium vivax burden in India is complicated by the potential threat of an emerging chloroquine (CQ) resistant parasite population . Plasmodium vivax resistance to CQ is rampant in some regions of Southeast Asia [2, 3] and Western Pacific , while the emergence of CQ resistance within other regions of the world, including Central and South America [5–7], Southern Asia  and even the Indian Subcontinent [9–13] is less evident and represented in the literature mainly in the form of medical case reports. Similar to other countries (reviewed in ), CQ resistant Plasmodium falciparum emerged in India in the 1970s [15, 16]. Decades of research have linked P. falciparum drug resistance with a variety of single nucleotide polymorphisms (SNPs) causing non-synonymous amino acid substitutions in the P. falciparum CQ resistance transporter gene (Pfcrt)  and copy number variation in the P. falciparum multidrug resistance gene 1 (Pfmdr1) . However, mutations within Pfcrt appear to be more associated with CQ resistant phenotypes than copy number variation in Pfmdr1[19–23]. Comparative studies between Pfcrt and the P. vivax ortholog, P. vivax candidate gene 10 (Pvcg10), have not revealed similar functionality . Yet, there is some evidence that specific SNPs causing non-synonymous amino acid substitutions within the P. vivax multidrug resistance protein 1 (Pvmdr1), the ortholog to Pfmdr1, may be associated with drug resistant phenotypes .
The Pvmdr1 gene has been characterized in several global genetic diversity studies in Thailand, Indonesia, Azerbaijan, Turkey, French Guyana, Brazil, Madagascar, Mauritania and most recently in China and India [26–36]. Over a dozen non-synonymous amino acid mutations have been reported [26–35]; of these, Y976F and F1076L are most frequently reported and correlated with CQ resistance [27, 34, 35, 37], although much work remains to irrefutably link these mutations with CQ resistance. In India, a recent study in Kolkata reported an absence of the Y976F mutation in 25 P. vivax samples determined to be CQ sensitive taken from patients enrolled in a 28-day in vivo efficacy study . Despite the fact that most drug resistant cases remain confined to specific regions in Southeast Asia and the Western Pacific, the impact of long-term CQ exposure on a parasite population is indisputable and has been studied extensively in P. falciparum. Recently, Mallick et al. showed that a mutant (SVMNT) SNP haplotype of Pfcrt was found to predominate during 2002-2006 in regions with high P. vivax transmission, in conjunction with high CQ exposure during this time. Due to the fact that CQ remains the first line of defense against P. vivax, this highlights the epidemiological impact on the selection of drug resistant parasites in mixed species infections and the maintenance of genetic diversity in a population through inbreeding . Though only a few cases of P. vivax CQ resistance have been observed in India [9–13], a selective sweep with a CQ resistant parasite population would have serious consequences for the control of vivax malaria [39–41].
Tamil Nadu, a state located in southern India, is one of several that are heavily burdened by P. vivax. The capital city, Chennai, has historically seen high P. vivax transmission. For these reasons, Chennai was selected as the sentinel site for investigating CQ sensitivity in P. vivax. In this study, the therapeutic efficacy of CQ was evaluated with an in vivo study, while CQ sensitivity was measured with an in vitro CQ drug assay. In both studies, isolates also underwent molecular genotyping to investigate if a correlation exists between parasite genetic diversity and drug susceptibility to CQ. Molecular genotyping included sequencing a 604 base pair (bp) fragment of Pvmdr1 for SNP haplotype generation and the amplification of eight microsatellite (MS) loci located across the genome on eight different chromosomes. No detectable CQ resistance was found in the patient samples tested in this study, and there was no evidence of previously correlated Pvmdr1 resistance alleles gaining in frequency in the parasite population. However, valuable population level genetic diversity information was obtained, and this pilot study paves the way for future studies on the impact of primaquine (PQ) on relapsing infections.
Field location and study design
Samples were collected from patients attending the Central Malaria Laboratory (CML) in George Town, Chennai, in the state of Tamil Nadu, southeast India. George Town, a predominantly residential/commercial area of the city, is densely populated (~51,000 people/km) and mainly comprised of business community members of a high socio-economic status. However, this area also attracts a large number of laborers mainly from the southern districts of Tamil Nadu, and the population is also composed of different ethnic groups from other states of India such as Rajasthan, Bihar, Gujarat, Haryana, Maharashtra, Uttar Pradesh and Nepal. Immigration combined with rapid urbanization within the catchment area of the CML (~2 km or 100,000 residents) maintains local endemicity in this region. The CML is under the administrative control of the Health Department of the Municipal Corporation of Chennai and serves individuals from all economic strata through passive case detection (PCD). Approximately 15,000 blood smear examinations are performed each year.
Samples collected for this study had ethical clearance from the Ethical Committee of the National Institute of Malaria Research (ICMR) and the Institutional Review Board of New York University Langone Medical Center. All participants provided informed consent and/or assent. Inclusion criteria were specific for each study and are detailed below. Patients were treated as per the National Vector Borne Disease Control Programme (NVBDCP) guidelines in India as given in the study methods.
In vivo therapeutic efficacy study
A total of 856 patients were screened for malaria at the CML Malaria Clinic between the start of enrollment in mid-July 2010 and the completion of this study in mid-November 2010. Of these patients, 125 were enrolled in the 28-day in vivo therapeutic efficacy study (Figure 1A). Inclusion criteria included: >10 and <70 years of age, non-gravid women, patients with no existing concurrent infection(s) (including mixed species determined by microscopy), parasitaemia >500 asexual parasites/μl blood (to ensure enough parasite DNA for downstream experiments), no administration of anti-malarials or antibiotics within three weeks of enrollment and consent for 28-day follow up. The CML utilizes convenience sampling, which often results in an over-representation of males in the study population; a higher frequency of women and children seek treatment at private clinics. Subjects were enrolled on Day 0 and administered 1,500 mg of CQ over 3 days administered by the CML clinic. Subjects were followed on Days 1, 2, 3, 7, 14, 21, 28 and by passive case surveillance for any clinic visits beyond Day 28 (Day of Recurrence), ranging from 44 to 148 days post-enrollment. PQ (0.25 mg/kg) for radical cure was given on Day 28 and administered for 14 days under direct observation of a medical officer. For each subject, thick and thin smears were collected on Days 0, 2, 7, 14, 21 and 28 and examined by expert microscopists. Blood spots for DNA extraction and isolate genotyping were collected before treatment on Day 0 and on any subsequent day when the subject presented with fever.
To investigate patterns in genetic diversity between infection pairs, DNA was extracted from blood spots and used for (i) amplification and sequencing of a 604 base pair (bp) fragment of Pvmdr1 to determine the presence of SNPs previously associated with CQ resistance (T958M, Y976F and F1076L)  and (ii) MS genotyping to describe the relatedness between the pairs , described below.
In vitro CQ drug assay
The in vitro drug study, carried out in August 2009 and August-September 2010, enrolled 96 subjects (Figure 1B). Inclusion criteria included: >1 and <70 years of age, non-gravid women, patients with no existing concurrent infection(s) (including mixed species determined by microscopy) and no administration of anti-malarials or antibiotics within 3 weeks of enrollment. Due to the quantity of material collected, no inclusion criterion was established for parasitaemia level. For each subject, thick and thin smears and 3-4 ml of intravenous blood were collected. For each 3-4 ml of patient venous blood, host white blood cells were removed by CF11 filtration within three hours of collection , and the packed red blood cells were divided as follows: 1 ml was cryopreserved in Glycerolyte 57 solution (Baxter, Deerfield, IL), 200 μl were spotted onto filter paper and 800 μl was used for the in vitro drug susceptibility assay as described [44, 45]. Parasite counts taken before and after CF11 filtration revealed a significant reduction in mature blood stage parasites; consequently, no additional synchronization was performed in an effort to maximize immature ring-stage count.
Drug plates were prepared by coating the wells with CQ di-phosphate (Sigma) dissolved in water, serially diluted from a maximum concentration of 514 ng/ml to a minimum of 8 ng/ml, dried in a non-humidified incubator overnight at 37°C, and stored in the dark at 4°C. Batches of plates were tested for efficacy of the drug by using two P. falciparum lab strains from the National Institute of Malaria Research (NIMR) Parasite Bank, MRC2 (CQ sensitive) and RKL9 (CQ resistant). Approximately 200 μl of a 2% haematocrit blood medium mixture consisting of McCoy’s 5A media (Gibco) and 20% AB + human serum was added to each well, including control wells, and patient samples were tested in duplicate. Plates were incubated at 37°C in a candle jar for 36 hours, a thick smear prepared for each of the wells, stained with Giemsa or Jaswant-Singh-Bhattacherji (JSB) stain and parasitaemia quantified by microscopy. The number of schizonts per 200 asexual stage parasites was determined for each slide and the result for each drug concentration normalized to the control well. Only healthy schizonts with six or more distinct chromatin dots were quantified. Dose-response data were then analyzed using nonlinear regression analysis (HN-NonLin v.1.1, USAMC-AFRIMS, Bangkok, Thailand), and the 50% inhibitory concentration (IC50) obtained.
A selection of isolates used in the in vitro drug assay was identified for further characterization by genotyping. Similar to the in vivo assay isolates, DNA was extracted from blood spots and used for (i) amplification and sequencing of a 604 base pair (bp) fragment of Pvmdr1 and (ii) MS genotyping to describe the relatedness between the pairs , as described below.
Genetic diversity analysis
DNA extraction and species-specific PCR
DNA extraction was performed using QIAamp® DNA Mini Kit (Qiagen Inc., Valencia, CA). Species-specific polymerase chain reaction (PCR)  was used to confirm P. vivax microscopy diagnosis and rule out mixed species infections.
Amplification and sequencing of Pvmdr1
Previously published primers and methods were used to amplify a 604 base pair (bp) fragment of Pvmdr1, capturing three single nucleotide polymorphisms that cause non-synonymous amino acid changes (T958M, Y976F and F1076L) associated with CQ resistance in P. vivax in some studies. PCR products were analyzed by gel electrophoresis and sequenced on an ABI 3730xl sequencer (Applied Biosystems, Foster City, CA) with ≥ 2x coverage using BigDye Terminator v 3.1 and using both forward and reverse primers. Sequences were viewed using Chromas 2.33 (Technelysium Pvt Ltd), and subsequently aligned using Seaview 4 .
Microsatellite amplification and analysis
Eight MS markers  were selected from the literature (Additional file 1) after consideration of recommended guidelines . Using previously published methods and PCR conditions, MS were amplified in 20 μl reactions with ~ 40 ng of extracted DNA using. All reactions were performed individually and not in multiplex. Forward (5′-3′) oligonucleotides were labeled with the phosphoramidite conjugate 6-FAM (Eurofins MWG Operon, Huntsville, AL) and the P. vivax laboratory strain Salvador I was used as a positive control. PCR amplicons were analyzed on an ABI 3730xl sequencer using GeneScan-500 LIZ size standard (Applied Biosystems, Foster City, CA) for size determination. Bands <200 relative fluorescence units (rfu) were excluded and used to define the background. Isolates with greater than one locus having multiple peaks were defined as multiclonal infections, and minor peaks less than 1/3 the height of the major peak were excluded .
Expected heterozygosity (H e ) was used to quantify the amount of genetic diversity within each MS locus in this study, but not across compiled haplotypes of loci due to low sample size. H e was calculated with the following standard formula:
where n is the number of infections sampled and p i is the frequency of the ith allele. For the maximum expected probability of allelic combinations ranging from one to eight loci, let P equal the maximum probability of A, highest frequency alleles, to the nth locus, defined as:
Due to the small sample size, it is not possible to accurately estimate the allele frequencies for the population. To account for this potential error, the standard deviation (SD) of the allele frequency was estimated by a parametric bootstrap procedure, which generates random samples from a population assumed to have allele frequencies equal to their maximum-likelihood values. The SD of each allele frequency was then estimated from the resulting distribution of allele frequencies [50, 51]. The resulting SD was used to calculate the standard error (SE) by: , where n is equal to the number of samples with a specific allele.
In vivo therapeutic efficacy study
To determine the in vivo therapeutic efficacy of CQ, 125 patients were enrolled in a 28-day follow up study at the CML in George Town, Chennai. Demographic and clinical data is provided in Table 1A. A total of 107 (85.6%) subjects completed the 28-day follow up, five subjects (4.0%) withdrew from the study and 13 (10.4%) were lost to follow up (Figure 1A). No fever was observed after 48 hours (Day 2) for all 107 subjects completing the study. Eighty-eight (82.2%) subjects had no sexual or asexual parasites by Day 2, and 17 subjects (15.9%) had no sexual or asexual parasites by Day 7. Two (1.9%) subjects remained positive for sexual stage parasites until a time between Day 7 and Day 14, at which point all patients remained negative for P. vivax for the remainder of the 28-day study (Figure 2). Though the study was completed on Day 28, passive case surveillance continued to monitor the enrolled subjects until four months post-enrollment. Twenty-two of the subjects who had completed the 28-day follow up presented with malaria-like symptoms at a time-point after Day 28. These 22 subjects were tested for malaria by microscopy and 17 subjects were found to have recurrent P. vivax infections on a day ranging from 44 to 148 post-enrollment (termed “Day of Recurrence”) (Figure 1A).
Genetic diversity and Pvmdr1 genotyping of in vivo isolates
Day 0 infections that were subsequently positive on Day 2 and/or Day of Recurrence were genotyped. Despite detecting two positive cases on Day 7, no blood was collected due to an absence of fever (see Methods), so these infections were excluded from genotyping experiments. In total, 69 isolates (33 Day 0, 19 Day 2 and 17 Day of Recurrence) were selected for genotyping in the following combinations: (1) Day 0 and 2 paired-infection (n = 15, 30 isolates); (2) Day 0, 2 and 7 paired-infection (n = 1, 2 isolates); (3) Day 0, 2, 7 and Day of Recurrence paired-infection (n = 1, 3 isolates); (4) Day 0, 2 and Day of Recurrence (n = 3, 6 isolates); and (5) 14 Day 0 and Day of Recurrence paired-infection (n = 14, 28 isolates) (Figure 1A). No Day 2 isolates successfully amplified for genotyping or sequencing, despite successful amplification by species-specific diagnostic PCR. Consequently, the only paired-infection suitable for genetic analysis were Day 0 (n = 33) and Day of Recurrence (n = 17) (Figure 1A).
A 604bp fragment of Pvmdr1 was successfully amplified and sequenced for all Day 0 isolates (n = 33); all isolates (100%) were identified as having the mutant 958M/Y976/1076L haplotype. For the Day of Recurrence isolates, 15 of the 17 isolates successfully amplified and were sequenced; 14 (93.3%) of these were identified as the mutant 958M/Y976/1076L haplotype and 1 (6.7%) as the mutant 958M/Y976/F1076 haplotype (Table 2). No new SNPs were detected and sequences have been deposited in GenBank (Accession nos: KC818367-KC818412).
Eight MS markers were used to determine the genetic diversity of the Day 0 (n = 33) and Day of Recurrence (n = 17) isolates (Figure 1A). MS haplotypes were found to be extremely diverse with no common haplotypes within a single group (within the Day 0 samples or within the Day of Recurrence samples). However, when comparing paired isolates collected from the same individual (on Day 0 paired with Day of Recurrence), there were distinct patterns of relatedness, likely indicating relapsing or recrudescing parasite infections. Of the 17 subjects that were parasite positive on Day of Recurrence, MS haplotypes revealed that 58.8% (n = 9) of the infection pairs were related to the Day 0 sample with ≥75% haplotype similarity (same allele at six or more loci) between samples. In fact, there were six instances where all eight (100%) MS alleles matched exactly, one instance where seven (87.5%) alleles matched and two instances with six (75.0%) alleles matched (Table 3). The probability of detecting identical allelic combinations was calculated for each paired sample based on the number of alleles repeated and the observed frequency of each specific allele present within the population (see Methods for calculation parameters and equation). Due to the extensive genetic diversity within the local P. vivax parasite population, it seems unlikely that any Day of Recurrence infection identified as having ≥75% relatedness between sample pairs was due to a new infection, with a probability ranging from 0.001 to 1.98E-05; rather it is more likely to be the result of a relapsing or recrudescing parasite infection.
Although all 17 subjects with recurrent infections began the 14-day PQ drug regimen during their enrollment in the 28-day in vivo study, only 10 were known to complete the recommended therapy. MS data from this study indicated that six of the nine subjects with highly related paired-infections did complete radical PQ therapy. This finding may support recurrence via recrudescence if the parasite population remains susceptible to the drug, or recurrence via relapse if the parasite population has become tolerant or resistant (Table 3).
Of these 17 paired-infections, 13 (76.5%) were single infections on both Day 0 and Day of Recurrence, while two (11.8%) were single on Day 0 and then multiclonal on Day of Recurrence, one (5.9%) was multiclonal on Day 0 and then single on Day of Recurrence, and one (5.9%) was multiclonal on both Day 0 and Day of Recurrence. In all four cases, all possible haplotype combinations were made and compared between paired-infections. Given the high allelic diversity, only two of the four multiclonal infections shared ≥75% haplotype similarity, with both cases sharing six of eight alleles (Table 3).
In vitro CQ drug assay
A total of 96 subjects were enrolled by PCD at the CML to assay for susceptibility to CQ in vitro. Demographic and clinical data is provided in Table 1B. An in vitro CQ drug assay completed on 81 isolates gave reliable IC50 values for 68 isolates (Figure 1B). IC50 values were low for all isolates, ranging from 7.8-30.2 nM. When plotted the distribution of the IC50 values is a non-linear curve, with a geometric mean of 16.3 nM (95% CI: 14.9-17.8 nM).
Unlike reports from other vivax malaria endemic regions, the percentage of ring-stage parasites upon collection and processing was low. The mean parasitaemia was 0.163% (95% CI: 0.122-0.203%) and was largely dominated by trophozoite-stage parasites (90.0%, 95% CI: 88.4-93.6%). This resulted in a reduced ring-stage inoculation of in vitro culture (Mean: 3.23%, 95% CI: 1.90-4.58%). After 36 hours of in vitro culture, no correlation between IC50 values and the percentage of ring-stage parasites prior to culture was observed (rsadj: -0.00063, p = 0.3307).
Genetic diversity and Pvmdr1 genotype of in vitro samples
The interquartile range of the IC50 values was used to group isolates into either Group 1 (n = 40, IQR <75% = 7.8 to 21.9) or Group 2 (n = 18, IQR ≥75% = 22.1 to 30.2), and a total of 20 isolates (10 from each group) were selected for further analysis (Figure 1B). Genotyping isolates across a range of IC50 values helped to gather evidence of Pvmdr1 resistance alleles gaining in frequency in the parasite population.
Pvmdr1 was successfully amplified and sequenced for nine of the 10 isolates in Group 1; eight of these (88.9%) were identified as having the mutant 958M/Y976/1076L haplotype and 1 (11.1%) as having the mutant 958M/Y976/F1076 haplotype. A similarly dominant pattern of the mutant 958M/Y976/1076L haplotype was observed in nine of the Group 2 isolates (100%). These sequences have been deposited in GenBank (Accession nos: KC818349-KC818366).
The same 20 isolates were genotyped using eight polymorphic MS to obtain background genetic diversity data in the form of the mean number of alleles per locus, the distribution of alleles per locus and the multiplicity of infection (MOI) (Table 2). Both groups had high levels of genetic diversity with no common MS haplotypes within or between the two groups. Though expected given the small sample size, there was no significant difference between the mean number of alleles per locus for Group 1 isolates was 5.89 (range = 4 to 9, SD = 1.55) and for Group 2 isolates 6.5 (range = 3 to 8, SD = 1.93). The distribution of specific alleles between the two groups also did not differ significantly (ranging from r2 = 0.17 and p = 0.18, to r2 = 0.06 and p = 0.60); that is, there did not appear to be a distribution of alleles that was specific to either group. Similarly, there was no significant difference between the total number of clones, multiclonal infections or MOI between Group 1 (n = ≥14, n = 4 and MOI = 1.40, respectively) and Group 2 (n = ≥12, n = 2 and MOI = 1.20, respectively) (Table 2).
A city-wide view of P. vivax genetic diversity in Chennai
Since the subjects for the in vivo and in vitro studies were collected from the same clinic, with sample collection occurring during the same period (2009-2010), MS data were collated into one database to provide a more comprehensive perspective of the parasite population-level diversity. In total 70 samples were genotyped, 20 from the in vitro CQ drug assay and 50 samples from the in vivo therapeutic efficacy study. For each locus the number of alleles, size range, frequency of the most common allele, expected heterozygosity (H e ) and the MOI was calculated (Table 4). On the whole, the mean number of alleles per locus was high (Mean = 13.3, SD = 4.71, range = 5 to 20), indicating high levels of P. vivax genetic diversity, which was confirmed by calculating H e (Mean = 0.856, SD = 0.063, range = 0.729 to 0.912). As expected, the number of alleles per locus was positively correlated with the H e (p = 0.0008, r2 = 0.86, ANOVA) and negatively correlated with the frequency of the most common allele per locus (p = 0.0029, r2 = 0.79, ANOVA). Consistent with this order, the frequency of the most common allele per locus was also negatively correlated with the H e (p = 0.0002, r2 = 0.92, ANOVA) and the MOI (p = 0.0206, r2 = 0.62, ANOVA). MOI was only significantly correlated with the frequency of the most common alleles for each locus, indicating that it is not necessarily the amount of diversity present within a high transmission population, but rather how the diversity is dispersed within the population that influences the accumulation of multiple clones within a single infection.
In this study, P. vivax CQ sensitivity was investigated among infected patients in Chennai, Tamil Nadu. A longitudinal in vivo efficacy study was used to monitor the clearance of parasites and recurrence of infection, and an in vitro CQ drug assay was used to monitor the dose-response effect on parasites isolated from patient samples (Table 1 and Figure 1). Molecular genotyping was used in both studies to investigate the genetic structure of the parasites and to determine if a correlation between parasite genetic diversity and drug susceptibility to CQ existed.
The in vivo efficacy study demonstrated P. vivax sensitivity to CQ in Chennai. Of the 125 subjects enrolled, 107 completed follow up and were microscopy negative for asexual parasites by Day 7; however, two subjects remained positive for sexual parasites until a time prior to the Day 14 collection. Seventeen subjects remained microscopy negative until a time point after Day 28, called Day of Recurrence, which ranged from 44 to 148 days after the initial enrollment (Day 0) (Figure 2). Isolates from Day 0 and Day of Recurrence were genotyped to determine relatedness between the infecting clones. Of the 17 paired-infections, nine were found to share ≥75% of the same alleles. Six (66.7%) of the nine subjects experiencing recurrent infections were known to take PQ radical treatment beginning on Day 28 of this study (Table 3). Though no distinction can be made on whether or not these recurrent infections were due to a relapse, recrudescence or reinfection, microscopy examination indicated parasite clearance by Day 7 and the time-to-recurrence for these infections was within the expected timeframe of a tropical relapse [52, 53]. Other studies in India have reported variable P. vivax relapse rates, between 2.2 and 40.1%, which may be linked with the geographic distribution of relapsing phenotypes [54–57]. However, as recently highlighted by White (2011), little is known about these relapsing phenotypes in India, which is a serious cause for concern given that India harbors a significant global P. vivax burden [58, 59]. In an effort to understand these phenotypes, cross-sectional, longitudinal and clinic epidemiological studies are underway at three sentinel sites in India, as part of the Center for the Study of Complex Malaria in India (CSCMi), an International Center of Excellence in Malaria Research .
Results from the in vitro CQ drug assay, which measures the inhibitory concentration (IC50) of drug needed to prevent parasite growth, also revealed P. vivax sensitivity to CQ in the parasite isolates collected from Chennai. IC50 values from this study ranged from 7.8-30.2 nM (N = 68), which is significantly lower than the 220 nM cutoff for CQ resistance define by Suwanarusk et al. in Indonesia and Thailand, as well as the 100 nM cutoff for CQ resistance defined by Druilhe et al.  in isolates from Myanmar. Regional differences in the above cutoff values are most likely due to methodological differences. Fragment sequencing of Pvmdr1 in both studies revealed that the mutant 958M/Y976/1076L was the dominant haplotype, while the mutant 958M/Y976/F1076 was the minor haplotype (Table 2). The majority of isolates containing the 1076L mutation, which is highly prevalent across the globe, is actually predominantly found in regions with P. vivax CQ-sensitivity [26–28, 35, 37, 62, 63]. The detection of these two haplotypes is consistent with a recent report in nearby Nepal, which may indicate relatedness between the parasite populations in India and Nepal . Similar to a genotyping study in Kolkata , the mutant 976F was not detected in any of the Chennai isolates tested in these two studies (Table 2). Previous studies have reported that this mutation is linked with drug resistance [35, 63]; however, subsequent studies have been unable to confirm this association [64–66]. MS genotyping of these P. vivax isolates revealed an unconstrained and highly diverse population of parasites, with no correlation to Pvmdr1 haplotypes or IC50 values. Rather, valuable information was gained about the population-level genetic diversity of extant P. vivax, which will help to direct future population diversity studies with larger sample sizes.
Considering the highly endemic nature of P. vivax in India and proximity to other countries such as Myanmar, Thailand and Indonesia, the emergence of P. vivax CQ resistance is a persistent threat. Although there are few case reports of P. vivax resistance to CQ in India, the potential for a similar devastating situation as happened with the emergence of CQ resistant P. vivax in Southeast Asia should encourage constant surveillance. Urban malaria epicenters, like Chennai, are of particular interest because the local incidence is not only shaped by common environmental factors, such as the seasonality of monsoons and the development of insecticide resistance, but also by population density, rapid urbanization, sanitation practices and inaccessibility to many sites for continuous monitoring. Though it is not completely understood why P. vivax drug resistance has yet to sweep across India, further investigation on this subject is currently underway as part of CSCMi. As India moves into the pre-elimination phase, the need for increased surveillance of drug resistance has never been more important.
Sneh Shalini, Saumyadripta Chaudhuri and Patrick L Sutton joint first author.
Dash AP, Valecha N, Anvikar AR, Kumar A: Malaria in India: challenges and opportunities. J Biosci. 2008, 33: 583-592. 10.1007/s12038-008-0076-x.
Baird JK, Basri H, Purnomo , Bangs MJ, Subianto B, Patchen LC, Hoffman SL: Resistance to chloroquine by Plasmodium vivax in Irian Jaya, Indonesia. Am J Trop Med Hyg. 1991, 44: 547-552.
Schwartz IK, Lackritz EM, Patchen LC: Chloroquine-resistant Plasmodium vivax from Indonesia. N Engl J Med. 1991, 324: 927-
Rieckmann KH, Davis DR, Hutton DC: Plasmodium vivax resistance to chloroquine?. Lancet. 1989, 2: 1183-1184.
Alecrim MG, Alecrim W, Macedo V: Plasmodium vivax resistance to chloroquine (R2) and mefloquine (R3) in Brazilian Amazon region. Rev Soc Bras Med Trop. 1999, 32: 67-68.
Phillips EJ, Keystone JS, Kain KC: Failure of combined chloroquine and high-dose primaquine therapy for Plasmodium vivax malaria acquired in Guyana, South America. Clin Infect Dis. 1996, 23: 1171-1173. 10.1093/clinids/23.5.1171.
Soto J, Toledo J, Gutierrez P, Luzz M, Llinas N, Cedeno N, Dunne M, Berman J: Plasmodium vivax clinically resistant to chloroquine in Colombia. Am J Trop Med Hyg. 2001, 65: 90-93.
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.
Dua VK, Kar PK, Sharma VP: Chloroquine resistant Plasmodium vivax malaria in India. Trop Med Int Health. 1996, 1: 816-819.
Garg M, Gopinathan N, Bodhe P, Kshirsagar NA: Vivax malaria resistant to chloroquine: case reports from Bombay. Trans R Soc Trop Med Hyg. 1995, 89: 656-657. 10.1016/0035-9203(95)90432-8.
Kshirsagar NA, Gogtay NJ, Rajgor D, Dalvi SS, Wakde M: An unusual case of multidrug-resistant Plasmodium vivax malaria in Mumbai (Bombay), India. Ann Trop Med Parasitol. 2000, 94: 189-190. 10.1080/00034980057536.
Nandy A, Addy M, Maji AK, Bandyopadhyay AK: Monitoring the chloroquine sensitivity of Plasmodium vivax from Calcutta and Orissa, India. Ann Trop Med Parasitol. 2003, 97: 215-220. 10.1179/000349803235001868.
Singh RK: Emergence of chloroquine-resistant vivax malaria in south Bihar (India). Trans R Soc Trop Med Hyg. 2000, 94: 327-10.1016/S0035-9203(00)90344-4.
Awasthi G, Satya Prasad GB, Das A: Pfcrt haplotypes and the evolutionary history of chloroquine-resistant Plasmodium falciparum. Mem Inst Oswaldo Cruz. 2012, 107: 129-134. 10.1590/S0074-02762012000100018.
Sharma VP: Battling the malaria iceberg with chloroquine in India. Malar J. 2007, 6: 105-10.1186/1475-2875-6-105.
Sehgal PN SM, Sharma SL, Gogai S: Resistance to chloroquine in falciparum malaria in Assam state, India. J Commun Dis. 1973, 5: 175-180.
Fidock DA, Nomura T, Talley AK, Cooper RA, Dzekunov SM, Ferdig MT, Ursos LM, Sidhu AB, Naude B, Deitsch KW, Su XZ, Wootton JC, Roepe PD, Wellems TE: Mutations in the P. falciparum digestive vacuole transmembrane protein PfCRT and evidence for their role in chloroquine resistance. Mol Cell. 2000, 6: 861-871. 10.1016/S1097-2765(05)00077-8.
Foote SJ, Thompson JK, Cowman AF, Kemp DJ: Amplification of the multidrug resistance gene in some chloroquine-resistant isolates of P. falciparum. Cell. 1989, 57: 921-930. 10.1016/0092-8674(89)90330-9.
Babiker HA, Pringle SJ, Abdel-Muhsin A, Mackinnon M, Hunt P, Walliker D: High-level chloroquine resistance in Sudanese isolates of Plasmodium falciparum is associated with mutations in the chloroquine resistance transporter gene pfcrt and the multidrug resistance gene pfmdr1. J Infect Dis. 2001, 183: 1535-1538. 10.1086/320195.
Chen N, Russell B, Staley J, Kotecka B, Nasveld P, Cheng Q: Sequence polymorphisms in pfcrt are strongly associated with chloroquine resistance in Plasmodium falciparum. J Infect Dis. 2001, 183: 1543-1545.
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.
Dorsey G, Kamya MR, Singh A, Rosenthal PJ: Polymorphisms in the Plasmodium falciparum pfcrt and pfmdr-1 genes and clinical response to chloroquine in Kampala, Uganda. J Infect Dis. 2001, 183: 1417-1420. 10.1086/319865.
Mayor AG, Gomez-Olive X, Aponte JJ, Casimiro S, Mabunda S, Dgedge M, Barreto A, Alonso PL: Prevalence of the K76T mutation in the putative Plasmodium falciparum chloroquine resistance transporter (pfcrt) gene and its relation to chloroquine resistance in Mozambique. J Infect Dis. 2001, 183: 1413-1416. 10.1086/319856.
Nomura T, Carlton JM, Baird JK, del Portillo HA, Fryauff DJ, Rathore D, Fidock DA, Su X, Collins WE, McCutchan TF, Wootton JC, Wellems TE: Evidence for different mechanisms of chloroquine resistance in 2 plasmodium species that cause human malaria. J Infect Dis. 2001, 183: 1653-1661. 10.1086/320707.
Brega S, de Monbrison F, Severini C, Udomsangpetch R, Sutanto I, Ruckert P, Peyron F, Picot S: Real-time PCR for dihydrofolate reductase gene single-nucleotide polymorphisms in Plasmodium vivax isolates. Antimicrob Agents Chemother. 2004, 48: 2581-2587. 10.1128/AAC.48.7.2581-2587.2004.
Barnadas C, Ratsimbasoa A, Tichit M, Bouchier C, Jahevitra M, Picot S, Menard D: Plasmodium vivax resistance to chloroquine in Madagascar: clinical efficacy and polymorphisms in pvmdr1 and pvcrt-o genes. Antimicrob Agents Chemother. 2008, 52: 4233-4240. 10.1128/AAC.00578-08.
Brega S, Meslin B, de Monbrison F, Severini C, Gradoni L, Udomsangpetch R, Sutanto I, Peyron F, Picot S: Identification of the Plasmodium vivax mdr-like gene (pvmdr1) and analysis of single-nucleotide polymorphisms among isolates from different areas of endemicity. J Infect Dis. 2005, 191: 272-277. 10.1086/426830.
Gama BE, Oliveira NK, Souza JM, Daniel-Ribeiro CT, Ferreira-da-Cruz Mde F: Characterisation of pvmdr1 and pvdhfr genes associated with chemoresistance in Brazilian Plasmodium vivax isolates. Mem Inst Oswaldo Cruz. 2009, 104: 1009-1011.
Ganguly S, Saha P, Guha SK, Das S, Bera DK, Biswas A, Kundu PK, Saha B, Ray K, Maji AK: In vivo therapeutic efficacy of chloroquine alone or in combination with primaquine in vivax malaria in Kolkata, West Bengal, India and polymorphism in pvmdr1 and pvcrt-0 genes. Antimicrob Agents Chemother. 2012, 57: 1246-1251.
Lu F, Gao Q, Chotivanich K, Xia H, Cao J, Udomsangpetch R, Cui L, Sattabongkot J: In vitro anti-malarial drug susceptibility of temperate Plasmodium vivax from central China. Am J Trop Med Hyg. 2011, 85: 197-201. 10.4269/ajtmh.2011.10-0070.
Paul RE, Brockman A, Price RN, Luxemburger C, White NJ, Looareesuwan S, Nosten F, Day KP: Genetic analysis of Plasmodium falciparum infections on the north-western border of Thailand. Trans R Soc Trop Med Hyg. 1999, 93: 587-593. 10.1016/S0035-9203(99)90057-3.
Mint Lekweiry K, Ould Mohamed Salem Boukhary A, Gaillard T, Wurtz N, Bogreau H, Hafid JE, Trape JF, Bouchiba H, Ould Ahmedou Salem MS, Pradines B, Rogier C, Basco LK, Briolant S: Molecular surveillance of drug-resistant Plasmodium vivax using pvdhfr, pvdhps and pvmdr1 markers in Nouakchott, Mauritania. J Antimicrob Chemother. 2012, 67: 367-374. 10.1093/jac/dkr464.
Price RN, Marfurt J, Chalfein F, Kenangalem E, Piera KA, Tjitra E, Anstey NM, Russell B: In vitro activity of pyronaridine against multidrug-resistant Plasmodium falciparum and Plasmodium vivax. Antimicrob Agents Chemother. 2010, 54: 5146-5150. 10.1128/AAC.00801-10.
Suwanarusk R, Chavchich M, Russell B, Jaidee A, Chalfein F, Barends M, Prasetyorini B, Kenangalem E, Piera KA, Lek-Uthai U, Anstey NM, Tjitra E, Nosten F, Cheng Q, Price RN: Amplification of pvmdr1 associated with multidrug-resistant Plasmodium vivax. J Infect Dis. 2008, 198: 1558-1564. 10.1086/592451.
Suwanarusk R, Russell B, Chavchich M, Chalfein F, Kenangalem E, Kosaisavee V, Prasetyorini B, Piera KA, Barends M, Brockman A, Lek-Uthai U, Anstey NM, Tjitra E, Nosten F, Cheng Q, Price RN: Chloroquine resistant Plasmodium vivax: in vitro characterisation and association with molecular polymorphisms. PLoS One. 2007, 2: e1089-10.1371/journal.pone.0001089.
Chehuan YF, Costa MR, Costa JS, Alecrim MG, Nogueira F, Silveira H, Brasil LW, Melo GC, Monteiro WM, Lacerda MV: In vitro chloroquine resistance for Plasmodium vivax isolates from the Western Brazilian Amazon. Malar J. 2013, 12: 226-10.1186/1475-2875-12-226.
Lu F, Lim CS, Nam DH, Kim K, Lin K, Kim TS, Lee HW, Chen JH, Wang Y, Sattabongkot J, Han ET: Genetic polymorphism in pvmdr1 and pvcrt-o genes in relation to in vitro drug susceptibility of Plasmodium vivax isolates from malaria-endemic countries. Acta Trop. 2011, 117: 69-75. 10.1016/j.actatropica.2010.08.011.
Mallick PK, Joshi H, Valecha N, Sharma SK, Eapen A, Bhatt RM, Srivastava HC, Sutton PL, Dash AP, Bhasin VK: Mutant pfcrt “SVMNT” haplotype and wild type pfmdr1 “N86” are endemic in Plasmodium vivax dominated areas of India under high chloroquine exposure. Malar J. 2012, 11: 16-10.1186/1475-2875-11-16.
Schneider KA, Kim Y: Genetic hitchhiking under heterogeneous spatial selection pressures. PLoS One. 2013, 8: e61742-10.1371/journal.pone.0061742.
Schneider KA, Kim Y: Approximations for the hitchhiking effect caused by the evolution of antimalarial-drug resistance. J Math Biol. 2011, 62: 789-832. 10.1007/s00285-010-0353-9.
Schneider KA, Kim Y: An analytical model for genetic hitchhiking in the evolution of antimalarial drug resistance. Theor Popul Biol. 2010, 78: 93-108. 10.1016/j.tpb.2010.06.005.
Karunaweera N, Wanasekara D, Chandrasekharan V, Mendis K, Carter R: Plasmodium vivax: paroxysm-associated lipids mediate leukocyte aggregation. Malar J. 2007, 6: 62-10.1186/1475-2875-6-62.
Sriprawat K, Kaewpongsri S, Suwanarusk R, Leimanis ML, Lek-Uthai U, Phyo AP, Snounou G, Russell B, Renia L, Nosten F: Effective and cheap removal of leukocytes and platelets from Plasmodium vivax infected blood. Malar J. 2009, 8: 115-10.1186/1475-2875-8-115.
Russell B, Chalfein F, Prasetyorini B, Kenangalem E, Piera K, Suwanarusk R, Brockman A, Prayoga P, Sugiarto P, Cheng Q, Tjitra E, Anstey NM, Price RN: Determinants of in vitro drug susceptibility testing of Plasmodium vivax. Antimicrob Agents Chemother. 2008, 52: 1040-1045. 10.1128/AAC.01334-07.
Russell BM, Udomsangpetch R, Rieckmann KH, Kotecka BM, Coleman RE, Sattabongkot J: Simple in vitro assay for determining the sensitivity of Plasmodium vivax isolates from fresh human blood to antimalarials in areas where P. vivax is endemic. Antimicrob Agents Chemother. 2003, 47: 170-173. 10.1128/AAC.47.1.170-173.2003.
Rubio JM, Post RJ, van Leeuwen WM, Henry MC, Lindergard G, Hommel M: Alternative polymerase chain reaction method to identify Plasmodium species in human blood samples: the semi-nested multiplex malaria PCR (SnM-PCR). Trans R Soc Trop Med Hyg. 2002, 96 (Suppl 1): S199-204.
Gouy M, Guindon S, Gascuel O: SeaView version 4: a multiplatform graphical user interface for sequence alignment and phylogenetic tree building. Mol Biol Evol. 2010, 27: 221-224. 10.1093/molbev/msp259.
Sutton PL: A call to arms: on refining Plasmodium vivax microsatellite marker panels for comparing global diversity. Malar J. 2013, 12: 447-10.1186/1475-2875-12-447.
Anderson TJ, Haubold B, Williams JT, Estrada-Franco JG, Richardson L, Mollinedo R, Bockarie M, Mokili J, Mharakurwa S, French N, Whitworth J, Velez ID, Brockman AH, Nosten F, Ferreira MU, Day KP: Microsatellite markers reveal a spectrum of population structures in the malaria parasite Plasmodium falciparum. Mol Biol Evol. 2000, 17: 1467-1482. 10.1093/oxfordjournals.molbev.a026247.
Excoffier L, Slatkin M: Maximum-likelihood estimation of molecular haplotype frequencies in a diploid population. Mol Biol Evol. 1995, 12: 921-927.
Rice JA: Mathematical statistics and data analysis. 1995, Belmont, CA: Duxburry Press, 2
Baird JK, Hoffman SL: Primaquine therapy for malaria. Clin Infect Dis. 2004, 39: 1336-1345. 10.1086/424663.
Kitchen SF: Malariology. 1939, London: W B Saunders
Adak T, Sharma VP, Orlov VS: Studies on the Plasmodium vivax relapse pattern in Delhi, India. Am J Trop Med Hyg. 1998, 59: 175-179.
Adak T, Valecha N, Sharma VP: Plasmodium vivax polymorphism in a clinical drug trial. Clin Diagn Lab Immunol. 2001, 8: 891-894.
Gogtay NJ, Desai S, Kadam VS, Kamtekar KD, Dalvi SS, Kshirsagar NA: Relapse pattern of Plasmodium vivax in Mumbai: a study of 283 cases of vivax malaria. J Assoc Physicians India. 2000, 48: 1085-1086.
Rajgor DD, Gogtay NJ, Kadam VS, Kamtekar KD, Dalvi SS, Chogle AR, Aigal U, Bichile LS, Kain KC, Kshirsagar NA: Efficacy of a 14-day primaquine regimen in preventing relapses in patients with Plasmodium vivax malaria in Mumbai, India. Trans R Soc Trop Med Hyg. 2003, 97: 438-440. 10.1016/S0035-9203(03)90082-4.
White NJ: Determinants of relapse periodicity in Plasmodium vivax malaria. Malar J. 2011, 10: 297-10.1186/1475-2875-10-297.
WHO: World malaria report. 2011, Geneva: World Health Organization,http://www.who.int/malaria/world_malaria_report_2011/9789241564403_eng.pdf,
Das A, Anvikar AR, Cator LJ, Dhiman RC, Eapen A, Mishra N, Nagpal BN, Nanda N, Raghavendra K, Read AF, Sharma SK, Singh OP, Singh V, Sinnis P, Srivastava HC, Sullivan SA, Sutton PL, Thomas MB, Carlton JM, Valecha N: Malaria in India: the center for the study of complex malaria in India. Acta Trop. 2012, 121: 267-273. 10.1016/j.actatropica.2011.11.008.
Druilhe P, Brasseur P, Blanc C, Makler M: Improved assessment of Plasmodium vivax response to antimalarial drugs by a colorimetric double-site plasmodium lactate dehydrogenase antigen capture enzyme-linked immunosorbent assay. Antimicrob Agents Chemother. 2007, 51: 2112-2116. 10.1128/AAC.01385-06.
Imwong M, Pukrittayakamee S, Pongtavornpinyo W, Nakeesathit S, Nair S, Newton P, Nosten F, Anderson TJ, Dondorp A, Day NP, White NJ: Gene amplification of the multidrug resistance 1 gene of Plasmodium vivax isolates from Thailand, Laos, and Myanmar. Antimicrob Agents Chemother. 2008, 52: 2657-2659. 10.1128/AAC.01459-07.
Ranjitkar S, Schousboe ML, Thomsen TT, Adhikari M, Kapel CM, Bygbjerg IC, Alifrangis M: Prevalence of molecular markers of anti-malarial drug resistance in Plasmodium vivax and Plasmodium falciparum in two districts of Nepal. Malar J. 2011, 10: 75-10.1186/1475-2875-10-75.
Orjuela-Sanchez P, de Santana Filho FS, Machado-Lima A, Chehuan YF, Costa MR, Alecrim M, del Portillo HA: Analysis of single-nucleotide polymorphisms in the crt-o and mdr1 genes of Plasmodium vivax among chloroquine-resistant isolates from the Brazilian Amazon region. Antimicrob Agents Chemother. 2009, 53: 3561-3564. 10.1128/AAC.00004-09.
Picot S, Brega S, Gerome P, Velut G, de Monbrison F, Cheminel V, Peyron F: Absence of nucleotide polymorphism in a Plasmodium vivax multidrug resistance gene after failure of mefloquine prophylaxis in French Guyana. Trans R Soc Trop Med Hyg. 2005, 99: 234-237. 10.1016/j.trstmh.2004.09.007.
Sa JM, Nomura T, Neves J, Baird JK, Wellems TE, del Portillo HA: Plasmodium vivax: allele variants of the mdr1 gene do not associate with chloroquine resistance among isolates from Brazil, Papua, and monkey-adapted strains. Exp Parasitol. 2005, 109: 256-259. 10.1016/j.exppara.2004.12.005.
We thank the Director, NIMR (ICMR) for providing the necessary facilities and support. The technical support of Mr. N. Bhaskar and Ms. B. Reeta, Technicians and the staff of NIMR field unit are gratefully acknowledged. This study was supported with funds from a Fogarty International Center Global Infectious Disease training grant (D43 TW007884). The content of this manuscript is solely the responsibility of the authors and does not necessarily represent the official views of the Fogarty International Center or the National Institutes of Health. We thank Drs. Bruce Russell and Simon Kang’a for technical support and Dr. Steven Sullivan for editing and proofreading the manuscript. This work would not have been possible without the support of the health staff of the Corporation of Chennai and the patients involved in the study. This paper bears the NIMR publication committee clearance number 006/2013.
The authors declare no competing interests.
AE and JMC designed the study protocol and provided training for the study; SS, SC undertook sample collection, in vitro drug tests and analysis with help from JC and AE. PLS, SC and SS performed molecular biology experiments; PLS undertook data and statistical analysis; AE, JKD and KJR organized logistics for students, provided access to laboratory space and equipment and provide expertise to facilitate the study protocol; NM and NS provided mentoring for SS; PLS, SS, SC, JC and AE wrote the manuscript. All authors read and approved the final manuscript.
Sneh Shalini, Saumyadripta Chaudhuri, Patrick L Sutton contributed equally to this work.