Skip to main content

Genetic analysis of the orthologous crt and mdr1 genes in Plasmodium malariae from Thailand and Myanmar

Abstract

Background

Plasmodium malariae is a widely spread but neglected human malaria parasite, which causes chronic infections. Studies on genetic polymorphisms of anti-malarial drug target genes in P. malariae are limited. Previous reports have shown polymorphisms in the P. malariae dihydrofolate reductase gene associated with pyrimethamine resistance and linked to pyrimethamine drug pressure. This study investigated polymorphisms of the P. malariae homologous genes, chloroquine resistant transporter and multidrug resistant 1, associated with chloroquine and mefloquine resistance in Plasmodium falciparum.

Methods

The orthologous P. malariae crt and mdr1 genes were studied in 95 patients with P. malariae infection between 2002 and 2016 from Thailand (N = 51) and Myanmar (N = 44). Gene sequences were analysed using BioEdit, MEGA7, and DnaSP programs. Mutations and gene amplifications were compared with P. falciparum and Plasmodium vivax orthologous genes. Protein topology models derived from the observed pmcrt and pmmdr1 haplotypes were constructed and analysed using Phyre2, SWISS MODEL and Discovery Studio Visualization V 17.2.

Results

Two non-synonymous mutations were observed in exon 2 (H53P, 40%) and exon 8 (E278D, 44%) of pmcrt. The topology model indicated that H53P and E278D were located outside of the transmembrane domain and were unlikely to affect protein function. Pmmdr1 was more diverse than pmcrt, with 10 non-synonymous and 3 synonymous mutations observed. Non-synonymous mutations were located in the parasite cytoplasmic site, transmembrane 11 and nucleotide binding domains 1 and 2. Polymorphisms conferring amino acid changes in the transmembrane and nucleotide binding domains were predicted to have some effect on PmMDR1 conformation, but were unlikely to affect protein function. All P. malariae parasites in this study contained a single copy of the mdr1 gene.

Conclusions

The observed polymorphisms in pmcrt and pmmdr1 genes are unlikely to affect protein function and unlikely related to chloroquine drug pressure. Similarly, the absence of pmmdr1 copy number variation suggests limited mefloquine drug pressure on the P. malariae parasite population, despite its long time use in Thailand for the treatment of falciparum malaria.

Background

Plasmodium malariae is one of the five important human Plasmodium species, but its genome sequence has just been revealed [1, 2]. Recently, developed highly sensitive molecular diagnosis methods have shown that the prevalence of P. malariae in many malaria-endemic regions is higher than previously assumed [3]. There is limited information on the biology and molecular genetics of P. malariae as well as on anti-malarial drug resistance in this species for which long-term in vitro culture methods are lacking. Small scale in vitro drug susceptibility testing and clinical efficacy [4] showed that P. malariae in Thailand responds well to pyronaridine [5] and artesunate [6] but that fever clearance time was delayed after treatment with chloroquine [7]. However, P. malariae exhibits a 72-hour erythrocytic stage, which suggests that it may require a longer time to clear parasites after treatment. Therefore, the extended clearance time may not be an indicator for chloroquine resistance in P. malariae [8]. In vitro drug testing in P. malariae is further complicated by the low parasitaemias present in patients and the high proportion of mixed infections with other Plasmodium species. An alternative approach for monitoring anti-malarial drug resistance in P. malariae is the assessment of polymorphisms in anti-malarial drug target genes. Since the molecular markers for drug resistance are not well characterized in P. malariae, the orthologous gene markers associated with drug resistance in Plasmodium falciparum and Plasmodium vivax were considered in P. malariae. Earlier studies have reported polymorphisms in the P. malariae dhps, dhfr and kelch orthologues, which relate to sulfadoxine, pyrimethamine and artemisinin resistance, respectively.

Plasmodium malariae infection rate is generally underestimated and overlooked because it is often asymptomatic and mostly mixed with other Plasmodium species. Drug-resistant P. falciparum and P. vivax have been reported in many regions, including Thailand [9, 10] and Myanmar [11, 12]. Chloroquine-resistant P. falciparum was first documented in 1957 in Thailand which is now widely spread around the globe [13]. Polymorphisms of pfcrt and pfmdr1 genes have been linked to chloroquine resistance with the key determinant of K76T mutation in pfcrt gene [10]. Several studies of pfmdr1 supported that amplification and polymorphism were useful for prediction of chloroquine and mefloquine resistance [14,15,16]. The pvmdr1 Y976F mutation has been found to correlate with reduced susceptibility to chloroquine [17, 18]. An increased copy number of pvmdr1 has been found in Thailand [17, 18], suggesting that there is a relationship between pvmdr1 copy number and mefloquine pressure. In the Greater Mekong Sub-region, previous drug pressure from chloroquine and mefloquine on the P. malariae parasite population is expected to be considerable, since these drugs have been used widely for the treatment of vivax and falciparum malaria, respectively. The current study focuses on the P. malariae crt and mdr1 orthologue genes, which in P. falciparum are involved in chloroquine and mefloquine resistance. These genes were evaluated in P. malariae isolated of patients from Thailand and Myanmar between 2002 and 2016.

