Molecular surveillance of pvdhfr, pvdhps, and pvmdr-1 mutations in Plasmodium vivax isolates from Yunnan and Anhui provinces of China

Surveillance of drug resistance molecular markers in Plasmodium vivax before and after introduction of dihydroartemisinin and piperaquine in Thailand: 2009-2019

Resistance to antimalarial drugs is a serious issue around the world. Widespread Plasmodium vivax and P. falciparum coinfections are commonly found in Thailand. Dihydroartemisinin and piperaquine (DHA-PPQ) have been used as first-line treatments for P. falciparum since 2015, and chloroquine (CQ) and primaquine (PQ) have remained first-line drugs for P. vivax for more than 60 years. Coinfections may lead parasites to evolve with regard to genetics under selective drug pressure. This study is aimed at investigating genes linked to antimalarial resistance in P. vivax before and after introduction of DHA-PPQ as a new drug regimen in Thailand. A total of 400 P. vivax isolates were collected from samples along the Thai-Myanmar and Thai-Malaysian borders before (2009-2015) and after (2016-2019) introduction of DHA-PPQ. Genomic DNA of P. vivax was obtained and subjected to analysis of five drug resistance-associated genes (Pvdhfr, Pvdhps, Pvmdr1, Pvcrt-o, and PvK12) by nested polymerase chain reaction (PCR), restriction fragment length polymorphism (RFLP), and nucleotide sequencing. A high prevalence of Pvdhfr was found in both endemic areas over the period. The quadruple (57I/58R/61M/117T) Pvdhfr haplotype was predominant in both periods in both endemic areas. Although the wild-type haplotype of Pvdhps was predominant in Thai-Malaysian isolates in both periods, a single mutant haplotype (383G) was dominant in Thai-Myanmar isolates during both periods. A low prevalence of the Pvmdr1 976F mutation was found in both periods among Thai-Myanmar isolates. A significant decrease in Pvmdr1 976F was identified in Thai-Malaysian isolates from the second period (p < 0.01). Only one nonsynonymous mutation of Pvcrt-o (193E) and one synonymous mutation of PvK12 (R584) were detected in four isolates (4.7%) and one isolate (0.5%) in the first period among Thai-Myanmar isolates, respectively. Thus, with limited clinical efficacy data, the low prevalence of drug-resistance markers may suggest that there is a low prevalence of P. vivax-resistant strains and that the current drug regimen for P. vivax is still effective for treating this P. vivax parasite population. Continued surveillance of antimalarial drug resistance markers and monitoring of clinical drug efficacy should be conducted for epidemiological and policy implications.

Global scenario of Plasmodium vivax occurrence and resistance pattern

Malaria caused by Plasmodium vivax is comparatively less virulent than Plasmodium falciparum, which can also lead to severe disease and death. It shows a wide geographical distribution. Chloroquine serves as a drug of choice, with primaquine as a radical cure. However, with the appearance of resistance to chloroquine and treatment has been shifted to artemisinin combination therapy followed by primaquine as a radical cure. Sulphadoxine-pyrimethamine, mefloquine, and atovaquone-proguanil are other drugs of choice in chloroquine-resistant areas, and later resistance was soon reported for these drugs also. The emergence of drug resistance serves as a major hurdle to controlling and eliminating malaria. The discovery of robust molecular markers and regular surveillance for the presence of mutations in malaria-endemic areas would serve as a helpful tool to combat drug resistance. Here, in this review, we will discuss the endemicity of P. vivax, a historical overview of antimalarial drugs, the appearance of drug resistance and molecular markers with their global distribution along with different measures taken to reduce malaria burden due to P. vivax infection and their resistance.

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

Prevalence of antifolate drug resistance markers in Plasmodium vivax in China

The dihydrofolate reductase (dhfr) and dihydropteroate synthetase (dhps) genes of Plasmodium vivax, as antifolate resistance-associated genes were used for drug resistance surveillance. A total of 375 P. vivax isolates collected from different geographical locations in China in 2009-2019 were used to sequence Pvdhfr and Pvdhps. The majority of the isolates harbored a mutant type allele for Pvdhfr (94.5%) and Pvdhps (68.2%). The most predominant point mutations were S117T/N (77.7%) in Pvdhfr and A383G (66.8%) in Pvdhps. Amino acid changes were identified at nine residues in Pvdhfr. A quadruple-mutant haplotype at 57, 58, 61, and 117 was the most frequent (57.4%) among 16 distinct Pvdhfr haplotypes. Mutations in Pvdhps were detected at six codons, and the double-mutant A383G/A553G was the most prevalent (39.3%). Pvdhfr exhibited a higher mutation prevalence and greater diversity than Pvdhps in China. Most isolates from Yunnan carried multiple mutant haplotypes, while the majority of samples from temperate regions and Hainan Island harbored the wild type or single mutant type. This study indicated that the antifolate resistance levels of P. vivax parasites were different across China and molecular markers could be used to rapidly monitor drug resistance. Results provided evidence for updating national drug policy and treatment guidelines.

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.

