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

Sterile protection against P. vivax malaria by repeated blood stage infection in the Aotus monkey model

The malaria parasite Plasmodium vivax remains a major global public health challenge, and no vaccine is approved for use in humans. Here, we assessed whether P. vivax strain-transcendent immunity can be achieved by repeated infection in Aotus monkeys. Sterile immunity was achieved after two homologous infections, whereas subsequent heterologous challenge provided only partial protection. IgG levels based on P. vivax lysate ELISA and protein microarray increased with repeated infections and correlated with the level of homologous protection. Parasite transcriptional profiles provided no evidence of major antigenic switching upon homologous or heterologous challenge. However, we observed significant sequence diversity and transcriptional differences in the P. vivax core gene repertoire between the two strains used in the study, suggesting that partial protection upon heterologous challenge is due to molecular differences between strains rather than immune evasion by antigenic switching. Our study demonstrates that sterile immunity against P. vivax can be achieved by repeated homologous blood stage infection in Aotus monkeys, thus providing a benchmark to test the efficacy of candidate blood stage P. vivax malaria vaccines.

Comparative transcriptomics reveal differential gene expression among Plasmodium vivax geographical isolates and implications on erythrocyte invasion mechanisms

The documentation of Plasmodium vivax malaria across Africa especially in regions where Duffy negatives are dominant suggests possibly alternative erythrocyte invasion mechanisms. While the transcriptomes of the Southeast Asian and South American P. vivax are well documented, the gene expression profile of P. vivax in Africa is unclear. In this study, we examined the expression of 4,404 gene transcripts belong to 12 functional groups and 43 erythrocyte binding gene candidates in Ethiopian isolates and compared them with the Cambodian and Brazilian P. vivax transcriptomes. Overall, there were 10-26% differences in the gene expression profile amongst geographical isolates, with the Ethiopian and Cambodian P. vivax being most similar. Majority of the gene transcripts involved in protein transportation, housekeeping, and host interaction were highly transcribed in the Ethiopian isolates. Members of the reticulocyte binding protein PvRBP2a and PvRBP3 expressed six-fold higher than Duffy binding protein PvDBP1 and 60-fold higher than PvEBP/DBP2 in the Ethiopian isolates. Other genes including PvMSP3.8, PvMSP3.9, PvTRAG2, PvTRAG14, and PvTRAG22 also showed relatively high expression. Differential expression patterns were observed among geographical isolates, e.g., PvDBP1 and PvEBP/DBP2 were highly expressed in the Cambodian but not the Brazilian and Ethiopian isolates, whereas PvRBP2a and PvRBP2b showed higher expression in the Ethiopian and Cambodian than the Brazilian isolates. Compared to Pvs25, gametocyte genes including PvAP2-G, PvGAP (female gametocytes), and Pvs47 (male gametocytes) were highly expressed across geographical samples.

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.

Genetic diversity of merozoite surface protein-5 (MSP-5) of Plasmodium vivax isolates from Malaria patients in Iran

Malaria has not yet been eradicated in Iran, and Plasmodium vivax (P. vivax) is the main cause of malaria in the country. This study aimed to investigate and analyze the amount of genetic diversity of Plasmodium vivax merozoite surface protein-5 (PvMSP-5) exon 1 gene in the southeast of Iran.Thirty-five patients with clinical symptoms of P. vivax malaria participated. The exon 1 of PvMSP-5 was amplified by PCR, and the PCR product of all isolates was sequenced, and genetic polymorphisms were determined using various genetic software.The analysis showed that studied isolates are different from one another in the DnaSP software version. Out of the 612 sites, 477 were monomorphic and 135 were segregated. The total number of mutations was 143. The singleton variable and the parsimony informative sites were 23 and 112, respectively. There were 17 specific haplotypes with haplotype diversity equal to 0.943. Nucleotide diversity was equal to 0.06766 in the isolates. The ratio of nonsynonymous (0.06446) to synonymous (0.07909) mutations was 0.815020. Tajima's D, which expressed coding, and non-coding regions, was 0.72403, which was not deemed significant (P > 0.10).The analysis of intrapopulation diversity revealed nucleotide and haplotype diversity in the msp-5 gene of Iranian P. vivax isolates. In addition to balancing or purifying selection, intragenic recombination also contributed to the variation observed in exon 1 of PvMSP-5, according to the findings.

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.

A new Plasmodium vivax reference genome for South American isolates

BACKGROUND

Plasmodium vivax is the second most important cause of human malaria worldwide, and accounts for the majority of malaria cases in South America. A high-quality reference genome exists for Papua Indonesia (PvP01) and Thailand (PvW1), but is lacking for South America. A reference genome specifically for South America would be beneficial though, as P. vivax is a genetically diverse parasite with geographical clustering.

RESULTS

