Prevalence of polymorphisms in antifolate drug resistance molecular marker genes pvdhfr and pvdhps in clinical isolates of Plasmodium vivax from Kolkata, India

Characterization of drug resistance genes in Indian Plasmodium falciparum and Plasmodium vivax field isolates

One of the major challenges for malaria control and elimination is the spread and emergence of antimalarial drug resistance. Mutations in Plasmodium falciparum (Pf) and Plasmodium vivax (Pv) field isolates for five drug resistance genes viz. crt, mdr1, dhps, dhfr and kelch known to confer resistance to choloroquine (CQ), sulfadoxine-pyrimethamine (SP) and artemisinin (ART) and its derivatives were analyzed. A total of 342 symptomatic isolates of P. falciparum (Pf) and P. vivax (Pv) from 1993 to 2014 were retrieved from malaria parasite repository at National Institute of Malaria Research (NIMR). Sample DNA was extracted from dried blood spots and various targeted single nucleotide polymorphisms (SNPs) associated with antimalarial drug resistance were analysed for these isolates. 72S (67.7%) and 76T (83.8%) mutations along with SVMNT haplotype (67.7%) predominated the study population for Pfcrt. The most prevalent SNPs were 108N (73.2%) and 437G (24.8%) and the most prevalent haplotypes were ACNRNI (51.9%) and SAKAA (74.5%) in Pfdhfr and Pfdhps respectively. Only two mutations in Pfmdr1, 86Y (26.31%) and 184F (56.26%), were seen frequently in our study population. No mutations associated with Pfk13 were observed. For Pv, all the studied isolates showed two Pvdhps mutations, 383G and 553G, and two Pfdhfr mutations, 58R and 117N. Similarly, three mutations, viz. 958M, 908L and 1076L were found in Pvmdr1. No variations were observed in Pvcrt-o and Pvk12 genes. Overall, our study demonstrates an increase in mutations associated with SP resistance in both Pf and Pv, however, no single nucleotide polymorphisms (SNPs) associated with ART resistance have been observed for either species. Various SNPs associated with CQ resistance were seen in Pf; whereas only Pvmdr1 associated resistant SNPs were observed in Pv. Therefore, molecular characterization of drug resistance genes is essential for timely monitoring and prevention of malaria by identifying the circulating drug resistant parasites in the country.

Nationwide spatiotemporal drug resistance genetic profiling from over three decades in Indian Plasmodium falciparum and Plasmodium vivax isolates

BACKGROUND

Drug resistance is a serious impediment to efficient control and elimination of malaria in endemic areas.

METHODS

This study aimed at analysing the genetic profile of molecular drug resistance in Plasmodium falciparum and Plasmodium vivax parasites from India over a ~ 30-year period (1993-2019). Blood samples of P. falciparum and/or P. vivax-infected patients were collected from 14 regions across India. Plasmodial genome was extracted and used for PCR amplification and sequencing of drug resistance genes in P. falciparum (crt, dhps, dhfr, mdr1, k13) and P. vivax (crt-o, dhps, dhfr, mdr1, k12) field isolates.

RESULTS

The double mutant pfcrt SVMNT was highly predominant across the country over three decades, with restricted presence of triple mutant CVIET from Maharashtra in 2012. High rates of pfdhfr-pfdhps quadruple mutants were observed with marginal presence of "fully resistant" quintuple mutant ACIRNI-ISGEAA. Also, resistant pfdhfr and pfdhps haplotype has significantly increased in Delhi between 1994 and 2010. For pfmdr1, only 86Y and 184F mutations were present while no pfk13 mutations associated with artemisinin resistance were observed. Regarding P. vivax isolates, the pvcrt-o K10 "AAG" insertion was absent in all samples collected from Delhi in 2017. Pvdhps double mutant SGNAV was found only in Goa samples of year 2008 for the first time. The pvmdr1 908L, 958M and 1076L mutations were highly prevalent in Delhi and Haryana between 2015 and 2019 at complete fixation. One nonsynonymous novel pvk12 polymorphism was identified (K264R) in Goa.

CONCLUSIONS

These findings support continuous surveillance and characterization of P. falciparum and P. vivax populations as proxy for effectiveness of anti-malarial drugs in India, especially for independent emergence of artemisinin drug resistance as recently seen in Africa.

Novel highly-multiplexed AmpliSeq targeted assay for Plasmodium vivax genetic surveillance use cases at multiple geographical scales