Molecular Surveillance and Ex Vivo Drug Susceptibilities of Plasmodium vivax Isolates From the China-Myanmar Border

Drug resistance in Plasmodium vivax may pose a challenge to malaria elimination. Previous studies have found that P. vivax has a decreased sensitivity to antimalarial drugs in some areas of the Greater Mekong Sub-region. This study aims to investigate the ex vivo drug susceptibilities of P. vivax isolates from the China–Myanmar border and genetic variations of resistance-related genes. A total of 46 P. vivax clinical isolates were assessed for ex vivo susceptibility to seven antimalarial drugs using the schizont maturation assay. The medians of IC₅₀ (half-maximum inhibitory concentrations) for chloroquine, artesunate, and dihydroartemisinin from 46 parasite isolates were 96.48, 1.95, and 1.63 nM, respectively, while the medians of IC₅₀ values for piperaquine, pyronaridine, mefloquine, and quinine from 39 parasite isolates were 19.60, 15.53, 16.38, and 26.04 nM, respectively. Sequence polymorphisms in pvmdr1 (P. vivax multidrug resistance-1), pvmrp1 (P. vivax multidrug resistance protein 1), pvdhfr (P. vivax dihydrofolate reductase), and pvdhps (P. vivax dihydropteroate synthase) were determined by PCR and sequencing. Pvmdr1 had 13 non-synonymous substitutions, of which, T908S and T958M were fixed, G698S (97.8%) and F1076L (93.5%) were highly prevalent, and other substitutions had relatively low prevalences. Pvmrp1 had three non-synonymous substitutions, with Y1393D being fixed, G1419A approaching fixation (97.8%), and V1478I being rare (2.2%). Several pvdhfr and pvdhps mutations were relatively frequent in the studied parasite population. The pvmdr1 G698S substitution was associated with a reduced sensitivity to chloroquine, artesunate, and dihydroartemisinin. This study suggests the possible emergence of P. vivax isolates resistant to certain antimalarial drugs at the China–Myanmar border, which demands continuous surveillance for drug resistance.

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.

Global assessment of genetic paradigms of Pvmdr1 mutations in chloroquine-resistant Plasmodium vivax isolates

BACKGROUND

Chloroquine (CQ) is generally prescribed as the front-line antimalarial drug of choice to treat Plasmodium vivax infections; however, some clinical CQ-resistant P. vivax isolates have been indigenously reported around the world during the last decade.

METHODS

In this study, P. vivax isolates (n=52) were obtained from autochthonous samples in southeast Iran during 2015-2017. The genomic DNA of samples was extracted, amplified (nested PCR) and sequenced by targeting the multidrug-resistance 1 gene. To verify the global genetic diversity of CQ-resistant P. vivax strains, the sequences of Pvmdr1 originating from Asia and the Americas were retrieved.

RESULTS

A total of 46 haplotypes were grouped into three distinct geographical haplogroups. The haplotype diversity and occurrence rates of Pvmdr1 976F/1076L mutations indicate that the efficacy of CQ is being compromised in Mexico, China, Nicaragua, Thailand, Brazil (2016), Ethiopia, Mauritania (2012) and southwest India in the near future. The cladistic phylogenetic tree showed that Pvmdr1 sequences isolated from the southeast Asian clade has a partial sister relationship with the American clade.

CONCLUSIONS

The current findings will serve as a basis to develop appropriate malaria control strategies and public health policies in symptomatic imported malaria cases or plausible CQ-resistant P. vivax strains.

Geographical spread and structural basis of sulfadoxine-pyrimethamine drug-resistant malaria parasites

