Genetic variability and population structure of Plasmodium falciparum parasite populations from different malaria ecological regions of Kenya

Molecular tools are crucial for malaria elimination

The eradication of Plasmodium parasites, responsible for malaria, is a daunting global public health task. It requires a comprehensive approach that addresses symptomatic, asymptomatic, and submicroscopic cases. Overcoming this challenge relies on harnessing the power of molecular diagnostic tools, as traditional methods like microscopy and rapid diagnostic tests fall short in detecting low parasitaemia, contributing to the persistence of malaria transmission. By precisely identifying patients of all types and effectively characterizing malaria parasites, molecular tools may emerge as indispensable allies in the pursuit of malaria elimination. Furthermore, molecular tools can also provide valuable insights into parasite diversity, drug resistance patterns, and transmission dynamics, aiding in the implementation of targeted interventions and surveillance strategies. In this review, we explore the significance of molecular tools in the pursuit of malaria elimination, shedding light on their key contributions and potential impact on public health.

Putative molecular markers of Plasmodium falciparum resistance to antimalarial drugs in malaria parasites from Ghana

Introduction

Antimalarial drugs including artemisinin-based combination therapy (ACT) regimens and sulphadoxine-pyrimethamine (SP) are used in Ghana for malaria therapeutics and prophylaxis respectively. The genetic basis of Plasmodium falciparum development of drug resistance involves single nucleotide polymorphisms in genes encoding proteins for multiple cellular and metabolic processes. The prevalence of single nucleotide polymorphisms in nine P. falciparum genes linked to ACT and SP resistance in the malaria parasite population was determined.

Methods

Archived filter paper blood blot samples from patients aged 9 years and below with uncomplicated malaria reporting at 10 sentinel sites located in three ecological zones for the Malaria Therapeutic Efficacy Studies were used. The samples used were collected from 2007-2018 malaria transmission seasons and mutations in the genes were detected using PCR and Sanger sequencing.

Results

In all 1,142 samples were used for the study. For falcipain-2 gene (pffp2), Sanger sequencing was successful for 872 samples and were further analysed. The prevalence of the mutants was 45% (392/872) with pffp2 markers V51I and S59F occurring in 15.0% (128/872) and 3.0% (26/872) of the samples respectively. Prevalence of other P. falciparum gene mutations: coronin (pfcoronin) was 44.8% (37/90); cysteine desulfurase (pfnfs) was 73.9% (68/92); apicoplast ribosomal protein S10 (pfarps10) was 36.8% (35/95); ferredoxin (pffd) was 8.8% (8/91); multidrug resistance protein-1 (pfmrp1) was 95.2.0% (80/84); multidrug resistance protein-2 (pfmrp2) was 91.4% (32/35); dihydrofolate reductase (pfdhfr) was 99.0% (84/85); dihydropteroate synthase (pfdhps) was 72% (68/95).

Discussion

The observation of numerous mutations in these genes of interest in the Ghanaian isolates, some of which have been implicated in delayed parasite clearance is of great interest. The presence of these genotypes may account for the decline in the efficacies of ACT regimens being used to treat uncomplicated malaria in the country. The need for continuous monitoring of these genetic markers to give first-hand information on parasite susceptibility to antimalarial drugs to inform policy makers and stakeholders in malaria elimination in the country is further discussed.

Implementing landscape genetics in molecular epidemiology to determine drivers of vector-borne disease: A malaria case study

