Genomic Analysis of Plasmodium vivax in Southern Ethiopia Reveals Selective Pressures in Multiple Parasite Mechanisms

PvGAP: Development of a globally-applicable, highly-multiplexed microhaplotype amplicon panel for Plasmodium vivax

Plasmodium vivax malaria research has yet to fully benefit from the advances in genomic surveillance that have revolutionized P. falciparum epidemiology. Closing this gap is critical because genomic tools are necessary to achieve certain malaria control program objectives: 1) they allow monitoring of the spread of drug resistance and thus selection of therapeutic drugs based on local prevalence of resistance; 2) they permit classification of infections as local or imported, enabling more precise targeting of control resources; and 3) they can distinguish reinfection, recrudescence, and relapse, a necessity for conducting therapeutic efficacy studies. To achieve these objectives, microhaplotype marker panels that allow powerful genotyping of polyclonal infections are needed. A handful of such panels have been published for P. vivax, but they may be limited to certain geographic areas. We here present a Globally-applicable Amplicon Panel for P. vivax (PvGAP) designed to maximize discriminatory capability between geographic regions. PvGAP has 80 high diversity targets suitable for population genomics and eight targets of specific epidemiological interest, such as putative markers of drug resistance. We demonstrate PvGAP achieves robust amplification with field data and that it clearly distinguishes samples from different locations both at a regional and global scale. PvGAP is ready for broad application that can support powerful and comprehensive studies of malaria genomic epidemiology.

Spatial distribution of Plasmodium vivax Duffy Binding Protein copy number variation and Duffy genotype, and their association with parasitemia in Ethiopia

BACKGROUND

Duffy Binding Protein (PvDBP) binding to the Duffy antigen receptor for chemokine (DARC) is essential for Plasmodium vivax invasion of human reticulocytes. PvDBP copy number variation (CNV) might increase parasite invasion and thus parasitemia. We examined the spatial distribution of PvDBP CNVs and DARC genotypes and their association with parasitemia in P. vivax endemic settings in Ethiopia.

METHODOLOGY/PRINCIPAL FINDINGS

P. vivax isolates (n = 435) collected from five P. vivax endemic settings in Ethiopia were genotyped by amplifying the GATA1 transcription factor-binding site of the Duffy blood group and the CNV of PvDBP was quantified. Parasitemia was determined using 18S-based qPCR. The majority of participants were Duffy positive (96.8%, 421/435). Of the few Duffy negative individuals, most (n = 8) were detected from one site (Gondar). Multiple copies of PvDBP were detected in 83% (363/435) isolates with significant differences between sites (range 60%-94%). Both heterozygous (p = 0.005) and homozygous (p = 0.006) patients were more likely to have been infected by parasites with multiple PvDBP copies than Duffy negatives. Parasitemia was higher among the Duffy positives (median 17,218 parasites/µL; interquartile range [IQR] 2,895-104,489) than Duffy negatives (170; 78-24,132, p = 0.004) as well as in infections with 2 to 3 PvDBP copies (20,468; 3,649-110,632, p = 0.001) and more than 3 PvDBP copies (17,139; 2,831-95,946, p = 0.004) than single copy (5,673; 249-76,605).

CONCLUSIONS/SIGNIFICANCE

A high proportion of P. vivax infection was observed in Duffy positives in this study, yet few Duffy negatives were found infected with P. vivax. The significant prevalence of multi-copy PvDBP observed among Ethiopian P. vivax isolates explains the high prevalence and parasitemia observed in clinical cases. This suggests that vivax malaria is a public health concern in the country where the Duffy positive population predominates. Investigating the relative contribution to the maintenance of the infectious reservoir of infections with different genotyping backgrounds (both host and parasite) might be required.

High-throughput genotyping of Plasmodium vivax in the Peruvian Amazon via molecular inversion probes

Plasmodium vivax transmission occurs throughout the tropics and is an emerging threat in areas of Plasmodium falciparum decline, causing relapse infections that complicate treatment and control. Targeted sequencing for P. falciparum has been widely deployed to detect population structure and the geographic spread of antimalarial and diagnostic resistance. However, there are fewer such tools for P. vivax. Leveraging global variation data, we designed four molecular inversion probe (MIP) genotyping panels targeting geographically differentiating SNPs, neutral SNPs, putative antimalarial resistance genes, and vaccine candidate genes. We deployed these MIP panels on 866 infections from the Peruvian Amazon and identified transmission networks with clonality (IBD[identity by descent]>0.99), copy number variation in Pvdbp and multiple Pvrbps, mutations in antimalarial resistance orthologs, and balancing selection in 13 vaccine candidate genes. Our MIP panels are the broadest genotyping panel currently available and are poised for successful deployment in other regions of P. vivax transmission.

Quest for malaria management using natural remedies

Malaria, transmitted through the bite of a Plasmodium-infected Anopheles mosquito, remains a significant global health concern. This review examines the complex life cycle of Plasmodium, emphasizing the role of humans and mosquitoes in its transmission and proliferation. Malarial parasites are transmitted as sporozoites to the human body by biting an infected female Anopheles mosquito. These sporozoites then invade liver cells, multiply, and release merozoites, which infect red blood cells, perpetuating the cycle. As this cycle continues, the affected person starts experiencing the clinical symptoms of the disease. The current treatments for malaria, including chloroquine, artemisinin-based combination therapy, and quinine, are discussed alongside the challenges of drug resistance and misdiagnosis. Although efforts have been made to develop a malarial vaccine, they have so far been unsuccessful. Additionally, the review explores the potential of medicinal plants as remedies for malaria, highlighting the efficacy of compounds derived from Artemisia annua, Cinchona species, and Helianthus annuus L., as well as exploration of plants and phytocompounds like cryptolepine, and isoliquiritigenin against drug-resistant Plasmodium species. Moreover, studies from Pakistan further highlight the diverse vegetal resources utilized in malaria treatment, emphasizing the need for further research into natural remedies. Despite the advantages of herbal medicines, including cost-effectiveness, and fewer side effects; their limitations must be taken into account, including variations in potency and potential drug interactions. The review concludes by advocating for a balanced approach to malaria treatment and prevention, emphasizing the importance of early detection, accurate diagnosis, and integrated efforts to combat the disease in the endemic regions.

Analytic Approaches in Genomic Epidemiological Studies of Parasitic Protozoa

Whole genome sequencing (WGS) plays an important role in the advanced characterization of pathogen transmission and is widely used in studies of major bacterial and viral diseases. Although protozoan parasites cause serious diseases in humans and animals, WGS data on them are relatively scarce due to the large genomes and lack of cultivation techniques for some. In this review, we have illustrated bioinformatic analyses of WGS data and their applications in studies of the genomic epidemiology of apicomplexan parasites. WGS has been used in outbreak detection and investigation, studies of pathogen transmission and evolution, and drug resistance surveillance and tracking. However, comparative analysis of parasite WGS data is still in its infancy, and available WGS data are mainly from a few genera of major public health importance, such as Plasmodium, Toxoplasma, and Cryptosporidium. In addition, the utility of third-generation sequencing technology for complete genome assembly at the chromosome level, studies of the biological significance of structural genomic variation, and molecular surveillance of pathogens has not been fully exploited. These issues require large-scale WGS of various protozoan parasites of public health and veterinary importance using both second- and third-generation sequencing technologies.

Malaria Species Positivity Rates Among Symptomatic Individuals Across Regions of Differing Transmission Intensities in Mainland Tanzania

BACKGROUND

