Skip to main content

Chloroquine resistance is associated to multi-copy pvcrt-o gene in Plasmodium vivax malaria in the Brazilian Amazon

Abstract

Background

The resistance of Plasmodium vivax to chloroquine has become an obstacle to control strategies based on the use of anti-malarials. The current study investigated the association between P. vivax CQ-resistance in vivo with copy number variation and mutations in the promoter region in pvcrt-o and pvmdr1 genes.

Methods

The study included patients with P. vivax that received supervised treatment with chloroquine and primaquine. Recurrences were actively recorded during this period.

Results

Among the 60 patients with P. vivax, 25 were CQ-resistant and 35 CQ-susceptible. A frequency of 7.1% of multi-copy pvcrt-o was observed in CQ-susceptible samples and 7.7% in CQ-resistant at D0 (P > 0.05) and 33.3% in CQ-resistant at DR (P < 0.05). For pvmdr1, 10.7% of the CQ-susceptible samples presented multiple copies compared to 11.1% in CQ-resistant at D0 and 0.0% in CQ-resistant at DR (P > 0.05). A deletion of 19 bp was found in 11/23 (47.6%) of the patients with CQ-susceptible P. vivax and 3/10 (23.1%) of the samples with in CQRPv at D0. At day DR, 55.5% of the samples with CQRPv had the 19 bp deletion. For the pvmdr-1 gene, was no variation in the analysed gene compared to the P. vivax reference Sal-1.

Conclusions

This was the first study with 42-day clinical follow-up to evaluate the variation of the number of copies and polymorphisms in the promoter region of the pvcrt-o and pvmdr1 genes in relation to treatment outcomes. Significantly higher frequency of multi-copy pvcrt-o was found in CQRPv samples at DR compared to CQ-susceptible, indicating parasite selection of this genotype after CQ treatment and its association with CQ-resistance in vivo.

Background

Malaria remains a major public health problem in many countries. Globally, 3.2 billion people are exposed to its transmission in 91 countries and territories around the world. In 2016, an estimated 216 million cases of malaria occurred worldwide. About 4% of estimated cases globally are caused by Plasmodium vivax [1]. In some countries, such as Brazil, the current treatment for P. vivax is based on two anti-malarials, chloroquine (CQ) and primaquine (PQ) [2]. One of the greatest threats to the control and elimination of malaria is the spread of parasites resistant to anti-malarial drugs. Plasmodium vivax resistance to chloroquine (CQRPv) is difficult to detect due to the low level of parasitaemia among parasite carriers and also to distinguish from relapses to reinfections [3].

The first reports of CQRPv were from 1989 [3, 4], 30 years after the resistance report for Plasmodium falciparum. In Brazil, the first reported case of CQRPv was from a patient treated in Manaus, in the Brazilian Amazon [5]. In this city, other subsequent studies assessed the efficiency of standard supervised therapy and the proportion of CQ treatment failures was 10.1% in 2007 [6] and 5.2% in 2014 [7].

pfcrt-o and pfmdr1 are molecular markers involved in CQ resistance in P. falciparum [8]. Several studies sought mutations in these genes orthologs in P. vivax samples and the majority did not show the involvement of these genes with CQ resistance [9, 10], suggesting that in this species the mechanism is probably different. In an in vivo study, in Manaus, P. vivax isolates from patients treated with CQ with treatment failure presented an increased pvcrt-o and pvmdr-1 gene expression on both D0 (day of admission) and DR (day of recrudescence) compared to patients that responded properly to treatment. No significant differences were observed between immature and mature asexual forms in CQRPv and susceptible groups [11]. Other study cited that pvcrt-o transcription is not a primary determinant of ex vivo drug susceptibility, but is rather related to stage of parasite development [12].

An AAG insertion in the exon I of pvcrt-o gene leading to an extra amino acid at position 10 (termed K10 insertion) was identified as a possible marker for CQRPv [13]. In Myanmar, more than half of the samples (63.9%) showed the K10 insertion [14]. During later years there have been an increasing number of whole genome sequences for P. vivax from Latin America where polymorphisms in potential drug resistance genes are mentioned [15, 16]. These studies showed that the genetic basis of drug resistance can be gained simply by examining the genomes of parasites.

This study aimed to investigate the association between CQRPv with copy number variation and mutations in the promoter region in pvcrt-o and pvmdr1.

Methods

Ethics statement

The study was approved by the Ethics Review Board of Fundação de Medicina Tropical Dr. Heitor Vieira Dourado (FMT-HVD) (approval number 426/2011). Participants were informed about the study and signed a consent form. In the case of under 18 years, the study participant signed an assent form and the guardians signed a consent form.