Methods

DNA extraction

Whole blood samples were collected from 95 P. malariae-infected patients from Thailand (51 samples) and Myanmar (44 samples) between 2002 and 2016 (Table 1). DNA was extracted using QIAamp DNA Mini Kit (Qiagen, Germany) and stored at −80 °C until use. Plasmodium species was confirmed by polymerase chain reaction based on 18 small-subunit ribosomal RNA [19, 20].

Table 1 The samples used in the study

Amplification of orthologous crt and mdr1 genes from Plasmodium malariae

Specific primers for amplification of pmcrt and pmmdr1 were designed based on the reference sequences (accession number LT594622.1 and LT594631.1). The primers and conditions used for amplification of pmcrt and pmmdr1 are listed in Additional files 1 and 2. Semi-nested and nested PCR were carried out in 20 μl, 2 mM of MgCl2, 250 µM of dNTPs, 250 nM of both forward and reverse primers, 0.5U BIOTAQ DNA polymerase (Bioline, UK), and 2.0 μl of DNA template. The approximated concentration of DNA templates were 80 to 150 ng/μl. The primary PCR products were used as a template in the secondary PCR. Cycling conditions were: initial denaturation at 94 °C for 5 min then denaturation at 94 °C for 1 min, annealing for 1 min, and extension at 72 °C for 1 min plus 30 s, using 30 cycles during the first round and 35 cycles during the second round, follow by a final extension of 72 °C for 10 min. The positive PCR products were purified and submitted for DNA sequencing in South Korea (Macrogen Inc., Korea).

DNA sequences analysis of pmcrt and pmmdr1 genes

DNA sequences were analysed using Clustal software in the Bioedit package [21] together with reference sequences of P. malariae accession numbers LT594622.1 (pmcrt) and LT594631.1 (pmmdr1). DNA sequence polymorphisms and haplotypes patterns were analysed using DnaSP version 6.10.04 [22] and MEGA7 [23].

Detection of pmmdr1 copy number variation

To assess pmmdr1 copy number variations, relative quantitative PCR was performed on the Applied Biosystems StepOnePlus™ (Applied Biosystems, USA). The primers for pmmdr1 and pmβtubulin were pmmdr1F (5′- CAGATGTGGGAAACGACAATG-3′), pmmdr1R (5′-TAGAAGCTCCCTCCCCGTTT-3′), pmβ-tubulinF (5′-TGAAGCAACTGGAGGAAGGT-3′), and pmβ-tubulinR (5′-GGACCTGCTCGGACACTATC-3′). SsoFast™ EvaGreen® Supermixes (Biorad, USA) was used as medium, and the thermal cycler profile was prepared according to the manufacturer’s instruction. As calibrator, a single copy control was constructed from a plasmid by insertion of pmmdr1 (nucleotides [nt] 1102 to 1993) and pmβtubulin (nucleotides [nt] 1132 to 1215) fragments in a ratio of 1:1 into pGEM®-T Easy Vector (Promega, USA). In addition, the two-copy pmmdr1 control was constructed by insertion of pmmdr1 and pmβtubulin fragments in a ratio of 2:1 into the pGEM®-T Easy Vector (Promega, USA). Both of the control plasmids were confirmed their insertions by DNA sequencing. The pmβtubulin served as an internal control to calculate the relative amount of pmmdr1 gene by comparing the Ct readings. Copy numbers were calculated as 2−ΔΔCt. All reactions were performed in triplicate. Samples with a copy number more than 1.5 were classified as multiple copies. Reactions were repeated at least twice in case the samples contained ΔΔCt spread > 1.5, a Ct value > 35, or a pmmdr1 copy number higher than 1.3.

Homology modelling of PmCRT and PmMDR1 proteins

The topology structure of PmCRT and PmMDR1 were predicted using Phyre2 [24] and SWISS-MODEL [25]. Homologous sequences of PmCRT and PmMDR1 proteins derived from different haplotypes were searched with PSI-Blast [26]. The secondary structure of PmCRT and PmMDR1 was predicted using Psi-pred 2.5 [27] and Disopred 2.4 softwares [28]. The membrane-spanning domains of PmCRT and PmMDR1 were predicted using Memsat_SVM program [29]. A multi-template approach was selected for prediction of the structural model of PmCRT and PmMDR1 proteins using SWISS-MODEL [25]. The model was evaluated with the VADAR tool and was visualized by Discovery Studio Visualizer V 17.2 [30].

Results

Genetic analysis of pmcrt