This study employs landscape genetics to investigate the environmental drivers of a deadly vector-borne disease, malaria caused by Plasmodium falciparum, in a more spatially comprehensive manner than any previous work. With 1804 samples from 44 sites collected in western Kenya in 2012 and 2013, we performed resistance surface analysis to show that Lake Victoria acts as a barrier to transmission between areas north and south of the Winam Gulf. In addition, Mantel correlograms clearly showed significant correlations between genetic and geographic distance over short distances (less than 70 km). In both cases, we used an identity-by-state measure of relatedness tailored to find highly related individual parasites in order to focus on recent gene flow that is more relevant to disease transmission. To supplement these results, we performed conventional population genetics analyses, including Bayesian clustering methods and spatial ordination techniques. These analyses revealed some differentiation on the basis of geography and elevation and a cluster of genetic similarity in the lowlands north of the Winam Gulf of Lake Victoria. Taken as a whole, these results indicate low overall genetic differentiation in the Lake Victoria region, but with some separation of parasite populations north and south of the Winam Gulf that is explained by the presence of the lake as a geographic barrier to gene flow. We recommend similar landscape genetics analyses in future molecular epidemiology studies of vector-borne diseases to extend and contextualize the results of traditional population genetics.

The Multistage Antimalarial Compound Calxinin Perturbates P. falciparum Ca2+ Homeostasis by Targeting a Unique Ion Channel

Malaria elimination urgently needs novel antimalarial therapies that transcend resistance, toxicity, and high costs. Our multicentric international collaborative team focuses on developing multistage antimalarials that exhibit novel mechanisms of action. Here, we describe the design, synthesis, and evaluation of a novel multistage antimalarial compound, 'Calxinin'. A compound that consists of hydroxyethylamine (HEA) and trifluoromethyl-benzyl-piperazine. Calxinin exhibits potent inhibitory activity in the nanomolar range against the asexual blood stages of drug-sensitive (3D7), multidrug-resistant (Dd2), artemisinin-resistant (IPC4912), and fresh Kenyan field isolated Plasmodium falciparum strains. Calxinin treatment resulted in diminished maturation of parasite sexual precursor cells (gametocytes) accompanied by distorted parasite morphology. Further, in vitro liver-stage testing with a mouse model showed reduced parasite load at an IC50 of 79 nM. A single dose (10 mg/kg) of Calxinin resulted in a 30% reduction in parasitemia in mice infected with a chloroquine-resistant strain of the rodent parasite P. berghei. The ex vivo ookinete inhibitory concentration within mosquito gut IC50 was 150 nM. Cellular in vitro toxicity assays in the primary and immortalized human cell lines did not show cytotoxicity. A computational protein target identification pipeline identified a putative P. falciparum membrane protein (Pf3D7_1313500) involved in parasite calcium (Ca2+) homeostasis as a potential Calxinin target. This highly conserved protein is related to the family of transient receptor potential cation channels (TRP-ML). Target validation experiments showed that exposure of parasitized RBCs (pRBCs) to Calxinin induces a rapid release of intracellular Ca2+ from pRBCs; leaving de-calcinated parasites trapped in RBCs. Overall, we demonstrated that Calxinin is a promising antimalarial lead compound with a novel mechanism of action and with potential therapeutic, prophylactic, and transmission-blocking properties against parasites resistant to current antimalarials.

Population structure and diversity of Plasmodium falciparum in children with asymptomatic malaria living in different ecological zones of Ghana

BACKGROUND

Genetic diversity in Plasmodium falciparum populations can be used to describe the resilience and spatial distribution of the parasite in the midst of intensified intervention efforts. This study used microsatellite analysis to evaluate the genetic diversity and population dynamics of P. falciparum parasites circulating in three ecological zones of Ghana.

METHODS

A total of 1168 afebrile children aged between 3 to 13 years were recruited from five (5) Primary schools in 3 different ecological zones (Sahel (Tamale and Kumbungu), Forest (Konongo) and Coastal (Ada and Dodowa)) of Ghana. Asymptomatic malaria parasite carriage was determined using microscopy and PCR, whilst fragment analysis of 6 microsatellite loci was used to determine the diversity and population structure of P. falciparum parasites.

RESULTS