Location of the study

The study was conducted from January 2012 to April 2013 at the FMT-HVD, an infectious disease referral centre located in Manaus, in the Western Brazilian Amazon.

Selection of patients

The study included patients with P. vivax malaria of both genders, aged 6 months–60 years, weight greater than 5 kg, with blood parasite density from 250 to 100,000 parasites/mL and axillary temperature ≥ 37.5 °C or history of fever in the last 48 h. Use of anti-malarials in the previous 30 days, impossibility to be followed up for 42 days and clinical complications were considered exclusion criteria [17]. Patients received supervised treatment with 25 mg/kg of CQ phosphate over a 3-day period (10 mg/kg on days 0 and 7.5 mg/kg on days 1 and 2. PQ was prescribed at the end of follow-up or if malaria recurred during follow-up, at the dosage of 0.5 mg/kg per day, during 7 days.

Patients were evaluated on days 0, 1, 2, 3, 7, 14, 28 and 42 and, if they felt ill, at any time during the follow-up period. Day 0 (D0) is defined as the admission day and day of recurrence (DR) as the day of recrudescence for the same patient. Both samples from the same patient were analysed whenever available. For some patients, only a one-time sample (D0 or DR) was available. Thick blood smears and full blood counts were collected from all patients on each visit. On D0 and in DR, whole blood was also collected for DNA storage. CQ and desethylchloroquine (DCQ) plasmatic levels were determined only in case of parasitological failure. Three aliquots of 100 μL of plasma from DR samples were spotted onto filter paper for later analysis by high performance liquid chromatography to determine the levels of CQ and DCQ, as previously described [2, 18]. CQRPv was defined as a parasite recurrence presenting plasma concentrations of CQ plus DCQ higher than 100 ng/mL when added up. Susceptible-CQ P. vivax group consisted of patients with no parasite recurrence during the follow-up period.

Plasmodium vivax diagnosis

Thick blood smear was performed by the Walker technique [19] and evaluated by an experienced microscopist. Parasite densities were calculated by counting the number of parasites per 500 leukocytes and the number of parasites/µL per patient was determined [20]. Real-time PCR was performed to confirm P. vivax monoinfection. The extraction of total DNA from whole blood was performed using the QIAamp DNA Blood Mini Kit (Qiagen, USA), according to the manufacturer’s protocol. The DNA was amplified in an Applied Biosystems 7500 Fast System using primers and TaqMan fluorescence labeled probes for real-time PCR [21].

Copy number variation of pvcrt-o and pvmdr-1

Quantitative PCR (qPCR) was performed using Taqman® Universal PCR Master Mix (Applied Biosystems) and 100–200 ng of template DNA were prepared from each sample. PCR products were amplified and identified on a 7500 FAST (Applied Biosystems). Sequences of primers used to determine copy number variation of pvcrt-o [22] and pvmdr-1 [23] genes were described in Additional file 1.

Cycling parameters for PCR were an initial denaturation step at 95 °C for 10 min, 40 cycles of 15 s at 95 °C and 1 min at 60 °C. The single-copy β-tubulin gene was used as a reference gene (normalizer), and a field sample with a single copy of the target gene was used as a calibrator. The ΔΔCt method was used to determine the number of copies of pvcrt-o and pvmdr-1 genes relative to a standard calibrator sample. Samples were considered to have a single copy of the gene when the estimated value was 0.5–1.5, while values above 1.5 were defined as gene amplification. Each DNA sample was assayed in triplicate.

To define the calibrators for copy number variation (CNV) analysis using real-time PCR assays, plasmids containing a single copy of each gene: pvcrt-o, pvmdr-1 and pvtubulin were constructed. The presence of a unique copy of each gene was confirmed by sequencing. Using the single-copy plasmids, a set of field samples to define the gene copy number using relative quantification in real-time PCR was evaluated. Thus, samples with a single copy of pvcrt-o and pvmdr-1 were selected to be used as calibrators in the CNV assays.

Nucleotide sequencing of the promoter region of pvcrt-o and pvmdr-1 genes

PCR primers and different reaction conditions used to amplify pvcrt-o and pvmdr-1 gene sequences are described in Additional file 2. Sequencing reaction was performed in triplicate. Nucleotide sequences were analysed using the NCBI BLAST algorithm and compared with reference sequences of GenBank [24]. Amino acid sequences were aligned and compared with sequences of the reference chloroquine transporter from P. vivax Sal1 (GeneBank ID: NC_009906) using Geneious® software (Biomatters, v6.0.5).

Polymorphism analysis of msp1F3, MS2 and MS8 molecular markers for patients with CQRPv

For the genotyping procedures for comparing day 0 and recurrences samples, 3 highly polymorphic genes were chosen according to diversity and discriminatory potential in the region, msp1F3, MS2, and MS8 [25, 26]. PCRs were performed in 20 µL reactions with 10 µM of each primer, 2 µL of 10 × Buffer B (50 mM KCl, 10 mM Tris, pH 8.3), 2.5 mM each dNTP, 25 mM MgCl2, 1.5 U Taq DNA polymerase and 5 µL genomic DNA, as described in Additional file 3. All reagents were used from Invitrogen. 3 µL of diluted primary PCR product was used as template for the nested PCR. PCRs were performed with conditions as follows: initial denaturation at 95 °C for 1 min, followed by 29 cycles (primary PCR) or 24 cycles (nested PCR) of denaturation at 95 °C for 30 s, annealing at 59 °C for 45 s, elongation at 72 °C for 1 min, followed by a final elongation at 72 °C for 5 min. PCR products were sent to Macrogen for capillary electrophoresis-based sequencing and analysed with Gene Marker version 2.6.0. To distinguish reinfection, recrudescence, and relapse by P. vivax, was used the classification method recommended by WHO for P. falciparum [27]. For each recurrence, samples were classified as homologous if at least 1 allele for each locus investigated was detected in both paired samples (D0 and DR) on two different markers and as heterologous if all alleles for a given marker were different on two or more markers [11].

Statistical analyses

Data were analysed using SPSS® version 21.0 for Windows (SPSS Inc.® Chicago, IL, USA). Normal distribution of data was evaluated with the Kolmogorov–Smirnov test. Asexual parasitaemia at D0, D1, D2 and D3 was compared in patients with CQRPv and CQ-susceptible P. vivax using T-Student test. Parasitaemia clearance rates at D1, D2 and D3 were compared in patients with CQRPv and CQ-susceptible P. vivax using Chi square or Fisher’s test. Chi square or Fisher’s test was used to test differences in proportions of polymorphism of the promoter region and copy number variation. A p value of 0.05 or lower was considered to be statistically significant. The genetic variation for each microsatellite locus was measured by calculating the expected heterozygosity (HE). HE was calculated using D0 and DR for each locus as = [n/(n − 1)] [1 − ∑p 2i ] where n is the number of isolates sampled and pi is the frequency of allele i.

Results

Sixty patients with P. vivax were included in the study. Twenty-five patients presented CQRPv and 35 presented CQ-susceptible P. vivax. Clinical and laboratory characteristics of the patients with CQRPv and CQ-susceptible P. vivax are presented in Tables 1, 2. Of the patients with CQRPv, 80% were male, with a mean age of 24.6 years. The mean haemoglobin was 12.8 g/dL at D0 and 12.6 g/dL at DR. The geometric mean parasitaemia was 2553.5 parasites/mm3 at D0 and 4672.1 parasites/mm3 at DR (Table 2). The mean blood levels of CQ plus DCQ were higher than 100 ng/mL in all patients with recurrent infection. Thirty-five patients with CQ-susceptible P. vivax were randomly selected and included as controls (Table 2). Of the patients with CQ-susceptible P. vivax, 77.2% were males, mean age of 33.1 years. The mean haemoglobin was 13.2 g/dL. The geometric mean parasitaemia was 2623.9/mm3.

Table 1 Demographic and clinical characteristics of the study CQRPv admitted to a tertiary health centre, Manaus, Amazon, Brazil
Table 2 Clinical characteristics of the study participants CQRPv and susceptible-chloroquine P. vivax admitted to a tertiary health centre, Manaus, Amazon, Brazil

Asexual parasitaemia mean was similar at D0, D1, D2 and D3 in patients with CQRPv and CQ-susceptible P. vivax (P > 0.05). At D1, parasitaemia cleared in 2.9 and 0% of the patients with CQRPv and CQ-susceptible P. vivax, respectively (P = 0.659); 47 and 59.1% of the patients with CQRPv and CQ-susceptible P. vivax cleared at D2, respectively (P = 0.274); and 85.3 and 81.8% of the patients with CQRPv and CQ-susceptible P. vivax cleared at D3, respectively (P = 0.225).

For the microsatellite analysis at MS1F3, MS2 and MS8 loci, reading of the alleles was successful in 21 of the 25 CQRPv samples analysed. In these samples, 86.0% presented at least one concordant allele in primary infection and recrudescence (Additional file 4). The expected heterozygosity (HE) for msp1F3 was: 0.880 at D0 and 0.780 at DR; for MS2, 0.960 at D0 and  0.940 at DR; and for MS8 was 1.00 at D0 and 0.970 at DR.

Analysis of copy number variation of the pvcrt-o and pvmdr1 genes

For pvcrt-o CNV analysis, with 28 patients with CQ-susceptible P. vivax and 25 patients with CQRPv being analysed (14 samples from D0 to 18 from DR). At D0, 7.1% of the CQ-susceptible P. vivax presented multiple copies of pvcrt-o, 7.7% of the CQRPv showed multiple copies of pvcrt-o at D0 while at DR 33.0% had multiple copies of pvcrt-o (Fig. 1a). There was a significant difference between CQ-susceptible P. vivax and CQRPv at DR (P = 0.02).

Fig. 1
figure 1

Determination of the copy number variation of genes in CQ-susceptible and patients resistant. a Determination of the number of copies for the pvcrt-o gene between the susceptible and resistant groups in D0 and DR. b Determination of the number of copies of the pvmdr-1 gene between the susceptible and resistant groups in D0 and DR. *P < 0.05

For pvmdr-1 CNV analysis, a total of 28 CQ-susceptible and 25 CQRPv (9 from D0 to 17 from DR) P. vivax samples were analysed. It was verified that 10.7% of the CQ-susceptible P. vivax samples, 11.1% of the CQRPv samples at D0 and 0.0% of the CQRPv samples at DR presented multiple copies (Fig. 1b). There was no significant difference between CQ-susceptible and CQRPv groups regarding pvmdr-1 CNV (P > 0.05). Samples that did not meet the criteria used in the methodology were excluded from the analysis of these genes. Overall, pvmdr1 and pvcrt-o CNV was not related to parasitaemia and age for individuals infected with P. vivax with either with multiple copies of pvmdr-1 or with one copy of the genes (P > 0.05).

Analysis of pvcrt-o and pvmdr1 genes polymorphisms in the promoter region

Nucleotide sequencing of a region flanking the promoter and exon I of the gene pvcrt-o showed a deletion of 19 bp in the promoter region at position 32,083–32,102 and an insertion of three bases in exon 1 leading to 5AAG repeats (AAG) after codon nine, leading to the K10 insertion compared to the reference sequence NC_009906-Chloroquine resistance transporter-P. vivax Sal1, with 4AAG repeats termed wild type here (Fig. 2).

Fig. 2
figure 2

Scheme of the pvcrt-o gene showing promoter region and exon 1. This figure show specimens with 19 kb deletion without K10 insertion (1), samples with 19 kb insertion without k10 insertion (2) and samples with 19 kb insert with k10 insertion (3)

A deletion of 19 bp was found in 11/23 (47.6%) of the patients with CQ-susceptible P. vivax and 3/10 (23.1%) of the samples with in CQRPv at D0. At day DR, 55.5% of the samples with CQRPv had the 19 bp deletion. There was no significant difference between groups (P ≥ 0.05, Fisher’s exact test) (Table 3).

All the samples with 19 bp deletion presented 4AAG repeats in the exon I (Table 2). Of the samples with wild type sequence, 8 (33.3%), 8 (60.0%) and 3 (25.0%) exhibited 5AAG (K10) repeats in the exon I of the gene in the samples with CQ-susceptible P. vivax, CQRPv at D0 and CQRPv at DR, respectively (Table 3).

For the pvmdr-1 gene, nucleotide sequencing was performed in 66 samples. Sixteen samples with CQ susceptible P. vivax and 23 CQRPv were sequenced, 14 in the D0 and 9 in the DR. There was no variation in the analysed gene compared to the P. vivax reference Sal-1 available in the NCBI gene bank.

Discussion

Chloroquine and PQ are still the drugs of choice for the treatment of P. vivax in some countries, especially in Latin America [28,29,30,31]. However, CQRPv has emerged in several parts of the world [32, 33], including Malaysia [30], Myanmar [34], India [33] and Brazil [7, 35]. The molecular mechanism of CQRPv is still not well defined. In this study, CNV of pvcrt-o and pvmdr1 was determined in patients with CQRPv from a well-characterized cohort under supervised treatment and this demonstrated a relationship between CQRPv in vivo and increased CNV of the pvcrt-o. These patients presented an effective concentration of CQ and DCQ and microsatellite revealed the presence of the same clonal nature at D0 and DR, indicating recrudescence. In Brazil, for the first time, CNV of the pvcrt-o was evaluated and the amplification of this gene was observed in P. vivax isolates in a study with treatment supervised for 42 days (Table 3).

Table 3 Frequency of the mutation and repetition AAG in the promoter region of the pvcrt-o gene in CQRPv and CQ susceptible P. vivax

In regions where mefloquine (MQ) is used to treatment of P. falciparum infection, drug pressure mediated by increase copy number variation of pvmdr1 may select MQ-resistant P. vivax [36]. Actually, some studies show that pvmdr1 amplification was associated to MQ-resistance [13, 23, 37, 38]. MQ was used in the Brazilian Amazon in the early 2000s, possibly contributing to the selection of this genotype in this region [5].

In this study, pvmdr-1 CNV was not associated to CQRPv. Although CQ was found to be effective in India [39, 40], amplification of the pvmdr-1 gene was found in 31.6% of the Indian isolates [33]. In Brazil, two studies described pvmdr1 copy number variation, with an amplification rate of 0.9% in isolates from Acre state [41] and 20% in isolates from the states of Mato Grosso and Rondônia [22]. However, these studies did not have clinical follow-up to confirm CQRPv. Pvmdr1 Y976F mutation was found in malaria-endemic regions where CQ is used as the first-line treatment [13, 38, 42].

In previous investigations, CQRPv isolates have shown high expression of pvcrt-o and/or pvmdr-1 genes in a malaria reference centre in the Brazilian Amazon [11]. The main factor involved in gene expression regulation is its promoter, since it contains the binding sites for the RNA polymerases responsible for gene transcription [43]. To better elucidate CQRPv, the presence of mutations in the promoter region of the pvcrt-o and pvmdr1 genes in CQ-susceptible P. vivax and CQRPv was investigated. For the pvcrt-o gene, an deletion 19 bp sequence was found in the promoter sequence, followed by 4AAG in the exon I. Although there was no significant difference in the frequency of mutations between these two groups, the CQRPv presented higher frequency of the deletion of 19 bp at D0 compared to at DR.

Here, we also described the occurrence of the wild type sequence in pvcrt-o, followed by 5AAG repeat (K10 insert) in CQ-susceptible P. vivax and CQRPv, with a frequency ranging from 25 to 60%. Similar frequency of insertion was reported by in vitro susceptibility to CQ studies using isolates from Thailand (46%) and Myanmar (56%) [13]. More recent studies have also demonstrated K10 insertion in isolates of P. vivax from India [32, 33] and Ethiopia [44], but the frequency was lower than that found in our study. Among the samples from India, a high frequency of polymorphisms in the pvcrt-o and pvmdr1 genes has been demonstrated also, indicating the circulation of strains that may acquire CQRPv phenotype in that region [32, 33].

A limitation of this study is the small number of samples, restricting the possibility of detecting differences. Polyclonal infections were present, challenging the discrimination of recrudescence from relapse to re-infection. However, it was shown that plasma levels of CQ plus DCQ were higher than the minimum effective limits at DR for all patients. Results from in vivo studies, as the presented in this study, are very useful since they provide the response to treatment in the real life, but are subject to host immune influence, drug metabolism variations and drug quality [42]. Ex vivo assays evaluate drug activity in the absence of host confounding factors [45]. Then, potential influence of host immunity can not be discarded here.

Conclusions

This was the first study with 42-day clinical follow-up to evaluate the variation of the number of copies and polymorphisms in the promoter region of the pvcrt-o and pvmdr1 genes in relation to treatment outcomes. Significantly higher frequency of multi-copy pvcrt-o was found in CQRPv samples at DR compared to CQ-susceptible, indicating parasite selection of this genotype after CQ treatment and its association with CQ-resistance in vivo.

References

  1. WHO. World mlaria report 2017. Geneva: World Health Organization; 2017.

    Google Scholar 

  2. Naing C, Aung K, Win DK, Wah MJ. Efficacy and safety of chloroquine for treatment in patients with uncomplicated Plasmodium vivax infections in endemic countries. Trans R Soc Trop Med Hyg. 2010;104:695–705.

    Article  CAS  Google Scholar 

  3. Rieckmann KH, Davis DR, Hutton DC. Plasmodium vivax resistance to chloroquine? Lancet. 1989;2:1183–4.

    Article  CAS  Google Scholar 

  4. Baird JK, Basri H, Purnomo Bangs MJ, Subianto B, Patchen LC, et al. Resistance to chloroquine by Plasmodium vivax in Irian Jaya, Indonesia. Am J Trop Med Hyg. 1991;44:547–52.

    Article  CAS  Google Scholar 

  5. Alecrim MdG, 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–8.

    Article  Google Scholar 

  6. de Santana Filho FS, Arcanjo AR, Chehuan YM, Costa MR, Martinez-Espinosa FE, Vieira JL, et al. Chloroquine-resistant Plasmodium vivax, Brazilian Amazon. Emerg Infect Dis. 2007;13:1125–6.

    Article  Google Scholar 

  7. Marques MM, Costa MR, Santana Filho FS, Vieira JL, Nascimento MT, Brasil LW, et al. Plasmodium vivax chloroquine resistance and anemia in the western Brazilian Amazon. Antimicrob Agents Chemother. 2014;58:342–7.

    Article  Google Scholar 

  8. Lo E, Hemming-Schroeder E, Yewhalaw D, Nguyen J, Kebede E, Zemene E, et al. Transmission dynamics of co-endemic Plasmodium vivax and P. falciparum in Ethiopia and prevalence of antimalarial resistant genotypes. PLoS Negl Trop Dis. 2017;11:e0005806.

    Article  Google Scholar 

  9. Nomura T, Carlton JM, Baird JK, del Portillo HA, Fryauff DJ, Rathore D, et al. Evidence for different mechanisms of chloroquine resistance in 2 Plasmodium species that cause human malaria. J Infect Dis. 2001;183:1653–61.

    Article  CAS  Google Scholar 

  10. Sa JM, Yamamoto MM, Fernandez-Becerra C, de Azevedo MF, Papakrivos J, Naude B, et al. Expression and function of pvcrt-o, a Plasmodium vivax ortholog of pfcrt, in Plasmodium falciparum and Dictyostelium discoideum. Mol Biochem Parasitol. 2006;150:219–28.

    Article  CAS  Google Scholar 

  11. Melo GC, Monteiro WM, Siqueira AM, Silva SR, Magalhaes BM, Alencar AC, et al. Expression levels of pvcrt-o and pvmdr-1 are associated with chloroquine resistance and severe Plasmodium vivax malaria in patients of the Brazilian Amazon. PLoS ONE. 2014;9:e105922.

    Article  Google Scholar 

  12. Pava Z, Handayuni I, Wirjanata G, To S, Trianty L, Noviyanti R, et al. Expression of Plasmodium vivax crt-o is related to parasite stage but not ex vivo chloroquine susceptibility. Antimicrob Agents Chemother. 2016;60:361–7.

    Article  CAS  Google Scholar 

  13. Lu F, Lim CS, Nam DH, Kim K, Lin K, Kim TS, et al. 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.

    Article  CAS  Google Scholar 

  14. Nyunt MH, Han JH, Wang B, Aye KM, Aye KH, Lee SK, et al. Clinical and molecular surveillance of drug resistant vivax malaria in Myanmar (2009–2016). Malar J. 2017;16:117.

    Article  Google Scholar 

  15. Winter DJ, Pacheco MA, Vallejo AF, Schwartz RS, Arevalo-Herrera M, Herrera S, et al. Whole genome sequencing of field isolates reveals extensive genetic diversity in Plasmodium vivax from Colombia. PLoS Negl Trop Dis. 2015;9:e0004252.

    Article  Google Scholar 

  16. Flannery EL, Wang T, Akbari A, Corey VC, Gunawan F, Bright AT, et al. Next-generation sequencing of Plasmodium vivax patient samples shows evidence of direct evolution in drug-resistance genes. ACS Infect Dis. 2015;1:367–79.

    Article  CAS  Google Scholar 

  17. WHO. Guidelines for the treatment of malaria. Geneva: World Health Organization; 2006.

    Google Scholar 

  18. Ruebush TK 2nd, Zegarra J, Cairo J, Andersen EM, Green M, Pillai DR, et al. Chloroquine-resistant Plasmodium vivax malaria in Peru. Am J Trop Med Hyg. 2003;69:548–52.

    Article  Google Scholar 

  19. WHO. World malaria report: 2010. Geneva: World Health Organization; 2010.

    Google Scholar 

  20. WHO. Malaria microscopy quality assurance manual. Geneva: World Health Organization; 2016.

    Google Scholar 

  21. Snounou G, Viriyakosol S, Zhu XP, Jarra W, Pinheiro L, do Rosario VE, et al. High sensitivity of detection of human malaria parasites by the use of nested polymerase chain reaction. Mol Biochem Parasitol. 1993;61:315–30.

    Article  CAS  Google Scholar 

  22. Costa GL, Amaral LC, Fontes CJF, Carvalho LH, de Brito CFA, de Sousa TN. Assessment of copy number variation in genes related to drug resistance in Plasmodium vivax and Plasmodium falciparum isolates from the Brazilian Amazon and a systematic review of the literature. Malar J. 2017;16:152.

    Article  Google Scholar 

  23. Imwong M, Pukrittayakamee S, Pongtavornpinyo W, Nakeesathit S, Nair S, Newton P, et al. Gene amplification of the multidrug resistance 1 gene of Plasmodium vivax isolates from Thailand, Laos, and Myanmar. Antimicrob Agents Chemother. 2008;52:2657–9.

    Article  CAS  Google Scholar 

  24. BLAST. Basic Local Alignment Search Tool (BLAST: Basic Local Alignment Search Tool).

  25. Beck HP, Wampfler R, Carter N, Koh G, Osorio L, Rueangweerayut R, et al. Estimation of the antirelapse efficacy of tafenoquine, using Plasmodium vivax genotyping. J Infect Dis. 2016;213:794–9.

    Article  CAS  Google Scholar 

  26. de Araujo FC, de Rezende AM, Fontes CJ, Carvalho LH, de Alves Brito CF. Multiple-clone activation of hypnozoites is the leading cause of relapse in Plasmodium vivax infection. PLoS ONE. 2012;7:e49871.

    Article  Google Scholar 

  27. WHO. Methods and techniques for clinical trials on antimalarial drug efficacy: genotyping to identify parasite populations: informal consultation organized by the Medicines for Malaria Venture and cosponsored by the World Health Organization. Amsterdam: World Health Organization; 2008.

    Google Scholar 

  28. Commons RJ, Thriemer K, Humphreys G, Suay I, Sibley CH, Guerin PJ, et al. The Vivax Surveyor: online mapping database for Plasmodium vivax clinical trials. Int J Parasitol Drugs Drug Resist. 2017;7:181–90.

    Article  CAS  Google Scholar 

  29. Negreiros S, Farias S, Viana GM, Okoth SA, Chenet SM, de Souza TM, et al. Efficacy of chloroquine and primaquine for the treatment of uncomplicated Plasmodium vivax malaria in Cruzeiro do Sul, Brazil. Am J Trop Med Hyg. 2016;95:1061–8.

    Article  CAS  Google Scholar 

  30. Grigg MJ, William T, Menon J, Barber BE, Wilkes CS, Rajahram GS, et al. Efficacy of artesunate-mefloquine for chloroquine-resistant Plasmodium vivax malaria in Malaysia: an open-label, randomized, controlled trial. Clin Infect Dis. 2016;62:1403–11.

    Article  CAS  Google Scholar 

  31. Price RN, von Seidlein L, Valecha N, Nosten F, Baird JK, White NJ. Global extent of chloroquine-resistant Plasmodium vivax: a systematic review and meta-analysis. Lancet Infect Dis. 2014;14:982–91.

    Article  Google Scholar 

  32. Bansal D, Acharya A, Bharti PK, Abdelraheem MH, Elmalik A, Abosalah S, et al. Distribution of mutations associated with antifolate and chloroquine resistance among imported Plasmodium vivax in the state of Qatar. Am J Trop Med Hyg. 2017;97:1797–803.

    Article  CAS  Google Scholar 

  33. Joy S, Mukhi B, Ghosh SK, Achur RN, Gowda DC, Surolia N. Drug resistance genes: pvcrt-o and pvmdr-1 polymorphism in patients from malaria endemic South Western Coastal Region of India. Malar J. 2018;17:40.

    Article  Google Scholar 

  34. Htun MW, Mon NCN, Aye KM, Hlaing CM, Kyaw MP, Handayuni I, et al. Chloroquine efficacy for Plasmodium vivax in Myanmar in populations with high genetic diversity and moderate parasite gene flow. Malar J. 2017;16:281.

    Article  Google Scholar 

  35. Siqueira AM, Alencar AC, Melo GC, Magalhaes BL, Machado K, Alencar Filho AC, et al. Fixed-dose artesunate-amodiaquine combination vs chloroquine for treatment of uncomplicated blood stage P. vivax infection in the Brazilian Amazon: an open-label randomized, controlled trial. Clin Infect Dis. 2017;64:166–74.

    Article  CAS  Google Scholar 

  36. Auburn S, Serre D, Pearson RD, Amato R, Sriprawat K, To S, et al. Genomic analysis reveals a common breakpoint in amplifications of the Plasmodium vivax multidrug resistance 1 locus in Thailand. J Infect Dis. 2016;214:1235–42.

    Article  CAS  Google Scholar 

  37. Lin JT, Patel JC, Kharabora O, Sattabongkot J, Muth S, Ubalee R, et al. Plasmodium vivax isolates from Cambodia and Thailand show high genetic complexity and distinct patterns of P. vivax multidrug resistance gene 1 (pvmdr1) polymorphisms. Am J Trop Med Hyg. 2013;88:1116–23.

    Article  Google Scholar 

  38. Suwanarusk R, Chavchich M, Russell B, Jaidee A, Chalfein F, Barends M, et al. Amplification of pvmdr1 associated with multidrug-resistant Plasmodium vivax. J Infect Dis. 2008;198:1558–64.

    Article  CAS  Google Scholar 

  39. Shalini S, Chaudhuri S, Sutton PL, Mishra N, Srivastava N, David JK, et al. Chloroquine efficacy studies confirm drug susceptibility of Plasmodium vivax in Chennai, India. Malar J. 2014;13:129.

    Article  Google Scholar 

  40. Saravu K, Kumar R, Ashok H, Kundapura P, Kamath V, Kamath A, et al. Therapeutic assessment of chloroquine-primaquine combined regimen in adult cohort of Plasmodium vivax malaria from primary care centres in Southwestern India. PLoS ONE. 2016;11:e0157666.

    Article  Google Scholar 

  41. del Vargas-Rodriguez RC, da Silva Bastos M, Menezes MJ, Orjuela-Sanchez P, Ferreira UM. Single-nucleotide polymorphism and copy number variation of the multidrug resistance-1 locus of Plasmodium vivax: local and global patterns. Am J Trop Med Hyg. 2012;87:813–21.

    Article  CAS  Google Scholar 

  42. Suwanarusk R, Russell B, Chavchich M, Chalfein F, Kenangalem E, Kosaisavee V, et al. Chloroquine resistant Plasmodium vivax: in vitro characterisation and association with molecular polymorphisms. PLoS ONE. 2007;2:e1089.

    Article  Google Scholar 

  43. Carvalho TL, Ribolla PE, Curi RA, Mota LS. Characterization and transcriptional analysis of the promoter region of the Duffy blood group, chemokine receptor (DARC) gene in cattle. Vet Immunol Immunopathol. 2009;132:153–9.

    Article  CAS  Google Scholar 

  44. Golassa L, Erko B, Baliraine FN, Aseffa A, Swedberg G. Polymorphisms in chloroquine resistance-associated genes in Plasmodium vivax in Ethiopia. Malar J. 2015;14:164.

    Article  Google Scholar 

  45. Russell B, Chalfein F, Prasetyorini B, Kenangalem E, Piera K, et al. Determinants of in vitro drug susceptibility testing of Plasmodium vivax. Antimicrob Agents Chemother. 2008;52:1040–5.

    Article  CAS  Google Scholar 

Download references

Authors’ contributions

Conceived and designed the experiments: SRS GCM RS SCPL GLC TNS JLFV. Analysed the data: SRS GCM. Wrote the paper: SRS GCM ACGA WMM MVGL. Acquisition of data: GCM AMS. Critical revision of the manuscript: GCM WMM AMS SCPL MVGL. All authors read and approved the final manuscript.

Acknowledgements

We acknowledge Nélida Teresa Sánchez Rúa, Maria Raimunda Costa for sample collection and Sequencing Platform of the Instituto Leônidas & Maria Deane.

Competing interests

The authors declare that they have no competing interests.

Availability of data and materials

The datasets during and/or analysed during the current study available from the corresponding author on reasonable request.

Ethics approval and consent to participate

The study was approved both by the Research Ethics Committee of the Fundação de Medicina Tropical Dr. Heitor Vieira Dourado (FMT-HVD) and by the National Commission for Research Ethics (CONEP) under opinion 426/2011.

Funding

Research Support Foundation of Amazonas (FAPEAM) through PPSUS (protocol 062.00650.2014) and Guyamazon (protocol 26873.UNI554.3814.29062015) project supported this study. MVGL are level 1 fellows from CNPq. SRS have fellowships from FAPEAM. ACGA have fellowships from CNPq

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Gisely Cardoso de Melo.

Additional files

Additional file 1.

Oligonucleotide primers used for copy number variation of P. vivax orthologs genes.

Additional file 2.

Oligonucleotide primers used for promotor region sequencing of P. vivax orthologs genes.

Additional file 3.

Oligonucleotide primers used for genotyping of P. vivax.

Additional file 4.

Infection haplotypes in patients with CQRPv.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Silva, S.R., Almeida, A.C.G., da Silva, G.A.V. et al. Chloroquine resistance is associated to multi-copy pvcrt-o gene in Plasmodium vivax malaria in the Brazilian Amazon. Malar J 17, 267 (2018). https://doi.org/10.1186/s12936-018-2411-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12936-018-2411-5

Keywords