Genome-wide diversity and differentiation in New World populations of the human malaria parasite Plasmodium vivax

Plasmodium vivax populations in the western Greater Mekong Subregion evaluated using a genetic barcode

An improved understanding of the Plasmodium vivax populations in the Great Mekong Subregion (GMS) is needed to monitor the progress of malaria elimination. This study aimed to use a P. vivax single nucleotide polymorphism (SNP) barcode to evaluate the population dynamics and explore the gene flow among P. vivax parasite populations in the western GMS (China, Myanmar and Thailand). A total of 315 P. vivax patient samples collected in 2011 and 2018 from four regions of the western GMS were genotyped for 42 SNPs using the high-throughput MassARRAY SNP genotyping technology. Population genetic analysis was conducted to estimate the genetic diversity, effective population size, and population structure among the P. vivax populations. Overall, 291 samples were successfully genotyped at 39 SNPs. A significant difference was observed in the proportion of polyclonal infections among the five P. vivax populations (P = 0.0012, Pearson Chi-square test, χ2 = 18.1), with western Myanmar having the highest proportion (96.2%, 50/52) in 2018. Likewise, the average complexity of infection was also highest in western Myanmar (1.31) and lowest in northeast Myanmar (1.01) in 2018. The older samples from western China in 2011 had the highest pairwise nucleotide diversity (π, 0.388 ± 0.046), expected heterozygosity (He, 0.363 ± 0.02), and the largest effective population size. In comparison, in the neighboring northeast Myanmar, the more recent samples in 2018 showed the lowest values (π, 0.224 ± 0.036; He, 0.220 ± 0.026). Furthermore, the 2018 northeast Myanmar parasites showed high and moderate genetic differentiation from other populations with FST values of 0.162-0.252, whereas genetic differentiation among other populations was relatively low (FST ≤ 0.059). Principal component analysis, phylogeny, and STRUCTURE analysis showed that the P. vivax population in northeast Myanmar in 2018 substantially diverged from other populations. Although the 42 SNP barcode is a valuable tool for tracking parasite origins of worldwide parasite populations, a more extended barcode with additional SNPs is needed to distinguish the more related parasite populations in the western GMS.

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.

New reference genomes to distinguish the sympatric malaria parasites, Plasmodium ovale curtisi and Plasmodium ovale wallikeri

Despite Plasmodium ovale curtisi (Poc) and wallikeri (Pow) being important human-infecting malaria parasites that are widespread across Africa and Asia, little is known about their genome diversity. Morphologically identical, Poc and Pow are indistinguishable and commonly misidentified. Recent rises in the incidence of Poc/Pow infections have renewed efforts to address fundamental knowledge gaps in their biology, and to develop diagnostic tools to understand their epidemiological dynamics and malaria burden. A major roadblock has been the incompleteness of available reference assemblies (PocGH01, PowCR01; ~ 33.5 Mbp). Here, we applied multiple sequencing platforms and advanced bioinformatics tools to generate new reference genomes, Poc221 (South Sudan; 36.0 Mbp) and Pow222 (Nigeria; 34.3 Mbp), with improved nuclear genome contiguity (> 4.2 Mbp), annotation and completeness (> 99% Plasmodium spp., single copy orthologs). Subsequent sequencing of 6 Poc and 15 Pow isolates from Africa revealed a total of 22,517 and 43,855 high-quality core genome SNPs, respectively. Genome-wide levels of nucleotide diversity were determined to be 2.98 × 10-4 (Poc) and 3.43 × 10-4 (Pow), comparable to estimates for other Plasmodium species. Overall, the new reference genomes provide a robust foundation for dissecting the biology of Poc/Pow, their population structure and evolution, and will contribute to uncovering the recombination barrier separating these species.

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.

The structure of a Plasmodium vivax Tryptophan Rich Antigen domain suggests a lipid binding function for a pan-Plasmodium multi-gene family

Tryptophan Rich Antigens (TRAgs) are encoded by a multi-gene family found in all Plasmodium species, but are significantly expanded in P. vivax and closely related parasites. We show that multiple P. vivax TRAgs are expressed on the merozoite surface and that one, PVP01_0000100 binds red blood cells with a strong preference for reticulocytes. Using X-ray crystallography, we solved the structure of the PVP01_0000100 C-terminal tryptophan rich domain, which defines the TRAg family, revealing a three-helical bundle that is conserved across Plasmodium and has structural homology with lipid-binding BAR domains involved in membrane remodelling. Biochemical assays confirm that the PVP01_0000100 C-terminal domain has lipid binding activity with preference for sulfatide, a glycosphingolipid present in the outer leaflet of plasma membranes. Deletion of the putative orthologue in P. knowlesi, PKNH_1300500, impacts invasion in reticulocytes, suggesting a role during this essential process. Together, this work defines an emerging molecular function for the Plasmodium TRAg family.

The need to increase antimicrobial resistance surveillance among forcibly displaced persons (FDPs)

Antimicrobial resistance (AMR) poses a significant threat to human health as 4.95 million deaths were associated with bacterial AMR in 2019 and is projected to reach 10 million by 2050. To mitigate AMR, surveillance is an essential tool for determining the burden of AMR and providing the necessary information for its control. However, the global AMR surveillance is inadequate and particularly limited among forcibly displaced persons (FDPs) despite having higher risks of harboring these pathogens. Predisposing factors among this group include poor living conditions, limited access to treatment and diagnostic tests, and inadequate trained health professionals in refugee camps. Strengthening AMR surveillance among FDPs would address the identified gaps and facilitate formulation and implementation of evidence-based policies on AMR control and prevention response. This article provides information on the growing population of FDPs, factors contributing to the AMR burden and AMR surveillance gaps in FDPs and highlighted recommendations for control.