Out of the 1168 samples examined, 16.1 and 39.5% tested positive for P. falciparum by microscopy and nested PCR respectively. The genetic diversity of parasites in the 3 ecological zones was generally high, with an average heterozygosity (He) of 0.804, 0.787 and 0.608 the rainy (peak) season for the Sahel, Forest and Coastal zones respectively. The mean He for the dry (off-peak) season were 0.562, 0.693 and 0.610 for the Sahel, Forest and Coastal zones respectively. Parasites from the Forest zone were more closely related to those from the Sahel than from the Coastal zone, despite the Coastal zone being closer in physical distance to the Forest zone. The fixation indexes among study sites ranged from 0.049 to 0.112 during the rainy season and 0.112 to 0.348 during the dry season.

CONCLUSION

A large asymptomatic parasite reservoir was found in the school children during both rainy and dry seasons, especially those in the Forest and Sahel savannah zones where parasites were also found to be related compared to those from the Coastal zone. Further studies are recommended to understand why despite the roll out of several malaria interventions in Ghana, high transmission still persist.

Within-household clustering of genetically related Plasmodium falciparum infections in a moderate transmission area of Uganda

BACKGROUND

Evaluation of genetic relatedness of malaria parasites is a useful tool for understanding transmission patterns, but patterns are not easily detectable in areas with moderate to high malaria transmission. To evaluate the feasibility of detecting genetic relatedness in a moderate malaria transmission setting, relatedness of Plasmodium falciparum infections was measured in cohort participants from randomly selected households in the Kihihi sub-county of Uganda (annual entomological inoculation rate of 27 infectious bites per person).

METHODS

All infections detected via microscopy or Plasmodium-specific loop mediated isothermal amplification from passive and active case detection during August 2011-March 2012 were genotyped at 26 microsatellite loci, providing data for 349 samples from 230 participants living in 80 households. Pairwise genetic relatedness was calculated using identity by state (IBS).

RESULTS

As expected, genetic diversity was high (mean heterozygosity [He] = 0.73), and the majority (76.5 %) of samples were polyclonal. Despite the high genetic diversity, fine-scale population structure was detectable, with significant spatiotemporal clustering of highly related infections. Although the difference in malaria incidence between households at higher (mean 1127 metres) versus lower elevation (mean 1015 metres) was modest (1.4 malaria cases per person-year vs. 1.9 per person-year, respectively), there was a significant difference in multiplicity of infection (2.2 vs. 2.6, p = 0.008) and, more strikingly, a higher proportion of highly related infections within households (6.3 % vs. 0.9 %, p = 0.0005) at higher elevation compared to lower elevation.

CONCLUSIONS

Genetic data from a relatively small number of diverse, multiallelic loci reflected fine scale patterns of malaria transmission. Given the increasing interest in applying genetic data to augment malaria surveillance, this study provides evidence that genetic data can be used to inform transmission patterns at local spatial scales even in moderate transmission areas.

Diversity and Multiplicity of P. falciparum infections among asymptomatic school children in Mbita, Western Kenya

Multiplicity of infection (MOI) and genetic diversity of P. falciparum infections are important surrogate indicators for assessing malaria transmission intensity in different regions of endemicity. Determination of MOI and diversity of P. falciparum among asymptomatic carriers will enhance our understanding of parasite biology and transmission to mosquito vectors. This study examined the MOI and genetic diversity of P. falciparum parasite populations circulating in Mbita, a region characterized as one of the malaria hotspots in Kenya. The genetic diversity and multiplicity of P. falciparum infections in 95 asymptomatic school children (age 5–15 yrs.) residing in Mbita, western Kenya were assessed using 10 polymorphic microsatellite markers. An average of 79.69% (Range: 54.84–95.74%) of the isolates analysed in this study were polyclonal infections as detected in at least one locus. A high mean MOI of 3.39 (Range: 2.24–4.72) and expected heterozygosity (He) of 0.81 (Range: 0.57–0.95) was reported in the study population. The analysed samples were extensively polyclonal infections leading to circulation of highly genetically diverse parasite populations in the study area. These findings correlated with the expectations of high malaria transmission intensity despite scaling up malaria interventions in the area thereby indicating the need for a robust malaria interventions particularly against asymptomatic carriers in order to attain elimination in the region.