Recent data indicate that non-Plasmodium falciparum species may be more prevalent than thought in sub-Saharan Africa. Although Plasmodium malariae, Plasmodium ovale spp., and Plasmodium vivax are less severe than P. falciparum, treatment and control are more challenging, and their geographic distributions are not well characterized.

METHODS

We randomly selected 3284 of 12 845 samples collected from cross-sectional surveys in 100 health facilities across 10 regions of Mainland Tanzania and performed quantitative real-time PCR to determine presence and parasitemia of each malaria species.

RESULTS

P. falciparum was most prevalent, but P. malariae and P. ovale were found in all but 1 region, with high levels (>5%) of P. ovale in 7 regions. The highest P. malariae positivity rate was 4.5% in Mara and 8 regions had positivity rates ≥1%. We only detected 3 P. vivax infections, all in Kilimanjaro. While most nonfalciparum malaria-positive samples were coinfected with P. falciparum, 23.6% (n = 13 of 55) of P. malariae and 14.7% (n = 24 of 163) of P. ovale spp. were monoinfections.

CONCLUSIONS

P. falciparum remains by far the largest threat, but our data indicate that malaria elimination efforts in Tanzania will require increased surveillance and improved understanding of the biology of nonfalciparum species.

Imported malaria into Australia: surveillance insights and opportunities

BACKGROUND

Malaria continues to pose a significant burden in endemic countries, many of which lack access to molecular surveillance. Insights from malaria cases in travellers returning to non-endemic areas can provide valuable data to inform endemic country programmes. To evaluate the potential for novel global insights into malaria, we examined epidemiological and molecular data from imported malaria cases to Australia.

METHODS

We analysed malaria cases reported in Australia from 2012 to 2022 using National Notifiable Disease Surveillance System data. Molecular data on imported malaria cases were obtained from literature searches.

RESULTS

Between 2012 and 2022, 3204 malaria cases were reported in Australia. Most cases (69%) were male and 44% occurred in young adults aged 20-39 years. Incidence rates initially declined between 2012 and 2015, then increased until 2019. During 2012-2019, the incidence in travellers ranged from 1.34 to 7.71 per 100 000 trips. Cases were primarily acquired in Sub-Saharan Africa (n = 1433; 45%), Oceania (n = 569; 18%) and Southern and Central Asia (n = 367; 12%). The most common countries of acquisition were Papua New Guinea (n = 474) and India (n = 277). Plasmodium falciparum accounted for 58% (1871/3204) of cases and was predominantly acquired in Sub-Saharan Africa, and Plasmodium vivax accounted for 32% (1016/3204), predominantly from Oceania and Asia. Molecular studies of imported malaria cases to Australia identified genetic mutations and deletions associated with drug resistance and false-negative rapid diagnostic test results, and led to the establishment of reference genomes for P. vivax and Plasmodium malariae.

CONCLUSIONS

Our analysis highlights the continuing burden of imported malaria into Australia. Molecular studies have offered valuable insights into drug resistance and diagnostic limitations, and established reference genomes. Integrating molecular data into national surveillance systems could provide important infectious disease intelligence to optimize treatment guidelines for returning travellers and support endemic country surveillance programmes.

Population genomic evidence of structured and connected Plasmodium vivax populations under host selection in Latin America

Pathogen genomic epidemiology has the potential to provide a deep understanding of population dynamics, facilitating strategic planning of interventions, monitoring their impact, and enabling timely responses, and thereby supporting control and elimination efforts of parasitic tropical diseases. Plasmodium vivax, responsible for most malaria cases outside Africa, shows high genetic diversity at the population level, driven by factors like sub-patent infections, a hidden reservoir of hypnozoites, and early transmission to mosquitoes. While Latin America has made significant progress in controlling Plasmodium falciparum, it faces challenges with residual P. vivax. To characterize genetic diversity and population structure and dynamics, we have analyzed the largest collection of P. vivax genomes to date, including 1474 high-quality genomes from 31 countries across Asia, Africa, Oceania, and America. While P. vivax shows high genetic diversity globally, Latin American isolates form a distinctive population, which is further divided into sub-populations and occasional clonal pockets. Genetic diversity within the continent was associated with the intensity of transmission. Population differentiation exists between Central America and the North Coast of South America, vs. the Amazon Basin, with significant gene flow within the Amazon Basin, but limited connectivity between the Northwest Coast and the Amazon Basin. Shared genomic regions in these parasite populations indicate adaptive evolution, particularly in genes related to DNA replication, RNA processing, invasion, and motility - crucial for the parasite's survival in diverse environments. Understanding these population-level adaptations is crucial for effective control efforts, offering insights into potential mechanisms behind drug resistance, immune evasion, and transmission dynamics.

Genomic insights into Plasmodium vivax population structure and diversity in central Africa

BACKGROUND

Though Plasmodium vivax is the second most common malaria species to infect humans, it has not traditionally been considered a major human health concern in central Africa given the high prevalence of the human Duffy-negative phenotype that is believed to prevent infection. Increasing reports of asymptomatic and symptomatic infections in Duffy-negative individuals throughout Africa raise the possibility that P. vivax is evolving to evade host resistance, but there are few parasite samples with genomic data available from this part of the world.

METHODS

Whole genome sequencing of one new P. vivax isolate from the Democratic Republic of the Congo (DRC) was performed and used in population genomics analyses to assess how this central African isolate fits into the global context of this species.

RESULTS

Plasmodium vivax from DRC is similar to other African populations and is not closely related to the non-human primate parasite P. vivax-like. Evidence is found for a duplication of the gene PvDBP and a single copy of PvDBP2.

CONCLUSION

These results suggest an endemic P. vivax population is present in central Africa. Intentional sampling of P. vivax across Africa would further contextualize this sample within African P. vivax diversity and shed light on the mechanisms of infection in Duffy negative individuals. These results are limited by the uncertainty of how representative this single sample is of the larger population of P. vivax in central Africa.

Genomic analysis of Plasmodium vivax describes patterns of connectivity and putative drivers of adaptation in Ethiopia

Ethiopia has the greatest burden of Plasmodium vivax in Africa, but little is known about the epidemiological landscape of parasites across the country. We analysed the genomic diversity of 137 P. vivax isolates collected nine Ethiopian districts from 2012 to 2016. Signatures of selection were detected by cross-country comparisons with isolates from Thailand (n = 104) and Indonesia (n = 111), representing regions with low and high chloroquine resistance respectively. 26% (35/137) of Ethiopian infections were polyclonal, and 48.5% (17/35) of these comprised highly related clones (within-host identity-by-descent > 25%), indicating frequent co-transmission and superinfection. Parasite gene flow between districts could not be explained entirely by geographic distance, with economic and cultural factors hypothesised to have an impact on connectivity. Amplification of the duffy binding protein gene (pvdbp1) was prevalent across all districts (16-75%). Cross-population haplotype homozygosity revealed positive selection in a region proximal to the putative chloroquine resistance transporter gene (pvcrt-o). An S25P variant in amino acid transporter 1 (pvaat1), whose homologue has recently been implicated in P. falciparum chloroquine resistance evolution, was prevalent in Ethiopia (96%) but not Thailand or Indonesia (35-53%). The genomic architecture in Ethiopia highlights circulating variants of potential public health concern in an endemic setting with evidence of stable transmission.

Genomics of Plasmodium vivax in Colombia reveals evidence of local bottle-necking and inter-country connectivity in the Americas