Purified PCR products covering the 14 exons of the pmcrt gene from all 95 P. malariae samples were sequenced and analysed. The assembled sequence covered the complete pmcrt coding sequence, with a size of 1272 nucleotides, encoding a protein of 423 amino acids. The multiple sequence alignments showed 9 nucleotide polymorphisms in exons 2, 8, 9, 10 and 12 at nucleotide positions A363C (40%), C346T (18.94%), C353T (53.68%), C394G (4.21%), A403T (28.42%), A2120T (44.2%), T2483T (2.10%), T2637A (1.05%), and A3151T (6.31%). There were 2 non-synonymous mutations in exons 2 and 8 at codons 53 (H53P) and 278 (E278D). The samples collected from Thailand showed more frequent mutations than samples collected from Myanmar. In Thailand, the prevalence of mutations in exons 2 (H53P) and 8 (E278D) were 66.66 and 60.78%, respectively, compared to Myanmar 9.09 and 25%, respectively (Table 2).

Table 2 Summary of pmcrt gene polymorphisms in Thailand and Myanmar

Four haplotypes were characterized from 95 PmCRT sequences (Table 3). All 4 haplotypes were presented in samples from both Thailand and Myanmar. Haplotype 1 (H53P + E278D) was the most prevalent in Thailand (41.18%), whereas in Myanmar haplotype 4 (H53 + E278) was the most prevalent (68.18%). Haplotype 2 (H53E + E278) was observed in higher frequency in Thailand (25.49%, 13/51) than in Myanmar (6.82%, 3/44). There were ten samples (19.61% for Thailand and 22.73% for Myanmar) from both countries showing haplotype 3 (H53 + E278D).

Table 3 Haplotype patterns of PmCRT in Thailand and Myanmar

For comparison, PmCRT haplotypes obtained from this study were aligned with the orthologous genes from human Plasmodium spp. The two-point mutations, H53P and E278D, found in P. malariae samples were not corresponded to the point mutations associated with chloroquine resistance described in P. falciparum. To predict the effect of these two-point mutations on PmCRT function, a topology model of PmCRT was constructed and then compared to PfCRT and PvCRT (Fig. 1). This showed that both H53P and E278D mutations are located outside the transmembrane domain (Fig. 1a), and are thus less likely to affect PmCRT function.

Fig. 1
figure1

Topology model of CRT from different Plasmodium spp a PmCRT. b. PfCRT and c PvCRT. The mutated amino acid residues of PmCRT are highlighted in red circle. The previous reported mutatations found in PfCRT and PvCRT are labelled in black

Genetic analysis of pmmdr1

DNA sequence analysis of P. malariae mdr1 gene (n = 95) revealed a single open reading frame of 4386 nucleotides, encoding a protein of 1461 amino acids. Results from multiple sequence alignments showed 13 single nucleotide polymorphisms. Of these, three were synonymous mutations at codons S77 (8.42%), S463 (42.10%), and L615 (31.57%). The most common polymorphism in pmmdr1 was in codon 1389 (S463), prevalent of 39.21% in the samples collected from Thailand and of 45.45% in the samples from Myanmar. The identified 10 non-synonymous mutations and their prevalence included N6I (2.10%), Y7C (5.26%), L490I (4.21%), L1063F (16.84%), N1248I (5.26%), T1266S (12.63%), R1361S (3.15%), T1406S (4.21%), A1460S (4.21%), and A1460T (2.10%) (Table 4). For comparison, the non-synonymous mutations observed in pmmdr1 gene were aligned to sequences of the orthologous pfmdr1 and pvmdr1 genes. This showed that the observed point mutations in pmmdr1 were not corresponded to the pfmdr1 or pvmdr1 mutations associated with drug resistance (Additional file 3).

Table 4 Summary of pmmdr1 gene polymorphisms in Thailand and Myanmar

The PmMDR1 sequences showed 16 haplotypes patterns (Table 5). Samples from Myanmar were classified into 13 haplotypes while the samples collected from Thailand were classified into 6 haplotypes. Fifty-five samples, 30 from Thailand and 25 from Myanmar, were wild type (haplotype 16). Haplotypes 3 (NCLLNTRTAA) and 4 (NYLLNSRTAA) were found in both countries. A total of 10 haplotypes (haplotypes 5–7 and 9–15) were only identified in Myanmar, whereas three other haplotypes (haplotypes 1, 2, and 8), were identified only in Thailand (Table 5).

Table 5 Haplotype pattern of PmMDR1 in Thailand and Myanmar

A topology model of PmMDR1 was constructed and compared to the models of PfMDR1 and PvMDR1 (Fig. 2). PmMDR1 contains two nucleotide-binding domains (NBD1 and NBD2), facing the cytoplasm. Ten PmMDR1 substitution residues were identified: N6I, Y7C, L490I, L1063F, N1248I, T1266S, R1361S, T1406S, A1460S, and A1460T. Among the 10 mutations, N6I, Y7C were located on the cytoplasmic side, L490I located in nucleotide binding domain 1, L1063F located in transmembrane domain 11, N1248I, T1266S, R1361S, T1406S, A1460S, and A1460T located in NBD2. Although the mutations observed in P. malariae were not corresponding to residues associated with drug resistance reported previously in P. falciparum, the residues located in TMD11 and NBD might have an effect on PmMDR1 function. For this, structural models for each haplotype carrying mutations in TMD11 and NBD were constructed and analysed. The structural models representing these haplotypes showed similar characteristics and the mutated residues were predicted to have only a moderate effect on the PmMDR1 structure (Fig. 3).