Development of a Plasmodium vivax biobank for functional ex vivo assays

BACKGROUND

Plasmodium vivax is the second most prevalent cause of malaria yet remains challenging to study due to the lack of a continuous in vitro culture system, highlighting the need to establish a biobank of clinical isolates with multiple freezes per sample for use in functional assays. Different methods for cryopreserving parasite isolates were compared and subsequently the most promising one was validated. Enrichment of early- and late-stage parasites and parasite maturation were quantified to facilitate assay planning.

METHODS

In order to compare cryopreservation protocols, nine clinical P. vivax isolates were frozen with four glycerolyte-based mixtures. Parasite recovery post thaw, post KCl-Percoll enrichment and in short-term in vitro culture was measured via slide microscopy. Enrichment of late-stage parasites by magnetic activated cell sorting (MACS) was measured. Short and long-term storage of parasites at either - 80 °C or liquid nitrogen were also compared.

RESULTS

Of the four cryopreservation mixtures, one mixture (glycerolyte:serum:RBC at a 2.5:1.5:1 ratio) resulted in improved parasite recovery and statistically significant (P < 0.05) enhancement in parasite survival in short-term in vitro culture. A parasite biobank was subsequently generated using this protocol resulting in a collection of 106 clinical isolates, each with 8 vials. The quality of the biobank was validated by measuring several factors from 47 thaws: the average reduction in parasitaemia post-thaw (25.3%); the average fold enrichment post KCl-Percoll (6.65-fold); and the average percent recovery of parasites (22.0%, measured from 30 isolates). During short-term in vitro culture, robust maturation of ring stage parasites to later stages (> 20% trophozoites, schizonts and gametocytes) was observed in 60.0% of isolates by 48 h. Enrichment of mature parasite stages via MACS showed good reproducibility, with an average of 30.0% post-MACS parasitaemia and an average of 5.30 × 105 parasites/vial. Finally, the effect of storage temperature was tested, and no large impacts from short-term (7 days) or long-term (7-10 years) storage at - 80 °C on parasite recovery, enrichment or viability was observed.

CONCLUSIONS

Here, an optimized freezing method for P. vivax clinical isolates is demonstrated as a template for the generation and validation of a parasite biobank for use in functional assays.

Development of a Plasmodium vivax biobank for functional ex vivo assays

Background

Plasmodium vivax is the second most prevalent cause of malaria yet remains challenging to study due to the lack of a continuous in vitro culture system, highlighting the need to establish a biobank of clinical isolates with multiple freezes per sample for use in functional assays. Different methods for cryopreserving parasite isolates were compared and subsequently the most promising one was validated. Enrichment of early- and late-stage parasites and parasite maturation were quantified to facilitate assay planning.

Methods

In order to compare cryopreservation protocols, nine clinical P. vivax isolates were frozen with four glycerolyte-based mixtures. Parasite recovery post thaw, post KCl-Percoll enrichment and in short-term in vitro culture was measured via slide microscopy. Enrichment of late-stage parasites by magnetic activated cell sorting (MACS) was measured. Short and long-term storage of parasites at either -80°C or liquid nitrogen were also compared.

Results

Of the four cryopreservation mixtures, one mixture (glycerolyte:serum:RBC at a 2.5:1.5:1 ratio) resulted in improved parasite recovery and statistically significant (P<0.05) enhancement in parasite survival in short-term in vitro culture. A parasite biobank was subsequently generated using this protocol resulting in a collection with 106 clinical isolates, each with 8 vials. The quality of the biobank was validated by measuring several factors from 47 thaws: the average reduction in parasitemia post-thaw (25.3%); the average fold enrichment post KCl-Percoll (6.65-fold); and the average percent recovery of parasites (22.0%, measured from 30 isolates). During short-term in vitro culture, robust maturation of ring stage parasites to later stages (>20% trophozoites, schizonts and gametocytes) was observed in 60.0% of isolates by 48 hours. Enrichment of mature parasite stages via MACS showed good reproducibility, with an average 30.0% post-MACS parasitemia and an average 5.30 × 10 5 parasites/vial. Finally, the effect of storage temperature was tested, and no large impacts from short-term (7 day) or long term (7 - 10 year) storage at -80°C on parasite recovery, enrichment or viability was observed.

Conclusions

Here, an optimized freezing method for P. vivax clinical isolates is demonstrated as a template for the generation and validation of a parasite biobank for use in functional assays.

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).

Drug resistance and population structure of Plasmodium falciparum and Plasmodium vivax in the Peruvian Amazon

Malaria is a major health problem in Peru despite substantial progress achieved by the ongoing malaria elimination program. This study explored the population genetics of 63 Plasmodium falciparum and 170 P. vivax cases collected in the Peruvian Amazon Basin between 2015 and 2019. Microscopy and PCR were used for malaria detection and positive samples were genotyped at neutral and drug resistance-associated regions. The P. falciparum population exhibited a low nucleotide diversity (π = 0.02) whereas the P. vivax population presented a higher genetic diversity (π = 0.34). All P. falciparum samples (n = 63) carried chloroquine (CQ) resistant mutations on Pfcrt. Most P. falciparum samples (53 out of 54) carried sulfadoxine (SD) resistant mutations on Pfdhfr and Pfdhps. No evidence was found of artemisinin resistance mutations on kelch13. Population structure showed that a single cluster accounted for 93.4% of the P. falciparum samples whereas three clusters were found for P. vivax. Our study shows a low genetic diversity for both species with significant differences in genetic sub-structuring. The high prevalence of CQ-resistance mutations could be a result of indirect selection pressures driven by the P. vivax treatment scheme. These results could be useful for public health authorities to safeguard the progress that Peru has achieved towards malaria elimination.

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.

