Progressive increase in point mutations associates chloroquine resistance: Even after withdrawal of chloroquine use in 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.

Unraveling Key Chloroquine Resistance-Associated Alleles Among Plasmodium falciparum Isolates in South Darfur State, Sudan Twelve Years After Drug Withdrawal

Background

Due to the increasing resistance of Plasmodium falciparum to chloroquine (CQ) in Sudan, a shift from CQ to artesunate combined with sulfadoxine/pyrimethamine as a first-line treatment for uncomplicated falciparum malaria was adopted in 2004. This study aimed to determine the frequency distribution of K76T and N86Y mutations in P. falciparum chloroquine resistance transporter (pfcrt) and P. falciparum multidrug resistance 1 (pfmdr1) genes as key markers of resistance to CQ among P. falciparum isolates from patients in Nyala district of South Darfur state, west of Sudan.

Methods

A descriptive, cross-sectional study was conducted among 75 P. falciparum isolates from Sudanese patients diagnosed with falciparum malaria mono-infection. Parasite DNA was extracted from dried blood spots and amplified using a nested polymerase chain reaction (PCR). Then, restriction fragment length polymorphism (RFLP) was used to detect the genetic polymorphisms in codons 76 of pfcrt and 86 of pfmdr1. PCR-RFLP products were analyzed using 1.5% gel electrophoresis to identify the genetic polymorphisms in the studied codons. The wild-type (pfcrt K76 and pfmdr1 N86), mutant (pfcrt 76T and pfmdr1 86Y) and mixed-type (pfcrt K76T and pfmdr1 N86Y) alleles were expressed as frequencies and proportions.

Results

The wild-type pfcrt K76 allele was observed among 34.7% of isolates and the mutant 76T allele among 20% of isolates, while the mixed-type K76T allele was observed among 45.3% of isolates. On the other hand, 54.7% of isolates harbored the wild-type pfmdr1 N86 allele and 5.3% of isolates had the mutant 86Y allele, while the mixed-type N86Y allele was observed among 40% of isolates.

Conclusion

The key molecular markers associated with CQ resistance (pfcrt 76T and pfmdr1 86Y) are still circulating in high frequency among P. falciparum isolates in South Darfur state, about twelve years after the official withdrawal of the drug as a treatment for uncomplicated falciparum malaria.

Molecular surveillance of Kelch-13 gene in Plasmodium falciparum field isolates from Mayurbhanj District, Odisha, India, and in silico artemisinin-Kelch-13 protein interaction study

The global malaria control and elimination program faces major threats due to the emergence and transmission of the anti-malarial drug-resistant strain of Plasmodium falciparum. Monitoring of artemisinin (ART) resistance marker Kelch-13 in the malaria-endemic region is essential in mitigating the disease's morbidity and mortality. The current study aimed to generate baseline information for further surveillance in the future. The current research was designed and conducted from July 2019 to June 2021 to monitor Pfkelch13 mutation at the molecular level in the eastern region of India. We also conducted an in silico study to understand the drug-protein interactions between ART and the protein crystal of PfKelch13 (KELCH) with PDB id:4ZGC. The kelch-13 gene was amplified by nested polymerase chain reaction (PCR) and sequenced through the Sanger sequencing method. Reference 3D7 clone (PF3D7_1343700) was used to align and probe all the sequences. The sequence analysis showed the absence of validated or associated mutation in the Kelch-13 propeller domain. The absence of natural selection in drug resistance was confirmed by the Tajima test. Further, in silico interaction studies between the drug ART and the Kelch propeller domain of P. falciparum were evaluated by structure predictions, molecular docking, molecular dynamics (MD) simulations, and estimations of binding free energies for the KELCH-ART complex. The results were compared with the apoprotein (KELCH-APO). The study confirmed the favorable binding of ART with the Kelch-13 propeller domain.

Artemisinin combination therapy fails even in the absence of Plasmodium falciparum kelch13 gene polymorphism in Central India

