Geographic distribution of amino acid mutations in DHFR and DHPS in Plasmodium vivax isolates from Lao PDR, India and Colombia

Genomics of Plasmodium vivax in Colombia reveals evidence of local bottle-necking and inter-country connectivity in the Americas

Colombia aims to eliminate malaria by 2030 but remains one of the highest burden countries in the Americas. Plasmodium vivax contributes half of all malaria cases, with its control challenged by relapsing parasitaemia, drug resistance and cross-border spread. Using 64 Colombian P. vivax genomes collected between 2013 and 2017, we explored diversity and selection in two major foci of transmission: Chocó and Córdoba. Open-access data from other countries were used for comparative assessment of drug resistance candidates and to assess cross-border spread. Across Colombia, polyclonal infections were infrequent (12%), and infection connectivity was relatively high (median IBD = 5%), consistent with low endemicity. Chocó exhibited a higher frequency of polyclonal infections (23%) than Córdoba (7%), although the difference was not significant (P = 0.300). Most Colombian infections carried double pvdhfr (95%) and single pvdhps (71%) mutants, but other drug resistance mutations were less prevalent (< 10%). There was no evidence of selection at the pvaat1 gene, whose P. falciparum orthologue has recently been implicated in chloroquine resistance. Global population comparisons identified other putative adaptations. Within the Americas, low-level connectivity was observed between Colombia and Peru, highlighting potential for cross-border spread. Our findings demonstrate the potential of molecular data to inform on infection spread and adaptation.

Global scenario of Plasmodium vivax occurrence and resistance pattern

Malaria caused by Plasmodium vivax is comparatively less virulent than Plasmodium falciparum, which can also lead to severe disease and death. It shows a wide geographical distribution. Chloroquine serves as a drug of choice, with primaquine as a radical cure. However, with the appearance of resistance to chloroquine and treatment has been shifted to artemisinin combination therapy followed by primaquine as a radical cure. Sulphadoxine-pyrimethamine, mefloquine, and atovaquone-proguanil are other drugs of choice in chloroquine-resistant areas, and later resistance was soon reported for these drugs also. The emergence of drug resistance serves as a major hurdle to controlling and eliminating malaria. The discovery of robust molecular markers and regular surveillance for the presence of mutations in malaria-endemic areas would serve as a helpful tool to combat drug resistance. Here, in this review, we will discuss the endemicity of P. vivax, a historical overview of antimalarial drugs, the appearance of drug resistance and molecular markers with their global distribution along with different measures taken to reduce malaria burden due to P. vivax infection and their resistance.

Geographical spread and structural basis of sulfadoxine-pyrimethamine drug-resistant malaria parasites

The global spread of sulfadoxine (Sdx, S) and pyrimethamine (Pyr, P) resistance is attributed to increasing number of mutations in DHPS and DHFR enzymes encoded by malaria parasites. The association between drug resistance mutations and SP efficacy is complex. Here we provide an overview of the geographical spread of SP resistance mutations in Plasmodium falciparum (Pf) and Plasmodium vivax (Pv) encoded dhps and dhfr genes. In addition, we have collated the mutation data and mapped it on to the three-dimensional structures of DHPS and DHFR which have become available. Data from genomic databases and 286 studies were collated to provide a comprehensive landscape of mutational data from 2005 to 2019. Our analyses show that the Pyr-resistant double mutations are widespread in Pf/PvDHFR (P. falciparum ∼61% in Asia and the Middle East, and in the Indian sub-continent; in P. vivax ∼33% globally) with triple mutations prevailing in Africa (∼66%) and South America (∼33%). For PfDHPS, triple mutations dominate South America (∼44%), Asia and the Middle East (∼34%) and the Indian sub-continent (∼27%), while single mutations are widespread in Africa (∼45%). Contrary to the status for P. falciparum, Sdx-resistant single point mutations in PvDHPS dominate globally. Alarmingly, highly resistant quintuple and sextuple mutations are rising in Africa (PfDHFR-DHPS) and Asia (Pf/PvDHFR-DHPS). Structural analyses of DHFR and DHPS proteins in complexes with substrates/drugs have revealed that resistance mutations map proximal to Sdx and Pyr binding sites. Thus new studies can focus on discovery of novel inhibitors that target the non-substrate binding grooves in these two validated malaria parasite drug targets.