Colombia aims to eliminate malaria by 2030 but remains one of the highest burden countries in the Americas. Plasmodium vivax contributes half of all malaria cases, with its control challenged by relapsing parasitaemia, drug resistance and cross-border spread. Using 64 Colombian P. vivax genomes collected between 2013 and 2017, we explored diversity and selection in two major foci of transmission: Chocó and Córdoba. Open-access data from other countries were used for comparative assessment of drug resistance candidates and to assess cross-border spread. Across Colombia, polyclonal infections were infrequent (12%), and infection connectivity was relatively high (median IBD = 5%), consistent with low endemicity. Chocó exhibited a higher frequency of polyclonal infections (23%) than Córdoba (7%), although the difference was not significant (P = 0.300). Most Colombian infections carried double pvdhfr (95%) and single pvdhps (71%) mutants, but other drug resistance mutations were less prevalent (< 10%). There was no evidence of selection at the pvaat1 gene, whose P. falciparum orthologue has recently been implicated in chloroquine resistance. Global population comparisons identified other putative adaptations. Within the Americas, low-level connectivity was observed between Colombia and Peru, highlighting potential for cross-border spread. Our findings demonstrate the potential of molecular data to inform on infection spread and adaptation.

No Association between the Plasmodium vivax crt-o MS334 or In9pvcrt Polymorphisms and Chloroquine Failure in a Pre-Elimination Clinical Cohort from Malaysia with a Large Clonal Expansion

Increasing reports of resistance to a frontline malaria blood-stage treatment, chloroquine (CQ), raises concerns for the elimination of Plasmodium vivax. The absence of an effective molecular marker of CQ resistance in P. vivax greatly constrains surveillance of this emerging threat. A recent genetic cross between CQ sensitive (CQS) and CQ resistant (CQR) NIH-1993 strains of P. vivax linked a moderate CQR phenotype with two candidate markers in P. vivax CQ resistance transporter gene (pvcrt-o): MS334 and In9pvcrt. Longer TGAAGH motif lengths at MS334 were associated with CQ resistance, as were shorter motifs at the In9pvcrt locus. In this study, high-grade CQR clinical isolates of P. vivax from a low endemic setting in Malaysia were used to investigate the association between the MS334 and In9pvcrt variants and treatment efficacy. Among a total of 49 independent monoclonal P. vivax isolates assessed, high-quality MS334 and In9pvcrt sequences could be derived from 30 (61%) and 23 (47%), respectively. Five MS334 and six In9pvcrt alleles were observed, with allele frequencies ranging from 2 to 76% and 3 to 71%, respectively. None of the clinical isolates had the same variant as the NIH-1993 CQR strain, and none of the variants were associated with CQ treatment failure (all P > 0.05). Multi-locus genotypes (MLGs) at 9 neutral microsatellites revealed a predominant P. vivax strain (MLG6) accounting for 52% of Day 0 infections. The MLG6 strain comprised equal proportions of CQS and CQR infections. Our study reveals complexity in the genetic basis of CQ resistance in the Malaysian P. vivax pre-elimination setting and suggests that the proposed pvcrt-o MS334 and In9pvcrt markers are not reliable markers of CQ treatment efficacy in this setting. Further studies are needed in other endemic settings, applying hypothesis-free genome-wide approaches, and functional approaches to understand the biological impact of the TGAAGH repeats linked to CQ response in a cross are warranted to comprehend and track CQR P. vivax.

Toxins from Animal Venoms as a Potential Source of Antimalarials: A Comprehensive Review

Malaria is an infectious disease caused by Plasmodium spp. and it is mainly transmitted to humans by female mosquitoes of the genus Anopheles. Malaria is an important global public health problem due to its high rates of morbidity and mortality. At present, drug therapies and vector control with insecticides are respectively the most commonly used methods for the treatment and control of malaria. However, several studies have shown the resistance of Plasmodium to drugs that are recommended for the treatment of malaria. In view of this, it is necessary to carry out studies to discover new antimalarial molecules as lead compounds for the development of new medicines. In this sense, in the last few decades, animal venoms have attracted attention as a potential source for new antimalarial molecules. Therefore, the aim of this review was to summarize animal venom toxins with antimalarial activity found in the literature. From this research, 50 isolated substances, 4 venom fractions and 7 venom extracts from animals such as anurans, spiders, scorpions, snakes, and bees were identified. These toxins act as inhibitors at different key points in the biological cycle of Plasmodium and may be important in the context of the resistance of Plasmodium to currently available antimalarial drugs.

Molecular investigation of malaria-infected patients in Djibouti city (2018-2021)

BACKGROUND

The Republic of Djibouti is a malaria endemic country that was in pre-elimination phase in 2006-2012. From 2013, however, malaria has re-emerged in the country, and its prevalence has been increasing every year. Given the co-circulation of several infectious agents in the country, the assessment of malaria infection based on microscopy or histidine-rich protein 2 (HRP2)-based rapid diagnostic tests (RDT) has shown its limitations. This study, therefore, aimed to assess the prevalence of malaria among febrile patients in Djibouti city using more robust molecular tools.

METHODS

All suspected malaria cases reported to be microscopy-positive were randomly sampled (n = 1113) and included in four health structures in Djibouti city over a 4-year period (2018-2021), mainly during the malaria transmission season (January-May). Socio-demographic information was collected, and RDT was performed in most of the included patients. The diagnosis was confirmed by species-specific nested polymerase chain reaction (PCR). Data were analysed using Fisher's exact test and kappa statistics.

RESULTS

In total, 1113 patients with suspected malaria and available blood samples were included. PCR confirmed that 788/1113 (70.8%) were positive for malaria. Among PCR-positive samples, 656 (83.2%) were due to Plasmodium falciparum, 88 (11.2%) Plasmodium vivax, and 44 (5.6%) P. falciparum/P. vivax mixed infections. In 2020, P. falciparum infections were confirmed by PCR in 50% (144/288) of negative RDTs. After the change of RDT in 2021, this percentage decreased to 17%. False negative RDT results were found more frequently (P < 0.05) in four districts of Djibouti city (Balbala, Quartier 7, Quartier 6, and Arhiba). Malaria occurred less frequently in regular bed net users than in non-users (odds ratio [OR]: 0.62, 95% confidence interval [CI]: 0.42-0.92).

CONCLUSIONS

The present study confirmed the high prevalence of falciparum malaria and, to a lesser extent, vivax malaria. Nevertheless, 29% of suspected malaria cases were misdiagnosed by microscopy and/or RDT. There is a need to strengthen the capacity for diagnosis by microscopy and to evaluate the possible role of P. falciparum hrp2 gene deletion, which leads to false negative cases of P. falciparum.

Population-based genomic study of Plasmodium vivax malaria in seven Brazilian states and across South America

Background

Brazil is a unique and understudied setting for malaria, with complex foci of transmission associated with human and environmental conditions. An understanding of the population genomic diversity of P. vivax parasites across Brazil can support malaria control strategies.

Methods

Through whole genome sequencing of P. vivax isolates across 7 Brazilian states, we use population genomic approaches to compare genetic diversity within country (n = 123), continent (6 countries, n = 315) and globally (26 countries, n = 885).

Findings

We confirm that South American isolates are distinct, have more ancestral populations than the other global regions, with differentiating mutations in genes under selective pressure linked to antimalarial drugs (pvmdr1, pvdhfr-ts) and mosquito vectors (pvcrmp3, pvP45/48, pvP47). We demonstrate Brazil as a distinct parasite population, with signals of selection including ABC transporter (PvABCI3) and PHIST exported proteins.

Interpretation

Brazil has a complex population structure, with evidence of P. simium infections and Amazonian parasites separating into multiple clusters. Overall, our work provides the first Brazil-wide analysis of P. vivax population structure and identifies important mutations, which can inform future research and control measures.

Funding