Fig. 2
figure2

Topology model of MDR1 from different Plasmodium spp a PmMDR1, b PfMDR1 and c PvMDR1. The mutated amino acid residues of PmMDR1 are highlight in red circle. The previous reported mutation found in PfMDR1 and PvMDR1 are labelled in black

Fig. 3
figure3

Three-dimensional structural model of PmMDR1 Transmembrane domains 1-12 are indicated as TMD1-TMD12. Two nucleotide binding domains are NBD1 and NBD2. Mutated residues in TMD11, NBD1, and NBD2 are labelled in red

Predicted changes in the physicochemical properties of each mutation were also addressed. For the L1063F mutant found in TMD11, leucine and phenylalanine share similar physicochemical properties, from which it can be inferred that this mutation has limited effect on the protein three-dimensional structure and function. A similar approach was followed for the four-point mutation identified in NBD, including L490I in NBD1 and T1266S, R1361S and A1460T/S in NBD2. To predict whether these point mutations affect the NBD structure and potentially its function, a structural model was constructed and residues involved in interactions between NBD and ATP were identified. For this, the highest conserved protein PGP1 from human was used as a template. Structural analysis of PmMDR1 model indicated that 20 and 19 residues are involved in ATP1 and ATP2 binding, respectively. Residue L490I was located outside of the ATP binding sites in NBD1. Residue T1266S, R1361S and A1460T/S in NBD2, were also located outside of the ATP binding sites.

A total of 95 P. malariae isolates were assessed for pmmdr1 gene amplification. The pmmdr1 copy number in all samples ranged from 0.75 to 1.25 (Fig. 4), which represent a single copy of the pmmdr1 gene.

Fig. 4
figure4

Distribution of pmmdr1 copy number in Thailand and Myanmar

Discussion

In this study, the anti-malarial drug target genes, pmcrt and pmmdr1, were characterized in P. malariae samples collected from Thailand and Myanmar and were compared to the orthologous mutations in pfcrt and pfmdr1 because polymorphisms in P. falciparum are well characterized for their association with chloroquine and mefloquine resistance [14,15,16, 31]. Since chloroquine and mefloquine have been used widely in the treatment of vivax and falciparum malaria in Thailand and Myanmar, it can be assumed that the P. malariae parasite population in these countries has also been exposed to these drugs, given the frequent co-infection with other human Plasmodium species in patients with P. malariae infection.

Nine nucleotide polymorphisms were identified in 5 out of the 14 exons of the pmcrt gene, with 2 non-synonymous mutations found in exon 2 (H53P) and exon 8 (E278D), both were not corresponding to the orthologous mutations in P. falciparum involved in chloroquine resistance. In P. falciparum, the majority of polymorphisms in PfCRT lay on exon 2, and mutations in PfCRT positions 72, 74, 75, and 76 are the main contributors to chloroquine resistance [32, 33]. The in silico topology model of PmCRT revealed that H53P and E278D are located in cytoplasmic region and food vacuole, which are outside of the transmembrane domain, and are thus considered unlikely to affect PmCRT function. Genetic polymorphisms in the pmcrt gene in both Thailand and Myanmar samples were limited, similar to previous studies describing very limited polymorphisms in the pvcrt gene in P. vivax populations in Thailand [34].

The number of polymorphisms in this study was higher in the pmmdr1 gene. There were 13 SNPs, including 10 non-synonymous and 3 synonymous mutations, and none of these showed an equivalent position in pfmdr1 associated with anti-malarial drug resistance in P. falciparum. Some of the P. malariae studied here were mixed infections with other species, which might have an impact on polymorphisms of pmcrt and pmmdr1. There were 12 and 11 mixed infections found in the samples from Thailand and Myanmar, respectively. The patterns of nucleotide polymorphisms in pmcrt and pmmdr1 between single and mixed infections were compared. The proportion of mixed infected samples carrying mutations in Thailand and Myanmar was accounted for 14.29-22.72% (Additional file 4), suggesting that the mixed infections are unlikely to affect pmcrt and pmmdr1 mutations in this study.