High-resolution micro-epidemiology of parasite spatial and temporal dynamics in a high malaria transmission setting in Kenya

Novel interventions that leverage the heterogeneity of parasite transmission are needed to achieve malaria elimination. To better understand spatial and temporal dynamics of transmission, we applied amplicon next-generation sequencing of two polymorphic gene regions (csp and ama1) to a cohort identified via reactive case detection in a high-transmission setting in western Kenya. From April 2013 to July 2014, we enrolled 442 symptomatic children with malaria, 442 matched controls, and all household members of both groups. Here, we evaluate genetic similarity between infected individuals using three indices: sharing of parasite haplotypes on binary and proportional scales and the L1 norm. Symptomatic children more commonly share haplotypes with their own household members. Furthermore, we observe robust temporal structuring of parasite genetic similarity and identify the unique molecular signature of an outbreak. These findings of both micro- and macro-scale organization of parasite populations might be harnessed to inform next-generation malaria control measures.

Here, Nelson et al. use amplicon next-generation sequencing of two P. falciparum polymorphic gene regions to investigate the genetic similarity of parasite populations across time and space in a pediatric cohort in Kenya. They identify both micro- and macro-scale structuring of malaria parasites in this high-transmission setting, which could inform future intervention strategies.

The Diversity, Multiplicity of Infection and Population Structure of P. falciparum Parasites Circulating in Asymptomatic Carriers Living in High and Low Malaria Transmission Settings of Ghana

Background: Diversity in Plasmodium falciparum poses a major threat to malaria control and elimination interventions. This study utilized 12 polymorphic microsatellite (MS) markers and the Msp2 marker to examine diversity, multiplicity of infection (MOI) as well as the population structure of parasites circulating in two sites separated by about 92 km and with varying malaria transmission intensities within the Greater Accra Region of Ghana. Methods: The diversity and MOI of P. falciparum parasites in 160 non-symptomatic volunteers living in Obom (high malaria transmission intensity) and Asutsuare (low malaria transmission intensity) aged between 8 and 60 years was determined using Msp2 genotyping and microsatellite analysis. Results: The prevalence of asymptomatic P. falciparum carriers as well as the parasite density of infections was significantly higher in Obom than in Asutsuare. Samples from Asutsuare and Obom were 100% and 65% clonal, respectively, by Msp2 genotyping but decreased to 50% and 5%, respectively, when determined by MS analysis. The genetic composition of parasites from Obom and Asutsuare were highly distinct, with parasites from Obom being more diverse than those from Asutsuare. Conclusion: Plasmodium falciparum parasites circulating in Obom are genetically more diverse and distinct from those circulating in Asutsuare. The MOI in samples from both Obom and Asutsuare increased when assessed by MS analysis relative to MSP2 genotyping. The TA40 and TA87 loci are useful markers for estimating MOI in high and low parasite prevalence settings.

Genetic diversity and population structure of Plasmodium falciparum in Kenyan-Ugandan border areas

Kenya has, in the last decade, made tremendous progress in the fight against malaria. Nevertheless, continued surveillance of the genetic diversity and population structure of Plasmodium falciparum is required to refine malaria control and to adapt and improve elimination strategies. Twelve neutral microsatellite loci were genotyped in 201 P. falciparum isolates obtained from the Kenyan-Ugandan border (Busia) and from two inland malaria-endemic sites situated in western (Nyando) and coastal (Msambweni) Kenya. Analyses were done to assess the genetic diversity (allelic richness and expected heterozygosity, [H_e ]), multilocus linkage disequilibrium ( I S A ) and population structure. A similarly high degree of genetic diversity was observed among the three parasite populations surveyed (mean H_e  = 0.76; P > 0.05). Except in Msambweni, random association of microsatellite loci was observed, indicating high parasite out-breeding. Low to moderate genetic structure (F_ST   = 0.022-0.076; P < 0.0001) was observed with only 5% variance in allele frequencies observed among the populations. This study shows that the genetic diversity of P. falciparum populations at the Kenyan-Ugandan border is comparable to the parasite populations from inland Kenya. In addition, high genetic diversity, panmixia and weak population structure in this study highlight the fitness of Kenyan P. falciparum populations to successfully withstand malaria control interventions.