Distribution pattern of amino acid mutations in chloroquine and antifolate drug resistance associated genes in complicated and uncomplicated Plasmodium vivax isolates from Chandigarh, North India

BACKGROUND

The increasing antimalarial drug resistance is a significant hindrance to malaria control and elimination programs. For the last six decades, chloroquine (CQ) plus pyrimethamine remains the first-line treatment for P. vivax malaria. Regions where both P. falciparum and P. vivax co-exist, P. vivax is exposed to antifolate drugs due to either misdiagnosis or improper treatment that causes selective drug pressure to evolve. Therefore, the present study aims to estimate antimalarial drug resistance among the complicated and uncomplicated P. vivax patients.

METHODS

A total of 143 P. vivax malaria positive patients were enrolled in this study, and DNA was isolated from their blood samples. Pvcrt-o, Pvmdr-1, Pvdhps, and Pvdhfr genes were PCRs amplified, and drug resistance-associated gene mutations were analyzed. Statistical analysis of the drug resistance genes and population diversity was performed using MEGA vs. 7.0.21 and DnaSP v software.

RESULTS

Among the CQ resistance marker gene Pvcrt-o, the prevalence of K10 insertion was 17.5% (7/40) and 9.5% (7/73) of complicated and uncomplicated P vivax group isolates respectively. In Pvmdr-1, double mutant haplotype (M₉₅₈/L₁₀₇₆) was found in 99% of the clinical isolates. Among the pyrimethamine resistance-associated gene Pvdhfr, the double mutant haplotype I₁₃P₃₃F₅₇R₅₈T₆₁N₁₁₇I₁₇₃ was detected in 23% (11/48) in complicated and 20% (17/85) in uncomplicated group isolates. In the sulphadoxine resistance-associated Pvdhps gene, limited polymorphism was observed with the presence of a single mutant (D459A) among 16 and 5% of the clinical isolates in the complicated and uncomplicated group respectively.

CONCLUSION

The study presents the situations of polymorphism in the antimalarial drug resistance-associated genes and emphasizes the need for regular surveillance. It is imperative for the development of suitable antimalarial drug policy in India.

Polymorphisms in Plasmodium vivax antifolate resistance markers in Afghanistan between 2007 and 2017

BACKGROUND

Plasmodium vivax is the predominant Plasmodium species in Afghanistan. National guidelines recommend the combination of chloroquine and primaquine (CQ-PQ) for radical treatment of P. vivax malaria. Artesunate in combination with the antifolates sulfadoxine-pyrimethamine (SP) has been first-line treatment for uncomplicated falciparum malaria until 2016. Although SP has been the recommended treatment for falciparum and not vivax malaria, exposure of the P. vivax parasite population to SP might still have been quite extensive because of community based management of malaria. The change in the P. vivax antifolate resistance markers between 2007 and 2017 were investigated.

METHODS

Dried blood spots were collected (n = 185) from confirmed P. vivax patients in five malaria-endemic areas of Afghanistan bordering Tajikistan, Turkmenistan and Pakistan, including Takhar, Faryab, Laghman, Nangarhar, and Kunar, in 2007, 2010 and 2017. Semi-nested PCR, RFLP and nucleotide sequencing were used to assess the pyrimethamine resistant related mutations in P. vivax dihydrofolate reductase (pvdhfr I13L, P33L, N50I, F57L, S58R, T61I, S93H, S117N, I173L) and the sulfonamide resistance related mutations in P. vivax dihydropteroate synthase (pvdhps A383G, A553G).

RESULTS

In the 185 samples genotyped for pvdhfr and pvdhps mutations, 11 distinct haplotypes were observed, which evolved over time. In 2007, wild type pvdhfr and pvdhps were the most frequent haplotype in all study sites (81%, 80/99). However, in 2017, the frequency of the wild-type was reduced to 36%, (21/58; p value ≤ 0.001), with an increase in frequency of the double mutant pvdhfr and pvdhps haplotype S58RS117N (21%, 12/58), and the single pvdhfr mutant haplotype S117N (14%, 8/58). Triple and quadruple mutations were not found. In addition, pvdhfr mutations at position N50I (7%, 13/185) and the novel mutation S93H (6%, 11/185) were observed. Based on in silico protein modelling and molecular docking, the pvdhfr N50I mutation is expected to affect only moderately pyrimethamine binding, whereas the S93H mutation does not.

