Molecular surveillance for drug resistance markers in Plasmodium vivax isolates from symptomatic and asymptomatic infections at the China-Myanmar border

Are pvcrt-o and pvmdr1 Gene Mutations Associated with Plasmodium vivax Chloroquine-Resistant Parasites?

(1) Background: Malaria remains a significant global public health issue. Since parasites quickly became resistant to most of the available antimalarial drugs, treatment effectiveness must be constantly monitored. In Brazil, up to 10% of cases of vivax malaria resistant to chloroquine (CQ) have been registered. Unlike P. falciparum, there are no definitive molecular markers for the chemoresistance of P. vivax to CQ. This work aimed to investigate whether polymorphisms in the pvcrt-o and pvmdr1 genes could be used as markers for assessing its resistance to CQ. (2) Methods: A total of 130 samples from P. vivax malaria cases with no clinical and/or parasitological evidence of CQ resistance were studied through polymerase chain reaction for gene amplification followed by target DNA sequencing. (3) Results: In the pvcrt-o exons, the K10 insert was present in 14% of the isolates. Regarding pvmdr1, T958M and F1076L haplotypes showed frequencies of 95% and 3%, respectively, while the SNP Y976F was not detected. (4) Conclusions: Since K10-pvcrt-o and F1076L/T958M-pvmdr1 polymorphisms were detected in samples from patients who responded well to CQ treatment, it can be concluded that mutations in these genes do not seem to have a potential for association with the phenotype of CQ resistance.

Molecular surveillance of chloroquine resistance in Plasmodium vivax isolates from malaria cases in Yunnan Province of China using pvcrt-o gene polymorphisms

BACKGROUND

The efficacy of chloroquine treatment for vivax malaria has been rarely evaluated due to a lack of an appropriate testing method. The objective of this study was to conduct molecular monitoring of chloroquine resistance in Plasmodium vivax strains from vivax malaria patients in Yunnan Province, focusing on the analysis of polymorphism in the P. vivax chloroquine resistance transporter protein orthologous gene (pvcrt-o).

METHODS

In accordance with the principles of a cohort study, blood samples were collected from malaria cases diagnosed with a P. vivax mono-infection in Yunnan Province from 2020 to 2022. Segmental PCR was used to amplify the whole pvcrt-o gene in the blood samples and their products were subsequently sequenced. The sequencing data were arranged to obtain the full coding DNA sequence (CDS) as well as the gene's promoter region sequences. The CDSs were aligned with the reference sequence (XM_001613407.1) of the P. vivax SalI isolate to identify the mutant loci.

RESULTS

From a total of 375 blood samples taken from vivax malaria cases, 272 both whole gene CDSs (1272-1275 bp) and promoter DNA sequences (707 bp) of pvcrt-o gene were obtained. Among the whole CDSs, there were 7 single nucleotide polymorphic sites in which c.7 A>G was the minor allele frequency (MAF) site with 4.4% (12/272) detection rate. The mutation detection rate showed a significant decrease from 9.8% (10/102) in 2020 to 1.1% (1/92) in 2021 and 1.3% (1/78) in 2022, indicating statistical significance (χ2 = 11.256, P < 0.05). Among the identified 12 haplotypes, the majority of which were wild type (75.7%; 206/272). These four mutant haplotypes (Hap_3, Hap_5, Hap_9, and Hap_10) were classified as "K10 insertion type" and accounted for 12.1% (33/272). The detection rate of Hap_3 increased from 1.0% (1/102) in 2020 to 13.0% (12/92) in 2021 and 14.1% (11/78) in 2022, indicating statistical significance. A total of 23.8% (65/272) of the samples exhibited 14 bp (bp) deletions in the promoter region, occurring most frequently in the wild type haplotype (Hap_1) samples at a rate of 28.6% (59/206).

CONCLUSIONS

In recent years in Yunnan Province, a notable proportion of vivax malaria patients are infected by P. vivax strains with a "K10 insertion" and partial sequence deletions in the promoter region of the pvcrt-o gene, necessitating vigilance.

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.

Characteristics of molecular markers associated with chloroquine resistance in Plasmodium vivax strains from vivax malaria cases in Yunnan Province, China

BACKGROUND

Chloroquine (CQ) has been the preferred clinical treatment for vivax malaria in Yunnan Province since 1958, with over 300,000 patients. This study aimed to help make trend predictions regarding variations the in anti-malarial drug susceptibility of Plasmodium vivax distributed in Yunnan Province and effectively implement monitoring measures on the efficacy of anti-malarial drugs for vivax malaria.