Population Genetics and Molecular Epidemiology of Eukaryotes

Molecular epidemiology uses the distribution and organization of a pathogen's DNA to understand the distribution and determinants of disease. Since the biology of DNA for eukaryotic pathogens differs substantially from that of bacteria, the analytic approach to their molecular epidemiology can also differ. While many of the genotyping techniques presented earlier in this series, "Advances in Molecular Epidemiology of Infectious Diseases," can be applied to eukaryotes, the output must be interpreted in the light of how DNA is distributed from one generation to the next. In some cases, parasite populations can be evaluated in ways reminiscent of bacteria. They differ, however, when analyzed as sexually reproducing organisms, where all individuals are unique but the genetic composition of the population does not change unless a limited set of events occurs. It is these events (migration, mutation, nonrandom mating, selection, and genetic drift) that are of interest. At a given time, not all of them are likely to be equally important, so the list can easily be narrowed down to understand the driving forces behind the population as it is now and even what it will look like in the future. The main population characteristics measured to assess these events are differentiation and diversity, interpreted in the light of what is known about the population from observation. The population genetics of eukaryotes is important for planning and evaluation of control measures, surveillance, outbreak investigation, and monitoring of the development and spread of drug resistance. *This article is part of a curated collection.

Selective sweeps and genetic lineages of Plasmodium falciparum multi-drug resistance (pfmdr1) gene in Kenya

BACKGROUND

There are concerns that resistance to artemisinin-based combination therapy might emerge in Kenya and sub-Saharan Africa (SSA) in the same pattern as was with chloroquine and sulfadoxine-pyrimethamine. Single nucleotide polymorphisms (SNPs) in critical alleles of pfmdr1 gene have been associated with resistance to artemisinin and its partner drugs. Microsatellite analysis of loci flanking genes associated with anti-malarial drug resistance has been used in defining the geographic origins, dissemination of resistant parasites and identifying regions in the genome that have been under selection.

METHODS

This study set out to investigate evidence of selective sweep and genetic lineages in pfmdr1 genotypes associated with the use of artemether-lumefantrine (AL), as the first-line treatment in Kenya. Parasites (n = 252) from different regions in Kenya were assayed for SNPs at codons 86, 184 and 1246 and typed for 7 neutral microsatellites and 13 microsatellites loci flanking (± 99 kb) pfmdr1 in Plasmodium falciparum infections.

RESULTS

The data showed differential site and region specific prevalence of SNPs associated with drug resistance in the pfmdr1 gene. The prevalence of pfmdr1 N86, 184F, and D1246 in western Kenya (Kisumu, Kericho and Kisii) compared to the coast of Kenya (Malindi) was 92.9% vs. 66.7%, 53.5% vs. to 24.2% and 96% vs. to 87.9%, respectively. The NFD haplotype which is consistent with AL selection was at 51% in western Kenya compared to 25% in coastal Kenya.

CONCLUSION

Selection pressures were observed to be different in different regions of Kenya, especially the western region compared to the coastal region. The data showed independent genetic lineages for all the pfmdr1 alleles. The evidence of soft sweeps in pfmdr1 observed varied in direction from one region to another. This is challenging for malaria control programs in SSA which clearly indicate effective malaria control policies should be based on the region and not at a country wide level.

Polymorphisms in the K13 Gene in Plasmodium falciparum from Different Malaria Transmission Areas of Kenya

