Molecular Epidemiology of P. vivax in Iran: High Diversity and Complex Sub-Structure Using Neutral Markers, but No Evidence of Y976F Mutation at pvmdr1

Molecular epidemiology of potential candidate markers for chloroquine resistance in imported Plasmodium vivax malaria cases in Iran

BACKGROUND

The spread of Plasmodium vivax strains resistant to chloroquine (CQ) has posed a challenge to control strategies aimed at eliminating malaria. Molecular analysis of candidate resistance markers is very important for monitoring the P. vivax resistance to CQ in different endemic regions. In the present study, the multidrug resistance 1 (pvmdr1) gene, a possible marker for CQ resistance in P. vivax, was evaluated by molecular methods.

METHODS

A simple PCR-RFLP method was developed for mutation analysis in pvmdr1 gene. A number of 120 blood spots were obtained from patients with P. vivax mono-infection in 2021. All of the samples were collected from Pakistani patients who travelled to Iran.

RESULTS

None of the samples had any mutation at codon 976 of pvmdr1, while the 1076 mutation was detected in 96.2% of the examined isolates. Only two pvmdr1 haplotypes were identified, including the single mutant (Y976/1076L) as the most prevalent haplotype (with 96.2% frequency) and the wild type (Y976/F1076; with 3.8% frequency).

CONCLUSIONS

In this study, the major CQ resistance-mediating mutation and multiple mutant haplotypes of the pvmdr1 gene was not detected. However, continuous monitoring of drug resistance markers and close supervision of the efficacy of CQ is essential to detect the potential emergence of CQ-resistant P. vivax isolates in Iran. This data is important for performing future epidemiological surveillance to monitor CQ resistance in this endemic area and the bordering regions.

Novel highly-multiplexed AmpliSeq targeted assay for Plasmodium vivax genetic surveillance use cases at multiple geographical scales