The mutations in pmmdr1 codons L490I and L1063F correspond to the previously reported mutations in P. vivax at pvmdr1 codons L493L and F1076L [17, 35]. In P. vivax, these polymorphisms have not been clearly linked to chloroquine resistance [36, 37]. The observed polymorphisms in pmmdr1 were translated into a topology model of PmMDR1, showing that the mutations resulted in predicted protein changes located in the parasite cytoplasmic side of the protein, in nucleotide binding domains and one change in TMD11. In P. falciparum, residues in TMD11 were suggested as part of an anti-malarial binding pocket [38,39,40]. An in silico homology model of PmMDR1 was constructed and analysed to predict the effect of amino acid changes in TMD11 and NBD, which showed that the observed mutations were unlikely to affect the tertiary structure of the protein. Additionally, structural analysis of PmMDR1 polymorphisms found in NBD1 and NBD2 suggested that the mutations (L490I and L1063F) had only a moderate effect on the conformation of these domains, unlikely to affect the NBD function. Overall, the described polymorphisms were predicted to have insignificant impact on PmMDR1 protein morphology and function.

The samples used in this study were from Thailand and Myanmar collected at different periods of time, which may potentially affect the pattern of gene mutations [34]. The samples from Myanmar were collected in year 2009 in which artemether-lumefantrine was used as the first line treatment for uncomplicated malaria [41]. For Thailand, P. malariae were collected from two periods of time, year 2002–2008 and year 2012–2016. During those time periods, artesunate-mefloquine was used as the first-line treatment before it was changed to dihydroartemisinin-piperaquine in 2015 [41]. The pmcrt in Thailand showed more mutations (haplotypes 1-3, 81.48-91.67%) when compared to Myanmar (31.82%) (Additional file 5). There was no difference in mutation pattern of pmmdr1 in both countries (Additional file 5). Although H53P and E278D mutations found in PmCRT were predicted that they were unlikely to have an impact on PmCRT function, the high proportion of those point mutations in Thailand might refer to geographical characteristics of the parasites. This will need to be confirmed in a larger sample size collected from different areas.

In P. falciparum, amplification of the pfmdr1 gene is strongly associated with mefloquine resistance [14,15,16]. Despite the presumed long-term mefloquine drug pressure on the P. malariae parasite population in the study areas, none of P. malariae samples carried amplification of pmmdr1. Possible explanations include firstly, that pmmdr1 amplification is not involved in mefloquine resistance in P. malariae and thus is not a good marker, and that there might be an alternative mechanism conferring mefloquine resistance in P. malariae other than mdr1 amplification. Secondly, the parasite loads during infection of P. malariae is rarely exceed 1000 parasites per µl of blood. Thus, the number of parasites under selective pressure would be low for P. malariae and the likelihood of selecting resistant is lower. However, some of the potential drug resistance markers that have been studied in P. malariae might be under selective pressure such as pmdhfr [42, 43], pmdhps [44], and pmkelch [45]. Thirdly, although clinical P. malariae infection often presents as co-infection with other human Plasmodium species, there might be a large P. malariae reservoir outside of these patients with co-infection or a reservoir in non-human primates, so that overall P. malariae population has been limited to the selective pressure. The study of chloroquine-resistant P. vivax revealed that increased expression of pvcrt and pvmdr1 are associated with chloroquine resistance [46, 47]. Moreover, the gene copy number of pvcrt was significantly higher in chloroquine-resistant P. vivax [48]. In addition to this study, expression level of pmcrt and pmmdr1, and the copy number variation of pmcrt should be evaluated.

Conclusions

Polymorphisms in pmmdr1 were more frequently observed than in pmcrt. The non-synonymous mutations found in both pmcrt and pmmdr1 were unlikely to affect protein function. No amplification of pmmdr1 was observed in this study. If the orthologous resistance genes in P. malariae are indeed associated with anti-malarial drug resistance in this Plasmodium species, the findings suggest limited chloroquine and mefloquine drug pressure on the P. malariae populations in the study regions. Alternatively, anti-malarial drug resistance in P. malariae could differ from that described in P. falciparum and P. vivax, which will require further investigation.

Availability of data and materials

The nucleotide sequences of pmcrt and pmmdr1 genes obtained from this study have been submitted to GenBank database under the Accession Numbers; MN623294-MN623342 (for pmcrt) and MN645870-MN645902 (for pmmdr1).

Abbreviations

pmcrt :

Plasmodium malariae chloroquine resistant transporter

pfcrt :

Plasmodium falciparum chloroquine resistant transporter

pvcrt :

Plasmodium vivax chloroquine resistant transporter

pmmdr1 :

Plasmodium malariae multidrug resistance 1

pfmdr1 :

Plasmodium falciparum multidrug resistance 1

pvmdr1 :

Plasmodium vivax multidrug resistance 1

SNPs:

Single nucleotide polymorphisms

TMD:

Transmembrane domain

NBD:

Nucleotide binding domain