The global spread of sulfadoxine (Sdx, S) and pyrimethamine (Pyr, P) resistance is attributed to increasing number of mutations in DHPS and DHFR enzymes encoded by malaria parasites. The association between drug resistance mutations and SP efficacy is complex. Here we provide an overview of the geographical spread of SP resistance mutations in Plasmodium falciparum (Pf) and Plasmodium vivax (Pv) encoded dhps and dhfr genes. In addition, we have collated the mutation data and mapped it on to the three-dimensional structures of DHPS and DHFR which have become available. Data from genomic databases and 286 studies were collated to provide a comprehensive landscape of mutational data from 2005 to 2019. Our analyses show that the Pyr-resistant double mutations are widespread in Pf/PvDHFR (P. falciparum ∼61% in Asia and the Middle East, and in the Indian sub-continent; in P. vivax ∼33% globally) with triple mutations prevailing in Africa (∼66%) and South America (∼33%). For PfDHPS, triple mutations dominate South America (∼44%), Asia and the Middle East (∼34%) and the Indian sub-continent (∼27%), while single mutations are widespread in Africa (∼45%). Contrary to the status for P. falciparum, Sdx-resistant single point mutations in PvDHPS dominate globally. Alarmingly, highly resistant quintuple and sextuple mutations are rising in Africa (PfDHFR-DHPS) and Asia (Pf/PvDHFR-DHPS). Structural analyses of DHFR and DHPS proteins in complexes with substrates/drugs have revealed that resistance mutations map proximal to Sdx and Pyr binding sites. Thus new studies can focus on discovery of novel inhibitors that target the non-substrate binding grooves in these two validated malaria parasite drug targets.

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.

Distribution pattern of amino acid mutations in chloroquine and antifolate drug resistance associated genes in complicated and uncomplicated Plasmodium vivax isolates from Chandigarh, North India

BACKGROUND

The increasing antimalarial drug resistance is a significant hindrance to malaria control and elimination programs. For the last six decades, chloroquine (CQ) plus pyrimethamine remains the first-line treatment for P. vivax malaria. Regions where both P. falciparum and P. vivax co-exist, P. vivax is exposed to antifolate drugs due to either misdiagnosis or improper treatment that causes selective drug pressure to evolve. Therefore, the present study aims to estimate antimalarial drug resistance among the complicated and uncomplicated P. vivax patients.

METHODS

A total of 143 P. vivax malaria positive patients were enrolled in this study, and DNA was isolated from their blood samples. Pvcrt-o, Pvmdr-1, Pvdhps, and Pvdhfr genes were PCRs amplified, and drug resistance-associated gene mutations were analyzed. Statistical analysis of the drug resistance genes and population diversity was performed using MEGA vs. 7.0.21 and DnaSP v software.

RESULTS

Among the CQ resistance marker gene Pvcrt-o, the prevalence of K10 insertion was 17.5% (7/40) and 9.5% (7/73) of complicated and uncomplicated P vivax group isolates respectively. In Pvmdr-1, double mutant haplotype (M₉₅₈/L₁₀₇₆) was found in 99% of the clinical isolates. Among the pyrimethamine resistance-associated gene Pvdhfr, the double mutant haplotype I₁₃P₃₃F₅₇R₅₈T₆₁N₁₁₇I₁₇₃ was detected in 23% (11/48) in complicated and 20% (17/85) in uncomplicated group isolates. In the sulphadoxine resistance-associated Pvdhps gene, limited polymorphism was observed with the presence of a single mutant (D459A) among 16 and 5% of the clinical isolates in the complicated and uncomplicated group respectively.

CONCLUSION

The study presents the situations of polymorphism in the antimalarial drug resistance-associated genes and emphasizes the need for regular surveillance. It is imperative for the development of suitable antimalarial drug policy in India.

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

BACKGROUND

In the Greater Mekong sub-region, Plasmodium vivax has become the predominant species and imposes a major challenge for regional malaria elimination. This study aimed to investigate the variations in genes potentially related to drug resistance in P. vivax populations from the China-Myanmar border area. In addition, this study also wanted to determine whether divergence existed between parasite populations associated with asymptomatic and acute infections.

METHODS

A total of 66 P. vivax isolates were obtained from patients with acute malaria who attended clinics at the Laiza area, Kachin State, Myanmar in 2015. In addition, 102 P. vivax isolates associated with asymptomatic infections were identified by screening of volunteers without signs or symptoms from surrounding villages. Slide-positive samples were verified with nested PCR detecting the 18S rRNA gene. Multiclonal infections were further excluded by genotyping at msp-3α and msp-3β genes. Parasite DNA from 60 symptomatic cases and 81 asymptomatic infections was used to amplify and sequence genes potentially associated with drug resistance, including pvmdr1, pvcrt-o, pvdhfr, pvdhps, and pvk12.

