Presence of novel triple mutations in the pvdhfr from Plasmodium vivax in Mangaluru city area in the southwestern coastal region of India

Genetic diversity of the PvMSP-3α gene in Plasmodium vivax isolates circulating in the National Capital Region (NCR) of India

Malaria is still a public health problem in tropical countries like India; major malaria parasite species are Plasmodium falciparum and P. vivax. Of which, P. vivax is responsible for ∼40% of the malaria burden at least in the Indian scenario. Unfortunately, there is limited data on the population structure and genetic diversity of P. vivax parasites in India. In this study, we investigated the genetic diversity of P. vivax strains in the South-west district, Delhi and, Nuh district, Haryana [National Capital Region (NCR)], using a polymorphic marker- P. vivax merozoite surface protein-3α (PvMSP-3α) gene. Dried blood spots from microscopically confirmed P. vivax patients were used for investigation of the PvMSP-3α gene. PCR-RFLP was performed on the PvMSP-3α gene to investigate the genotypes and allelic variability with HhaI and AluI restriction enzymes. In total, 40 successfully PCR amplified PvMSP-3α gene segments were subjected to RFLP analysis. Amplified products showed three different base pair size variations viz. genotype A in 31(77.5%), genotype B in 4(10%) and genotype C in 5(12.5%) P. vivax specimens. RFLP with HhaI and AluI revealed 17 (H1-H17) and 25 (A1-A25) allelic variants, respectively. Interestingly, two similar sub-allelic variants, ie. H8 (with HhaI), and A4 (with AluI) clustered within the rural area of Nuh district, Haryana in two samples. With this study, we propose to commission such type of genetic diversity analysis of P. vivax to investigate the circulating genotypes of the parasites from distinct geographical locations across India, that can have significant implications in understanding the population structures of P. vivax.

Drug resistance markers in Plasmodium vivax isolates from a Kanchanaburi province, Thailand between January to May 2023

BACKGROUND

Plasmodium vivax has become the predominant species in the border regions of Thailand. The emergence and spread of antimalarial drug resistance in P. vivax is one of the significant challenges for malaria control. Continuous surveillance of drug resistance is therefore necessary for monitoring the development of drug resistance in the region. This study aims to investigate the prevalence of the mutation in the P. vivax multidrug resistant 1 (Pvmdr1), dihydrofolate reductase (Pvdhfr), and dihydropteroate synthetase (Pvdhps) genes conferred resistance to chloroquine (CQ), pyrimethamine (P) and sulfadoxine (S), respectively.

METHOD

100 P. vivax isolates were obtained between January to May 2023 from a Kanchanaburi province, western Thailand. Nucleotide sequences of Pvmdr1, Pvdhfr, and Pvdhps genes were amplified and sequenced. The frequency of single nucleotide polymorphisms (SNPs)-haplotypes of drug-resistant alleles was assessed. The linkage disequilibrium (LD) tests were also analyzed.

RESULTS

In Pvmdr1, T958M, Y976F, and F1076L, mutations were detected in 100%, 21%, and 23% of the isolates, respectively. In Pvdhfr, the quadruple mutant allele (I57R58M61T117) prevailed in 84% of the samples, followed by (L57R58M61T117) in 11%. For Pvdhps, the double mutant allele (G383G553) was detected (48%), followed by the triple mutant allele (G383M512G553) (47%) of the isolates. The most prevalent combination of Pvdhfr (I57R58M61T117) and Pvdhps (G383G553) alleles was sextuple mutated haplotypes (48%). For LD analysis, the association in the SNPs pairs was found between the intragenic and intergenic regions of the Pvdhfr and Pvdhps genes.

CONCLUSION

The study has recently updated the high prevalence of three gene mutations associated with CQ and SP resistance. Genetic monitoring is therefore important to intensify in the regions to further assess the spread of drug resistant. Our data also provide evidence on the distribution of drug resistance for the early warning system, thereby threatening P. vivax malaria treatment policy decisions at the national level.

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.

Nutrient Limitation Mimics Artemisinin Tolerance in Malaria

Mounting evidence demonstrates that nutritional environment can alter pathogen drug sensitivity. While the rich media used for in vitro culture contains supraphysiological nutrient concentrations, pathogens encounter a relatively restrictive environment in vivo. We assessed the effect of nutrient limitation on the protozoan parasite that causes malaria and demonstrated that short-term growth under physiologically relevant mild nutrient stress (or "metabolic priming") triggers increased tolerance of a potent antimalarial drug. We observed beneficial effects using both short-term survival assays and longer-term proliferation studies, where metabolic priming increases parasite survival to a level previously defined as resistant (>1% survival). We performed these assessments by either decreasing single nutrients that have distinct roles in metabolism or using a media formulation that simulates the human plasma environment. We determined that priming-induced tolerance was restricted to parasites that had newly invaded the host red blood cell, but the effect was not dependent on genetic background. The molecular mechanisms of this intrinsic effect mimic aspects of genetic tolerance, including translational repression and protein export. This finding suggests that regardless of the impact on survival rates, environmental stress could stimulate changes that ultimately directly contribute to drug tolerance. Because metabolic stress is likely to occur more frequently in vivo compared to the stable in vitro environment, priming-induced drug tolerance has ramifications for how in vitro results translate to in vivo studies. Improving our understanding of how pathogens adjust their metabolism to impact survival of current and future drugs is an important avenue of research to slow the evolution of resistance. IMPORTANCE There is a dire need for effective treatments against microbial pathogens. Yet, the continuing emergence of drug resistance necessitates a deeper knowledge of how pathogens respond to treatments. We have long appreciated the contribution of genetic evolution to drug resistance, but transient metabolic changes that arise in response to environmental factors are less recognized. Here, we demonstrate that short-term growth of malaria parasites in a nutrient-limiting environment triggers cellular changes that lead to better survival of drug treatment. We found that these strategies are similar to those employed by drug-tolerant parasites, which suggests that starvation "primes" parasites to survive and potentially evolve resistance. Since the environment of the human host is relatively nutrient restrictive compared to growth conditions in standard laboratory culture, this discovery highlights the important connections among nutrient levels, protective cellular pathways, and resistance evolution.

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.