Humanized mice for investigating sustained Plasmodium vivax blood-stage infections and transmission

Plasmodium vivax is the most widespread human malaria parasite. Due to the presence of extravascular reservoirs and relapsing infections from dormant liver stages, P. vivax is particularly difficult to control and eliminate. Experimental research is hampered by the inability to maintain P. vivax cultures in vitro, due to its tropism for immature red blood cells (RBCs). Here, we describe a new humanized mice model that can support efficient human erythropoiesis and maintain long-lasting multiplication of inoculated cryopreserved P. vivax parasites and their sexual differentiation, including in bone marrow. Mature gametocytes were transmitted to Anopheles mosquitoes, which led to the formation of salivary gland sporozoites. Importantly, blood-stage P. vivax parasites were maintained after the secondary transfer of fresh or frozen infected bone marrow cells to naïve chimeras. This model provides a unique tool for investigating, in vivo, the biology of intraerythrocytic P. vivax.

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.

Cohort profile: the Mâncio Lima cohort study of urban malaria in Amazonian Brazil

PURPOSE

This population-based open cohort study aims to investigate biological and sociodemographic drivers of malaria transmission in the main urban hotspot of Amazonian Brazil.

PARTICIPANTS

Nearly 20% of the households in the northwestern town of Mâncio Lima were randomly selected and 2690 participants were enrolled since April 2018. Sociodemographic, housing quality, occupational, behavioural and morbidity information and travel histories were collected during consecutive study visits. Blood samples from participants>3 months old were used for malaria diagnosis and human genetic studies; samples from participants with laboratory-confirmed malaria have been cryopreserved for genetic and phenotypic characterisation of parasites. Serology was introduced in 2020 to measure the prevalence and longevity of SARS-CoV-2 IgG antibodies.

FINDINGS TO DATE

Malaria prevalence rates were low (up to 1.0% for Plasmodium vivax and 0.6% for P. falciparum) during five consecutive cross-sectional surveys between April-May 2018 and October-November 2020; 63% of infections diagnosed by microscopy were asymptomatic. Malaria risk is heterogeneously distributed, with 20% study participants contributing 86% of the overall burden of P. vivax infection. Adult males are at greatest risk of infection and human mobility across the urban-rural interface may contribute to sustained malaria transmission. Local P. vivax parasites are genetically diverse and fragmented into discrete inbred lineages that remain stable across space and time.

FUTURE PLANS

Two follow-up visits, with similar study protocols, are planned in 2021. We aim to identify high-risk individuals that fuel onwards malaria transmission and represent a priority target for more intensive and effective control interventions.

TRIAL REGISTRATION NUMBER

NCT03689036.

Ongoing host-shift speciation in Plasmodium simium

Plasmodium simium, a malaria parasite that infects platyrrhine monkeys and humans in the New World, is nearly identical to Plasmodium vivax. Recent genomic comparative analyses of these sister species have identified elevated divergence in a gene that may underlie P. simium adaptation to non-human primates during its gradual speciation process.

Increasing proportions of relapsing parasite species among imported malaria in China's Guangxi Province from Western and Central Africa

BACKGROUND

Travel-related malaria in non-endemic areas returning from endemic areas presents important challenges to diagnosis and treatment. Imported malaria to newly malaria-free countries poses further threats of malaria re-introduction and potential resurgence. For those traveling to places with high Plasmodium falciparum prevalence, prophylaxis against this parasite is recommended, whereas causal prophylaxis against relapsing malaria is often overlooked.

METHODS

We analyzed a cluster of imported malaria among febrile patients in Shanglin County, Guangxi Province, China, who had recent travel histories to Western and Central Africa. Malaria was diagnosed by microscopy and subsequently confirmed by species- and subspecies-specific PCR. Plasmodium vivax was genotyped using a barcode consisting of 42 single nucleotide polymorphisms.

RESULTS

Investigations of 344 PCR-confirmed malaria cases revealed that in addition to Plasmodium falciparum being the major parasite species, the relapsing parasites Plasmodium ovale and P. vivax accounted for ~40% of these imported cases. Of the 114 P. ovale infections, 65.8% and 34.2% were P. ovale curtisi and P. ovale wallikeri, respectively, with the two subspecies having a ~2:1 ratio in both Western and Central Africa. Phylogenetic analysis of 14 P. vivax isolates using a genetic barcode demonstrated that 11 formed a distinct clade from P. vivax populations from Eastern Africa.

CONCLUSION

This study provides support for active P. vivax transmission in areas with the predominant Duffy-negative blood group. With relapsing malaria making a substantial proportion of the imported malaria, causal prophylaxis should be advocated to travelers with a travel destination to Western and Central Africa.

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

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

Plasmodium simium: population genomics reveals the origin of a reverse zoonosis

The population history of Plasmodium simium, which causes malaria in sylvatic Neotropical monkeys and humans along the Atlantic Coast of Brazil, remains disputed. Genetically diverse P. vivax populations from various sources, including the lineages that founded the species P. simium, are thought to have arrived in the Americas in separate migratory waves. However, here we find a minimal genome-level differentiation between P. simium and present-day New World P. vivax isolates, consistent with their common geographic origin and subsequent divergence on this continent. The meagre genetic diversity in P. simium samples from humans and monkeys implies a recent transfer from humans to non-human primates - a unique example of malaria as a reverse zoonosis of public health significance. Likely genomic signatures of P. simium adaptation to new hosts include the deletion of >40% of a key erythrocyte invasion ligand, PvRBP2a, which may have favored more efficient simian host cell infection.

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.

Monitoring Plasmodium vivax resistance to antimalarials: Persisting challenges and future directions