CONCLUSIONS

In the course of ten years, there has been a strong increase in the frequency pyrimethamine resistance related mutations in pvdhfr in the P. vivax population in Afghanistan, although triple and quadruple mutations conferring high grade resistance were not observed. This suggests relatively low drug pressure from SP on the P. vivax parasite population in the study areas. The impact of two newly identified mutations in the pvdhfr gene on pyrimethamine resistance needs further investigation.

Plasmodium vivax liver stage assay platforms using Indian clinical isolates

Background

Vivax malaria is associated with significant morbidity and economic loss, and constitutes the bulk of malaria cases in large parts of Asia and South America as well as recent case reports in Africa. The widespread prevalence of vivax is a challenge to global malaria elimination programmes. Vivax malaria control is particularly challenged by existence of dormant liver stage forms that are difficult to treat and are responsible for multiple relapses, growing drug resistance to the asexual blood stages and host-genetic factors that preclude use of specific drugs like primaquine capable of targeting Plasmodium vivax liver stages. Despite an obligatory liver-stage in the Plasmodium life cycle, both the difficulty in obtaining P. vivax sporozoites and the limited availability of robust host cell models permissive to P. vivax infection are responsible for the limited knowledge of hypnozoite formation biology and relapse mechanisms, as well as the limited capability to do drug screening. Although India accounts for about half of vivax malaria cases world-wide, very little is known about the vivax liver stage forms in the context of Indian clinical isolates.

Methods

To address this, methods were established to obtain infective P. vivax sporozoites from an endemic region in India and multiple assay platforms set up to detect and characterize vivax liver stage forms. Different hepatoma cell lines, including the widely used HCO4 cells, primary human hepatocytes as well as hepatocytes obtained from iPSC’s generated from vivax patients and healthy donors were tested for infectivity with P. vivax sporozoites.

Results

Both large and small forms of vivax liver stage are detected in these assays, although the infectivity obtained in these platforms are low.

Conclusions

This study provides a proof of concept for detecting liver stage P. vivax and provide the first characterization of P. vivax liver stage forms from an endemic region in India.

Lack of quadruple and quintuple mutant alleles associated with sulfadoxine-pyrimethamine resistance in Plasmodium vivax isolates from Brazilian endemic areas

BACKGROUND AND OBJECTIVE

Brazil is responsible for a large number of Plasmodium vivax cases in America. Given the emergence of P. vivax parasites resistant to chloroquine and the effectiveness of antifolates in vivax malaria treatment together with a correlation between mutations in P. vivax dhfr and dhps genes and SP treatment failure, the point mutations in these genes were investigated.

METHODS

Blood samples from 54 patients experiencing vivax malaria symptomatic episodes in the Amazonian Region were investigated. Genomic DNA was extracted using a DNA extraction kit (QIAGENTM). Nested polymerase chain reaction (PCR) amplification was carried out followed by Sanger sequencing to detect single nucleotide polymorphisms (SNPs).

FINDINGS

All tested isolates showed non-synonymous mutations in pvdhfr gene: 117N (54/54, 100%) and 58R (25/54, 46%). Double mutant allele 58R/117N (FRTNI, 28%) was the most frequent followed by triple mutant alleles (58R/117N/173L, FRTNL, 11%; 58R/61M/117N, FRMNI, 5% 117N/173L, FSTNL, 4%) and quadruple mutant allele (58R/61M/117N/173L, FRMNL, 2%). A single mutation was observed at codon C383G in pvdhps gene (SGKAV, 48%).

CONCLUSION

No evidence of molecular signatures associated with P. vivax resistance to SP was observed in the Brazilian samples.

Molecular monitoring of dihydrofolatereductase (dhfr) and dihydropteroatesynthetase (dhps) associated with sulfadoxine-pyrimethamine resistance in Plasmodium vivax isolates of Palawan, Philippines