References

  1. 1.

    Ansari HR, Templeton TJ, Subudhi AK, Ramaprasad A, Tang JX, Lu F, et al. Genome-scale comparison of expanded gene families in Plasmodium ovale wallikeri and Plasmodium ovale curtisi with Plasmodium malariae and with other Plasmodium species. Int J Parasitol. 2016;46:685–96.

    CAS  PubMed  Google Scholar 

  2. 2.

    Rutledge GG, Bohme U, Sanders M, Rrid AJ, Cotton JA, Maiga-Ascofare O, et al. Plasmodium malariae and P. ovale genomes provide insights into malaria parasite evolution. Nature. 2017;542:101–4.

    CAS  PubMed  PubMed Central  Google Scholar 

  3. 3.

    Sutherland CJ. Persistent parasitism: the adaptive biology of malariae and ovale malaria. Trends Parasitol. 2016;32:808–19.

    PubMed  Google Scholar 

  4. 4.

    Tan-ariya PPS. First report of in vitro susceptibility of Plasmodium malariae Thai isolates to chloroquine. Southeast Asian J Trop Med Public Health. 1994;25:784–7.

    CAS  PubMed  Google Scholar 

  5. 5.

    Ringwald PBJ, Same-Ekobo A, Basco LK. Pyronaridine for treatment of Plasmodium ovale and Plasmodium malariae infections. Antimicrob Agents Chemother. 1997;41:2317–9.

    CAS  PubMed  PubMed Central  Google Scholar 

  6. 6.

    Borrmann S, Szlezak N, Binder RK, Missinou MA, Lell B, Kremsner PG. Evidence for the efficacy of artesunate in asymptomatic Plasmodium malariae infections. J Antimicrob Chemoth. 2002;50:751–4.

    CAS  Google Scholar 

  7. 7.

    Maguire JD, Sumawinata IW, Masbar S, Laksana B, Prodjodipuro P, Susanti I, et al. Chloroquine-resistant Plasmodium malariae in south Sumatra Indonesia. Lancet. 2002;360:58–60.

    CAS  PubMed  Google Scholar 

  8. 8.

    Collins WE, Jeffery GM. Extended clearance time after treatment of infections with Plasmodium malariae may not be indicative of resistance to chloroquine. Am J Trop Med Hyg. 2002;67:406–10.

    PubMed  Google Scholar 

  9. 9.

    Noisang C, Prosser C, Meyer W, Chemoh W, Ellis J, Sawangjaroen N, et al. Molecular detection of drug resistant malaria in Southern Thailand. Malar J. 2019;18:275.

    PubMed  PubMed Central  Google Scholar 

  10. 10.

    Setthaudom C, Tan-ariya P, Sitthichot N, Khositnithikul R, Suwandittakul N, Leelayoova S, et al. Role of Plasmodium falciparum chloroquine resistance transporter and multidrug resistance 1 genes on in vitro chloroquine resistance in isolates of Plasmodium falciparum from Thailand. Am J Trop Med Hyg. 2011;85:606–11.

    CAS  PubMed  PubMed Central  Google Scholar 

  11. 11.

    Li JY, Zhang J, Li Q, Hu Y, Ruan YH, Tao ZY, et al. Ex vivo susceptibilities of Plasmodium vivax isolates from the China-Myanmar border to antimalarial drugs and association with polymorphisms in Pvmdr1 and Pvcrt-o genes. PLoS Negl Trop Dis. 2020;14:e0008255.

    CAS  PubMed  PubMed Central  Google Scholar 

  12. 12.

    Zhao Y, Liu ZL, Soe MT, Wang L, Soe TN, Wei HP, et al. Genetic variations associated with drug resistance markers in asymptomatic Plasmodium falciparum infections in Myanmar. Genes. 2019;10:692.

    CAS  PubMed Central  Google Scholar 

  13. 13.

    Packard RM. The origins of antimalarial-drug resistance. N Engl J Med. 2014;371:397–9.

    CAS  PubMed  Google Scholar 

  14. 14.

    Cowman AFGD, Thompson JK. Selection for mefloquine resistance in Plasmodium falciparum is linked to amplification of the pfmdr1 gene and cross-resistance to halofantrine and quinine. Proc Natl Acad Sci USA. 1994;91:1143–7.

    CAS  PubMed  Google Scholar 

  15. 15.

    Preechapornkul P, Imwong M, Chotivanich K, Pongtavornpinyo W, Dondorp AM, Day NPJ, et al. Plasmodium falciparum pfmdr1 amplification, mefloquine resistance, and parasite fitness. Antimicrob Agents Chemother. 2009;53:1509–15.

    CAS  PubMed  PubMed Central  Google Scholar 

  16. 16.

    Price RN, Uhlemann AC, Brockman A, McGready R, Ashley E, Phaipun L, et al. Mefloquine resistance in Plasmodium falciparum and increased pfmdr1 gene copy number. Lancet. 2004;364:438–47.

    CAS  PubMed  PubMed Central  Google Scholar 

  17. 17.

    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.

    CAS  PubMed  PubMed Central  Google Scholar 

  18. 18.

    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.

    CAS  PubMed  PubMed Central  Google Scholar 

  19. 19.

    Snounou GPL, Gonçalves A, Fonseca L, Dias F, Brown KN, do Rosario VE. The importance of sensitive detection of malaria parasites in the human and insect hosts in epidemiological studies, as shown by the analysis of field samples from Guinea Bissau. Trans R Soc Trop Med Hyg. 1993;87:649–53.

    CAS  PubMed  Google Scholar 

  20. 20.

    Snounou G, Singh B. Nested PCR analysis of Plasmodium parasites. Methods Mol Med. 2002;72:189–203.

    CAS  PubMed  Google Scholar 

  21. 21.

    Hall TA. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucl Acids Symp Ser. 1999;41:95–8.

    CAS  Google Scholar 

  22. 22.

    Rozas J, Ferrer-Mata A, Sanchez-DelBarrio JC, Guirao-Rico S, Librado P, Ramos-Onsins SE, et al. DnaSP 6: DNA sequence polymorphism analysis of large data sets. Mol Biol Evol. 2017;34:3299–302.

    CAS  PubMed  Google Scholar 

  23. 23.

    Kumar S, Stecher G, Tamura K. Molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016;33:1870–4.

    CAS  PubMed  PubMed Central  Google Scholar 

  24. 24.

    Kelley LA, Mezulis S, Yates CM, Wass MN, Sternberg MJE. The Phyre2 web portal for protein modeling, prediction and analysis. Nat Protoc. 2015;10:845–58.

    CAS  PubMed  PubMed Central  Google Scholar 

  25. 25.

    Waterhouse A, Bertoni M, Bienert S, Studer G, Tauriello G, Gumienny R, et al. SWISS-MODEL: homology modelling of protein structures and complexes. Nucleic Acids Res. 2018;46:W296-03.

    Google Scholar 

  26. 26.

    Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997;25:3389–402.

    CAS  PubMed  PubMed Central  Google Scholar 

  27. 27.

    McGuffin LJ, Bryson K, Jones DT. The PSIPRED protein structure prediction server. Bioinformatics. 2000;16:404–5.

    CAS  PubMed  Google Scholar 

  28. 28.

    Ward JJ, McGuffin LJ, Bryson K, Buxton BF, Jones DT. The DISOPRED server for the prediction of protein disorder. Bioinformatics. 2004;20:2138–9.

    CAS  PubMed  Google Scholar 

  29. 29.

    Nugent T, Jones DT. Predicting transmembrane helix packing arrangements using residue contacts and a force-directed algorithm. PLoS Comput Biol. 2010;6:e1000714.

    PubMed  PubMed Central  Google Scholar 

  30. 30.

    BIOVIA. Discovery StudioVisualizer. 2017.

  31. 31.

    Happi CT, Gbotosho GO, Folarin OA, Bolaji OM, Sowunmi A, Kyle DE, et al. Association between mutations in Plasmodium falciparum chloroquine resistance transporter and P. falciparum multidrug resistance 1 genes and in vivo amodiaquine resistance in P. falciparum malaria-infected children in Nigeria. Am J Trop Med Hyg. 2006;75:155–61.

    CAS  PubMed  Google Scholar 

  32. 32.

    Fidock DA, Nomura T, Talley AK, Cooper RA, Dzekunov SM, Ferdig MT, et al. Mutations in the P. falciparum digestive vacuole transmembrane protein PfCRT and evidence for their role in chloroquine resistance. Mol Cell. 2000;6:861–71.

    CAS  PubMed  PubMed Central  Google Scholar 

  33. 33.

    Wootton JC, Feng XR, Ferdig MT, Cooper RA, Mu JB, Baruch DI, et al. Genetic diversity and chloroquine selective sweeps in Plasmodium falciparum. Nature. 2002;418:320–3.

    CAS  PubMed  Google Scholar 

  34. 34.

    Tantiamornkul K, Pumpaibool T, Piriyapongsa J, Culleton R, Lek-Uthai U. The prevalence of molecular markers of drug resistance in Plasmodium vivax from the border regions of Thailand in 2008 and 2014. Int J Parasitol Drugs Drug Resist. 2018;8:229–37.

    PubMed  PubMed Central  Google Scholar 

  35. 35.

    Kittichai V, Nguitragool W, Mbenda HGN, Sattabongkot J, Cui LW. Genetic diversity of the Plasmodium vivax multidrug resistance 1 gene in Thai parasite populations. Infect Genet Evol. 2018;64:168–77.

    CAS  PubMed  PubMed Central  Google Scholar 

  36. 36.

    Nomura T, Carlton JMR, 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.

    CAS  PubMed  Google Scholar 

  37. 37.

    Brega S, Meslin B, de Monbrison R, Severini C, Gradoni L, Udomsangpetch R, et al. 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–7.

    CAS  PubMed  Google Scholar 

  38. 38.

    Ferreira PE, Holmgren G, Veiga MI, Uhlen P, Kaneko A, Gil JP. PfMDR1: Mechanisms of transport modulation by functional polymorphisms. PLoS ONE. 2011;6:e23875.

    CAS  PubMed  PubMed Central  Google Scholar 

  39. 39.

    Lekostaj JK, Amoah LE, Roepe PD. A single S1034C mutation confers altered drug sensitivity to PfMDR1 ATPase activity that is characteristic of the 7G8 isoform. Mol Biochem Parasit. 2008;157:107–11.

    CAS  Google Scholar 

  40. 40.

    Sidhu ABS, Valderramos SG, Fidock DA. pfmdr1 mutations contribute to quinine resistance and enhance mefloquine and artemisinin sensitivity in Plasmodium falciparum. Mol Microbiol. 2005;57:913–26.

    CAS  PubMed  Google Scholar 

  41. 41.

    WHO. World malaria report 2016. Geneva: World Health Organization; 2016.

  42. 42.

    Khim N, Kim S, Bouchier C, Tichit M, Ariey F, Fandeur T, et al. Reduced impact of pyrimethamine drug pressure on Plasmodium malariae dihydrofolate reductase gene. Antimicrob Agents Chemother. 2012;56:863–8.

    CAS  PubMed  PubMed Central  Google Scholar 

  43. 43.

    Tanomsing N, Imwong M, Pukrittayakamee S, Chotivanich K, Looareesuwan S, Mayxay M, et al. Genetic analysis of the dihydrofolate reductase-thymidylate synthase gene from geographically diverse isolates of Plasmodium malariae. Antimicrob Agents Chemother. 2007;51:3523–30.

    CAS  PubMed  PubMed Central  Google Scholar 

  44. 44.

    Tanomsing N, Mayxay M, Newton PN, Nosten F, Dolecek C, Hien TT, et al. Genetic Variability of Plasmodium malariae dihydropteroate synthase (dhps) in Four Asian Countries. PLoS ONE. 2014;9:e93942.

    PubMed  PubMed Central  Google Scholar 

  45. 45.

    Nakeesathit S, Saralamba N, Pukrittayakamee S, Dondorp A, Nosten F, White NJ, et al. Limited polymorphism of the kelch propeller domain in Plasmodium malariae and P. ovale isolates from Thailand. Antimicrob Agents Chemother. 2016;60:4055–62.

    CAS  PubMed  PubMed Central  Google Scholar 

  46. 46.

    Fernandez-Becerra C, Pinazo MJ, Gonzalez A, Alonso PL, del Portillo HA, Gascon J. Increased expression levels of the pvcrt-o and pvmdr1 genes in a patient with severe Plasmodium vivax malaria. Malar J. 2009;8:55.

    PubMed  PubMed Central  Google Scholar 

  47. 47.

    Melo GC, Monteiro WM, Siqueira AM, Silva SR, Magalhaes BML, Alencar ACC, 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.

    PubMed  PubMed Central  Google Scholar 

  48. 48.

    Silva SR, Almeida ACG, da Silva GAV, Ramasawmy R, Lopes SCP, Siqueira AM, et al. Chloroquine resistance is associated to multi-copy pvcrt-o gene in Plasmodium vivax malaria in the Brazilian Amazon. Malar J. 2018;17:267.

    PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

