DNA sequence polymorphisms of the pfmdr1 gene and association of mutations with the pfcrt gene in Indian Plasmodium falciparum isolates

Genetic diversity and evolution pattern of Plasmodium falciparum multidrug resistance 1 gene (Pfmdr1) in Indian and global populations

BACKGROUND

Among the crucial molecular markers contributing to multidrug resistance, the Plasmodium falciparum multidrug resistance-1 gene (Pfmdr1) remains understudied as compared to other drug-resistant genes in terms of its genetic diversity and evolution pattern. This study presents a comprehensive analysis of Pfmdr1 gene's genetic diversity aiming to discern its dynamics, distribution and evolutionary trends especially in Indian and global populations.

METHODS AND RESULTS

The Pfmdr1 gene was amplified and sequenced from 256 Plasmodium falciparum mono-infected samples collected from 14 Indian states during the years 1993-2023. Analysis revealed six non-synonymous (N86Y, N86F, S137C, D144F, F157L and Y184F) and one synonymous mutation (G182G) in N-terminal fragment. Among these, N86F, S137C, D144F and F157L were novel findings. The most prevalent mutations were N86Y (18.91%), Y184F (64.71%) and G182G (GGT > GGG) (59.24%; exclusive to India), with Y184F showing increasing trend when compared to N86Y over time. The mutation GGT > GGG is experiencing a hitchhiking by Y184F mutation which is likely undergoing a selective sweep. High haplotype and nucleotide diversity were observed in most Indian states, particularly in Odisha and Delhi. However, a decrease in diversity was noted in samples from 2020 onwards throughout India. Globally Pfmdr1 showed tendency of negative selection, except for populations from Liberia, Nigeria, Sudan and Central African Republic. Notably, samples from Sudan depicted a distinct haplotype and population structure compared to other countries.

CONCLUSIONS

These findings contribute significantly to our understanding of the genetic structure and evolutionary trends of Pfmdr1, which can help to strengthen the current malaria control policies for emergence of drug resistance.

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.

Single-Nucleotide Polymorphism Associates' β-Tubulin Isotype-1 Gene in Onchocerca volvulus Populations in Ivermectin-Treated Communities in Taraba State, Nigeria

PURPOSE

The occurrence of Single-Nucleotide Polymorphisms (SNPs) associated with repeated ivermectin treatment and sub-optimal responses reported by previous findings is of great concern in Onchocerciasis endemic areas. This study investigated SNPs' occurrence after 15 years of ivermectin intervention in Onchocerciasis endemic communities in two Local Government Areas of Taraba State, Nigeria.

METHODS

Microfilariae samples were collected by skin snip from individuals treated with ivermectin for 10-15 years of annual distribution and preserved in RNAlater^® in a 1.5 ml micro-centrifuge tube. Genomic DNA was extracted from microfilariae and residual skin, amplification in two regions within the β-tubulin gene, sequenced and analyzed for SNPs using Bioinformatics tools.

RESULTS

Three distinct SNP positions: 1183 (T/G), 1188 (T/C) and 1308 (C/T) on the β-tubulin gene on the targeted 1083-1568 bp fragment, associate's with the ivermectin-treated population. Furthermore, SNPs positions detected in this study are 1730 (A/G) and 1794 (T/G) in the β-tub gene in the 1557-1857 (bp) region. The 1794 (T/G) SNP position (Phe243Val) in the exon within the β-tubulin gene region were observed in this study.

CONCLUSION

The present study indicates that SNPs are observed in Onchocerca volvulus, thus strengthening the warning that genetic changes could occur in some parasite populations in some ivermectin-treated areas.

Characterization of drug resistance and genetic diversity of Plasmodium falciparum parasites from Tripura, Northeast India

Monitoring of anti-malarial drug resistance is vital in Northeast India as this region shares its international border with Southeast Asia. Genetic diversity of Plasmodium parasites regulates transmission dynamics, disease severity and vaccine efficacy. P. falciparum chloroquine resistance transporter (Pfcrt), multidrug resistance-1 (Pfmdr-1) and kelch 13 propeller (PfK-13) genes which govern antimalarial drug resistance and three genetic diversity markers, merozoite surface protein 1 and 2 (Pfmsp-1, Pfmsp-2) and glutamate rich protein (Pfglurp) were evaluated from Tripura, Northeast India using molecular tools. In the Pfcrt gene, 87% isolates showed triple mutations at codons M74I, N75E and K76T. 12.5% isolates in Pfmdr-1 gene showed mutation at N86Y. No polymorphism in PfK-13 propeller was found. Polyclonal infections were observed in 53.85% isolates and more commonly in adults (p = 0.0494). In the Pfmsp-1 locus, the K1 allelic family was predominant (71.2%) followed by the 3D7/IC family (69.2%) in the Pfmsp-2 locus. RII region of Pfglurp exhibited nine alleles with expected heterozygosity of 0.85. The multiplicity of infection for Pfmsp-1, Pfmsp-2 and Pfglurp were 1.56, 1.31 and 1.06 respectively. Overall, the study demonstrated a high level of chloroquine resistance and extensive parasite diversity in the region, necessitating regular surveillance in this population group.