The emergence of drug-resistant Plasmodium vivax poses problems for malaria control and elimination in some parts of the world, especially in developing countries where individuals are routinely exposed to the infection. The aim of this study was to determine the single nucleotide polymorphisms (SNPs) in dihydropteroate synthase (pvdhps) and dihydrofolate reductase (pvdhfr) genes associated with sulfadoxine-pyrimethamine (SP) drug resistance among P. vivax isolates collected in Palawan, Philippines. Genetic polymorphisms of pvdhps and pvdhfr were analysed by nested PCR. Analysis at specific codons I₁₃P₃₃F₅₇S₅₈T₆₁S₁₁₇I₁₇₃ associated with pyrimethamine resistance in the pvdhfr gene revealed that most of the samples (66/87, 75.9%) carried double mutation at positions I₁₃P₃₃F₅₇R₅₈T₆₁N₁₁₇I₁₇₃, while only 18.4% (16/87) of the isolates carried the wild-type haplotype (I₁₃P₃₃F₅₇S₅₈T₆₁S₁₁₇I₁₇₃). For the pvdhps gene, the codons involved in sulfadoxine resistance S₃₈₂A₃₈₃K₅₁₂A₅₅₃V₅₈₅ were investigated. Single mutation at position S₃₈₂G₃₈₃K₅₁₂A₅₅₃V₅₈₅ was most observed in 68.0% (68/100) of the samples, whereas wild-type haplotype was found in 26.0% (26/100) of samples. The pvdhps and pvdhfr combination S₃₈₂A₃₈₃K₅₁₂A₅₅₃V₅₈₅/I₁₃P₃₃F₅₇S₅₈T₆₁S₁₁₇I₁₇₃ (wild-type), S₃₈₂G₃₈₃K₅₁₂A₅₅₃V₅₈₅/I₁₃P₃₃F₅₇R₅₈T₆₁N₁₁₇I₁₇₃, and S₃₈₂A₃₈₃K₅₁₂A₅₅₃V₅₈₅-I₁₃P₃₃F₅₇R₅₈T₆₁N₁₁₇I₁₇₃ were the most frequently observed combination haplotypes from the three study sites. The information on molecular markers associated with antifolate drug-resistance could help better understanding ofthe molecular epidemiology and situation of SP resistant P. vivax malaria in the country. Continuous surveillance of these genetic markers is necessary to monitor the evolution of SP resistance in the Philippines.

Simultaneous detection of Plasmodium vivax dhfr, dhps, mdr1 and crt-o resistance-associated mutations in the Colombian Amazonian region

Background

Malaria continues being a public health problem worldwide. Plasmodium vivax is the species causing the largest number of cases of malaria in Asia and South America. Due to the lack of a completely effective anti-malarial vaccine, controlling this disease has been based on transmission vector management, rapid diagnosis and suitable treatment. However, parasite resistance to anti-malarial drugs has become a major yet-to-be-overcome challenge. This study was thus aimed at determining pvmdr1, pvdhfr, pvdhps and pvcrt-o gene mutations and haplotypes from field samples obtained from an endemic area in the Colombian Amazonian region.

Methods

Fifty samples of parasite DNA infected by a single P. vivax strain from symptomatic patients from the Amazonas department in Colombia were analysed by PCR and the pvdhfr, pvdhps, pvmdr1 and pvcrt-o genes were sequenced. Diversity estimators were calculated from the sequences and the haplotypes circulating in the Colombian Amazonian region were obtained.

Conclusion

pvdhfr, pvdhps, pvmdr1 and pvcrt-o genes in the Colombian Amazonian region are characterized by low genetic diversity. Some resistance-associated mutations were found circulating in this population. New variants are also being reported. A selective sweep signal was located in pvdhfr and pvmdr1 genes, suggesting that these mutations (or some of them) could be providing an adaptive advantage.

Distribution of Mutations Associated with Antifolate and Chloroquine Resistance among Imported Plasmodium vivax in the State of Qatar