Emerging antimalarial drug resistance may undermine current efforts to control and eliminate Plasmodium vivax, the most geographically widespread yet neglected human malaria parasite. Endemic countries are expected to assess regularly the therapeutic efficacy of antimalarial drugs in use in order to adjust their malaria treatment policies, but proper funding and trained human resources are often lacking to execute relatively complex and expensive clinical studies, ideally complemented by ex vivo assays of drug resistance. Here we review the challenges for assessing in vivo P. vivax responses to commonly used antimalarials, especially chloroquine and primaquine, in the presence of confounding factors such as variable drug absorption, metabolism and interaction, and the risk of new infections following successful radical cure. We introduce a simple modeling approach to quantify the relative contribution of relapses and new infections to recurring parasitemias in clinical studies of hypnozoitocides. Finally, we examine recent methodological advances that may render ex vivo assays more practical and widely used to confirm P. vivax drug resistance phenotypes in endemic settings and review current approaches to the development of robust genetic markers for monitoring chloroquine resistance in P. vivax populations.

Barriers to implementing antimicrobial stewardship programs in three low- and middle-income country tertiary care settings: findings from a multi-site qualitative study

BACKGROUND

Antimicrobial resistance has been named as one of the top ten threats to public health in the world. Hospital-based antimicrobial stewardship programs (ASPs) can help reduce antimicrobial resistance. The purpose of this study was to determine perceived barriers to the development and implementation of ASPs in tertiary care centers in three low- and middle-income countries (LMICs).

METHODS

Interviews were conducted with 45 physicians at tertiary care hospitals in Sri Lanka (n = 22), Kenya (12), and Tanzania (11). Interviews assessed knowledge of antimicrobial resistance and ASPs, current antimicrobial prescribing practices, access to diagnostics that inform antimicrobial use, receptiveness to ASPs, and perceived barriers to implementing ASPs. Two independent reviewers coded the interviews using principles of applied thematic analysis, and comparisons of themes were made across the three sites.

RESULTS

Barriers to improving antimicrobial prescribing included prohibitively expensive antimicrobials, limited antimicrobial availability, resistance to changing current practices regarding antimicrobial prescribing, and limited diagnostic capabilities. The most frequent of these barriers in all three locations was limited drug availability. Many physicians in all three sites had not heard of ASPs before the interviews. Improved education was a suggested component of ASPs at all three sites. The creation of guidelines was also recommended, without prompting, by interviewees at all three sites. Although most participants felt microbiological results were helpful in tailoring antibiotic courses, some expressed distrust of laboratory culture results. Biomarkers like erythrocyte sedimentation rate and c-reactive protein were not felt to be specific enough to guide antimicrobial therapy. Despite limited or no prior knowledge of ASPs, most interviewees were receptive to implementing protocols that would include documentation and consultation with ASPs regarding antimicrobial prescribing.

CONCLUSIONS

Our study highlighted several important barriers to implementing ASPs that were shared between three tertiary care centers in LMICs. Improving drug availability, enhancing availability of and trust in microbiologic data, creating local guidelines, and providing education to physicians regarding antimicrobial prescribing are important steps that could be taken by ASPs in these facilities.

Heterogeneity in response to serological exposure markers of recent Plasmodium vivax infections in contrasting epidemiological contexts

BACKGROUND

Antibody responses as serological markers of Plasmodium vivax infection have been shown to correlate with exposure, but little is known about the other factors that affect antibody responses in naturally infected people from endemic settings. To address this question, we studied IgG responses to novel serological exposure markers (SEMs) of P. vivax in three settings with different transmission intensity.

METHODOLOGY

We validated a panel of 34 SEMs in a Peruvian cohort with up to three years' longitudinal follow-up using a multiplex platform and compared results to data from cohorts in Thailand and Brazil. Linear regression models were used to characterize the association between antibody responses and age, the number of detected blood-stage infections during follow-up, and time since previous infection. Receiver Operating Characteristic (ROC) analysis was used to test the performance of SEMs to identify P. vivax infections in the previous 9 months.

PRINCIPAL FINDINGS

Antibody titers were associated with age, the number of blood-stage infections, and time since previous P. vivax infection in all three study sites. The association between antibody titers and time since previous P. vivax infection was stronger in the low transmission settings of Thailand and Brazil compared to the higher transmission setting in Peru. Of the SEMs tested, antibody responses to RBP2b had the highest performance for classifying recent exposure in all sites, with area under the ROC curve (AUC) = 0.83 in Thailand, AUC = 0.79 in Brazil, and AUC = 0.68 in Peru.

CONCLUSIONS

In low transmission settings, P. vivax SEMs can accurately identify individuals with recent blood-stage infections. In higher transmission settings, the accuracy of this approach diminishes substantially. We recommend using P. vivax SEMs in low transmission settings pursuing malaria elimination, but they are likely to be less effective in high transmission settings focused on malaria control.

Molecular surveillance of the Plasmodium vivax multidrug resistance 1 gene in Peru between 2006 and 2015

BACKGROUND

The high incidence of Plasmodium vivax infections associated with clinical severity and the emergence of chloroquine (CQ) resistance has posed a challenge to control efforts aimed at eliminating this disease. Despite conflicting evidence regarding the role of mutations of P. vivax multidrug resistance 1 gene (pvmdr1) in drug resistance, this gene can be a tool for molecular surveillance due to its variability and spatial patterns.

METHODS

Blood samples were collected from studies conducted between 2006 and 2015 in the Northern and Southern Amazon Basin and the North Coast of Peru. Thick and thin blood smears were prepared for malaria diagnosis by microscopy and PCR was performed for detection of P. vivax monoinfections. The pvmdr1 gene was subsequently sequenced and the genetic data was used for haplotype and diversity analysis.