Artemisinin is the frontline fast-acting anti-malarial against P. falciparum. Emergence and spread of resistant parasite in eastern-India poses a threat to national malaria control programs. Therefore, the objective of our study is to evaluate the artesunate-sulfadoxine-pyrimethamine efficacy in Central India. 180 monoclonal P. falciparum-infected patients received standard ASSP therapy during August 2015–January 2017, soon after diagnosis and monitored over next 42-days. Artemisinin-resistance was assessed through in-vivo parasite clearance half-life (PC_1/2), ex-vivo ring-stage survivability (RSA), and genome analysis of kelch13 and other candidate gene (pfcrt, pfmdr1, pfatpase 6, pfdhfr and pfdhps). Of 180 P. falciparum positive patients, 9.5% showed increased PC_1/2 (> 5.5 h), among them eleven isolates (6.1%) showed reduced sensitivity to RSA. In 4.4% of cases, parasites were not cleared by 72 h and showed prolonged PC_1/2(5.6 h) (P < 0.005) along with significantly higher RSA (2.2%) than cured patients (0.4%). None of day-3 positive isolates contained the pfkelch13 mutation implicated in artemisinin resistance. Parasite recrudescence was observed in 5.6% patients, which was associated with triple dhfr–dhps (A₁₆I₅₁R₅₉N₁₀₈I₁₆₄–S₄₃₆G₄₃₇K₅₄₀G₅₈₁T₆₁₃) combination mutation. Emergence of reduced sensitivity to artesunate-sulfadoxine-pyrimethamine, in central India highlighted the risk toward spread of resistant parasite across different parts of India. Day-3 positive parasite, featuring the phenotype of artemisinin-resistance without pfkelch13 mutation, suggested kelch13-independent artemisinin-resistance.

Plasmodium falciparum multidrug resistance gene-1 polymorphisms in Northern Nigeria: implications for the continued use of artemether-lumefantrine in the region

Background

The analysis of single nucleotide polymorphism (SNPs) in drug-resistance associated genes is a commonly used strategy for the surveillance of anti-malarial drug resistance in populations of parasites. The present study was designed and performed to provide genetic epidemiological data of the prevalence of N86Y-Y184F-D1246Y SNPs in Plasmodium falciparum multidrug resistance 1 (pfmdr1) in the malaria hotspot of Northern Nigeria.

Methods

Plasmodium falciparum-positive blood samples on Whatman-3MM filter papers were collected from 750 symptomatic patients from four states (Kano, Kaduna, Yobe and Adamawa) in Northern Nigeria, and genotyped via BigDye (v3.1) terminator cycle sequencing for the presence of three SNPs in pfmdr1. SNPs in pfmdr1 were used to construct NYD, NYY, NFY, NFD, YYY, YYD, YFD and YFY haplotypes, and all data were analysed using Pearson Chi square and Fisher’s exact (FE) tests.

Results

The prevalence of the pfmdr1 86Y allele was highest in Kaduna (12.50%, ² = 10.50, P = 0.02), whilst the 184F allele was highest in Kano (73.10%, ² = 13.20, P = 0.00), and the pfmdr1 1246Y allele was highest in Yobe (5.26%, ² = 9.20, P = 0.03). The NFD haplotype had the highest prevalence of 69.81% in Kano (² = 36.10, P = 0.00), followed by NYD with a prevalence of 49.00% in Adamawa, then YFD with prevalence of 11.46% in Kaduna. The YYY haplotype was not observed in any of the studied states.

Conclusion

The present study suggests that strains of P. falciparum with reduced sensitivity to the lumefantrine component of AL exist in Northern Nigeria and predominate in the North-West region.

Discovery of FNDR-20123, a histone deacetylase inhibitor for the treatment of Plasmodium falciparum malaria

Background

Emergence of anti-malarial drug resistance and perpetual increase in malaria incidence necessitates the development of novel anti-malarials. Histone deacetylases (HDAC) has been shown to be a promising target for malaria, despite this, there are no HDAC inhibitors in clinical trials for malaria treatment. This can be attributed to the poor pharmacokinetics, bioavailability and selectivity of the HDAC inhibitors.

Methods

A collection of HDAC inhibitors were screened for anti-malarial activity, and the best candidate was profiled in parasite-killing kinetics, growth inhibition of sensitive and multi-drug resistant (MDR) strains and against gametocytes. Absorption, distribution, metabolism and excretion pharmacokinetics (ADME-PK) parameters of FNDR-20123 were determined, and in vivo efficacy was studied in a mouse model for Plasmodium falciparum infection.

Results

A compound library of HDAC inhibitors (180 in number) was screened for anti-malarial activity, of which FNDR-20123 was the most potent candidate. The compound had been shown to inhibit Plasmodium HDAC with IC₅₀ of 31 nM and human HDAC with IC₅₀ of 3 nM. The IC₅₀ obtained for P. falciparum in asexual blood-stage assay was 42 nM. When compared to atovaquone and pyrimethamine, the killing profiles of FNDR-20123 were better than atovaquone and comparable to pyrimethamine. The IC₅₀ values for the growth inhibition of sensitive and MDR strains were similar, indicating that there is no cross-resistance and a low risk of resistance development. The selected compound was also active against gametocytes, indicating a potential for transmission control: IC₅₀ values being 190 nM for male and > 5 µM for female gametocytes. FNDR-20123 is a stable candidate in human/mouse/rat liver microsomes (> 75% remaining post 2-h incubation), exhibits low plasma protein binding (57% in humans) with no human Ether-à-go–go-Related Gene (hERG) liability (> 100 µM), and does not inhibit any of the cytochrome P450 (CYP) isoforms tested (IC₅₀ > 25 µM). It also shows negligible cytotoxicity to HepG-2 and THP-1 cell lines. The oral pharmacokinetics in rats at 100 mg/kg body weight shows good exposures (C_max = 1.1 µM) and half-life (T_1/2 = 5.5 h). Furthermore, a 14-day toxicokinetic study at 100 mg/kg daily dose did not show any abnormality in body weight or gross organ pathology. FNDR-20123 is also able to reduce parasitaemia significantly in a mouse model for P. falciparum infection when dosed orally and subcutaneously.