Plasmodium vivax is the most prevalent parasite worldwide, escalating by spread of drug resistance. Currently, in Qatar, chloroquine (CQ) plus primaquine are recommended for the treatment of P. vivax malaria. The present study examined the prevalence of mutations in dihydrofolate reductase (dhfr), dihydropteroate synthase (dhps) genes and CQ resistance transporter (crt-o) genes, associated with sulphadoxine-pyrimethamine (SP) and chloroquine resistance, among imported P. vivax cases in Qatar. Blood samples were collected from patients positive for P. vivax and seeking medical treatment at Hamad General Hospital, Doha, during 2013-2016. The Sanger sequencing method was performed to examine the single nucleotide polymorphisms in Pvdhfr, Pvdhps, and Pvcrt-o genes. Of 314 examined P. vivax isolates, 247 (78.7%), 294 (93.6%) and 261 (83.1%) were successfully amplified and sequenced for Pvdhfr, Pvdhps, and Pvcrt-o, respectively. Overall, 53.8% (N = 133) carried mutant alleles (58R/117N) in Pvdhfr, whereas 77.2% (N = 227) and 90% (N = 235) isolates possessed wild type allele in Pvdhps and Pvcrt-o genes, respectively. In addition, a total of eleven distinct haplotypes were detected in Pvdhfr/Pvdhps genes. Interestingly, K10 insertion in the Pvcrt-o gene was observed only in patients originating from the Indian subcontinent. The results suggested that CQ remains an acceptable treatment regimen but further clinical data are required to assess the effectiveness of CQ and SP in Qatar to support the current national treatment guidelines. In addition, limited distribution of genetic polymorphisms associated with CQ and SP resistance observed in imported P. vivax infections, necessitates regular monitoring of drug resistant P. vivax malaria in Qatar.

Genomic variation in Plasmodium vivax malaria reveals regions under selective pressure

BACKGROUND

Although Plasmodium vivax contributes to almost half of all malaria cases outside Africa, it has been relatively neglected compared to the more deadly P. falciparum. It is known that P. vivax populations possess high genetic diversity, differing geographically potentially due to different vector species, host genetics and environmental factors.

RESULTS

We analysed the high-quality genomic data for 46 P. vivax isolates spanning 10 countries across 4 continents. Using population genetic methods we identified hotspots of selection pressure, including the previously reported MRP1 and DHPS genes, both putative drug resistance loci. Extra copies and deletions in the promoter region of another drug resistance candidate, MDR1 gene, and duplications in the Duffy binding protein gene (PvDBP) potentially involved in erythrocyte invasion, were also identified. For surveillance applications, continental-informative markers were found in putative drug resistance loci, and we show that organellar polymorphisms could classify P. vivax populations across continents and differentiate between Plasmodia spp.

CONCLUSIONS

This study has shown that genomic diversity that lies within and between P. vivax populations can be used to elucidate potential drug resistance and invasion mechanisms, as well as facilitate the molecular barcoding of the parasite for surveillance applications.

Clinical and molecular surveillance of drug resistant vivax malaria in Myanmar (2009-2016)

BACKGROUND

One of the major challenges for control and elimination of malaria is ongoing spread and emergence of drug resistance. While epidemiology and surveillance of the drug resistance in falciparum malaria is being explored globally, there are few studies on drug resistance vivax malaria.

METHODS

To assess the spread of drug-resistant vivax malaria in Myanmar, a multisite, prospective, longitudinal study with retrospective analysis of previous therapeutic efficacy studies, was conducted. A total of 906 from nine study sites were included in retrospective analysis and 208 from three study sites in prospective study. Uncomplicated vivax mono-infected patients were recruited and monitored with longitudinal follow-up until day 28 after treatment with chloroquine. Amplification and sequence analysis of molecular markers, such as mutations in pvcrt-O, pvmdr1, pvdhps and pvdhfr, were done in day-0 samples in prospective study.

RESULTS

Clinical failure cases were found only in Kawthaung, southern Myanmar and western Myanmar sites within 2009-2016. Chloroquine resistance markers, pvcrt-O 'AAG' insertion and pvmdr1 mutation (Y976F) showed higher mutant rate in southern and central Myanmar than western site: 66.7, 72.7 vs 48.3% and 26.7, 17.0 vs 1.7%, respectively. A similar pattern of significantly higher mutant rate of antifolate resistance markers, pvdhps (S382A, K512M, A553G) and pvdhfr (F57L/I, S58R, T61M, S117T/N) were noted.

CONCLUSIONS

Although clinical failure rate was low, widespread distribution of chloroquine and antifolate resistance molecular makers alert to the emergence and spread of drug resistance vivax malaria in Myanmar. Proper strategy and action plan to eliminate and contain the resistant strain strengthened together with clinical and molecular surveillance on drug resistance vivax is recommended.