Genomic analysis of Indian isolates of Plasmodium falciparum: Implications for drug resistance and virulence factors

The emergence of drug resistance to frontline treatments such as Artemisinin-based combination therapy (ACT) is a major obstacle to the control and eradication of malaria. This problem is compounded by the inherent genetic variability of the parasites, as many established markers of resistance do not accurately predict the drug-resistant status. There have been reports of declining effectiveness of ACT in the West Bengal and Northeast regions of India, which have traditionally been areas of drug resistance emergence in the country. Monitoring the genetic makeup of a population can help to identify the potential for drug resistance markers associated with it and evaluate the effectiveness of interventions aimed at reducing the spread of malaria. In this study, we performed whole genome sequencing of 53 isolates of Plasmodium falciparum from West Bengal and compared their genetic makeup to isolates from Southeast Asia (SEA) and Africa. We found that the Indian isolates had a distinct genetic makeup compared to those from SEA and Africa, and were more similar to African isolates, with a high prevalence of mutations associated with antigenic variation genes. The Indian isolates also showed a high prevalence of markers of chloroquine resistance (mutations in Pfcrt) and multidrug resistance (mutations in Pfmdr1), but no known mutations associated with artemisinin resistance in the PfKelch13 gene. Interestingly, we observed a novel L152V mutation in PfKelch13 gene and other novel mutations in genes involved in ubiquitination and vesicular transport that have been reported to support artemisinin resistance in the early stages of ACT resistance in the absence of PfKelch13 polymorphisms. Thus, our study highlights the importance of region-specific genomic surveillance for artemisinin resistance and the need for continued monitoring of resistance to artemisinin and its partner drugs.