Prevalence of antifolate drug resistance markers in Plasmodium vivax in China

The dihydrofolate reductase (dhfr) and dihydropteroate synthetase (dhps) genes of Plasmodium vivax, as antifolate resistance-associated genes were used for drug resistance surveillance. A total of 375 P. vivax isolates collected from different geographical locations in China in 2009-2019 were used to sequence Pvdhfr and Pvdhps. The majority of the isolates harbored a mutant type allele for Pvdhfr (94.5%) and Pvdhps (68.2%). The most predominant point mutations were S117T/N (77.7%) in Pvdhfr and A383G (66.8%) in Pvdhps. Amino acid changes were identified at nine residues in Pvdhfr. A quadruple-mutant haplotype at 57, 58, 61, and 117 was the most frequent (57.4%) among 16 distinct Pvdhfr haplotypes. Mutations in Pvdhps were detected at six codons, and the double-mutant A383G/A553G was the most prevalent (39.3%). Pvdhfr exhibited a higher mutation prevalence and greater diversity than Pvdhps in China. Most isolates from Yunnan carried multiple mutant haplotypes, while the majority of samples from temperate regions and Hainan Island harbored the wild type or single mutant type. This study indicated that the antifolate resistance levels of P. vivax parasites were different across China and molecular markers could be used to rapidly monitor drug resistance. Results provided evidence for updating national drug policy and treatment guidelines.

Surveillance of molecular markers of antimalarial drug resistance in Plasmodium falciparum and Plasmodium vivax in Federally Administered Tribal Area (FATA), Pakistan

This molecular epidemiological study was designed to determine the antimalarial drug resistance pattern, and the genetic diversity of malaria isolates collected from a war-altered Federally Administered Tribal Area (FATA), in Pakistan. Clinical isolates were collected from Bajaur, Mohmand, Khyber, Orakzai and Kurram agencies of FATA region between May 2017 and May 2018, and they underwent DNA extraction and amplification. The investigation of gene polymorphisms in drug resistance genes (dhfr, dhps, crt, and mdr1) of Plasmodium falciparum and Plasmodium vivax was carried out by pyrosequencing and Sanger sequencing, respectively. Out of 679 PCR-confirmed malaria samples, 523 (77%) were P. vivax, 121 (18%) P. falciparum, and 35 (5%) had mixed-species infections. All P. falciparum isolates had pfdhfr double mutants (C59R+S108N), while pfdhfr/pfdhps triple mutants (C59R+S108N+A437G) were detected in 11.5% of the samples. About 97.4% of P. falciparum isolates contained pfcrt K76T mutation, while pfmdr1 N86Y and Y184F mutations were present in 18.2% and 10.2% of the samples. P. vivax pvdhfr S58R mutation was present in 24.9% of isolates and the S117N mutation in 36.2%, while no mutation in the pvdhps gene was found. Pvmdr1 F1076L mutation was found in nearly all samples, as it was observed in 98.9% of isolates. No significant anti-folate and chloroquine resistance was observed in P. vivax; however, mutations associated with antifolate-resistance were found, and the chloroquine-resistant gene has been observed in 100% of P. falciparum isolates. Chloroquine and sulphadoxine-pyrimethamine resistance were found to be high in P. falciparum and low in P. vivax. Chloroquine could still be used for P. vivax infection but need to be tested in vivo, whereas a replacement of the artemisinin combination therapy for P. falciparum appears to be justified.

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.

Assessing the in vitro sensitivity with associated drug resistance polymorphisms in Plasmodium vivax clinical isolates from Delhi, India

The drug resistance of Plasmodium vivax in clinical cases remains largely unknown till date because of the difficulty in diagnosing the resistant P. vivax strains. The present study was undertaken to determine the prevalence of mutant alleles in drug resistance genes viz P. vivax multi-drug resistance (pvmdr-1), chloroquine resistance transporter (pvcrt-o), dihydrofolate reductase (pvdhfr) and dihydropteroate synthase (pvdhps) along with in vitro chloroquine (CQ) sensitivity in P. vivax clinical isolates. During August-October 2017 a total of 86 samples of the febrile patients were screened and 31 samples were found to be positive for P. vivax in Safdarjung hospital, New Delhi. Sequence genotyping of the drug resistance genes was carried out in these P. vivax samples and in vitro CQ susceptibility for 23 isolates was determined by the schizont maturation assay (SMA). The CQ inhibitory concentrations (IC₅₀) for the clinical isolates was found to be in the range of 25.6-176.7 nM. All the 31 clinical isolates analyzed for pvmdr-1 gene, showed mutant alleles and in only two isolates novel mutations at 861 and 898 codons were observed. Sequence analysis of pvcrt-o, pvdhfr and pvdhps genes revealed wild type genotypes in all the 31 studied isolates. The presence of mutations in pvmdr-1 gene and the increase in the CQ IC₅₀ value indicates the possibility of shift in drug tolerance where CQ with primaquine (PQ) is still the first line of treatment for P. vivax malaria in the country. The regular molecular surveillance in P. vivax would provide useful information for the policy makers of the malaria control programme.

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.