RESULTS

The pvmdr1 Y976F and F1076L mutations were present in 3/113 (2.7%) and 97/113 (85.5%) P. vivax isolates, respectively. The K10 insertion in pvcrt-o gene was found in 28.2% of the parasites. Four mutations in the two antifolate resistance genes reached relatively high levels of prevalence: pvdhfr S58R (53.4%), S117N/T (50.8%), pvdhps A383G (75.0%), and A553G (36.3%). Haplotypes with wild-type pvmdr1 (976Y/997K/1076F) and quadruple mutations in pvdhfr (13I/57L/58R/61M/99H/117T/173I) were significantly more prevalent in symptomatic than asymptomatic infections, whereas the pvmdr1 mutant haplotype 976Y/997K/1076L was significantly more prevalent in asymptomatic than symptomatic infections. In addition, quadruple mutations at codons 57, 58, 61 and 117 of pvdhfr and double mutations at codons 383 and 553 of pvdhps were found both in asymptomatic and symptomatic infections with similar frequencies. No mutations were found in the pvk12 gene.

CONCLUSIONS

Mutations in pvdhfr and pvdhps were prevalent in both symptomatic and asymptomatic P. vivax infections, suggestive of resistance to antifolate drugs. Asymptomatic carriers may act as a silent reservoir sustaining drug-resistant parasite transmission necessitating a rational strategy for malaria elimination in this region.

Prevalence of molecular markers associated with drug resistance of Plasmodium vivax isolates in Western Yunnan Province, China

BACKGROUND

Plasmodium vivax is the most widely distributed malaria parasite, and its drug resistance poses unique challenges to malaria elimination. The Greater Mekong Subregion (GMS) is known as the global epicenter of multidrug resistance. Surveillance of molecular markers associated with drug resistance in this area will help to inform drug policy.

METHODS

Dry blood spots from 58 patients out of 109 with P. vivax infection between 2017, December and 2019, March were obtained from Yingjiang County, Yunnan Province, along the China-Myanmar border. Pvdhfr, Pvdhps, Pvmdr1 and Pvcrt-o were amplified and sequenced to assess gene mutations. The polymorphism and prevalence of these molecular markers were analyzed.

RESULTS

Mutations in Pvdhfr at codons 57, 58, 61, 99 and 117 were detected in 27.59, 48.28, 27.59, 32.76 and 48.28% of the isolates, respectively. Single mutant haplotype (I13F57S58T61S99S117I173) was the most frequent (29.31%, 17/58), followed by double mutant haplotype (20.69%, 12/58). Of three types of tandem repeat variations of Pvdhfr, deletion type was the most common. Pvdhps showed a lower prevalence among mutation genotypes. Single mutant was dominant and accounted for 34.48% (20/58). Prevalence of Pvmdr1 mutations at codons 958 and 1076 were 100.00% and 84.48%, respectively. The proportion of double and single mutant types was 84.48% (49/58) and 15.52% (9/58), respectively. Eleven samples (18.97%, 11/58) showed K10 "AAG" insertion in chloroquine resistance transporter gene Pvcrt-o.

CONCLUSIONS

There was moderate diversity of molecular patterns of resistance markers of Pvdhfr, Pvdhps, Pvmdr1 and Pvcrt-o in imported P. vivax cases to Yingjiang county in Western Yunnan, along the China-Myanmar border. Prevalence and molecular pattern of candidate drug resistance markers Pvdhfr, Pvdhps, Pvmdr1 and Pvcrt-o were demonstrated in this current study, which would help to update drug policy.

Plasmodium vivax Malaria Viewed through the Lens of an Eradicated European Strain

The protozoan Plasmodium vivax is responsible for 42% of all cases of malaria outside Africa. The parasite is currently largely restricted to tropical and subtropical latitudes in Asia, Oceania, and the Americas. Though, it was historically present in most of Europe before being finally eradicated during the second half of the 20th century. The lack of genomic information on the extinct European lineage has prevented a clear understanding of historical population structuring and past migrations of P. vivax. We used medical microscope slides prepared in 1944 from malaria-affected patients from the Ebro Delta in Spain, one of the last footholds of malaria in Europe, to generate a genome of a European P. vivax strain. Population genetics and phylogenetic analyses placed this strain basal to a cluster including samples from the Americas. This genome allowed us to calibrate a genomic mutation rate for P. vivax, and to estimate the mean age of the last common ancestor between European and American strains to the 15th century. This date points to an introduction of the parasite during the European colonization of the Americas. In addition, we found that some known variants for resistance to antimalarial drugs, including Chloroquine and Sulfadoxine, were already present in this European strain, predating their use. Our results shed light on the evolution of an important human pathogen and illustrate the value of antique medical collections as a resource for retrieving genomic information on pathogens from the past.