Conclusion

FNDR-20123 may be a suitable candidate for the treatment of malaria, which can be further developed.

Epidemiology of malaria and chloroquine resistance in Mizoram, northeastern India, a malaria-endemic region bordering Myanmar

Background

Mizoram, a northeastern state in India, shares international borders with Myanmar and Bangladesh and is considered to be one of the key routes through which drug-resistant parasites of Southeast Asia enter mainland India. Despite its strategic location and importance, malaria epidemiology and molecular status of chloroquine resistance had not been well documented, and since chloroquine (CQ), as the first-line treatment in Plasmodium falciparum infection was discontinued since 2008, it was expected that CQ-sensitive haplotype would be more abundant.

Methods

Malaria epidemiology data for the period 2010 to 2018 was collected from the office of State Vector Disease Control Programme. Plasmodium falciparum-positive blood samples were collected from government district hospitals, community health centres, primary health centres, sub-centres, and diagnostic centres from six malaria-prone districts. The samples were processed and analysed using genes–P. falciparum chloroquine-resistant transporter (pfcrt) and P. falciparum multidrug resistance 1 (pfmdr1) via sequencing of PCR amplicon from 2015 to 2017.

Results

Malaria occurred throughout the year and P. falciparum accounted for > 89% of total malaria cases. During 2010–2018, the highest number of malaria incidence was recorded in Lawngtlai (36% of total malaria cases; average API_2010–2018 of 34.8) while Champhai remained consistently low (0.4%; average API_2010–2018 of 0.04). Males of ≥ 15 years old contributed maximum (35.7%) among gender and age malarial distribution recorded during 2014–2018. Death due to malaria gradually decreased over the years. A higher abundance of mutated pfcrt (58.5% of the total sample analysed) and a lower prevalence of mutated pfmdr1 (48.7%) were observed. All mutations identified for pfcrt belong to the Southeast Asian CVIET haplotype. Only a single point mutation was observed at 86 (N → Y) position in pfmdr1 (48.7%). The key N86Y mutation in pfmdr1 that had been shown to modulate CQR was found in 67.1% of the samples positive for the CVIET haplotype.

Conclusions

This is the first report that details malaria epidemiology and also the molecular status of CQ-resistance in P. falciparum population of the region. The efforts of the State Vector Borne Disease Control Programme have proved to be quite effective in controlling the malaria burden in the state. Despite the discontinuation of CQ for a decade, local P. falciparum is observed with decreased CQ-sensitive haplotype. It is believed that the present findings will form a basis for further studies on genetic diversity in P. falciparum, which could confer better understanding of the complexity of the disease in Southeast Asia.

Prevalence of Plasmodium falciparum field isolates with deletions in histidine-rich protein 2 and 3 genes in context with sub-Saharan Africa and India: a systematic review and meta-analysis

BACKGROUND

In 2017, nearly 80% of malaria morbidity and mortality occurred in sub-Saharan African (SSA) countries and India. Rapid diagnostic tests (RDTs), especially those targeting histidine-rich protein 2 (PfHRP2) of Plasmodium falciparum, have become an important diagnostic tool in these malaria-endemic areas. However, the chances of RDT-oriented successful treatment are increasingly jeopardized by the appearance of mutants with deletions in pfhrp2 and pfhrp3 genes. This systematic review and meta-analysis determines the prevalence of field P. falciparum isolates with deletion in pfhrp2 and/or pfhrp3 genes and their proportion among false-negative results in the PfHRP2-based RDTs in SSA and India.

METHODS

Eight electronic databases were used for searching potentially relevant publications for the systematic review analysis, wherein the main methodological aspects of included studies were analysed and some missing links in the included studies were identified.

RESULTS

A total of 19 studies were included, 16 from SSA and 3 from India. The pooled prevalence of pfhrp2 deletions was 8 and 5% while 16 and 4% for pfhrp3 gene deletions in Africa and India, respectively. The pooled proportion of pfhrp2 gene deletions found among false negative PfHRP2-based RDTs results was about 27.0 and 69.0% in Africa and India, respectively.