Although the power of genetic surveillance tools has been acknowledged widely, there is an urgent need in malaria endemic countries for feasible and cost-effective tools to implement in national malaria control programs (NMCPs) that can generate evidence to guide malaria control and elimination strategies, especially in the case of Plasmodium vivax. Several genetic surveillance applications (‘use cases’) have been identified to align research, technology development, and public health efforts, requiring different types of molecular markers. Here we present a new highly-multiplexed deep sequencing assay (Pv AmpliSeq). The assay targets the 33-SNP vivaxGEN-geo panel for country-level classification, and a newly designed 42-SNP within-country barcode for analysis of parasite dynamics in Vietnam and 11 putative drug resistance genes in a highly multiplexed NGS protocol with easy workflow, applicable for many different genetic surveillance use cases. The Pv AmpliSeq assay was validated using: 1) isolates from travelers and migrants in Belgium, and 2) routine collections of the national malaria control program at sentinel sites in Vietnam. The assay targets 229 amplicons and achieved a high depth of coverage (mean 595.7 ± 481) and high accuracy (mean error-rate of 0.013 ± 0.007). P. vivax parasites could be characterized from dried blood spots with a minimum of 5 parasites/µL and 10% of minority-clones. The assay achieved good spatial specificity for between-country prediction of origin using the 33-SNP vivaxGEN-geo panel that targets rare alleles specific for certain countries and regions. A high resolution for within-country diversity in Vietnam was achieved using the designed 42-SNP within-country barcode that targets common alleles (median MAF 0.34, range 0.01-0.49. Many variants were detected in (putative) drug resistance genes, with different predominant haplotypes in the pvmdr1 and pvcrt genes in different provinces in Vietnam. The capacity of the assay for high resolution identity-by-descent (IBD) analysis was demonstrated and identified a high rate of shared ancestry within Gia Lai Province in the Central Highlands of Vietnam, as well as between the coastal province of Binh Thuan and Lam Dong. Our approach performed well in geographically differentiating isolates at multiple spatial scales, detecting variants in putative resistance genes, and can be easily adjusted to suit the needs in other settings in a country or region. We prioritize making this tool available to researchers and NMCPs in endemic countries to increase ownership and ensure data usage for decision-making and malaria policy.

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 Malaria Viewed through the Lens of an Eradicated European Strain

The protozoan Plasmodium vivax is responsible for 42% of all cases of malaria outside Africa. The parasite is currently largely restricted to tropical and subtropical latitudes in Asia, Oceania, and the Americas. Though, it was historically present in most of Europe before being finally eradicated during the second half of the 20th century. The lack of genomic information on the extinct European lineage has prevented a clear understanding of historical population structuring and past migrations of P. vivax. We used medical microscope slides prepared in 1944 from malaria-affected patients from the Ebro Delta in Spain, one of the last footholds of malaria in Europe, to generate a genome of a European P. vivax strain. Population genetics and phylogenetic analyses placed this strain basal to a cluster including samples from the Americas. This genome allowed us to calibrate a genomic mutation rate for P. vivax, and to estimate the mean age of the last common ancestor between European and American strains to the 15th century. This date points to an introduction of the parasite during the European colonization of the Americas. In addition, we found that some known variants for resistance to antimalarial drugs, including Chloroquine and Sulfadoxine, were already present in this European strain, predating their use. Our results shed light on the evolution of an important human pathogen and illustrate the value of antique medical collections as a resource for retrieving genomic information on pathogens from the past.

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.

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

Background

Non-synonymous mutations in dhfr and dhps genes in Plasmodium vivax are associated with sulfadoxine–pyrimethamine (SP) resistance. The present study aimed to assess the prevalence of point mutations in P. vivax dhfr (pvdhfr) and P. vivax dhps (pvdhps) genes in three countries: Lao PDR, India and Colombia.

Methods

Samples from 203 microscopically diagnosed vivax malaria were collected from the three countries. Five codons at positions 13, 57, 58, 61, and 117 of pvdhfr and two codons at positions 383 and 553 of pvdhps were examined by polymerase chain reaction-restriction fragment length polymorphism methodology.

Results

The largest number of 58R/117 N double mutations in pvdhfr was observed in Colombia (94.3 %), while the corresponding wild-type amino acids were found at high frequencies in Lao PDR during 2001–2004 (57.8 %). Size polymorphism analysis of the tandem repeats within pvdhfr revealed that 74.3 % of all the isolates carried the type B variant. Eighty-nine per cent of all the isolates examined carried wild-type pvdhps A383 and A553.

Conclusions

Although SP is not generally used to treat P. vivax infections, mutations in dhfr and dhps that confer antifolate resistance in P. vivax are common. The data strongly suggest that, when used primarily to treat falciparum malaria, SP can exert a substantial selective pressure on P. vivax populations, and this can lead to point mutations in dhfr and dhps. Accurate data on the global geographic distribution of dhfr and dhps genotypes should help to inform anti-malarial drug-use policies.

Low prevalence of dihydro folate reductase (dhfr) and dihydropteroate synthase (dhps) quadruple and quintuple mutant alleles associated with SP resistance in Plasmodium vivax isolates of West Bengal, India

Background

Emergence of chloroquine resistant Plasmodium vivax is a serious obstacle towards malaria control in India. This study elucidates the temporal pattern of antifolate [sulfadoxine–pyrimethamine (SP)] resistance in P. vivax infection by means of genetic polymorphisms, especially analysing the single nucleotide polymorphisms of dihydrofolate reductase (pvdhfr) and dihydropteroate synthase (pvdhps) gene among the field isolates of urban Kolkata Municipal Corporation and rural Purulia region of West Bengal, India.

Methods

Blood samples were collected from 99 microscopically diagnosed P. vivax patients (52 from Kolkata Municipal Corporation and 47 from Purulia). Parasitic DNA was extracted followed by polymerase chain reaction and sequencing of different codons of pvdhfr gene (15, 33, 50, 57, 58, 61, 64, 117, and 173 codons) and pvdhps gene (373, 380, 382, 383, 384, 512, 553, 585, and 601 codons) were performed to identify the mutations.

Results

Prevalence of double mutant dhfr A₁₅P₃₃N₅₀F₅₇R₅₈T₆₁V₆₄N₁₁₇I₁₇₃ allele (53.85 %) was observed in Kolkata Municipal Corporation (KMC) whereas in Purulia, wild dhfr A₁₅P₃₃N₅₀F₅₇S₅₈T₆₁V₆₄S₁₁₇I₁₇₃ allele was predominated (48.94 %). In pvdhps gene a significant number of isolates (17.31 %) in KMC contained the double mutant S₃₇₃E₃₈₀S₃₈₂G₃₈₃P₃₈₄K₅₁₂G₅₅₃V₅₈₅M₆₀₁ allele. pvdhfr and pvdhps combination haplotype revealed the emergence of quadruple (13.46 %) and quintuple (3.84 %) mutant allele in KMC, which might result in poor clinical response against antifolate drugs.

Conclusion

The study reveals that P. vivax parasites in rural Purulia may still be susceptible to SP but additional caution should be taken for treatment of vivax malaria in KMC to limit the blooming of quadruple and quintuple mutant allele in the remainder of the West Bengal, India.

Advances in genetics and genomics: use and limitations in achieving malaria elimination goals

Success of the global research agenda towards eradication of malaria will depend on the development of new tools, including drugs, vaccines, insecticides and diagnostics. Genetic and genomic information now available for the malaria parasites, their mosquito vectors and human host, can be harnessed to both develop these tools and monitor their effectiveness. Here we review and provide specific examples of current technological advances and how these genetic and genomic tools have increased our knowledge of host, parasite and vector biology in relation to malaria elimination and in turn enhanced the potential to reach that goal. We then discuss limitations of these tools and future prospects for the successful achievement of global malaria elimination goals.

Next-Generation Sequencing of Plasmodium vivax Patient Samples Shows Evidence of Direct Evolution in Drug-Resistance Genes

Understanding the mechanisms of drug resistance in Plasmodium vivax, the parasite that causes the most widespread form of human malaria, is complicated by the lack of a suitable long-term cell culture system for this parasite. In contrast to P. falciparum, which can be more readily manipulated in the laboratory, insights about parasite biology need to be inferred from human studies. Here we analyze the genomes of parasites within 10 human P. vivax infections from the Peruvian Amazon. Using next-generation sequencing we show that some P. vivax infections analyzed from the region are likely polyclonal. Despite their polyclonality we observe limited parasite genetic diversity by showing that three or fewer haplotypes comprise 94% of the examined genomes, suggesting the recent introduction of parasites into this geographic region. In contrast we find more than three haplotypes in putative drug-resistance genes, including the gene encoding dihydrofolate reductase-thymidylate synthase and the P. vivax multidrug resistance associated transporter, suggesting that resistance mutations have arisen independently. Additionally, several drug-resistance genes are located in genomic regions with evidence of increased copy number. Our data suggest that whole genome sequencing of malaria parasites from patients may provide more insight about the evolution of drug resistance than genetic linkage or association studies, especially in geographical regions with limited parasite genetic diversity.

Antimalarial Drug Resistance: Literature Review and Activities and Findings of the ICEMR Network

Antimalarial drugs are key tools for the control and elimination of malaria. Recent decreases in the global malaria burden are likely due, in part, to the deployment of artemisinin-based combination therapies. Therefore, the emergence and potential spread of artemisinin-resistant parasites in southeast Asia and changes in sensitivities to artemisinin partner drugs have raised concerns. In recognition of this urgent threat, the International Centers of Excellence for Malaria Research (ICEMRs) are closely monitoring antimalarial drug efficacy and studying the mechanisms underlying drug resistance. At multiple sentinel sites of the global ICEMR network, research activities include clinical studies to track the efficacies of antimalarial drugs, ex vivo/in vitro assays to measure drug susceptibilities of parasite isolates, and characterization of resistance-mediating parasite polymorphisms. Taken together, these efforts offer an increasingly comprehensive assessment of the efficacies of antimalarial therapies, and enable us to predict the emergence of drug resistance and to guide local antimalarial drug policies. Here we briefly review worldwide antimalarial drug resistance concerns, summarize research activities of the ICEMRs related to drug resistance, and assess the global impacts of the ICEMR programs.