Evolution of the Plasmodium vivax multidrug resistance 1 gene in the Greater Mekong Subregion during malaria elimination

Background

The malaria elimination plan of the Greater Mekong Subregion (GMS) is jeopardized by the increasing number of Plasmodium vivax infections and emergence of parasite strains with reduced susceptibility to the frontline drug treatment chloroquine/primaquine. This study aimed to determine the evolution of the P. vivax multidrug resistance 1 (Pvmdr1) gene in P. vivax parasites isolated from the China–Myanmar border area during the major phase of elimination.

Methods

Clinical isolates were collected from 275 P. vivax patients in 2008, 2012–2013 and 2015 in the China–Myanmar border area and from 55 patients in central China. Comparison was made with parasites from three border regions of Thailand.

Results

Overall, genetic diversity of the Pvmdr1 was relatively high in all border regions, and over the seven years in the China–Myanmar border, though slight temporal fluctuation was observed. Single nucleotide polymorphisms previously implicated in reduced chloroquine sensitivity were detected. In particular, M908L approached fixation in the China–Myanmar border area. The Y976F mutation sharply decreased from 18.5% in 2008 to 1.5% in 2012–2013 and disappeared in 2015, whereas F1076L steadily increased from 33.3% in 2008 to 77.8% in 2015. While neutrality tests suggested the action of purifying selection on the pvmdr1 gene, several likelihood-based algorithms detected positive as well as purifying selections operating on specific amino acids including M908L, T958M and F1076L. Fixation and selection of the nonsynonymous mutations are differently distributed across the three border regions and central China. Comparison with the global P. vivax populations clearly indicated clustering of haplotypes according to geographic locations. It is noteworthy that the temperate-zone parasites from central China were completely separated from the parasites from other parts of the GMS.

Conclusions

This study showed that P. vivax populations in the China–Myanmar border has experienced major changes in the Pvmdr1 residues proposed to be associated with chloroquine resistance, suggesting that drug selection may play an important role in the evolution of this gene in the parasite populations.

Molecular detection of antimalarial drug resistance in Plasmodium vivax from returned travellers to NSW, Australia during 2008-2018

To monitor drug resistance in Plasmodium vivax, a multidrug resistance 1 (Pvmdr1) gene and a putative transporter protein (Pvcrt-o) gene were used as molecular markers for chloroquine resistance. The biomarkers, the dihydrofolate reductase (Pvdhfr) gene and the dihydropteroate synthetase (Pvdhps) gene, were also used for the detection of resistance to sulphadoxine-pyrimethamine (SP); this drug is often accidentally used to treat P. vivax infections. Clinical blood samples (n = 120) were collected from patients who had been to one of eight malaria-endemic countries and diagnosed with P. vivax infection. The chloroquine resistance marker, the Pvmdr1 gene, showed F976:L1076 mutations and L1076 mutation. A K10 insertion in the Pvcrt-o gene was also found among the samples successfully sequenced. A combination of L/I57:R58:M61:T117 mutations in the Pvdhfr gene and G383:G553 mutations in the Pvdhps gene were also observed. Mutations found in these genes indicate that drug resistance is present in these eight countries. Whether or not countries are using chloroquine to treat P. vivax, there appears to be an increase in mutation numbers in resistance gene markers. The detected changes in mutation rates of these genes do suggest that there is still a trend towards increasing P. vivax resistance to chloroquine. The presence of the mutations associated with SP resistance indicates that P. vivax has had exposure to SP and this may be a consequence of either misdiagnosis or coinfections with P. falciparum in the past.

Molecular detection of drug resistant malaria in Southern Thailand

BACKGROUND

Drug resistance within the major malaria parasites Plasmodium vivax and Plasmodium falciparum threatens malaria control and elimination in Southeast Asia. Plasmodium vivax first-line treatment drug is chloroquine together with primaquine, and the first-line treatment for P. falciparum malaria is artemisinin in combination with a partner drug. Plasmodium vivax and P. falciparum parasites resistant to their respective first-line therapies are now found within Southeast Asia. The resistance perimeters may include high transmission regions of Southern Thailand which are underrepresented in surveillance efforts.