Although the power of genetic surveillance tools has been acknowledged widely, there is an urgent need in malaria endemic countries for feasible and cost-effective tools to implement in national malaria control programs (NMCPs) that can generate evidence to guide malaria control and elimination strategies, especially in the case of Plasmodium vivax. Several genetic surveillance applications (‘use cases’) have been identified to align research, technology development, and public health efforts, requiring different types of molecular markers. Here we present a new highly-multiplexed deep sequencing assay (Pv AmpliSeq). The assay targets the 33-SNP vivaxGEN-geo panel for country-level classification, and a newly designed 42-SNP within-country barcode for analysis of parasite dynamics in Vietnam and 11 putative drug resistance genes in a highly multiplexed NGS protocol with easy workflow, applicable for many different genetic surveillance use cases. The Pv AmpliSeq assay was validated using: 1) isolates from travelers and migrants in Belgium, and 2) routine collections of the national malaria control program at sentinel sites in Vietnam. The assay targets 229 amplicons and achieved a high depth of coverage (mean 595.7 ± 481) and high accuracy (mean error-rate of 0.013 ± 0.007). P. vivax parasites could be characterized from dried blood spots with a minimum of 5 parasites/µL and 10% of minority-clones. The assay achieved good spatial specificity for between-country prediction of origin using the 33-SNP vivaxGEN-geo panel that targets rare alleles specific for certain countries and regions. A high resolution for within-country diversity in Vietnam was achieved using the designed 42-SNP within-country barcode that targets common alleles (median MAF 0.34, range 0.01-0.49. Many variants were detected in (putative) drug resistance genes, with different predominant haplotypes in the pvmdr1 and pvcrt genes in different provinces in Vietnam. The capacity of the assay for high resolution identity-by-descent (IBD) analysis was demonstrated and identified a high rate of shared ancestry within Gia Lai Province in the Central Highlands of Vietnam, as well as between the coastal province of Binh Thuan and Lam Dong. Our approach performed well in geographically differentiating isolates at multiple spatial scales, detecting variants in putative resistance genes, and can be easily adjusted to suit the needs in other settings in a country or region. We prioritize making this tool available to researchers and NMCPs in endemic countries to increase ownership and ensure data usage for decision-making and malaria policy.

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.

The changing epidemiology of Plasmodium vivax: Insights from conventional and novel surveillance tools

Sarah Auburn and co-authors discuss the unique biology and epidemiology of P. vivax and current evidence on conventional and new approaches to surveillance.

Implementing parasite genotyping into national surveillance frameworks: feedback from control programmes and researchers in the Asia-Pacific region

The Asia-Pacific region faces formidable challenges in achieving malaria elimination by the proposed target in 2030. Molecular surveillance of Plasmodium parasites can provide important information on malaria transmission and adaptation, which can inform national malaria control programmes (NMCPs) in decision-making processes. In November 2019 a parasite genotyping workshop was held in Jakarta, Indonesia, to review molecular approaches for parasite surveillance and explore ways in which these tools can be integrated into public health systems and inform policy. The meeting was attended by 70 participants from 8 malaria-endemic countries and partners of the Asia Pacific Malaria Elimination Network. The participants acknowledged the utility of multiple use cases for parasite genotyping including: quantifying the prevalence of drug resistant parasites, predicting risks of treatment failure, identifying major routes and reservoirs of infection, monitoring imported malaria and its contribution to local transmission, characterizing the origins and dynamics of malaria outbreaks, and estimating the frequency of Plasmodium vivax relapses. However, the priority of each use case varies with different endemic settings. Although a one-size-fits-all approach to molecular surveillance is unlikely to be applicable across the Asia-Pacific region, consensus on the spectrum of added-value activities will help support data sharing across national boundaries. Knowledge exchange is needed to establish local expertise in different laboratory-based methodologies and bioinformatics processes. Collaborative research involving local and international teams will help maximize the impact of analytical outputs on the operational needs of NMCPs. Research is also needed to explore the cost-effectiveness of genetic epidemiology for different use cases to help to leverage funding for wide-scale implementation. Engagement between NMCPs and local researchers will be critical throughout this process.

Molecular surveillance over 14 years confirms reduction of Plasmodium vivax and falciparum transmission after implementation of Artemisinin-based combination therapy in Papua, Indonesia

Genetic epidemiology can provide important insights into parasite transmission that can inform public health interventions. The current study compared long-term changes in the genetic diversity and structure of co-endemic Plasmodium falciparum and P. vivax populations. The study was conducted in Papua Indonesia, where high-grade chloroquine resistance in P. falciparum and P. vivax led to a universal policy of Artemisinin-based Combination Therapy (ACT) in 2006. Microsatellite typing and population genetic analyses were undertaken on available isolates collected between 2004 and 2017 from patients with uncomplicated malaria (n = 666 P. falciparum and n = 615 P. vivax). The proportion of polyclonal P. falciparum infections fell from 28% (38/135) before policy change (2004-2006) to 18% (22/125) at the end of the study (2015-2017); p<0.001. Over the same period, polyclonal P. vivax infections fell from 67% (80/119) to 35% (33/93); p<0.001. P. falciparum strains persisted for up to 9 years compared to 3 months for P. vivax, reflecting higher rates of outbreeding in the latter. Sub-structure was observed in the P. falciparum population, but not in P. vivax, confirming different patterns of outbreeding. The P. falciparum population exhibited 4 subpopulations that changed in frequency over time. Notably, a sharp rise was observed in the frequency of a minor subpopulation (K2) in the late post-ACT period, accounting for 100% of infections in late 2016-2017. The results confirm epidemiological evidence of reduced P. falciparum and P. vivax transmission over time. The smaller change in P. vivax population structure is consistent with greater outbreeding associated with relapsing infections and highlights the need for radical cure to reduce recurrent infections. The study emphasizes the challenge in disrupting P. vivax transmission and demonstrates the potential of molecular data to inform on the impact of public health interventions.

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.

Genetic Diversity of Plasmodium vivax Causing Epidemic Malaria in the Republic of Korea

Plasmodium vivax is more challenging to control and eliminate than P. falciparum due to its more asymptomatic infections with low parasite densities making diagnosis more difficult, in addition to its unique biological characteristics. The potential re-introduction of incidence cases, either through borders or via human migrations, is another major hurdle to sustained control and elimination. The Republic of Korea has experienced re-emergence of vivax malaria in 1993 but is one of the 32 malaria-eliminating countries to-date. Despite achieving successful nationwide control and elimination of vivax malaria, the evolutionary characteristics of vivax malaria isolates in the Republic of Korea have not been fully understood. In this review, we present an overview of the genetic variability of such isolates to increase understanding of the epidemiology, diversity, and dynamics of vivax populations in the Republic of Korea.

Genomic analysis of a pre-elimination Malaysian Plasmodium vivax population reveals selective pressures and changing transmission dynamics

The incidence of Plasmodium vivax infection has declined markedly in Malaysia over the past decade despite evidence of high-grade chloroquine resistance. Here we investigate the genetic changes in a P. vivax population approaching elimination in 51 isolates from Sabah, Malaysia and compare these with data from 104 isolates from Thailand and 104 isolates from Indonesia. Sabah displays extensive population structure, mirroring that previously seen with the emergence of artemisinin-resistant P. falciparum founder populations in Cambodia. Fifty-four percent of the Sabah isolates have identical genomes, consistent with a rapid clonal expansion. Across Sabah, there is a high prevalence of loci known to be associated with antimalarial drug resistance. Measures of differentiation between the three countries reveal several gene regions under putative selection in Sabah. Our findings highlight important factors pertinent to parasite resurgence and molecular cues that can be used to monitor low-endemic populations at the end stages of P. vivax 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.

Nationwide genetic surveillance of Plasmodium vivax in Papua New Guinea reveals heterogeneous transmission dynamics and routes of migration amongst subdivided populations

The Asia Pacific Leaders in Malaria Alliance (APLMA) have committed to eliminate malaria from the region by 2030. Papua New Guinea (PNG) has the highest malaria burden in the Asia-Pacific region but with the intensification of control efforts since 2005, transmission has been dramatically reduced and Plasmodium vivax is now the dominant malaria infection in some parts of the country. To gain a better understanding of the transmission dynamics and migration patterns of P. vivax in PNG, here we investigate population structure in eight geographically and ecologically distinct regions of the country. A total of 219 P. vivax isolates (16-30 per population) were successfully haplotyped using 10 microsatellite markers. A wide range of genetic diversity (H_e=0.37-0.87, R_s=3.60-7.58) and significant multilocus linkage disequilibrium (LD) was observed in six of the eight populations (I_A^S=0.08-0.15 p-value<0.05) reflecting a spectrum of transmission intensities across the country. Genetic differentiation between regions was evident (Jost's D=0.07-0.72), with increasing divergence of populations with geographic distance. Overall, P. vivax isolates clustered into three major genetic populations subdividing the Mainland lowland and coastal regions, the Islands and the Highlands. P. vivax gene flow follows major human migration routes, and there was higher gene flow amongst Mainland parasite populations than among Island populations. The Central Province (samples collected in villages close to the capital city, Port Moresby), acts as a sink for imported infections from the three major endemic areas. These insights into P. vivax transmission dynamics and population networks will inform targeted strategies to contain malaria infections and to prevent the spread of drug resistance in PNG.

Malaria Epidemiology at the Clone Level

Genotyping to distinguish between parasite clones is nowadays a standard in many molecular epidemiological studies of malaria. It has become crucial in drug trials and to follow individual clones in epidemiological studies, and to understand how drug resistance emerges and spreads. Here, we review the applications of the increasingly available genotyping tools and whole-genome sequencing data, and argue for a better integration of population genetics findings into malaria-control strategies.

Chloroquine efficacy for Plasmodium vivax in Myanmar in populations with high genetic diversity and moderate parasite gene flow

BACKGROUND

Plasmodium vivax malaria remains a major public health burden in Myanmar. Resistance to chloroquine (CQ), the first-line treatment for P. vivax, has been reported in the country and has potential to undermine local control efforts.

METHODS

Patients over 6 years of age with uncomplicated P. vivax mono-infection were enrolled into clinical efficacy studies in Myawaddy in 2014 and Kawthoung in 2012. Study participants received a standard dose of CQ (25 mg/kg over 3 days) followed by weekly review until day 28. Pvmdr1 copy number (CN) and microsatellite diversity were assessed on samples from the patients enrolled in the clinical study and additional cross-sectional surveys undertaken in Myawaddy and Shwegyin in 2012.

RESULTS

A total of 85 patients were enrolled in the CQ clinical studies, 25 in Myawaddy and 60 in Kawthoung. One patient in Myawaddy (1.2%) had an early treatment failure and two patients (2.3%) in Kawthoung presented with late treatment failures on day 28. The day 28 efficacy was 92.0% (95% CI 71.6-97.9) in Myawaddy and 98.3% (95% CI 88.7-99.8) in Kawthoung. By day 2, 92.2% (23/25) in Myawaddy and 85.0% (51/60) in Kawthoung were aparasitaemic. Genotyping and pvmdr1 CN assessment was undertaken on 43, 52 and 46 clinical isolates from Myawaddy, Kawthoung and Shwegyin respectively. Pvmdr1 amplification was observed in 3.2% (1/31) of isolates in Myawaddy, 0% (0/49) in Kawthoung and 2.5% (1/40) in Shwegyin. Diversity was high in all sites (H _E 0.855-0.876), with low inter-population differentiation (F _ST 0.016-0.026, P < 0.05).

CONCLUSIONS

Treatment failures after chloroquine were observed following chloroquine monotherapy, with pvmdr1 amplification present in both Myawaddy and Shwegyin. The results emphasize the importance of ongoing P. vivax drug resistance surveillance in Myanmar, particularly given the potential connectivity between parasite population at different sites.

Genetic micro-epidemiology of malaria in Papua Indonesia: Extensive P. vivax diversity and a distinct subpopulation of asymptomatic P. falciparum infections

BACKGROUND

Genetic analyses of Plasmodium have potential to inform on transmission dynamics, but few studies have evaluated this on a local spatial scale. We used microsatellite genotyping to characterise the micro-epidemiology of P. vivax and P. falciparum diversity to inform malaria control strategies in Timika, Papua Indonesia.

METHODS

Genotyping was undertaken on 713 sympatric P. falciparum and P. vivax isolates from a cross-sectional household survey and clinical studies conducted in Timika. Standard population genetic measures were applied, and the data was compared to published data from Kalimantan, Bangka, Sumba and West Timor.

RESULTS

Higher diversity (HE = 0.847 vs 0.625; p = 0.017) and polyclonality (46.2% vs 16.5%, p<0.001) were observed in P. vivax versus P. falciparum. Distinct P. falciparum substructure was observed, with two subpopulations, K1 and K2. K1 was comprised solely of asymptomatic infections and displayed greater relatedness to isolates from Sumba than to K2, possibly reflecting imported infections.

CONCLUSIONS

The results demonstrate the greater refractoriness of P. vivax versus P. falciparum to control measures, and risk of distinct parasite subpopulations persisting in the community undetected by passive surveillance. These findings highlight the need for complimentary new surveillance strategies to identify transmission patterns that cannot be detected with traditional malariometric methods.

VivaxGEN: An open access platform for comparative analysis of short tandem repeat genotyping data in Plasmodium vivax populations

BACKGROUND

The control and elimination of Plasmodium vivax will require a better understanding of its transmission dynamics, through the application of genotyping and population genetics analyses. This paper describes VivaxGEN (http://vivaxgen.menzies.edu.au), a web-based platform that has been developed to support P. vivax short tandem repeat data sharing and comparative analyses.

RESULTS

The VivaxGEN platform provides a repository for raw data generated by capillary electrophoresis (FSA files), with fragment analysis and standardized allele calling tools. The query system of the platform enables users to filter, select and differentiate samples and alleles based on their specified criteria. Key population genetic analyses are supported including measures of population differentiation (FST), expected heterozygosity (HE), linkage disequilibrium (IAS), neighbor-joining analysis and Principal Coordinate Analysis. Datasets can also be formatted and exported for application in commonly used population genetic software including GENEPOP, Arlequin and STRUCTURE. To date, data from 10 countries, including 5 publicly available data sets have been shared with VivaxGEN.

CONCLUSIONS

VivaxGEN is well placed to facilitate regional overviews of P. vivax transmission dynamics in different endemic settings and capable to be adapted for similar genetic studies of P. falciparum and other organisms.