AI is funded by an MRC LiD PhD studentship. TGC is funded by the Medical Research Council (Grant no. MR/M01360X/1, MR/N010469/1, MR/R025576/1, MR/R020973/1 and MR/X005895/1). SC is funded by Medical Research Council UK grants (MR/M01360X/1, MR/R025576/1, MR/R020973/1 and MR/X005895/1) and Bloomsbury SET (ref. CCF17-7779). FN is funded by The Shloklo Malaria Research Unit - part of the Mahidol Oxford Research Unit, supported by the Wellcome Trust (Grant no. 220211). ARSB is funded by São Paulo Research Foundation - FAPESP (Grant no. 2002/09546-1). RLDM is funded by Brazilian National Council for Scientific and Technological Development - CNPq (Grant no. 302353/2003-8 and 471605/2011-5); CRFM is funded by FAPESP (Grant no. 2020/06747-4) and CNPq (Grant no. 302917/2019-5 and 408636/2018-1); JGD is funded by FAPESP fellowships (2016/13465-0 and 2019/12068-5) and CNPq (Grant no. 409216/2018-6).

A molecular barcode and web-based data analysis tool to identify imported Plasmodium vivax malaria

Traditionally, patient travel history has been used to distinguish imported from autochthonous malaria cases, but the dormant liver stages of Plasmodium vivax confound this approach. Molecular tools offer an alternative method to identify, and map imported cases. Using machine learning approaches incorporating hierarchical fixation index and decision tree analyses applied to 799 P. vivax genomes from 21 countries, we identified 33-SNP, 50-SNP and 55-SNP barcodes (GEO33, GEO50 and GEO55), with high capacity to predict the infection's country of origin. The Matthews correlation coefficient (MCC) for an existing, commonly applied 38-SNP barcode (BR38) exceeded 0.80 in 62% countries. The GEO panels outperformed BR38, with median MCCs > 0.80 in 90% countries at GEO33, and 95% at GEO50 and GEO55. An online, open-access, likelihood-based classifier framework was established to support data analysis (vivaxGEN-geo). The SNP selection and classifier methods can be readily amended for other use cases to support malaria control programs.

Plasmodium vivax: the potential obstacles it presents to malaria elimination and eradication

Initiatives to eradicate malaria have a good impact on P. falciparum malaria worldwide. P. vivax, however, still presents significant difficulties. This is due to its unique biological traits, which, in comparison to P. falciparum, pose serious challenges for malaria elimination approaches. P. vivax's numerous distinctive characteristics and its ability to live for weeks to years in liver cells in its hypnozoite form, which may elude the human immune system and blood-stage therapy and offer protection during mosquito-free seasons. Many malaria patients are not fully treated because of contraindications to primaquine use in pregnant and nursing women and are still vulnerable to P. vivax relapses, although there are medications that could radical cure P. vivax. Additionally, due to CYP2D6's highly variable genetic polymorphism, the pharmacokinetics of primaquine may be impacted. Due to their inability to metabolize PQ, some CYP2D6 polymorphism alleles can cause patients to not respond to treatment. Tafenoquine offers a radical treatment in a single dose that overcomes the potentially serious problem of poor adherence to daily primaquine. Despite this benefit, hemolysis of the early erythrocytes continues in individuals with G6PD deficiency until all susceptible cells have been eliminated. Field techniques such as microscopy or rapid diagnostic tests (RDTs) miss the large number of submicroscopic and/or asymptomatic infections brought on by reticulocyte tropism and the low parasitemia levels that accompany it. Moreover, P. vivax gametocytes grow more quickly and are much more prevalent in the bloodstream. P. vivax populations also have a great deal of genetic variation throughout their genome, which ensures evolutionary fitness and boosts adaptation potential. Furthermore, P. vivax fully develops in the mosquito faster than P. falciparum. These characteristics contribute to parasite reservoirs in the human population and facilitate faster transmission. Overall, no genuine chance of eradication is predicted in the next few years unless new tools for lowering malaria transmission are developed (i.e., malaria elimination and eradication). The challenging characteristics of P. vivax that impede the elimination and eradication of malaria are thus discussed in this article.

Plasmodium vivax blood stage invasion pathways: Contribution of omics technologies in deciphering molecular and cellular mechanisms

Vivax malaria is an infectious disease caused by Plasmodium vivax, a parasitic protozoan transmitted by female Anopheline mosquitoes. Historically, vivax malaria has often been regarded as a benign self-limiting infection due to the observation of low parasitemia in Duffy-positive patients in endemic transmission areas and the virtual absence of infections in Duffy-negative individuals in Sub Saharan Africa. However, the latest estimates show that the burden of the disease is not decreasing in many countries and cases of vivax infections in Duffy-negative individuals are increasingly reported throughout Africa. This raised questions about the accuracy of diagnostics and the evolution of interactions between humans and parasites. For a long time, our knowledge on P. vivax biology has been hampered due to the limited access to biological material and the lack of robust in vitro culture methods. Consequently, little is currently known about P. vivax blood stage invasion mechanisms. The introduction of omics technologies with novel and accessible techniques such as third generation sequencing and RNA sequencing at single cell level, two-dimensional electrophoresis, liquid chromatography, and mass spectrometry, has progressively improved our understanding of P. vivax genetics, transcripts, and proteins. This review aims to provide broad insights into P. vivax invasion mechanisms generated by genomics, transcriptomics, and proteomics and to illustrate the importance of integrated multi-omics studies.

Population genomics in neglected malaria parasites

Malaria elimination includes neglected human malaria parasites Plasmodium vivax, Plasmodium ovale spp., and Plasmodium malariae. Biological features such as association with low-density infection and the formation of hypnozoites responsible for relapse make their elimination challenging. Studies on these parasites rely primarily on clinical samples due to the lack of long-term culture techniques. With improved methods to enrich parasite DNA from clinical samples, whole-genome sequencing of the neglected malaria parasites has gained increasing popularity. Population genomics of more than 2200 P. vivax global isolates has improved our knowledge of parasite biology and host-parasite interactions, identified vaccine targets and potential drug resistance markers, and provided a new way to track parasite migration and introduction and monitor the evolutionary response of local populations to elimination efforts. Here, we review advances in population genomics for neglected malaria parasites, discuss how the rich genomic information is being used to understand parasite biology and epidemiology, and explore opportunities for the applications of malaria genomic data in malaria elimination practice.

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.

Geographical distribution and genetic diversity of Plasmodium vivax reticulocyte binding protein 1a correlates with patient antigenicity

Plasmodium vivax is the most widespread cause of human malaria. Recent reports of drug resistant vivax malaria and the challenge of eradicating the dormant liver forms increase the importance of vaccine development against this relapsing disease. P. vivax reticulocyte binding protein 1a (PvRBP1a) is a potential vaccine candidate, which is involved in red cell tropism, a crucial step in the merozoite invasion of host reticulocytes. As part of the initial evaluation of the PvRBP1a vaccine candidate, we investigated its genetic diversity and antigenicity using geographically diverse clinical isolates. We analysed pvrbp1a genetic polymorphisms using 202 vivax clinical isolates from six countries. Pvrbp1a was separated into six regions based on specific domain features, sequence conserved/polymorphic regions, and the reticulocyte binding like (RBL) domains. In the fragmented gene sequence analysis, PvRBP1a region II (RII) and RIII (head and tail structure homolog, 152-625 aa.) showed extensive polymorphism caused by random point mutations. The haplotype network of these polymorphic regions was classified into three clusters that converged to independent populations. Antigenicity screening was performed using recombinant proteins PvRBP1a-N (157-560 aa.) and PvRBP1a-C (606-962 aa.), which contained head and tail structure region and sequence conserved region, respectively. Sensitivity against PvRBP1a-N (46.7%) was higher than PvRBP1a-C (17.8%). PvRBP1a-N was reported as a reticulocyte binding domain and this study identified a linear epitope with moderate antigenicity, thus an attractive domain for merozoite invasion-blocking vaccine development. However, our study highlights that a global PvRBP1a-based vaccine design needs to overcome several difficulties due to three distinct genotypes and low antigenicity levels.

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.