METHODS

Blood samples collected from patients with mono-P. vivax infections were employed in this study based on the principle of cluster sampling. The whole gene of P. vivax multidrug resistance 1 protein gene (pvmdr1) was amplified by nested-PCR techniques and the PCR amplification produce were sequenced by Sanger bidirectional sequencing. The mutant loci and haplotypes of coding DNA sequence (CDS) were identified by comparison with the reference sequence (NC_009915.1) of the P. vivax Sal I isolate. Parameters such as Ka/Ks ratio were calculated using MEGA 5.04 software.

RESULTS

A total of 753 blood samples from patients infected with mono-P. vivax were collected, of which 624 blood samples yielded the full gene sequence (4392 bp) of the pvmdr1 gene, with 283, 140, 119, and 82 sequences from 2014, 2020, 2021 and 2022, respectively. A total of 52 single nucleotide polymorphic (SNP) loci were detected for the 624 CDSs, of which 92.3% (48/52), 34.6% (18/52), 42.3% (22/52), and 36.5% (19/52) SNPs were detected in 2014, 2020, 2021 and 2022, respectively. All of 624 CDSs were defined for a total of 105 mutant haplotypes, with CDSs of 2014, 2020, 2021, and 2022 containing 88, 15, 21, and 13 haplotypes, respectively. Of the 105 haplotypes, the threefold mutant haplotype (Hap_87) was the starting point for stepwise evolution, and the most drastic tenfold mutations were Hap_14 and Hap_78, and the fivefold, sixfold, sevenfold, and eightfold mutations.

CONCLUSIONS

In the majority of vivax malaria cases in Yunnan Province, most of them were infected with strains carrying demonstrating highly mutated in pvmdr1 genes. However, the dominant mutation strains types varied from year to year, which warrants further exploration in order to confirm the correlation between with phenotypic changes in P. vivax strains and their susceptibility to anti-malarial drugs such as chloroquine.

Risk of hemolysis in Plasmodium vivax malaria patients receiving standard primaquine treatment in a population with high prevalence of G6PD deficiency

BACKGROUND

Primaquine is essential for the radical cure of Plasmodium vivax malaria, but it poses a potential danger of severe hemolysis in G6PD-deficient (G6PDd) patients. This study aimed to determine whether primaquine is safe in a population with high G6PD prevalence but lacking G6PD diagnosis capacity.

METHODS

In Myanmar, 152 vivax patients were gender- and age-matched at 1:3 for G6PDd versus G6PD-normal (G6PDn). Their risk of acute hemolysis was followed for 28 days after treatment with the standard chloroquine and 14-day primaquine (0.25 mg/kg/day) regimen.

RESULTS

Patients anemic and non-anemic at enrollment showed a rising and declining trend in the mean hemoglobin level, respectively. In males, the G6PDd group showed substantially larger magnitudes of hemoglobin reduction and lower hemoglobin nadir levels than the G6PDn group, but this trend was not evident in females. Almost 1/3 of the patients experienced clinically concerning declines in hemoglobin, with five requiring blood transfusion.

CONCLUSIONS

The standard 14-day primaquine regimen carries a significant risk of acute hemolytic anemia (AHA) in vivax patients without G6PD testing in a population with a high prevalence of G6PD deficiency and anemia. G6PD testing would avoid most of the clinically significant Hb reductions and AHA in male patients.

Spatio-temporal trends of malaria incidence from 2011 to 2017 and environmental predictors of malaria transmission in Myanmar

BACKGROUND

Myanmar bears the heaviest malaria burden in the Greater Mekong Subregion (GMS). This study assessed the spatio-temporal dynamics and environmental predictors of Plasmodium falciparum and Plasmodium vivax malaria in Myanmar.

METHODS

Monthly reports of malaria cases at primary health centers during 2011-2017 were analyzed to describe malaria distribution across Myanmar at the township and state/region levels by spatial autocorrelation (Moran index) and spatio-temporal clustering. Negative binomial generalized additive models identified environmental predictors for falciparum and vivax malaria, respectively.

RESULTS