METHODS

This study investigated blood samples from malaria centres in Southern Thailand. Genetic loci associated with drug resistance were amplified and sequenced. Drug resistance associated genes Pvmdr1, Pvcrt-o, Pvdhfr, and Pvdhps were characterized for 145 cases of P. vivax malaria, as well as the artemisinin resistance-associated Pfkelch13 gene from 91 cases of P. falciparum malaria.

RESULTS

Plasmodium vivax samples from Southern Thai provinces showed numerous chloroquine and antifolate resistance-associated mutations, including SNP and Pvcrt-o K10-insertion combinations suggestive of chloroquine resistant P. vivax phenotypes. A high proportion of the C580Y coding mutation (conferring artemisinin resistance) was detected in P. falciparum samples originating from Ranong and Yala (where the mutation was previously unreported).

CONCLUSIONS

The results demonstrate a risk of chloroquine and antifolate resistant P. vivax phenotypes in Southern Thailand, and artemisinin resistant P. falciparum observed as far south as the Thai-Malaysian border region. Ongoing surveillance of antimalarial drug resistance markers is called for in Southern Thailand to inform case management.

Lack of quadruple and quintuple mutant alleles associated with sulfadoxine-pyrimethamine resistance in Plasmodium vivax isolates from Brazilian endemic areas

BACKGROUND AND OBJECTIVE

Brazil is responsible for a large number of Plasmodium vivax cases in America. Given the emergence of P. vivax parasites resistant to chloroquine and the effectiveness of antifolates in vivax malaria treatment together with a correlation between mutations in P. vivax dhfr and dhps genes and SP treatment failure, the point mutations in these genes were investigated.

METHODS

Blood samples from 54 patients experiencing vivax malaria symptomatic episodes in the Amazonian Region were investigated. Genomic DNA was extracted using a DNA extraction kit (QIAGENTM). Nested polymerase chain reaction (PCR) amplification was carried out followed by Sanger sequencing to detect single nucleotide polymorphisms (SNPs).

FINDINGS

All tested isolates showed non-synonymous mutations in pvdhfr gene: 117N (54/54, 100%) and 58R (25/54, 46%). Double mutant allele 58R/117N (FRTNI, 28%) was the most frequent followed by triple mutant alleles (58R/117N/173L, FRTNL, 11%; 58R/61M/117N, FRMNI, 5% 117N/173L, FSTNL, 4%) and quadruple mutant allele (58R/61M/117N/173L, FRMNL, 2%). A single mutation was observed at codon C383G in pvdhps gene (SGKAV, 48%).

CONCLUSION

No evidence of molecular signatures associated with P. vivax resistance to SP was observed in the Brazilian samples.

Plasmodium genomics: an approach for learning about and ending human malaria

Malaria causes high levels of morbidity and mortality in human beings worldwide. According to the World Health Organization (WHO), about half a million people die of this disease each year. Malaria is caused by six species of parasites belonging to the Plasmodium genus: P. falciparum, P. knowlesi, P. vivax, P. malariae, P. ovale curtisi, and P. ovale wallikeri. Currently, malaria is being kept under control with varying levels of elimination success in different countries. The development of new molecular tools as well as the use of next-generation sequencing (NGS) technologies and novel bioinformatic approaches has improved our knowledge of malarial epidemiology, diagnosis, treatment, vaccine development, and surveillance strategies. In this work, the genetics and genomics of human malarias have been analyzed. Since the first P. falciparum genome was sequenced in 2002, various population-level genetic and genomic surveys, together with transcriptomic and proteomic studies, have shown the importance of molecular approaches in supporting malaria elimination.

Genetic diversity of the Plasmodium vivax multidrug resistance 1 gene in Thai parasite populations

Plasmodium vivax resistance to chloroquine (CQ) was first reported over 60 years ago. Here we analyzed sequence variations in the multidrug resistance 1 gene (Pvmdr1), a putative molecular marker for P. vivax CQ resistance, in field isolates collected from three sites in Thailand during 2013-2016. Several single nucleotide polymorphisms previously implicated in reduced CQ sensitivity were found. These genetic variations encode amino acids in the two nucleotide-binding domains as well as the transmembrane domains of the protein. The high level of genetic diversity of Pvmdr1 provides insights into the evolutionary history of this gene. Specifically, there was little evidence of positive selection at amino acid F1076L in global isolates to be promoted as a possible marker for CQ resistance. Population genetic analysis clearly divided the parasites into eastern and western populations, which is consistent with their geographical separation by the central malaria-free area of Thailand. With CQ-primaquine remaining as the frontline treatment for vivax malaria in all regions of Thailand, such a population subdivision could be shaped and affected by the current drugs for P. falciparum since mixed P. falciparum/P. vivax infections often occur in this region.