Gene Polymorphisms Among Plasmodium vivax Geographical Isolates and the Potential as New Biomarkers for Gametocyte Detection

The unique biological features of Plasmodium vivax not only make it difficult to control but also to eliminate. For the transmission of the malaria parasite from infected human to the vector, gametocytes play a major role. The transmission potential of a malarial infection is inferred based on microscopic detection of gametocytes and molecular screening of genes in the female gametocytes. Microscopy-based detection methods could grossly underestimate the reservoirs of infection as gametocytes may occur as submicroscopic or as micro- or macro-gametocytes. The identification of genes that are highly expressed and polymorphic in male and female gametocytes is critical for monitoring changes not only in their relative proportions but also the composition of gametocyte clones contributing to transmission over time. Recent transcriptomic study revealed two distinct clusters of highly correlated genes expressed in the P. vivax gametocytes, indicating that the male and female terminal gametocytogeneses are independently regulated. However, the detective power of these genes is unclear. In this study, we compared genetic variations of 15 and 11 genes expressed, respectively, in the female and male gametocytes among P. vivax isolates from Southeast Asia, Africa, and South America. Further, we constructed phylogenetic trees to determine the resolution power and clustering patterns of gametocyte clones. As expected, Pvs25 (PVP01_0616100) and Pvs16 (PVP01_0305600) expressed in the female gametocytes were highly conserved in all geographical isolates. In contrast, genes including 6-cysteine protein Pvs230 (PVP01_0415800) and upregulated in late gametocytes ULG8 (PVP01_1452800) expressed in the female gametocytes, as well as two CPW-WPC family proteins (PVP01_1215900 and PVP01_1320100) expressed in the male gametocytes indicated considerably high nucleotide and haplotype diversity among isolates. Parasite samples expressed in male and female gametocyte genes were observed in separate phylogenetic clusters and likely represented distinct gametocyte clones. Compared to Pvs25, Pvs230 (PVP01_0415800) and a CPW-WPC family protein (PVP01_0904300) showed higher expression in a subset of Ethiopian P. vivax samples. Thus, Pvs230, ULG8, and CPW-WPC family proteins including PVP01_0904300, PVP01_1215900, and PVP01_1320100 could potentially be used as novel biomarkers for detecting both sexes of P. vivax gametocytes in low-density infections and estimating transmission reservoirs.

In vivo efficacy of anti-malarial drugs against clinical Plasmodium vivax malaria in Ethiopia: a systematic review and meta-analysis

Background

Ethiopia is one of the few countries in Africa where Plasmodium vivax commonly co-exists with Plasmodium falciparum, and which accounts for ~ 40% of the total number of malaria infections in the country. Regardless of the growing evidence over many decades of decreasing sensitivity of this parasite to different anti-malarial drugs, there has been no comprehensive attempt made to systematically review and meta-analyse the efficacy of different anti-malarial drugs against P. vivax in the country. However, outlining the efficacy of available anti-malarial drugs against this parasite is essential to guide recommendations for the optimal therapeutic strategy to use in clinical practice. The aim of this study was to synthesize evidence on the efficacy of anti-malarial drugs against clinical P. vivax malaria in Ethiopia.

Methods

All potentially relevant, peer-reviewed articles accessible in PubMed, Scopus, Web of Science, and Clinical Trial.gov electronic databases were retrieved using a search strategy combining keywords and related database-specific subject terms. Randomized controlled trials (RCTs) and non-randomized trials aiming to investigate the efficacy of anti-malarial drugs against P. vivax were included in the review. Data were analysed using Review Manager Software. Cochrane Q (χ²) and the I² tests were used to assess heterogeneity. The funnel plot and Egger’s test were used to examine risk of publication bias.

Results

Out of 1294 identified citations, 14 articles that presented data on 29 treatment options were included in the analysis. These studies enrolled 2144 clinical vivax malaria patients. The pooled estimate of in vivo efficacy of anti-malarial drugs against vivax malaria in Ethiopia was 97.91% (95% CI: 97.29–98.52%), with significant heterogeneity (I² = 86%, p < 0.0001) and publication bias (Egger’s test = -12.86, p < 0.001). Different anti-malarial drugs showed varied efficacies against vivax malaria. The duration of follow-up significantly affected the calculated efficacy of any given anti-malarial drug, with longer duration of the follow-up (42 days) associated with significantly lower efficacy than efficacy reported on day 28. Also, pooled PCR-corrected efficacy and efficacy estimated from altitudinally lower transmission settings were significantly higher than PCR-uncorrected efficacy that estimated for moderate transmission settings, respectively.

Conclusion

The overall efficacy of anti-malarial drugs evaluated for the treatment of vivax malaria in Ethiopia was generally high, although there was wide-ranging degree of efficacy, which was affected by the treatment options, duration of follow-up, transmission intensity, and the confirmation procedures for recurrent parasitaemia. Regardless of evidence of sporadic efficacy reduction reported in the country, chloroquine (CQ), the first-line regimen in Ethiopia, remained highly efficacious, supporting its continuous utilization for confirmed P. vivax mono-infections. The addition of primaquine (PQ) to CQ is recommended, as this is the only approved way to provide radical cure, and thus ensure sustained efficacy and longer protection against P. vivax. Continuous surveillance of the efficacy of anti-malarial drugs and clinical trials to allow robust conclusions remains necessary to proactively act against possible emergence and spread of drug-resistant P. vivax in Ethiopia.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12936-021-04016-2.

Longitudinal ex vivo and molecular trends of chloroquine and piperaquine activity against Plasmodium falciparum and P. vivax before and after introduction of artemisinin-based combination therapy in Papua, Indonesia

Drug resistant Plasmodium parasites are a major threat to malaria control and elimination. After reports of high levels of multidrug resistant P. falciparum and P. vivax in Indonesia, in 2005, the national first-line treatment policy for uncomplicated malaria was changed in March 2006, to dihydroartemisinin-piperaquine against all species. This study assessed the temporal trends in ex vivo drug susceptibility to chloroquine (CQ) and piperaquine (PIP) for both P. falciparum and P. vivax clinical isolates collected between 2004 and 2018, by using schizont maturation assays, and genotyped a subset of isolates for known and putative molecular markers of CQ and PIP resistance by using Sanger and next generation whole genome sequencing. The median CQ IC50 values varied significantly between years in both Plasmodium species, but there was no significant trend over time. In contrast, there was a significant trend for increasing PIP IC50s in both Plasmodium species from 2010 onwards. Whereas the South American CQ resistant 7G8 pfcrt SVMNT isoform has been fixed since 2005 in the study area, the pfmdr1 86Y allele frequencies decreased and became fixed at the wild-type allele in 2015. In P. vivax isolates, putative markers of CQ resistance (no pvcrt-o AAG (K10) insertion and pvmdr1 Y967F and F1076L) were fixed at the mutant alleles since 2005. None of the putative PIP resistance markers were detected in P. falciparum. The ex vivo drug susceptibility and molecular analysis of CQ and PIP efficacy for P. falciparum and P. vivax after 12 years of intense drug pressure with DHP suggests that whilst the degree of CQ resistance appears to have been sustained, there has been a slight decline in PIP susceptibility, although this does not appear to have reached clinically significant levels. The observed decreasing trend in ex vivo PIP susceptibility highlights the importance of ongoing surveillance.