RESULTS

A total of 550 positive P. vivax samples were sequenced; 445 from the Northern Amazon Basin, 48 from the Southern Amazon Basin and 57 from the North Coast. Eight non-synonymous mutations and three synonymous mutations were analysed in 4,395 bp of pvmdr1. Amino acid changes at positions 976F and 1076L were detected in the Northern Amazon Basin (12.8%) and the Southern Amazon Basin (4.2%) with fluctuations in the prevalence of both mutations in the Northern Amazon Basin during the course of the study that seemed to correspond with a malaria control programme implemented in the region. A total of 13 pvmdr1 haplotypes with non-synonymous mutations were estimated in Peru and an overall nucleotide diversity of π = 0.00054. The Northern Amazon Basin was the most diverse region (π = 0.00055) followed by the Southern Amazon and the North Coast (π = 0.00035 and π = 0.00014, respectively).

CONCLUSION

This study showed a high variability in the frequencies of the 976F and 1076L polymorphisms in the Northern Amazon Basin between 2006 and 2015. The low and heterogeneous diversity of pvmdr1 found in this study underscores the need for additional research that can elucidate the role of this gene on P. vivax drug resistance as well as in vitro and clinical data that can clarify the extend of CQ resistance in Peru.

Balancing selection and high genetic diversity of Plasmodium vivax circumsporozoite central region in parasites from Brazilian Amazon and Rio de Janeiro Atlantic Forest

Circumsporozoite protein (CSP) is the primary pre-erythrocytic vaccine target in Plasmodium species. Knowledge about their genetic diversity can help predict vaccine efficacy and the spread of novel parasite variants. Thus, we investigated pvcsp gene polymorphisms in 219 isolates (136 from Brazilian Amazon [BA], 71 from Rio de Janeiro Atlantic Forest [AF], and 12 from non-Brazilian countries [NB]). Forty-eight polymorphic sites were detected, 46 in the central repeat region (CR), and two in the C-terminal region. Also, the CR presents InDels and a variable number of repeats. All samples correspond to the VK210 variant, and 24 VK210 subtypes based on CR. Nucleotide diversity (π = 0.0135) generated a significant number of haplotypes (168) with low genetic differentiation between the Brazilian regions (Fst = 0.208). The haplotype network revealed similar distances among the BA and AF regions. The linkage disequilibrium indicates that recombination does not seem to be acting in diversity, reinforcing natural selection's role in accelerating adaptive evolution. The high diversity (low Fst) and polymorphism frequencies could be indicators of balancing selection. Although malaria in BA and AF have distinct vector species and different host immune pressures, consistent genetic signature was found in two regions. The immunodominant B-cell epitope mapped in the CR varies from seven to 19 repeats. The CR T-cell epitope is conserved only in 39 samples. Concerning to C-terminal region, the Th2R epitope presented nonsynonymous SNP only in 6% of Brazilian samples, and the Th3R epitope remained conserved in all studied regions. We conclude that, although the uneven distribution of alleles may jeopardize the deployment of vaccines directed to a specific variable locus, a unique vaccine formulation could protect populations in all Brazilian regions.

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.

Plasmodium vivax is the most geographically widespread parasite species that causes malaria in humans. Although it occurs in Africa as a member of a mix of Plasmodium species, P. vivax is dominant in other parts of the world outside of Africa (e.g., Brazil). It was previously thought that most African populations were immune to P. vivax infections due to the absence of Duffy antigen chemokine receptor (DARC) gene expression required for erythrocyte invasion. However, several recent reports have indicated the emergence and potential spread of P. vivax across human populations in Africa. Compared to Southeast Asia and South America where P. vivax is highly endemic, data on polymorphisms in erythrocyte binding gene candidates of P. vivax from Africa is limited. Filling this knowlege gap is critical for identifying functional genes in erythrocyte invasion, biomarkers for tracking the P. vivax isolates from Africa, as well as potential gene targets for vaccine development. This paper examined the level of genetic polymorphisms in a panel of 43 potential erythrocyte binding protein genes based on whole genome sequences and described transmission patterns of P. vivax infections from different study sites in Ethiopia based on the genetic variants. Our analyses showed that chromosomes 9 and 10 of the P. vivax genomes isolated in Ethiopia had the most high-quality genetic polymorphisms. Among all erythrocyte binding protein gene candidates, the merozoite surface proteins 1 and merozoite surface protein 3 showed high levels of polymorphism. MAEBL and MSP3.8 related to immunogenicity and erythrocyte binding function were detected with significant signals of positive selection. The expansion of the Duffy binding protein and merozoite surface protein 3 gene copies was an independent process among the P. vivax lineages in Ethiopia. Various levels of gene flow were observed even at a smaller geographical scale. Our study provided baseline data for future comparison with P. vivax in Duffy negative individuals and help develop a panel of genetic markers that are informative 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.

New Plasmodium vivax Genomes From the China-Myanmar Border

Plasmodium vivax is increasingly the dominant species of malaria in the Greater Mekong Subregion (GMS), which is pursuing regional malaria elimination. P. vivax lineages in the GMS are poorly characterized. Currently, P. vivax reference genomes are scarce due to difficulties in culturing the parasite and lack of high-quality samples. In addition, P. vivax is incredibly diverse, necessitating the procurement of reference genomes from different geographical regions. Here we present four new P. vivax draft genomes assembled de novo from clinical samples collected in the China-Myanmar border area. We demonstrate comparable length and content to existing genomes, with the majority of structural variation occurring around subtelomeric regions and exported proteins, which we corroborated with detection of copy number variations in these regions. We predicted peptides from all PIR gene subfamilies, except for PIR D. We confirmed that proteins classically labeled as PIR D family members are not identifiable by PIR motifs, and actually bear stronger resemblance to DUF (domain of unknown function) family DUF3671, potentially pointing to a new, closely related gene family. Further, phylogenetic analyses of MSP7 genes showed high variability within the MSP7-B family compared to MSP7-A and -C families, and the result was comparable to that from whole genome analyses. The new genome assemblies serve as a resource for studying P. vivax within the GMS.