Evolutionary interplay of single nucleotide polymorphisms at the promoter region of TNF-α gene in different clinical outcomes of malaria in India

Host genetic factors are frequently ascribed to differential malaria outcomes as a by-product of evolutionary adaptation. To this respect, Tumor Necrosis factor alpha (TNF-α), a human cytokine, is known to be associated with malaria through its differential regulation in diverse malaria manifestations. Since diversity in differential malaria outcome is uncommon in every endemic settings, possible association of TNF-α and malaria is not commonly established. In order to check for association between the occurrence of Single Nucleotide Polymorphisms (SNPs) in the TNF-α gene with different malaria manifestations, we have sequenced a 4011 bp region constituting the promoter and the whole gene of human TNF-α in 61 patients [(16 cerebral plus severe (SCM), 21 severe (SM) and 24 uncomplicated (UM)] samples in a highly malaria endemic state (Odisha) of India. Multiple sequence alignment revealed presence of six SNPs (-1031 T > C, -863C > A, -857C > T, -308G > A, -806C > T, +787C > A), out of which the -806C > T and +787C > A are novel in malaria patients in general and the +787C > A was detected for the first time in humans. Although alleles due to six different SNPs segregate differentially in the three groups of malaria (SCM, SM and UM) in the present study, interestingly, for the -1031 T > C position, the frequency of individuals possessing the homozygous rare allele was higher in the SCM group with a higher number of heterozygotes in the UM group. The Tajima's D values considering all the SNPs in a defined group were positive and statistically insignificant conforming no evolutionary constraint. However, statistically significant deviation from expectation under Hardy-Weinberg equilibrium for -1031 T > C SNP in the UM group points towards the probable role of natural selection providing some kind of protection to malaria in Odisha, India.

Prevalence of mutations linked to antimalarial resistance in Plasmodium falciparum from Chhattisgarh, Central India: A malaria elimination point of view

Antimalarial drug resistance is a major global challenge in malaria control and elimination. Mutations in six different genes of Plasmodium falciparum (crt, mdr1, dhfr, dhps, ATPase6 and K-13 propeller) that confer resistance to chloroquine, sulphadoxine-pyrimethamine and artemisinin-based combination therapy were analyzed in samples from Chhattisgarh. Seventy-eight percent of the samples were found to have a pfcrt mutation (53% double, 24% triple and 1% single mutant), and 59% of pfmdr1 genes were found to have an N86Y mutation. Double mutations were recorded in pfdhfr gene among 76% of the samples while only 6% of the samples harbored mutant genotypes in pfdhps. No mutation was found in the K-13 propeller gene, while only one sample showed a mutant genotype for the PfATPase6 gene. The Tajima test confirmed that there is no role of evolutionary natural selection in drug resistance, and gene pairwise linkage of disequilibrium showed significant intragenic association. The high level of pfcrt mutations suggests that parasite resistance to chloroquine is almost at a fixed level, whereas resistance to SP is evolving in the population and parasites remain sensitive to artemisinin derivatives. These findings provide potential information and understanding of the evolution and spread of different drug resistance alleles in Chhattisgarh.

Polymorphisms in Pfcrt and Pfmdr-1 genes after five years withdrawal of chloroquine for the treatment of Plasmodium falciparum malaria in West Bengal, India

BACKGROUND

The emergence of resistant power against different antimalarial agents particularly by Plasmodium falciparum is a challenge to combat malaria. Regular monitoring is essential not only to determine the efficacy and development of resistance by the parasite but also to detect early sign of regaining sensitivity to any anti-malarial agent that has been withdrawn for a long period. Studies on molecular markers associated with antimalarial drug resistance of prevailing Plasmodium population play an important role in this aspect. The present protocol was designed to study the polymorphisms in pfcrt and pfmdr-1 gene to determine any sign of regaining sensitivity to chloroquine among P. falciparum after five years of artemisinin combination therapy (ACT) implementation.

METHODS

Clinical isolates were collected from P. falciparum positive patients attending the malaria clinic of Calcutta School of Tropical Medicine during December 2014 to December 2015. Genomic parasitic DNA was extracted and subjected to sequencing of pfcrt and pfmdr-1 gene directly from purified PCR products.

RESULTS

A total of 89 isolates were sequenced for pfcrt and 73 isolates for pfmdr-1 genes. In pfcrt gene mutant K76T was detected in all isolates and all were SVMNT haplotype. Out of three important polymorphisms in pfmdr-1 gene mutant Y184F was detected among all isolates. One synonymous G182G and one non-synonymous S232F/Y, mutation were detected in 99% isolates.

CONCLUSION

All isolates carrying mutant K76T in pfcrt gene, considered as hall mark for CQ resistance, indicate that there is no sign of regaining CQ sensitivity among the prevailing P. falciparum population of the study area after five years of ACT implementation.

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.

Sequence analysis of pfcrt and pfmdr1 genes and its association with chloroquine resistance in Southeast Indian Plasmodium falciparum isolates

BACKGROUND

Due to the widespread resistance of Plasmodium falciparum to chloroquine drug, artemisinin-based combination therapy (ACT) has been recommended as the first-line treatment. This study aims to evaluate the extent of chloroquine resistance in P. falciparum infection after the introduction of ACT. This study was carried out based on the mutation analysis in P. falciparum chloroquine resistant transporter (pfcrt) and P. falciparum multidrug resistance 1 (pfmdr1) genes. Identification of these molecular markers plays a significant role in monitoring and assessment of drug resistance as well as in designing an effective antimalarial drug policy in India.

METHODS

Sixty blood samples were collected from patients infected with P. falciparum from JIPMER, Puducherry and MKCG Medical College, Odisha. Polymerase chain reaction-restriction fragment length polymorphism was performed, targeting the point mutation of K76T in pfcrt and N86Y in pfmdr1 gene. The PCR products were sequenced, genotyped and further analysed for amino acid changes in these codons.

RESULTS

The frequency of pfcrt mutation at 76th position was dominant for mutant T allele with 56.7% and wild type K, 43.3%. Majority of pfmdr1 86 allele were wild type, with N (90%) and mutant, Y (10%). Additionally, we found three haplotypes for CQ resistance, SVMNT, CVIET and CVIKT in association with the pfcrt gene. However, a poorly studied SNP in pfmdr1 gene (Y184F) associated with CQ resistance showed high frequency (70%) in P. falciparum isolates.

CONCLUSIONS

The point mutation K76T of pfcrt is high in P. falciparum suggesting a sustained high CQ resistance even after five years of the introduction of ACTs for antimalarial therapy. The present study suggests a strong association of CQ resistance with pfcrt T76, but not with the pfmdr1 Y86 mutation. However, sequence analysis showed that Y184F mutation on pfmdr1 gene was found to be associated with high resistance. Also, a new pfcrt haplotype 'CVIKT' associated with CQ resistance was found to be present in Indian strains of P. falciparum. The data obtained from this study helps in continuous monitoring of drug resistance in malaria and also suggests the need for careful usage of CQ in Plasmodium vivax malarial treatment.

Combinatorial Genetic Modeling of pfcrt-Mediated Drug Resistance Evolution in Plasmodium falciparum

The emergence of drug resistance continuously threatens global control of infectious diseases, including malaria caused by the protozoan parasite Plasmodium falciparum. A critical parasite determinant is the P. falciparum chloroquine resistance transporter (PfCRT), the primary mediator of chloroquine (CQ) resistance (CQR), and a pleiotropic modulator of susceptibility to several first-line artemisinin-based combination therapy partner drugs. Aside from the validated CQR molecular marker K76T, P. falciparum parasites have acquired at least three additional pfcrt mutations, whose contributions to resistance and fitness have been heretofore unclear. Focusing on the quadruple-mutant Ecuadorian PfCRT haplotype Ecu1110 (K76T/A220S/N326D/I356L), we genetically modified the pfcrt locus of isogenic, asexual blood stage P. falciparum parasites using zinc-finger nucleases, producing all possible combinations of intermediate pfcrt alleles. Our analysis included the related quintuple-mutant PfCRT haplotype 7G8 (Ecu1110 + C72S) that is widespread throughout South America and the Western Pacific. Drug susceptibilities and in vitro growth profiles of our combinatorial pfcrt-modified parasites were used to simulate the mutational trajectories accessible to parasites as they evolved CQR. Our results uncover unique contributions to parasite drug resistance and growth for mutations beyond K76T and predict critical roles for the CQ metabolite monodesethyl-CQ and the related quinoline-type drug amodiaquine in driving mutant pfcrt evolution. Modeling outputs further highlight the influence of parasite proliferation rates alongside gains in drug resistance in dictating successful trajectories. Our findings suggest that P. falciparum parasites have navigated constrained pfcrt adaptive landscapes by means of probabilistically rare mutational bursts that led to the infrequent emergence of pfcrt alleles in the field.