Single-genome sequencing reveals within-host evolution of human malaria parasites

Population genomics of bulk malaria infections is unable to examine intrahost evolution; therefore, most work has focused on the role of recombination in generating genetic variation. We used single-cell sequencing protocol for low-parasitaemia infections to generate 406 near-complete single Plasmodium vivax genomes from 11 patients sampled during sequential febrile episodes. Parasite genomes contain hundreds of de novo mutations, showing strong signatures of selection, which are enriched in the ApiAP2 family of transcription factors, known targets of adaptation. Comparing 315 P. falciparum single-cell genomes from 15 patients with our P. vivax data, we find broad complementary patterns of de novo mutation at the gene and pathway level, revealing the importance of within-host evolution during malaria infections.

A population genetic perspective on the origin, spread and adaptation of the human malaria agents Plasmodium falciparum and Plasmodium vivax

Malaria is considered one of the most important scourges that humanity has faced during its history, being responsible every year for numerous deaths worldwide. The disease is caused by protozoan parasites, among which two species are responsible of the majority of the burden, Plasmodium falciparum and Plasmodium vivax. For these two parasite species, the questions of their origin (how and when they appeared in humans), of their spread throughout the world, as well as how they have adapted to humans have long been of interest to the scientific community. Here, we review the current knowledge that has accumulated on these different questions, thanks in particular to the analysis of the genetic and genomic variability of these parasites and comparison with related Plasmodium species infecting other host species (like non-human primates). In this paper we review the existing body of knowledge, including current research dealing with these questions, focusing particularly on genetic analysis and genomic variability of these parasites and comparison with related Plasmodium species infecting other species of host (such as non-human primates).

Advances and opportunities in malaria population genomics

Almost 20 years have passed since the first reference genome assemblies were published for Plasmodium falciparum, the deadliest malaria parasite, and Anopheles gambiae, the most important mosquito vector of malaria in sub-Saharan Africa. Reference genomes now exist for all human malaria parasites and nearly half of the ~40 important vectors around the world. As a foundation for genetic diversity studies, these reference genomes have helped advance our understanding of basic disease biology and drug and insecticide resistance, and have informed vaccine development efforts. Population genomic data are increasingly being used to guide our understanding of malaria epidemiology, for example by assessing connectivity between populations and the efficacy of parasite and vector interventions. The potential value of these applications to malaria control strategies, together with the increasing diversity of genomic data types and contexts in which data are being generated, raise both opportunities and challenges in the field. This Review discusses advances in malaria genomics and explores how population genomic data could be harnessed to further support global disease control efforts.

High Proportion of Genome-Wide Homology and Increased Pretreatment pvcrt Levels in Plasmodium vivax Late Recurrences: a Chloroquine Therapeutic Efficacy Study

Chloroquine (CQ) is the first-line treatment for Plasmodium vivax malaria in most countries where malaria is endemic. Monitoring P. vivax CQ resistance (CQR) is critical but remains challenged by the difficulty to distinguish real treatment failure from reinfection or liver relapse. The therapeutic efficacy of CQ against uncomplicated P. vivax malaria was evaluated in Gia Lai Province, Vietnam. Sixty-seven patients were enrolled and followed for 42 days using microscopy and quantitative PCR. Adequate clinical and parasitological response (ACPR) was 100% (66/66) on day 28 but 75.4% (49/65) on day 42. Eighteen recurrences (27.7%) were detected, with a median time to recurrence of 42 days (interquartile range [IQR], 35 to 42) and blood CQ concentration of <100 ng/ml. Primary infections leading to recurrence occurred in younger individuals (median age for ACPR = 25 years [IQR, 20 to 28]; recurrences = 18 [16 to 21]; P = 0.002) had a longer parasite clearance time (PCT for ACPR = 47.5 h [IQR, 36.2 to 59.8 h]; recurrences = 54.2 [48.4 to 62.0]; P = 0.035) and higher pvcrt gene expression (median relative expression ratio for ACPR = 0.09 [IQR, 0.05 to 0.22]; recurrences = 0.20 [0.15 to 0.56]; P = 0.002), but showed no differences in ex vivo CQ sensitivity. Parasite genotyping by microsatellites, single nucleotide polymorphism (SNP) barcoding, and whole-genome sequencing (WGS) identified a majority of homologous recurrences, with 80% (8/10) showing >98% identity by descent to paired day 0 samples. This study shows that CQ remained largely efficacious to treat P. vivax in Gia Lai; i.e., recurrences occurred late (>day 28) and in the presence of low blood CQ concentrations. However, the combination of both WGS and gene expression analysis (pvcrt) data with clinical data (PCT) allowed us to identify potential emergence of low-grade CQR, which should be closely monitored. (This study has been registered at ClinicalTrials.gov under identifier NCT02610686.).

Malaria in the 'Omics Era'

Genomics has revolutionised the study of the biology of parasitic diseases. The first Eukaryotic parasite to have its genome sequenced was the malaria parasite Plasmodium falciparum. Since then, Plasmodium genomics has continued to lead the way in the study of the genome biology of parasites, both in breadth-the number of Plasmodium species' genomes sequenced-and in depth-massive-scale genome re-sequencing of several key species. Here, we review some of the insights into the biology, evolution and population genetics of Plasmodium gained from genome sequencing, and look at potential new avenues in the future genome-scale study of its biology.

Distinctive genetic structure and selection patterns in Plasmodium vivax from South Asia and East Africa

Despite the high burden of Plasmodium vivax malaria in South Asian countries, the genetic diversity of circulating parasite populations is not well described. Determinants of antimalarial drug susceptibility for P. vivax in the region have not been characterised. Our genomic analysis of global P. vivax (n = 558) establishes South Asian isolates (n = 92) as a distinct subpopulation, which shares ancestry with some East African and South East Asian parasites. Signals of positive selection are linked to drug resistance-associated loci including pvkelch10, pvmrp1, pvdhfr and pvdhps, and two loci linked to P. vivax invasion of reticulocytes, pvrbp1a and pvrbp1b. Significant identity-by-descent was found in extended chromosome regions common to P. vivax from India and Ethiopia, including the pvdbp gene associated with Duffy blood group binding. Our investigation provides new understanding of global P. vivax population structure and genomic diversity, and genetic evidence of recent directional selection in this important human pathogen.

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.

Population genomic evidence of Plasmodium vivax Southeast Asian origin

Plasmodium vivax is the most common and widespread human malaria parasite. It was recently proposed that P. vivax originates from sub-Saharan Africa based on the circulation of its closest genetic relatives (P. vivax-like) among African great apes. However, the limited number of genetic markers and samples investigated questions the robustness of this hypothesis. Here, we extensively characterized the genomic variations of 447 human P. vivax strains and 19 ape P. vivax-like strains collected worldwide. Phylogenetic relationships between human and ape Plasmodium strains revealed that P. vivax is a sister clade of P. vivax-like, not included within the radiation of P. vivax-like By investigating various aspects of P. vivax genetic variation, we identified several notable geographical patterns in summary statistics in function of the increasing geographic distance from Southeast Asia, suggesting that P. vivax may have derived from a single area in Asia through serial founder effects.

Reconstruction of Plasmodium vivax outbreaks in a low malaria endemic setting utilizing conventional restriction fragment length polymorphism