From 2011 to 2017, there was an apparent reduction in malaria incidence in Myanmar. Malaria incidence peaked in June each year. There were significant spatial autocorrelation and clustering with extreme spatial heterogeneity in malaria cases and test positivity across the nation (P < 0.05). Areas with higher malaria incidence were concentrated along international borders. Primary clusters of P. falciparum persisted in western townships, while clusters of P. vivax shifted geographically over the study period. The primary cluster was detected from January 2011 to December 2013 and covered two states (Sagaing and Kachin). Annual malaria incidence was highest in townships with a mean elevation of 500‒600 m and a high variance in elevation (states with both high and low elevation). There was an apparent linear relationship between the mean normalized difference vegetative index and annual P. falciparum incidence (P < 0.05).

CONCLUSION

The decreasing trends reflect the significant achievement of malaria control efforts in Myanmar. Prioritizing the allocation of resources to high-risk areas identified in this study can achieve effective disease control.

Molecular Profiles of Multiple Antimalarial Drug Resistance Markers in Plasmodium falciparum and Plasmodium vivax in the Mandalay Region, Myanmar

Emergence and spreading of antimalarial drug resistant malaria parasites are great hurdles to combating malaria. Although approaches to investigate antimalarial drug resistance status in Myanmar malaria parasites have been made, more expanded studies are necessary to understand the nationwide aspect of antimalarial drug resistance. In the present study, molecular epidemiological analysis for antimalarial drug resistance genes in Plasmodium falciparum and P. vivax from the Mandalay region of Myanmar was performed. Blood samples were collected from patients infected with P. falciparum and P. vivax in four townships around the Mandalay region, Myanmar in 2015. Partial regions flanking major mutations in 11 antimalarial drug resistance genes, including seven genes (pfdhfr, pfdhps, pfmdr-1, pfcrt, pfk13, pfubp-1, and pfcytb) of P. falciparum and four genes (pvdhfr, pvdhps, pvmdr-1, and pvk12) of P. vivax were amplified, sequenced, and overall mutation patterns in these genes were analyzed. Substantial levels of mutations conferring antimalarial drug resistance were detected in both P. falciparum and P. vivax isolated in Mandalay region of Myanmar. Mutations associated with sulfadoxine-pyrimethamine resistance were found in pfdhfr, pfdhps, pvdhfr, and pvdhps of Myanmar P. falciparum and P. vivax with very high frequencies up to 90%. High or moderate levels of mutations were detected in genes such as pfmdr-1, pfcrt, and pvmdr-1 associated with chloroquine resistance. Meanwhile, low frequency mutations or none were found in pfk13, pfubp-1, pfcytb, and pvk12 of the parasites. Overall molecular profiles for antimalarial drug resistance genes in malaria parasites in the Mandalay region suggest that parasite populations in the region have substantial levels of mutations conferring antimalarial drug resistance. Continuous monitoring of mutations linked with antimalarial drug resistance is necessary to provide useful information for policymakers to plan for proper antimalarial drug regimens to control and eliminate malaria in the country.

Multidisciplinary Investigations of Sustained Malaria Transmission in the Greater Mekong Subregion

In the course of malaria elimination in the Greater Mekong Subregion (GMS), malaria epidemiology has experienced drastic spatiotemporal changes with residual transmission concentrated along international borders and the rising predominance of Plasmodium vivax. The emergence of Plasmodium falciparum parasites resistant to artemisinin and partner drugs renders artemisinin-based combination therapies less effective while the potential spread of multidrug-resistant parasites elicits concern. Vector behavioral changes and insecticide resistance have reduced the effectiveness of core vector control measures. In recognition of these problems, the Southeast Asian International Center of Excellence for Malaria Research (ICEMR) has been conducting multidisciplinary research to determine how human migration, antimalarial drug resistance, vector behavior, and insecticide resistance sustain malaria transmission at international borders. These efforts allow us to comprehensively understand the ecology of border malaria transmission and develop population genomics tools to identify and track parasite introduction. In addition to employing in vivo, in vitro, and molecular approaches to monitor the emergence and spread of drug-resistant parasites, we also use genomic and genetic methods to reveal novel mechanisms of antimalarial drug resistance of parasites. We also use omics and population genetics approaches to study insecticide resistance in malaria vectors and identify changes in mosquito community structure, vectorial potential, and seasonal dynamics. Collectively, the scientific findings from the ICEMR research activities offer a systematic view of the factors sustaining residual malaria transmission and identify potential solutions to these problems to accelerate malaria elimination in the GMS.

Genetic diversity of Plasmodium vivax reticulocyte binding protein 2b in global parasite populations

Background

Plasmodium vivax reticulocyte binding protein 2b (PvRBP2b) plays a critical role in parasite invasion of reticulocytes by binding the transferrin receptor 1. PvRBP2b is a vaccine candidate based on the negative correlation between antibody titers against PvRBP2b recombinant proteins and parasitemia and risk of vivax malaria. The aim of this study was to analyze the genetic diversity of the PvRBP2b gene in the global P. vivax populations.