This study presents a new high-quality assembly of a South American P. vivax isolate, referred to as PvPAM (P. vivax Peruvian AMazon). The genome was obtained from a low input patient sample from the Peruvian Amazon and sequenced using PacBio technology, resulting in a highly complete assembly with 6497 functional genes. Telomeric ends were present in 17 out of 28 chromosomal ends, and additional (sub)telomeric regions are present in 12 unassigned contigs. A comparison of multigene families between PvPAM and the PvP01 genome revealed remarkable variation in vir genes, and the presence of merozoite surface proteins (MSP) 3.6 and 3.7. Three dhfr and dhps drug resistance associated mutations are present in PvPAM, similar to those found in other Peruvian isolates. Mapping of publicly available South American whole genome sequencing (WGS) data to PvPAM resulted in significantly fewer variants and truncated reads compared to the use of PvP01 or PvW1 as reference genomes. To minimize the number of core genome variants in non-South American samples, PvW1 is most suited for Southeast Asian isolates, both PvPAM and PvW1 are suited for South Asian isolates, and PvPAM is recommended for African isolates. Interestingly, non-South American samples still contained the least subtelomeric variants when mapped to PvPAM, indicating high quality of the PvPAM subtelomeric regions.

CONCLUSIONS

Our findings show that the PvPAM reference genome more accurately represents South American P. vivax isolates in comparison to PvP01 and PvW1. In addition, PvPAM has a high level of completeness, and contains a similar number of annotated genes as PvP01 or PvW1. The PvPAM genome therefore will be a valuable resource to improve future genomic analyses on P. vivax isolates from the South American continent.

Comparative transcriptomics reveal differential gene expression in Plasmodium vivax geographical isolates and implications on erythrocyte invasion mechanisms

Plasmodium vivax uses Duffy binding protein (PvDBP1) to bind to the Duffy Antigen-Chemokine Receptor (DARC) to invade human erythrocytes. Individuals who lack DARC expression (Duffy-negative) are thought to be resistance to P. vivax. In recent years, P. vivax malaria is becoming more prevalent in Africa with a portion of these cases detected in Duffy-negatives. Apart from DBP1, members of the reticulocyte binding protein (RBP) and tryptophan-rich antigen (TRAg) families may also play a role in erythrocyte invasion. While the transcriptomes of the Southeast Asian and South American P. vivax are well documented, the gene expression profile of P. vivax in Africa and more specifically the expression level of several erythrocyte binding gene candidates as compared to DBP1 are largely unknown. This paper characterized the first P. vivax transcriptome in Africa and compared with those from the Southeast Asian and South American isolates. The expression of 4,404 gene transcripts belong to 12 functional groups including 43 specific erythrocyte binding gene candidates were examined. Overall, there were 10-26% differences in the gene expression profile amongst the geographical isolates, with the Ethiopian and Cambodian P. vivax being most similar. Majority of the gene transcripts involved in protein transportation, housekeeping, and host interaction were highly transcribed in the Ethiopian P. vivax. Erythrocyte binding genes including PvRBP2a and PvRBP3 expressed six-fold higher than PvDBP1and 60-fold higher than PvEBP/DBP2. Other genes including PvRBP1a, PvMSP3.8, PvMSP3.9, PvTRAG2, PvTRAG14, and PvTRAG22 also showed relatively high expression. Differential expression was observed among geographical isolates, e.g., PvDBP1 and PvEBP/DBP2 were highly expressed in the Cambodian but not the Brazilian and Ethiopian isolates, whereas PvRBP2a and PvRBP2b showed higher expression in the Ethiopian and Cambodian than the Brazilian isolates. Compared to Pvs25, the standard biomarker for detecting female gametocytes, PvAP2-G (PVP01_1440800), GAP (PVP01_1403000), and Pvs47 (PVP01_1208000) were highly expressed across geographical samples. These findings provide an important baseline for future comparisons of P. vivax transcriptomes from Duffy-negative infections and highlight potential biomarkers for improved gametocyte detection.

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.

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.

Systems biology of malaria explored with nonhuman primates

"The Primate Malarias" book has been a uniquely important resource for multiple generations of scientists, since its debut in 1971, and remains pertinent to the present day. Indeed, nonhuman primates (NHPs) have been instrumental for major breakthroughs in basic and pre-clinical research on malaria for over 50 years. Research involving NHPs have provided critical insights and data that have been essential for malaria research on many parasite species, drugs, vaccines, pathogenesis, and transmission, leading to improved clinical care and advancing research goals for malaria control, elimination, and eradication. Whilst most malaria scientists over the decades have been studying Plasmodium falciparum, with NHP infections, in clinical studies with humans, or using in vitro culture or rodent model systems, others have been dedicated to advancing research on Plasmodium vivax, as well as on phylogenetically related simian species, including Plasmodium cynomolgi, Plasmodium coatneyi, and Plasmodium knowlesi. In-depth study of these four phylogenetically related species over the years has spawned the design of NHP longitudinal infection strategies for gathering information about ongoing infections, which can be related to human infections. These Plasmodium-NHP infection model systems are reviewed here, with emphasis on modern systems biological approaches to studying longitudinal infections, pathogenesis, immunity, and vaccines. Recent discoveries capitalizing on NHP longitudinal infections include an advanced understanding of chronic infections, relapses, anaemia, and immune memory. With quickly emerging new technological advances, more in-depth research and mechanistic discoveries can be anticipated on these and additional critical topics, including hypnozoite biology, antigenic variation, gametocyte transmission, bone marrow dysfunction, and loss of uninfected RBCs. New strategies and insights published by the Malaria Host-Pathogen Interaction Center (MaHPIC) are recapped here along with a vision that stresses the importance of educating future experts well trained in utilizing NHP infection model systems for the pursuit of innovative, effective interventions against malaria.