Molecular evaluation of pvdhfr and pvmdr-1 mutants in Plasmodium vivax isolates after treatment with sulfadoxine/pyrimethamine and chloroquine in Iran during 2001-2016

The rising use of sulfadoxine/pyrimethamine (SP) in the treatment of chloroquine (CQ)-resistant Plasmodium falciparum has resulted in increased exposure to P. vivax isolates in Iran, where both species are being circulated. In this investigation, the frequency of pvdhfr and pvmdr-1 mutants was assessed in P. vivax strains during 2001-2016 after the introduction of SP/CQ in malarious areas of Iran. The P. vivax isolates (n, 52) were obtained from autochthonous samples in Southeast Iran during 2015-2016. The genomic DNA was extracted and examined using nested polymerase chain reaction-(PCR) and sequencing. Mutations were detected in pvdhfr codons P33L (21.2%), T61 M (25%), S93H (3.9%), and S117 T (1.9%) and 5 isolates showed double mutations (33 L/61 M, 7.7%; 33 L/117 T, 1.9%). No mutation was identified in pvdhfr codons F57 and S58. The pvmdr-1 1076 L mutation was detected in 93.3% of P. vivax isolates. The findings indicated that the frequency of three codons of pvdhfr F57/S58/S117 has decreased from 2001 (1.05%/7.0%/16.9%) to 2016 (0%/0%/1.9%). Genomic analysis of pvmdr-1 showed that the frequency of 1076 L has gradually increased from 2013 (93%) to 2016 (93.3%) (P > .05). The results demonstrated that P. vivax isolates are probably being exited under SP pressure, which reflects the appropriate level of training for field microscopists, as established by Iranian policymakers. Emergent pvdhfr codons 33L, 61M, and 93H should be noticed in plausible drug tolerance and treatment plans. The high prevalence of pvmdr-1 1076L mutation shows that efficacy of CQ combination with primaquine may be in danger of being compromised, however further investigations are needed to evaluate the clinical importance of CQ-resistant P. vivax isolates.

The prevalence of molecular markers of drug resistance in Plasmodium vivax from the border regions of Thailand in 2008 and 2014

The prevalence of Plasmodium vivax is increasing in the border regions of Thailand; one potential problem confounding the control of malaria in these regions is the emergence and spread of drug resistance. The aim of this study was to determine the genetic diversity in genes potentially linked to drug resistance in P. vivax parasites isolated from four different border regions of Thailand; Thai-Myanmar (Tak, Mae Hong Son and Prachuap Khiri Khan Provinces), and Thai-Cambodian borders (Chanthaburi Province). Isolates were collected from 345 P. vivax patients in 2008 and 2014, and parasite DNA extracted and subjected to nucleotide sequencing at five putative drug-resistance loci (Pvdhfr, Pvdhps, Pvmdr1, Pvcrt-o and Pvk12). The prevalence of mutations in Pvdhfr, Pvdhps and Pvmdr1 were markedly different between the Thai-Myanmar and Thai-Cambodian border areas and also varied between sampling times. All isolates carried the Pvdhfr (58R and 117N/T) mutation, however, whereas the quadruple mutant allele (I₅₇R₅₈M₆₁T₁₁₇) was the most prevalent (69.6%) in the Thai-Myanmar border region, the double mutant allele (F₅₇R₅₈T₆₁N₁₁₇) was at fixation on the Thai-Cambodian border (100%). The most prevalent genotypes of Pvdhps and Pvmdr1 were the double mutant (S₃₈₂G₃₈₃K₅₁₂G₅₅₃) (65.1%) and single mutant (M₉₅₈Y₉₇₆F₁₀₇₆) (46.5%) alleles, respectively on the Thai-Myanmar border while the single Pvdhps mutant (S₃₈₂G₃₈₃K₅₁₂A₅₅₃) (52.7%) and the triple Pvmdr1 mutant (M₉₅₈F₉₇₆L₁₀₇₆) (81%) alleles were dominant on the Thai-Cambodian border. No mutations were observed in the Pvcrt-o gene in either region. Novel mutations in the Pvk12 gene, the P. vivax orthologue of PfK13, linked to artemisinin resistance in Plasmodium falciparum, were observed with three nonsynonymous and three synonymous mutations in six isolates (3.3%).