Methods

Near full-length PvRBP2b nucleotide sequences (190–8349 bp) were obtained from 88 P. vivax isolates collected from the China–Myanmar border (n = 44) and Thailand (n = 44). An additional 224 PvRBP2b sequences were retrieved from genome sequences from parasite populations worldwide. The genetic diversity, neutral selection, haplotype distribution and genetic differentiation of PvRBP2b were examined.

Results

The genetic diversity of PvRBP2b was distributed unevenly, with peak diversity found in the reticulocyte binding region in the N-terminus. Neutrality analysis suggested that this region is subjected to balancing selection or population bottlenecks. Several amino acid variants were found in all or nearly all P. vivax endemic regions. However, the critical residues responsible for reticulocyte binding were highly conserved. There was substantial population differentiation according to the geographical separation. The distribution of haplotypes in the reticulocyte binding region varied among regions; even the two major haplotypes Hap_6 and Hap_8 were found in only five populations.

Conclusions

Our data show considerable genetic variations of PvRBPb in global parasite populations. The geographic divergence may pose a challenge to PvRBP2b-based vaccine development.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-022-05296-6.

Polymorphisms of potential drug resistant molecular markers in Plasmodium vivax from China–Myanmar border during 2008‒2017

Background

Plasmodium vivax remains the predominant species at the China–Myanmar border, imposing a major challenge to the recent gains in regional malaria elimination. To closely supervise the emerging of drug resistance in this area, we surveyed the variations in genes potentially correlated with drug resistance in P. vivax parasite and the possible drug selection with time.

Methods

A total of 235 P. vivax samples were collected from patients suffering uncomplicated malaria at Yingjiang, Tengchong, and Longling counties, and Nabang port in China, Yunnan province, and Laiza sub-township in Myanmar, from 2008 to 2017. Five potential drug resistance genes were amplified utilizing nested-PCR and analyzed, including pvdhfr, pvdhps, pvmdr1, pvcrt-o, and pvk12. The Pearson’s Chi-squared test or Fisher’s exact test were applied to determine the statistical frequency differences of mutations between categorical data.

Results

The pvdhfr F57I/L, S58R, T61M and S117T/N presented in 40.6%, 56.7%, 40.1%, and 56.0% of the sequenced P. vivax isolates, and these mutations significantly decreased with years. The haplotype formed by these quadruple mutations predominated in Yingjiang, Tengchong, Longling and Nabang. While a mutation H99S/R (56.6%) dominated in Laiza and increased with time. In pvdhps, the A383G prevailed in 69.2% of the samples, which remained the most prevalent haplotype. However, a significant decrease of its occurrence was also noticed over the time. The S382A/C and A553G existed in 8.4% and 30.8% of the isolates, respectively. In pvmdr1, the mutation Y976F occurred at a low frequency in 5/232 (2.2%), while T958M was fixed and F1076L was approaching fixed (72.4%). The K10 insertion was detected at an occurrence of 33.2% in pvcrt-o, whereas there was no significant difference among the sites or over the time. No mutation was identified in pvk12.

Conclusions

Mutations related with resistance to antifolate drugs are prevalent in this area, while their frequencies decrease significantly with time, suggestive of increased susceptibility of P. vivax parasite to antifolate drugs. Resistance to chloroquine (CQ) is possibly emerging. However, since the molecular mechanisms underneath CQ resistance is yet to be better understood, close supervision of clinical drug efficiency and continuous function investigation is urgently needed to alarm drug resistance.

Graphical abstract

Genetic polymorphisms in genes associated with drug resistance in Plasmodium vivax parasites from northeastern Myanmar

Background

Anti-malarial drug resistance is still a major threat to malaria elimination in the Great Mekong Sub-region. Plasmodium vivax parasites resistant to anti-malarial drugs are now found in Myanmar. Molecular surveillance on drug resistance genes in P. vivax parasites from northeastern Myanmar was aimed at estimating the underlying drug resistance in this region.

Methods

Blood samples from patients with vivax malaria were collected from Laiza city in northeastern Myanmar in 2020. Drug resistance genes including Pvcrt-o, Pvmdr1, Pvdhfr and Pvdhps were amplified and sequenced. Genetic polymorphisms and haplotypes were analysed to evaluate the prevalence of mutant alleles associated with drug resistance.