Effectiveness of pyronaridine-artesunate against Plasmodium malariae, Plasmodium ovale spp, and mixed-Plasmodium infections: a post-hoc analysis of the CANTAM-Pyramax trial

Background

Methods

Findings

Interpretation

Funding

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.

GIP: an open-source computational pipeline for mapping genomic instability from protists to cancer cells

Genome instability has been recognized as a key driver for microbial and cancer adaptation and thus plays a central role in many diseases. Genome instability encompasses different types of genomic alterations, yet most available genome analysis software are limited to just one type of mutation. To overcome this limitation and better understand the role of genetic changes in enhancing pathogenicity we established GIP, a novel, powerful bioinformatic pipeline for comparative genome analysis. Here, we show its application to whole genome sequencing datasets of Leishmania, Plasmodium, Candida and cancer. Applying GIP on available data sets validated our pipeline and demonstrated the power of our tool to drive biological discovery. Applied to Plasmodium vivax genomes, our pipeline uncovered the convergent amplification of erythrocyte binding proteins and identified a nullisomic strain. Re-analyzing genomes of drug adapted Candida albicans strains revealed correlated copy number variations of functionally related genes, strongly supporting a mechanism of epistatic adaptation through interacting gene-dosage changes. Our results illustrate how GIP can be used for the identification of aneuploidy, gene copy number variations, changes in nucleic acid sequences, and chromosomal rearrangements. Altogether, GIP can shed light on the genetic bases of cell adaptation and drive disease biomarker discovery.

Plasmodium vivax From Duffy-Negative and Duffy-Positive Individuals Share Similar Gene Pools in East Africa

Plasmodium vivax malaria was thought to be rare in Africa, but an increasing number of P. vivax cases reported across Africa and in Duffy-negative individuals challenges this dogma. The genetic characteristics of P. vivax in Duffy-negative infections, the transmission of P. vivax in East Africa, and the impact of environments on transmission remain largely unknown. This study examined genetic and transmission features of P. vivax from 107 Duffy-negative and 305 Duffy-positive individuals in Ethiopia and Sudan. No clear genetic differentiation was found in P. vivax between the 2 Duffy groups, indicating between-host transmission. P. vivax from Ethiopia and Sudan showed similar genetic clusters, except samples from Khartoum, possibly due to distance and road density that inhibited parasite gene flow. This study is the first to show that P. vivax can transmit to and from Duffy-negative individuals and provides critical insights into the spread of P. vivax in sub-Saharan Africa.

Basic Research of Plasmodium vivax Biology Enabling Its Management as a Clinical and Public Health Problem

The emerging understanding of Plasmodium vivax as an infection seated in extravascular spaces of its human host carries fundamentally important implications for its management as a complex clinical and public health problem. This progress begins to reverse decades of neglected research borne of the false dogma of P. vivax as an intrinsically benign and inconsequential parasite. This Review provides real world context for the on-going laboratory explorations of the molecular and cellular events in the life of this parasite. Chemotherapies against the latent reservoir impose extraordinarily complex and difficult problems of science and medicine, but great strides in studies of the biology of hepatic P. vivax promise solutions. Fundamental assumptions regarding the interpretation of parasitaemia in epidemiology, clinical medicine, and public health are being revisited and reassessed in light of new studies of P. vivax cellular/molecular biology and pathogenesis. By examining these long overlooked complexities of P. vivax malaria, we open multiple new avenues to vaccination, chemoprevention, countermeasures against transmission, epidemiology, diagnosis, chemotherapy, and clinical management. This Review expresses how clarity of vision of biology and pathogenesis may rationally and radically transform the multiple means by which we may combat this insidiously harmful infection.

Alternative Invasion Mechanisms and Host Immune Response to Plasmodium vivax Malaria: Trends and Future Directions

Plasmodium vivax malaria is a neglected tropical disease, despite being more geographically widespread than any other form of malaria. The documentation of P. vivax infections in different parts of Africa where Duffy-negative individuals are predominant suggested that there are alternative pathways for P. vivax to invade human erythrocytes. Duffy-negative individuals may be just as fit as Duffy-positive individuals and are no longer resistant to P.vivax malaria. In this review, we describe the complexity of P. vivax malaria, characterize pathogenesis and candidate invasion genes of P. vivax, and host immune responses to P. vivax infections. We provide a comprehensive review on parasite ligands in several Plasmodium species that further justify candidate genes in P. vivax. We also summarize previous genomic and transcriptomic studies related to the identification of ligand and receptor proteins in P. vivax erythrocyte invasion. Finally, we identify topics that remain unclear and propose future studies that will greatly contribute to our knowledge of P. vivax.

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.