Suriname is on track to eliminate local malaria transmission. P. vivax malaria reemerged in March and September 2019 in the Amerindian village Palumeu, free of malaria for two years and concurrently, a case was reported in another village Alalaparoe. The outbreaks were contained through targeted interventions including Mass Drug Administration (MDA). Molecular outbreak analysis was performed on 23 dried blood spots (DBS) using combined polymerase chain reaction/restriction fragment length polymorphism (PCR-RFLP) with Pvmsp-1 F2 and Pvmsp-3α as polymorphic marker genes. Independent controls substantiated the discriminating capacities of the utilized PCR-RFLP method. All isolates from the first and second Palumeu outbreak shared a distinctive haplotype presuming single clonal lineage. An imported case probably triggered the first outbreak, while a delayed episode, prompted by withdrawal of drug pressure at the end of the prophylactic MDA, was suggested as source of the second outbreak. A diverging variant was demonstrated in Alalaparoe, implicating an infection from a different source. PCR-RFLP proved to be a useful molecular tool for P. vivax outbreak management in low endemic malaria settings.

Multi-locus genotyping reveals established endemicity of a geographically distinct Plasmodium vivax population in Mauritania, West Africa

BACKGROUND

Plasmodium vivax has been recently discovered as a significant cause of malaria in Mauritania, although very rare elsewhere in West Africa. It has not been known if this is a recently introduced or locally remnant parasite population, nor whether the genetic structure reflects epidemic or endemic transmission.

METHODOLOGY/PRINCIPAL FINDINGS

To investigate the P. vivax population genetic structure in Mauritania and compare with populations previously analysed elsewhere, multi-locus genotyping was undertaken on 100 clinical isolates, using a genome-wide panel of 38 single nucleotide polymorphisms (SNPs), plus seven SNPs in drug resistance genes. The Mauritanian P. vivax population is shown to be genetically diverse and divergent from populations elsewhere, indicated consistently by genetic distance matrix analysis, principal components analyses, and fixation indices. Only one isolate had a genotype clearly indicating recent importation, from a southeast Asian source. There was no linkage disequilibrium in the local parasite population, and only a small number of infections appeared to be closely genetically related, indicating that there is ongoing genetic recombination consistent with endemic transmission. The P. vivax diversity in a remote mining town was similar to that in the capital Nouakchott, with no indication of local substructure or of epidemic population structure. Drug resistance alleles were virtually absent in Mauritania, in contrast with P. vivax in other areas of the world.

CONCLUSIONS/SIGNIFICANCE

The molecular epidemiology indicates that there is long-standing endemic transmission that will be very challenging to eliminate. The virtual absence of drug resistance alleles suggests that most infections have been untreated, and that this endemic infection has been more neglected in comparison to P. vivax elsewhere.

Whole genome sequencing of Plasmodium vivax isolates reveals frequent sequence and structural polymorphisms in erythrocyte binding genes

Plasmodium vivax malaria is much less common in Africa than the rest of the world because the parasite relies primarily on the Duffy antigen/chemokine receptor (DARC) to invade human erythrocytes, and the majority of Africans are Duffy negative. Recently, there has been a dramatic increase in the reporting of P. vivax cases in Africa, with a high number of them being in Duffy negative individuals, potentially indicating P. vivax has evolved an alternative invasion mechanism that can overcome Duffy negativity. Here, we analyzed single nucleotide polymorphism (SNP) and copy number variation (CNV) in Whole Genome Sequence (WGS) data from 44 P. vivax samples isolated from symptomatic malaria patients in southwestern Ethiopia, where both Duffy positive and Duffy negative individuals are found. A total of 123,711 SNPs were detected, of which 22.7% were nonsynonymous and 77.3% were synonymous mutations. The largest number of SNPs were detected on chromosomes 9 (24,007 SNPs; 19.4% of total) and 10 (16,852 SNPs, 13.6% of total). There were particularly high levels of polymorphism in erythrocyte binding gene candidates including merozoite surface protein 1 (MSP1) and merozoite surface protein 3 (MSP3.5, MSP3.85 and MSP3.9). Two genes, MAEBL and MSP3.8 related to immunogenicity and erythrocyte binding function were detected with significant signals of positive selection. Variation in gene copy number was also concentrated in genes involved in host-parasite interactions, including the expansion of the Duffy binding protein gene (PvDBP) on chromosome 6 and MSP3.11 on chromosome 10. Based on the phylogeny constructed from the whole genome sequences, the expansion of these genes was an independent process among the P. vivax lineages in Ethiopia. We further inferred transmission patterns of P. vivax infections among study sites and showed various levels of gene flow at a small geographical scale. The genomic features of P. vivax provided baseline data for future comparison with those in Duffy-negative individuals and allowed us to develop a panel of informative Single Nucleotide Polymorphic markers diagnostic at a micro-geographical scale.

Population genomics reveals the expansion of highly inbred Plasmodium vivax lineages in the main malaria hotspot of Brazil

BACKGROUND

Plasmodium vivax is a neglected human malaria parasite that causes significant morbidity in the Americas, the Middle East, Asia, and the Western Pacific. Population genomic approaches remain little explored to map local and regional transmission pathways of P. vivax across the main endemic sites in the Americas, where great progress has been made towards malaria elimination over the past decades.

METHODOLOGY/PRINCIPAL FINDINGS

We analyze 38 patient-derived P. vivax genome sequences from Mâncio Lima (ML)-the Amazonian malaria hotspot next to the Brazil-Peru border-and 24 sequences from two other sites in Acre State, Brazil, a country that contributes 23% of malaria cases in the Americas. We show that the P. vivax population of ML is genetically diverse (π = 4.7 × 10-4), with a high polymorphism particularly in genes encoding proteins putatively involved in red blood cell invasion. Paradoxically, however, parasites display strong genome-wide linkage disequilibrium, being fragmented into discrete lineages that are remarkably stable across time and space, with only occasional recombination between them. Using identity-by-descent approaches, we identified a large cluster of closely related sequences that comprises 16 of 38 genomes sampled in ML over 26 months. Importantly, we found significant ancestry sharing between parasites at a large geographic distance, consistent with substantial gene flow between regional P. vivax populations.

CONCLUSIONS/SIGNIFICANCE

We have characterized the sustained expansion of highly inbred P. vivax lineages in a malaria hotspot that can seed regional transmission. Potential source populations in hotspots represent a priority target for malaria elimination in the Amazon.

Hotspots in Plasmodium and RBC Receptor-Ligand Interactions: Key Pieces for Inhibiting Malarial Parasite Invasion

Protein-protein interactions (IPP) play an essential role in practically all biological processes, including those related to microorganism invasion of their host cells. It has been found that a broad repertoire of receptor-ligand interactions takes place in the binding interphase with host cells in malaria, these being vital interactions for successful parasite invasion. Several trials have been conducted for elucidating the molecular interface of interactions between some Plasmodium falciparum and Plasmodium vivax antigens with receptors on erythrocytes and/or reticulocytes. Structural information concerning these complexes is available; however, deeper analysis is required for correlating structural, functional (binding, invasion, and inhibition), and polymorphism data for elucidating new interaction hotspots to which malaria control methods can be directed. This review describes and discusses recent structural and functional details regarding three relevant interactions during erythrocyte invasion: Duffy-binding protein 1 (DBP1)–Duffy antigen receptor for chemokines (DARC); reticulocyte-binding protein homolog 5 (PfRh5)-basigin, and erythrocyte binding antigen 175 (EBA175)-glycophorin A (GPA).

Genetic diversity and neutral selection in Plasmodium vivax erythrocyte binding protein correlates with patient antigenicity