Results

A total of 149 blood samples from P. vivax patients were collected. The prevalence of Pvmdr1 mutations at codons 958 and 1076 was 100.0% and 52.0%, respectively, whereas no single nucleotide polymorphism was present at codon 976. The proportions of single and double mutant types were 48.0% and 52.0%, respectively. A K10 “AAG” insertion in the Pvcrt-o gene was not detected. Mutations in Pvdhfr at codons 57, 58, 61, 99 and 117 were detected in 29.9%, 54.3%, 27.6%, 44.9% and 55.1% of the samples, respectively. Wild type was predominant (46.3%), followed by quadruple and double mutant haplotypes. Of three types of tandem repeat variations of Pvdhfr, Type B, with three copies of GGDN repeats, was the most common. Pvdhps mutations were only detected at codons 383 and 553 and the wild type Pvdhps was dominant (78.0%). Eleven haplotypes were identified when combining the mutations of Pvdhfr and Pvdhps, among which the predominant one was the wild type (33.9%), followed by double mutant alleles S58R/S117N /WT (24.6%).

Conclusions

This study demonstrated resistant P. vivax phenotypes exists in northeastern Myanmar. Continued surveillance of drug resistance markers is needed to update treatment guidelines in this region.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12936-022-04084-y.

Polymorphism of Antifolate Drug Resistance in Plasmodium vivax From Local Residents and Migrant Workers Returned From the China-Myanmar Border

Drug-resistant Plasmodium vivax malaria impedes efforts to control, eliminate, and ultimately eradicate malaria in Southeast Asia. P. vivax resistance to antifolate drugs derives from point mutations in specific parasite genes, including the dihydropteroate synthase (pvdhps), dihydrofolate reductase (pvdhfr), and GTP cyclohydrolase I (pvgch1) genes. This study aims to investigate the prevalence and spread of drug resistance markers in P. vivax populating the China-Myanmar border. Blood samples were collected from symptomatic patients with acute P. vivax infection. Samples with single-clone P. vivax infections were sequenced for pvdhps and pvdhfr genes and genotyped for 6 flanking microsatellite markers. Copy number variation in the pvgch1 gene was also examined. Polymorphisms were observed in six different codons of the pvdhps gene (382, 383, 512, 549, 553, and 571) and six different codons of the pvdhfr gene (13, 57, 58, 61, 99, 117) in two study sites. The quadruple mutant haplotypes 57I/L/58R/61M/117T of pvdhfr gene were the most common (comprising 76% of cases in Myitsone and 43.7% of case in Laiza). The double mutant haplotype 383G/553G of pvdhps gene was also prevalent at each site (40.8% and 31%). Microsatellites flanking the pvdhfr gene differentiated clinical samples from wild type and quadruple mutant genotypes (F_ST= 0.259-0.3036), as would be expected for a locus undergoing positive selection. The lack of copy number variation of pvgch1 suggests that SP-resistant P. vivax may harbor alternative mechanisms to secure sufficient folate.

The molecular basis of antimalarial drug resistance in Plasmodium vivax

Plasmodium vivax is the most geographically widespread cause of human malaria and is responsible for the majority of cases outside of the African continent. While great progress has been made towards eliminating human malaria, drug resistant parasite strains pose a threat towards continued progress. Resistance has arisen to multiple antimalarials in P. vivax, including to chloroquine, which is currently the first line therapy for P. vivax in most regions. Despite its importance, an understanding of the molecular mechanisms of drug resistance in this species remains elusive, in large part due to the complex biology of P. vivax and the lack of in vitro culture. In this review, we will cover the extent and challenges of measuring clinical and in vitro drug resistance in P. vivax. We will consider the roles of candidate drug resistance genes. We will highlight the development of molecular approaches for studying P. vivax biology that provide the opportunity to validate the role of putative drug resistance mutations as well as identify novel mechanisms of drug resistance in this understudied parasite. Validated molecular determinants and markers of drug resistance are essential for the rapid and cost-effective monitoring of drug resistance in P. vivax, and will be useful for optimizing drug regimens and for informing drug policy in control and elimination settings.

•

Drug resistance is emerging in Plasmodium vivax, an important cause of malaria.

•

The complex biology of P. vivax and the limited range of research tools make it difficult to identify drug resistance.

•

The molecular mechanisms of drug resistance in P. vivax remain elusive.

•

This review highlights the extent of drug resistance, the putative mechanisms of resistance and new technologies for the study of P. vivax drug resistance.