Human Plasmodium vivax diversity, population structure and evolutionary origin

More than 200 million malaria clinical cases are reported each year due to Plasmodium vivax, the most widespread Plasmodium species in the world. This species has been neglected and understudied for a long time, due to its lower mortality in comparison with Plasmodium falciparum. A renewed interest has emerged in the past decade with the discovery of antimalarial drug resistance and of severe and even fatal human cases. Nonetheless, today there are still significant gaps in our understanding of the population genetics and evolutionary history of P. vivax, particularly because of a lack of genetic data from Africa. To address these gaps, we genotyped 14 microsatellite loci in 834 samples obtained from 28 locations in 20 countries from around the world. We discuss the worldwide population genetic structure and diversity and the evolutionary origin of P. vivax in the world and its introduction into the Americas. This study demonstrates the importance of conducting genome-wide analyses of P. vivax in order to unravel its complex evolutionary history.

Among the five Plasmodium species infecting humans, P. vivax is the most prevalent parasite outside Africa. To date, there has been less research on this species than for Plasmodium falciparum, a more lethal species, principally because of the lack of an in vitro culture system and also because P. vivax is considered relatively benign. Nevertheless, P. vivax is responsible for severe and incapacitating clinical symptoms with significant effects on human health. The emergence of new drug resistance and the discovery of severe and even fatal cases due to P. vivax question the benign status of P. vivax malaria. In recent years, there has been increased interest in characterizing the distribution of genetic variation in P. vivax. However, these studies either generated genetic information from a regional geographic scale or combine genetic datasets generated in different molecular platforms, which is known to generate biased results. In this study, we used a single genotyping platform to genotype 14 microsatellite markers in 834 samples of P. vivax obtained from 28 locations in 20 countries from around the world, including several populations from East and West Africa. We discuss the worldwide population genetic structure and the evolutionary origins of P. vivax, as well as its introduction into the Americas.

Plasmodium vivax transcriptional profiling of low input cryopreserved isolates through the intraerythrocytic development cycle

Approximately one-third of the global population is at risk of Plasmodium vivax infection, and an estimated 7.51 million cases were reported in 2017. Although, P. vivax research is currently limited by the lack of a robust continuous in vitro culture system for this parasite, recent work optimizing short-term ex vivo culture of P. vivax from cryopreserved isolates has facilitated quantitative assays on synchronous parasites. Pairing this improved culture system with low-input Smart-seq2 RNAseq library preparation, we sought to determine whether transcriptional profiling of P. vivax would provide insight into the differential survival of parasites in different culture media. To this end we probed the transcriptional signature of three different ex vivo P. vivax samples in four different culture media using only 1000 cells for each time point taken during the course of the intraerythrocytic development cycle (IDC). Using this strategy, we achieved similar quality transcriptional data to previously reported P. vivax transcriptomes. We found little effect with varying culture media on parasite transcriptional signatures, identified many novel gametocyte-specific genes from transcriptomes of FACS-isolated gametocytes, and determined invasion ligand expression in schizonts in biological isolates and across the IDC. In total, these data demonstrate the feasibility and utility of P. vivax RNAseq-based transcriptomic studies using minimal biomass input to maximize experimental capacity.

Plasmodium vivax is the most prevalent malaria-causing parasite species outside of Sub-Saharan Africa and has many unique and poorly understood biological characteristics that make it particularly challenging to study and combat. Transcriptomic profiling of P. vivax under various conditions has the potential to unlock new experimental abilities and aid in elucidating biology and the development of clinical interventions. However, a lack of a robust in vitro culture system for this parasite has restricted transcriptomic studies to researchers with timely access to fresh human isolates from clinics, which often are in resource-poor settings, as well as nearby, well-equipped laboratories for sample processing. This study aimed to gain insight into the differential survival of P. vivax in various culture media from the parasites transcriptional signature in each media. By implementing low-input RNA library preparation strategies, this study obtains robust transcriptomic data at various parasite development stages and in different culture conditions from just 1000 FACS-purified, P. vivax-infected erythrocytes from viable cryopreserved patient isolates. With these data, we find culture media has little effect on transcriptional profile, we characterize invasion ligand expression across intraerythrocytic development and between clinical isolates, and we define the transcriptome of sexual, transmissible stages of the P. vivax parasite. These results highlight the establishment and utility of a powerful platform for studying the transcriptomic biology of this particularly challenging parasite.

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

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

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.

Plasmodium Genomics and Genetics: New Insights into Malaria Pathogenesis, Drug Resistance, Epidemiology, and Evolution

Protozoan Plasmodium parasites are the causative agents of malaria, a deadly disease that continues to afflict hundreds of millions of people every year. Infections with malaria parasites can be asymptomatic, with mild or severe symptoms, or fatal, depending on many factors such as parasite virulence and host immune status. Malaria can be treated with various drugs, with artemisinin-based combination therapies (ACTs) being the first-line choice. Recent advances in genetics and genomics of malaria parasites have contributed greatly to our understanding of parasite population dynamics, transmission, drug responses, and pathogenesis. However, knowledge gaps in parasite biology and host-parasite interactions still remain. Parasites resistant to multiple antimalarial drugs have emerged, while advanced clinical trials have shown partial efficacy for one available vaccine. Here we discuss genetic and genomic studies of Plasmodium biology, host-parasite interactions, population structures, mosquito infectivity, antigenic variation, and targets for treatment and immunization. Knowledge from these studies will advance our understanding of malaria pathogenesis, epidemiology, and evolution and will support work to discover and develop new medicines and vaccines.