Plasmodium vivax is the most widespread and difficult to treat cause of human malaria. The development of vaccines against the blood stages of P. vivax remains a key objective for the control and elimination of vivax malaria. Erythrocyte binding-like (EBL) protein family members such as Duffy binding protein (PvDBP) are of critical importance to erythrocyte invasion and have been the major target for vivax malaria vaccine development. In this study, we focus on another member of EBL protein family, P. vivax erythrocyte binding protein (PvEBP). PvEBP was first identified in Cambodian (C127) field isolates and has subsequently been showed its preferences for binding reticulocytes which is directly inhibited by antibodies. We analysed PvEBP sequence from 316 vivax clinical isolates from eight countries including China (n = 4), Ethiopia (n = 24), Malaysia (n = 53), Myanmar (n = 10), Papua New Guinea (n = 16), Republic of Korea (n = 10), Thailand (n = 174), and Vietnam (n = 25). PvEBP gene exhibited four different phenotypic clusters based on the insertion/deletion (indels) variation. PvEBP-RII (179-479 aa.) showed highest polymorphism similar to other EBL family proteins in various Plasmodium species. Whereas even though PvEBP-RIII-V (480-690 aa.) was the most conserved domain, that showed strong neutral selection pressure for gene purifying with significant population expansion. Antigenicity of both of PvEBP-RII (16.1%) and PvEBP-RIII-V (21.5%) domains were comparatively lower than other P. vivax antigen which expected antigens associated with merozoite invasion. Total IgG recognition level of PvEBP-RII was stronger than PvEBP-RIII-V domain, whereas total IgG inducing level was stronger in PvEBP-RIII-V domain. These results suggest that PvEBP-RII is mainly recognized by natural IgG for innate protection, whereas PvEBP-RIII-V stimulates IgG production activity by B-cell for acquired immunity. Overall, the low antigenicity of both regions in patients with vivax malaria likely reflects genetic polymorphism for strong positive selection in PvEBP-RII and purifying selection in PvEBP-RIII-V domain. These observations pose challenging questions to the selection of EBP and point out the importance of immune pressure and polymorphism required for inclusion of PvEBP as a vaccine candidate.

A molecular barcode to inform the geographical origin and transmission dynamics of Plasmodium vivax malaria

Although Plasmodium vivax parasites are the predominant cause of malaria outside of sub-Saharan Africa, they not always prioritised by elimination programmes. P. vivax is resilient and poses challenges through its ability to re-emerge from dormancy in the human liver. With observed growing drug-resistance and the increasing reports of life-threatening infections, new tools to inform elimination efforts are needed. In order to halt transmission, we need to better understand the dynamics of transmission, the movement of parasites, and the reservoirs of infection in order to design targeted interventions. The use of molecular genetics and epidemiology for tracking and studying malaria parasite populations has been applied successfully in P. falciparum species and here we sought to develop a molecular genetic tool for P. vivax. By assembling the largest set of P. vivax whole genome sequences (n = 433) spanning 17 countries, and applying a machine learning approach, we created a 71 SNP barcode with high predictive ability to identify geographic origin (91.4%). Further, due to the inclusion of markers for within population variability, the barcode may also distinguish local transmission networks. By using P. vivax data from a low-transmission setting in Malaysia, we demonstrate the potential ability to infer outbreak events. By characterising the barcoding SNP genotypes in P. vivax DNA sourced from UK travellers (n = 132) to ten malaria endemic countries predominantly not used in the barcode construction, we correctly predicted the geographic region of infection origin. Overall, the 71 SNP barcode outperforms previously published genotyping methods and when rolled-out within new portable platforms, is likely to be an invaluable tool for informing targeted interventions towards elimination of this resilient human malaria.

The enigmatic mechanisms by which Plasmodium vivax infects Duffy-negative individuals

The absence of the Duffy protein at the surface of erythrocytes was considered for decades to confer full protection against Plasmodium vivax as this blood group is the receptor for the key parasite ligand P. vivax Duffy binding protein (PvDBP). However, it is now clear that the parasite is able to break through this protection and induce clinical malaria in Duffy-negative people, although the underlying mechanisms are still not understood. Here, we briefly review the evidence of Duffy-negative infections by P. vivax and summarize the current hypothesis at the basis of this invasion process. We discuss those in the perspective of malaria-elimination challenges, notably in African countries.

Plasmodium vivax chloroquine resistance links to pvcrt transcription in a genetic cross

Mainstay treatment for Plasmodium vivax malaria has long relied on chloroquine (CQ) against blood-stage parasites plus primaquine against dormant liver-stage forms (hypnozoites), however drug resistance confronts this regimen and threatens malaria control programs. Understanding the basis of P. vivax chloroquine resistance (CQR) will inform drug discovery and malaria control. Here we investigate the genetics of P. vivax CQR by a cross of parasites differing in drug response. Gametocytogenesis, mosquito infection, and progeny production are performed with mixed parasite populations in nonhuman primates, as methods for P. vivax cloning and in vitro cultivation remain unavailable. Linkage mapping of progeny surviving >15 mg/kg CQ identifies a 76 kb region in chromosome 1 including pvcrt, an ortholog of the Plasmodium falciparum CQR transporter gene. Transcriptional analysis supports upregulated pvcrt expression as a mechanism of CQR.

Frequent expansion of Plasmodium vivax Duffy Binding Protein in Ethiopia and its epidemiological significance

Plasmodium vivax invasion of human erythrocytes depends on the Duffy Binding Protein (PvDBP) which interacts with the Duffy antigen. PvDBP copy number has been recently shown to vary between P. vivax isolates in Sub-Saharan Africa. However, the extent of PvDBP copy number variation, the type of PvDBP multiplications, as well as its significance across broad samples are still unclear. We determined the prevalence and type of PvDBP duplications, as well as PvDBP copy number variation among 178 Ethiopian P. vivax isolates using a PCR-based diagnostic method, a novel quantitative real-time PCR assay and whole genome sequencing. For the 145 symptomatic samples, PvDBP duplications were detected in 95 isolates, of which 81 had the Cambodian and 14 Malagasy-type PvDBP duplications. PvDBP varied from 1 to >4 copies. Isolates with multiple PvDBP copies were found to be higher in symptomatic than asymptomatic infections. For the 33 asymptomatic samples, PvDBP was detected with two copies in two of the isolates, and both were the Cambodian-type PvDBP duplication. PvDBP copy number in Duffy-negative heterozygotes was not significantly different from that in Duffy-positives, providing no support for the hypothesis that increased copy number is a specific association with Duffy-negativity, although the number of Duffy-negatives was small and further sampling is required to test this association thoroughly.

Plasmodium vivax invasion of human erythrocytes relies on interaction between the Duffy antigen and P. vivax Duffy Binding Protein (PvDBP). Whole genome sequences from P. vivax field isolates in Madagascar identified a duplication of the PvDBP gene and PvDBP duplication has also been detected in non-African P. vivax-endemic countries. Two types of PvDBP duplications have been reported, termed Cambodian and Malagasy-type duplications. Our study used a combination of PCR-based diagnostic method, a novel quantitative real-time PCR assay, and whole genome sequencing to determine the prevalence and type of PvDBP duplications, as well as PvDBP copy number on a broad number of P. vivax samples in Ethiopia. We found that over 65% of P. vivax isolated from the symptomatic infections were detected with PvDBP duplications and PvDBP varied from 1 to >4 copies. The majority of PvDBP duplications belongs to the Cambodian-type while the Malagasy-type duplications was also detected. For the asymptomatic infections, despite a small sample size, the majority of P. vivax were detected with a single-copy based on both PCR and qPCR assays. There was no significant difference in PvDBP copy number between Duffy-null heterozygote and Duffy-positive homozygote/heterozygote. Further investigation is needed with expanded Duffy-null homozygotes to examine the functional significance of PvDBP expansion.