Molecular monitoring of dihydrofolatereductase (dhfr) and dihydropteroatesynthetase (dhps) associated with sulfadoxine-pyrimethamine resistance in Plasmodium vivax isolates of Palawan, Philippines

The emergence of drug-resistant Plasmodium vivax poses problems for malaria control and elimination in some parts of the world, especially in developing countries where individuals are routinely exposed to the infection. The aim of this study was to determine the single nucleotide polymorphisms (SNPs) in dihydropteroate synthase (pvdhps) and dihydrofolate reductase (pvdhfr) genes associated with sulfadoxine-pyrimethamine (SP) drug resistance among P. vivax isolates collected in Palawan, Philippines. Genetic polymorphisms of pvdhps and pvdhfr were analysed by nested PCR. Analysis at specific codons I₁₃P₃₃F₅₇S₅₈T₆₁S₁₁₇I₁₇₃ associated with pyrimethamine resistance in the pvdhfr gene revealed that most of the samples (66/87, 75.9%) carried double mutation at positions I₁₃P₃₃F₅₇R₅₈T₆₁N₁₁₇I₁₇₃, while only 18.4% (16/87) of the isolates carried the wild-type haplotype (I₁₃P₃₃F₅₇S₅₈T₆₁S₁₁₇I₁₇₃). For the pvdhps gene, the codons involved in sulfadoxine resistance S₃₈₂A₃₈₃K₅₁₂A₅₅₃V₅₈₅ were investigated. Single mutation at position S₃₈₂G₃₈₃K₅₁₂A₅₅₃V₅₈₅ was most observed in 68.0% (68/100) of the samples, whereas wild-type haplotype was found in 26.0% (26/100) of samples. The pvdhps and pvdhfr combination S₃₈₂A₃₈₃K₅₁₂A₅₅₃V₅₈₅/I₁₃P₃₃F₅₇S₅₈T₆₁S₁₁₇I₁₇₃ (wild-type), S₃₈₂G₃₈₃K₅₁₂A₅₅₃V₅₈₅/I₁₃P₃₃F₅₇R₅₈T₆₁N₁₁₇I₁₇₃, and S₃₈₂A₃₈₃K₅₁₂A₅₅₃V₅₈₅-I₁₃P₃₃F₅₇R₅₈T₆₁N₁₁₇I₁₇₃ were the most frequently observed combination haplotypes from the three study sites. The information on molecular markers associated with antifolate drug-resistance could help better understanding ofthe molecular epidemiology and situation of SP resistant P. vivax malaria in the country. Continuous surveillance of these genetic markers is necessary to monitor the evolution of SP resistance in the Philippines.

Amplification of pfmdr1, pfcrt, pvmdr1, and K13 propeller polymorphisms associated with Plasmodium falciparum and Plasmodium vivax isolates from the China-Myanmar border

Malaria in the China-Myanmar border region is still severe; local transmission of both falciparum and vivax malaria persists, and there is a risk of geographically expanding antimalarial resistance. In this research, the pfmdr1, pfcrt, pvmdr1, and K13-propeller genotypes were determined in 26 Plasmodium falciparum and 64 Plasmodium vivax isolates from Yingjiang county of Yunnan province. The pfmdr1 (11.5%), pfcrt (34.6%), and pvmdr1 (3.1%) mutations were prevalent at the China-Myanmar border. The indigenous samples exhibited prevalences of 14.3%, 28.6%, and 14.3% for pfmdr1 N86Y, pfcrt K76T, and pfcrt M74I, respectively, whereas the samples from Myanmar showed prevalences of 10.5%, 21.1%, and 5.3%, respectively. The most prevalent genotypes of pfmdr1 and pfcrt were Y86Y184 and M74N75T76, respectively. No pvmdr1 mutation occurred in the indigenous samples but was observed in two cases coming from Myanmar. In addition, we are the first to report on 10 patients (38.5%) with five different K13 point mutations. The F446I allele is predominant (19.2%), and its prevalence was 28.6% in the indigenous samples of Yingjiang county and 15.8% in samples from Myanmar. The present data might be helpful for enrichment of the molecular surveillance of antimalarial resistance and useful for developing and updating guidance for the use of antimalarials in this region.