Plasmodium knowlesi as a model system for characterising Plasmodium vivax drug resistance candidate genes

Plasmodium vivax causes the majority of malaria outside Africa, but is poorly understood at a cellular level partly due to technical difficulties in maintaining it in in vitro culture conditions. In the past decades, drug resistant P. vivax parasites have emerged, mainly in Southeast Asia, but while some molecular markers of resistance have been identified, none have so far been confirmed experimentally, which limits interpretation of the markers, and hence our ability to monitor and control the spread of resistance. Some of these potential markers have been identified through P. vivax genome-wide population genetic analyses, which highlighted genes under recent evolutionary selection in Southeast Asia, where chloroquine resistance is most prevalent. These genes could be involved in drug resistance, but no experimental proof currently exists to support this hypothesis. In this study, we used Plasmodium knowlesi, the most closely related species to P. vivax that can be cultured in human erythrocytes, as a model system to express P. vivax genes and test for their role in drug resistance. We adopted a strategy of episomal expression, and were able to express fourteen P. vivax genes, including two allelic variants of several hypothetical resistance genes. Their expression level and localisation were assessed, confirming cellular locations conjectured from orthologous species, and suggesting locations for several previously unlocalised proteins, including an apical location for PVX_101445. These findings establish P. knowlesi as a suitable model for P. vivax protein expression. We performed chloroquine and mefloquine drug assays, finding no significant differences in drug sensitivity: these results could be due to technical issues, or could indicate that these genes are not actually involved in drug resistance, despite being under positive selection pressure in Southeast Asia. These data confirm that in vitro P. knowlesi is a useful tool for studying P. vivax biology. Its close evolutionary relationship to P. vivax, high transfection efficiency, and the availability of markers for colocalisation, all make it a powerful model system. Our study is the first of its kind using P. knowlesi to study unknown P. vivax proteins and investigate drug resistance mechanisms.

Monitoring the Efficacy of Chloroquine-Primaquine Therapy for Uncomplicated Plasmodium vivax Malaria in the Main Transmission Hot Spot of Brazil

Emerging Plasmodium vivax resistance to chloroquine (CQ) may undermine malaria elimination efforts in South America. CQ-resistant P. vivax has been found in the major port city of Manaus but not in the main malaria hot spots across the Amazon Basin of Brazil, where CQ is routinely coadministered with primaquine (PQ) for radical cure of vivax malaria. Here we randomly assigned 204 uncomplicated vivax malaria patients from Juruá Valley, northwestern Brazil, to receive either sequential (arm 1) or concomitant (arm 2) CQ-PQ treatment. Because PQ may synergize the blood schizontocidal effect of CQ and mask low-level CQ resistance, we monitored CQ-only efficacy in arm 1 subjects, who had PQ administered only at the end of the 28-day follow-up. We found adequate clinical and parasitological responses in all subjects assigned to arm 2. However, 2.2% of arm 1 patients had microscopy-detected parasite recrudescences at day 28. When PCR-detected parasitemias at day 28 were considered, response rates decreased to 92.1% and 98.8% in arms 1 and 2, respectively. Therapeutic CQ levels were documented in 6 of 8 recurrences, consistent with true CQ resistance in vivo In contrast, ex vivo assays provided no evidence of CQ resistance in 49 local P. vivax isolates analyzed. CQ-PQ coadministration was not found to potentiate the antirelapse efficacy of PQ over 180 days of surveillance; however, we suggest that larger studies are needed to examine whether and how CQ-PQ interactions, e.g., CQ-mediated inhibition of PQ metabolism, modulate radical cure efficacy in different P. vivax-infected populations. (This study has been registered at ClinicalTrials.gov under identifier NCT02691910.).

Malaria Situation in Latin America and the Caribbean: Residual and Resurgent Transmission and Challenges for Control and Elimination

Despite recent progress toward malaria elimination in Latin America and the Caribbean, with an overall 62% decrease in incidence between 2000 and 2015, malaria remains endemic to 21 countries and territories in the region, where 120 million people are exposed to some risk of infection. Here we review recent epidemiologic trends, highlight current challenges, and briefly discuss the relative role of traditional and novel strategies for better malaria control and elimination across the continent.

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

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

Bone Marrow Is a Major Parasite Reservoir in Plasmodium vivax Infection

Plasmodium vivax causes heavy burdens of disease across malarious regions worldwide. Mature P. vivax asexual and transmissive gametocyte stages occur in the blood circulation, and it is often assumed that accumulation/sequestration in tissues is not an important phase in their development. Here, we present a systematic study of P. vivax stage distributions in infected tissues of nonhuman primate (NHP) malaria models as well as in blood from human infections. In a comparative analysis of the transcriptomes of P. vivax and Plasmodium falciparum blood-stage parasites, we found a conserved cascade of stage-specific gene expression despite the greatly different gametocyte maturity times of these two species. Using this knowledge, we validated a set of conserved asexual- and gametocyte-stage markers both by quantitative real-time PCR and by antibody assays of peripheral blood samples from infected patients and NHP (Aotus sp.). Histological analyses of P. vivax parasites in organs of 13 infected NHP (Aotus and Saimiri species) demonstrated a major fraction of immature gametocytes in the parenchyma of the bone marrow, while asexual schizont forms were enriched to a somewhat lesser extent in this region of the bone marrow as well as in sinusoids of the liver. These findings suggest that the bone marrow is an important reservoir for gametocyte development and proliferation of malaria parasites.IMPORTANCEPlasmodium vivax malaria continues to cause major public health burdens worldwide. Yet, significant knowledge gaps in the basic biology and epidemiology of P. vivax malaria remain, largely due to limited available tools for research and diagnostics. Here, we present a systematic examination of tissue sequestration during P. vivax infection. Studies of nonhuman primates and malaria patients revealed enrichment of developing sexual stages (gametocytes) and mature replicative stages (schizonts) in the bone marrow and liver, relative to those present in peripheral blood. Identification of the bone marrow as a major P. vivax tissue reservoir has important implications for parasite diagnosis and treatment.