The authors would like to thank all staff involving in sample collection across all study sites in Myanmar and Thailand.

Funding

This study was financed in part by Mahidol University, Thailand, and was part of the Wellcome Trust Mahidol University-Oxford Tropical Medicine Research Programme supported by the Wellcome Trust of Great Britain, and Thailand Science Research and Innovation (TSRI), RTA6280006. NS was supported by Research Fund 2017, Faculty of Tropical Medicine, Mahidol University.

Author information

Affiliations

Authors

Contributions

YP, NS, AD, and MI contributed to study design. SP, FN, and FS collected samples. YP and NS carried out laboratory work. YP, NS, UB, ND, AD and MI analysed the data. YP, NS, UB, AD and MI drafted the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Naowarat Saralamba.

Ethics declarations

Ethics approval and consent to participate

Ethical approval was obtained and reviewed by the Ethics Committee from Faculty of Tropical Medicine, Mahidol University, Bangkok, Thailand (MUTM 2018-013-01).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

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

Supplementary information

Additional file 1.

The PCR primers and condition for amplification of pmcrt gene.

Additional file 2.

The PCR primers and condition for amplification of pmmdr1 gene.

Additional file 3.

The equivalence position of pmmdr1 gene compared with pfmdr1 and pvmdr1.

Additional file 4.

Summary of point mutations in pmcrt and pmmdr1 and mixed P. malariae infections.

Additional file 5.

Summary of point mutations in pmcrt and pmmdr1 and sampling times.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. 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 in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Pimpat, Y., Saralamba, N., Boonyuen, U. et al. Genetic analysis of the orthologous crt and mdr1 genes in Plasmodium malariae from Thailand and Myanmar. Malar J 19, 315 (2020). https://doi.org/10.1186/s12936-020-03391-6

Download citation

Keywords

  • Malaria
  • Plasmodium malariae
  • Chloroquine resistant transporter
  • Multidrug resistance proteins 1