The development of artemisinin (ART)-resistant parasites in Southeast Asia (SEA) threatens malaria control globally. Mutations in the Kelch 13 (K13)-propeller domain have been useful in identifying ART resistance in SEA. ART combination therapy (ACT) remains highly efficacious in the treatment of uncomplicated malaria in Sub-Saharan Africa (SSA). However, it is crucial that the efficacy of ACT is closely monitored. Toward this effort, this study profiled the prevalence of K13 nonsynonymous mutations in different malaria ecological zones of Kenya and in different time periods, before (pre) and after (post) the introduction of ACT as the first-line treatment of malaria. Nineteen nonsynonymous mutations were present in the pre-ACT samples (N = 64) compared with 22 in the post-ACT samples (N = 251). Eight of these mutations were present in both pre- and post-ACT parasites. Interestingly, seven of the shared single-nucleotide polymorphisms were at higher frequencies in the pre-ACT than the post-ACT parasites. The A578S mutation reported in SSA and the V568G mutation reported in SEA were found in both pre- and post-ACT parasites, with their frequencies declining post-ACT. D584Y and R539K mutations were found only in post-ACT parasites; changes in these codons have also been reported in SEA with different amino acids. The N585K mutation described for the first time in this study was present only in post-ACT parasites, and it was the most prevalent mutation at a frequency of 5.2%. This study showed the type, prevalence, and frequency of K13 mutations that varied based on the malaria ecological zones and also between the pre- and post-ACT time periods.

Population Genetic Analysis of Theileria annulata from Six Geographical Regions in China, Determined on the Basis of Micro- and Mini-satellite Markers

Theileria annulata, a tick-borne apicomplexan protozoan, causes a lymphoproliferative disease of cattle with high prevalence in tropical and sub-tropical regions. Understanding the genetic diversity and structure of local populations will provide more fundamental knowledge for the population genetics and epidemics of protozoa. In this study, 78 samples of T. annulata collected from cattle/yaks representing 6 different geographic populations in China were genotyped using eight micro- and mini-satellite markers. High genetic variation within population, moderate genetic differentiation, and high level of diversity co-occurring with significant linkage disequilibrium were observed, which indicates there is gene flow between these populations in spite of the existence of reproductive and geographical barriers among populations. Furthermore, some degree of genetic differentiation was also found between samples from China and Oman. These findings provide a first glimpse of the genetic diversity of the T. annulata populations in China, and might contribute to the knowledge of distribution, dynamics, and epidemiology of T. annulata populations and optimize the management strategies for control.

Geographic-genetic analysis of Plasmodium falciparum parasite populations from surveys of primary school children in Western Kenya

Background. Malaria control, and finally malaria elimination, requires the identification and targeting of residual foci or hotspots of transmission. However, the level of parasite mixing within and between geographical locations is likely to impact the effectiveness and durability of control interventions and thus should be taken into consideration when developing control programs.

Methods. In order to determine the geographic-genetic patterns of Plasmodium falciparum parasite populations at a sub-national level in Kenya, we used the Sequenom platform to genotype 111 genome-wide distributed single nucleotide polymorphic (SNP) positions in 2486 isolates collected from children in 95 primary schools in western Kenya. We analysed these parasite genotypes for genetic structure using principal component analysis and assessed local and global clustering using statistical measures of spatial autocorrelation. We further examined the region for spatial barriers to parasite movement as well as directionality in the patterns of parasite movement.

Results. We found no evidence of population structure and little evidence of spatial autocorrelation of parasite genotypes (correlation coefficients <0.03 among parasite pairs in distance classes of 1km, 2km and 5km; p value<0.01). An analysis of the geographical distribution of allele frequencies showed weak evidence of variation in distribution of alleles, with clusters representing a higher than expected number of samples with the major allele being identified for 5 SNPs. Furthermore, we found no evidence of the existence of spatial barriers to parasite movement within the region, but observed directional movement of parasites among schools in two separate sections of the region studied.

Conclusions. Our findings illustrate a pattern of high parasite mixing within the study region. If this mixing is due to rapid gene flow, then “one-off” targeted interventions may not be currently effective at the sub-national scale in Western Kenya, due to the high parasite movement that is likely to lead to re-introduction of infection from surrounding regions. However repeated targeted interventions may reduce transmission in the surrounding regions.