Enhanced Ex Vivo Plasmodium vivax Intraerythrocytic Enrichment and Maturation for Rapid and Sensitive Parasite Growth Assays

Plasmodium vivax chloroquine resistance has been documented in nearly every region where this malaria-causing parasite is endemic. Unfortunately, P. vivax resistance surveillance and drug discovery are challenging due to the low parasitemias of patient isolates and poor parasite survival through ex vivo maturation that reduce the sensitivity and scalability of current P. vivax antimalarial assays. Using cryopreserved patient isolates from Brazil and fresh patient isolates from India, we established a robust enrichment method for P. vivax parasites. We next performed a medium screen for formulations that enhance ex vivo survival. Finally, we optimized an isotopic metabolic labeling assay for measuring P. vivax maturation and its sensitivity to antimalarials. A KCl Percoll density gradient enrichment method increased parasitemias from small-volume ex vivo isolates by an average of >40-fold. The use of Iscove's modified Dulbecco's medium for P. vivax ex vivo culture approximately doubled the parasite survival through maturation. Coupling these with [³H]hypoxanthine metabolic labeling permitted sensitive and robust measurements of parasite maturation, which was used to measure the sensitivities of Brazilian P. vivax isolates to chloroquine and several novel antimalarials. These techniques can be applied to rapidly and robustly assess the P. vivax isolate sensitivities to antimalarials for resistance surveillance and drug discovery.

CD39 and immune regulation in a chronic helminth infection: The puzzling case of Mansonella ozzardi

Background

Chronic helminth infections typically induce an immunoregulatory environment, with markedly reduced immune responses to both parasite-specific and unrelated bystander antigens. Here we tested whether these changes are also observed in human infections with Mansonella ozzardi, a neglected filarial nematode widely distributed across Latin America.

Methods

CD4⁺ T cell populations from microfilaremic (Fil+) and uninfected (Fil-) inhabitants in M. ozzardi-endemic riverine communities in Brazil were characterized by flow cytometry analysis. Plasma concentrations of a wide range of cytokines and chemokines were measured. We examined whether M. ozzardi infection is associated with suppressed in vitro lymphoproliferative and inflammatory cytokine responses upon stimulation with filarial antigen, unrelated antigens or mitogens.

Principal findings/Conclusions

Fil+ subjects had lower plasma levels of selected inflammatory cytokines, such as TNF-α, IL-8, and IL-6, than their Fil- counterparts. However, we found no evidence for attenuated T-cell responses to filarial antigens or co-endemic pathogens, such as malaria parasites and Toxoplasma gondii. CD4⁺ T cells expressing CD39, an ectonucleosidase involved in the generation of the anti-inflammatory molecule adenosine, were increased in frequency in Fil+ subjects, compared to uninfected controls. Significantly, such an expansion was directly proportional to microfilarial loads. Surprisingly, CD39 blocking with a neutralizing antibody suppressed antigen-driven lymphoproliferation in vitro, while decreasing inflammatory cytokine responses, in Fil+ and Fil- individuals. These findings suggest that circulating CD4⁺ CD39⁺ T cells comprise subsets with both regulatory and stimulatory roles that contribute to the immune homeostasis in chronic M. ozzardi infection.

Helminth infections downregulate immunity and reduce host’s inflammatory responses, but the filarial nematode Mansonella ozzardi, which is widely distributed across Latin America, appears to represent an exception to this rule. We found similar lymphoproliferative responses to filarial and unrelated antigens and comparable regulatory cytokine responses in subjects harboring M. ozzardi microfilariae, compared to local uninfected controls. The proportion of CD4⁺ T cell subtypes expressing CD39 was significantly increased in infected subjects and correlated positively with their microfilarial density. However, antibody blocking of CD39, an ectonucleosidase involved in the synthesis of the immunosuppresive molecule adenosine, paradoxically reduced, rather than promoted, antigen-driven lymphoproliferation in vitro. We suggest that CD39⁺ CD4⁺ T cells circulating in microfilaremics comprise both regulatory and stimulatory cell subsets that are concomitantly expanded. The balance between these cell subsets with opposing regulatory functions may be crucial to maintain immune homeostasis during chronic M. ozzardi infections.

Human migration and the spread of malaria parasites to the New World

We examined the mitogenomes of a large global collection of human malaria parasites to explore how and when Plasmodium falciparum and P. vivax entered the Americas. We found evidence of a significant contribution of African and South Asian lineages to present-day New World malaria parasites with additional P. vivax lineages appearing to originate from Melanesia that were putatively carried by the Australasian peoples who contributed genes to Native Americans. Importantly, mitochondrial lineages of the P. vivax-like species P. simium are shared by platyrrhine monkeys and humans in the Atlantic Forest ecosystem, but not across the Amazon, which most likely resulted from one or a few recent human-to-monkey transfers. While enslaved Africans were likely the main carriers of P. falciparum mitochondrial lineages into the Americas after the conquest, additional parasites carried by Australasian peoples in pre-Columbian times may have contributed to the extensive diversity of extant local populations of P. vivax.