CONCLUSIONS

This review study indicates a relatively high proportion of both pfhrp2/3 genes deletions in P. falciparum isolates and among false-negative malaria cases using PfHRP2-based RDT results in SSA and India. Recently the deletions in pfhrp2/3 genes have also been reported from two African countries (Nigeria and Sudan). This review emphasizes the importance of more extensive studies and standardization of studies addressing the pfhrp2/3 gene deletions in malarious areas.

Persistence of chloroquine resistance alleles in malaria endemic countries: a systematic review of burden and risk factors

Background

Chloroquine, a previous highly efficacious, easy to use and affordable anti-malarial agent was withdrawn from malaria endemic regions due to high levels of resistance. This review collated evidence from published-reviewed articles to establish prevalence of Pfcrt 76T and Pfmdr-1 86Y alleles in malaria affected countries following official discontinuation of chloroquine use.

Methods

A review protocol was developed, registered in PROSPERO (#CRD42018083957) and published in a peer-reviewed journal. Article search was done in PubMed, Scopus, Lilacs/Vhl and Embase databases by two experienced librarians (AK, RS) for the period 1990-to-Febuary 2018. Mesh terms and Boolean operators (AND, OR) were used. Data extraction form was designed in Excel spread sheet 2007. Data extraction was done by three reviewers (NL, BB and MO), discrepancies were resolved by discussion. Random effects analysis was done in Open Meta Analyst software. Heterogeneity was established using I²-statistic.

Results

A total of 4721 citations were retrieved from article search (Pubmed = 361, Lilac/vhl = 28, Science Direct = 944, Scopus = 3388). Additional targeted search resulted in three (03) eligible articles. After removal of duplicates (n = 523) and screening, 38 articles were included in the final review. Average genotyping success rate was 63.6% (18,343/28,820) for Pfcrt K76T and 93.5% (16,232/17,365) for Pfmdr-1 86Y mutations. Prevalence of Pfcrt 76T was as follows; East Africa 48.9% (2528/5242), Southern Africa 18.6% (373/2163), West Africa 58.3% (3321/6608), Asia 80.2% (1951/2436). Prevalence of Pfmdr-1 86Y was; East Africa 32.4% (1447/5722), Southern Africa 36.1% (544/1640), West Africa 52.2% (1986/4200), Asia 46.4% (1276/2217). Over half, 52.6% (20/38) of included studies reported continued unofficial chloroquine use following policy change. Studies done in Madagascar and Kenya reported re-emergence of chloroquine sensitive parasites (IC₅₀ < 30.9 nM). The average time (years) since discontinuation of chloroquine use to data collection was 8.7 ± 7.4. There was high heterogeneity (I² > 95%).

Conclusion

The prevalence of chloroquine resistance alleles among Plasmodium falciparum parasites have steadily declined since discontinuation of chloroquine use. However, Pfcrt K76T and Pfmdr-1 N86Y mutations still persist at moderate frequencies in most malaria affected countries.

Lack of significant recovery of chloroquine sensitivity in Plasmodium falciparum parasites following discontinuance of chloroquine use in Papua New Guinea

Background

Chloroquine treatment for Plasmodium falciparum has been discontinued in almost all endemic regions due to the spread of resistant isolates. Reversal of chloroquine susceptibility after chloroquine discontinuation has been reported in dozens of endemic regions. However, this phenomenon has been mostly observed in Africa and is not well documented in other malaria endemic regions. To investigate this, an ex vivo study on susceptibility to chloroquine and lumefantrine was conducted during 2016–2018 in Wewak, Papua New Guinea where chloroquine had been removed from the official malaria treatment regimen in 2010. Genotyping of pfcrt and pfmdr1 was also performed.

Results

In total, 368 patients were enrolled in this study. Average IC₅₀ values for chloroquine were 106.6, 80.5, and 87.6 nM in 2016, 2017, and 2018, respectively. These values were not significantly changed from those obtained in 2002/2003 (108 nM). The majority of parasites harboured a pfcrt K76T the mutation responsible for chloroquine resistance. However, a significant upward trend was observed in the frequency of the K76 (wild) allele from 2.3% in 2016 to 11.7% in 2018 (P = 0.008; Cochran–Armitage trend test).

Conclusions

Eight years of chloroquine withdrawal has not induced a significant recovery of susceptibility in Papua New Guinea. However, an increasing tendency of parasites harbouring chloroquine-susceptible K76 suggests a possibility of resurgence of chloroquine susceptibility in the future.

Electronic supplementary material

The online version of this article (10.1186/s12936-018-2585-x) contains supplementary material, which is available to authorized users.