Geographic-genetic analysis of Plasmodium falciparum parasite populations from surveys of primary school children in Western Kenya

Background. Malaria control, and finally malaria elimination, requires the identification and targeting of residual foci or hotspots of transmission. However, the level of parasite mixing within and between geographical locations is likely to impact the effectiveness and durability of control interventions and thus should be taken into consideration when developing control programs. Methods. In order to determine the geographic-genetic patterns of Plasmodium falciparum parasite populations at a sub-national level in Kenya, we used the Sequenom platform to genotype 111 genome-wide distributed single nucleotide polymorphic (SNP) positions in 2486 isolates collected from children in 95 primary schools in western Kenya. We analysed these parasite genotypes for genetic structure using principal component analysis and assessed local and global clustering using statistical measures of spatial autocorrelation. We further examined the region for spatial barriers to parasite movement as well as directionality in the patterns of parasite movement. Results. We found no evidence of population structure and little evidence of spatial autocorrelation of parasite genotypes (correlation coefficients <0.03 among parasite pairs in distance classes of 1km, 2km and 5km; p value<0.01). An analysis of the geographical distribution of allele frequencies showed weak evidence of variation in distribution of alleles, with clusters representing a higher than expected number of samples with the major allele being identified for 5 SNPs. Furthermore, we found no evidence of the existence of spatial barriers to parasite movement within the region, but observed directional movement of parasites among schools in two separate sections of the region studied. Conclusions. Our findings illustrate a pattern of high parasite mixing within the study region. If this mixing is due to rapid gene flow, then “one-off” targeted interventions may not be currently effective at the sub-national scale in Western Kenya, due to the high parasite movement that is likely to lead to re-introduction of infection from surrounding regions. However repeated targeted interventions may reduce transmission in the surrounding regions.

Comparative assessment on the prevalence of mutations in the Plasmodium falciparum drug-resistant genes in two different ecotypes of Odisha state, India

Considering malaria as a local and focal disease, epidemiological understanding of different ecotypes of malaria can help in devising novel control measures. One of the major hurdles in malaria control lies on the evolution and dispersal of the drug-resistant malaria parasite, Plasmodium falciparum. We herewith present data on genetic variation at the Single Nucleotide Polymorphism (SNP) level in four different genes of P. falciparum (Pfcrt, Pfmdr1, Pfdhfr, and Pfdhps) that confer resistance to different antimalarials in two different eco-epidemiological settings, i.e. Hilly-Forest (HF) and Riverine-Plain (RP), in a high malaria endemic district of Odisha state, India. Greater frequency of antimalarial resistance conferring SNPs and haplotypes was observed in all four genes in P. falciparum, and Pfdhps was the most variable gene among the four. No significant genetic differentiation could be observed in isolates from HF and RP ecotypes. Twelve novel, hitherto unreported nucleotide mutations could be observed in the Pfmdr1 and Pfdhps genes. While the Pfdhps gene presented highest haplotype diversity, the Pfcrt gene displayed the highest nucleotide diversity. When the data on all the four genes were complied, the isolates from HF ecotype were found to harbour higher average nucleotide diversity than those coming from RP ecotype. High and positive Tajima's D values were obtained for the Pfcrt and Pfdhfr genes in isolates from both the HF and RP ecotypes, with statistically significant deviation from neutrality in the RP ecotype. Different patterns of Linkage Disequilibrium (LD) among SNPs located in different drug-resistant genes were found in the isolates collected from HF and RP ecotypes. Whereas in the HF ecotype, SNPs in the Pfmdr1 and Pfdhfr were significantly associated, in the RP ecotype, SNPs located in Pfcrt were associated with Pfmdr1, Pfdhfr and Pfdhps. These findings provide a baseline understanding on how different micro eco-epidemiological settings influence evolution and spread of different drug resistance alleles. Our findings further suggest that drug resistance to chloroquine and sulfadoxine-pyrimethamine is approaching fixation level, which requires urgent attention of malaria control programme in India.